粗粒土细观组构分析的影响因素岩土工程论文

来源：核心期刊咨询网时间：2018-08-29 15:4912

摘要：摘要：基于离散元方法，采用PFC2D模拟分析了土样在双轴等向压缩和双轴剪切两种不同加载方式下的细观组构及接触演化规律.通过采取不同的颗粒级配、不同的试样尺寸、不同的制样方法和不同的加载方式等数值模拟条件，从试样颗粒的接触力网络、平均接触数、组构

　　摘要：基于离散元方法，采用PFC2D模拟分析了土样在双轴等向压缩和双轴剪切两种不同加载方式下的细观组构及接触演化规律.通过采取不同的颗粒级配、不同的试样尺寸、不同的制样方法和不同的加载方式等数值模拟条件，从试样颗粒的接触力网络、平均接触数、组构偏量和应力应变关系变化等方面研究各种条件变化对粗粒土结构及细观组构演化规律的影响.研究表明：试样的颗粒级配和制样方法对粗粒土的细观组构和各向异性组构分析有一定的影响;试样的几何尺寸对组构各向异性分析的影响较小，但是试样尺寸大小对研究土样的应力—应变关系会有影响;加载方式是土样诱发各向异性和导致土样力学行为差异性的根本原因.

　　关键词：粗粒土;细观组构;颗粒接触;各向异性;离散元方法

　　土的力学性质的复杂性，包括应力路径相关性、各向异性和时间效应等，是与其细观组构特征密切相关的，因此对土的细观组构的形成和演化规律进行研究是认识土的力学行为机理的重要途径.在土的力学特性和细观组构联合研究方面，Oda等[1-3]通过研究发现散粒材料的初始各向异性主要表现在接触应力的各向异性分布、孔隙空间的优势方向和非圆球体颗粒的优势方向3个方面，基于图像分析法提出了使用二阶组构张量来描述诱发各向异性，并发现了颗粒材料的孔隙分布与外力加载的方向具有一致性[4].

　　Zhang等[5]研究了在循环剪切荷载条件下二维光弹颗粒结构中应力发展，利用光弹技术检测到接触网格与接触力，并且研究了组构张量、应力张量以及应力在法向和切向的分布规律.杨仲轩等[6-9]利用三轴试验研究砂土的各向异性并给出了定量方法，发现砂土力学特性受到材料各向异性和应力各向异性的耦合影响;周健等[10-14]利用了离散元程序对颗粒材料的力学变形特性进行了数值模拟研究，并且在颗粒流理论的基础上，还引入了不同颗粒的接触本构模型，探索适用于砂土和黏性土等不同颗粒材料的颗粒模型，这为深入分析土的细观力学机理提供了可行的途径和研究基础.Hu等[15]采用PFC2D模拟了循环荷载作用下砂土的诱发各向异性，研究发现试样的初始不等向固结方式、颗粒材料的法向与切向刚度比及循环荷载峰值等对砂土的诱发各向异性都会产生一定程度的影响.

　　研究表明：土的组构各向异性是颗粒材料最重要的组构特征，而采用离散元细观数值模拟的方法可以对土的各向异性组构和受力特性进行联合分析，是有效的方法.但是在现有的离散元细观组构研究中，采用的试样颗粒级配、试样尺寸和数值制样方法等都不统一，这对各种组构分析研究结果之间的相互参考性提出了质疑.为此，通过采用离散元数值模拟，对不同的颗粒级配、不同的试样尺寸、不同的制样方法以及不同的加载方式等各种因素下粒状土的细观参数变化及各向异性组构特征等进行分析，比较各种不同条件下土样的颗粒接触分布、平均配位数和组构偏量的变化规律，并且结合试样的应力应变特征分析了应用离散元方法研究粗粒土各向异性组构特征时需考虑的前期影响因素，可为室内试验和理论研究提供参考.

　　1数值模拟试验设计及细观组构参数

　　1.1颗粒材料参数及接触模型选择

　　为了研究制样方法和试样的物理及几何特征对试样的组构及初始接触的影响，采用颗粒流程序PFC2D对数值试样进行了双轴等向压缩和双轴剪切试验.

　　在数值试验中合理选择颗粒间的接触模型及模型参数对正确分析颗粒材料组构特征是至关重要的.颗粒流程序中提供了线性接触刚度模型、粘结模型以及滑动等模型，其中线性接触刚度模型能较好的反映颗粒材料在接触点处的挤压与摩擦作用，是一种可用于模拟散体颗粒单元间作用的简单且有效的模型[16]，为此，采用线性接触刚度模型来模拟砂土的力学性质.

　　由于颗粒材料的细观模拟参数很难由试验精确确定，为此，参考文献[15-18]中的取值方法，经过反复试算调整颗粒间法、切向刚度比.选取土样的细观模拟参数如表1所示

　　

　　1.2试样制备及数值模拟方案

　　为了研究试样的颗粒级配、制样方法、试样尺寸和加载方式等条件的影响，采用的具体数值模拟方案如表2所示.其中各个影响因素的控制方法如下：

　　1)颗粒级配.采用均匀分布和高斯分布(不均匀分布)方法生成了两种颗粒试样，分别表示为均匀试样(U)和不均匀试样(G)，其中前者控制颗粒粒径分布在较小的区间内(0.3～0.399mm);后者的颗粒粒径的分布范围较大(0.05～0.926mm).

　　2)几何尺寸.采用相同形状的2种不同尺寸，分别为80mm×30mm(B)，40mm×15mm(L).

　　3)制样方法.以孔隙比为标准，采用压实法(P)和振动法(Z)2种不同方式制备试样.前者是通过上边界墙体以一定速度向下运动压实试样，当边界应力达到一定值时停止运动，稳定一定时间后再继续压实，如此循环若干次;后者是控制4个边界以相同的速度向同一个方向运动，待其位移达到一定的值时，改变其速度的方向，达到原位置时继续反向以相同速度压实试样，如此反复若干次.

　　4)加载方式.采用双轴等向压缩(C)和双轴剪切(S)两种加荷方式.前者控制试样的上下边界和左右边界以相同的应变率向内移动压缩试样，即压缩的速率与运动方向上的试样尺寸成正比，实现试样的等向压缩;后者是通过控制上下边界以一定的速率相对移动进行压缩，而左右边界同时以一定的速率相对扩张，通过伺服控制以维持均应力稳定，实现试样的纯剪切加载

　　

　　1.3细观组构参数

　　颗粒流方法提供了一种从细观层面上观察和分析散体颗粒材料在荷载作用下的力学变形规律的研究手段.为了分析砂样的接触及其细观组构的变化规律，采用如下统计方法：

　　1)接触力网络又可称颗粒间的接触力链，可以直观地观察试样内部颗粒间接触力的分布以及接触力的优势方向.

　　2)接触力方向分布图可以反映颗粒间接触力在不同方向上概率分布情况，也可以反映在不同方向上平均接触力的大小.

　　2PFC2D数值模拟验证

　　为了验证离散元方法分析颗粒材料组构及接触特征的合理性，对Majmudar等[[21]进行的等向压缩光弹试验进行数值模拟.采用线性接触模型，并根据光弹试验模型确定PFC2D模型中试样几何尺寸、颗粒粒径及细观参数，具体如表3所示.

　　

　　通过统计等向压缩条件下，颗粒法向和切向接触力大小的概率分布如图1所示，横坐标为无量纲化的法向或切向接触力，即法向/切向接触力与平均法向/平均切向接触力的比值(F/Fmean);纵坐标表示法向或切向接触力的概率密度P(F/Fmean).从图1中可以看出：光弹试验[21]和数值试验中颗粒间的法向接触力大小的概率密度在平均法向接触力附近达到峰值，大于平均法向接触力时，法向接触力的概率密度分布逐渐减小，而切向接触力大小的概率密度分布一直随着切向接触力的增大而减小.

　　图2给出了试样在等向压缩条件下大于平均接触力的接触方向分布.对比图2(a，b)的结果可以发现：在等向压缩条件下，大于平均接触力的接触，其方向分布规律也具有明显的相似性.

　　

　　

　　综上分析可知：两类试验得到的接触力大小的概率密度分布规律是一致的，且接触方向的分布规律也是相似的.因此，利用PFC2D模拟砂样的细观组构及结构变化过程是可行且可靠的.

　　3数值试验结果分析与讨论为了分析各种条件对PFC2D方法研究砂土组构各向异性特征的影响，开展了表2中16个试样的数值加载试验，通过统计加载后试样内部的平均接触数和组构偏量的变化等细观参数，明确各种因素对模拟结果的影响.

　　3.1颗粒级配的影响

　　图3给出了2种颗粒级配试样在等向压缩和剪切加载条件下的应力—应变关系，从图3可以看出：等向压缩条件下的应力应变曲线基本相同，在剪切加载条件下的偏应力轴应变关系曲线的规律一致，在后期略有区别，表明不同颗粒级配对试样的应力应变关系的影响不大.

　　

　　图4为2种不同颗粒级配试样采用压实法制样后加载前的接触力网络，图4中线段粗细表示接触力的大小，把大于平均接触力接触定义为主要接触，通过对比图4(a，b)中在加载前的接触力网络可以发现不同颗粒级配试样的接触力网络略有不同，前者主要接触力较大，分布不均匀，而后者主要接触力相对较小，分布比较均匀.原因可能是不均匀分布粒径分布范围较大，起骨架作用的颗粒数较少，而均匀分布恰好相反，起到主要骨架作用的颗粒数较多.

　　

　　4结论

　　在验证离散元数值模拟研究粗粒土细观组构合理性的基础上，研究试样颗粒级配、试样尺寸、制样方式和加载方式等条件对进行砂土细观组构分析的影响，结果表明：试样颗粒级配对相同加载下土样的应力应变响应影响不大，但是对土样的细观组构分析有一些影响，主要表现在颗粒的接触数和组构偏量取值有所不同.因此，分析不同土类在特定应力路径下的细观组构特征时，应注意土样的来源、选取有代表性的颗粒级配;试样尺寸对土样的应力-应变响应有一定影响，但是对分析试样细观组构及其演化影响不明显.

　　因此，从砂土细观组构分析的角度，采用离散元分析可以采用较小的试样尺寸来提高数值模拟的效率;试样制备方法对试样的初始组构及加载后的组构演化有较显著的影响.因此当以土样的细观组构及各向异性演化为研究内容时，须考虑土体的生成特征及沉积环境的影响，采取相应的方法制备土样;加载方式对土样的组构特征和各向异性演化有重要影响，因此为了分析土样的组构特征和演化，须考虑土样的应力历史及加载路径的影响;采用等向压缩的方法可以消除土样由于制样不同形成的各向异性的影响，因此研究正常固结土和超固结土在不同加载路径下的组构变化时，可以通过先施加等向压力再施加固结压力的方式对试样进行预压.

转载请注明来自：http://www.qikan2017.com/lunwen/lig/12461.html