RESEARCH ARTICLE

A family of non-conforming crack front elements of quadrilateral and triangular types for 3D crack problems using the boundary element method

  • Guizhong XIE 1 ,
  • Fenglin ZHOU , 2 ,
  • Hao LI 1 ,
  • Xiaoyu WEN 1 ,
  • Fannian MENG 1
Expand
  • 1. Henan Provincial Key Laboratory of Intelligent Manufacturing of Mechanical, Mechanical and Electrical Engineering Institute, Zhengzhou University of Light Industry, Zhengzhou 450002, China
  • 2. College of Mechanical Engineering, Hunan University of Technology, Zhuzhou 412007, China

Received date: 03 Feb 2019

Accepted date: 30 Mar 2019

Published date: 15 Sep 2019

Copyright

2019 Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature

Abstract

This study focuses on establishing non- conforming crack front elements of quadrilateral and triangular types for 3D crack problems when the dual boundary element method is applied. The asymptotic behavior of the physical variables in the area near the crack front is fully considered in the construction of the shape function. In the developed quadrilateral and triangular crack front elements, the asymptotic term, which captures the asymptotic behavior of the physical variable, is multiplied directly by the conventional Lagrange shape function to form a new crack front shape function. Several benchmark numerical examples that consider penny-shaped cracks and straight-edge crack problems are presented to illustrate the validity and efficiency of the developed crack front elements.

Cite this article

Guizhong XIE , Fenglin ZHOU , Hao LI , Xiaoyu WEN , Fannian MENG . A family of non-conforming crack front elements of quadrilateral and triangular types for 3D crack problems using the boundary element method[J]. Frontiers of Mechanical Engineering, 2019 , 14(3) : 332 -341 . DOI: 10.1007/s11465-019-0540-3

Introduction

Experimental [14] and numerical [59] are two of the most popular methods to solve engineering problems, such as solid mechanical, heat transfer, acoustic, and electro-magnetic problems. In structural mechanical problems with cracks or V-notch interfaces, physical quantities vary singularly in the areas near the crack front and interface [59]. Accurate approximation of such singular variations is important but difficult in numerical methods for these problems. Many additional methods have been developed to circumvent these difficulties [59]. The finite element method (FEM) [1012] and the boundary element method (BEM) [1315] are two of the most extensive methods for structural mechanical problems. Special elements with specially defined shape functions are typically utilized to approximate the variations of physical quantities in both methods and eventually achieve high accuracy and efficiency. Many specially defined crack elements have been developed in FEM and BEM implementations to compute the stress intensity factors (SIFs) along the crack front accurately and efficiently. The most efficient and valuable crack front elements are quarter-point triangular and quadrilateral elements, which can be integrated into standard FEM [16,17] and BEM [1821] codes. In BEM application in crack problems, the dual boundary integral equation method, in which the traction boundary integral equation is involved, is the most widely used scheme [5,1924].
Non-conforming elements, inside which interpolation nodes are arranged, are frequently used due to the existence of hyper-singular integrals. The calculation of hyper-singular integrals depends largely on the transformation of traditional conforming quarter-point elements into nonconforming ones. However, quarter-point triangular and quadrilateral elements are usually regarded as conforming elements, which are employed in the multidomain BEM. The shape functions of non-conforming crack front elements are obtained by solving linear equations. The derivation process is tedious, and the special shape functions are not closely related with traditional shape functions [20,25,26]. This study proposes a family of non-conforming crack front elements of quadrilateral and triangular types for 3D crack problems.
In the deduction of the presented crack front elements, the shape functions are constructed to capture the asymptotic behavior of displacement fields in parametric space. The function of the analytical distribution of displacement is expressed by Taylor series up to the second order to capture the asymptotic behavior. Variable transformation is performed to generate a square-root term that captures the variation of the physical quantities. The shape functions are constructed in consideration of the square-root term. Then, the conforming element is transformed into a nonconforming one by a multiplication of the transforming matrix. Several benchmark examples are presented to illustrate the accuracy of the new crack front elements. The advantages of the non-conforming crack front elements are as follows. The special shape functions are obtained using the traditional shape functions and displacement approximate terms in the intrinsic coordinate system. Unlike the traditional construction method for crack front elements [20,25,26], the method proposed in this study can be easily extended to the construction of special shape functions for V-notch interfaces, whose displacements may possess a property of order rα (1>α>0).
The rest of this paper is outlined as follows. In Section 2, crack front elements within a quadrilateral mesh are proposed. Then, crack front elements within a triangular mesh are obtained in the same manner in Section 3. The process of obtaining SIFs is described in Section 4, and several numerical examples are presented in Section 5.

Crack front elements within a quadrilateral mesh

We provide a detailed derivation of a three-node, non-conforming crack tip element. The method for the three-node, non-conforming crack tip element is then extended to the 3D boundary surface element.

Construction of a three-node, non-conforming crack tip element

As shown in Fig. 1, the coordinates of the distance vector from arbitrary point P(x, y) in the element to the crack tip can be expressed as
x x0=N0x0+N1x1+ N2 x2x0= 12( ξ+1)[(ξ 2) x0+2(1ξ)x2+ξx1] ,
y y0=N0y0+N1y1+ N2 y2y0= 12( ξ+1)[(ξ 2) y0+2(1ξ)y2+ξy1] ,
where ξ is the parametric coordinate of the point (x,y), (x,y ) is the coordinates of any point in the element, (x0,y0), (x 1,y1), and (x 2,y2) are the coordinates of the three nodes in Fig. 1 (Points 0, 1 and 2), respectively, N0(ξ )= 12ξ (ξ1), N1( ξ)=12ξ(1+ξ) and N2( ξ)=(1+ξ)(1 ξ) are shape functions that are defined on these nodes, and (ξ+1)/2 represents the distance property in the parametric space, which also reflects the distance in the Cartesian space. The source points in the parametric space are λ, 0, and λ.
Fig.1 Three-node, non-conforming quadratic element

Full size|PPT slide

The displacement u and traction fields t near the crack tip behave as square root singularity or inverse square root singularity, which can be written as [27]
u=a0+a1r+ a 2r+,
t=b0/r+b1+b2r+.
where ai and bi (i=0, 1, …) represent the coefficient vectors, respectively.
For the three-node, non-conforming quadratic element, the shape functions that are defined based on these collocation nodes can be written as
Ncoll0(ξ )= ξ(ξ λ)2λ2 ,Ncoll1(ξ)=( ξ+λ)ξ2 λ2,
Ncoll2(ξ )=(ξ +λ)(ξλ)λ2.
where λ is the offset parameter, and Ncoll0(ξ), Ncoll1(ξ), and Ncoll2(ξ) are the non-conforming shape functions of the three collocation nodes in Fig. 1 (Points 0, 1 and 2). Using a two-order Taylor expansion of u near the crack front, we can obtain the following:
u=u0 0+u0 1(ξ+1)+ 12! u 02(ξ+1)2+,
where
u00= m =02 d (0)N collm(ξ)d ξ(0)ξ=1um
=( 12λ + 12λ2)u0+ (12 λ2 12λ) u 1+(1 1 λ2)u2,
u01=m=02dNcollm(ξ) dξξ=1um
=( 1 λ2 12λ) u 0+( 1 2λ1λ 2)u1 2λ2 u2,
u02= m =02 d 2 Ncollm(ξ)dξ 2ξ=1um= 1 λ2u0+ 1λ2 u1+ 2 λ2u2,
and u0, u1, and u2 are the displacements at three collocation nodes. u00, u01, and u02 represent the zero- first- and second-order derivatives at Node 0.
Compared with Eq. (2), Eqs. (4) and (5) show that the shape functions of the non-conforming quadratic element cannot reflect the square-root asymptotic behavior of the displacements and tractions near the crack front. To introduce the asymptotic behavior of the displacements, a substitution of ( ξ+1)/2 by (ξ+1)/2 is introduced in the Taylor expansion. Thus, the Taylor expansion in Eq. (5) can be transformed into
u=u00+u01 2(ξ+1)+u02(ξ+1)+.
In Eq. (6), u00, u01, and u02 are unchanged. After substituting expressions for u00, u01, and u02 into Eq. (6), the following equation can be obtained:
u =(12λ+ 12λ2) u 0+ ( 1 2λ212λ)u1+( 1 1λ2)u2 + 2 ( (1λ 2+ 1 2λ)u0+ (12λ 1λ2)u1 2 λ2u2 ) ξ+1+( 1 λ2 u0+ 1 λ2u1+ 2λ2 u2)(ξ+ 1) =[ ( 1 2λ+12 λ2) 2 (1λ 2+ 1 2λ)ξ+1+1λ 2(ξ+1 )] u 0 +[ ( 1 2λ212λ)+2( 1 2λ1λ 2)ξ+1+1λ 2(ξ+1 )] u 1 +[ (1 1 λ2) 22λ2 ξ+1+ 2 λ2(ξ +1)] u2.
New shape functions for collocation points can be defined as
u =(12λ+ 12λ2) u 0+ ( 1 2λ212λ)u1+( 1 1λ2)u2 + 2 ( (1λ 2+ 1 2λ)u0+ (12λ 1λ2)u1 2 λ2u2 ) × ξ+1+ ( 1 λ2u0+ 1λ2 u1+ 2 λ2u2 )( ξ+1) =[( 12λ + 12λ2)2( 1 λ2+ 1 2λ)ξ+1+1λ 2(ξ+1 )] u 0 +[ ( 1 2λ212λ)+2( 1 2λ1λ 2)ξ+1+1λ 2(ξ+1 )] u 1 +[ (1 1 λ2) 22λ2 ξ+1+ 2 λ2(ξ +1)] u2.
where N scoll0(ξ ), Nscoll1(ξ), and Nscoll2(ξ) are the new shape functions at collocation points for crack tip elements.
A three-node non-conforming crack tip element can be successfully constructed with the shape functions of Eq. (8). In the next section, we extend this crack tip element to a nine-node quadrilateral crack front element.

Construction of a nine-node non-conforming crack front element

A nine-node non-conforming quadratic element is illustrated in Fig. 2. The relations between the three-node and nine-node non-conforming quadratic elements can be established as follows. In the nine-node non-conforming quadratic element, interpolation nodes do not coincide with geometric nodes, a condition that is similar to that in the three-node non-conforming quadratic element. The geometric shape functions in the nine-node non-conforming quadratic element are constructed through the products of shape functions along two different directions:
N0( ξ,η)=N1(ξ)N1( η), N1(ξ ,η)= N0( ξ)N1(η ),N2(ξ,η)=N0(ξ )N0(η),N 3(ξ ,η)= N1( ξ)N0(η ),N4(ξ,η)=N2(ξ )N1(η),N 5(ξ ,η)= N0( ξ)N2(η ),N6(ξ,η)=N2(ξ )N0(η),N 7(ξ ,η)= N1( ξ)N2(η ),N8(ξ,η)=N2(ξ )N2(η),
where Ni(ξ,η) ( i=0,1,... ,8) are the shape functions that are defined on the eight geometric nodes (from 0 to 8), N0(ξ), N1( ξ) and N2(ξ) have been defined before, and ξ and η are the parametric coordinates.
Fig.2 Nine-node non-conforming quadratic element

Full size|PPT slide

Similarly, the non-conforming shape functions of the nine-node non-conforming quadratic element can be obtained through a similar method, which can be expressed as follows:
Ncoll0(ξ ,η)= Ncoll 1(ξ)Ncoll1(η) , Ncoll1( ξ,η)=Ncoll0 (ξ)N coll1(η ),Ncoll2(ξ,η)=N coll0(ξ )Ncoll0 (η),N coll3(ξ ,η)= Ncoll 1(ξ)Ncoll0(η) , Ncoll4( ξ,η)=Ncoll2 (ξ)N coll1(η ),Ncoll5(ξ,η)=N coll0(ξ )Ncoll2 (η),N coll6(ξ ,η)= Ncoll 2(ξ)Ncoll0(η) , Ncoll7( ξ,η)=Ncoll1 (ξ)N coll2(η ),Ncoll8(ξ,η)=N coll2(ξ )Ncoll2 (η),
where Ncolli (ξ,η) ( i=0,1,... ,8) are the non-conforming shape functions that are defined on the eight collocation nodes (from 0 to 8).
In Fig. 2, given that ξ =1 is the crack front, the special shape functions can be expressed as
NNscoll0(ξ ,η)= Nscoll 1(ξ)Ncoll1(η) , NNscoll1( ξ,η)=Nscoll0 (ξ)N coll1(η ),NNscoll2(ξ,η)=N scoll0(ξ )Ncoll0 (η),N Nscoll3(ξ ,η)= Nscoll 1(ξ)Ncoll0(η) ,
NNscoll4(ξ ,η)= Nscoll 2(ξ)Ncoll1(η) , NNscoll5( ξ,η)=Nscoll0 (ξ)N coll2(η ),NNscoll6(ξ,η)=N scoll2(ξ )Ncoll0 (η),N Nscoll7(ξ ,η)= Nscoll 1(ξ)Ncoll2(η) , NNscoll8( ξ,η)=Nscoll2 (ξ)N coll2(η ),
where N Nscolli(ξ,η ) ( i=0,1,... ,8) are the new shape functions on the eight collocation nodes for crack front elements (from 0 to 8).
The displacement distribution near the crack front can be accurately approximated using Eq. (11). The results are demonstrated in the following numerical examples.

Crack front elements within a triangular mesh

In this section, we construct a triangular crack front element based on a six-node conforming quadratic element in the first step. The six-node conforming crack front element is then converted into a six-node non-conforming crack front element in the second step. The six-node non-conforming quadratic element is shown in Fig. 3, where k is the position parameter. From the geometric meaning of the area coordinate, ( 1ξ η) represents the distance in the parametric space, which also reflects the distance in the Cartesian space. The geometric shape functions of the element are
N0= ξ(2ξ1 ),N1= η(2 η1 ),N2=(1ξη)[ 2(1ξη) 1] ,N3=4ξη,N4= 4η( 1ξ η),N5=4 ξ(1ξη ).
Fig.3 Six-node non-conforming quadratic element

Full size|PPT slide

ƒTo perform Taylor expansion of the displacement fields near the neighborhood of the crack front (ξ+ η=1), the following coordinate transformation is applied:
ξ1= 1ξ η,η1=ξη.
Then, Taylor expansion of the displacement fields in the neighborhood of ξ+η =1 is performed similar to that in the neighborhood of ξ1=0. We can derive the following expression:
u=u00+u01( ξ10)+12!u02(ξ 10)2+.
where u00= m=05 (0) Nmξ1 (0)ξ 1= 0 um, u01=m=05 Nmξ1 ξ1=0um, and u02=m=052 Nmξ12ξ 1= 0 um. um(m=0, 1,...,5) denotes the displacements at the six interpolation nodes in the element. From Eq. (13), we obtain
ξ =1ξ 1+η1 2,η= 1ξ1η1 2.
Considering Eq. (15), the following relations are obtained:
Nm(ξ)ξ 1= Nm(ξ) ξξξ1 +Nm( ξ) ηηξ1 ,
2N ξ12= 2N ξ2 ξξ1ξξ1 +22N ξη ξξ1ηξ1 + 2Nη 2 ηξ1 η ξ1.
To introduce the asymptotic behavior of the displacements, we correct the Taylor expansion expression by replacing ξ1 with ξ1. The Taylor expansion in Eq. (14) can be transformed into
u=u00+u01 ξ1+12! u 02ξ1+,
where
{ u 00=ξ(2ξ1)u0+(1ξ)(12ξ)u1+4ξ (1 ξ)u3,u01=0.5[(4ξ1 )u0 +(34ξ )u1 +2u2+4u3 8(1ξ) u4 8ξ u5], u 02=u0+u1+4 u2+2u3 4u4 4u5.
From Eq. (17), we obtain the shape functions for the conforming crack front element with the crack front at edge ξ+η =1 as follows:
{Nscoll0(ξ,η)=ξ( 2ξ 1) 0.5(4ξ1)1ξη +0.5(1ξη),Nscoll1(ξ,η)=(1ξ)(1 2ξ)0.5(34ξ) 1ξη +0.5(1ξη),N scoll2(ξ ,η)=1 ξη+2 (1 ξη) , Nscoll3( ξ,η)=4ξ(1ξ)21ξη+2 (1 ξη) , Nscoll4( ξ,η)=4(1ξ)1ξη2(1 ξη) , Nscoll5( ξ,η)=4ξ 1ξη 2( 1ξ η).
In the application of the dual BEM for crack problems, due to the existence of hyper-singular integrals, the conforming crack front element should be transformed into a non-conforming one to retain Hölder’s continuity condition through the following linear transformation:
[ N Nscoll0 NNscoll 1NNscoll2 NNscoll3 NNscoll4 NNscoll 5] T=[ Nscoll0 Nscoll 1Nscoll2 Nscoll3 Nscoll4 Nscoll 5] T [N scoll0( P0)N scoll1( P0)N scoll2( P0)N scoll3( P0)N scoll4( P0)N scoll5( P0) Nscoll0(P 1) Nscoll1(P1 ) Nscoll2( P1)N scoll3( P1)N scoll4( P1)N scoll5( P1) Nscoll0(P 2) Nscoll1(P2 ) Nscoll2( P2)N scoll3( P2)N scoll4( P2)N scoll5( P2) Nscoll0(P 3) Nscoll1(P3 ) Nscoll2( P3)N scoll3( P3)N scoll4( P3)N scoll5( P3) Nscoll0(P 4) Nscoll1(P4 ) Nscoll2( P4)N scoll3( P4)N scoll4( P4)N scoll5( P4) Nscoll0(P 5) Nscoll1(P5 ) Nscoll2( P5)N scoll3( P5)N scoll4( P5)N scoll5( P5)] 1,
where NNscollm(m=0,1 ,..., 5) denotes the shape functions of the non-conforming crack front element, P m (m=0, 1,...,5) denotes the parameter coordinate of the interpolation nodes in ( ξ,η)coordinate systems, which can be obtained by the position parameter. In this work, Eqs. (11) and (19) are applied to approximate the distribution of displacements near the crack front.
The crack front element may lie on another edge. In this case, we can modify the order of the nodes or construct another type of shape function in the same manner.

Calculation of SIFs

Quadrilateral and triangular non-conforming crack front elements can be conveniently integrated into the dual boundary element method (DBEM). To compute the SIFs along the crack front, we construct a local coordinate system for Point O on the crack front, as shown in Fig. 4. The classic one-point formula is used to compute the SIFs, which is similar to that in Eq. (20) [5]:
{ KI= E4(1ν 2) π2rcosφ( ub P+ ubP), KII= E4(1ν 2) π2rcosφ( un P+ unP), KIII= E4(1+ν)π2rcosφ (u tP+ u tP),
where E is Young’s modulus, ν is Poisson’s ratio, and φ is the angle between P+O (or P O) and n. P is the point on the crack surface, which may be P1, P2, or P3. By using φ, the distance rP2=P2O, rP3=P3O, K p2 and Kp3 in Eq. (21) can be obtained. Then, a two-point formula considering the whole element can be derived after performing linear extrapolation between ( rP2,K p2) and (r P3,Kp 3) at location r=0. We can compute the SIFs at the crack front by using Eq. (21):
KO=r p3 Kp2r p2Kp3 rp3r p2.
Fig.4 Local coordinate system at Point O

Full size|PPT slide

Numerical examples

Several numerical examples are presented in this section to verify the accuracy and efficiency of the newly developed crack front elements. In the first example, the two constructed non-conforming crack front elements are tested in the DBEM analysis of a penny-shaped crack under uniform or linear polynomial tractions. The crack open displacements, Δu, are obtained for a comparison with existing results. Efficiency is also verified. In the second example, a penny-shaped crack under uniform inclined traction is analyzed with DBEM. Three types of SIFs, which contain KI, KII, and KIII, are computed, and their definitions can be found in Ref. [27]. The numerical results obtained by the two presented elements are compared with the exact results. In the final example, straight edge crack problems are compared with existing finite element analysis and other results.

Penny-shaped crack under uniform or linear polynomial tractions

This example concerns a penny-shaped crack embedded in an infinite space, as illustrated in Fig. 5. The boundary of the penny-shaped crack is described as
x =acosω, y=asinω, z=0,
where ω varies from 0 to 2π, and a is the radius. The crack is opened by imposing traction σz z(x,y) along the z-direction symmetrically to the upper and lower surfaces of the penny-shaped crack. In this example, a=1.0, Young’s modulus E=1.0, and Poisson’s ratio ν=0.25 are evaluated. Three traction boundary conditions, namely, σz z(x,y)=σ 0, σzz (x,y) =σ0x, and σzz (x,y) =σ0y with constant σ 0, are imposed on the penny-shaped crack. The reference solutions to the considered problems are 2σ0/ π, 4σ0x/(3aπ), or 4σ0y/(3aπ), respectively [28]. In this example, 96 nine-node non-conforming elements are used inside the crack surface and 32 nine-node non-conforming crack front elements are employed near the crack front, as shown in Fig. 5(a).
Fig.5 Mesh of the crack surface. (a) Non-conforming quadrilateral element; (b) non-conforming triangular element

Full size|PPT slide

Figure 6 shows that the KI calculated by non-conforming crack front elements of the quadrilateral type agrees well with the reference solutions. The largest errors are approximately 0.2%, 1.7%, and 2%, respectively.
Fig.6 Results of the normalized SIF KI calculated by non-conforming crack front elements of the quadrilateral type

Full size|PPT slide

Then, we perform an analysis using the triangular elements. As shown in Fig. 5(b), on the crack surface, 162 six-node non-conforming elements, including 26 six-node non-conforming crack front elements, are involved in the analysis.
Figure 7 indicates that the KI computed by employing non-conforming crack front elements of the triangular type agrees well with the reference solutions. The largest errors are about 0.9%, 2%, and 2%, respectively.
Fig.7 Results of the normalized SIF KI calculated by non-conforming crack front elements of the triangular type

Full size|PPT slide

Figure 8 indicates that the Δu calculated by non-conforming crack front elements of the quadrilateral and triangular type agrees well with the reference solutions [29]. The largest errors are about 1.5% in the two cases.
Fig.8 Results of Δu over the crack surface obtained using our method

Full size|PPT slide

Figure 9 shows that even when position parameter k varies from 0.1 to 0.3, the SIFs still retain considerable accuracy with the analytical solution. Notably, k=1λ for crack front elements of the quadrilateral type. The largest errors are about 1.2% in the two cases.
Fig.9 Influence of position parameter k

Full size|PPT slide

When the two types of crack front elements are used, high accuracy can be achieved within 1% by using fewer elements than Liu’s work [30]. In Liu’s work, 4704 elements are employed and the error is about 1.66%, while less than 200 elements are employed in our work, and the errors are within 1% stably.

Penny-shaped crack under uniform inclined traction

This example concerns the same crack as that in the first example. The load in this example, however, is different from that in the first numerical example. In this example, the crack is under a uniform inclined traction, as illustrated in Fig. 10. The parameters illustrated in Fig. 10 are evaluated by a=1.0, σ=1.0, E=1.0, γ=π /4, and ν=0.25. Here, γ is the angle between the uniform inclined traction and coordinate plane xoy. The same meshes as those in the first example are applied. The exact SIFs for this problem can be found in Tada et al.’s work [31]. The numerical and exact results are illustrated and compared in Figs. 11 and 12. The results obtained by our method agree well with the exact solutions. The largest errors are about 0.3%, 0.7%, 0.7%, 1%, 2%, and 2%, respectively.
Fig.10 Penny-shaped crack under uniform inclined traction

Full size|PPT slide

Fig.11 Results of three normalized SIFs KI, KII, and KIII calculated by non-conforming crack front elements of the quadrilateral type

Full size|PPT slide

Fig.12 Results of three normalized SIFs KI, KII, and KIII calculated by non-conforming crack front elements of the triangular type

Full size|PPT slide

Straight-edge crack compared with existing finite element analysis

In this example, a single edge crack, as illustrated in Fig. 13, is considered. The geometry parameters are c/b=1, t/c=3,and h/b=6. Uniform tension is imposed on two ends. The results are compared with those of the plane strain solution obtained by Mi and Aliabadi [5], Murakami and Keer [32], and Raju and Newman [33]. A total of 96 and 88 elements are used for the crack faces. The crack front elements proposed in this study are employed. In the work of Mi and Aliabadi [5], BEM was employed, whereas in the work of Raju and Newman [33], FEM was employed. The plane strain solution is 2.827, which can be found in the work of Murakami and Keer [32]. In our work, when z=0, the result is 2.2834 when crack front elements of the quadrilateral type are used; the result is 2.828 when crack front elements of the triangular type are employed. In Raju and Newman’s work [33], the SIF result is 2.776. Figure 14 shows that the results in the two cases agree well with those of Mi and Aliabadi [5] when z/ω changes from 0 to 0.4. When z/ω approaches 0.5, the results in the two cases agree well with those of Raju and Newman [33]. This example demonstrates the credibility of our method.
Fig.13 Geometry model of the singular edge crack and mesh

Full size|PPT slide

Fig.14 Normalized stress intensity factors along the crack front for the single edge crack

Full size|PPT slide

Conclusions

This study developed a family of non-conforming crack front elements of quadrilateral and triangular types in DBEM for 3D crack problems. In the deduction of the special shape functions, the asymptotic behavior of the distribution of the displacements and tractions near the crack front is considered. In the neighborhood of the crack front, the traditional shape function is expanded into Taylor expansion. The expanded shape function is then modified by considering the asymptotic behavior of the physical quantities. For the non-conforming quadrilateral crack front elements, the relations between the 2D and 3D Lagrange-type elements are determined. For the non-conforming triangular crack front elements, the transformation matrix between the non-conforming and conforming elements is introduced in consideration of the positions of the collocation nodes. Then, numerical examples using these non-conforming crack front elements are provided to verify the efficiency of our method. The numerical results are in good agreement with existing results. The proposed method may serve as a reference for addressing V-notch problems.

Acknowledgements

This work was supported in part by the National Natural Science Foundation of China (Grant Nos. 11602229 and 11602082), Hunan Provincial Natural Science Foundation of China (Grant No. 2017JJ3061), and Key Scientific and Technological Project of Henan Province (Grant No. 192102210227).
1
Zhang D H, Li Y Q, Xie G Z, Digital image correlation method for measuring deformations of vinyl chloride-coated metal multilayer sheets. Modern Physics Letters B, 2019, 33(5): 1950050

DOI

2
Zhang D H, Li Y Q, Liu J X, A novel 3D optical method for measuring and evaluating springback in sheet metal forming process. Measurement, 2016, 92: 303–317

DOI

3
Zhang D H, Xie G Z, Li Y Q, Strain and mechanical properties of the VCM multilayer sheet and their composites using digital speckle correlation method. Applied Optics, 2015, 54(25): 7534–7541

DOI

4
Zhang D H, Bai D P, Liu J B, Formability behaviors of 2A12 thin-wall part based on DYNAFORM and stamping experiment. Composites Part B: Engineering, 2013, 55: 591–598

DOI

5
Mi Y, Aliabadi M H. Dual boundary element method for three-dimensional fracture mechanics analysis. Engineering Analysis with Boundary Elements, 1992, 10(2): 161–171

DOI

6
Sladek V, Sladek J, Tanaka M. Nonsingular BEM formulations for thin-walled structures and elastostatic crack problems. Acta Mechanica, 1993, 99(1–4): 173–190

DOI

7
Zhang Y M, Gu Y, Chen J T. Boundary element analysis of the thermal behaviour in thin-coated cutting tools. Engineering Analysis with Boundary Elements, 2010, 34(9): 775–784

DOI

8
Yao Z H, Wang H T. Some benchmark problems and basic ideas on the accuracy of boundary element analysis. Engineering Analysis with Boundary Elements, 2013, 37(12): 1674–1692

DOI

9
Zhang J M, Lin W C, Dong Y, A double-layer interpolation method for implementation of BEM analysis of problems in potential theory. Applied Mathematical Modelling, 2017, 51: 250–269

DOI

10
Feng S Z, Han X, Wang G. An efficient on-line algorithm for the optimal design of multi-material structures under thermal loads. International Journal of Thermal Sciences, 2018, 132: 567–577

DOI

11
Feng S Z, Cheng Y H. An element decomposition method for heat transfer analysis. International Journal of Heat and Mass Transfer, 2018, 123: 437–444

DOI

12
Feng S Z, Bordas S P A, Han X, A gradient weighted extended finite element method (GW-XFEM) for fracture mechanics. Acta Mechanica, 2019 (in press)

DOI

13
Cheng C, Niu Z, Recho N. Analysis of the stress singularity for a bi-material V-notch by the boundary element method. Applied Mathematical Modelling, 2013, 37(22): 9398–9408

DOI

14
Cheng C, Ge S, Yao S, et al. Thermal stress singularity analysis for V-notches by natural boundary element method. Applied Mathematical Modelling, 2016, 40(19–20): 8552–8563

DOI

15
Zhu B J, Qin T Y. Application of hypersingular integral equation method to three-dimensional crack in electromagnetothermoelastic multiphase composites. International Journal of Solids and Structures, 2007, 44(18–19): 5994–6012

DOI

16
Henshell R D, Shaw K G. Crack tip finite elements are unnecessary. International Journal for Numerical Methods in Engineering, 1975, 9(3): 495–507

DOI

17
Lv J H, Jiao Y Y, Wriggers P, Efficient integration of crack singularities in the extended finite element method: Duffy-distance transformation and conformal preconditioning strategy. Computer Methods in Applied Mechanics and Engineering, 2018, 340: 559–576

DOI

18
Sáez A, Gallego R, Dominguez J. Hypersingular quarter-point boundary elements for crack problems. International Journal for Numerical Methods in Engineering, 1995, 38(10): 1681–1701

DOI

19
Hong H K, Chen J T. Derivations of integral equations of elasticity. Journal of Engineering Mechanics, 1988, 114(6): 1028–1044

DOI

20
Mi Y, Aliabadi M H. Discontinuous crack-tip elements: Application to 3D boundary element method. International Journal of Fracture, 1994, 67(3): R67–R71

DOI

21
Pan E, Yuan F G. Boundary element analysis of three-dimensional cracks in anisotropic solids. International Journal for Numerical Methods in Engineering, 2000, 48(2): 211–237

DOI

22
Xie G Z, Zhang J, Huang C, A direct traction boundary integral equation method for three-dimension crack problems in infinite and finite domains. Computational Mechanics, 2014, 53(4): 575–586

DOI

23
Li J, Feng W Z, Gao X W. A dual boundary integral equation method based on direct evaluation of higher order singular integral for crack problems. Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(2): 387–398 (in Chinese)

24
Xie G Z, Zhang D, Meng F, Calculation of stress intensity factor along the 3D crack front by dual BIE with new crack front elements. Acta Mechanica, 2017, 228(9): 3135–3153

DOI

25
Ariza M P, Saez A, Dominguez J. A singular element for three-dimensional fracture mechanics analysis. Engineering Analysis with Boundary Elements, 1997, 20(4): 275–285

DOI

26
Ariza M P, Dominguez J. Boundary element formulation for 3D transversely isotropic cracked bodies. International Journal for Numerical Methods in Engineering, 2004, 60(4): 719–753

DOI

27
Aliabadi M H, Rooke D P. Numerical Fracture Mechanics. Boston: Kluwer Academic Publishers, 1991

28
Shah R C, Kobayashi A S. Stress intensity factor for an elliptical crack under arbitrary normal loading. Engineering Fracture Mechanics, 1971, 3(1): 71–96

DOI

29
Chen L S, Kuang J H. A displacement extrapolation method for determining the stress intensity factors along flaw border. International Journal of Fracture, 1992, 57(4): R51–R58

DOI

30
Liu Y. On the displacement discontinuity method and the boundary element method for solving 3-D crack problems. Engineering Fracture Mechanics, 2016, 164: 35–45

DOI

31
Tada H, Paris P C, Irwin G R. The Stress Analysis of Cracks Handbook. Vol. 130. New York: ASME, 2000

32
Murakami Y, Keer L M. Stress intensity factors handbook, Vol. 3. Journal of Applied Mechanics, 1993, 60(4): 1063

DOI

33
Raju I S, Newman J C Jr. Three Dimensional Finite-Element Analysis of Finite-Thickness Fracture Specimens. NASA Technical Note, NASA TN D-8414. 1977

Outlines

/