The mechanical behavior near crack tip for periodic parallel cracks in an orthotropic composite plate subjected to the uniformly distributed load within the cracks surface is studied. The mechanical problem is turned into the boundary value problem of partial differential equation. By using the periodicity of the hyperbolic function in the complex domain and constructing proper Westergaard stress function, the periodicity of parallel cracks can be removed. Using the complex variable function method and the undetermined coefficients method, the boundary value problem of partial differential equation can be solved with the help of boundary conditions. The analytic expressions for stress intensity factor, stress, and displacement near the crack tip of periodical parallel cracks are obtained. When the vertical distance of cracks tends to infinity, the stress intensity factor degenerates into a single central crack situation. The stress intensity factor around the crack tip of periodic parallel cracks in an orthotropic composite plate depends on the shape factor. The interaction happens between the cracks. Finally, a numerical analysis of the stress and displacement changed with the polar angle is done.
1. Introduction
Defects in the materials will cause singular stress and cracks. Cracks in the interface, in particular, are the main reason that lessens the structural strength. As we all know, in engineering practice, one certain crack is rare; Crack is always gathered there. It is difficult to deal with a body containing agminate cracks. Therefore, one simple way to model a body containing agminate cracks is to assume that the cracks are arranged in a regular pattern. For simplicity, some of the agminate cracks can be considered ideally as periodic cracks. Periodic crack is the important mechanical model to study the interaction of multiple cracks. Consequently, the research on periodic cracks problem contributes to making an intensive understanding of failure mechanism of materials. Mechanics analysis of crack tip field is very important for engineering practice.
In recent years, the problem of collinear periodic cracks was investigated by many researchers. Hwu, Hu, Guo, et al. [1–4] studied the cracks-tip field problem on collinear periodic cracks in infinite homogeneous materials by means of complex variable function method. The expressions for stress intensity factor, stress, and displacements near crack tip were derived. Erdogan, Ozturk, Chen, and Ding [5–8] studied the antiplane problem of periodical collinear cracks in functionally graded materials by using Fourier transforms method. Sih and Zuo, Hao, Gao et al., Zhao and Meguid [9–13] studied the problem of collinear periodic cracks in piezoelectric material plane. However, to the authors’ knowledge, few papers considered the solutions for the problem on periodic parallel cracks. Pak and Goloubeva [14] studied the anti-plane problem of periodic parallel cracks in piezoelectric materials by using distributed dislocation method. The stress and the electric displacement intensity factors were obtained. Using the method of conformal mapping, Hao and Wu [15, 16] considered the anti-plane problem on parallel periodical cracks of finite length starting from the interface of two half planes. The stress intensity factor was obtained. Sanada et al. [17, 18] studied the stress intensity factors for glass-fiber reinforced plastics with an infinite row of parallel cracks at low temperatures under tension in generalized plane strain condition. By using the Fourier transforms to solve a pair of dual integral equations, the expression for the stress intensity factor was obtained. Chen and Liu [19] studied the dynamic anti-plane problem for a functionally graded piezoelectric strip containing a periodic array of parallel cracks, which were perpendicular to the boundary. By using Laplace transform and Fourier transform, Wang and Mai [20] analyzed the dynamic anti-plane problem in an infinite functionally graded material containing a periodic array of parallel cracks. Using eigenfunction expansion variational method, Chen [21, 22] analyzed the infinite strip problem of periodic parallel cracks subjected to the uniformly distributed load at infinity. The stress intensity factor at the cracks tip and the T-stress were evaluated. Tong, Jiang, Lo and Cheung [23] studied the anti-plane problem of doubly periodic cracks of unequal size in piezoelectric materials. A closed form solution of stress intensity factor was obtained by using complex variable function method. Zhou, Zhang and Li [24] studied the interactions of multiple parallel symmetric and permeable finite length cracks in a piezoelectric material plane subjected to anti-plane shear stress loading by the Schmidt method. Xiao and Jiang [25, 26] used the mapping technique to obtain a closed form solution of stress intensity factor to the problem of periodic open type parallel cracks in an infinite orthotropic elastic body. Bogdanov [27] investigated the axisymmetric problem of fracture of a prestressed composite material with a periodic system of parallel coaxial normal tensile cracks, using the harmonic potential functions and the technique of Hankel integral transformations. Rizk [28, 29] studied two periodic edge cracks in an elastic infinite strip located symmetrically along the free boundaries under thermal shock.
The mechanical behavior around the crack tip of periodic parallel cracks in an orthotropic composite plate subjected to the uniformly distributed load within the cracks surface is studied. The mechanical problem is turned into the boundary value problem of partial differential equation. Hyperbolic function is a periodic function in the complex domain. By constructing proper Westergaard stress function and using the periodicity of the hyperbolic function, the problem of periodic parallel cracks subjected to the uniformly distributed load within the cracks surface is ultimately turned into algebra problem. The analytic expressions for stress intensity factor, stress, and displacement near the crack tip of periodical parallel cracks are obtained.
2. Mechanical Model
Consider an infinite linear elastic orthotropic composite plate with periodic parallel cracks of length 2a as shown in Figure 1. Suppose that the two coordinate axes are parallel to the principal directions of material elasticity. All cracks are parallel to the x-axis and equally spaced apart from a distance ω along the y-axis. The surfaces of the cracks are subjected to the uniformly distributed stress load σ.
Orthotropic plate with periodic parallel cracks.
In the plane stress condition, the compatibility equation of the two-dimensional linear elastic body is as follows [30]:
(1)∂2εx∂y2+∂2εy∂x2=∂2γxy∂x∂y.
Assume that U=U(x,y) is the stress function which is defined as
(2)σx=∂2U∂y2,σy=∂2U∂x2,τxy=-∂2U∂x∂y.
The relations between the strain and the stress are as follows:
(3)(εxεyγxy)=(b11b120b12b22000b66)(σxσyτxy),
where b11, b12, b22, and b66 are the flexibility coefficients in the principal directions of elasticity. From the elasticity theory [31], substituting (2) and (3) into compatibility equation (1), the governing equation of the plane problem in an orthotropic composite plate can be obtained as follows:
(4)b22∂4U∂x4+(2b12+b66)∂4U∂x2∂y2+b11∂4U∂y4=0.
As shown in Figure 1, the periodic parallel cracks are subjected to the uniformly distributed stress load within the cracks surface. The boundary conditions are as follows [25, 26]:
(5)|x|⟶∞,|y|⟶∞:τxy=σy=0,-a<x<a,y=nω(n=0,±1,±2⋯):τxy=0,σy=σ.
An analysis of fracture problem near cracks tip for periodic parallel cracks subjected to the uniformly distributed load within the cracks surface can be reduced to finding the solution of the boundary value problem of partial differential equations (4) and (5).
Let
(6)U=U[x+s(y-nω)+inω](n=0,±1,±2⋯).
Substituting (6) into (4), we can obtain the characteristic equation as follows:
(7)b11s4+(2b12+b66)s2+b22=0.
It is a biquadratic equation, and its discriminant is written as
(8)Δ=(2b12+b66b11)2-4b22b11.
When Δ>0, the solutions of the characteristic equation (7) are as follows [32]:
(9)s1=iβ1,s2=iβ2,s3=s1¯,s4=s2¯,
where β2>β1>0 and
(10)β12=2b12+b662b11-(2b12+b662b11)2-b22b11β22=2b12+b662b11+(2b12+b662b11)2-b22b11.
When Δ<0, the solutions of the characteristic equation (7) are as follows [32]:
(11)s1=α+iβ,s2=-α+iβ,s3=s1¯,s4=s2¯,
where β>α>0 and
(12)2α2=b22b11-2b12+b662b11,2β2=b22b11+2b12+b662b11.
Let
(13)zj=x+sj(y-nω)+inω=xj+iyj,j=1,2.
Then, for Δ>0,
(14)x1=x,x2=x,y1=β1(y-nω)+nω,y2=β2(y-nω)+nω;
for Δ<0,
(15)x1=x+α(y-nω),x2=x-α(y-nω),y1=β(y-nω)+nω,y2=β(y-nω)+nω.
Using formulas (14) and (15), the governing equation (4) can be rewritten as a generalized biharmonic equation [32]:
(16)∇12∇22U=(∂2∂x12+∂2∂y12)(∂2∂x22+∂2∂y22)U=0.
By the theory of complex variable, the real part and the imaginary part of the analytic function are the solutions of the governing equation (4). The solution for partial differential equation (4) may be chosen as
(17)U=∑j=12[ajRe(U̿j)+bjIm(U̿j)],
where aj,bj(j=1,2) are undetermined real parameters, Uj (j=1,2) is an analytic function of zj, and
(18)dU̿jdzj=Uj¯,dUj¯dzj=Uj,j=1,2.
Substituting (17) and (18) into (2), the stress expressions can be written as
(19)σx=∂2U∂y2=∑j=12[ajRe(sj2Uj)+bjIm(sj2Uj)],σy=∂2U∂x2=∑j=12[ajRe(Uj)+bjIm(Uj)],τxy=-∂2U∂x∂y=-∑j=12[ajRe(sjUj)+bjIm(sjUj)].
3. Westergaard Stress Function
In order to solve the boundary value problem of partial differential equations (4) and (5), the Westergaard stress function is selected as follows [25]:
(20)Uj=σ[sinh(2πzj/ω)sinh2(2πzj/ω)-sinh2(2π(a+nωi)/ω)-1]j=1,2,
when |x|→∞,|y|→∞:Uj=0.
When -a<x<a,y=nω(n=0,±1,±2⋯):
(21)Uj=σ[-sinh(2πx/ω)isinh2(2πa/ω)-sinh2(2πx/ω)-1].
When Δ>0, substituting (9), (19), and (21) into boundary conditions (5), we can obtain a system of linear equations in 4 unknowns about coefficients aj,bj (j=1,2):
(22)a1β1+a2β2=0,a1+a2=-1b1+b2=0,b1β1+b2β2=0.
Solving nonhomogeneous linear equations, the unique solution of equations can be derived as follows:
(23)a1=-β2β2-β1,a2=β1β2-β1,b1=0,b2=0.
When Δ<0, substituting (11), (19), and (21) into boundary conditions (5), we can obtain a system of linear equations in 4 unknowns about coefficients aj,bj (j=1,2):
(24)a1+a2=-1,b1+b2=0,a1β+a2β-b1α+b2α=0,-a1α+a2α-b1β-b2β=0.
Solving the non-homogeneous linear equations, the unique solution of equations can be derived as follows:
(25)a1=-12,a2=-12,b1=-β2α,b2=β2α.
Substituting (23) and (25) into (17), the real analytic solution U is obtained, which meets the governing equation (4) and the boundary conditions (5).
4. Stress Intensity Factor
Fully considering the loadings and geometry size of orthotropic composite plate, we introduce the following stress intensity factor:
(26)KI=limzj→a+nωi[2π(zj-(a+nωi))]1/2Uj(zj).
Substituting (20) into (26), we obtain
(27)KI=σωtanh(2πa/ω)2=Yσπa,
where Y=(ω/2πa)·tanh(2πa/ω) is called the shape factor.
Labeling KI′=σπa, KI′ is the stress intensity factor of a single central crack in an infinite linear elastic orthotropic composite plate subjected to the uniformly distributed load within the crack surface. It can be concluded from (27) that stress intensity factor around the tip of periodic parallel cracks depends on the shape factor Y. When ω→∞, Y=(ω/2πa)·tanh(2πa/ω)→1 then KI→KI′; namely, the stress intensity factor degenerates into a single central crack situation when the vertical distance of periodic parallel cracks tends to infinity. The variation curves of KI′, Y, and KI with the cracks spacing ω are given as shown in Figure 2. It can be seen from Figure 2 that the stress intensity factor and the shape factor increase rapidly with the increase of the distance between cracks and then reach a steady state, that is, when ω→∞, Y→1, and KI→KI′. In other words, the periodic parallel cracks problem degenerates into a single central crack situation when the vertical distance of periodic parallel cracks tends to infinity, and it is entirely consistent with the previous results.
Variation curves of KI′, Y, and KI with crack spacing ω.
In order to research the interaction between cracks, the variation curve of KI/KI′ with ω/a is given as shown in Figure 3. It can be seen from Figure 3, when ω/a<10, KI/KI′ increases quickly with the increase of ω/a, and the distance between cracks is the main reason that influenced the interaction between cracks; when ω/a>10, KI/KI′ increases slowly with the increase of ω/a, and KI/KI′→1; the reason is that the interaction between the cracks is small with the increase of the distance between cracks, and so, it is considered that the periodic parallel cracks problem degenerates into a single crack problem.
Interaction of periodic parallel cracks.
5. Stress Field and Displacement Field
Substituting (13) and (26) into (20), in the vicinity of the cracks tip, we can obtain
(28)Uj(zj)=KI[2π(zj-(a+nωi))]1/2,(n=0,±1,±2⋯),aszj⟶a.
Let
(29)zj-(a+nωi)=r(cosθ+sjsinθ)
so,
(30)Uj(zj)=KI[2πr(cosθ+sjsinθ)]1/2,asr⟶0,
where polar radius r is the distance to the crack tip which is the shortest distance away from the point.
Substituting (9), (23), (30), (11), (25), and (30) into (19), respectively, and according to the relationships of the stress-strain and strain-displacement, the unified analytic expressions for stress and displacement of the periodical parallel cracks tip are achieved as follows. Consider
(31)σx=KI2πrRe{-s2(cosθ+s2sinθ)1/2]s1s2s1-s2[s1(cosθ+s1sinθ)1/2-s2(cosθ+s2sinθ)1/2]}σy=KI2πrRe{1s1-s2[s2(cosθ+s1sinθ)1/2-s1(cosθ+s2sinθ)1/2]}τxy=KI2πrRe{s1s2s1-s2[1(cosθ+s2sinθ)1/2-1(cosθ+s1sinθ)1/2]}u=KIr2πRe{1s1-s2[(cosθ+s2sinθ)1∕2s2(s12b11+b12)×(cosθ+s1sinθ)1/2-s1(s22b11+b12)×(cosθ+s2sinθ)1/2]1s1-s2}v=KIr2πRe{1s1-s2[s2s1(s12b12+b22)×(cosθ+s1sinθ)1/2-s1s2(s22b12+b22)×(cosθ+s2sinθ)1/2s2s1]}.
For an orthotropic composite plate, the material parameters are as follows:
(32)E1=70GPa,E2=11.4GPav12=0.23,G12=5.3GPaa=0.01m,r=0.001m.
For the different a/ω, the variations curves of the stress and displacement with the polar angle are given as shown in Figures 4 and 5. It can be seen from Figure 4 that stress can reach maximum value and minimum value in the range of [−90° to 90°]. For the different a/ω, the angles in which stress reaches the maximum (or minimum) value are the same, and this property is very important for S-fracture criterion and Z-fracture criterion. As seen from Figure 5, the displacement increases with the increase of a/ω. For the different a/ω, the angles in which displacement reach maximum (or minimum) value are the same, and this theory is useful to study the fracture criterion.
Variation curve of the stress with the polar angle.
Variation curve of the displacement with the polar angle.
6. Conclusions
In this paper the mechanical behaviour around the periodic parallel cracks in orthotropic composite plate is studied. The mechanical problem is turned into the boundary value problem of partial differential equation. By constructing proper Westergaard stress function and using the complex variable function method, the analytic expressions for stress intensity factor, stress, and displacement are obtained with the help of boundary conditions.
The stress intensity factor around the crack tip of periodic parallel cracks depends on the shape factor Y. The interaction happens between the cracks. When ω/a is small, strong interaction between the cracks can be found. With the increase of ω/a, the mutual influence between the cracks decreases; that is, the interaction between the cracks decreases.
The analytic expressions for stress and displacement of the periodical parallel cracks tip are obtained. The variations curves of the stress and displacement with the polar angles are given. For the different a/ω, stress and displacement can reach maximum value and minimum value, and the angles in which stress and displacement reach maximum (or minimum) value are the same. This property is very important for S-fracture criterion and Z-fracture criterion.
Conflict of Interests
The authors declare no direct financial relation with any commercial entities mentioned in the paper that might lead to a conflict of interests.
Acknowledgments
The authors acknowledge the financial support provided for this work by the Science and Technology Major Project of the Ministry of Education of China (no. 208022), the Doctoral Initiating Project of the Taiyuan University of Science and Technology (20122005), and the Graduate Scientific and Technological Innovation Project of the Taiyuan University of Science and Technology (20125027).
HwuC.Collinear cracks in anisotropic bodies19915242392562-s2.0-002640199310.1007/BF00034082HuY. T.ZhaoX. H.Collinear periodic cracks in an anisotropic medium19967632072192-s2.0-2342626232GuoJ.LuZ.Analysis of stress fields for plane problem of periodic cracks in orthotropic composites20102711621662-s2.0-77950886004ZhangX. X.LiC.CuiX. C.ZhaoW. B.Analysis of mode III collinear periodic cracks-tip stress field of an infinite orthotropic plate2012446-449208020842-s2.0-8485716578010.4028/www.scientific.net/AMR.446-449.2080ErdoganF.OzturkM.Periodic cracking of functionally graded coatings1995331521792195MR1365999ZBL0899.734212-s2.0-0001215219ChenY. Z.LinX. Y.Collinear crack problem for a strip of functionally graded materials in anti-plane elasticity2006271713ChenY. Z.Anti-plane problem of periodic crack in a strip of functionally graded materials2004364501506DingS. H.LiX.Anti-plane problem of periodic interface cracks in a functionally graded coating-substrate structure2008153153622-s2.0-5894908844210.1007/s10704-008-9302-7SihG. C.ZuoJ. Z.Multiscale behavior of crack initiation and growth in piezoelectric ceramics20003421231412-s2.0-0034300391HaoT. H.Periodical collinear air containing cracks in a piezoelectric material200111231972042-s2.0-003569979110.1023/A:1013588221105HaoT. H.Multiple collinear cracks in a piezoelectric material20013850-51920192082-s2.0-003597666510.1016/S0020-7683(01)00069-5GaoC. F.HäuslerC.BalkeH.Periodic permeable interface cracks in piezoelectric materials20044123233352-s2.0-034621694110.1016/j.ijsolstr.2003.09.044ZhaoX.MeguidS. A.On the dynamic behaviour of a piezolectric laminate with multiple interfacial collinear cracks2002399247724942-s2.0-003719806910.1016/S0020-7683(02)00154-3PakY. E.GoloubevaE.Electroelastic properties of cracked piezoelectric materials under longitudinal shear19962442873032-s2.0-003040655110.1016/S0167-6636(96)00038-5HaoT. H.An exact solution of the anti-plane parallel periodical transverse crack field in a bimaterial infinite plane1991473R49R512-s2.0-0026106865HaoT. H.WuY. C.Elastic plane problem of collinear periodical rigid lines19893369799812-s2.0-0024877265SanadaK.ShindoY.UedaS.Stress intensity factors for glass-fiber reinforced plastics with an infinite row of parallel cracks at low temperatures19982831831962-s2.0-0012914404TakedaT.ShindoY.NaritaF.SanadaK.Stress intensity factors for woven glass/epoxy laminates with cracks at cryogenic temperatures20041121091322-s2.0-164230393910.1080/15376490490277312ChenJ.LiuZ.On the dynamic behavior of a functionally graded piezoelectric strip with periodic cracks vertical to the boundary20054211-12313331462-s2.0-1344427165210.1016/j.ijsolstr.2004.10.021WangB. L.MaiY. W.A periodic array of cracks in functionally graded materials subjected to transient loading2006445-63513642-s2.0-3364669524910.1016/j.ijengsci.2005.12.007ChenY. L.Stress analysis for an infinite strip weakened by periodic cracks20042511118911942-s2.0-12444260681ChenY. Z.LinX. Y.Periodic group crack problems in an infinite plate2005429-10283728502-s2.0-1224430909810.1016/j.ijsolstr.2004.09.046TongZ. H.JiangC. P.LoS. H.CheungY. K.A closed form solution to the antiplane problem of doubly periodic cracks of unequal size in piezoelectric materials20063842692862-s2.0-3034448513610.1016/j.mechmat.2005.06.008ZhouZ. G.ZhangP. W.Interactions of multiple parallel symmetric permeable mode-III cracks in a piezoelectric material plane20092847287372-s2.0-6504908467910.1016/j.euromechsol.2009.02.004XiaoJ.JiangC.Study on the problem of periodic open type parallel cracks20073922782822-s2.0-34247369834XiaoJ.JiangC.Exact solution for orthotropic materials weakened by doubly periodic cracks of unequal size under antiplane shear200922153632-s2.0-6314916234610.1016/S0894-9166(09)60090-2BogdanovV. L.Influence of initial stresses on the stressed state of a composite with a periodic system of parallel coaxial normal tensile cracks20121861113RizkA. E. F. A.An elastic strip with periodic surface cracks under thermal shock20064411-128078182-s2.0-3374734150910.1016/j.ijengsci.2006.04.011RizkA. E. F. A.Periodic array of cracks in a strip subjected to surface heating20044116-17468546962-s2.0-304253244710.1016/j.ijsolstr.2004.02.065LekhnitskiiS. G.1963Beijing, ChinaScience PressMR0180018SihG. C.1973Groningen, The NetherlandsNoordhoff International PublishingMR0449095YangW. Y.LiJ. L.ZhangX. X.2005Beijing, ChinaScience Press