颗粒流方法在砂土双轴试验中的应用.docxVIP

颗粒流方法在砂土双轴试验中的应用 1 pfc2c颗粒流土体性能模型 分散元法由cukal在70年代提出。作为一种分散元，二维流程序（pbc2d-fs2）是模拟固体力学变形和流动问题的计算方法。它模拟了圆形离散单元的运动和颗粒介质的相互作用。由平面内的平动和转动运动方程来确定每一时刻颗粒的位置和速度。作为研究颗粒介质特性的一种工具, 它采用有代表性的数百个至上万个颗粒单元, 通过数值模拟实验可以得到颗粒介质本构模型。PFC2D潜在的应用很多, 如岩石与土体的开挖问题;模拟颗粒间的相互作用问题、大变形问题和断裂问题等;槽、管、料斗和筒仓中的松散物体流动问题;动态冲击问题;以及介质基本特性研究, 例如屈服、流动、体变等。 本文充分利用颗粒流能模拟颗粒聚合体变形性能的特点, 构造相应的颗粒流砂土试样模型, 对砂土的双轴试验进行数值模拟。通过改变颗粒的性质以及颗粒单元的大小和分布, 得到一系列不同的试样宏观参数。对这些数值结果进行分析处理, 发现颗粒流程序能较好地模拟砂土的室内双轴试验, 尤其对土体变形的宏观和微观性能关系能给出很好的规律性。通过室内试样的颗粒流数值模拟, 调整颗粒的各种参数, 得到和试验结果相似的应力-应变关系曲线, 有一定的理论分析和实际应用价值。 2 试验条件及程序 为了进行颗粒流模拟试验, 先生成颗粒聚合试样, 通过移动墙来模拟加载过程和保持试样的围压 (见图1) , 最后可以得到试样二维应力-应变关系曲线。通过一系列试验可以得出试样的弹塑性关系曲线及破坏特性曲线, 如轴向偏应力-轴向应变试验曲线;加卸荷体变-轴向应变曲线;近峰值荷载时颗粒单元破坏图;应变软化与卸荷再加荷时的破坏图等。本文通过颗粒流模型模拟砂土试样在室内双轴的试验条件, 侧重于探讨颗粒单元的微观输入参数与模拟试样宏观应力-应变关系及强度之间的关系。 采用程序进行颗粒流模拟, 当只研究圆形颗粒的运动与相互作用问题时可以直接采用PFC2D来模拟, 不需增加单元的组合形式。为了简化, 本文将砂土颗粒理想化成圆形颗粒。构造长方形试样, 用模型的顶、底部墙模拟加载, 左、右侧墙模拟围压约束。给定顶、底部墙的移动速度模拟应变控制加载方式, 两侧墙的速度由程序自动控制, 使在整个试验过程中围压保持恒定。记录整个试验过程中墙体的位置, 颗粒的运动状况和颗粒之间的作用力, 通过后处理, 得到颗粒试样的宏观变形过程数据。 2.1 样品的生成 (1) 颗粒相互作用 土体试样通过一系列颗粒来模拟。为了和文献进行对比, 土体试样的尺寸为105 mm×40 mm。为消除采用均匀颗粒可能出现的各向异性, 颗粒单元半径R的分布采用从Rmin到Rmax均匀分布。设r为所有颗粒单元的平均半径, 由于颗粒半径均匀分布, 故r= (Rmin+Rmax) /2。综合考虑计算量和计算精度因素, 下文未特别说明时, 均取Rmin=0.3 mm, Rmax=0.6mm, 相应产生的颗粒数目为2000左右。 由于程序中采用的颗粒为圆球形颗粒, 每个颗粒应与其周围颗粒接触, 以免颗粒悬浮于试样内部。这样, 对均匀圆形颗粒来说, 其孔隙比的变化范围可以容易地计算出来为0.35～0.90, 对于不均匀颗粒的组合, 其孔隙比变化范围将有所差别。因此, 在计算中, 笔者取孔隙比的变化范围为0.35～0.90。 为模拟砂土颗粒间的相互作用, 颗粒在接触处有法向接触力Fnij、切向接触力Fsij, 摩擦力Ffij, 下标表示力由第i个颗粒单元通过接触作用于第j个颗粒单元上。分别通过法向刚度knij, 切向刚度ksij和摩擦系数μij和法向相对位移Unij和切向相对位移Usij按式 (1) 计算。 法向接触力沿两颗粒单元圆心的连线, 切向接触力与摩擦力则与之垂直。在本文的计算中法向刚度、切向刚度和摩擦系数取为knij=ksij=k和μij=μ, 颗粒的比重在试验中取为2.63。 (2) 应力-应变关系曲线 先定义墙体, 共4道, 其包围的矩形为105 mm×40 mm, 见图1 (a) 。注意, 墙体要长一些, 否则, 可能会出现颗粒溢出的现象, 即颗粒由于没有墙体的约束而自由运动。考虑到PFC2D在颗粒生成时, 把给定半径的颗粒往区域内填充, 如果没有已生成颗粒与之重叠, 则生成此颗粒, 否则, 改变颗粒的位置重试。为了保证颗粒的生成以及效率, 先生成小直径颗粒, 然后再把半径复原, 最后通过循环来消除试样内部非均匀应力。颗粒的数目由孔隙比e来控制, e可近似地用式 (2) 计算: 式中b为试样的宽度;h为试样的高度;N为颗粒单元总数。 图1给出了用颗粒流程序模拟室内双轴试验整个试样的颗粒分布情况。图示试样的颗粒半径变化范围 为0.3～0.6 mm, 共有2043个颗粒。图1 (b) 为生成小粒径颗粒后颗粒