A plane problem for a crack between two anisotropic semi-infinite spaces under remote tensile-shear loading is considered. In the framework of the assumption that the crack faces are free of stresses an exact analytical solution of the problem is given on basis of the complex potentials approach. This solution possesses oscillating square root singularities in stresses and in the derivatives of the displacement jumps at the crack tips. To remove these singularities a new model founded on the introduction of the shear yield zones at the crack tips is suggested. This model is appropriate for the cases where interface adhesive layer is softer than the surrounding matrixes. Under this assumption the problem is reduced to the nonhomogeneous combined Dirichlet-Riemann boundary value problem with the conditions at infinity. An exact analytical solution of this problem is presented for the case of a single yield zone. The length of this zone is found from the finiteness of the shear stress at the end point of the zone. Due to such simulation the shear stress becomes finite at any point and the normal stress possesses only square root singularity at the crack tip. Therefore, the conventional stress intensity factor of the normal stress at the crack tip is used. The numerical illustration of the obtained solution is given.

1. Introduction

Interface cracks in many cases are the main reason of composite constructions failure. Therefore, much attention has been paid to the investigation of interface fracture in the framework of the open crack model. Such model is connected with oscillating singularity, which was the most clearly investigated in the paper [1] for a crack between two isotropic materials. The contact zone model, assuming a partial closing of the crack faces and eliminating a physically unreal oscillation, has been suggested in [2] and analytically developed in the papers [35].

The problem of an interface crack between two different anisotropic materials is much more complicated and therefore less studied compared to an isotropic case. An analytical analysis of this problem in the framework of an open interface crack model was carried out in [6] and continued in [711]. The phenomenon of oscillation was confirmed in these papers for an anisotropic case. Accounting of the crack faces contact for anisotropic bimaterials has been done in papers [12, 13] in an analytically numerical way and in [14, 15] analytically.

In many cases the interface is much softer than the adhered matrixes; therefore, thin yield zones develop at the crack tips along the interface. Accounting of such zones was performed in [16] with use of Coulomb law of dry friction and in [17, 18] by means of applying Dugdale [19] model for the simulation of the mechanical fields in these zones. However to the author’s knowledge a thin yield zone based upon the Treska’s-Saint-Venant’s theory has been never used for an interface crack in neither orthotropic nor isotropic bimaterial.

In this paper an interface crack with a yield zone based upon maximum shear stress criterion is studied. An exact analytical solution for the associated mathematical model is derived. This solution is free from an oscillation and, therefore, the conventional form of stress fracture criterion can be used.

2. Formulation of the Basic Relations

The constitutive relations of elasticity for a linear anisotropic material in the absence of body forces in a fixed rectangular coordinate system can be presented in the form [20]where , are stress and displacements components; are the elastic moduli, for which hold true and Einstein’s summation convention from 1 to 3 for Latin suffixes has been used.

Substituting (1) in (2) one gets the following Lame equations:

Assuming that all fields are independent of the coordinate , the solution of (3) according to the method suggested in [21] and developed in [6] can be presented in the form:where is an arbitrary analytic function of the complex variable ; and a1, a2, a3 are eigenvalue and eigenvector components of the following system:

The elements of the 3×3 matrices Q, R and T are defined asHere and afterwards the superscript T stands for the transposed matrix. A nontrivial solution of (5) exists if is a root of the equation

Since (7) has no real roots [9] we denote the roots of this equation with positive imaginary parts as and the associated eigenvectors of (5) as (subscript here and afterwards takes the numerals 1, 2, 3). The most general real solution of (3) can be presented as [9]where , is a matrix composed of eigenvectors, is an arbitrary vector function, , and the overbar stands for the complex conjugate. Introducing the vectorand using (2), this vector can be presented in the formwhere the components of 3×3 matrix B are defined as

Further a bimaterial composed of two different anisotropic semi-infinite spaces and with mechanical properties defined by the matrices and , respectively, is considered. We assume that the vector t is continuous across the whole bimaterial interface and the part of the interface , are bounded. That is, the boundary conditions at the interface are as follows:

Satisfying (12), using for each material the presentations (8), (10), and applying the method written in detail in the paper [14] the following expressions at the interface are obtained:where is the jump of the function across the material interface, , , .

It should be noted that the vector function W(z)= is analytic in the whole plane with a cut along (-∞,∞)L. The relations (14), (15) play an important role for the formulation of various problems of linear relationship for anisotropic bimaterials with cuts at the material interfaces under the action of mechanical loadings.

The main attention in the following will be paid to the consideration of orthotropic materials with the axis of the material symmetry parallel to as the most important class of anisotropic materials. In this case the matrix G has the following structure:where all are real. It is clear that the plane and out of plane problems can be decoupled. Because of the simplicity of the out of plane problem solution our attention will be focused on the plane problem for the displacement components (). In this case similarly to the contracted notations of the anisotropic elasticity [20] the following relations for the elastic coefficients related to the ()-plane can be introduced: and the combination of each of the first two equations of (14), (15) can be presented in the form (see for details [14])where

It is worthy to note that the function is analytic in the whole plane with a cut along (-∞,∞)L.

3. A Crack at the Interface of Two Materials

Let us assume further that a crack takes place at the section [c, a] of the material interface. The half-spaces are subjected to uniformly distributed normal stress and shear stress at infinity, which do not depend on the coordinate . The crack faces are free of loading. This kind of external fields initiates plane deformation state; therefore, only the cross-section orthogonal to (Figure 1) can be considered and the relations (17), (18) with are valid. We do not pay attention for a while to the section (a, b) of the interface, which will need consideration later.

The boundary conditions for the formulated problem are as follows:

Satisfying conditions (21) and (22) with the use of (17), (18) provides the continuity of the function over the segments of the material interface and also leads to the following equation:

The conditions at infinity follow from (17) and can be written in the formwhere , , .

The solution of (23) under the condition at infinity (24) has been obtained with use of [22] in the formwhere ,. .

The stresses at the interface are found from (17), (25) as follows:

The derivative of the displacement jumps is obtained by use of the formula (18) in the form

After integrating the last relation, we obtain

4. Formulation of the Problem and Development of the Interface Crack Model Free from Oscillation

It is clearly seen from (26), (27) that the mechanical stresses and the derivative of displacement jumps are singular at the crack tips. Moreover this singularity is oscillating. To remove this oscillation different models were suggested, for example, [2, 17, 18]. In the present paper the model based upon the introduction of the shear yield zones is suggested. As it will be shown later, either both of these zones are very short or one zone is substantially shorter than the other one; therefore their mutual influence can be neglected. Taking into account this circumstance we will consider only the longer zone (a, b) for simplicity, assuming that it occurs at the right crack tip (Figure 1). If the longer zone occurs at another crack tip then this zone can be considered by simple transformation of half-spaces.

Thus the boundary conditions for the considered model can be written as follows:where is the shear yield limit of the adhesive layer.

Satisfying the interface conditions (29) and using (17) one gets (23). Additionally the first and third conditions (30) lead to the equations

Satisfaction of the boundary conditions (31) provides the analyticity of the function outside of the interval and the last relations lead to the equationwhere .

Equations (23) and (33) present the nonhomogeneous combined Dirichlet-Riemann boundary value problem. The conditions at infinity (24) are valid for this problem also. By using the results of the paper by Nahmein and Nuller [23] the general solution of the homogeneous problem corresponding to (23), (33) can be presented in the formwhereThe functionspresent the canonical solutions of the homogeneous problem (23), (33), where , , and are arbitrary real coefficients.

A particular solution of the nonhomogeneous combined Dirichlet-Riemann boundary value problem for certain right sides of (23) and (33) was analyzed in [24]. In the present case of the problem (23), (33) its particular solution can be found in the formwhere is assumed to be analytic in the whole complex plane with a cut along the -axis. It should be mentioned that satisfies (33). Substituting (37) into (33) and taking into account that on one has the following equation:where .

A solution of the Dirichlet problem (38) has the following form [25, formula (46.25)]:where and , .

Usingthe formula (39) takes the formwhere , .

The general solution of the problem (23), (33) can be found by summing the solutions (34) and (37). Arbitrary constants can be found from the condition at infinity (24) together with condition of the displacement uniqueness which due to (18) can be written in the form

Considering that for the validity of last equation the coefficient before in the expansion of at infinity should be equal to zero [24] and also one gets the following expressions for the unknown coefficientswhere , ,

Using the obtained solutions the general solution of the nonhomogeneous combined Dirichlet-Riemann boundary value problem (23), (33) can be presented in the formFrom this solution all required quantities at the material interface can be found.

5. Stress Intensity Factor and Yield Zone Length

According to (17) the stress field on the right side from the yield zone can be presented in the form:

For an arbitrary position of the point , which defines the yield zone length, the right hand side of (47) is singular for and besides is real and is pure imaginary. Therefore, for any the normal stress is finite for whilst is singular. For removing of this singularity the equationshould be valid. After some transformation this equation can be written as follows:where .

This equation should be solved with respect to and after that the position of the point can be found from (45). Usually (49) can be solved numerically and that allows finding the largest root of this equation from the interval which we denote .

The normal stress at the interval according to (17), (46) can be found in the form

Substituting the formula (46), taking into account that [24]and using Plemeli formulas [22] one gets the following expression:where the integrals and should be considered here in sense of principal value on Cauchy [25].

Consider further the stress intensity factor (SIF) of the normal stress at the point

Considering that has a square root singularity for and has the logarithmic singularity at this point we get the following formula:

This formula after some transformations takes the form

The derivative of the crack faces displacement jump at the interval (crack opening) can be found due to (18) in the form

Substituting the expression (46) one getswhere .

The crack opening can be found by using the following formula:

6. Numerical Results and Discussion

The calculations were performed for a bimaterial composed of boron-epoxy orthotropic material (upper one) with , , , and isotropic material (lower one) having the following characteristics: , . , , were chosen and different values of the external mechanical loadings were considered.

In Table 1 the yield zone lengths and the SIF are presented for and different values of . It can be seen from these results that growing of leads to increasing of and . Such conclusion completely agrees with the finding which follows from simple physical inferences.

The dependence of the same variables as in Table 1 with respect to the intensity of the shear stress at infinity is given in Table 2. It can be seen that the yield zone length essentially grows with growing of the mentioned stress; however the SIF remain slightly dependent on this parameter. It is also explainable that the applied shear stress does not greatly influence the SIF of the normal stress.

The crack opening for ,  ,   is presented in Figure 2 for different values of . Lines I, II, III, and IV correspond to , , and , respectively. It follows from these results that the crack opening is almost symmetrical for because the yield zone length is extremely small in this case and the problem is almost symmetrical. However, increasing of leads to some contortion of the curve and its deviation from the symmetrical state.

The variation of the shear stress in the right hand side of the yield area for , the same crack geometry, mechanical loading and equals to 0 (line I), (II), (III) and (IV) is presented in Figure 3. The values of for the lines I, II, III, and IV are found with use of Table 1 and are equal to 10.13mm, 10.84mm, 12.48mm, and 15.38mm, respectively. As it can be seen the considered yield model eliminated the singularities of the shear stress and the oscillation of the normal stress at the right crack tip. Thereby such approach allows getting shear stress field at the point free from singularities and transforming the oscillating singularity of the normal stress into conventional square root singularity at the point . The last circumstance gives the possibility to introduce the SIF of the normal stress in a commonly used form.

7. Conclusion

An interface crack in a bimaterial anisotropic space under tensile-shear loading at infinity is considered. Plane strain conditions are assumed and the crack faces are free of stresses. Due to complex potentials approach an exact analytical solution (25) of the formulated problem is given. The stresses and the derivatives of the displacement jumps obtained from this solution have the oscillating square root singularities at the crack tips. The new model based upon the introduction of the shear yield zones at the crack tips is suggested for removing of these singularities. This model is used under the assumption that the interface adhesive layer is softer than the surrounding matrixes. The problem is reduced to the nonhomogeneous combined Dirichlet-Riemann boundary value problem (23), (33) with the conditions at infinity (24). An exact analytical solution of this problem is presented for the case of a single yield zone. The assumption concerning consideration of a single yield zone is approved by the fact that another zone is extremely short and does not influence the longer zone. The length of yield zone is found from the finiteness of the shear stress at its end. This gives the simple transcendental equation (49) for the determination of the yield zone length.

As a result of the performed modeling the shear stress becomes restricted in the whole region and the normal stress has only square root singularity at the crack tip. It means that the stress intensity factor of the normal stress can be introduced and calculated in the conventional way. The dependencies of the mentioned stress intensity factor, yield zone length, crack opening, and the shear stress on the applied loading are illustrated in the tables and graphical forms. It is particularly shown that the stress intensity factor of the normal stress essentially depends on the normal applied stress and moderately on the external shear stress. Finally, it is worthy to mention that the results of the paper are obtained in a simple analytical form which is convenient for the engineering applications.

Data Availability

The data obtained in the published article can be used for future investigation on the fracture mechanics of composite materials related to an interface crack in anisotropic and isotropic materials. The data is available upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.