含裂隙岩体单轴压缩裂纹扩展机制离散元分析EI北大核心CSCD

将由室内试验总结得到的岩石微观胶结模型嵌入离散元软件,对Lac du Bonnet花岗岩石进行预制单裂隙单轴压缩试验DEM数值模拟,分析了压缩过程中裂隙试样中应力的分布,并与理论计算结果进行对比分析,同时对各种断裂判据中裂纹起裂角的预测值进行了适用性的对比分析。结果表明,离散元模拟试样破坏形态与试验结果相近;离散元分析得到的应力分布与理论解在定性上相似;当预制角度较小时,侧向应力都处于拉压状态;由于裂隙左右两端压应变的集中造成了裂隙上下面拉应变的产生,造成了裂隙周围特殊的应力分布;当裂隙角度较大时,应力集中现象已不明显,因而,理论值与试验值有偏差;在断裂判据中最大周应力准则和最大能量释放率准则得到的裂纹扩展角与室内试验与DEM结果中的数值较为吻合。     

简化接触模型的月壤离散元数值分析

根据月壤其颗粒级配可归类于粉质砂土。针对真实月壤所处的环境（无水、低重力场、低气压等）,将Perko等2001年提出的月壤颗粒间的范德华力植入离散元分析软件PFC2D中,模拟了刚性边界下加入该模型与未加该模型试样的双轴压缩试验,研究了颗粒间范德华力对试样的宏观力学特性与微观颗粒接触的影响。结果表明,颗粒间的范德华力对试样的抗剪强度、体应变以及颗粒平均配位数都有显著的影响     

离散元微观参数对砂土宏观参数的影响

采用离散元商业软件PFC2D对密砂和松砂试样共进行了66组双轴压缩试验,试验中改变颗粒接触刚度和粒间摩擦系数以探究离散元微观参数对试样宏观参数(弹性模量和内摩擦角)的影响,其中颗粒法向刚度变化范围为1.5×108～1.5×109N.m-1,粒间摩擦系数变化范围为0.1～0.7 N.m-1。结果表明:颗粒接触刚度对弹性模量有一定程度影响且对密样影响程度大于松样,而对内摩擦角基本无影响;粒间摩擦系数对试样弹性模量有明显影响且对松样的影响程度更大,对密样和松样的内摩擦角都有非常显著的影响。     

含裂隙岩体单轴压缩裂纹扩展机制离散元分析

将由室内试验总结得到的岩石微观胶结模型嵌入离散元软件,对Lac du Bonnet花岗岩石进行预制单裂隙单轴压缩试验DEM数值模拟,分析了压缩过程中裂隙试样中应力的分布,并与理论计算结果进行对比分析,同时对各种断裂判据中裂纹起裂角的预测值进行了适用性的对比分析。结果表明,离散元模拟试样破坏形态与试验结果相近;离散元分析得到的应力分布与理论解在定性上相似;当预制角度较小时,侧向应力都处于拉压状态;由于裂隙左右两端压应变的集中造成了裂隙上下面拉应变的产生,造成了裂隙周围特殊的应力分布;当裂隙角度较大时,应力集中现象已不明显,因而,理论值与试验值有偏差;在断裂判据中最大周应力准则和最大能量释放率准则得到的裂纹扩展角与室内试验与DEM结果中的数值较为吻合。     

不同开采方法下深海能源土离散元模拟

天然气水合物以胶结及孔隙填充等形式存在于深海能源土中,开采时因其分解会劣化地层力学特性进而引发海底事故,使得人们对能源土开采过程进行中力学特性的变化愈发重视。在前期室内试验的基础上,将一个温度-水压-力学二维微观胶结模型引入离散元商业软件PFC2D中,通过对排气、排水性较好的土体进行升温及降压法开采进行数值模拟,并将模拟结果与相同条件下的室内试验结果对比,验证了该胶结模型的适用性。进一步分析了颗粒接触分布与颗粒平均纯转动率（averaged pure rotation rate, APR）在水合物分解时的变化情况。升温分解过程中随温度升高,颗粒总接触分布各向异性程度增大;胶结接触逐渐减少并始终保持主方向为水平方向,无胶结接触增多并始终保持主方向为竖直方向;APR值逐渐增大且正负值分布逐渐趋于集中。降压分解过程中随反（水）压降低,颗粒总接触由各向同性分布逐渐发展为主方向为竖直方向的各向异性,APR值较小且分布均匀;恢复反压后,试样进一步破坏,颗粒总接触各向异性更加明显,APR值增大且正负值呈集中分布。     

应力路径试验前后不同黄土的孔隙分形特征

为了解不同黄土孔隙形状复杂程度的差异和应力路径对孔隙形状复杂程度的影响,对两种黄土应力路径前后孔隙分形特征进行了研究。首先比较了3种分维模型所得孔隙分形维数的可靠性,然后选用热力学关系模型,由进汞、退汞试验得到地裂缝区黄土、充填黄土初始样和三轴应力路径试验后的孔隙分布,据此分析了两种土体进汞孔隙和退汞孔隙分形维数的差异和受载后的分形维数变化,根据退汞过程仅管形孔内的汞流出及孔隙由管形孔和球形孔组成的假定,得到了应力路径试验前后孔隙形状的改变。结果表明,热力学关系模型得到的孔隙分形维数合理、可靠。初始状态,原状和重塑充填黄土的孔隙形状比原状和重塑黄土复杂;常规三轴试验后试样孔隙比减少,孔隙分形维数增加,球形孔向管形孔转变;相对于常规三轴压缩试验,减围压三轴压缩试验后试样的孔隙分形维数较小,管形孔占总孔隙的体积比例较少;总体上管形孔的分形维数比球形孔大,且基本不受应力路径的影响,其占总孔隙的体积比例随着试样宏观孔隙比的减少而增加。