|
|
|
| APPLICATION OF THE LINEARLY INDEPENDENT HIGH-ORDER NUMERICAL MANIFOLD METHOD IN FRACTURE MECHANICS |
| XU Dongdong1,YANG Yongtao2,ZHENG Hong2,WU Aiqing1 |
| (1. Key Laboratory of Geotechnical Mechanics and Engineering of Ministry of Water Resources,Yangtze River Scientific Research Institute,Hubei,Wuhan 430010,China;2. State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,Hubei,Wuhan 430071,China) |
|
|
|
|
Abstract The numerical manifold method(NMM) has succeeded in providing a unified solution to continuum and discontinuum problems and therefore it is highly suitable for solving fracture mechanics problems. However,the conventional high-order NMM using the first order polynomial as the local displacement function has the problem of linear dependence,which restricts to a certain degree its further development and application. A new NMM framework was established in this research by introducing a new localized displacement function,as well as a special displacement function for modeling the stress singularity around crack tips. A new paradigm that eliminates the problem of linear dependence is then derived to solve linear elastic fracture mechanics problems. The numerical examples show that:(1) The proposed method successfully eliminates the problem of linear dependence;(2) For classic linear elastic fracture problems,the stress intensity factors at the crack tip can be calculated accurately even if the mesh is relatively sparse;(3) The stress function at interpolation points inside the physical domain is continuous;(4) All the degrees of freedom defined on non-singular physical patches are physically meaningful,with the third to the fifth being the strain components at the interpolation point of the patch. As a result,the stress components at the interpolation point can be directly obtained. Finally,the proposed approach can be easily extended to other methods based on the theory of the partition of unity.
|
|
|
|
|
|
| [1] 张楚汉. 论岩石、混凝土离散–接触–断裂分析[J]. 岩石力学与工程学报,2008,27(2):217–235.(ZHANG Chuhan. Discrete-contact- fracture analysis of rock and concrete[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(2):217–235.(in Chinese))
[2] GERRARD C M. Joint compliances as a basis for rock mass properties and the design of supports[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts,1982,19(6):285–305.
[3] GOODMAN R E,JOHN C ST. Finite element analysis for discontinuous rocks[J]. Numerical Methods in Geotechnical Engineering,1977:148–175.
[4] BABUŠKA I,MELENK J M. The partition of unity method[J]. International Journal for Numerical Methods in Engineering,1997,40(4):727–758.
[5] 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.
[6] DOLBOW J,BELYTSCHKO T. A finite element method for crack growth without remeshing[J]. International Journal for Numerical Methods in Engineering,1999,46(1):131–150.
[7] SUKUMAR N,MOËS N,MORAN B,et al. Extended finite element method for three-dimensional crack modelling[J]. International Journal for Numerical Methods in Engineering,2000,48(11):1 549–1 570.
[8] MOËS N,GRAVOUIL A,BELYTSCHKO T. Non-planar 3D crack growth by the extended finite element and level sets-part I:Mechanical model[J]. International Journal for Numerical Methods in Engineering,2002,53(11):2 549–2 568.
[9] GRAVOUIL A,MOËS N,BELYTSCHKO T. Non-planar 3D crack growth by the extended finite element and level sets-Part II: Level set update[J]. International Journal for Numerical Methods in Engineering,2002,53(11):2 569–2 586.
[10] STROUBOULIS T,BABUŠKA I,COPPS K. The design and analysis of the generalized finite element method[J]. Computer Methods in Applied Mechanics and Engineering,2000,181(1):43–69.
[11] STROUBOULIS T,COPPS K,BABUŠKA I. The generalized finite element method:an example of its implementation and illustration of its performance[J]. International Journal for Numerical Methods in Engineering,2000,47(8):1 401–1 417.
[12] DUARTE C A,HAMZEH O N,LISZKA T J,et al. A generalized finite element method for the simulation of three-dimensional dynamic crack propagation[J]. Computer Methods in Applied Mechanics and Engineering,2001,190(15):2 227–2 262.
[13] DUARTE C A,RENO L G,SIMONE A. A high-order generalized FEM for through-the-thickness branched cracks[J]. International Journal for Numerical Methods in Engineering,2007,72(3):325–351.
[14] BELYTSCHKO T,KRONGAUZ Y,ORGAN D,et al. Meshless methods:an overview and recent developments[J]. Computer Methods in Applied Mechanics and Engineering,1996,139(1):3–47.
[15] BELYTSCHKO T,GU L,LU Y Y. Fracture and crack growth by element free Galerkin methods[J]. Modelling and Simulation in Materials Science and Engineering,1994,2(3A):519–534.
[16] BELYTSCHKO T,LU Y Y,GU L. Element-free Galerkin methods[J]. International Journal for Numerical Methods in Engineering,1994,37(2):229–256.
[17] CUNDALL P A. A computer model for simulating progressive large-scale movements in blocky rock systems[C]// Proceedings of the International Symposium Rock Fracture. Nancy,France:[s.n.],1971:No. II–8.
[18] CUNDALL P A. The measurement and analysis of acceleration on rock slopes[Ph. D. Thesis][D]. London:University of London,Imperial College of Science and Technology,1971.
[19] SHI G H,GOODMAN R E. Two-dimensional discontinuous deformation analysis[J]. International Journal for Numerical and Analytical Methods in Geomechanics,1985,9(6):541–556.
[20] SHI G H,GOODMAN R E. Generalization of two-dimensional discontinuous deformation analysis for forward modelling[J]. International Journal for Numerical and Analytical Methods in Geomechanics,1989,13(4):359–380.
[21] CAMONES L A M,VARGAS JR E A,DE FIGUEIREDO R P,et al. Application of the discrete element method for modeling of rock crack propagation and coalescence in the step-path failure mechanism[J]. Engineering Geology,2013,153:80–94.
[22] 焦玉勇,张秀丽,刘泉声,等. 用非连续变形分析方法模拟岩石裂纹扩展[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))
[23] NING Y,YANG J,MA G,et al. Modelling rock blasting considering explosion gas penetration using discontinuous deformation analysis[J]. Rock Mechanics and Rock Engineering,2011,44(4):483–490.
[24] SHI G H. Manifold method of material analysis[C]// Transactions of the Ninth Army Conference on Applied Mathematics and Computing. Minneapolish,Minncsoda,USA:[s. n.],1992:51–76.
[25] 姜清辉,周创兵. 四面体有限单元覆盖的三维数值流形方法[J]. 岩石力学与工程学报,2005,24(24):4 455–4 460.(JIANG Qinghui,ZHOU Chuangbing. Three-dimensional numerical manifold method with tetrahedron finite element covers[J]. Chinese Journal of Rock Mechanics and Engineering,2005,24(24):4 455–4 460.(in Chinese))
[26] 姜清辉,邓书申,周创兵. 三维高阶数值流形方法研究[J]. 岩土力学,2006,27(9):1 471–1 474.(JIANG Qinghui,DENG Shushen,ZHOU Chuangbing. Study of three-dimensional high-order numerical manifold method[J]. Rock and Soil Mechanics,2006,27(9):1 471– 1 474.(in Chinese))
[27] ZHENG H,LIU Z J,GE X R. Numerical manifold space of Hermitian form and application to Kirchhoff?s thin plate problems[J]. International Journal for Numerical Methods in Engineering,2013,95(9):721–739.
[28] 王水林,葛修润. 流形元方法在模拟裂纹扩展中的应用[J]. 岩石力学与工程学报,1997,16(5):7–12.(WANG Shuilin,GE Xiurun. Application of manifold method in simulating crack propagation[J]. Chinese Journal of Rock Mechanics and Engineering,1997,16(5):7–12.(in Chinese))
[29] AN X M. Extended numerical manifold method for engineering failure analysis[Ph. D. Thesis][D]. Singapore:Nanyang Technology University Singapore,2010.
[30] 张慧华,祝晶晶. 复杂裂纹问题的多边形数值流形方法求解[J]. 固体力学学报,2013,34(1):38–46.(ZHANG Huihua,ZHU Jingjing. Numerical manifold analysis of complex crack problems on polygonal elements[J]. Chinese Journal of Solid Mechanics,2013,34(1):38–46.(in Chinese))
[31] 苏海东,祁勇峰,龚亚琦. 裂纹尖端解析解与周边数值解联合求解应力强度因子[J]. 长江科学院院报,2013,30(6):83–89.(SU Haidong,QI Yongfeng,GONG Yaqi. Compute stress intensity factors via combining analytical solutions around crack tips with surrounding numerical solutions[J]. Journal of Yangtze River Scientific Research Institute,2013,30(6):83–89.(in Chinese))
[32] LI S C,CHENG Y M. Enriched meshless manifold method for two-dimensional crack modeling[J]. Theoretical and Applied Fracture Mechanics,2005,44(3):234–248.
[33] GAO H,CHENG Y. A complex variable meshless manifold method for fracture problems[J]. International Journal of Computational Methods,2010,7(1):55–81.
[34] ZHU H,ZHUANG X,CAI Y,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.
[35] ZHANG G X,ZHU B F,LU Z C. Cracking simulation of the Wuqiangxi ship lock by manifold method[C]// Proceedings of the Sixth Internatonsl Conference on Analysis of Discontinuous Deformation (ICADD–6). [S. l.]:[s. n.],2003:133–140.
[36] WU Z,WONG L N Y. Frictional crack initiation and propagation analysis using the numerical manifold method[J]. Computers and Geotechnics,2012,39:38–53.
[37] ZHENG H,XU D D. New strategies for some issues of numerical manifold method in simulation of crack propagation[J]. International Journal for Numerical Methods in Engineering,2014,97(13):986–1 010.
[38] ZHENG H,LIU F,LI C G. The MLS-based numerical manifold method with applications to crack analysis[J]. International Journal of Fracture,2014,190(2):47–166.
[39] ZHENG H,LIU F,LI C G. Primal mixed solution to unconfined seepage flow in porous media with numerical manifold method[J]. Applied Mathematical Modelling,2015,39(2):794–808.
[40] 杨永涛,徐栋栋,郑 宏. 动载下裂纹应力强度因子计算的数值流形元法[J]. 力学学报,2014,46(5):730–738.(YANG Yongtao,XU Dongdong,ZHENG Hong. Evaluation on stress intensity factor of crack under dynamic load using numerical manifold method[J]. Chinese Journal of Theoretical and Applied Mechanics,2014,46(5):730–738.(in Chinese))
[41] 徐栋栋,郑 宏,杨永涛. 线性无关高阶数值流形法[J]. 岩土工程学报,2014,36(3):482–488.(XU Dongdong,ZHENG Hong,YANG Yongtao. Linearly independent higher-order numerical manifold method[J]. Chinese Journal of Geotechnical Engineering,2014,36(3):482–488.(in Chinese))
[42] 郭朝旭,郑 宏. 高阶数值流形方法中的线性相关问题研究[J]. 工程力学,2012,29(12):228–232.(GUO Chaoxu,ZHENG Hong. Study on linear dependence problem in high-order numerical manifold method[J]. Engineering Mechanics,2012,29(12):228–232.(in Chinese))
[43] 徐栋栋,郑 宏. 数值流形法在处理强奇异性问题时的网格无关性[J]. 岩土力学,2014,35(8):2 385–2 394.(XU Dongdong,ZHENG Hong. Mesh independence of numerical manifold method in treating strong singularity[J]. Rock and Soil Mechanics,2014,35(8):2 385–2 394.(in Chinese))
[44] 徐栋栋,郑 宏,夏开文,等. 高阶扩展数值流形法在裂纹扩展中的应用[J]. 岩石力学与工程学报,2014,33(7):1 375–1 387.(XU Dongdong,ZHENG Hong,XIA Kaiwen,et al. Application of higher-order enriched numerical manifold method to crack propagation[J]. Chinese Journal of Rock Mechanics and Engineering,2014,33(7):1 375–1 387.(in Chinese))
[45] XU D D,ZHENG H,XIA K W. New Numerical Quadrature of Integrand with Singularity of 1/r and its Application[J]. Applied Mechanics and Materials,2014,444:641–649.
[46] XU D D,ZHENG H. Mesh independence test of numerical manifold method in treating strong singularity[C]// Proceedings of 11th International Conference on Analysis of Discontinuous Deformation. Japan,Fukuoka:[s. n.],2013,91–96.
[47] ZHANG H H,LI L X,AN X M,et al. Numerical analysis of 2D crack propagation problems using the numerical manifold method[J]. Engineering Analysis with Boundary Elements,2010,34(1):41–50.
[48] 中国航空研究院主编. 应力强度因子手册(增订版)[M]. 北京:科学出版社,1993:356–357.(Chinese Aeronautical Establishment. Handbook of the stress intensity factor(revised edition)[M]. Beijing:Science Press,1993:356–357.(in Chinese))
[49] SIH G C. Energy-density concept in fracture mechanics[J]. Engineering Fracture Mechanics,1973,5(4):1 037–1 040 |
|
|
|