SIMULATION FOR 2D COMPRESSION-SHEAR JOINT DYNAMIC PROPAGATION PROCESS USING MESHLESS METHODS BASED ON
LEVEL SETS COORDINATES
ZHUANG Xiaoying1，2，3，HUANG Runqiu3，ZHU Hehua1，2
(1. Department of Geotechnical Engineering，Tongji University，Shanghai 200092，China；2. Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education，Tongji University，Shanghai 200092，China；3. State Key Laboratory of Geohazard Prevention and Geoenvironment Protection，Chengdu University of Technology，Chengdu，Sichuan 610059，China)
Abstract：The contract and frictional algorithms are developed in the meshless method，element-free Galerkin method for modelling 2D crack propagation in rock mass subjected to compressive and tensile loading conditions. Here，the level set method is used to describe the crack geometry. To capture the displacement jump across the joint when the join is open，the diffraction rule is used to modify the weight function；and hence construction of the discontinuous meshless approximation. The joint closure is modeled by a new meshless frictional sliding and contact algorithm，which is similar to the well known Goodman element in the family of finite elements to model the sliding and contact between joint surfaces. Different from Goodman element and former meshless methods，the construction of present method is based on the material points instead of the nodes，and thereby it performs better and is more flexible in model generation of a joint of complex geometry. The material points and nodes generation after joint propagation is based on the level sets functions which have a general form for arbitrary crack in 2D. A number of examples are tested to model the crack propagation in rock samples subjected to compression and tension showing the correctness of the method.
[1] LAJTAI E Z. Shear strength of weakness planes in rock[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts，1969，6(5)：499–515.
[2] OLSON J E. Joint pattern development：effects of subcritical crack growth and mechanical crack interaction[J]. Journal of Geophysical Research，1993，98(B7)：12 251–12 265.
[3] 李术才，陈卫忠，朱维申，等. 加锚节理岩体裂纹扩展失稳的突变模型研究[J]. 岩石力学与工程学报，2003，22(10)：1 661–1 666.(LI Shucai，CHEN Weizhong，ZHU Weishen，et al. Catastrophe model of crack propagating in bolt-supported jointed rock mass[J]. Chinese Journal of Rock Mechanics and Engineering，2003，22(10)：1 661– 1 666.(in Chinese))
[4] 周小平，王建华，哈秋舲. 压剪应力作用下断续节理岩体的破坏分析[J]. 岩石力学与工程学报，2003，22(9)：1 437–1 440.(ZHOU Xiaoping，WANG Jianhua，HA Qiuling. Failure analysis of intermittent jointed rocks under compression-shear stress[J]. Chinese Journal of Rock Mechanics and Engineering，2003，22(9)：1 437– 1 440.(in Chinese))
[5] LI Y P，CHEN L Z，WANG Y H. Experimental research on pre-cracked marble under compression[J]. International Journal of Solids and Structures，2005，42(9/10)：2 505–2 516.
[6] 刘远明，夏才初，李宏哲. 节理研究进展及在非贯通节理岩体研究的应用[J]. 地下空间，2007，3(4)：682–693.(LIU Yuanming，XIA Caichu，LI Hongzhe. Development in joints research and its application to rock mass containing discontinuous joints[J]. Chinese Journal of Underground Space and Engineering，2007，3(4)：682–693.(in Chinese))
[7] JIANG L. A review of techniques advances and outstanding issues in numerical modelling for rock mechanics and rock engineering[J]. International Journal of Rock Mechanics and Mining Sciences，2003，40(3)：283–353.
[8] BELYTSCHKO T，BLACK T. Elastic crack growth in finite elements with minimal remeshing[J]. International Journal for Numerical Methods in Engineering，1999，45(5)：601–620.
[9] NGUYEN V，RABCZUK T，BORDAS S，et al. Meshless methods：a review and computer implementation aspects[J]. Mathematics and Computers in Simulation，2008，79(3)：763–813.
[10] 庄晓莹，蔡永昌，朱合华. 锚固边坡稳定性分析有限元模型及锚固效应的探讨[J]. 岩土工程学报，2008，30(7)：1 099–1 104. (ZHUANG Xiaoying，CAI Yongchang，ZHU Hehua. FEM modeling method for anchored slope stability analysis and discussions on anchoring effect[J]. Chinese Journal of Geotechnical Engineering，2008，30(7)：1 099–1 104.(in Chinese))
[11] 雷远见，王水林. 基于离散元的强度折减法分析岩体边坡稳定性[J]. 岩土力学 2006，27(10)：1 693–1 698.(LEI Yuanjian，WANG Shuilin. Slope stability analysis of jointed rock slope by strength reduction method based on UDEC[J]. Rock and Soil Mechanics，2006，27(10)：1 693–1 698.(in Chinese))
[12] 焦玉勇，张秀丽，刘泉声，等. 用非连续变形分析方法模拟岩石裂纹扩展[J]. 岩石力学与工程学报，2007，26(4)：682–691.(JIAO Yuyong，ZHANG Xiuli，LIU Quansheng，et al. Simulation of rock crack propagation using discontinuous deformation analysis method[J]. Chinese Journal of Rock Mechanics and Engineering，2007，26(4)：682–691.(in Chinese))
[13] 甯尤军，杨 军，陈鹏万. 节理岩体爆破的DDA方法模拟[J]. 岩土力学，2010，31(7)：2 259–2 263.(NING Youjun，YANG Jun，CHEN Pengwan. Numerical simulation of rock blasting in jointed rock mass by DDA method[J]. Rock and Soil Mechanics，2010，31(7)：2 259–2 263.(in Chinese))
[14] 庄晓莹，朱合华，蔡永昌，等. 基于无网格–图论方法的岩体边坡稳定性分析研究[J]. 力学季刊，2008，29(4)：537–543.(ZHUANG Xiaoying，ZHU Hehua，CAI Yongchang，et al. Rock slope stability analysis based on meshless method and shortest path algorithm[J]. Chinese Quarterly of Mechanics，2008，29(4)：537–543.(in Chinese))
[15] 陈 炜，王均星. 节理岩质边坡的块体元塑性极限分析下限法[J]. 岩土工程学报，2008，30(2)：272–277.(CHEN Wei，WANG Junxing. Plastic lower bound limit analysis using block element method for jointed rock slope[J]. Chinese Journal of Geotechnical Engineering，2008，30(2)：272–277.(in Chinese))
[16] 赵尚毅，郑颖人，邓卫东. 用有限元强度折减法进行节理岩质边坡稳定性分析[J]. 岩石力学与工程学报，2003，22(2)：254–260. (ZHAO Shangyi，ZHENG Yingren，DENG Weidong. Stability analysis of jointed rock slope by strength reduction FEM[J]. Chinese Journal of Rock Mechanics and Engineering，2003，22(2)：254–260.(in Chinese))
[17] ZHU H H，ZHUANG X Y，CAI Y C，et al. High rock slope stability analysis using the enriched meshless shepard and least squares method[J]. International Journal of Computational Methods，2011，8(2)：209–228.
[18] ZHUANG X Y，AUGARDE C，BORDAS S. Accurate fracture modelling using meshless methods and level sets：formulation and 2D modelling[J]. International Journal for Numerical Methods in Engineering，2011，86(2)：249–268.
[19] OSHER S，FEDKIW R. Level set methods and dynamic implicit surfaces[M]. New York：Springer，2003：63–74.
[20] DUFLOT M，NGUYEN-DANG H. A meshless method with enriched weight functions for fatigue crack growth[J]. International Journal for Numerical Methods in Engineering，2004，59(14)：1 945–1 961.
[21] RIGBY R，ALIABADI M. Decomposition of the mixed-mode J-integral revisited[J]. International Journal of Solids and Structures，1998，35(17)：2 073–2 099.
[22] MURAKAMI Y. The stress intensity factors handbook[M]. Oxford：Pergamon Press，1987：1 255–1 256.
[23] 朱维申，李术才，陈卫忠. 节理岩体破坏机制和锚固效应及工程应用[M]. 北京：科学出版社，2002：116–146.(ZHU Weishen，LI Shucai，CHEN Weizhong. The failure mechanism and anchorage effect of jointed rock mass and their application to engineering[M]. Beijing：Science Press，2002：116–146.(in Chinese))
[24] MUGHIEDA O，ALZO′UBI A K. Fracture mechanisms of offset rock joints—a laboratory investigation[J]. Geotechnical and Geological Engineering，2004，22(4)：545–562.
2012, Vol. 31 
