考虑颗粒破碎的粗粒土剪胀性统一本构模型

粗粒土作为无黏性散粒状材料具有状态依赖特性,土体的剪切特性受密度和应力水平影响。易破碎是粗粒土的另一个特点,颗粒破碎影响粗粒土的剪胀、内摩擦角、峰值强度和渗透系数。为了能够准确地描述粗粒土的应力-应变关系,采用初始状态参量描述粗粒土的内部状态,根据三轴试验数据建立考虑颗粒破碎耗能的应力-应变关系,采用相关联流动法则推导考虑颗粒破碎的粗粒土剪胀性“统一本构模型”,并建立初始状态参量与模型参数之间的关系。所建立的统一本构模型既考虑了颗粒破碎对剪胀、内摩擦角的影响,又考虑了剪切特性对土体初始状态的依赖。采用变异粒子群算法拟合试验曲线,确定模型参数。模型计算结果能够很好地拟合试验曲线。采用同一组参数对假定的初始状态进行模拟计算,计算结果表明,模型能够模拟不同初始密度和应力水平下粗粒土变形的一般规律。     

考虑颗粒破碎效应的堆石料静动力本构模型

为合理反映颗粒破碎对堆石料力学特性的影响,基于试验结果分析,得出了堆石料在压缩和剪切作用下的颗粒破碎特性规律。通过引入压缩破碎和剪切破碎的相关参数,借鉴已有本构模型的合理定义,吸收临界状态理论和边界面理论的优点,发展了考虑颗粒破碎和状态相关的堆石料静动力统一弹塑性本构模型,并阐述了模型参数的确定方法。该模型不仅能够反映堆石料在静力荷载作用下的低压剪胀、高压剪缩、应变软化和硬化等特性,还能够反映在循环荷载作用下应力-应变的滞回特性和残余变形的累积效应。最后为验证模型的合理性,分别对堆石料的静力三轴和循环三轴试验进行了数值模拟预测,结果表明：模型预测与试验数据吻合良好,所提出的本构模型能够合理地描述颗粒破碎对堆石料静动力变形特性的影响。     

考虑颗粒破碎对特征孔隙比影响的堆石体亚塑性本构模型

颗粒破碎是影响堆石体强度和变形特性的主要问题之一。相比于砂土,堆石料在较低的应力水平下就会发生严重的颗粒破碎,因此,在进行堆石体力学特性及本构模型研究时必须考虑颗粒破碎的影响。同时,堆石体在受力过程中孔隙比是变化的,而传统本构模型不能使用一组参数模拟不同孔隙比的同种材料。因此,以能够考虑应力水平和土体孔隙比影响的Gudehus-Bauer亚塑性本构模型为基础,考虑堆石体有别于砂土的孔隙比变化特征,提出了考虑堆石破碎的亚塑性本构模型。亚塑性理论是目前可最大限度地减少人为假定的一种本构理论,颗粒材料在不同特征应力路径下,破碎造成的过度变形量不同;但相同应力水平、不同特征应力路径下孔隙比已不满足Gudehus-Bauer亚塑性本构模型中提出的等比例变化规律。据此,结合考虑颗粒破碎的临界状态理论和堆石体常规三轴试验和循环加载试验结果,提出了考虑颗粒破碎堆石体特征孔隙比的表达式,并将其引入到Gudehus-Bauer亚塑性本构模型中,建立了考虑颗粒破碎的堆石体亚塑性本构模型,提出了模型参数的确定方法。与堆石体试验结果对比表明,该本构模型可以较好地模拟其力学与变形特性。     

堆石料的临界状态与考虑颗粒破碎的本构模型

高应力水平时堆石料的颗粒破碎对其强度和变形机制有重要影响。临界状态土力学理论对重塑土的应力-应变关系的描述较为成功,但目前颗粒破碎对堆石料的临界状态的影响及其数学描述鲜有研究。对堆石料进行了固结应力从0.4 MPa到4 MPa的18组固结排水和固结不排水常规三轴压缩试验,以及6组等向压缩试验。试验结果表明：在排水条件和不排水条件下,不同的固结应力试样都趋于临界状态;堆石料的临界状态在q-p′平面和e-lgp′平面均为非线性变化。基于此试验结果,通过引入状态参数,在广义塑性力学的理论框架下,建立了考虑颗粒破碎的堆石料本构模型,并给出了模型参数的确定方法。与长河坝料的试验进行了对比,结果表明所建议的本构模型可以较好地模拟堆石料从低围压0.4 MPa到高围压 3.5 MPa下的应力-应变特性     

考虑颗粒破碎影响的粗粒土本构模型

颗粒破碎直接改变了粗粒土本身结构,对粗粒土的剪胀、内摩擦角、峰值强度、渗透系数都会产生影响。为了能够准确地描述粗粒土的应力-应变关系,特别是高应力条件下出现显著颗粒破碎时的应力-应变关系,亟待建立考虑颗粒破碎的粗粒土本构模型。根据三轴试验数据,建立考虑颗粒破碎耗能的应力-应变关系,采用相关联流动法则导出考虑颗粒破碎的粗粒土本构模型。所建立的本构模型考虑了颗粒破碎对粗粒土剪胀、内摩擦角的影响。通过变异粒子群优化算法确定模型参数,拟合试验曲线。模型计算结果与试验曲线拟合较好,能够很好地描述粗粒土在不同围压下的体积剪胀、剪缩和应力硬化、软化现象。     

考虑颗粒破碎的冻结砂土非线性本构模型研究

颗粒破碎是粒状材料在高应力状态下的一种基本现象。为了研究冻结砂土中颗粒破碎对应力应变关系的影响,将冻结砂土视为复合颗粒材料,忽略冰的压融,考虑内摩擦角随应力状态的变化,构建一个适用于冻结砂土的考虑颗粒破碎的非线性本构模型。构建过程分为三步,首先是基于三轴剪切前后颗粒分析对冻结砂土颗粒破碎模式和产生机理进行探讨;其次是基于考虑颗粒破碎的能量平衡方程,对冻土在三轴剪切试验过程中的颗粒破碎耗能进行分析,结果表明颗粒破碎耗能随轴向应变呈双曲线变化趋势;最后应用考虑颗粒破碎的剪胀方程修正沈珠江三参数非线性模型中的体积切线模量ν_t,得到一个考虑颗粒破碎的非线性本构模型,模型参数可以通过单轴压缩试验和常规三轴试验确定。将原模型和修正后模型的计算结果与控制温度为-6℃,围压为1 MPa、4 MPa、6 MPa、8 MPa和10 MPa时冻结砂土的试验结果进行对比,结果表明该模型能够较好的模拟冻结砂土从低围压到高围压的应变软化特征与剪胀特征。      

颗粒破碎的三维离散元模拟研究

采用PFC等离散元方法研究岩土材料的颗粒破碎已经成为热点。采用考虑局部应力集中的点荷载破碎准则,利用阿波罗填充和膨胀法保证破碎前、后颗粒之间的种群平衡,并引入尺寸因子来表征不同粒径的颗粒强度。在此基础上,开展了石英砂、钙质砂和萨克拉门托河砂3种不同破碎难易程度材料的数值试验,并与室内试验结果进行对比。结果表明:建立的三维颗粒破碎模型能够很好地描述破碎难易程度不同的颗粒材料的压缩特性;考虑应力集中效应的点荷载破碎准则比基于平均应力Mohr-Coulomb理论的颗粒破碎准则更能真实地反应颗粒材料的破碎现象。同时,所建立的模型能够揭示破碎对颗粒材料各向异性消散和级配曲线演化的影响规律。     

考虑颗粒破碎引起级配演变的道砟边界面本构模型

道砟即使在低围压下也会产生颗粒破碎,从而改变其原始级配。级配作为道砟的一个重要物理特性,显著地影响着道砟本身物理、力学特性。建立可以反映颗粒破碎引起级配演变的本构模型至关重要。基于Einav的分形破碎指标和土体边界面塑性理论,通过分析道砟三轴试验结果,尝试建立了低围压下可以反映颗粒破碎引起级配演变的道砟边界面本构模型。与试验结果对比发现:模型可以较好地反映道砟在低围压下的应力-应变特性;同时,可以对加载过程中颗粒破碎引起的级配演变进行合理的预测。     

颗粒形状对粒状材料破碎演化规律及强度准则影响

为探讨颗粒形状对粒状材料的颗粒破碎演化规律及强度特征的影响,提出了一个新的粒状材料颗粒形状量化参数,设计了一种考虑三维颗粒形状的人工试样制备方法,随即进行了常规三轴压缩试验,并分析了颗粒破碎和强度特征,最终建立了一个二元介质强度准则,具体的研究成果为：建立颗粒形状量化参数——球形模量GM,在此基础上制备了5种不同形状的可破碎粒状材料三轴试样,并发现球形模量影响着粒状材料的三轴压缩强度特征;通过筛分确定试样的颗粒破碎情况,对试样的颗粒破碎演化规律和临界状态进行探讨,发现颗粒形状通过影响颗粒破碎规律而控制着宏观强度的非线性演化特征;以二元介质理论为基础,建立了考虑颗粒形状的可破碎粒状材料强度准则,并通过试验对其适用性进行了验证。     

亚塑性临界状态的本构模型

给出了无粘性粒状材料的三维非线性应力应变和体积变化特点的本构模型。基于临界状态理论,提出了一种极限面的特殊类型,证明和它与试验数据逼近具有更大灵活性。标准本构方程因子表达式允许本构常量独立校准。通过对材料压力水平和密度表现出的相关性的统一描述,使昨这个基     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.     

Machine learning reveals the influences of grain morphology on grain crushing strength

Acta Geotechnica - Grain morphology has significant impacts on the mechanical behaviors of granular materials. However, its influences on grain breakage are still poorly understood due to the...     

Effect of grain crushing on shear localization in granular bodies during plane strain compression

This article deals with the effect of grain crushing on shear localization in granular materials during plane strain monotonic compression tests under constant lateral pressure. The grain diameter and the initial void ratio were stochastically distributed using a spatial correlation. To describe the mechanical behavior of cohesionless granular materials during a monotonic deformation path in plane strain compression, we used a micropolar hypoplastic constitutive model that is able to describe the salient properties of granular bodies including shear localization. The model was extended by introducing changes to the grain diameter with varying pressure using formulae from breakage mechanics proposed for crushable granulates. The initial void ratios and grain diameters took the form of correlated random spatial fields described by both symmetric and nonsymmetric random distributions using a homogeneous correlation function. The field realizations were generated with the help of an original conditional rejection method. A few representative samples of the random fields selected from the generated set were taken into account in numerical calculations. Copyright © 2011 John Wiley & Sons, Ltd.     

脆性颗粒材料的动态多尺度模型研究

脆性颗粒材料的多尺度模型一般包含微观尺度的基本粒子、细观尺度的颗粒和宏观尺度的颗粒堆积体3个尺度。基于离散元方法（DEM）构建多尺度模型,并将该模型应用于动态加载。首先,对多尺度模型所涉及的两种接触模型和两种黏结模型的参数进行分析,详细讨论微细观模型参数与宏观材料常数之间的联系。然后,选用Hertz-Mindlin接触模型[1]和平行键黏结模型,建造石英砂的动态多尺度模型。通过选择合适的强度和局部阻尼参数发现,模型宏细观尺度上的动态压缩响应与对石英砂的相关试验结果吻合很好。利用多尺度模型和选定的参数,探讨与动态加载密切相关的局部阻尼机制对多尺度模型各个尺度上力学响应的影响。结果表明,阻尼越大则颗粒材料对波的衰减能力越强,但过高的阻尼会使团簇强度和模型的宏观压缩曲线都表现出异常的加载速度效应（后者实际是阻尼引起的微惯性效应）。另外,高阻尼会过度衰减颗粒破碎过程产生的应力波,从而阻碍颗粒破碎。最后,应用改进的动态多尺度模型,对脆性颗粒材料的动态破碎特性进行研究,发现该模型不但能给出与试验相吻合的颗粒级配曲线,还能揭示出颗粒破碎过程中微裂纹分布的空间不均匀性,即颗粒破碎过程中波的产生机制和衰减机制相互作用导致的微裂纹聚团分布的现象。     

DEM investigation of the effect of intermediate principle stress on particle breakage of granular materials

This study investigates the influence of the intermediate principle stress on the particle breakage of granular materials. The crushable agglomerate method is applied to model soil particles and numerical true triaxial tests were carried out. The results show that particle breakage increases with increasing b value, the relationship of which follows an exponential function and agrees well with previous experimental results. More importantly, the study found that the relationship between particle breakage and total energy input is independent of the intermediate principle stress, which provides a good basis for the constitutive modeling of granular materials.     

Modeling for critical state line of granular soil with evolution of grain size distribution due to particle breakage

Determination of the critical state line(CSL)is important to characterize engineering properties of granular soils.Grain size distribution(GSD)has a significant influence on the location of CSL.The influence of particle breakage on the CSL is mainly attributed to the change in GSD due to particle breakage.However,GSD has not been properly considered in modeling the CSL with influence of particle breakage.This study aims to propose a quantitative model to determine the CSL considering the effect of GSD.We hypothesize that the change of critical state void ratio with respect to GSD is caused by the same mechanism that influences of the change of minimum void ratio with respect to GSD.Consequently,the particle packing model for minimum void ratio proposed by Chang et al.(2017)is extended to predict critical state void ratio.The developed model is validated by experimental results of CSLs for several types of granular materials.Then the evolution of GSD due to particle breakage is incorporated into the model.The model is further evaluated using the experimental results on rockfill material,which illustrates the applicability of the model in predicting CSL for granular material with particle breakage.     

颗粒级配对粒状材料不排水力学特性的影响

为了研究颗粒级配对粒状材料不排水力学特性的影响,采用两种形状自相似度很高的材料（玻璃球和Hostun砂）进行了一系列传统三轴不排水压缩试验。对于每种材料,保持相近的相对密实度制作了6个不同颗粒级配（不均匀系数 1.1～20）的试样,然后进行剪切。试验结果表明,颗粒级配对粒状材料的不排水力学性能存在明显影响：两种材料试样的偏应力水平随着不均匀系数的增加而降低,抗剪强度随着材料的不均匀系数 的增加而降低,并在 5后呈现出趋稳势态。基于不排水剪切试验中颗粒集合体二阶功的演化规律,对材料受荷载过程中出现不稳定性的条件进行分析。结果表明,在相同密实度下,粒状材料的静态液化潜能随着材料不均匀系数 的提高而增强,高不均匀系数的材料变得不稳定。     

Modelling size effect on rock aggregates strength using a DEM bonded-cell model

Empirical evidence has shown that particle breakage affects the mechanical behaviour of granular materials. The source of this mechanism takes place at the particle scale, and the main consequence on the macromechanical behaviour is increasing compressibility. Due to the inverse correlation between particle size and particle crushing strength, coarse rockfill materials are particularly vulnerable to mechanical degradation due to particle breakage. However, such coarse materials do not fit in standard laboratory devices, and the alternative of large sample testing is usually unavailable or too expensive. Alternatively, recent works have proposed multi-scale approaches using the discrete element method (DEM) to carry out numerical testing of coarse crushable materials, although few studies have focused on size effects. This article presents the application of a DEM bonded-cell model to study particle size-strength correlation on angular rock aggregates. Each particle is modelled by a cluster of perfectly rigid polyhedral cells with Mohr–Coulomb contact law. Constant cell density within particles implies that the presence of potential fragmentation planes increases with size. Therefore, particle strength decreases with size. A comprehensive sensitivity analysis was carried out through 1477 particle crushing simulations in a given particle size. Based on published experimental data on calcareous rock aggregates, part of the simulations were used for calibration, and 97 additional simulations of a coarser size fraction were performed for validation. The results show a good agreement with the empirical data in terms of size effect and data scatter through Weibull statistics.

     

Effect of micro and macro parameters in 3D modeling of grain crushing

Acta Geotechnica - Coarse granular materials exhibit important grain breakage in civil engineering structures, making it more complicated to predict the settlement and collapse of structures. A...     

Suction and time effects in rockfill deformation

Ambient relative humidity controls the stress-strain-time behaviour of granular media exhibiting particle breakage. Breakage is the result of crack propagation inside loaded particles. The DEM model developed uses the results of subcritical crack propagation in brittle rock to predict breakage and the subsequent structural rearrangement. Two closed-form solutions at grain level, stress distribution under diametrically opposed forces and the mode I crack propagation velocity, were included in the DEM formulation. Crack propagation velocity depends explicitly on relative humidity. The resulting model incorporates comminution and splitting modes of particle breakage as well as arbitrary grain shapes by means of aggregation of spherical microparticles. The model was first validated against a large diameter suction-controlled oedometer test on hard limestone crushed gravel. Direct tests on contact properties (local stiffness and friction) helped to validate the model. The model was then used to predict the response of large diameter suction-controlled triaxial tests on gravels having different (uniform) initial sizes. Time delayed deformations and wetting-induced collapse deformations under maintained load are also a natural outcome of the model. They were shown to follow quite accurately the experimental results.