A generalized cover renewal strategy for multiple crack propagation in two-dimensional numerical manifold method

Abstract

Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system, i.e., the mathematical cover and physical cover. However, renewal of the topology of the two-cover system poses a challenge for multiple crack propagation problems and there are few references. In this study, a robust and efficient strategy is proposed to update the cover system of the numerical manifold method in simulation of multiple crack propagation problems. The proposed algorithm updates the cover system with a bottom-up process: 1) identification of fractured manifold elements according to the previous and latest crack tip position; and 2) local topological update of the manifold elements, physical patches, block boundary loops, and non-persistent joint loops according to the scenario classification of the propagating crack. The proposed crack tracking strategy and classification of the renewal cases promote a robust and efficient cover renewal algorithm for multiple crack propagation analysis. Three crack propagation examples show that the proposed algorithm performs well in updating the cover system. This cover renewal methodology can be extended for numerical manifold method with polygonal mathematical covers.

Keywords

numerical manifold method / multiple crack propagation / physical cover renewal / polygonal mathematical cover

References

[1] ZhengF, ZhuangX-Y, ZhengH, JiaoY-Y, RabczukT. Kinetic analysis of polyhedral block system using an improved potential-based penalty function approach for explicit discontinuous deformation analysis. Applied Mathematical Modelling, 2020, 82: 314-335

[2] ZhengF, LeungA, ZhuJ-B, JiaoY-Y. Modified predictor-corrector solution approach for efficient discontinuous deformation analysis of jointed rock masses. International Journal for Numerical and Analytical Methods in Geomechanics, 2019, 43(2): 599-624

[3] ScheldtTComparison of continuous and discontinuous modelling for computational rock mechanics, 2003, Norway, Department of Geology and Mineral Rersources Engineering, Noraegion University of Science and Technology

[4] JingL. A review of techniques, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering. International Journal of Rock Mechanics and Mining Sciences, 2003, 40(3): 283-353

[5] WipplerK, KunaM. Crack analyses in three-dimensional piezoelectric structures by the BEM. Computational Materials Science, 2007, 39(1): 261-266

[6] CordeiroS G F, LeonelE D. Cohesive crack propagation modelling in wood structures using BEM and the tangent operator technique. Engineering Analysis with Boundary Elements, 2016, 64: 111-121

[7] NguyenV P, RabczukT, BordasS, DuflotM. Meshless methods: A review and computer implementation aspects. Mathematics and Computers in Simulation, 2008, 79(3): 763-813

[8] LuY Y, BelytschkoT, GuL. A new implementation of the element free Galerkin method. Computer Methods in Applied Mechanics and Engineering, 1994, 113(3): 397-414

[9] CaiY-C, HanL, TianL-G, ZhangL-Y. Meshless method based on Shepard function and partition of unity for two-dimensional crack problems. Engineering Analysis with Boundary Elements, 2016, 65: 126-135

[10] FRIES T P, MATTHIES H. Classification and overview of meshfree methods [R]. Technical Report, TU Brunswick, Germany Nr. 2003-03.

[11] BouchardP O, BayF, ChastelY, TovenaT. Crack propagation modelling using an advanced remeshing technique. Computer Methods in Applied Mechanics and Engineering, 2000, 189(3): 723-742

[12] BouchardP O, BayF, ChastelY. Numerical modelling of crack propagation: Automatic remeshing and comparison of different criteria. Computer Methods in Applied Mechanics and Engineering, 2003, 192(35): 3887-390836

[13] LiuG R, DaiK Y, NguyenT T. A smoothed finite element method for mechanics problems. Computational Mechanics, 2006, 39(6): 859-877

[14] LiuG R, NguyenT T, DaiK Y, LamK Y. Theoretical aspects of the smoothed finite element method (SFEM). Int J Numer Methods Eng, 2007, 71(8): 902-930

[15] OliverJ. Modelling strong discontinuities in solid mechanics via strain softening constitutive equations. Part 1: Fundamentals. Int J Numer Methods Eng, 1996, 39(21): 3575-3600

[16] RabczukT, StéphaneB, ZiG. On three-dimensional modelling of crack growth using partition of unity methods. Computers & Structures, 2010, 88(23): 1391-141124

[17] DuarteC A, HamzehO N, LiszkaT J, TworzydloW W. A generalized finite element method for the simulation of three-dimensional dynamic crack propagation. Computer Methods in Applied Mechanics and Engineering, 2001, 190(15–17): 2227-2262

[18] SalimzadehS, KhaliliN. A three-phase XFEM model for hydraulic fracturing with cohesive crack propagation. Computers and Geotechnics, 2015, 69: 82-92

[19] XieY-S, CaoP, LiuJ, DongL-W. Influence of crack surface friction on crack initiation and propagation: A numerical investigation based on extended finite element method. Computers and Geotechnics, 2016, 74: 1-14

[20] BERGARA A, DORADO J I, MARTíN-MEIZOSO A, MARTINEZ-ESNADLA J M. Fatigue crack propagation in complex stress fields: experiments and numerical simulations using the extended finite element method (Xfem) [J]. International Journal of Fatigue, 2017. DOI: https://doi.org/10.1016/j.ijfatigue.2017.05.026.

[21] BabuskaI, CalozG, OsbornJ E. Special finite element methods for a class of second order elliptic problems with rough coefficients. SIAM Journal on Numerical Analysis, 1994, 31(4): 945-981

[22] VesgaL, VallejoL, Lobo-GuerreroS. DEM analysis of the crack propagation in brittle clays under uniaxial compression tests. International Journal for Numerical and Analytical Methods in Geomechanics, 2008, 32: 1405-1415

[23] YangD-M, ShengY, YeJ-Q, TanY-Q. Dynamic simulation of crack initiation and propagation in cross-ply laminates by DEM. Composites Science and Technology, 2011, 71(11): 1410-1418

[24] JiaoY-Y, ZhangX-L, ZhaoJ. Two-dimensional DDA contact constitutive model for simulating rock fragmentation. Journal of Engineering Mechanics, 2012, 138: 199-209

[25] ZhangX-L, JiaoY-Y, ZhaoJ. Simulation of failure process of jointed rock. Journal of Central South University of Technology, 2008, 15: 888-894

[26] ShiG-HDiscontinuous deformation analysis: A new numerical model for the statics and dynamics of block systems, 1988, Berkeley, University of California

[27] NingY-J, YangJ, AnX-M, MaG-W. Modelling rock fracturing and blast-induced rock mass failure via advanced discretisation within the discontinuous deformation analysis framework. Computers and Geotechnics, 2011, 38(1): 40-49

[28] TsayR J, ChiouY J, ChuangW L. Crack growth prediction by manifold method. Journal of Engineering Mechanics, 1999, 125: 884

[29] ZhangH H, LiL X, AnX M, MaG W. Numerical analysis of 2-D crack propagation problems using the numerical manifold method. Engineering Analysis with Boundary Elements, 2010, 34(1): 41-50

[30] WuZ-J, WongL N Y. Frictional crack initiation and propagation analysis using the numerical manifold method. Computers and Geotechnics, 2012, 3938-53

[31] HeJ, LiuQ-S, MaG-W, ZengB. An improved numerical manifold method incorporating hybrid crack element for crack propagation simulation. International Journal of Fracture, 2016, 199(1): 21-38

[32] KourepinisD, PearceC, BicanicN. Higher-order discontinuous modeling of fracturing in concrete using the numerical manifold method. International Journal of Computational Methods, 2010, 7(1): 83-106

[33] ZhaoG-F, MaG-W, ZhangH-H, ZhaoJ. A numerical manifold method for plane micropolar elasticity. International Journal of Computational Methods, 2010, 7(1): 151-166

[34] LiuF, ZhengH, DuX-L. Hybrid analytical and MLS-based NMM for the determination of generalized stress intensity factors. Mathematical Problems in Engineering, 2015, 2015: 1-9

[35] ZhengH, XuD-D. New strategies for some issues of numerical manifold method in simulation of crack propagation. Int J Numer Methods Eng, 2014, 97(13): 986-1010

[36] ZhengH, LiuF, DuX-L. Complementarity problem arising from static growth of multiple cracks and MLS-based numerical manifold method. Computer Methods in Applied Mechanics and Engineering, 2015, 295: 150-171

[37] YangS-K, MaG-W, RenX-H, RenF. Cover refinement of numerical manifold method for crack propagation simulation. Engineering Analysis with Boundary Elements, 2014, 43: 37-49

[38] YangY-T, ZhengH, SivaselvanM V. A rigorous and unified mass lumping scheme for higher-order elements. Computer Methods in Applied Mechanics & Engineering, 2017, 319(1): 491-514

[39] CaiY-C, ZhuangX-Y, ZhuH-H. A generalized and efficient method for finite cover generation in the numerical manifold method. International Journal of Computational Methods, 2013, 10(5): 1350028

[40] ZhangH H, MaG W. Fracture modeling of isotropic functionally graded materials by the numerical manifold method. Engineering Analysis with Boundary Elements, 2014, 38: 61-71

[41] ZhangH-H, YanJ-X. Investigation of the accuracy of the numerical manifold method on n-sided regular elements for crack problems. Applied Mechanics and Materials, 2012, 157–1581093-1096

[42] MaG W, AnX M, ZhangH H, LiL X. Modeling complex crack problems using the numerical manifold method. International Journal of Fracture, 2009, 156(1): 21-35

[43] ChenG Q, OhnishiY, ItoT. Development of high-order manifold method. Int J Numer Methods Eng, 1998, 43(4): 685-712

[44] LinJ S, KuC Y. Two-scale modeling of jointed rock masses. International Journal of Rock Mechanics and Mining Sciences, 2006, 43(3): 426-436

[45] MikiS, SasakiT, KoyamaT, NishiyamaS, OhnishiY ODevelopment of coupled discontinuous deformation analysis and numerical manifold method (NMM-DDA) and its application to dynamic problems, 2010, Singapore, Research Publishing Services

[46] ZhangG X, ZhaoY, PengX C. Simulation of toppling failure of rock slope by numerical manifold method. International Journal of Computational Methods, 2010, 7(1): 167-189

[47] CHEN G Q, JIANG Z S, WU Y Q. A new approach for numerical manifold method [J]. IEIT Journal of Adaptive & Dynamic Computing, 2012: 23–34. DOI: https://doi.org/10.5813/www.ieit-web.org/IJADC/2012.2.5.

[48] ZhangH H, ChenY L, LiL X, et al.. Accuracy comparison of rectangular and triangular mathematical elements in the numerical manifold method. Analysis of Discontinuous of Deformation New Developments and Applications, 2010, Singapore, Research Publishing Services, 297303Vol.1

[49] ZHANG Z, ZHANG X. Direct simulation of low-re flow around a square cylinder by numerical manifold method for navier-stokes equations [J]. Journal of Applied Mathematics, 2012: 465972. DOI: https://doi.org/10.1155/2012/465972.

[50] ZhangZ R, ZhangX W, LuW G. Numerical method based on compatible manifold element for thin plate bending. Chin J Mech Eng, 2010, 23(1): 100-109

[51] WuZ-J, WongL N Y. Modeling cracking behavior of rock mass containing inclusions using the enriched numerical manifold method. Engineering Geology, 2013, 162: 1-13

[52] WuZ-J, WongL N Y. Elastic-plastic cracking analysis for brittle-ductile rocks using manifold method. International Journal of Fracture, 2013, 180(1): 71-91

[53] WuZ-J, WongL N Y, FanL-F. Dynamic study on fracture problems in viscoelastic sedimentary rocks using the numerical manifold method. Rock Mechanics and Rock Engineering, 2013, 46(6): 1415-1427

[54] NingY J, AnX M, MaG W. Footwall slope stability analysis with the numerical manifold method. International Journal of Rock Mechanics and Mining Sciences, 2011, 48(6): 964-975

[55] AnX-M, NingY-J, MaG-W, HeL. Modeling progressive failures in rock slopes with non-persistent joints using the numerical manifold method. International Journal for Numerical and Analytical Methods in Geomechanics, 2014, 38(7): 679-701

[56] LimI L, JohnstonI W, ChoiS K. Stress intensity factors for semi-circular specimens under three-point bending. Engineering Fracture Mechanics, 1993, 44(3): 363-382

[57] LimI L, JohnstonI W, ChoiS K, BolandJ N. Fracture testing of a soft rock with semi-circular specimens under three-point bending. Part 1—mode I. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1994, 31(3): 185-197

[58] LimI L, JohnstonI W, ChoiS K, BolandJ N. Fracture testing of a soft rock with semi-circular specimens under three-point bending. Part 2—mixed-mode. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1994, 31(3): 199-212

[59] AdamsonR M, DempseyJ P, MulmuleS V. Fracture analysis of semi-circular and semi-circular-bend geometries. International Journal of Fracture, 1996, 77(3): 213-222

[60] AyatollahiM R, AlihaM R M. On determination of mode II fracture toughness using semi-circular bend specimen. International Journal of Solids and Structures, 2006, 43(17): 5217-5227

[61] ChenR, XiaK, DaiF, LuF, LuoS N. Determination of dynamic fracture parameters using a semi-circular bend technique in split Hopkinson pressure bar testing. Engineering Fracture Mechanics, 2009, 76(9): 1268-1276

[62] AyatollahiM R, AlihaM R M, HassaniM M. Mixed mode brittle fracture in PMMA—An experimental study using SCB specimens. Materials Science and Engineering A, 2006, 417(1): 348-3562

[63] ZhangQ B, ZhaoJ. Determination of mechanical properties and full-field strain measurements of rock material under dynamic loads. International Journal of Rock Mechanics and Mining Sciences, 2013, 60: 423-439

[64] AlihaM R M, AyatollahiM R, SmithD J, PavierM J. Geometry and size effects on fracture trajectory in a limestone rock under mixed mode loading. Engineering Fracture Mechanics, 2010, 77(11): 2200-2212

[65] XuY, DaiF, XuN W, ZhaoT. Numerical investigation of dynamic rock fracture toughness determination using a semi-circular bend specimen in split hopkinson pressure bar testing. Rock Mechanics and Rock Engineering, 2015, 49(3): 731-745

[66] XieY, CaoP, JinJ, WangM. Mixed mode fracture analysis of semi-circular bend (SCB) specimen: A numerical study based on extended finite element method. Computers and Geotechnics, 2017, 82: 157-172

[67] BoccaP, CarpinteriA, ValenteS. Size effects in the mixed mode crack propagation: Softening and snap-back analysis. Engineering Fracture Mechanics, 1990, 35(1): 159-170

[68] SchlangenEExperimental and numerical analyses of fracture processes in concrete, 1993, The Netherlands, Department of Civil Engineering, Delft University of Technology

[69] GeersM, BorstR, PeerlingsR H J. Damage and crack modeling in single-edge and double-edge notched concrete beams. Engineering Fracture Mechanics, 2000, 65: 247-261

[70] ZhuW-C, TangC-A. Numerical simulation on shear fracture process of concrete using mesoscopic mechanical model. Construction and Building Materials, 2002, 16: 453-463

[71] OliverJ, HuespeA, SamaniegoE, ChavesE W. Continuum approach to the numerical simulation of material failure in concrete. International Journal for Numerical and Analytical Methods in Geomechanics, 2004, 28: 609-632