arxiv: v2 [math-ph] 9 Jul 2012

Transcription

1 Stroh formalism in analysis of skew-symmetric and symmetric weight functions for interfacial cracks L. Morini ), E. Radi ), A.B. Movchan ), and N.V. Movchan ) ) Dipartimento di Scienze e Metodi dell Ingegneria, Universitá di Modena e Reggio Emilia, Via Amendola, 400, Reggio Emilia, Italy ) Department of Mathematical Sciences, University of Liverpool, Liverpool L69 3BX, U.K. arxiv:0.548v [math-ph] 9 Jul 0 Abstract The focus of the article is on analysis of skew-symmetric weight matrix functions for interfacial cracks in two dimensional anisotropic solids. It is shown that the Stroh formalism proves to be an efficient approach to this challenging task. Conventionally, the weight functions, both symmetric and skew-symmetric, can be identified as a non-trivial singular solutions of the homogeneous boundary value problem for a solid with a crack. For a semi-infinite crack, the problem can be reduced to solving a matrix Wiener-Hopf functional equation. Instead, the Stroh matrix representation of displacements and tractions, combined with a Riemann-Hilbert formulation, is used to obtain an algebraic eigenvalue problem, that is solved in a closed form. The proposed general method is applied to the case of a quasistatic semi-infinite crack propagation between two dissimilar orthotropic media: explicit expressions for the weight matrix functions are evaluated and then used in the computation of complex stress intensity factor corresponding to an asymmetric load acting on the crack faces. Keywords: Interfacial crack, Riemann-Hilbert problem, Stroh formalism, Weight functions, Stress intensity factor. Introduction Evaluation of coefficients in asymptotic representations of displacements and stress fields represents an important issue for addressing vector problems of crack propagation in elastic materials Bercial-Velez et al., 005; Mishuris & Kuhn, 00). The explicit derivation of weight functions is fundamental for the evaluation of stress intensity factors corresponding to a general distribution of forces acting on thecrack faces, as well as for thecalculation of higher order coefficients in the asymptotic expressions of the fields. The latter are used in asymptotic models of incremental crack growth, and hence are essential for evaluation of the crack path and analysis of fracture stability. In the work by Bueckner Bueckner, 985, 989), weight functions for several types of cracks in homogeneous elastic media, both in two and three dimensions, have been defined. In this paper, the term symmetric load is associated to forces of the same magnitude applied on both crack faces in opposite directions, while the load generated by forces acting on both crack faces in the same direction is called skew-symmetric or anti-symmetric. For cracks in homogeneous elastic materials, in the two-dimensional setting the skew-symmetric loading does not contribute to stress intensity factors, whereas it becomes relevant and it must to be accounted for in three-dimensional solids Bueckner, 985; Meade & Keer, 984). The situation is different when the crack is placed at the interface between two dissimilar elastic materials:

2 even for two-dimensional problems the skew-symmetric loads generate a non-zero contribution to stress intensity factors Bercial-Velez et al., 005; Lazarus & Leblond, 998; Piccolroaz et al., 009). In particular, for Mode III interfacial cracks the stress components do not oscillate, but a non-vanishing skew-symmetric component of the weight function still has to be accounted Piccolroaz et al., 009, 00). In the case of isotropic media, weight functions for semi-infinite cracks can be defined as singular non-trivial solutions of the homogeneous boundary problem with zero tractions on the crack faces but unbounded elastic energy Willis & Movchan, 995). For interfacial cracks between dissimilar isotropic elastic media, the weight functions are well discussed in literature Lazarus & Leblond, 998; Piccolroaz et al., 009, 00). The problem is generally reduced to a functional equation of Wiener-Hopf type, and its solution gives the symmetric weight function matrix Antipov, 999), while the skew-symmetric component is obtained by the construction of the corresponding full-field singular solution of the elasticity boundary value problem discussed in Piccolroaz et al. 009). For interfacial cracks between anisotropic dissimilar elastic media, although weight functions have been derived by Gao 99, 99) and Ma & Chen 004), the results on skew-symmetric weight functions are not readily available. In this article we illustrate a general procedure for the calculation of symmetric and skewsymmetric weight functions matrices for semi-infinite two-dimensional interfacial crack problems in anisotropic elastic bi-materials. It is shown that the challenging analysis of the matrix functional Wiener-Hopf equation can be replaced by solving a matrix eigenvalue problem deduced via an equivalent formulation, based on Stroh representation of displacement and stress fields Stroh, 96). By means of this new approach, general expressions for weight functions matrices, valid for plane interfacial cracks problems between any anisotropic media, are obtained. This general result is reported and discussed in details in Section 3). Section 4) illustrates the proposed method for the case of a semi-infinite two-dimensional stationary crack between two dissimilar orthotropic materials under plane stress deformation Suo, 990b), corresponding to a general non-symmetric load. Stroh representations for displacements and stress corresponding to this problem are explicitly calculated and used for deriving symmetric and skew-symmetric weight functions matrices. In Section 5) both symmetric and skew-symmetric weight functions are utilized together with Betti formula in order to evaluate complex stress intensity factor for an interfacial crack in orthotropic bi-material subject to an asymmetric loading configuration. In the particular case of isotropic media, the obtained result is consistent with the stress intensity factor derived for the non-symmetric distribution of forces obtained by Piccolroaz et al. 009). Finally, in Appendix A, the evaluated skew-symmetric weight function is compared to those calculated by the full field singular solution of the plane elasticity problem, following the approach illustrated in Piccolroaz et al. 009). The perfect agreement detected between the expressions derived by means of two alternative formulations represents an important benchmarking for the correctness of the obtained results. Interfacial cracks: preliminary results Here we introduce the main notations and the mathematical framework of the model. Consider a quasi-static advance of a semi-infinite plane crack between two dissimilar anisotropic elastic

3 materials with asymmetric loading applied to the crack faces, the geometry of the system is shown in Fig.). The loading is defined via tractions acting upon the crack faces. Considering a Cartesian coordinate system with the origin at the crack tip, see Fig.)), traction components behind the crack front are then defined as follows: σ ± j x,0 ± ) = p ± j x ) for x < 0 j =,, ) where p ± j x ) are given functions. Since the load is assumed to be self-balanced, its resultant force and moment vectors are equal to zero. Moreover, we assume that forces are applied outside a neighborhood of the crack tip and vanish at the infinity. The body forces are assumed to be zero. Symmetric and skew-symmetric parts of the loading are given by the following expressions: p j x ) = ) p + j x )+p j x ), [p j ]x ) = p + j x ) p j x ), j =,, ) The solutions in form of functions which vanish at the infinity and possess finite elastic energy are sought. Expressions for the stress field and displacements for a semi-infinite interfacial crack in anisotropic bi-materials have been obtained by means of Stroh formalism Stroh, 96) by Suo 990b). This approach to the physical problem of the crack, that is used for the derivation of weight function matrix in the paper, is reported in Section.). In Section.) the weight functions are defined as non-trivial singular solutions of the homogeneous elasticity problem for interfacial cracks with zero tractions on the crack faces but unbounded elastic energy, following the approach of Piccolroaz et al. 009). Section.3) reports the application of the Betti formula to physical fields and weight functions and the derivation of the fundamental identity in the Fourier space, discussed in details in Willis & Movchan 995); Piccolroaz & Mishuris 0) and Piccolroaz et al. 007). This integral identity will be used in the paper for the evaluation of coefficients in the asymptotics of the stress fields near the crack tip by means of a procedure based on weight functions theory Piccolroaz et al., 009). The asymptotic representations of physical displacements and stress fields near the crack tip, including the stress intensity factor as the coefficient of the first term Piccolroaz & Mishuris, 0; Piccolroaz et al., 009), are introduced in Section.4).. Stroh formalism in analysis of interfacial cracks Physical displacements and stress fields for an intefacial crack between two different anisotropic materials can be derived by means of the Lekhnitskii or Stroh approaches Lekhnitskii, 963; Stroh, 96; Suo, 990b). Following the article by Suo Suo, 990b), we introduce the stress vectors t j = σ j,σ j ) T together with the displacement u = u,u ) T. Theconstitutive relations for both the elastic media occupying respectively the upper and lower half-planes can be written using the Stroh formulation Stroh, 96): t = Qu, +Ru, 3) t = R T u, +Tu,, 4) 3

4 Figure : Geometry of the model. where the matrices Q, R and T depend on the material constants. A semi-infinite static crack placed at the interface between the two materials is considered, as it is illustrated in Fig.). Accordingly, thederivativeofthedisplacementu, x,x )andthetraction τx,x ) = t x,x ) can be written as τx,x ) = Bgz) +Bgz) 5) and u, x,x ) = Agz)+Agz), 6) where A and B are constant matrices, gz) is an analytic vector function with components g j x + µ j x ), and µ j are complex numbers with positive imaginary parts. According to Suo 990b), if g j z j ) is an analytic function of z j = x + µ j x in the upper half-plane or in the lower half-plane) for one µ j, where µ j is a complex number with positive imaginary parts, it is analytic for any µ j. On the basis of this property, here and in the text that follows, we reduce the analysis to a single complex variable. The connection between the elements of A, B and Q,R,T is given by the following relations see Ting 996), pages 70, 7): Q ik +R ik +R ki )µ j +T ik µ j)a kj = 0 7) B ij = R ki +µ j T ik )A kj 8) Thus, each column of A is a non-trivial solution of the eigensystem 7), while the eigenvalues µ j are roots of the characteristic equation: Q ik +R ik +R ki )µ j +T ik µ j = 0 9) In turn, for each of the two phases we introduce the Hermitian matrix Y = iab, which will be used in the text below. Further in the text, we shall use the superscripts ) and ) to denote the quantities related to the upper and lower half-planes, respectively. The expression 5), in the limit x 0 ±, leads to a non-homogeneous Riemann-Hilbert problem: h + x ) h x ) = τx ), x R, 0) 4

5 A crack advancing along the negative semi-axis x < 0 is considered, the traction-free condition is imposed for x < 0, while the continuity of the tractions at the interface is assumed ahead the crack front. The following equations are satisfied on the real axis Suo, 990b): h + x )+H Hh x ) = τx ) for x > 0 ) h + x )+H Hh x ) = 0 for x < 0 ) The detailed analysis of this problem is included in Suo 990b); Suo et al. 99), and it shows that the stress and displacement fields do not have oscillation near the crack tip for the case when the matrix H = Y ) + Y ) is real, otherwise the stress and displacement are characterized by the oscillatory behavior near the crack tip. The branch cut for the function hz) is assumed to bealong the crack linex < 0. Assumingthat thestress vanish at theinfinity, the solution of the homogeneous Riemann-Hilbert problem ) is hz) = wz +iε, where w is solution of the following eigenvalues problem Gao et al., 99; Suo, 990b; Suo et al., 99): Hw = e πε Hw. 3) The traction τ ahead of the crack has the form: ) τx ) = Re Kx iε w. 4) πx where K = K I +ik II is the complex stress intensity factor including both mode I and mode II contributions to the traction, ε is the bi-material parameter a real non-dimensional number measuring an aspect of elastic dissimilarity of the two materials), and w is the eigenvector obtained from 3). For the case of plane strain load, it has been shown in Suo Suo, 990b; Suo et al., 99) that the displacement jump [u] across a semi-infinite crack running along the negative semi-axis x < 0 is given in the form. Weight functions definition ) /H+H) x ) [u]x ) = π coshπε Re K x ) iε ) w +iε Following the theory developed by Willis and Movchan Willis & Movchan, 995), we define a vector functionu = U,U ) T asthesingularsolutionoftheelasticity problemwithzerotraction on the faces where the crack is placed along the positive semi-axis x > 0. The traces of these functions on the plane containing the crack are known as the weight functions Piccolroaz et al., 009), and notations [U] and U will be used in the paper to denote symmetric and skewsymmetric weight functions respectively: 5) [U]x ) = Ux,x = 0 + ) Ux,x = 0 ) 6) U x ) = Ux,x = 0 + )+Ux,x = 0 )) 7) 5

6 The traction vector Σ = Σ,Σ ) T associated to the singular solutions is continuous on the plane containing the crack see Piccolroaz et al. 007, 009)) and vanishes for x < 0 homogeneous boundary conditions are imposed). In practice, for singular solutions we impose traction-free condition for x > 0 and traction continuity at the interface for x < 0: h + x )+H Hh x ) = 0 for x > 0 8) h + x )+H Hh x ) = Σx ) for x < 0 9) It is important to note that the domain of singular solutions and then of weight functions is different from the domain of the physical solution, defined in previous Section, where the crack is placed along the negative semi-axis. For the singular solution, the branch cut for hx ) is chosen to be along the line x > 0. Using the representation hz) = vz 3 +iε we reduce 9) to the eigenvalues problem: Hv = e πε Hv. 0) The bi-material matrix H is Hermitian, and the eigenvector corresponding to the eigenvalue ε is v = w Suo, 990b). The singular traction Σ is sought in the form: Σ = x ) 3 Re Cw x ) iε), ) π where C = C I +ic II is a complex constant representing both mode I and mode II contributions to the traction ), as it is for the complex stress intensity factor K with regard to the physical traction see eq.4)). The shear and normal opening modes for plane stress and plane strain problems are coupled. Two linearly independent vectors Σ and associated weight functions can be identified similar to Piccolroaz & Mishuris 0); Piccolroaz et al. 009)..3 Fundamental Betti identity As discussed above, U is discontinuous along the positive semi-axis x > 0, whereas u is discontinuous along the positive semi-axis x < 0. Σ is zero for x > 0, whereas τ is zero for x < 0. Asymmetric loading are applied to the crack faces see Fig.)), thus the physical traction acting at the interface on the entire x axis can be written as: σx,x = 0 + ) = p + x )+τx ), σx,x = 0 ) = p x )+τx ), ) According to the approach illustrated in Willis & Movchan 995) and Piccolroaz et al. 007), in order to derive explicit formulas for calculating the coefficients for the asymptotic of the stress fields near to the crack tip, we apply the Betti formula to the physical fields and to weight functions. The following relations, respectively in the upper and in the lower-midplane, are obtained: x =0 + x =0 } U T x x,0 + )Rσx,0 + ) Σ T x x,0 + )Rux,0 + ) dx = 0 3) } U T x x,0 )Rσx,0 ) Σ T x x,0 )Rux,0 ) dx = 0 4) 6

7 where R is the rotation matrix: R = ) 0, 0 Subtracting the 4) from the 3) and using ), we derive the following integral identity: x =0 U T x x,0 + )Rp + x )+U T x x,0 + )Rτx ) U T x x,0 )Rp x ) } U T x x,0 )Rτx ) [Σ T x x,0 + )Rux,0 + )+Σ T x x,0 )Rux,0 )] dx = 0 Referring to 6) and 7) we deduce = x =0 x =0 } [U] T x x,0 + )Rτx ) Σ T x x,0 + )R[u]x ) = } [U] T x x,0 + )R p x ) Σ T x x,0 + )R[p]x ). 5) Let us define the Fourier transform of the skew-symmetric weight function and associated traction Σ respect to x : [Û]+ ξ) = + 0 [U]x )e iξx dx, ˆΣ ξ) = 0 Σx )e iξx dx, 6) Where the superscript + indicates that [Û] ξ) is analytical in the upper half-plane Imξ 0,+ )) and the superscript indicates that ˆΣ ξ) is analytical in the lower half-plane Imξ,0)). The physical traction and the jump function are defined in such a way that the transforms ˆτ + ξ) and [û] ξ) are analytic in the upper and in the lower half-plane, respectively. Applying the Fourier transform to the 5), by means of convolutions properties, the following relation, valid for ξ R, is derived Piccolroaz & Mishuris, 0; Piccolroaz et al., 007): [Û]+T Rˆτ + ˆΣ T R[û] = [Û]+T R ˆp Û +T R[ˆp], ξ R, 7) This identity, obtained by Willis & Movchan 995) and Piccolroaz et al. 007), relates transforms of the physical solutions to transforms of the weight functions, and it is used for the evaluation of the stress identity factors..4 Asymptotic representation of the fields and stress intensity factors The physical traction 4) and the displacement jump across the crack face 5) for x 0, can be written in the matrix form, as follows: τx ) = x π T x )K+ x π T x )J+ x π T x )L+Ox 5 ) 8) [u]x ) = x ) π Ux )K+ x ) 3 π Ux )J+ x ) 5 π Ux )L+O x ) 7 ) 9) 3 7

8 Where K = K,K),J = J,J),L = L,L), and K and L are higher order coefficients, defined in the same way of the stress intensity factor: J = J I + ij II,L = L I + il II. Matrices T x ) and Ux ) are given by: ) T x ) = wx iε,wx iε, Ux ) = H+H) w x ) iε cosh πε +iε, w x ) ) iε, iε The evaluation of coefficients in the asymptotic formulas8) and9) is included in Piccolroaz & Mishuris 0); Piccolroaz et al. 009) for interfacial cracks between dissimilar isotropic media, by means of a general integral formula involving symmetric and skew-symmetric weight functions matrices. Here this approach is extended to anisotropic bi-materials: in the next Section general expressions for Fourier transforms of symmetric and skew-symmetric weight functions are derived using the Stroh formulation, as defined in Section.), while in Section 4) these expressions will be specialized to the case of a crack between two orthotropic materials Suo, 990b) and the general formula obtained by Piccolroaz et al. Piccolroaz & Mishuris, 0; Piccolroaz et al., 009) will be applied for calculating stress intensity factor for this case. 3 Symmetric and skew-symmetric weight functions for anisotropic bi-materials Here we derive general expressions for symmetric and skew-symmetric weight function defined in Section.) as the jump [U] and the average U of the singular solution for a semi-infinite interfacial crack with free-traction conditions. Let us introduce the vector function hz), solution of the Riemann-Hilbert problem for the physical traction 0), such that: hz) = Bgz), Im z 0, Bgz), Im z < 0. According to the Plemelj formula, hz) can be written as follows: hz) = τη) πi η z dη, where the integral is understood in the principal value sense. Taking the Fourier transform respect to x, for x = 0 ±, at the interface, we obtain: ĥξ,x = 0 ± ) = ĥ± ξ) = ± e iξη H±ξ)τη)dη = ±H ξ)ˆτ + ξ), ξ R, where H is the Heaviside function, and the traction τ is defined in such a way that its transform is analytic in the upper half-plane. To derive the above representation we used the following relations e iξx ±i0) η x ±i0) dx = ±πie iξη H ξ). 8

9 and Hence the Fourier transforms of g and g at the interface are obtained in the form ĝξ) = B ĥ + ξ) = H ξ)b ˆτ + ξ), ξ R, ĝξ) = B ĥ ξ) = Hξ)B ˆτ + ξ), ξ R. By applying the Fourier transform to 6) for the derivative of the physical displacements, we deduce: iξûξ,x = 0 ± ) = Aĝξ)+Aĝξ), ξ R Hence, the Fourier transform of the displacements on the boundary of the upper half-plane is given by: ûξ,x = 0 + ) = i H ξ)ab +Hξ)A B } ˆτ + ξ) ξ = ξ By expressing the Heaviside function in the form: H ξ)y ) Hξ)Y )}ˆτ + ξ). H±ξ) = ) ±sign ξ we obtain: ûξ,x = 0 + ) = ξ Y) Y ) ) ξ Y) +Y ) } ) ˆτ + ξ), ξ R. 30) Following the same pattern of the derivation as above, we derive the Fourier transform of the physical displacement on the boundary of the lower half-plane: ûξ,x = 0 ) = ξ Y) Y ) )+ ξ Y) +Y ) } ) ˆτ + ξ), ξ R. 3) Next, we replace u in the above text by the singular solution U for a semi-infinite interfacial crack. Correspondingly, the vector of tractions τ is replaced by Σ, which is the traction vector corresponding to the singular solution U. The Fourier transforms of the singular displacements on the boundary can then be derived: Ûξ,x = 0 + ) = Ûξ,x = 0 ) = ξ Y) Y ) ) ξ Y) +Y ) ) }ˆΣ ξ), ξ R. 3) ξ Y) Y ) )+ ξ Y) +Y ) ) }ˆΣ ξ), ξ R. 33) According to definitions 6) and 7), we take the difference and the average of 3) and 33), and we derive that: [Û]+ ξ) = i signξ) ImY ) Y ) ) ReY ) +Y ) )}ˆΣ ξ), 34) ξ 9

10 and Û ξ) = i signξ) ImY ) +Y ) ) ReY ) Y ) )}ˆΣ ξ), ξ R. 35) ξ We note that 34) is the functional equation of the Wiener-Hopf type, similar to the one studied in Piccolroaz et al. 009, 007) and Piccolroaz & Mishuris0) for the case of isotropic media. Using the continuity of tractions Σ across the interface, together with 34) and 35), we express the Fourier transform of the skew-symmetric weight function U in the form: Û ξ) = A[Û]+ ξ)+ i ξ BˆΣ ξ), ξ R, 36) where the diagonal matrix A and the off-diagonal matrix B are A =, ReY) Y ) ) ReY ) +Y )) ) and B = ImY) +Y ) ) A ImY ) Y ) ). The formulas 34) and 35) give expressions for the Fourier transform of the symmetric and skew-symmetric weight functions for interfacial cracks in anisotropic bi-materials in terms of the transformed singular traction ˆΣ. These expressions are compared to those obtained by the construction of the full-field singular solution for the elasticity problem in an half-plane Piccolroaz et al., 009) in Appendix A. As expected, the weight functions derived via two different approaches agree. Making the inversion of the34) and the35) the explicit representation of the weight functions can be obtained and the complex stress intensity factors associated to an arbitrary system of forces can be evaluated by means of integral formulas derived by the application of the Betty identity Piccolroaz et al., 009). 4 Interfacial crack in orthotropic bi-materials In this section explicit expressions for the Fourier transforms of the symmetric and skewsymmetric weight functions for an interfacial crack in orthotropic bi-materials are derived using equations 34) and 35). 4. Stroh representation for orthotropic bi-materials For the case of orthotropic two-dimensional media, the matrices Q,R and T introduced in section.) are given by: ) ) ) c 0 0 c c66 0 Q =,R =,T =. 0 c 66 c c Where c ij are the elements of the materials stiffness matrix, which can be expressed in function of the elements of the compliance matrix s ij Suo, 990b): s s c = s s c = s s s 37) s s c = s s c 66 = 38) s s 66 0

11 The characteristic equation 0) becomes: s µ 4 +s +s 66 )µ +s = 0 39) The matrices A and B, defined by relations 8) and 8), are equivalent to those provided by alternative Lekhnitskii formulation Lekhnitskii, 963), more precisely the Lekhnitskii approach gives a specially normalized eigenvector matrix A, and the characteristics equation derived using the Lekhnitskii formalism is identical to the 39) Suo, 990b; Ting, 996). The relationships between the two alternative formulations are derived and reported in details in reference Hwu 993), where the elements of A and B in the Stroh representation are reported in function of the coefficient proposed by Lekhnitskii. Here we write the Stroh matrices A and B using this particular normalization, reported in Hwu 993): A = B = s µ +s µ s s µ 4 ) s µ + s µ µ s s µ 4 ) µ s µ s µ 4 ) s µ s µ 4 ) s µ +s µ s s µ 4 ) s µ + s µ µ s s µ 4 ) µ µ s s µ 4 ) µ s s µ 4 ) 40) 4) Where µ and µ are the two roots of the characteristics equation with positive imaginary part. The Hermitian matrix Y evaluated using 40) and 4) is: ) s Y = iab Imµ +µ ) is µ µ s ) = is µ µ s ) s Im + 4) µ µ) Following the notation introduced by Suo, Suo, 990b; Gupta et al., 99), we define two adimentional parameters measuring the material anisotropy: λ = s ρ = s +s 66 s s s If λ = the material has transversely cubic symmetry, while if λ = ρ = the material is transversely isotropic. The positive definetess of the strain energy density requires that: λ > 0, < ρ < +, The characteristics equation 39) then becomes: λµ 4 +ρλ µ + = 0 43) The roots with positive imaginary parts are: µ = iλ 4 n+m), µ = iλ 4 n m), for < ρ < +, µ = λ 4in+m), µ = λ 4in m), for < ρ <, µ = µ = iλ 4, for ρ =,

12 Where: n = ) +ρ), m = ρ) Using this notation the matrix Y becomes: ) Y = iab nλ 4s = s ) is s ) +s ) is s ) +s ) nλ 4s s ), 44) The bi-material matrix H for two orthotropic materials Suo, 990b; Gupta et al., 99) is: Where: H = Y ) +Y ) H = iβ ) H H iβ H H H H = [nλ 4s s ) ] ) +[nλ 4s s ) ] ), H = [nλ 4s s ) ] ) +[nλ 4s s ) ] ), β H H = [s s ) +s )] ) [s s ) +s )] ), β is the generalization of one of the Dundurs parameters Dundurs, 969), and is connected to the bi-material oscillatory index ε by the relation: ε = ) β π ln +β The eigenvector w of the eigensystem 3), associated with the component of displacement jump 5) and of the traction ahead of the crack 4), is assumed in the normalized form: 4. Weight functions for in-plane deformations 45) w = i, ) T H, 46) H In order to derive explicit expressions for the Fourier transforms of the weight functions matrices from the relations 34) and 35), we need to evaluate the Fourier transform of the singular traction Σ. The singular traction on the boundary Σ for an interfacial crack between two generic anisotropic materials is given by ). Here we report this expression: Σ = x ) 3 Re Cw x ) iε) 47) π For studying a crack between two orthotropic media, we assume for w the expression 46). As just anticipated, mode I and mode II are coupled, and associated to one single complex constant C. Therefore the spaces of singular tractions and of singular solutions U are twodimensional linear spaces, and two linearly independent vectors Σ and weight functions must be defined Piccolroaz et al., 009; Piccolroaz & Mishuris, 0). Assuming C I =,C II = 0,

13 and C I = 0,C II =, and substituing explicit components of w, the two independent traction vectors are defined: Σ x ) = x i x ) 3 ) iε x ) iε ) ; 48) π H H x ) iε + x ) iε ) Σ x ) = x ) 3 π Applying the Fourier transform to 48) and 49), the result is: ˆΣ ξ) = ξ +4ε ˆΣ ξ) = ξ +4ε x ) iε + x ) iε ) ; 49) i H H x ) iε x ) iε ) e 0ξ iε ) c + iε i H e0 ξ iε H c + iε ) e0 ξ i iε ) c + iε H H Where Imξ,0), e 0 = e επ and: e 0ξ iε c + iε ) e 0 = e επ, c ± = +i) π [ ], Γ ±iε ) + ξiε e 0 c +iε + ξiε e 0 c ) } +iε ; 50) ) } + ξiε e 0 c +iε + ξiε e 0 c ) } ; 5) +iε Substituing the 50) and 50) into 34) and 35), and expressing the matrices in term of elements of H, we derive the Fourier transforms of the two independent symmetric and skew symmetrix weight functions: [Û ]+ [Û ]+ [Û ]+ [Û ]+ = H H ξ H H iβsignξ) iβsignξ) [ˆΣ ] [ˆΣ ] ; 5) [ˆΣ ] [ˆΣ ] H H Û Û δ H H H H iγsignξ) [ˆΣ = ] [ˆΣ ] ; 53) Û Û ξ iγsignξ) δ H H [ˆΣ ] [ˆΣ ] Where the following Dundurs-like parameters have been defined: δ = [nλ 4s s ) ] ) [nλ 4s s ) ] ) H, δ = [nλ 4s s ) ] ) [nλ 4s s ) ] ) H, γ = [s s ) +s )] ) +[s s ) +s )] ) H H, 3

14 The superscript + indicates that Û+ ξ) is analytic in the upper half-plane Imξ 0,+ )). The expression 5) for the transform of the symmetric weight functions, as just anticipated in previous session, is a Wiener-Hopf type equation similar to the one solved for interfacial crack in isotropic materials in Piccolroaz & Mishuris 0); Piccolroaz et al. 007, 009). The skew-symmetric part of the weight function can also be decomposed by means of relation 53): Û = A[Û]+ + i ξ BˆΣ, 54) Where the matrices A and B for orthotropic bi-materials using our notation are given by: A = ) ) δ 0 H H 0 γ +βδ ; B = ; 0 δ γ +βδ ) 0 The explicit expressions for the weight functions are obtained by the inversion of the derived transforms. The symmetric weight function [U]x ) is equal to zero for x < 0, while for x > 0 it is given by: [ ] U x ) = [ ] U x ) = [ ] U x ) = [ ] U x ) = H x c + c β ) ) π+4ε ) +iε x iε +β +) ) iε i H H x c + c π+4ε ) ih x c + c π+4ε ) H H x c + c π+4ε ) β ) ) +iε x iε β +) ) iε β ) ) +iε x iε β +) ) iε β ) ) +iε x iε +β +) ) iε } x iε } x iε } x iε } x iε The skew-symmetric weight function U x ), is equal to A[U]x) for x > 0, while for x < 0 it is given by: U x ) = ih γ +βδ ) x ) 4c + c ) π+4ε ) +iε x ) iε e 0 U x H H γ +βδ ) x ) ) = 4c + c ) π+4ε ) +iε x ) iε e 0 U x ) = H γ +βδ ) x ) 4c + c ) π+4ε ) +iε x ) iε U x ) = i H H γ +βδ ) x ) 4c + c π+4ε ) e 0 +iε + } )e iε 0 x ) iε } )e iε 0 x ) iε } )e iε 0 x ) iε ) x ) iε e 0 55) + } )e iε 0 x ) iε 56) 4

15 4.3 Stress-intensity factor for orthotropic bi-materials The symmetric and skew-symmetric weight functions are now used for evaluating the complex stress factor for a loading correspondent to an arbitrary system of forces, following the procedure outlined in Piccolroaz et al. 007, 009) and Piccolroaz et al. 00). In Section.4, we have shown that complex stress identity factor represents the coefficient of the first order asymptotic term of physical traction field near the crack tip. Remembering equation 8), for x 0 + the traction becomes: τx ) = x π T x )K+ x π T x )J+ x π T x )L+Ox 5 ) 57) For orthotropic bi-materials, considering the components of w given by 46), the matrix T x ) is: ix iε ix iε T x ) =, H H x iε H H x iε The Fourier transform of the 57), as ξ and Imξ 0,+ ), is: + ˆτ + ξ) = ξ + 4 ˆT ξ + )K+ ξ 4ξ ˆT ξ + )J+ ξ 4ξ ˆT ξ + )L+Oξ 7 ) 58) + 3 Where: ˆT ξ + ) = ξ iε + c + e 0 H iξ iε + H c + e 0 e 0ξ iε + c H ie 0 ξ iε + H c Using this expression, the explicit transforms of the traction ˆΣ and of the symmetric weight functions matrix [Û]+ derived in previous Section, and evaluating [û], the asymptotic expansions for the members of Betti identity 7) are derived: [Û]+T Rˆτ + = ξ M K+ξ M J+Oξ 3 ) for Imξ 0,+ ) 59) ˆΣ T R[û] = ξ M K+ξ M J+Oξ 3 ) for Imξ,0) 60) The explicit form for the matrix M is: H M = 4c + c +4ε 4 ) β ) iε) e 0 iβ ) iε) e 0 e 0 β +)+iε) 6) β +)+iε) ie 0 Now, we rewrite the Fourier transform of the Betti identity 7) in terms of a Riemann- Hilbert problem: Ψ + ξ) Ψ ξ) = [Û]+T ξ)r ˆp ξ) Û +T ξ)r[ˆp]ξ), ξ R 6) 5

16 Where Ψξ), according to the Plemelj formula, is: Then the solution of the 6) is: Ψ ± ξ) = πi Ψη) dη, 63) η ξ [Û]+T Rˆτ + = Ψ +, Imξ 0,+ ); ˆΣ T R[û] = Ψ, Imξ,0); From these expressions asymptotic estimates can be extracted. For ξ we can expand Plemelj s formula, considering only the first term we have: Ψ ± ξ) = πi Ψη) dη = ξ η ξ πi Ψη)dη +Oξ ) 64) Substituing this expression and the expansions 59) and 60) in 6), comparing the corresponding terms, and considering only the first order, the following general formula for the complex stress intensity factor is obtained: K = πi M } [Û]+Tη)R ˆp η)+ Û +T η)r[ˆp]η) dη 65) Since we have explicit expressions for matrix M and for the Fourier transforms of both symmetric and skew-symmetric weight functions, now using the formula 65) we can evaluate the complex stress intensity factor associated to an arbitrary loading for interfacial cracks in orthotropic bi-materials. In the next section an illustrative example of computation of K by means of the 65) is reported. 5 An illustrative example In this section we present an illustrative example of computation of complex stress intensity factor K for an interfacial crack loaded by a simple asymmetric force system in orthotropic bi-materials. The considered force system is illustrated in Fig.): the loading consists in a point force F acting upon the upper crack face at a distance a behind the crack tip and two point forces F/ acting upon the lower crack face at a distance a b/ and a+b/, respectively, behind the crack tip. The loading can be expressed in terms of the Dirac delta function Piccolroaz et al., 009): p + x ) = Fδx +a), p x ) = F δx +a+b) F δx +a b), 66) As it is shown in Fig.), the loading can be decomposed into symmetric and skew-symmetric part: p x ) = F δx +a) F 4 δx +a+b) F 4 δx +a b) [p]x ) = Fδx +a)+ F δx +a+b)+ F δx +a b) 67) 6

17 Figure : Three point loading of an interfacial crack. The Fourier transform of the symmetric and skew symmetric part of the loading are given by: ˆp ξ) = F e iξa F 4 e iξa+b) F 4 e iξa b) [ˆp]ξ) = Fe iξa + F e iξa+b) + F e iξa b) 68) Using these expressions and the explicit transforms of the symmetric and skew-symmetric weight functions in orthotropic media, 5) and 54), both symmetric and anti-symmetric part of the complex stress intensity factor K = K S +K A corresponding to this loading system have been evaluated by means of the integral formula 65): K S = e 0 β F H H π a iε K A = e 0 γ +δ ) β) F H H π a iε + 4 b/a) iε + } 4 +b/a) iε 4 b/a) iε } 4 +b/a) iε In order to study the behavior of these symmetric and anti-symmetric contributions to the stress intensity factors in function of b/a, the following non-dimensional parameters have been defined Suo, 990b): Φ = [s s ) ] ) [ s = [s s ) ] ),Θ) s s ) ] ),Θ ) = [ s s s ) Then we can exprime H,H,δ and γ in function of these parameters: H = [nλ 4] ) +[nλ 4] ) Φ H [nλ4] ) +[nλ 4] ) Φ ] ), 7

18 δ = [nλ 4] ) [nλ 4] ) Φ [nλ 4] ) +[nλ 4] ) Φ γ = β [+Θ]) [+Θ] ) Φ Φ[+Θ] ) [+Θ] ) It is important to note that the Dundurs parameter γ, associated to the skew-symmetric part of the loading, depends on β. Moreover, oscillations of the stress and displacement fields are excluded for β = 0, ε = 0 Suo, 990b). In this case also γ = 0. The complex stress intensity factor has been computed for ρ ) = 0.74,/λ ) =,Θ ) = /,ρ ) = 4.9,Φ = 6.4 here we have considered that material ) is alluminium and material ) is boron Suo, 990a)), Θ ) = /,Θ ) = and five different values of the Dundurs oscillation parameter β = /4, /, 0, /4, /}. The values have been normalized multiplying by a F and plotted in Fig.3) in function of the ratio b/a. The symmetric stress intensity factor is reported on the left of the figure, while the skew-symmetric is on the right, the real part is reported on the top while the imaginary is on the bottom. Observing the figure, we note that both the real and the imaginary part of the skew-symmetric stress intensity factor are zero for b/a = 0, as expected, since for b = 0 the loading is symmetric. As we increase b/a, the skew-symmetric contribution to the loading become more relevant, and the skew symmetric stress intensity factor correspondingly increases. As we can see, in the case without oscillation, corresponding to β = 0, both the imaginary parts KII S and KA II vanish and the stress intensity factor becomes real. Both the symmetric and the skew-symmetric stress intensity factors diverge for b/a 0, because a point force is approaching the crack tip. In order to characterize the magnitude of the skew-symmetric stress intensity factor respect to the symmetric stress intensity factor, the ratio KI A/KS I is evaluated as a function of b/a in Fig.4). We observe that as b/a increases, KI A may reach the 40% of KI S, as a consequence we can say that the contribution of the skew-symmetric part of the loading is not negligible, and needs to be taken into account in perturbative analysis of interfacial cracks between two dissimilar anisotropic elastic materials subject to asymmetric forces systems applied on the crack faces Piccolroaz et al., 00, 009). 8

19 $! a #! "! b %a KI S /F "! #! & ='/ & = '/4 & = 0 %a KI A /F "! & =)/ & = )/4 #! & = 0 & = /4 #! & = /!! 0!" b/a! 0!" b/a!*!+ c!*!+ d %a KII S /F!"! & =)/ )!" & = )/4 & = 0 )!+ & = /4 & = / )!* 0!" b/a %a KII A /F!"! )!" )!+ )!* 0 & = )/ & = )/4 & = 0 & = /4 & = /!" b/a Figure 3: Symmetric and anti-symmetric stress intensity factors as functions of b/a computed for ρ ) = 0.74,/λ ) =,Θ ) = /,ρ ) = 4.9,/λ ) = 4.3,Θ ) =,Φ = 6.4 and different values of the Dundurs parameter β: β = /, β = /4, β = 0, β = /4, β = /. 9

20 KI A /KI S!&!%!$!# ' = / ' = /4 ' = 0 ' = /4 ' = /!"! 0!# Figure 4: Ratio K A I /K S I computed for ρ ) = 0.74,/λ ) =,Θ ) = /,ρ ) = 4.9,/λ ) = 4.3,Θ ) =,Φ = 6.4 and different values of the Dundurs parameter β: β = /, β = /4, β = 0, β = /4, β = / reported in function of b/a. 6 Conclusions The developed general approach for the derivation of symmetric and skew-symmetric weight function matrices for interfacial plane cracks between dissimilar anisotropic materials, based on Stroh formulation of displacements and stress fields, have been discussed in details and tested by means of the application to the case of a crack placed at the interface between two orthotropic materials under plane stress. The skew-symmetric weight functions obtained by means of the proposed method have been compared to those obtained for the same problem by the construction of the full-field singular solution of the elasticity problem in a half-plane Piccolroaz et al., 009), the comparison between the two different solutions is reported in Appendix A. The perfect equality detected between the expressions derived by means of two distinct approaches is an important ulterior proof for the obtained results. Since the proposed Stroh representation is valid for stationary and steady-state elasticity problems in many anisotropic media Stroh, 96; Suo, 990b; Ting, 996), it can be utilized for evaluating explicit weight functions for plane interfacial cracks in several kind of materials Stroh analysis have been proposed for example in quasi-crystals Radi & Mariano, 00), piezoelectrics Suo et al., 99) and poroelastics media Gautier et al., 0). The derived weight functions can be used in many important applications: in perturbative expansions for growing cracks or wavy cracks problems Piccolroaz et al., 007), and in the computation of the stress intensity factor for non-symmetric self-balanced load generated by a system of point-forces applied on the crack faces. An example of stress intensity factor evaluation for an asymmetric loading b/a 0

21 is reported in Section 5): weight functions matrices obtained for orthotropic bi-materials have been used in the computation, and the results show that the contribution of the skew-symmetric part of the loading is not negligible and must be considered in the asymptotic expressions of the stress near the crack tip, as it has already been demonstrated for the case of isotropic media Piccolroaz et al., 009, 00). 7 Acknowledgements The authors are grateful to Prof. G. Mishuris and Dr. A. Piccolroaz for the fruitful discussions about weight functions theory and their valuable suggestions. L. M. and E. R. gratefully acknowledge financial support from the Cassa di Risparmio di Modena in the framework of the International Research Project Modelling of crack propagation in complex materials. Appendix A In this appendix, the Fourier transforms of the singular solution of the interfacial crack problem between two dissimilar orthotropic materials are derived by solving a boundary value problem for a semi-infinite half-plane subjected to traction boundary conditions at its boundary, following the procedure illustrated in Piccolroaz et al. 009). The derived expressions for the singular displacements and for the symmetric and skew-symmetric weight functions are compared to those obtained in Section 3) by means of the direct solution of the Riemann-Hilbert problem 0). The perfect agreement detected between the expressions derived using two different approaches is an important test for the obtained results. Initially, we consider the lower half-plane, denoted in the article by the superscript ). Introducing the Fourier transform of the stresses respect to the variable x, we consider the component ˆσ as the primary unknown function, so that the plane strain elasticity problem for the orthotropic material ) is reduce to the following ordinary differential equation: s ) ˆσ ξ s ) 66 +s) )ˆσ +ξ4 s ) ˆσ = 0 69) Where a prime denotes the derivatives respect to x. The characteristic equation associated to 69) is: [ω ) ] 4 s ) ξ s ) ) 66 +s) [ω ) ] +ξ 4 s ) = 0, ξ R Introducing ν ) = ω ) / ξ and using the same notation of section 4, this characteristic equation becomes: λ ) [ν ) ] 4 ρ ) λ ) ) [ν ) ] + = 0 70) Assuming that this equation possesses four distinct roots, ρ ) ), the general solution of the plane strain elasticity problem in the lower half-plane is: ˆσ ξ,x ) = A ) e ξ ν) x +A ) e ξ ν) x, ˆσ ξ,x ) = ξ ˆσ, ˆσ ξ,x ) = i ξˆσ, 7)

22 The Fourier transform on the displacements components are: û = i ξ s ) ) ˆσ +s) ˆσ, û = ξ s ) ) ˆσ +s) ˆσ, 7) Whereonly thetwo eigenvalues withpositive real part Re ) ν > 0), suchthat thestresses vanish at the infinity ˆσ ij 0 for x ), has been accounted. Remembering the conditions for having positive definetess of the strain energy density introduced in section 4, the two eigenvalues with positive real part become: ν ) [ = ν ) = [ = ] ), λ ) 4n+m) ν [ ] ), λ ) 4n+im) ν = ), λ 4n m)] for < ρ ) < +, [ ), λ 4n im)] for < ρ ) <, From this form it is straightforward to note that these eigenvalues can be expressed in function of Stroh eigenvalues introduced in second section by means of the relation: ν ) = iµ ), ν) = iµ ), 73) In order to derive the weight functions, we need to evaluate explicit expressions for the singular displacements utilizing the 7). The boundary conditions along the boundary x = 0 are defined as follows: ˆσ ξ,x = 0 ) = ˆΣ ξ), ˆσ ξ,x = 0 ) = ˆΣ ξ), ξ R where Σ,Σ are the components of the singular traction defined in section 4, equations 48) and 49)). It follows that: and thus: ˆσ ξ,x = 0 ) = A ) = ˆΣ ξ), ˆσ ξ,x = 0 ) = isignξ) A ) ν) +A ) ) ν) = ˆΣ ξ), +A ) A ) = ν) ˆΣ isignξ)ˆσ ν ) ν ) = isignξ)ˆσ ν) A ) ν ) ν ) The Fourier transforms of the singular displacements fields are then: = û ξ,x ) = ξν ) ν ) + ) [ signξ) [ signξ) s ) [ s) s ) [ ν ) ] ) ) s ν ) ] ) ˆΣ +iν) ˆΣ +iν) s ) [ ˆΣ, 74), 75) s ) [ s) ν ) ν ) ] ) ˆΣ ] ) ] } ) s ˆΣ e ξ ν) x ] e ξ ν) x + 76)

23 û ξ,x ) = = ξν ) +signξ) s ) +s) 66 s) +signξ) ν ) ) [ ν ) s ) [ i ] ) ν ) [ s ) [ ν) ˆΣ ν ) ν ) ] ) ) s s) 66 ν ) ˆΣ + ] [ e ξ ν) x + i s ) [ +s) 66 s) } ] ) s s) 66 ) ν ) ν) ˆΣ ] e ξ ν) x ν ) ] ) ν ) ˆΣ + For the upper half-plane, we find the same expressions, subject to replacing ξ with ξ and the superscript ) with ) Piccolroaz et al. 009). From equations 76) and 77) and their corresponding expressions on the upper half-plane, we can derive the traces of the singular displacements transforms on the plane containing the crack: Û + ξ) = û+ ξ,x = 0 + ) = [ν +ν )s ] ),i [s +s ν ν ] ) ) ˆΣ ξ ξ Û + ξ) = û+ ξ,x = 0 + ) = i [s +s 66 ) ν +ν +ν ν )s ] ), [s ν ν ν +ν )] ) ) ˆΣ ξ ξ Û ξ) = û ξ,x = 0 [ν +ν )s ] ) ) =,i [s ) +s ν ν ] ) ˆΣ ξ ξ Û ξ) = û ξ,x = 0 ) = i [s +s 66 ) ν +ν +ν ν )s ] ), [s ) ν ν ν +ν )] ) ˆΣ ξ ξ Using the relation 73) between η, and the Stroh eigenvalues: and considering the following relations ν = iµ, ν = iµ, ν +ν = µ +µ ) = s +s 66 s, ν ν = µ µ = s s ; we deduce the expressions for singular displacements along the axes of propagation of the crack x = 0), as follows: Û + ξ) = û+ ξ,x = 0 + ) = [Imµ +µ )s ] ), i [s ) µ µ s ] ) ˆΣ ξ ξ Û + ξ) = û+ ξ,x = 0 + ) = i [s µ µ s ] ), [ Im + )] ) ) ˆΣ ξ ξ µ µ Û ξ) = û ξ,x = 0 [Imµ +µ )s ] ) ) =, i [s µ µ s ] ) ) ˆΣ ξ ξ Û ξ) = û ξ,x = 0 ) = i [s µ µ s ] ), [ Im + )] ) ) ˆΣ ξ ξ µ µ 3 77) 78)

24 These expressions can be written in the same form of the physical displacements 30) and 3): Û + ξ) = Û ξ) = ξ Y) Y ) ) ξ Y) +Y ) ) }ˆΣ 79) ξ Y) Y ) )+ ξ Y) +Y ) ) }ˆΣ 80) where the hermitian matrices Y ) and Y ) possess exactly the same form 4), evaluated by specializing the Stroh formalism to the case of a two-dimensional orthotropic material Suo, 990b): [s Y ),) Imµ +µ )] ),) i[s µ µ s ] ),) ) = [ )] i[s µ µ s ] ),) ),) s Im + 8) µ µ The Fourier transforms of the symmetric and skew-symmetric weight functions are defined respectively as the jump and the average of the singular displacements across the plane containing the crack: [Û]ξ) = Û+ ξ) Û ξ) = ξ Û ξ) = Û+ ) ξ)+û ξ) = ξ i signξ) ImY ) Y ) ) ReY ) +Y ) )}ˆΣ i signξ) ImY ) +Y ) ) ReY ) Y ) )}ˆΣ We have finally recovered the expressions 34) and 35), previously derived from the direct solution of the Riemann-Hilbert problem for the interfacial crack by means of the Stroh formalism, consequently, we can say that the two alternative formulations are perfectly equivalent, and that our result is proved by this further test. 8) 4

25 References Antipov, Y. A. 999). An exact solution of the 3-D problem of an interface semi-infinite plane crack. J. Mech. Phys. Solids, 47, Bercial-Velez, J. P., Antipov, Y. A., & Movchan, A. B. 005). High-order asymptotics and perturbation problems for 3D interfacial cracks. J. Mech. Phys. Solids, 53, 8 6. Bueckner, H. F. 985). Weight functions and fundamental fields for the penny-shaped and the half plane crack in three-space. Int. J. Solids Struct., 3, Bueckner, H. F. 989). Observations on weight functions. Eng. Anal. Bound. Elem., 6, 3 8. Dundurs, J. 969). Discussion of a paper by D. B. Bogy. J. Appl. Mech., 36, Gao, H. 99). Weight function analysis for interface cracks: mis-mach versus oscillation. J. Appl. Mech., 58, Gao, H. 99). Weight function method for interfacial cracks in anisotropic bimaterials. Int. J. Fract., 56, Gao, H., Abbudi, M., & Barnett, D. M. 99). Interfacial crack-tip field in anisotropic elastic solids. J. Mech. Phys. Solids, 40, Gautier, G., Kelders, L., Groby, J. P., Dazel, O., Ryck, L. D., & Leclaire, P. 0). Propagation of acoustic waves in a one-dimensional macroscopically inhomogeneous poroelastic material. J. Acoust. Soc. Am., 30, Gupta, V., Argon, A. S., & Suo, Z. 99). Crack deflection at an interface between two orthotropic media. J. Appl. Mech., 59, S79 S87. Hwu, C. 993). Relations among Stroh, Lekhnitskii and Muskhelishvili formulations. In Proc. of the 0th National Conference on Mechanical Engineering, pp ). Taipei, Taiwan, R.O.C:. Lazarus, V., & Leblond, J. B. 998). Three-dimensional crack-face weight functions for the semi-infinite interface crack-i: variation of the stress intensity factors due to some small perturbation of the crack front. J. Mech. Phys. Solids, 46, Lekhnitskii, S. G. 963). Theory of elasticity of an anisotropic mody. Holden-Day, San Francisco. Ma, L., & Chen, Y. 004). Weight functions for interface cracks in dissimilar anisotropic materials. Acta Mech. Sin., 0, Meade, K. P., & Keer, L. M. 984). On the problem of a pair of point forces applied to the faces of a semi-infinite plane crack. J. Elasticity, 4, 3 4. Mishuris, G. S., & Kuhn, G. 00). Asymptotic behaviour of the elastic solution near the tip of a crack situated at a nonideal interface. Z. Angew. Math. Mech., 8,

26 Piccolroaz, A., & Mishuris, G. 0). Integral identities for a semi-infinite interfacial crack in D and 3D elasticity. ArXiv:0.36v. Piccolroaz, A., Mishuris, G., & Movchan, A. B. 007). Evaluation of the Lazarus-Leblond constants in the asymptotic model for the interfacial wavy crack. J. Mech. Phys. Solids, 55, Piccolroaz, A., Mishuris, G., & Movchan, A. B. 009). Symmetric and skew symmetric weight functions in D perturbation models for semi-infinite interfacial cracks. J. Mech. Phys. Solids, 57, Piccolroaz, A., Mishuris, G., & Movchan, A. B. 00). Perturbation of mode III interfacial cracks. Int. J. Fract., 66, 4 5. Radi, E., & Mariano, P. M. 00). Stationary straight cracks in quasicrystals. Int. J. Fract., 66, Stroh, A. N. 96). Steady state problems in anisotropic elasticity. J. Math. Phys., 4, Suo, Z. 990a). Delamination specimens for orthotropic materials. J. Appl. Mech., 57, Suo, Z. 990b). Singularities, interfaces and cracks in dissimilar anisotropic media. Proc. R. Soc. Lond. A, 47, Suo, Z., Kuo, C.-M., Barnett, D. M., & Willis, J. R. 99). Fracture mechanics for piezoelectric ceramics. J. Mech. Phys. Solids, 40, Ting, T. C. T. 996). Anisotropic elasticity: theory and applications. Oxford University Press. Willis, J. R., & Movchan, A. B. 995). Dynamic weight function for a moving crack. I. mode I loading. J. Mech. Phys. Solids, 43,