Processing math: 0%

基于离散元法的砂石混合体直剪试验结果分析

苏永华, 王栋. 基于离散元法的砂石混合体直剪试验结果分析[J]. 水文地质工程地质, 2021, 48(6): 97-104. doi: 10.16030/j.cnki.issn.1000-3665.202012029
引用本文: 苏永华, 王栋. 基于离散元法的砂石混合体直剪试验结果分析[J]. 水文地质工程地质, 2021, 48(6): 97-104. doi: 10.16030/j.cnki.issn.1000-3665.202012029
SU Yonghua, WANG Dong. An analysis of direct shear test results of sand - gravel mixture based on the discrete element method[J]. Hydrogeology & Engineering Geology, 2021, 48(6): 97-104. doi: 10.16030/j.cnki.issn.1000-3665.202012029
Citation: SU Yonghua, WANG Dong. An analysis of direct shear test results of sand - gravel mixture based on the discrete element method[J]. Hydrogeology & Engineering Geology, 2021, 48(6): 97-104. doi: 10.16030/j.cnki.issn.1000-3665.202012029

基于离散元法的砂石混合体直剪试验结果分析

  • 基金项目: 国家自然科学基金项目(51878266)
详细信息
    作者简介: 苏永华(1965-),男,教授,博士,主要从事地下结构和岩土工程研究。E-mail:yong_su1965@126.com
    通讯作者: 王栋(1995-),男,硕士研究生,主要从事地表塌陷方面研究。E-mail:975662835@qq.com
  • 中图分类号: TU43

An analysis of direct shear test results of sand - gravel mixture based on the discrete element method

More Information
  • 砂石混合体由力学性质以及结构相差极大的材料组成,其组成的重塑地层易发生塌陷等问题,因此对砂石混合体力学特性的研究具有重要的工程意义。砾石形状是砂砾石力学特性研究的重要属性参数,但采用规则图形对砾石进行描述不能反映出其真实的力学性质,采用数字图像处理技术构建的砾石数据库能反映砾石真实形状并可对特定形状参数进行具体分析。由于砂石混合体的粒径分布较广,采用特征粒径等无法描述整体粒度分布,故本文结合分形理论构建砂石混合体的二重分形结构模型,通过粒度分维值反演出完整的级配分布曲线。考虑到砂石混合体离散型的特点,采用离散元软件进行直剪试验数值模拟并对细观结构进行分析,研究结果表明,砂石混合体一般具有2个粒度分维值:砂粒度分维值和砾石粒度分维值,砂、砾石粒度分维值越接近,抗剪强度和内摩擦角越大;当两者相等时,砂石混合体具有一重分维,此时均一性最好,抗剪强度和内摩擦角最大;轴向系数是形容砾石形状的一个重要参数,随着轴向系数的增加,砾石显示出明显的条状性,在直剪试验中抗转动能力增强、周围接触数量增加,导致抗剪强度和内摩擦角不断增加。

  • 对于地表高程较大甚至悬殊区域,在建、构筑物施工前都要进行一定的或较大、甚至大规模的回填处理,这种状况下建、构筑物下的地层就称为重塑地层,砂石地层是常见的重塑地层之一,在基坑、路基等工程中广泛存在。砂石混合体由力学性质以及结构相差极大的材料组成,是典型的土石混合体。砂石混合体构成的地层由于颗粒粒度变化较大、分布不均匀、块石形态各异等原因,在实际工程中易出现沉陷或塌陷等问题[1-2],因此,砂石混合体力学性质研究是亟需进行的研究课题之一。

    室内试验是研究土石混合体力学性质的基本方法,考虑不同影响因素对宏观力学性能的影响是研究的重点。但室内试验往往受尺寸效应和试验条件等因素的影响,存在可重复性差、耗时较久、不易对颗粒形状进行分析的限制。随着计算机的发展,颗粒流模拟(PFC)成为研究土石混合体力学性质的一种方法,它能够模拟微观层面的结构,充分考虑砂石混合体离散性的特点,从不同角度探讨力学性质,如:颗粒实际形状、颗粒转动、力链演化和孔隙变化等[3]

    砾石形状是砂砾石力学特性研究的重要属性参数,目前许多学者采用规则图像对砾石进行描述,如圆形、正多边形,没有考虑砾石的真实形状,很难反映出砾石真实的力学性质。许多学者利用数字图像处理技术获取真实的细观结构模型,如徐文杰[4]、廖秋林等[5]通过实际拍摄的数码图像生成土石混合体结构模型。还有学者通过识别自然块石的轮廓构建块石数据库,实现块石二维[6]、三维[7]随机模型建立,数字图像处理技术可以在不忽略砾石真实形状的前提下,对形状参数进行细致化分析,探究其影响作用。

    粒度分布影响着土石混合体的变形、强度等物理力学特性。在现有的岩土工程散体材料的粒度分布描述中,大多是采用特征粒径、不均匀系数和曲率系数等参数来描述粒度分布,这些指标只能表示某些关键点的信息,缺乏对粒度整体分布的描述。针对上述情况,杜修力等[8]通过对多组试样的研究,发现土石混合体的粒度分布均满足分形结构,可以利用粒度分维值定量描述。舒志乐等[9]建立了土石混合体的二重分形结构模型,对土石混合体的粒度分维值与抗剪强度关系进行分析。

    本文通过数字图像处理技术获取砾石边界信息,建立砾石形状数据库。结合分形理论建立了砂石混合体二重分形结构模型。通过颗粒流软件对不同粒度分布、颗粒形状的砂石混合体进行直剪试验,并对结果进行分析,为进一步研究砂石混合体重塑地层塌陷机理奠定基础。

    通过数字图像处理技术获取砾石的边界坐标和形状参数,并导入PFC中生成颗粒簇模板以供模型构建时调用,详细步骤如下:

    (1)图像获取及预处理:通过高精度数码相机获取大量砾石图像,但是砾石因光线等环境因素图像会出现大量噪音,需要利用图像处理软件对图像进行消噪,并使每张图片仅显示1个砾石(图1a)。

    图 1.  数字图像处理技术获取边界信息流程示意图
    Figure 1.  Flow chart showing the digital image processing technology to obtain boundary information

    (2)边缘提取:预处理后的图像需进行二值化处理,可通过设定合适阈值,使图像变得简单,凸显感兴趣目标轮廓,处理后图像呈现明显的黑白效果(图1b)。二值化图像可通过边缘检测获取砾石边界,本文利用MATLAB软件中内嵌的Canny算法提取砾石边界(图1c)。

    (3)边界拟合:通过上述方法提取的边界点数目过多,会对后续工作产生较大负担,需要对边界进行拟合,拟合方法如下[6]:选定初始点,然后顺时针获取所有边界点,依次计算边界点和初始点的距离,如果该距离满足式(1)则保留该点为拟合后的边界点,并作为下一轮的初始点,依此循环直至回到初始点。

    d (1)

    式中:d——拟合后两个相邻边界点的距离;

    D——与砾石颗粒面积相等的等效圆直径;

    \alpha ——控制边界拟合精度的参数。

    通过调整\alpha 可以控制砾石边界的拟合精度。\alpha 越大所需边界点越少,但精度越低。本文选取\alpha 为0.1,此时砾石的形状和原图像基本一致(图1d),坐标点的数量也满足后续的数值模拟需求。

    (4)图像标准化处理:为了方便后续对砾石大小、旋转角度等参数的控制,需进行标准化处理,如粒径统一化、形心归零化、长轴水平化。

    (5)PFC颗粒簇模板生成:获取的砾石边界信息不能直接被PFC识别,需要将其转换为矢量数据格式文件(dxf)。本文利用Python第三方库dxfwrite完成批量转换,导入PFC并利用内置的颗粒填充算法,就可生成砾石的颗粒簇模板。

    颗粒形状是研究砂石混合体宏观力学行为的一个重要参数。根据相关研究[10],颗粒可以通过轴向系数和圆度进行描述,其中轴向系数S为:

    S = \frac{{L_{{\rm{A}},\max }}} {{L_{{\rm{A}},\min }}} (2)

    式中:LA,maxLA,min——砾石颗粒等效椭圆的长轴和短轴。

    砾石颗粒的圆度公式定义为R

    R = \frac{L^2}{4{\text{π}} {\rm{A}}} (3)

    式中:LA——颗粒的周长和面积。

    通过砾石形状数据库可导出形状参数,并利用数据分析软件得出圆度和轴向系数的频率分布直方图(图2),圆度分布在1.0~1.35之间,均值为1.15,圆度指标可有效地表征颗粒整体形状接近圆的程度。圆度越大,其形状越偏离圆形。轴向系数分布在1.0~1.8之间,均值为1.28,轴向系数越大,砾石的长轴与短轴比值越大,条状性越明显。

    图 2.  圆度频率(a)和轴向系数频率(b)分布直方图
    Figure 2.  Distribution histogram of (a) roundness frequency and (b) axial coefficient frequency

    Tyler等[11]提出了基于分形理论的分形结构模型:

    \frac{{M\left( {r \lt R} \right)}}{{{M_{\rm{T}}}}} = {\left( {\frac{R}{{{R_{\rm{L}}}}}} \right)^{3 - D}} (4)

    式中:r——颗粒粒径;

    R——某一颗粒粒径;

    MT——总质量;

    RL——最大颗粒尺寸;

    M(r<R)——粒径小于R颗粒累积百分含量;

    D——粒度分维值。

    根据式(4)得出,如果在双对数坐标下lg[M(r<R)- MT]-lg(R)存在直线段,就表明粒度分布满足分形结构,其中直线斜率为3-D

    砂石混合体的粒度分布具有良好的分形结构,但是由于其复杂的颗粒组成,通常具有多重分维。本文结合砂石混合体的组成以及土石阈值构建二重分形结构模型:

    \frac{{M\left( {r \lt R} \right)}}{{{M_{\rm{T}}}}} = \left\{ \begin{aligned} &{{{\left( {\frac{R}{{{R_{\rm{L}}}}}} \right)}^{3 - {D_{\rm{s}}}}},\;R \lt {R_{{\rm{S/RT}}}}}\\ &{{{\left( {\frac{R}{{{R_{\rm{L}}}}}} \right)}^{3 - {D_{\rm{r}}}}},\;R \geqslant {R_{{\rm{S/RT}}}}} \end{aligned} \right. (5)

    式中:Ds——砂粒度分维值;

    Dr——砾石粒度分维值;

    RS/RT——土石阈值。

    土石阈值,即“土”和“石”的界限粒径,本文结合以前学者的研究结果[12-13],选取5 mm作为砂和砾石的界限粒径。令式(5)中R= RS/RT即可确定砾石的最大粒径,公式如下[14]

    {R_{\rm{L}}} = {R_{{\rm{S/RT}}}}{\left( {\frac{{100}}{{100 - P}}} \right)^{\frac{1}{{3 - {D_{\rm{r}}}}}}} (6)

    式中:P——含石量。

    徐文杰等[15]指出含石量在25%~75%之间时,岩土体的强度取决于其中的“土体”与“块石”,其力学强度是两者相互作用的共同反映。含石量在这范围之外,整体的强度会主要取决于“土体”或者“块石”。为了分别研究砂和砾石级配变化对土石混合体强度的影响,需考虑两者共同作用,因此本文选取含石量为50%。

    根据式(6)得出最大粒径、含石量和砾石粒度分维值的关系(图3),可以看出:随着含石量和砾石粒度分维值的增加,最大粒径变大;最大粒径保持不变时,砾石粒度分维值的增加,会导致含石量减少,两者关系曲线呈近似抛物线形状。

    图 3.  最大粒径、含石量和砾石粒度分维值的关系
    Figure 3.  Relationship among the maximum grain size, rock content and fractal dimension of gravel size

    在确定含石量和土石阈值后即可根据不同粒度分维值确定砂石混合体的完整级配分布。首先需要确定每组粒径区间,取区间的粒径上限为该组粒径,并且粒组按粒径由大到小的顺序,即第1个粒组的粒径为砾石最大粒径,则第i个粒组含量占总试样的百分数:

    \Delta R\left( {{i}} \right){\rm{ = 100}}\left[ {{{\left( {\frac{{R\left( 1 \right)}}{{R\left( i \right)}}} \right)}^{ - \left( {3 - D} \right)}} - {{\left( {\frac{{R\left( 1 \right)}}{{R\left( {i + 1} \right)}}} \right)}^{ - \left( {3 - D} \right)}}} \right] (7)

    式中:R(i)——砾石第i组粒径上限。

    将粒度分维值DrDs代入式(7)就可确定不同粒度分维值下砂石混合体的试样级配。

    (1)剪切盒建立:砂和砾石颗粒单元均需要在剪切盒内生成。剪切盒尺寸和室内大型直剪试验截面一致,为300 mm×150 mm。模拟过程中剪切盒被认为是刚性体,剪切盒的刚度需大于砂和砾石的刚度,保证颗粒不飞出,并且避免模型内部产生较大初始接触力。通过试算,本文剪切盒刚度选取为1×1010 N/m 。

    (2)试样生成:砂采用ball模拟,砾石采用前述构建的颗粒簇模板生成,接触模型采取线性接触模型。颗粒级配利用确定的二重分形结构模型,计算每个粒组需要生成颗粒的数量,并采用粒径扩大法生成特定级配的试样。需要注意的是:如果颗粒大小按照1∶1的比例进行模拟,会对计算机提出较高要求,造成效率的降低。本文通过试算,设置颗粒放大系数为1.6,并忽略0.5 mm以下颗粒,减少生成颗粒的数量,从而平衡PFC模拟的精度和效率。

    (3)施加竖向荷载:颗粒流软件中墙体无法直接施加力的作用,需要利用墙体的伺服机制,即在每一时步调整墙体的速度来控制墙体的受力大小。通过施加竖向荷载形成均匀试样,待稳定后,便可进行试样剪切。本文选取的主要正应力为100,200,300 kPa。

    (4)试样剪切:剪切过程中维持上剪切盒不动,下剪切盒以一定速度向右移动,当剪切位移达到30 mm时停止剪切。剪切过程中监测剪应力和剪切位移的变化,记录数据。

    PFC中颗粒的细观参数不能直接通过宏观参数获得,需要通过室内试验获取应力-应变曲线,并通过不断改变细观参数,使得数值模拟的结果与试验结果相近,完成细观参数的标定。本文采用Geotest S2450型摩擦剪切仪进行直剪试验,砾石和粗砂作为试验材料。

    由于砂石混合体参数较多,因此本文首先通过对单一砂进行直剪试验,获取砂的细观参数。并以此为基准进行砂石混合体直剪试验,完成其余参数的标定。经过大量参数试算,最终模拟曲线和实际试验曲线良好(图4),具体参数见表1

    图 4.  砂(a)和砂石混合体(b)剪应力-剪切位移模拟和实际试验结果对比曲线
    Figure 4.  Shear stress-shear displacement comparison curve between the simulation and actual test results of (a) sand and (b) sand - gravel mixture
    表 1.  基于离散元模拟的直剪试验主要计算参数
    Table 1.  Main computational parameters in DEM simulation of direct shear
    计算参数砾石
    法向刚度kn/(5×1085×109
    切向刚度ks/(2×1082×109
    摩擦系数0.50.9
    粒径/mm0.5~5.0>5.0
     | Show Table
    DownLoad: CSV

    砂石混合体的抗剪强度明显高于单一砂时的抗剪强度。随着正应力从100 kPa增加至300 kPa,颗粒间受到挤压排列更加紧密,剪应力逐渐增加。不同正应力下剪应力-剪切位移曲线呈明显的非线性特征,可分为3个阶段[16]:线弹性阶段、初始屈服阶段以及应变硬化阶段,剪切过程中没有出现明显的峰值。

    砂石混合体的抗剪强度分别受砂粒度分维值和砾石粒度分维值影响。为单独研究砾石粒度分维值的影响,采用控制变量法进行研究,首先保持砂的粒度分维值为2.55不变,砾石的粒度分维值分别取:2.45,2.50,2.55,2.60,2.65。

    砂石混合体的粒度分布曲线按照前文构建的二重分形结构模型计算,图5为砂石混合体级配示意图,为了方便观察仅列出3条典型曲线(下同),可以看出随着Dr增大,最大粒径也增大,在Dr=2.65时,最大粒径达到32.23 mm。

    图 5.  不同砾石粒度分维值的级配曲线
    Figure 5.  Grading curves of fractal dimension of different gravel sizes

    图6是相应的剪应力-剪切位移曲线,采用了逐步平均法进行平滑,可以发现3条曲线都近似呈应变硬化型,在Dr=2.55时,曲线末尾其剪应力上升趋势还十分明显。抗剪强度与砾石的粒度分维值不是简单的线性变化,Dr为2.55时抗剪强度最大,其值为220.05 kPa,而在Dr=2.45和2.65时,其值分别为196.07 kPa和195.48 kPa,均明显小于Dr=2.55的抗剪强度。

    图 6.  300 kPa正应力下不同Dr剪应力-剪切位移曲线
    Figure 6.  Shear stress-shear displacement curves of different Dr under the normal stress of 300 kPa

    依据摩尔库伦准则τ=σtanφ+c,可以进行抗剪强度参数的计算。砂石混合体由砂和砾石组成,两者的黏聚力非常小,所以本文不考虑黏聚力,取值为0 MPa,此时抗剪强度主要由内摩擦角产生(表2)。

    表 2.  不同砾石粒度分维值在直剪试验中获取的相关参数
    Table 2.  Correlation parameters of fractal dimension of different gravel sizes obtained in the direct shear test
    砾石粒度分维值100 kPa200 kPa300 kPa内摩擦角/(°)
    2.4568.57144.84196.0734.06
    2.5065.22147.33209.9735.26
    2.5569.56150.10220.0536.34
    2.6065.22123.23210.7833.99
    2.6564.84112.112195.4832.02
     | Show Table
    DownLoad: CSV

    采用同样方法研究砂粒度分维值的影响,保持砾石的粒度分维值为2.55不变,改变砂的粒度分维值,最终得出不同正应力下抗剪强度及内摩擦角的变化(表3)。

    表 3.  不同砂粒度分维值在直剪试验中获取的相关参数
    Table 3.  Relevant parameters obtained in the direct shear test of the fractal dimension of different sand particle sizes
    砂粒度分维值100 kPa200 kPa300 kPa内摩擦角/(°)
    2.4565.22143.75175.3132.11
    2.5066.09150.45202.7534.86
    2.5569.57150.10220.0536.34
    2.6068.70135.24203.5334.15
    2.6564.44134.01170.9931.13
     | Show Table
    DownLoad: CSV

    为了方便观察,以砂、砾石粒度分维差值作为指标进行分析(图7),明显看出内摩擦角与DrDs不是简单的线性变化,而是呈现先增大后减小的趋势。随着砂、砾石粒度分维值的接近,内摩擦角逐渐增加,当DrDs=0即砂石混合体具有单一分维时,内摩擦角达到最高值36.34°。

    图 7.  粒度分维差值与内摩擦角关系
    Figure 7.  Relationship between the difference in particle size fractal dimension and the angle of internal friction

    为了探究影响原因,通过PFC绘制出剪切结束时颗粒位移图(图8),颗粒的位移主要由两部分造成:第一部分是由下剪切盒的移动造成的,剪切过程中,上剪切盒保持不动,下剪切盒向右移动了30 mm;第二部分是由于剪切过程中颗粒的相互错动造成的位移。剪切盒的右下角颗粒位移最大,有的超过了30 mm,这主要是由于剪切过程中颗粒的移动并不能与剪切盒同步,在短时间里颗粒右侧形成一段空隙,使其变得较为松散,颗粒接触减少,易产生错动位移。由于上剪切盒在剪切过程中没有移动,位移主要由颗粒之间相互作用影响,因此具体分析上剪切盒颗粒位移。分析10~15 mm剪切带附近的颗粒位移可以看出,随着砂、砾石粒度分维值的接近,10~15 mm位移的颗粒范围在不断扩大,未受影响的颗粒范围在减少。这说明随着砂、砾石粒度分维值的接近,可观察到颗粒排列的越紧密,因剪切盒移动产生的相互错动越明显,颗粒位移增加,导致抗剪强度增加,内摩擦角也增大。

    图 8.  不同粒度分维差值的位移变化图
    Figure 8.  Displacement change diagram of difference of different fractal dimensions

    综上,当砂、砾石粒度分维值越接近,其抗剪强度越大,内摩擦角也相应越大。当两者相等时,即砂石混合体只具有单一分维,此时试样的均一性达到最好,具有最大的抗剪强度和内摩擦角。

    由于本文在边界平滑和颗粒填充算法中对砾石细部边界信息有所缺失,因此只对轴向系数进行具体分析。从砾石数据库选取颗粒圆度大致相同,轴向系数分别为1.0,1.4,1.8的砾石进行模拟(图9)。为了减少其他方面影响,模拟时选取单一颗粒形状。由图9可以看出,轴向系数越大,颗粒条状性越明显,反之,越接近于圆形。

    图 9.  不同轴向系数的砾石
    Figure 9.  Gravel with different axial coefficients

    对不同正应力得到的剪应力峰值进行线性拟合(图10),可以发现拟合相关系数R2都在0.99以上(表4),表明通过摩尔库伦强度准则得到的内摩擦角比较可靠。随着轴向系数的增加,抗剪强度和内摩擦角逐渐增加,S=1.8时内摩擦角达到最高值35.92°。产生这种情况主要是因为,砾石轴向系数越大,颗粒条状性越明显,从而有更多机会接触不同颗粒,颗粒的抗转动能力也得到了提高[17]

    图 10.  不同轴向系数下剪应力峰值-正应力曲线
    Figure 10.  Shear stress peak-normal stress curve under different axial coefficients
    表 4.  不同轴向系数在直剪试验中获取的强度参数
    Table 4.  Strength parameters obtained in direct shear test with different axial coefficients
    主要参数轴向系数
    S=1.0S =1.4S =1.8
    100 kPa65.3168.9183.95
    200 kPa127.39128.25143.04
    300 kPa174.61184.58211.15
    400 kPa235.34251.42292.50
    内摩擦角/(°)30.7632.1435.92
    相关系数R20.99190.99680.9925
     | Show Table
    DownLoad: CSV

    由图11可以发现,颗粒在剪切面附近旋转角度偏高,并且越靠近左右边缘旋转角度越大。由于颗粒移动的延迟性以及不均匀导致剪切盒右下方砾石周围的接触数量也明显降低。

    图 11.  砾石颗粒旋转角度及周围接触示意图
    Figure 11.  Schematic diagram of rotation angle and surrounding contact of gravel particles

    为具体了解相关变化,通过Fish语言获取300 kPa正应力下砾石和周围颗粒的接触数量以及旋转角度(表5)。由表5可以看出,随着轴向系数从1.0至1.8变化,接触数量从1 581增加至1 742,砾石的平均旋转角度从9.43°减少至6.42°。这说明随着轴向系数增加,砾石颗粒的条状性越强,和其它颗粒的接触可能性越多,也导致其在旋转时遇到的阻力越多,抗转动能力较强,在剪切过程中旋转的角度总体上偏低,从而导致抗剪强度提高。

    表 5.  不同轴向系数下砾石接触数量和旋转角度
    Table 5.  Gravel contact number and rotation angle under different axial coefficients
    轴向系数砾石接触数量砾石平均旋转角度/(°)
    1.015819.43
    1.416568.70
    1.817426.42
     | Show Table
    DownLoad: CSV

    (1)利用数字图像处理技术能够自动化完成砾石真实图像至颗粒簇模板的转换,构建砾石的形状数据库,并可对轴向系数、圆度等形状参数进行分析。

    (2)基于分形理论建立了砂石混合体的二重分形结构模型,得出完整的粒度曲线。土石阈值一定时,砂石混合体的最大粒径随含石量、砾石粒度分维值的增加而变大。最大粒径保持不变时,砾石粒度分维值的增加,会导致含石量减少,两者关系曲线呈近似抛物线形状。

    (3)砂、砾石粒度分维值越接近,抗剪强度越大,内摩擦角也相应越大,当两者相等时,即砂石混合体只具有单一分维时,此时试样的均一性达到最好,具有最大的抗剪强度和内摩擦角。

    (4)轴向系数S是形容砾石形状的一个重要参数,其分布在1.0~1.8之间。随着轴向系数S的增加,砾石条状性越强,在直剪试验中抗转动能力越强、与周围接触可能性越大,导致抗剪强度和内摩擦角不断增加。

  • 图 1  数字图像处理技术获取边界信息流程示意图

    Figure 1. 

    图 2  圆度频率(a)和轴向系数频率(b)分布直方图

    Figure 2. 

    图 3  最大粒径、含石量和砾石粒度分维值的关系

    Figure 3. 

    图 4  砂(a)和砂石混合体(b)剪应力-剪切位移模拟和实际试验结果对比曲线

    Figure 4. 

    图 5  不同砾石粒度分维值的级配曲线

    Figure 5. 

    图 6  300 kPa正应力下不同Dr剪应力-剪切位移曲线

    Figure 6. 

    图 7  粒度分维差值与内摩擦角关系

    Figure 7. 

    图 8  不同粒度分维差值的位移变化图

    Figure 8. 

    图 9  不同轴向系数的砾石

    Figure 9. 

    图 10  不同轴向系数下剪应力峰值-正应力曲线

    Figure 10. 

    图 11  砾石颗粒旋转角度及周围接触示意图

    Figure 11. 

    表 1  基于离散元模拟的直剪试验主要计算参数

    Table 1.  Main computational parameters in DEM simulation of direct shear

    计算参数砾石
    法向刚度kn/(5×1085×109
    切向刚度ks/(2×1082×109
    摩擦系数0.50.9
    粒径/mm0.5~5.0>5.0
    下载: 导出CSV

    表 2  不同砾石粒度分维值在直剪试验中获取的相关参数

    Table 2.  Correlation parameters of fractal dimension of different gravel sizes obtained in the direct shear test

    砾石粒度分维值100 kPa200 kPa300 kPa内摩擦角/(°)
    2.4568.57144.84196.0734.06
    2.5065.22147.33209.9735.26
    2.5569.56150.10220.0536.34
    2.6065.22123.23210.7833.99
    2.6564.84112.112195.4832.02
    下载: 导出CSV

    表 3  不同砂粒度分维值在直剪试验中获取的相关参数

    Table 3.  Relevant parameters obtained in the direct shear test of the fractal dimension of different sand particle sizes

    砂粒度分维值100 kPa200 kPa300 kPa内摩擦角/(°)
    2.4565.22143.75175.3132.11
    2.5066.09150.45202.7534.86
    2.5569.57150.10220.0536.34
    2.6068.70135.24203.5334.15
    2.6564.44134.01170.9931.13
    下载: 导出CSV

    表 4  不同轴向系数在直剪试验中获取的强度参数

    Table 4.  Strength parameters obtained in direct shear test with different axial coefficients

    主要参数轴向系数
    S=1.0S =1.4S =1.8
    100 kPa65.3168.9183.95
    200 kPa127.39128.25143.04
    300 kPa174.61184.58211.15
    400 kPa235.34251.42292.50
    内摩擦角/(°)30.7632.1435.92
    相关系数R20.99190.99680.9925
    下载: 导出CSV

    表 5  不同轴向系数下砾石接触数量和旋转角度

    Table 5.  Gravel contact number and rotation angle under different axial coefficients

    轴向系数砾石接触数量砾石平均旋转角度/(°)
    1.015819.43
    1.416568.70
    1.817426.42
    下载: 导出CSV
  • [1]

    柴贺军, 陈谦应, 孔祥臣, 等. 土石混填路基修筑技术研究综述[J]. 岩土力学,2004,25(6):1005 − 1010. [CHAI Hejun, CHEN Qianying, KONG Xiangchen, et al. Overview of soil-stone high embankment construction study[J]. Rock and Soil Mechanics,2004,25(6):1005 − 1010. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-7598.2004.06.037

    [2]

    廖秋林, 李晓, 郝钊, 等. 土石混合体的研究现状及研究展望[J]. 工程地质学报,2006,14(6):800 − 807. [LIAO Qiulin, LI Xiao, HAO Zhao, et al. Current status and future trends of studies on rock and soil aggregates (RSA)[J]. Journal of Engineering Geology,2006,14(6):800 − 807. (in Chinese with English abstract) doi: 10.3969/j.issn.1004-9665.2006.06.014

    [3]

    马昊, 黄达, 肖衡林, 等. 江北机场高填方夯后碎块石土剪切力学性质研究[J]. 水文地质工程地质,2019,46(3):88 − 94. [MA Hao, HUANG Da, XIAO Henglin, et al. A study of the shear mechanical properties of high-filled gravel-block soil after dynamic compaction near the Jiangbei airport[J]. Hydrogeology & Engineering Geology,2019,46(3):88 − 94. (in Chinese with English abstract)

    [4]

    徐文杰, 胡瑞林, 岳中琦, 等. 土石混合体细观结构及力学特性数值模拟研究[J]. 岩石力学与工程学报,2007,26(2):300 − 311. [XU Wenjie, HU Ruilin, YUE Zhongqi, et al. Mesostructural character and numerical simulation of mechanical properties of soil-rock mixtures[J]. Chinese Journal of Rock Mechanics and Engineering,2007,26(2):300 − 311. (in Chinese with English abstract) doi: 10.3321/j.issn:1000-6915.2007.02.011

    [5]

    廖秋林, 李晓, 朱万成, 等. 基于数码图像土石混合体结构建模及其力学结构效应的数值分析[J]. 岩石力学与工程学报,2010,29(1):155 − 162. [LIAO Qiulin, LI Xiao, ZHU Wancheng, et al. Structure model construction of rock and soil aggregate based on digital image technology and its numerical simulation on mechanical structure effects[J]. Chinese Journal of Rock Mechanics and Engineering,2010,29(1):155 − 162. (in Chinese with English abstract)

    [6]

    赵鑫曜, 陈建功, 张海权, 等. 基于块石形状数据库的土石混合体模型随机生成方法[J]. 岩土力学,2018,39(12):4691 − 4697. [ZHAO Xinyao, CHEN Jiangong, ZHANG Haiquan, et al. Random generation of soil-rock mixture models by rock shape database using digital imaging technology[J]. Rock and Soil Mechanics,2018,39(12):4691 − 4697. (in Chinese with English abstract)

    [7]

    袁斌, 霍宇翔, 巨能攀, 等. 颗粒棱度指标的改进及其对剪切特性的影响[J]. 水文地质工程地质,2020,47(4):167 − 173. [YUAN Bin, HUO Yuxiang, JU Nengpan, et al. Improvement of grain edge index and its effect on shear characteristics[J]. Hydrogeology & Engineering Geology,2020,47(4):167 − 173. (in Chinese with English abstract)

    [8]

    杜修力, 张佩, 金浏, 等. 基于分形理论的北京地区砂砾石地层细观建模[J]. 岩石力学与工程学报,2017,36(2):437 − 445. [DU Xiuli, ZHANG Pei, JIN Liu, et al. A mesoscopic model of sand-gravel stratum in Beijing based on fractal theory[J]. Chinese Journal of Rock Mechanics and Engineering,2017,36(2):437 − 445. (in Chinese with English abstract)

    [9]

    舒志乐, 刘新荣, 刘保县, 等. 土石混合体粒度分形特性及其与含石量和强度的关系[J]. 中南大学学报(自然科学版),2010,41(3):1096 − 1101. [SHU Zhile, LIU Xinrong, LIU Baoxian, et al. Granule fractal properties of earth-rock aggregate and relationship between its gravel content and strength[J]. Journal of Central South University (Science and Technology),2010,41(3):1096 − 1101. (in Chinese with English abstract)

    [10]

    汪海年, 郝培文. 粗集料二维形状特征的图像描述[J]. 建筑材料学报,2009,12(6):747 − 751. [WANG Hainian, HAO Peiwen. Digital description of two-dimensional shape characteristics of coarse aggregate[J]. Journal of Building Materials,2009,12(6):747 − 751. (in Chinese with English abstract) doi: 10.3969/j.issn.1007-9629.2009.06.025

    [11]

    TYLER S W, WHEATCRAFT S W. Fractal scaling of soil particle-size distributions: analysis and limitations[J]. Soil Science Society of America Journal,1992,56(2):362 − 369. doi: 10.2136/sssaj1992.03615995005600020005x

    [12]

    王江营, 曹文贵, 张超, 等. 基于正交设计的复杂环境下土石混填体大型直剪试验研究[J]. 岩土工程学报,2013,35(10):1849 − 1856. [WANG Jiangying, CAO Wengui, ZHANG Chao, et al. Large-scale direct shear tests on soil-rock aggregate mixture under complicated environment based on orthogonal design[J]. Chinese Journal of Geotechnical Engineering,2013,35(10):1849 − 1856. (in Chinese with English abstract)

    [13]

    董云. 土石混合料强度特性的试验研究[J]. 岩土力学,2007,28(6):1269 − 1274. [DONG Yun. Experimental study on intensity character of rock-soil aggregate mixture[J]. Rock and Soil Mechanics,2007,28(6):1269 − 1274. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-7598.2007.06.040

    [14]

    徐文杰, 胡瑞林, 岳中崎. 土–石混合体随机细观结构生成系统的研发及其细观结构力学数值试验研究[J]. 岩石力学与工程学报,2009,28(8):1652 − 1665. [XU Wenjie, HU Ruilin, YUE Zhongqi. Development of random mesostructure generating system of soil-rock mixture and study of its mesostructural mechanics based on numerical test[J]. Chinese Journal of Rock Mechanics and Engineering,2009,28(8):1652 − 1665. (in Chinese with English abstract) doi: 10.3321/j.issn:1000-6915.2009.08.017

    [15]

    徐文杰, 胡瑞林. 土石混合体概念、分类及意义[J]. 水文地质工程地质,2009,36(4):50 − 56. [XU Wenjie, HU Ruilin. Conception, classification and significations of soil-rock mixture[J]. Hydrogeology & Engineering Geology,2009,36(4):50 − 56. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-3665.2009.04.012

    [16]

    赵金凤, 严颖, 季顺迎. 基于离散元模型的土石混合体直剪试验分析[J]. 固体力学学报,2014,35(2):124 − 134. [ZHAO Jinfeng, YAN Ying, JI Shunying. Analysis of direct shear test of soil-rock mixture based on discrete element model[J]. Chinese Journal of Solid Mechanics,2014,35(2):124 − 134. (in Chinese with English abstract)

    [17]

    张昭, 魏星, 王刚. 颗粒长短轴比对饱和砂土液化影响的离散元分析[J]. 路基工程,2018(5):40 − 44. [ZHANG Zhao, WEI Xing, WANG Gang. DEM study on effect of particle long-short-axes ratio on liquefaction of saturated sand[J]. Subgrade Engineering,2018(5):40 − 44. (in Chinese with English abstract)

  • 期刊类型引用(4)

    1.  刘鑫,董广阳,史旦达. 考虑颗粒破碎的砂土中鱼雷锚贯入离散元分析. 水文地质工程地质. 2024(01): 91-101 . 本站查看
    2.  洪勇,姜奕辰,倪嘉楠,于超. 基于离散元方法的混合砂土剪切特性细观分析. 河南理工大学学报(自然科学版). 2024(01): 172-179 . 百度学术
    3.  杨志勇,朱俊伟,杨星,邵小康,江玉生,王允龙. 基于改进LCPC测试方法的盾构隧道土体磨蚀性试验. 水文地质工程地质. 2023(06): 90-98 . 本站查看
    4.  匡俊,陈芋如,朱仁民,曹建林,周俊,唐强,孙俊芳. 基于离散元法的边坡填土剪切性能分析. 地基处理. 2022(05): 380-386 . 百度学术

    其他类型引用(1)

  • 加载中
    Created with Highcharts 5.0.7访问量Chart context menu近一年内文章摘要浏览量、PDF下载量统计信息摘要浏览量PDF下载量2024-052024-062024-072024-082024-092024-102024-112024-122025-012025-022025-032025-04012345Highcharts.com
    Created with Highcharts 5.0.7Chart context menu访问类别分布DOWNLOAD: 3.6 %DOWNLOAD: 3.6 %摘要: 96.4 %摘要: 96.4 %DOWNLOAD摘要Highcharts.com
    Created with Highcharts 5.0.7Chart context menu访问地区分布其他: 5.6 %其他: 5.6 %其他: 0.8 %其他: 0.8 %Aurora: 2.1 %Aurora: 2.1 %Beijing: 0.6 %Beijing: 0.6 %Çanakkale: 0.8 %Çanakkale: 0.8 %Chicago: 0.4 %Chicago: 0.4 %Colombo: 0.6 %Colombo: 0.6 %Ghaziabad: 0.2 %Ghaziabad: 0.2 %Haidian: 2.5 %Haidian: 2.5 %Hanoi: 0.6 %Hanoi: 0.6 %Mountain View: 17.7 %Mountain View: 17.7 %Phoenix: 2.5 %Phoenix: 2.5 %Seattle: 0.2 %Seattle: 0.2 %[]: 0.6 %[]: 0.6 %上海: 0.8 %上海: 0.8 %东京都: 0.4 %东京都: 0.4 %北京: 0.4 %北京: 0.4 %台州: 0.2 %台州: 0.2 %呼和浩特: 0.2 %呼和浩特: 0.2 %宣城: 0.2 %宣城: 0.2 %延安: 0.4 %延安: 0.4 %张家口: 0.2 %张家口: 0.2 %徐州: 0.4 %徐州: 0.4 %成都: 0.6 %成都: 0.6 %扬州: 0.2 %扬州: 0.2 %洛杉矶: 0.4 %洛杉矶: 0.4 %湖州: 0.2 %湖州: 0.2 %漯河: 0.2 %漯河: 0.2 %石家庄: 0.6 %石家庄: 0.6 %芒廷维尤: 52.8 %芒廷维尤: 52.8 %苏州: 0.6 %苏州: 0.6 %莫斯科: 1.0 %莫斯科: 1.0 %衢州: 0.2 %衢州: 0.2 %西宁: 2.1 %西宁: 2.1 %运城: 0.4 %运城: 0.4 %郑州: 0.8 %郑州: 0.8 %长沙: 0.2 %长沙: 0.2 %阿斯图里亚斯自治区: 0.6 %阿斯图里亚斯自治区: 0.6 %其他其他AuroraBeijingÇanakkaleChicagoColomboGhaziabadHaidianHanoiMountain ViewPhoenixSeattle[]上海东京都北京台州呼和浩特宣城延安张家口徐州成都扬州洛杉矶湖州漯河石家庄芒廷维尤苏州莫斯科衢州西宁运城郑州长沙阿斯图里亚斯自治区Highcharts.com

(11)

(5)

计量
  • 文章访问数:  919
  • PDF下载数:  77
  • 施引文献:  5
出版历程
收稿日期:  2020-12-14
修回日期:  2021-02-04
刊出日期:  2021-11-15

目录

  • 表 1.  基于离散元模拟的直剪试验主要计算参数
    Table 1.  Main computational parameters in DEM simulation of direct shear
    计算参数砾石
    法向刚度kn/(5×1085×109
    切向刚度ks/(2×1082×109
    摩擦系数0.50.9
    粒径/mm0.5~5.0>5.0
     | Show Table
    DownLoad: CSV
  • 表 2.  不同砾石粒度分维值在直剪试验中获取的相关参数
    Table 2.  Correlation parameters of fractal dimension of different gravel sizes obtained in the direct shear test
    砾石粒度分维值100 kPa200 kPa300 kPa内摩擦角/(°)
    2.4568.57144.84196.0734.06
    2.5065.22147.33209.9735.26
    2.5569.56150.10220.0536.34
    2.6065.22123.23210.7833.99
    2.6564.84112.112195.4832.02
     | Show Table
    DownLoad: CSV
  • 表 3.  不同砂粒度分维值在直剪试验中获取的相关参数
    Table 3.  Relevant parameters obtained in the direct shear test of the fractal dimension of different sand particle sizes
    砂粒度分维值100 kPa200 kPa300 kPa内摩擦角/(°)
    2.4565.22143.75175.3132.11
    2.5066.09150.45202.7534.86
    2.5569.57150.10220.0536.34
    2.6068.70135.24203.5334.15
    2.6564.44134.01170.9931.13
     | Show Table
    DownLoad: CSV
  • 表 4.  不同轴向系数在直剪试验中获取的强度参数
    Table 4.  Strength parameters obtained in direct shear test with different axial coefficients
    主要参数轴向系数
    S=1.0S =1.4S =1.8
    100 kPa65.3168.9183.95
    200 kPa127.39128.25143.04
    300 kPa174.61184.58211.15
    400 kPa235.34251.42292.50
    内摩擦角/(°)30.7632.1435.92
    相关系数R20.99190.99680.9925
     | Show Table
    DownLoad: CSV
  • 表 5.  不同轴向系数下砾石接触数量和旋转角度
    Table 5.  Gravel contact number and rotation angle under different axial coefficients
    轴向系数砾石接触数量砾石平均旋转角度/(°)
    1.015819.43
    1.416568.70
    1.817426.42
     | Show Table
    DownLoad: CSV