Open access peer-reviewed chapter

Determination of Stresses in Composite Plates with Holes and Cracks Based on Singular Integral Equations

By Olesia Maksymovych and Adam Podhorecki

Submitted: May 8th 2019Reviewed: June 6th 2019Published: August 23rd 2019

DOI: 10.5772/intechopen.87718

Downloaded: 223


The problems of determination of stresses at crack in bounded plates with holes of different shapes under the action of concentrated forces or distributed forces at its boundary are considered. The study is performed by the singular integral equation method. They were determined based on the established interdependences between the Lekhnitskii potentials and the stress and strain. The numerical method for solving integral equations is developed based on the quadrature method for the systems of holes and cracks. The eigen solutions of the problem were taken into account in this method. The research of stresses at cracks in samples which are used in experimental studies of crack fracture resistance was performed.


  • stress intensity factors (SIF)
  • composite plates
  • holes
  • cracks
  • crack fracture resistance
  • BIEM
  • stress-strain state (SSS)

1. Introduction

The boundary integral equation method (BIEM) is widely used to study the stress-strain state (SSS) of anisotropic plates with holes [1, 2, 3] and cracks [4, 5, 6, 7, 8, 9]. The integral equations for anisotropic plates are usually determined based on the Somigliana identity. Such equations for plates with given stresses at the boundaries of the plate are hypersingular. At the same time, the same problem for isotropic plates is reduced to singular integral equations [10, 11], for which simple numerical algorithms for solving with given precision are obtained.

In [12, 13], the simple dependencies between the Lekhnitskii complex potentials and stress and strain are obtained. In a simple form based on them and the Cauchy theorem, the integral equations are written for anisotropic plates with holes [12, 13] and cracks [14, 15, 16]. We will use the established dependencies for the construction and regularization of integral equations for anisotropic plates with holes and cracks.

For conducting experimental studies of crack fracture resistance on experimental samples in relation to isotropic materials, theoretical estimates for stresses at cracks are performed.

For such materials, the stresses in samples of different shapes with cracks under the action of stretching or compressing concentrated forces are studied in detail [10]. The experimental samples for the experimental determination of the characteristics of crack fracture resistance of various types of materials are made based on performed studies. We perform similar studies for composite samples.

2. The integral representations for anisotropic plates with holes and cracks

We consider a plate which is weakened with a system hole with boundaries L1,,LJ(j = 1, , J), and cracks are placed along curves Γkk=1K. The L0is the outer boundary of plates. Assume ( Figure 1 ) that a plate is loaded with concentrated forces (Xj , Yj ), j = 1, ..., M acting at the points (aj , bj ); tractions XTYTare applied to the crack edges, which are accepted the same on its opposite edges; and tractions XLYLare applied to the boundaries of the holes and plate.

Figure 1.

Scheme of the problem.

2.1 Governing equations

Let us start from the Lekhnitskii complex potentials Φz1,Ψz2, where zj=x+sjyand sj,j=1,2are roots with positive imaginary part of the characteristic equation Δs=0[10]:



αijare elastic compliances which are included in the Hooke’s law [10]:


where εx,εy,γxyare strains and σx,σy,τxyare stresses.

Consider an arbitrary path Γ, which belongs to the domain D occupied by the plate, and select a positive direction of traversal ( Figure 2 ).

Figure 2.

q Γ → is the stress vector at plane AB .

Then introduce in consideration the stress vectors qΓat the plane tangent to the curve. The normal to it is located right relative to the selected direction of traversal. The projections XΓYΓof stress vectors qΓand derivatives of displacements uvwith respect to an arc coordinate at the curve through Lekhnitskii complex potentials are determined by the formula [17]:


where u=du/ds,v=dv/dsand zj=dx/ds+sjdy/ds, where ds is a differential of arc at Γ.

The stress vectors qΓz=XΓ+iYΓat path Γare determined using the formulas (2) by the formula:


Assume that the vectors XYand uvare known at path Γ. Then based on Eqs. (2) and (3) at Γone has [12, 15]


where Δj=Δsj,j=1,2.

2.2 Integral equations for anisotropic bounded plate with holes and cracks

Let us write a general solution of the problem based on [12, 15] through the Lekhnitskii potentials in the form


where L=L0+L1++LJ, Γ=Γ1+Γ2++ΓK, s is an arc coordinate, and ΦΔz1and ΨΔz2are the known functions, which are determined by the following formulas:


Here, u,vare the values of the derivatives of the displacements with respect to the arc coordinate at the boundary of the plate and holes, g1=u+u,g2=v+vare the displacements discontinuity at the cracks, u±,v±are limit values of displacements in the approach to the section at the left and the right relative to the selected direction, and the potentials ΦS,ΨScorrespond to the concentrated forces and have the form [12]:


where in zkj=aj+skbj, j = 1, 2, and k = 1, 2.

Note that when the boundary is traction free, then ΦΔ=ΨΔ=0.

Let us substitute the potentials (6) into the formulas (4) for projections of stress vectors determined at the boundaries path Land Γ. Using Plemelj-Sokhotski formula, we obtain a system of integral equations [12, 15]:


where QjZTare stress vectors qLat point Zwith coordinates xy. LΓ, the stress vector is determined by the formula (4) accordingly through complex potentials Φjz1t1,Ψjz2t2,j=1,2; Tis a point with coordinates ξη, which belongs to the contour LΓ; QZ=QLZQSZQΔZwith ZLand QZ=QTZQSZQΔZwith ZΓ; QL=XL+iYL; and Qm=Xm+iYm, where Xm=2Res1Φmz1z1+s2Ψmz2z2,Ym=2ReΦmz1z1+Ψmz2z2and


Using the results [12], we obtained that the unknown functions u,vat the boundary of each of the holes Lj,j=0,1,,Jin representation (6) are defined up to a summand u˜=ωjdy/ds,v˜=ωjdx/ds, where ωjare arbitrary constants. At numerical solution of the problem, the constants ωj,j=0,,Jare to be necessarily fixed. In addition, to ensure the displacement continuity condition, it is necessary to impose the following conditions on unknown functions:


Let us consider a problem-solving equation (10) for the case of one hole and a crack. Let us assume that the contour on which the crack is placed is described parametrically in the form x=αΓτ, y=βΓτ,1τ1, and the equation of the boundary hole is described in the form x=αLθ, y=βLθ,0θ<2π.

Let us assume the representation for the displacement discontinuity at the cracks:


Let us replace the integrals with Lobatto-type quadrature formulas [15], and the integrals at the boundaries of the holes replaced by the quadrature of a rectangle, which, for periodic functions, are Gauss quadrature-type formulas [12]. Then we obtain the system of equations:


where qνkj=QjZνTkL,pνmj=QjZνTmΓ,qν=QLZνQSZνQΔZν,

UkL=uxkLykLsk,VkL=vxkLykLsk, UmΓ=g1xmΓymΓsm,VmΓ=g2xmΓymΓsm,


xmT=αΓτm,ymT=βΓτm,τm=cosπNm1,m=1,,NΓ;Cm=πNatm1and mNΓ;C1=CNΓ=0,5πN; πN=πNΓ1;TkLis a point with coordinates xkLykL, TmΓis a point with coordinates xmΓymΓ, Zνis a point with coordinates x˜νLy˜νLwith 1νNL, and is a point with coordinates x˜νNLΓy˜νNLΓwith NL<νNL+NΓ1, where x˜kL=αLθ˜k,y˜kL=βLθ˜k,θ˜n=θnH/2,x˜mΓ=αΓτ˜m,y˜mΓ=βΓτ˜m,τ˜m=cosπNm05.

We obtain the additional equation of system (12) from condition (11)


Analogously to [12], we should remove three equations with 1νNLand add the following three equations to the received incomplete system:


The first two equations follow from the displacements continuity conditions (9). The last equation is obtained when fixing an arbitrary constant ω.

The system of Eqs. (12)(14) is generalized in the case of hole and crack system in the same way as it was done in [12].

3. Stresses in circular samples with cracks under the action of concentrated forces

Let us consider the circular composite plate with radius a, which is weakened by a central crack with the half-length L. The plate is stretched by the concentrated forces ±Papplied at points 0±y0. Destruction of the plate is happening when the stress intensity factors (SIFs) reach a certain limit value. Therefore, when we calculated the limiting loads, we considered the SIF which explicitly takes into account the length of the cracks. Due to this, we performed calculations of the relative SIF Ka=KIaPπwith different relative distances α=y0/a, depending on the half-length of the crack, which is divided into a.

Calculations are made for the composite plates with elastic constants shown in Table 1 .


Table 1.

Elastic constants of LU and EF materials.

The results of the calculations for the plate made of an EF material (with a small degree of anisotropy) with the maximum stiffness in the direction of the OX axis are shown in Figure 3 .

Figure 3.

Relative SIF for a circular plate made of an EF material: a direction with maximum stiffness is parallel to the OX axis.

In general, the character of the distribution for an EF material is not significantly different from that of an isotropic material. It is necessary to increase monotonically the load for a stable growth of the crack when the distances of forces to a crack are smaller than 0,2a. With y0/a=0,3, the stable growth of cracks (without jumping) will occur at L/a>0,2; with y0/a=0,4if L/a>0,25; and with y0/a=0,5if L/a>0,35. At greater distances to forces, after reaching the corresponding level of values of traction, the circle fractures. For a case where the crack is perpendicular to the direction with the maximum stiffness of the material, the SIF is slightly increasing, especially at greater distances to forces.

For a weakly anisotropic material, the incline of the crack to the main axis of orthotropy had little effect on the SIF KI, and the SIF KIIis practically absent.

The calculations have shown that for the case of placing the crack in parallel to the direction with the maximum stiffness of the material, the above set of specifics of the SIF remain unchanged for substantially anisotropic LU-1 material. When the crack is placed perpendicular to the direction with maximum stiffness, stable crack growth occurs of small cracks (L/a0,1) when the relative distance to forces is y0/a<0,2. Moreover, for these cases, the fracture is spasmodic. In all other cases, the SIF increases monotonically with increasing crack length. Hence, the plate fractures completely after reaching forces of critical value.

Testing the developed algorithm is conducted for the case of isotropic plate with ν=0,4488and y0/a=0,15;0,16;0,18;0,2. The calculation results of the relative SIF Ka=KIaPare shown in Figure 4 .

Figure 4.

Relative SIF for isotropic material.

On the right are shown figures from the book [10]. Such calculations were also performed for the same material, and here the corresponding relative SIF is represented by dashed lines. It is seen that the results obtained by different methods coincide.

4. Determination of working intervals of crack lengths at circular samples

Two types of samples are used in experimental studies of crack fracture resistance [10]. The first is a sample for which the SIF grows monotonically with the growth of the crack. In the second type, the range of cracks’ lengths is selected in such a way that the SIF KIis practically constant. Hence, it is varied in this range from the mean value to the small value (∼2–4%). This range of crack lengths is called working. The samples of the second type are particularly suitable for conducting experimental studies including a wide range of problems in the area of destruction. In particular, with the constant force factor (the SIF is constant) in such samples, the possibility of an effective study of the rate of growth of fatigue cracks with cyclic loads, the study of crack fracture resistance depending on the influence of working environment, etc. arises.

Based on the studies in the literature for isotropic material, it has been established [10] that the range of working lengths is most favorable with α=0,18, although allowed, and α=0,16;0,2, where α=y0/a. With 0<α<0,2237the SIF KIincreases from zero to a certain maximum (depending on α), then falls to a minimum, and then increases monotonously. With α>0,2237the SIF increases monotonically with the increasing length of the crack. We note that such conclusions are valid for material with a Poisson ratio ν=0,4488.

Since at big lengths of cracks L/a>0,5the SIF depends little on α, then for the first type of samples α=0,65is taken.

Based on the obtained results, let us perform a similar study of samples of two types of composite materials.

For samples of the second type, we perform calculations only for small ratios (at α0,2). The results of calculations of the SIF in a circular isotropic sample with ν=1/3are shown in Figure 5 . Here and further is assumed that the relative value of the SIF is equal to Ka=KIaPπand the parameter value αis indicated near the curves.

Figure 5.

Relative SIF for isotropic material.

Figure 5 shows that the range with few changed SIF is necessary to determine in the vicinity of the lengths of the cracks with L/a0,25.

Similar results for the plate made of an EF material are shown in Figures 6 and 7 . Here two cases are considered: the crack is parallel or perpendicular to the direction in which the stiffness of the material is maximal.

Figure 6.

Relative SIF for an EF material: the crack is parallel to the direction of greater stiffness of the material.

Figure 7.

Relative SIF for an EF material: the crack is perpendicular to the direction of greater stiffness of the material.

From the given data, it is seen that the SIF for an EF material is bigger in the case when the crack is perpendicular to the direction of maximum stiffness of the material.

The following conclusions are made based on the data of the calculations: for the crack that is parallel to the direction of maximum stiffness of the material, the minimum deviations from the constant SIF (with an error of not more than 2%) are achieved on the ranges of lengths of cracks with relative dimensions with Δ<0,21with α=0,240,26, where Δ=L2L1/a. The biggest range (Δ=0,27) with SIF values close to constants with an error 2,1%is achieved with α=0,24. For the crack that is perpendicular to the direction of maximum stiffness of the material, the range with SIF values close to constants is reduced. Moreover, the distance of force application needs to be increased. The biggest range (Δ=0,27) with SIF values close to constants with an error 3,2%is achieved with α=0,24.

Let us consider the case of samples of the first type 1. For them, the forces are selected that are distant from the crack. For isotropic materials, as a rule, α=0,65is taken. The above-obtained results of the calculations show that the same distance can be chosen for the composite materials with a crack parallel to the principal axes of the orthotropy.

To compare the effect of a sample shape, similar calculations are made for a square plate with a crack. The results of calculations for such an isotropic sample that are similar to the results of calculations for a circle are shown in Figure 8 .

Figure 8.

Relative SIF for a square isotropic sample.

The conclusion is made based on the comparison of Figures 5 and 8 that with small distances of forces from cracks, the shape of the sample has little effect on the SIF.

Similar results of calculations for LU material are shown in Figure 9 for a horizontal crack and in Figure 10 for a diagonal crack.

Figure 9.

LU material, square sample, horizontal crack.

Figure 10.

LU material, square sample, diagonal crack.

5. Determination of the SSS of samples under the action of the tractions applied to the hole’s boundary

Let us apply the developed algorithm to study a square plate with a half-side a; weakened by a central crack with a half-length L ( Figure 11 ), the edges are not loaded.

Figure 11.

Scheme of the sample.

Two identical circular holes of radius R, the centers of which are located at points 0±c, are created for stretching in a plate. It was assumed that the load was applied to the boundary of the circular holes. Using [10], we accept that the forces act normally on the domain θθc<γand are given on it in the form


where θis angle coordinate on each of the holes, θcis angle coordinate of the middle of the domain, and Р is the principal vector applied to the domain of forces, which is directed from the center of the hole at an angle θc.

At first, for the purpose of testing the algorithm, the calculations are performed for the case of a localized load at γ=π/32and α=c/a=0,2;0,4;0,6, with θc=π/2for the upper hole and θc=π/2for the lower hole (thus the stretching of the plate in the direction of the OY axis is considered). The relative SIF Y=KIPπLfor an isotropic material with a Poisson ratio ν=1/3at different crack lengths at R/a=0,1is calculated and given in Table 2 .

α= 0.2α= 0.4α= 0.6

Table 2.

Relative SIFs for a square sample, isotropy.

The values of relative SIF YSobtained by another method in [10] for the case of stretching by concentrated forces (i.e., with γ0) are given in the same table. As you can see, the results for these cases were close, with the exception of the values L/a=0,8. Some differences in them are due to different localizations of applied tractions.

The results of the calculations of SIF are given in Table 3 for the case R0,γ=0, that is, the case is considered when forces are applied at the center of the holes and more smoothly applied efforts (at γ=π/8).


Table 3.

Relative SIFs for a square and circular sample, isotropy.

The results for a circular sample with γ=π/8are shown ( Figure 12 ) in the same table.

Figure 12.

Circular sample with holes and a crack.

The following conclusions are made based on Tables 2 and 3 : SIF does not differ significantly in the case of distributed loads with different degrees of localization, SIF increases somewhat with the growth of the domain of action of tractions, and SIFs are bigger at small crack lengths at point action of tractions (at R0) and at all lengths with α>0,6. SIFs in a circular sample are bigger than in a square one under the same load conditions.

Similar results for a square sample made from a LU material are given in Table 4 . Here, the relative SIFs Kain which the crack length is explicitly taken into account and the case where the crack is parallel to the direction of bigger stiffness (data on the left) and is perpendicular to it (data on the right) are given.

LU (Ох)LU (Оу)

Table 4.

Relative SIFs Kafor square sample, LU material.

Similar results for a circular sample for the same material are given in Table 5 .

LU (Ох)LU (Оу)

Table 5.

Relative SIFs for a circular sample, LU material.

6. Determination of the SSS of samples loaded with concentrated forces at the boundary (compression test)

The methods of studied crack fracture resistance based on sample compression in experimental practice are widely used.

The direct application of the abovementioned variant of the method of integral equations for the case when the concentrated forces act on the boundary of the domain is associated with significant errors, because unknown functions in the vicinity of the points of application of forces have a singularity. Due to this, it is necessary to separate a singular part in the solution for a more precise solution of this type of task.

6.1 Determining a singular part of the solution of this problem

Let us consider a point of the plate boundary z0, in which the concentrated force (X, Y) is applied. A singular part of the solution (of Lekhnitskii potentials) will be the same as in the half-plane, whose boundary is tangent to the plate at the point of action of the concentrated force. Let us mark this angle through φand potentials for the half-plane through Φ0z1,Ψ0z2.

Let us consider at first a half-plane y<0, which is loaded with force (X, Y) at an arbitrary point x0on the boundary. The Lekhnitskii potentials for this half-plane will be [17]




It can be shown [17] that the half-plane also corresponds to the Φ0,Ψ0potentials, whose boundary passes through the point x00and is inclined at an arbitrary angle under the action of the same force. Let us consider a semicircle with the center at point x00of radius ρ0, which belongs to the half-plane. It is easy to show that the principal vector of all forces applied to the arc of the semicircle is equal to XY. This proves that the case of loading by the concentrated force of the half-plane corresponds to the potentials Φ0,Ψ0.

Let us now consider a bounded plate occupying the domain D. The self-balanced concentrated forces XjYj,j=1J.are applied to the boundary of this domain at points zj=xj+iyj. Let us represent the complex potentials in the form


where z1j=xj+s1yj,z2j=xj+s2yj. Here the coefficients Aj,Bjare determined based on expressions for A,Bby the substitution of X and Y on Xjand Yj, respectively. By substituting formulas (16) into boundary conditions, we obtain the boundary problem for obtaining the introduced complex potentials at xyL:


where L is the boundary of domain D:


It is easy to show that the right-hand side of formula (17) is a continuous and limited function, and therefore the introduced complex potentials with an index Δare continuous and limited in the vicinity of the points of application of forces. In this regard, the above-developed numerical algorithm based on BIEM can be used to determine these potentials.

6.2 Calculation of the SIF for a rectangular sample with compression

Let us consider a square plate with half-side a, which contains a diagonal central vertical crack with half-length L. The crack fracture resistance of such a sample is determined based on compression by force R applied in vertical direction. The tips of the angles, in the vicinity of which forces are applied, can be cut off. Due to this, a sample is considered whose tops have coordinates


where in c=2a,zp=hih,zm=hih,h is the height of the cut triangle and а is the half-side of the square.

The calculations are performed at h=c/8; moreover, all tops are rounded by the arcs of the circle of the radius a/10(the shape of the sample— Figure 13 ).

Figure 13.

Sample view.

Calculation of the relative SIF Ka=KIaPπfor isotopic material and composites of EF and LU with different directions of the orthotropic axis are given in Table 6 . The angle between the crack and the direction with the maximum stiffness of the material is indicated in brackets.

L/aIsotr.EF (90°)EF (0°)LU (90°)LU (0°)

Table 6.

Relative SIF Kawhen compressing a sample with a diagonal crack.

The table shows a significant effect on the SIF of the placement of the crack relative to the axis with the maximum stiffness of the material. In particular, for cracks parallel to the maximum stiffness direction, the SIFs appeared to be significantly larger than those returned for 90°. The difference between the SIFs for these two directions is increasing for a substantially anisotropic LU material.

Table 7 shows the results of calculations for the case of stretching the same sample with a horizontal crack ( Figure 14 ).

L/aIsotr.EF (0°)EF (90°)LU (0°)LU (90°)

Table 7.

Relative SIFs when stretching a sample with a diagonal crack.

Figure 14.

Sample view.

Based on the comparison of data from Tables 6 and 7 , it follows that, unlike the case of compression, the SIF with stretching is larger for cracks that are perpendicular to the direction with maximum stiffness of the material.

7. Conclusions

An algorithm for calculating stresses at cracks in bounded plate with holes of various shapes due to concentrated forces or distributed forces at its boundary has developed. The solution of integral equations is performed by quadrature Gauss-type formulas for regular and singular integrals.

The research of stresses at cracks in the samples which are used in experimental studies of crack fracture resistance was performed.

The calculation of the stresses at cracks in samples of various forms is performed, in which ones’ experimental research are performed. To study the crack fracture resistance of composite samples, the optimal distances from the central crack to the forces at which the SIF increases monotonously with increasing crack length are determined. In particular, for square samples with a half-side a, forces should be placed at a distance of 0.6a–0.7a from the crack. For the experimental study of the growth rate of fatigue cracks, there are definite ranges of lengths of cracks for which SIFs are practically constant values. At the same time, the distances are determined at which it is expedient to apply forces. The problem for studying samples with cracks with compression is considered.

List of symbols with explanations

s is an arc coordinate.

νis a Poisson ratio.

εx,εy,γxyare strains.

σx,σy,τxyare stresses.

αijare elastic compliances which are included in the Hooke’s law.

Φz1,Ψz2are Lekhnitskii complex potentials.

ΦS,ΨSare the potentials which correspond to the concentrated forces.

u,vare the values of the derivatives of the displacements with respect to the arc coordinate at the boundary of the plate and holes.

u±,v±are limit values of displacements in the approach to the section at the left and the right relative to the selected direction.

qΓis the stress vector qΓz=XΓ+iYΓat path Γ.

KI,KIIare the stress intensity factors (SIFs).

Ka, Y,YSare the relative SIFs.

Р is the principal vector.

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Olesia Maksymovych and Adam Podhorecki (August 23rd 2019). Determination of Stresses in Composite Plates with Holes and Cracks Based on Singular Integral Equations, Dynamical Systems Theory, Jan Awrejcewicz and Dariusz Grzelczyk, IntechOpen, DOI: 10.5772/intechopen.87718. Available from:

chapter statistics

223total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Distinctive Characteristics of Cosserat Plate Free Vibrations

By Lev Steinberg and Roman Kvasov

Related Book

First chapter

Applications of 2D Padé Approximants in Nonlinear Shell Theory: Stability Calculation and Experimental Justification

By Igor Andrianov, Jan Awrejcewicz and Victor Olevs’kyy

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us