应力波在各向异性介质中传播及损伤数值模拟.pdfVIP

第５０卷第１期山西建筑Ｖｏｌ．５０Ｎｏ．１

２０２４年１月Ｊａｎ．２０２４·７７·

ＳＨＡＮＸＩＡＲＣＨＩＴＥＣＴＵＲＥ

１０．１３７１９／ｊ．ｃｎｋｉ．１００９－６８２５．２０２４．０１．０２０

ＤＯＩ：

应力波在各向异性介质中传播及损伤数值模拟★

１１１２２

鞠冬冬，王清，陈鹏，胡少斌，蔡余康

（１．南京市江北新区公共工程建设中心，江苏南京２１００９８；２．河海大学土木与交通学院，江苏南京２１００９８）

摘要：基于ＣＯＭＳＯＬ的相场损伤对应力波在各向异性岩层中的传播和所造成的损伤进行数值模拟研究，使用无限元边界来

避免固定边界下应力波反射造成的影响，通过对爆破孔孔间距以及孔心连线的转角进行参数化扫描以研究二者对应力波传

播和损伤演化的影响。最终得到应力波传播和损伤演化云图，以及应力时程曲线和破坏区面积（损伤＞０．９）变化曲线并得到

以下结论：应力波在弹性模量更小的方向应力幅值较高；轴线转角减小会延迟中心处应力波第一个波峰的时间，且峰值应力

减小；孔间距减小会导致更大的峰值应力和更明显的叠加效应，但会导致破碎区面积的减小。

关键词：有限元；应力波；各向异性；无限元；损伤；相场法

中图分类号：ＴＵ４５９．４文献标识码：Ａ１００９－６８２５（２０２４）０１－００７７－０５

文章编号：

在人类的生活中经常有诸如冲击、爆炸、碰撞等强荷下爆破花岗岩层的数值模拟模型，并给出本模型的相关

载作用下材料发生高应变的现象。比如地震、洪水、交通控制方程、几何建模、网格划分、边界条件，以及模型所涉

事故这样的灾难，还有矿山开采、钻爆法开挖、结构抗震及的材料参数。

等造福民生的技术方法。这些现象中诸如地震、矿山开１．１控制方程

采以及钻爆法开挖都可纳入岩石力学范畴，更具体地讲，本模型采用相场损伤法来描述线弹性岩石的损伤应

－１４－１

属于岩石动力学范畴，即应变率在１×１０～１×１０ｓ变软化过程，其本构方程采用胡克定律：

［１］

之间的情况。而岩石动力学以应力波作为其理论基σ＝Ｄ×ε（１）

础，主要研究涵盖应力波在岩体中的传播与衰减［２－４］，应其中，为弹性应力；Ｄ为考虑损伤的弹性矩阵；为

σε

力波与节理裂隙的相互作用［５－７］，应力波与结构面的透、弹性应变。

折、反射关系［８－１０］，应力波入射角［１１－１２］等多个方面。除本文研究应力波在各向异性介质中的传播，故弹性

此之外，岩石动力学所研究的另一个基本问题是岩石或矩阵Ｄ自定义为：

岩体受冲击后的动态响应，包括变形、损伤、断裂的结果以１６５５０５０

及演化过程。目前，对应力波传播及岩石的动态响应的研５０６２５００

究多种多样，从一开始常规的应力波的传播和衰减规５０５０６２

［１３］［１４－１５］