Numerical methods such as boundary element methods are widely used for the stress analysis in solid mechanics. These methods are also used for crack analysis in rock fracture mechanics. There are singularities for the stresses and displacements at the crack tips in fracture mechanics problem, which decrease the accuracy of the numerical results in areas very close to the crack ends. To overcome this, higher order elements and isoperimetric higher order elements have been used. Recently, special crack tip elements have been proposed and used in most of the numerical fracture mechanics models. These elements can drastically increase the accuracy of the results near the crack tips, but in most of the models only one special crack tip element has been used for each crack end. In this study the uses of higher order crack tip elements are discussed and a higher order displacement discontinuity method is used to investigate the effect of these elements on the accuracy of the results in some crack problems. The useful shape functions for two special crack tip elements, are derived and given in the text and appendix for both infinite and semi-infinite plane problems. In this analysis both Mode I and Mode II stress intensity factors are computed . Some example problems are solved and the computed results are compared with the results given in the literature. The numerical results obtained here are in good agreement with those cited in the literature. For the curved crack problem, the strain energy release rate, G can be calculated accurately in the vicinity of the crack tips by using the higher order displacement discontinuity method with a quadratic variation of displacement discontinuity elements and with two special crack tip elements at each crack end.
This paper describes progress on the development of theoretical models required for studying failure mechanism, crack initiation and growth around the boreholes driven by hydrofracturing processes in Hot Dry Rock (HDR) reservoirs of geothermal energy. Due to the importance of the stress intensity factor concept (K) in Fracture Mechanics, some advanced modeling techniques for accurate and fast determination of K for relevant problems are proposed. Alternative tools to deal with stress intensity factor determination are developed and assessed from the points of view of accuracy and computational cost.We concentrate on residual strength, crack initiation and crack growth as a means to model and understand experimentally observed behaviors. Several modeling methods such as compounding and weight function techniques, and boundary and finite element modeling for stress intensity factor calculation are discussed. Further to reviews of those techniques, work performed included (i) developing alternative solutions to deal with boundary-to-boundary interaction when using the compounding technique, (ii) relating the precision of K calculations with the level of precision of the crack opening displacement of a reference solution, in order to assess the precision of weight function technique, (iii) modeling relevant geometries using the finite element method (FEM), (iv) working on the implementation of direct stress intensity factor K determination in the Higher Order Displacement Discontinuity Method (HODDM), and (v) developing tools to deal with residual stress fields around the boundary of the hydraulically pressurized boreholes.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.