Mesh-independent equivalent domain integral method for J-integral evaluation. (October 2016)

• Record Type: Journal Article
• Authors: Nikishkov, G.P.; Vershinin, A.V.; Nikishkov, Y.G.
• Is Part Of: Advances in engineering software. Volume 100(2016)
• Page Start: 308
• Page End: 318
• Publication Date: 2016-10
• Subjects: Equivalent domain integral method -- J-Integral -- Stress intensity factors -- Moving least squares -- Finite element method -- Hexahedral elements -- Tetrahedral elements; Computer-aided engineering -- Periodicals; Engineering -- Computer programs -- Periodicals; Engineering -- Software -- Periodicals; Periodicals; 620.0028553
• DOI: 10.1016/j.advengsoft.2016.08.006
• Languages: English
• ISSNs: 0965-9978

Abstract

Highlights: Further development of the domain integral method for J-integral calculation in 3D. Moving least squares approximation provides independence of finite element mesh. Domain integration in global coordinates. Using vector weight function to obtain the J -integral in crack front coordinates. Special integration rule with double coordinate change employed in J 2 calculation. Abstract: The equivalent domain integral method is a reliable tool for J -integral computation in two- and three-dimensional elastic and elastic-plastic fracture mechanics problems. A variant of this method that is independent of finite element mesh is presented. Finite element solution of a boundary value problem is performed on a mesh composed of arbitrary elements. Nodal results are approximated by the moving least squares method that does not require knowledge of mesh topology. Domain integrals are evaluated on a background mesh of hexahedral elements. The mesh has the polar structure with the refinement towards the crack front. Elements of the background mesh are generated in the coordinate system associated with the crack front and then transformed to the global system. Domain integration for each background element is performed once during computations. Evaluation of J -integral for multiple domains is achieved by multiplication of an element domain integral with multiple domain weight functions. Performance of the proposed algorithm is demonstrated by the examples of three-dimensional cracks using meshes of both hexahedral and tetrahedral elements.