Elastic constants of LU and EF materials.
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
- crack fracture resistance
- stress-strain state (SSS)
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 . 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 (j = 1, …, J), and cracks are placed along curves . The is 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 are applied to the crack edges, which are accepted the same on its opposite edges; and tractions are applied to the boundaries of the holes and plate.
2.1 Governing equations
Let us start from the Lekhnitskii complex potentials , where and are roots with positive imaginary part of the characteristic equation :
are elastic compliances which are included in the Hooke’s law :
where are strains and are stresses.
Consider an arbitrary path , which belongs to the domain D occupied by the plate, and select a positive direction of traversal ( Figure 2 ).
Then introduce in consideration the stress vectors at the plane tangent to the curve. The normal to it is located right relative to the selected direction of traversal. The projections of stress vectors and derivatives of displacements with respect to an arc coordinate at the curve through Lekhnitskii complex potentials are determined by the formula :
where and , where ds is a differential of arc at .
The stress vectors at path are determined using the formulas (2) by the formula:
2.2 Integral equations for anisotropic bounded plate with holes and cracks
where , , s is an arc coordinate, and and are the known functions, which are determined by the following formulas:
Here, are the values of the derivatives of the displacements with respect to the arc coordinate at the boundary of the plate and holes, are the displacements discontinuity at the cracks, 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 correspond to the concentrated forces and have the form :
where in , j = 1, 2, and k = 1, 2.
Note that when the boundary is traction free, then .
Let us substitute the potentials (6) into the formulas (4) for projections of stress vectors determined at the boundaries path and . Using Plemelj-Sokhotski formula, we obtain a system of integral equations [12, 15]:
where are stress vectors at point with coordinates . , the stress vector is determined by the formula (4) accordingly through complex potentials ; is a point with coordinates , which belongs to the contour ; with and with ; ; and , where and
Using the results , we obtained that the unknown functions at the boundary of each of the holes in representation (6) are defined up to a summand , where are arbitrary constants. At numerical solution of the problem, the constants are 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 , , and the equation of the boundary hole is described in the form , .
Let us assume the representation for the displacement discontinuity at the cracks:
Let us replace the integrals with Lobatto-type quadrature formulas , 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 . Then we obtain the system of equations:
and ; is a point with coordinates , is a point with coordinates , is a point with coordinates with , and is a point with coordinates with , where .
Analogously to , we should remove three equations with and 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 .
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 applied at points . 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 with different relative distances , 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 .
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 .
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 . With , the stable growth of cracks (without jumping) will occur at ; with if ; and with if . 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 , and the SIF is 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 () when the relative distance to forces is . 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 and . The calculation results of the relative SIF are shown in Figure 4 .
On the right are shown figures from the book . 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 . 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 is 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  that the range of working lengths is most favorable with , although allowed, and , where . With the SIF increases from zero to a certain maximum (depending on ), then falls to a minimum, and then increases monotonously. With the SIF increases monotonically with the increasing length of the crack. We note that such conclusions are valid for material with a Poisson ratio .
Since at big lengths of cracks the SIF depends little on , then for the first type of samples is 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 ). The results of calculations of the SIF in a circular isotropic sample with are shown in Figure 5 . Here and further is assumed that the relative value of the SIF is equal to and the parameter value is indicated near the curves.
Figure 5 shows that the range with few changed SIF is necessary to determine in the vicinity of the lengths of the cracks with .
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.
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 with , where . The biggest range () with SIF values close to constants with an error is achieved with . 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 () with SIF values close to constants with an error is achieved with .
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, is 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 .
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.
Two identical circular holes of radius R, the centers of which are located at points , are created for stretching in a plate. It was assumed that the load was applied to the boundary of the circular holes. Using , we accept that the forces act normally on the domain and are given on it in the form
where is angle coordinate on each of the holes, is 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 .
At first, for the purpose of testing the algorithm, the calculations are performed for the case of a localized load at and , with for the upper hole and for the lower hole (thus the stretching of the plate in the direction of the OY axis is considered). The relative SIF for an isotropic material with a Poisson ratio at different crack lengths at is calculated and given in Table 2 .
|= 0.2||= 0.4||= 0.6|
The values of relative SIF obtained by another method in  for the case of stretching by concentrated forces (i.e., with ) are given in the same table. As you can see, the results for these cases were close, with the exception of the values . 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 , that is, the case is considered when forces are applied at the center of the holes and more smoothly applied efforts (at ).
|= 0.2||0.4||0.6||= 0.2||0.4||0.6||= 0.2||0.4||0.6|
The results for a circular sample with are shown ( Figure 12 ) in the same table.
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 ) and at all lengths with . 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 in 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 (Оу)|
|= 0.2||0.4||0.6||= 0.2||0.4||0.6|
Similar results for a circular sample for the same material are given in Table 5 .
|LU (Ох)||LU (Оу)|
|= 0.2||0.4||0.6||= 0.2||0.4||0.6|
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 , 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 .
Let us consider at first a half-plane , which is loaded with force (X, Y) at an arbitrary point on the boundary. The Lekhnitskii potentials for this half-plane will be 
It can be shown  that the half-plane also corresponds to the potentials, whose boundary passes through the point and is inclined at an arbitrary angle under the action of the same force. Let us consider a semicircle with the center at point of radius , 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 . This proves that the case of loading by the concentrated force of the half-plane corresponds to the potentials .
Let us now consider a bounded plate occupying the domain D. The self-balanced concentrated forces are applied to the boundary of this domain at points . Let us represent the complex potentials in the form
where . Here the coefficients are determined based on expressions for by the substitution of X and Y on and , respectively. By substituting formulas (16) into boundary conditions, we obtain the boundary problem for obtaining the introduced complex potentials at :
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 . 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 h is the height of the cut triangle and а is the half-side of the square.
The calculations are performed at ; moreover, all tops are rounded by the arcs of the circle of the radius (the shape of the sample— Figure 13 ).
Calculation of the relative SIF 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.
|Isotr.||EF (90°)||EF (0°)||LU (90°)||LU (0°)|
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.
|Isotr.||EF (0°)||EF (90°)||LU (0°)||LU (90°)|
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.
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.
are elastic compliances which are included in the Hooke’s law.
are Lekhnitskii complex potentials.
are the potentials which correspond to the concentrated forces.
are the values of the derivatives of the displacements with respect to the arc coordinate at the boundary of the plate and holes.
are limit values of displacements in the approach to the section at the left and the right relative to the selected direction.
is the stress vector at path .
are the stress intensity factors (SIFs).
, ,are the relative SIFs.
Р is the principal vector.