On the accuracy of higher order displacement discontinuity method (HODDM) in the solution of linear elastic fracture mechanics problems

Abstract

The higher order displacement discontinuity method (HODDM) utilizing special crack tip elements has been used in the solution of linear elastic fracture mechanics (LEFM) problems. The paper has selected several example problems from the fracture mechanics literature (with available analytical solutions) including center slant crack in an infinite and finite body, single and double edge cracks, cracks emanating from a circular hole. The numerical values of Mode I and Mode II SIFs for these problems using HODDM are in excellent agreement with analytical results (reaching up to 0.001% deviation from their analytical results). The HODDM is also compared with the XFEM and a modified XFEM results. The results show that the HODDM needs a considerably lower computational effort (with less than 400 nodes) than the XFEM and the modified XFEM (which needs more than 10000 nodes) to reach a much higher accuracy. The proposed HODDM offers higher accuracy and lower computation effort for a wide range of problems in LEFM.

Keywords

accuracy / higher order displacement discontinuity method / LEFM / mode I stress intensity factor / mode II stress intensity factor

References

[1] GalvisA F, SolleroP. 2D analysis of intergranular dynamic crack propagation in polycrystalline materials a multiscale cohesive zone model and dual reciprocity boundary elements [J]. Comput Struct, 2016, 164: 1-14

[2] BenedettiI, AliabadiM H, DaviG. A fast 3D dual boundary element method based on hierarchical matrices [J]. Int J Solids Struct, 2008, 45: 922-933

[3] DI PisaC, AliabadiM H. Fatigue crack growth analysis of assembled plate structures with dual boundary element method [J]. Eng Fract Mech, 2013, 98: 200-213

[4] Thanh TuB, PopovV. Boundary element dual reciprocity method with overlapping sub-domains [J]. Bound Elem Other Mesh Reduct Methods XXXI, 2008179-187

[5] ChuangS H, YuehC Y, HuangL H. Dual boundary element model coupled with the dual reciprocity method to determine wave scattering by a concentric cylindrical system mounted on a conical shoal [J]. Eng Anal Bound Elem, 2015, 56: 30-38

[6] TanM, WangP, LiJ, LiuQ, YangQ. Numerical simulation and fracture evaluation method of dual laterolog in organic shale [J]. J Appl Geophys, 2014, 100: 1-13

[7] YunB I, AngW T. A dual-reciprocity boundary element approach for axisymmetric nonlinear time-dependent heat conduction in a nonhomogeneous solid [J]. Eng Anal Bound Elem, 2010, 34: 697-706

[8] Kamali YazdiA, OmidvarB, RahimianM. Improving the stability of time domain dual boundary element method for three dimensional fracture problems: A time weighting approach [J]. Eng Anal Bound Elem, 2011, 35: 1142-1148

[9] RomlayF, RM, OuyangH, AriffinA K, MohamedN, AN. Modeling of fatigue crack propagation using dual boundary element method and Gaussian Monte Carlo method [J]. Eng Anal Bound Elem, 2010, 34: 297-305

[10] NatarajanS, MahapatraD R, BordasS, PA. Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework [J]. Int J Numer Meth Eng, 2010, 83: 269-294

[11] SimpsonR, TrevelyanJ. A partition of unity enriched dual boundary element method for accurate computations in fracture mechanics [J]. Comput Methods Appl Mech Eng, 2011, 200: 1-10

[12] CrouchS L. Solution of plane elasticity problems by the displacement discontinuity method I. Infinite body solution [J]. Int J Numer Methods Eng, 1976, 10: 301-343

[13] BehniaM, GoshtasbiK, MarjiM F, GolshaniA. The effect of layers elastic parameters on hydraulic fracturing propagation utilizing displacement discontinuity method [J]. Anal Numer Methods Min Eng, 2012, 3: 1-14

[14] HAERI H, SHAHRIAR K, FATEHI MARJI M, MOAREF VAND P. Simulating the bluntness of TBM disc cutters in rocks using displacement discontinuity method [C// 13th Int Conf Fract. Beijing: International Congresson Fracture, 2013.

[15] Fatehi MarjiM, GholamnejadJ, EghbalM. on the crack propagation mechanism of brittle substances under various loading conditions [C]. 11th Int Multidiscip Sci Geo-Conference. Albena, Bulgaria: SGEM, 2011

[16] HaeriH, ShahriarK, Fatehi MarjiM, MoarefvandP. An experimental and numerical study of crack propagation and cracks coalescence in the pre-cracked rock-like disc specimens under compression [J]. Int J Rock Mech Min Sci Geomech Abstr, 2013, 67: 20-28

[17] AbdollahipourA, Fatehi MarjiM, YarahmadibafghiA, GholamnejadJ. Simulating the propagation of hydraulic fractures from a circular wellbore using the displacement discontinuity method [J]. Int J Rock Mech Min Sci, 2015, 80: 281-291

[18] CourtesenD L. Cavities and gas penetration from blasts in stressed rock with flooded joints [J]. Acta Astron, 1979, 6: 341-363

[19] HaeriH, ShahriarK, Fatehi MarjiM, MoarefvandP. A boundary element analysis of the crack propagation mechanism of random micro cracks in rock-like specimens under uniform normal tension [J]. J Min Environ, 2014, 6: 73-93

[20] CrawfordA M, CurranJ H. Higher order functional varriation displacement discontinuity elements [J]. Int J Rock Mech Min Sci Geomech Abstr, 1982, 19: 143-148

[21] Fatehi MarjiMHigher order displacement discontinuity method in rock fracture mechanics [D], 2015YazdYazd University

[22] AbdollahipourACrack propagation mechanism in hydraulic fracturing procedure in oil reservoirs [D], 2015YazdYazd University

[23] CrouchS L, StarfieldA MBoundary element methods in solid mechanics [M], 1983LondonGeorge allen & Unwin

[24] Fatehi MarjiMModelling of cracks in rock fragmentation with a higher order displacement discontinuity method [D], 1997AnkaraMiddle East Technical University

[25] Fatehi MarjiM, Hosseini-NasabH. Application of higher order displacement discontinuity method using special crack tip elements in rock fracture mechanics [C]. 20th World Min Congr Expo. Tehran: Tehran university, 2005699-704

[26] ItouS. Effect of couple-stresses on the Mode I dynamic stress intensity factors for two equal collinear cracks in an infinite elastic medium during passage of time-harmonic stress waves [J]. Int J Solids Struct, 2013, 50: 1597-1604

[27] AbdollahipourA, Fatehi MarjiM, YarahmadibafghiA, GholamnejadJ. Numerical investigation on the effect of crack geometrical parameters in hydraulic fracturing process of hydrocarbon reservoirs [J]. J Min Environ, 2016, 7: 205-214

[28] Fatehi MarjiM, Hosseini NasabH, KohsaryA H. On the uses of special crack tip elements in numerical rock fracture mechanics [J]. Int J Solids Struct, 2006, 43: 1669-1692

[29] Fatehi MarjiM, Hosseini-NasabH, MorshediA H. Numerical modeling of the mechanism of crack propagation in rocks under TBM disc cutters [J]. Mech Mater Struct, 2009, 4: 605-627

[30] SanfordR JPrinciples of fracture mechanics [M], 2003USAPrentice Hall

[31] ShouK J, CrouchS L. A higher order displacement discontinuity method for analysis of crack problems [J]. Int J Rock Mech Min Sci Geomech Abstr, 1995, 32: 49-55

[32] GuoH, AzizN I, SchmittL C. Linear elastic crack tip modeling by displacement discontinuity method [J]. Engin Fract Mech, 1990, 36: 933-943

[33] WhittakerB N, SinghR N, SunGRock fracture mechanics, principles design and applications [M], 1992NetherlandElsevier Science Ltd

[34] TadaH, ParisP C, IrwinG RThe stress analysis of cracks handbook [M], 2000New YorkASME

[35] MurakamiYStress intensity factors handbook [M], 1987EnglandPergamon Press