单位文秘网 2021-08-30 08:57:11 点击: 次
打开文本图片集
引言
离散元方法是当今非连续介质力学数值方法中应用于岩石和混凝土等脆性材料的最广泛而且最有效的方法之一。其最大特点是可以反映岩石块体之间的接触面滑移、分离以及倾覆旋转,表征岩体内部节理、裂隙和断层等非连续特性,同时还能够计算块体内部的变形和应力,尤其对于描述岩石破坏过程的细观机制具有内在的优越性。已经广泛应用于岩石边坡稳定性分析、地下厂房围岩稳定性和破坏过程模拟、混凝土结构的裂缝扩展、破坏与倒塌过程分析、盾构机滚刀破岩过程分析和混凝土大坝及其基础破坏机制等岩石工程领域。Gu[1]提出了刚体-弹簧模型模拟混凝土材料初始缺陷和骨料分布特性,所建立的离散元模型包括粗骨料、细骨料、砂浆和粘结界面,并且考虑了这些材料的随机分布特性。James[2]研究了非均质岩石细观裂纹的分布特性,提出了颗粒断裂的破坏准则。Gray[3]回顾了地质结构体的有限元和离散元建模问题,分析了有限元、离散元和有限差分法的在表征非连续地质体力学特性的适用性。Shmulevich[4]应用离散元方法(PFC)模拟了切刀与土体之间的相互作用特性,研究表明该方法模拟的结果与实验结果基本一致,其目的在于通过优化切刀的形状,提高切刀破土的能量利用率。Ono[5]采用三维离散元方法模拟了切土工具与土体之间的作用机制,讨论了6种离散元形状对数值模拟结果的影响,并将实验结果与模拟结果进行了对比分析。Ucgul[6]采用三维离散元模型模拟了切刀切土的破坏过程,重点分析了无内聚力散体砂土模型参数估计问题,讨论了线性和非线性接触模型的差异性和适用性。Arena[7]提出了一种新方法用于定量计算实验过程中岩石细观裂纹,该方法以二维数字图像为基础,对于定量分析岩石试件表面裂纹的演化过程具有借鉴和参考价值。Mak[8]介绍了用于表征土体与切刀之间相互作用的离散元模型,讨论了细观模型中某些参数的简化估计方法,并且进行了切刀切土的实验和模拟分析研究。Asaf[9]提出了如何根据切刀破土实验数据进行离散元模型参数校准问题。Innaurato[10]采用实验和数值模拟方法研究了TBM滚刀的破岩过程,建立了滚刀贯入深度与贯入力之间的关系,研究了滚刀作用下岩石的剪切破坏问题。从以上相关研究工作中可以看出,采用离散元方法研究土体材料破坏过程所进行的工作较多,而研究岩石材料细观破坏过程的工作相对较少。与此同时,如何较为准确估计岩石材料的细观模型参数是制约离散元方法解决工程实际问题的瓶颈。岩石的力学行为由其内在的细观结构和模型参数决定,尤其是岩石细观裂纹及其演化和发展过程。笔者的目的在于在岩石颗粒尺度上研究压头作用下的岩石破坏过程,揭示岩石破碎漏斗的细观形成机制。
1岩石材料平行粘结模型细观参数估计
岩石是多种矿物和胶结物组成的一种高度不均匀材料,在细观层次上表现在颗粒、粘结物和结构面等力学特性的随机分布特点。颗粒流(Particle Flow Code)方法是由Cundall和Strack提出的模拟球形颗粒运动与相互作用的细观离散元方法。该方法已经应用到岩石、土体和混凝土材料的力学特性研究中,着重从细观力学角度揭示材料的损伤和断裂机理。颗粒接触刚度模型包括线性颗粒接触刚度模型和HerzMindlin非线性接触模型等。如图1所示,颗粒接触刚度模型参数包括法向刚度(kn)、切向刚度(ks)和摩擦系数(f)以及阻尼(η)等。
用于模拟岩石和混凝土材料(非散体材料)的开裂问题,Azevedo提出了其简化模式[11],如图2所示,岩石的细观模型参数还包括:颗粒之间胶结材料的法向刚度(kpn)和切向刚度(kps)、胶结材料切向剪切强度(τs)和法向抗拉强度(σt),以及胶结圆环的半径因子(λ)。需要说明的是,颗粒之间胶结材料的切向剪切强度(τs)和法向抗拉强度(σt)与岩石宏观实验的剪切强度和抗拉强度并不是一个概念,但二者之间应该具有某种相关性,因为平行粘结模型将颗粒之间的粘结简化成了一个等效的圆盘。
平行粘结模型法向刚度kn为[12]
kn=2tEc,(1)
式中t为宽度,对于2维模型,t=1;Ec为颗粒接触变形模量,近似取为岩石宏观实验的弹性模量E,取Ec=43 GPa.切向刚度ks与法向刚度的关系为
ks=Rckn,(2)
式中Rc为颗粒切向刚度与法向刚度的比,研究表明,这两个值近似相等[13],取Rc=1.0.颗粒之间的细观摩擦系数f与宏观内摩擦角的近似关系为[14]
f=2πtg,(3)
颗粒之间胶结材料的法向刚度kpn为
kpn=Epc(RA+RB),(4)
式中Epc为平行粘结变形模量;RA和RB分别为颗粒A和B的半径。平行粘结颗粒之间胶结材料的切向刚度kps为
kps=Rpkpn,(5)
式中Rp为平行粘结模型切向刚度与法向刚度的比,取Rp=1.平行粘结模型细观抗拉强度与宏观内聚力c和内摩擦角的关系为[8]
σt=c×cot,(6)
式中σt为平行粘结模型细观法向抗拉强度。平行粘结细观抗剪强度与宏观内聚力的关系为[8]
τs=c,(7)
式中τs为颗粒之间胶结材料的切向剪切强度;颗粒之间固定平行粘结半径因子取为λ=0.5.颗粒之间粘结环半径Rb为
Rb=λmin(RA,RB).(8)
岩石是由多种矿物和胶结物组成的非均匀体,并且含有大量的天然节理和微裂隙。天然岩体中的断层、节理和裂隙的随机分布,其空间几何和力学参数的分布具有统计特征。在长期的地质作用和水文环境下岩体产生的微裂纹、缺陷和孔隙,在外界荷载作用下,在岩石矿物颗粒和胶结物结合面上存在初始微裂纹损伤,在破坏机制上存在局部化现象,并导致应变软化的本构特性[14]。
在单轴压缩情况下,岩石试样的破坏过程是随机分布强度参数的岩石颗粒的刚度逐渐劣化的结果。从连续介质力学观点分析,岩石试样的破坏方式就其力学本质可以归结为压缩破坏、拉伸破坏和剪切破坏。为了实现表征颗粒连接力学破坏模式,应用FISH语言编写一个伺服程序,检测每个颗粒连接的受力状态,依据所建立的颗粒连接力学破坏准则,并对每个颗粒连接进行破断与否的评估,进而实现描述颗粒连接的破坏功能。
under 4 mm penetration dept在平行粘结力链图中,深色部分表示颗粒之间处于受拉状态,黑色部分表示颗粒之间处于受压状态。从图5中可以看出,在压头垂直压力作用下,在岩石试件内,形成了3条主环向拉力链,它与宏观的小主应力方向一致,环向拉力链最大影响深度约为30 mm.同时,出现了几条径向压力链,它与宏观的大主应力方向基本一致。
从图6中可以看出,随着压头侵入深度的增加,环向拉力链最大影响深度扩展到约为42 mm;主环向拉力链发展得更加明显。
10 mm penetration depth从图9中可以发现,在压头正下方形成了一个明显的破碎漏斗(浅色区域),破碎漏斗的角度为70左右。岩石是一种典型的非均质和准脆性材料[16],岩石晶粒和初始缺陷的分布具有随机性和不确定性,在荷载作用下,各组分成分的刚度差异很大,使得对力的传递响应和自身变形的差异,导致在岩石内部应力场和应变场分布极不均匀,产生局部应力集中,在初始缺陷或弱结合面区域产生微裂纹,且裂纹的扩展路径具有明显的非对称性。岩石中裂纹扩展路径强烈依赖于岩石材料的非均匀性,材料的非均匀性对裂纹的贯通和抑制有着重要的影响。在细观尺度上,岩石内部所含有的晶粒和胶结物的力学性质差异较大,岩石材料具有明显的非均匀性。
3结论
1)岩石材料在弹塑性变形直至破坏过程中,其宏观体积响应只是各种微结构之间相互作用的平均结果,变形的主要机制是微裂纹成核和增长及其聚集而生成细观裂纹、扩展成宏观裂纹、直至宏观破坏。由于岩石(岩体)内含有大量的初始缺陷(节理、层理、微裂纹等),连续介质力学方法对于表征岩石的这些非连续特性存在着固有的缺点,尤其是分析裂纹群扩展过程时遇到的困难更加突出,离散元方法为解决这类难题提供了新的途径;
2)采用非连续介质力学数值方法,发展岩石的细观力学模型有助于从材料破坏机理的层面上对宏观岩石破坏过程进行理解和解释。研究表明,在压头垂直压力作用下,在岩石试件内,形成了3条主环向拉力链,它与宏观的小主应力方向一致。环向拉力链最大影响深度约随着压头压入深度的增加而增加。与此同时,出现了几条径向压力链,它与宏观的大主应力方向基本一致;
3)岩石试件裂纹是由环向拉应力和剪应力大于颗粒平行连接胶结材料的强度引起的。压头作用下岩石的破坏特性细观模拟分析成果,为TBM滚刀破岩分析、滚刀优化布置和滚刀形状设计提供了一种可以借鉴的方法;
4)如何精确定量确定岩石细观模型中的参数是制约离散元方法发展和应用的瓶颈问题之一。由于岩石颗粒大多在mm量级,在该尺度上直接测量力学参数是非常困难的。如果能够在岩石试样宏观尺度上进行实验,并且基于宏观实验数据进行力学细观模型参数反演,将是一项十分有意义的工作。
参考文献References
[1]GU Xianglin,LI Hong,WANG Zhuolin,et al.A modified rigidbodyspring concrete model for prediction of initial defects and aggregates distribution effect on behavior of concrete[J].Computational Materials Science,2013,77(3):355-365.
[2]James P M,Ford M,Jivkov A P.A novel particle failure criterion for cleavage fracture modeling allowing measured brittle particle distributions[J].Engineering Fracture Mechanics,2014,121-122(1):98-115.
[3]Gray G G,Morgan J K,Sanz P F.Overview of continuum and particle dynamics methods for mechanical modeling of contractional geologic structures[J].Journal of Structural Geology,2014,59(2):19-36.
[4]Shmulevich I,Asaf Z,Rubinstein D,Interaction between soil and a wide cutting blade using the discrete element method[J].Soil & Tillage Research,2007,97(1):37-50.
[5]Ono I,Nakashima H,Shimizu H,et al.Investigation of elemental shape for 3D DEM modeling of interaction between soil and a narrow cutting tool[J].Journal of Terramechanics,2013,50(4):265-276.
[6]Ucgul M,Fielke J M,Saunders.Threedimensional discrete element modeling of tillage:Determination of a suitable contact model and parameters for a cohesionless soil[J].Biosystems Engineering,2014,121(2):105-117.
[7]Arena A,Piane C D,Sarout J.A new computational approach to cracks quantification from 2D image analysis:Application to microcracks description in rocks[J].Computers & Geoscience,2014,66(1):106-120.
[8]Mak J,Chen Y,Sadek M A.Determining parameters of a discrete element model for soiltool interaction[J].Soil & Tillage Research,2012,118(5):117-122.
[9]Asaf Z,Rubinstein D,Shmulevich I.Determination of discrete element model parameters required for soil tillage[J].Soil & Tillage Research,2007,92(1-2):227-242.
[10]Innaurato N,Oggeri C,Oreste P P,et al.Experimental and numerical studies on rock breaking with TBM tools under high stress confinement rock mech[J].Rock Mechanics and Rock Engineering,2007,40(5):429-451
[11]Azevedo N M,Lemos J V,Almeida J R.Influence of aggregate deformation and contact behavior on discrete particle modeling of fracture of concrete[J].Engineering Fracture Mechanics,2008,75(6):1 569-1 586.
[12]Potyondy D O,Cundall P A.A bondedparticle model for rock[J].International Journal of Rock Mechanics and Mining Sciences,2004,41(8):1 329-1 364.
[13]李守巨,李德,于申.基于宏观实验数据堆石料细观本构模型参数反演[J].山东科技大学学报,2015,34(5):20-26.
LI Shouju,LI De,YU Shen.Mesoparameters Inversion of constitutive model for rockfill materials based on micro experimental data[J].Journal of Shandong University of Science and Technology,2015,34(5):20-26.
[14]孙其诚,厚美瑛,金峰.颗粒物质物理与力学[M].北京:科学出版社,2011.
SUN Qicheng,HOU Meiying,JIN Feng.Physical and mechanicals of granular material[M].Beijing:Science Press,2011.
[15]乔卫国,杨麟,李大勇.高频振动沉桩的离散元模拟分析[J].西安科技大学学报,2012,32(5):604-609.
QIAO Weiguo,YANG Lin,LI Dayong.Numerical analysis of piles driving by high frequency vibration using discrete element method[J].Journal of Xi’an University of Science and Technology,2012,32(5):604-609.
[16]马文伟,赵光明,孟祥瑞.岩石SHPB 劈裂分析理论研究[J].西安科技大学学报,2013,33(5):283-290.
(责任编辑:单位文秘网) )地址:https://www.kgf8887.com/show-244-90136-1.html
上一篇:认知行为护理干预对改善乳腺癌术后患者焦虑抑郁情绪的影响
下一篇:流体运动的几类模型
版权声明:
本站由单位文秘网原创策划制作,欢迎订阅或转载,但请注明出处。违者必究。单位文秘网独家运营 版权所有 未经许可不得转载使用