颗粒形状对双轴试验的影响.docxVIP

颗粒形状对双轴试验的影响 0 基于pfcdd的模拟双轴试验 传统的三轴压缩试验是一种相对完善的方法，可以测量土壤或其他颗粒的弯曲强度。我们可以得到一些重要的力学性质，如抗切割强度、粘度、内摩擦角等。但土是由颗粒组成的集合体，颗粒的形状和大小、颗粒表面的粗糙程度、颗粒的级配、颗粒的含水率等这些细观因素的变化导致了颗粒体宏观力学性能的差异。常规的三轴压缩试验仅能得到土的宏观力学表现，而无法探明土在颗粒层面上的细观力学本质。离散单元法是一种广泛使用的分析离散颗粒体的数值方法，基于该方法的二维颗粒流程序PFC2D可方便的建立颗粒体试样进行模拟双轴压缩试验，研究不同颗粒特性对试样宏观力学性能的影响。还可采用三维颗粒流程序PFC3D更准确的模拟三轴压缩试验。采用离散元程序PFC2D进行颗粒体的模拟双轴试验还有一个重要的目的，因为在用离散元程序分析涉及颗粒体的具体问题时（如筒仓内的储料流动问题），首先需要得到具有一定宏观力学性能（如内摩擦角）的颗粒体，而PFC2D是通过指定颗粒的细观参数来得到颗粒体的宏观特性的，这两者之间并没有明确的对应，需要不断改变颗粒的细观参数，并进行模拟双轴试验来尝试得到具有期望宏观力学性能的颗粒体。因此采用离散元程序PFC2D进行模拟双轴试验非常重要，通常是采用PFC2D分析实际问题的第一步。 目前利用PFC2D模拟双轴压缩试验主要研究了颗粒刚度，颗粒摩擦系数和颗粒大小等细观因素的影响，以及颗粒接触的非线性本构关系。但上述研究中颗粒形状仅限于圆形颗粒，主要原因是PFC2D程序的基本颗粒形状只能是圆形，目前虽然有椭圆形等其它形状颗粒的离散元程序，但由于这些非圆形颗粒的引入，颗粒之间接触关系的判断变得非常复杂，使程序的计算速度大大降低，分析颗粒数量较多的试样需耗费大量的时间。 本文基于PFC2D程序的簇单元（clump），通过圆形颗粒组成三种异形颗粒，研究了颗粒形状对颗粒试样宏观特性的影响。簇单元是通过若干圆形的从属颗粒形成异形颗粒，簇单元内部从属颗粒之间的接触被忽略，与外界颗粒之间的接触才被考虑。因此PFC2D程序中簇单元形成的异形颗粒与真正的非圆形颗粒还是有较大区别的，但在颗粒数量相同的情况下，可较大的缩短计算时间，是研究非圆形颗粒试样的较好的方法。 1 数值模拟方法 1.1 接触刚度模型 本文生成的试样尺寸为85 mm×40 mm，采用了4种颗粒形状，其中圆形颗粒的最小颗粒半径为0.3mm，最大颗粒半径为0.6 mm，所生成的颗粒半径在0.3～0.6 mm随机分布。生成颗粒体的孔隙率为0.15，孔隙率的定义为 式中，ρ为孔隙率，Ap为颗粒体的面积，A为试样面积，总的颗粒数量约为4200个左右，所生成的四种试样的孔隙率和颗粒数量均基本相同。程序中提供了两种接触刚度模型，线性接触刚度模型和Hertz-Mindlin接触刚度模型，为了更好的考虑不同围压下颗粒接触刚度的变化，本文在下面的分析中均采用了Hertz-Mindlin接触刚度模型，该模型与线性接触刚度模型不同，是一种非线性接触模型，在程序中由颗粒的剪切模量G和泊松比ν两个参数决定。其颗粒间接触的法向割线刚度和切向切线刚度为 式中，Un为颗粒的重叠量，iFn为颗粒间的法向接触力，G和ν分别为两接触颗粒剪切模量以及泊松比的平均值，R～为两接触颗粒半径乘积与半径平均值的比值，颗粒接触刚度及其它细观参数的定义参考了文献的结果，颗粒密度为2650 kg/m3，孔隙率为0.15，剪切模量500 GPa，泊松比0.35，摩擦系数为0.5。 1.2 距离式异形颗粒的空间定位 本文通过簇单元由圆形从属颗粒建立异形颗粒，建立的4种试样的颗粒形状如图1所示。第一种为最基本的半径为r1的圆形颗粒；第二种为类长形颗粒，由3个半径为r2的圆形颗粒形成，3个圆形颗粒圆心之间的连线为一直线，且圆心之间距离为r2；第三种为类三角形颗粒，由3个半径为r3的圆形颗粒形成，3个圆形颗粒圆心的连线为边长为r3的正三角形；第四种为类正方形颗粒，由4个半径为r4的圆形颗粒形成，4个颗粒圆心的连线为边长为r4的正方形。为了保证三种异形颗粒外轮廓线所围成的面积均与半径为r1的圆形颗粒的面积相同，组成3种异形颗粒的从属颗粒半径如表1所示。 3种异形颗粒的空间定位除了需要其形心点的坐标之外，还需要一个表示其倾斜角度的方位角，如图2所示。类长形颗粒的方位角在0°～180°间变化，类三角形颗粒的方位角在0°～120°间变化，类四边形颗粒的方位角在0°～90°间变化。在生成3种异形颗粒时，其方位角是随机指定的，3种异形颗粒试样的初始方位角分布情况如图3所示。 1.3 宏观特性分析 双轴压缩试验的模拟环境如图4所示，颗粒体的双轴压缩是在上、下、左、右四面墙体所围成的区域内进行的。通过上部墙体向下移动和下部墙体向上移动对颗