胶结砂土强度特性的真三轴试验离散元分析

采用三维离散元法从宏微观尺度研究胶结砂土在真三轴试验条件下的强度特性。首先,将一个已有的三维胶结接触模型引入到离散元程序中,对胶结砂土数值试样进行真三轴试验离散元模拟(平均应力p保持不变,包括结构屈服前p=100 kPa和屈服后p=800 kPa);然后,对比分析离散元模拟结果和室内试验结果,并验证已有的胶结材料破坏准则;最后,从微观尺度对强度特征给予机制性分析。结果表明:当p小于结构屈服应力py时,应力–应变关系表现为显著的应变软化,当ppy时,应力–应变关系则呈现为应变硬化现象,且均受中主应力系数b的影响。归一化的峰值强度q_(peak)随b增加而逐渐减小,该变化规律与已有室内试验结果符合较好。Modified Lade-Duncan准则可以较好地预测p=100,800 kPa时的峰值强度,Extended SMP准则能较好地预测p=100 kPa时的峰值强度,但高估了p=800 kPa时的峰值强度。在微观尺度上,胶结试样的宏观力学特性主要受强接触点控制;随b增加,颗粒间法向接触力分布逐渐不均匀,导致胶结材料强度降低。     

不同沉积方向各向异性结构性砂土离散元力学特性分析

为探究沉积方向对各向异性结构性砂土力学特性的影响,首先采用椭圆颗粒生成两种不同沉积方向(水平与竖直)的各向异性净砂样,其次,引入一个考虑胶结厚度影响的微观胶接触模型从而生成各向异性结构性砂土。最后,对两种各向异性结构性砂土试样进行双轴压缩试验,并将水平沉积试样试验结果与室内试验结果对比验证该模型的可行性。同时,将两种不同方向试样的试验结果进行对比以探究沉积方向的影响。结果表明:两种试样应力–应变关系均呈软化及剪胀现象,水平沉积试样峰值偏应力较竖直沉积试样大,而二者残余阶段偏应力无明显差别;水平沉积试样临界孔隙比较竖直沉积试样大;胶结破坏速率及胶结破坏率变化与宏观力学特性变化相对应,且变化规律基本相同;水平方向沉积试样的胶结接触主方向始终保持竖直而竖直沉积试样的胶结接触主方向始终保持水平;在水平沉积试样中始终为水平分布颗粒长轴分布主方向,而竖直沉积试样中颗粒长轴逐渐向各向同性分布靠近。     

考虑流固耦合作用的深基坑有限元分析

基于比奥固结理论和修正剑桥模型,利用流固耦合模型对复杂地质条件下深基坑降水开挖过程中基坑的时间效应进行研究,并以地铁车站深基坑为例,对基坑降水开挖过程中的渗流固结进行了有限元计算。计算结果揭示了基坑渗流场分布及基坑水土压力和地下连续墙变形、周围地面沉降及基坑底部回弹的分布规律;并将计算结果与现场实测数据进行了对比,发现考虑流固耦合的计算结果与实测数据较为接近,能够比较真实地反映基坑渗流与变形特征。     

理想胶结砂土力学特性及剪切带形成的离散元分析

根据近期胶结铝棒接触力学特性的实测结果,提炼出用于模拟胶结砂土粒间胶结作用的胶结接触模型,并将该模型引入二维离散元商业软件PFC^2D。通过对不同胶结强度和不同围压下胶结砂土的平面应变双轴压缩试验的离散元模拟,分析了理想胶结砂土的宏观力学特性及其剪切带的形成规律。结果表明：相比同一孔隙比的无胶结试样,胶结试样具有更高的峰值强度、显著的应变软化和剪胀现象以及明显的剪切带,宏观力学特性与其胶结接触微观力学机理密切相关,模拟结果与已有室内试验结果具有规律上的一致性;由胶结试样内部的微观信息统计可知,胶结试样剪切带的形成一般在其峰值强度之后,且剪切带的形成是试样变形、胶结破坏、孔隙比、平均纯转动率和位移场等微观参量局部化的综合表现。     

不同胶结砂土力学特性及胶结破坏的离散元模拟

采用离散元法研究不同胶结砂土的宏观力学特性以及微观胶结破坏,其中一种胶结砂土胶结特性由离散元商业软件PFC2D中的胶结接触模型(Contact bond model)进行控制,另一种则采用蒋明镜等提出的无厚度改进胶结接触模型(改进的蒋氏模型)控制.首先,将改进的蒋氏模型引入PFC2D;其次,使用PFC2D对不同胶结强度和不同围压下的上述两种胶结砂土进行双轴压缩试验模拟,对比分析了两种不同胶结砂土的模拟结果.结果表明蒋氏试样(胶结特性由改进的蒋氏模型控制)应变软化和体积剪胀显著,峰值内摩擦角基本随胶结强度增大而增大;而PFC试样(胶结特性由Contact bond model控制)在高胶结强度时应变软化和体积剪胀性比较显著,低胶结强度时表现为应变硬化和体积剪缩,峰值内摩擦角随胶结强度增加而减小.离散元模拟结果在一定程度上与已有室内实测结果相符,蒋氏试样更能反映胶结砂土的主要力学特性.由胶结破坏微观信息统计可知,蒋氏试样中胶结点受拉破坏率远大于受剪破坏率,而PFC试样中两者相当,且蒋氏试样中的总胶结点破坏速率峰值要大于PFC试样.通过对改进的蒋氏模型参数分析可知,低围压条件下,试样宏观力学特性与胶结破坏形式对切向胶结强度与法向胶结强度的比值较为敏感.     

基于微观破损机理的胶结砂土三维本构模型研究

在岩土破损力学和临界状态土力学框架内,遵循宏微观土力学的研究思路,建立了胶结砂土三维本构模型。定义与重塑砂土屈服面几何相似但尺寸扩大的胶结砂土屈服面;采用经三维离散元验证的Lade-Duncan强度准则作为临界状态强度面;基于胶结材料微观力学理论并结合三维离散元模拟结果,获得具有微观力学机制的胶结破损规律;将胶结破损规律引入到重塑砂土的硬化规律和流动法则,得到胶结砂土的硬化规律和流动法则。将该本构模型应用于人工制备胶结砂土室内常规三轴压缩试验和等平均应力真三轴试验的模拟,初步验证了该模型的适用性。     

结构性砂土力学特性三维离散元分析

利用离散元法对结构性砂土的三轴试验进行了三维数值模拟并对其宏观特性进行了分析。首先将考虑胶结尺寸(宽度和厚度)的三维胶结接触模型导入离散元软件PFC3D中,对结构性砂土数值试样进行三轴试验数值模拟;然后对比分析离散元模拟与室内试验结果;最后从宏观力学角度对试验结果进行了分析。离散元模拟结果表明:结构性砂土与无胶结松散砂土表现不同,其在低围压时表现出应变软化和体积剪胀特征,并随胶结含量的增加或围压的减少而愈发显著,在高围压时则呈应变硬化和体积剪缩现象;低平均应力时,随胶结含量的增加,试样峰值内摩擦角、黏聚力以及内摩擦角均增加,其中黏聚力增加较为明显,随着平均应力的增加,峰值强度包线逐渐趋向于无胶结土。     

3种胶结岩土材料本构模型中破损规律的离散元验证

采用离散元法（DEM）对胶结岩土材料本构模型中的破损规律进行验证。首先,根据胶结颗粒间胶结特性,建立适用于胶结颗粒的微观理论;其次,基于该微观理论,给出3种经典胶结岩土材料本构模型破损参数的微观解释;最后,对数值试样进行等向压缩、等应力比压缩及双轴压缩试验的离散元数值模拟,并对比验证数值试样破损规律的模拟结果与本构模型的预测结果。结果表明：采用DEM对胶结岩土材料本构模型的破损规律进行验证是可行的。3种加载条件下,扰动状态模型假定的破损规律同模拟结果均较为一致;Nova模型在双轴压缩试验条件下的预测结果同模拟结果存在很大差异;而上负荷面剑桥模型假设的破损参数则不能很好描述上述加载条件下数值试样的破损规律。     

不排水循环荷载条件下胶结砂土宏微观力学性质离散元模拟研究

天然或人工胶结的存在能够提高砂土的抗液化能力,从宏微观尺度对其动力学性质进行研究具有重大意义。将已有的三维完整胶结接触模型引入到三维离散元程序中,对胶结砂土不排水循环三轴剪切试验进行三维离散元模拟,研究颗粒间胶结、循环应力比对离散元试样宏微观力学性质的影响。研究结果表明,胶结的存在能够抑制轴应变和孔压的发展,提高砂土的抗液化强度,循环应力比与液化振次之间具有指数函数关系,证实了本文离散元模拟能够反映胶结砂土的宏观动力学性质。在微观尺度上,当循环应力比较小时,胶结试样内部仅有极少量胶结发生破坏,力学配位数基本不变,外界输入功主要用于增加颗粒和胶结弹性能。对于特定胶结程度的试样,在初始液化发生之前,随循环应力比增加,试样内部胶结破坏更为剧烈,力学配位数下降速率更快,颗粒和胶结弹性能更快地趋向于0,颗粒摩擦耗能、弯转耗能、扭转耗能更快地达到最大值,而破坏胶结接触点、胶结接触点和无胶结接触点法方向的空间分布更快地趋向于各向同性性质。