1. Introduction
The construction of the homogenized strength criteria (HPC) of composite materials attracts the attention of both researchers and engineers [
1,
2,
3,
4,
5,
6,
7,
8,
9,
10]. HSC is a criterion for the strength of composite components, written in terms of macroscopic SSS. In [
2,
3], a theoretical scheme for constructing an HSC was proposed. The restriction of the approach [
2,
3] is that it is a framework scheme that cannot be directly applied to a particular composite. We applied in [
1] the theoretical scheme [
2,
3] to fibrous composites widely used in engineering [
10,
11,
12,
13,
14,
15,
16]. The approach developed in [
1] is briefly presented below (see Formulas (1)–(6) and computation scheme in
Section 4). The calculations in [
1] are based on the effect of SSS localization [
16,
17,
18,
19,
20,
21] in a composite with densely packed fibers. The dense packing of fibers in composites is a common but still special case. In addition, calculations based on the SSS localization effect take into account only the principal terms, which gives approximate results. If one wants to have precise strength criteria suitable for general cases, it is necessary to implement the general procedure described in [
2,
3]. Below, we will explain that the procedure from [
2,
3] cannot be implemented explicitly, but it is possible to construct a HSC in the form of computer programs. We develop relevant programs and present an example of HSC.
Although the idea of HPC was formulated several decades ago, it has not been fully implemented, primarily due to computational problems. The theory of homogenization in the earliest publications [
22,
23,
24,
25] discussed both homogenization (also called microscopic, general, and global) characteristics and the relationship between local and homogenized stress–strain state. It has been stated, see [
25], for example, that the aforementioned relationship, in perspective, leads to HSC. The calculation of both homogenization characteristics and local stress–strain states in a composite involves solving the so-called periodicity cell problem. To calculate all the homogenized characteristics, it is necessary to solve six periodicity cell problems (according to the number of elements of the stress or strain tensor). When using the theory of homogenization to directly construct a fracture surface (or similar objects), it is necessary to solve six periodicity cell problems for each value of the homogenized stress/strain. The more or less accurate design of the fracture surface involves taking into account the thousands of values of the homogenized stresses/strains. In the 1960s–1980s, the power of computers was sufficient to solve (often with low accuracy) problems with periodicity cells with simple geometry. A very popular object for numerical analysis in 1960s–1980s were fiber-reinforced composites reinforced with a unidirectional system of fibers. This theme remained actual until now, see, e.g., [
26,
27]. Frequently, the numerical analysis results for the unidirectional layer were used for theoretical reasoning in order to construct the strength criteria of composites with multidirectional reinforcement.
For a composite with complex geometry, computational difficulties (both related to the power of the computer and problems of theoretical nature) increase drastically. This likely led to a “gap” between research on unidirectional fiber-reinforced composites and complex structure composites. With increasing computer power, researchers are turning their attention to composites with a complex structure. Considerable attention is being paid to the numerical modeling of textile composites, such as knitted and woven, see [
28,
29,
30,
31,
32,
33] and references in this publications. As a rule, the authors developed numerical methods for calculating a composite subject to a specific SSS. Although the construction of the HSC (the composite calculation taking into account the set of SSSs) is the logical next step, this next step has not been taken in practice. Once again, this step, if it follows the homogenization theory directly, involves a much larger calculation than that required in the solution to specific SSS.
According to the results of a bibliographic search on the Internet, it should be stated that the number of publications devoted to the use of averaging theory for calculating strength is significantly (by an order with guarantee, perhaps two) less than the number of publications devoted to the use of homogenization theory to calculate the average characteristics of a composite. From an engineering point of view, the importance of material characteristics such as Young’s modulus and tensile strength seems to be equivalent.
Recently, new papers on the HSCs have appeared, see [
34,
35,
36,
37]. The recent works are restricted via the construction of certain 2D fragments or sections of HRS failure surfaces. By using our computer program, we computed a “safety zone” and “HSC failure surface” for the planar overall stress–strain state and presented corresponding 3D pictures of the aforementioned objects.
2. Homogenization Method as Applied to Composite Reinforced by Systems of Fibers
The homogenization method [
22,
23,
24,
25] has solved the principal problem of the mechanics of composite materials—the computation of the macroscopic (homogenized) characteristics of composite material. It gave rise to new scientific directions related to homogenized characteristics, for example, the topology design theory [
38,
39,
40] and auxetic materials [
41,
42,
43,
44,
45]. In another fundamental problem—the calculation of the strength properties of a composite material—progress was faring worse than in the previously mentioned issue. Especially small progress was made in the construction of HSC for specific composites. As the theoretical scheme of the construction of the HSC of composite material has been known since the 1990s (see [
2,
3]), we propose that the limited progress in the field is related to computational difficulties. We refer to computational difficulties as both problems with computational resources (hardware and software) and computational methodology.
The problem of composite strength has continuously attracted the attention of the researchers and engineers since the first papers on the strength [
2,
3,
4,
5,
6,
7,
8] of the composite [
46,
47,
48]. Many strength criteria have been developed for composite materials. Note that in all the articles mentioned above we find the same basic scheme: typical modes of local SSS are introduced; a connection between local SSS modes and macroscopic SSS is established; and the strength criteria of the composite components are written in terms of macroscopic SSS. This means that the purpose of these articles was the construction of the HSC, even if the term “homogenized strength criterion” was not used at all.
Today, computing resources (hardware and software) are sufficient to implement the methodology from [
2,
3] with the accuracy required in engineering. So, the problem is methodology. In this work, we construct HSC for a composite reinforced with layers of orthogonal fibers. The constructed HSCs resemble certain well-known strength criteria but do not fully match any of them.
Let us consider a composite reinforced with layers of orthogonal fibers. We will consider this simple reinforcement scheme in order not to clutter up our presentation with unnecessary details. We assume that the layers are parallel to the
-plane and the fibers in the layers are parallel to
or
-axes (
Figure 1). The characteristic dimension
of the microstructure of the composite (the radius of the fibers and the distances between the fibers) is assumed to be small compared to the size of the material sample:
.
To describe this two-scale material, the “fast” (microscopic) variables
are introduced in addition to the “slow” (macroscopic) variables
[
15,
16,
17].
Figure 1b displays the periodicity cell (PC)
of the composite in the “fast” variables. Denote
and
(
and
) as the opposite lateral faces of PC.
We consider the elasticity problem of composite (it means the fibers and matrix are assumed to be elastic and the connection between the fibers and the matrix is “ideal” [
9,
10,
11]). The basis for the construction of the strength criterion is the analysis of the local SSS in the composite. A complete analysis of the SSS in the components of the fiber-reinforced composite was carried out in [
1] by using the two-scale method, see [
15,
16,
17]. The main result of [
1] is that the local stresses
in the composite, subjected to the macroscopic strains
, are computed by the formula:
where elastic constants
in the fibers,
in the matrix, and the functions
(
,
) are solutions to the boundary-value problem (4) from [
1] (often referred to PCP).
By using the local stresses
corresponding to PCP (4) from [
1], we can write (1) as the following:
After the functions
have been computed, the local stresses
in the composite are computed following (1) or (2). The effective (homogenized) constants
of the composite are given by the formulas [
15,
16,
17]
Hereafter, means volume in 3D cases and square in 2D cases.
The homogenized stress and strains are connected by the homogenized Hook’s low .
3. The Strength of the Composite
Our goal is to obtain the HSC of the composite. We call the homogenized strength criterion (HSC) the strength criterion of composite components (fibers and matrix), written in terms of homogenized strains
or homogenized stresses
. We assume that the strength criterion of the material of the fibers and the binder may be written in the form
, where
In (4), is the strength limit of the fibers and is the strength limit of the matrix.
The following condition
ensures no damage at all the points of PC
. For this reason, we call the strains {
:
} the safety zone
.
The destruction of composite starts when the condition
is satisfied and occurs at point(s)
, at which the maximum in (5) reaches values 1:
The destruction of composites is a complex multistage process. This thesis was formulated in the initial works on composite materials [
18,
19] and is still accepted, see, for example, [
20,
21,
47,
48]. Therefore, when speaking about the strength or fracture of a composite, it is necessary to indicate which stage of fracture is being discussed. HSC (6)
is the “first crack” criterion and
indicates the “weakest element” of the composite. The fulfillment of Equation (6)
does not mean that the composition breaks up into separate parts. It is not excluded that the composite will not even lose its load-bearing capacity. In any case, if the condition
is met, the “first cracks” appear in numerous PCs, which means that damage to the composite will be massive, and this stage should be singled out as the specific stage of the destruction process. We say that deformations {
:
} form a fracture surface
.
The condition (5) guarantees that the composite does not fail when strains are applied to the composite. For this reason, we say that deformations {: } form a safety zone .
The safety zone and the destruction surface belong to . The destruction surface is the boundary of the safety zone .
Formulas (5) and (6) can be written in terms of homogenized stresses if the homogenized strains are expressed through by using the homogenized Hook’s low and then this expression is substituted to (5) and (6).
4. Construction of HSC in the Form of a Computer Program
The construction of HSC in [
1] was based on the SSS localization effect [
34,
35,
36,
37,
38] in the high-contrast composite. In [
1], we use the SSSs localized in different domains for the different deformation modes. As a rule, this is not the case. In addition, calculations based on the SSS localization effect are approximate, since only leading terms are taken into account. To construct a universal and accurate HSC, it is necessary to implement the general procedure described above. The construction of HSC, in general or special cases, assumes the computation of the local stresses
in the PC. In the general cases, these computations can be performed only numerically. Let us describe our proposed computational scheme. By virtue of Formula (1), the construction of the GSP can be divided into two computational blocks. The first is for PCP solution. This block is responsible for the microstructure of the composite. The other is for calculating HSC value. The detailed description of the blocks is given below.
Block of micro-structural analysis: solve PCP (4) for the unit value strains , compute the local stresses in all finite elements, and save into files. |
|
Block HSC: compute the local stresses in all finite elements for the given macroscopic strains following Formula (1) and check the condition (5). |
In order to determine the failure surface, it is necessary to replace “check condition (5)” in the HSC block with “select macroscopic strains that satisfy the equation (and draw these points if one wants for a visualization of the fracture surface)”. The calculations mentioned in the first block can be performed using the ANSYS (or similar) FEM software. The calculations mentioned in the second block can be performed using a program developed for this case (the authors developed program in C). The main component of our computer program, except the file exchange procedure, is the procedure for numerical solution to Equation (6). The program also includes a graphics procedure, which can be used for illustrative purposes only. Another version of the second program computes the HSC value (5) and decides whether the homogenized stresses/strains belong to the safety zone. The integration of HSC with ANSYS (or similar) FEM software is a separate problem.
In the present configuration, the computational programs are independent. The first block program accounts for the local (material and geometrical) characteristics of the composite. It does not depend on the macroscopic strains . The second block program—HSC itself—computes the strength criterion value for prescribed macroscopic . This program is the same for any composite material. The information about the specific structure of the composite plate is input into this program in the form of the -files.
The computations mentioned in the first block are time-consuming (from a few minutes to an hour) due to numerical solution to the 3D elasticity problem. The HCS procedure is fast (a few seconds). The failure surface construction procedure is time-consuming due to the many repetitions of the fast HCS procedure.
In our calculations, we used the ANSYS FEM software to solve the linear elasticity problem. For this reason, we did not pay much attention to the convergence and stability of numerical procedures (ANSYS FEM computational procedures satisfy the necessary requirements for an accurate solution of a linear problem of elasticity theory). We adopted the FE mesh, which twice refined the changes solution to less than 5%.
CPU 3 GHz and 1 GB of computer RAM was enough to solve problems with periodicity cells in the “Microstructural Analysis Unit” in 1–10 min. The main problem in the “Block of microstructural analysis” was the generation of a periodic FE grid (the insufficiency of ANSYS in this area is known [
49]). Executing a C program in the “HSR Block” in HSR mode takes a few seconds. The execution of program C in the “build security zone” mode was time-consuming and took between 10 and 30 min. This was due to the repeated repetition of the calculations of the “VSM block” in the “Construction of a security zone” mode.
We do not discuss specific issues here, such as interaction between Windows and DOS simulators and the like, which are specific problems associated with the software used by the authors.
5. An Example. Construction of the Failure Surfaces and the Safety Zones
We construct the HSC for the fiber-reinforced material considered in [
1]. The material parameters of the fibers and the matrix used in our computations are presented in
Table 1. These elastic characteristics correspond to the carbon/epoxy composite.
The dimensions of the components of the composite are the following (
Figure 1):
,
,
,
0.1,
, and
0.45, see
Figure 1. These values are included in the non-dimensional “fast” variables
. The corresponding dimensional values are computed by multiplying by
. For carbon fibers,
varies from 5 to 20 μm (
m).
The solution to PCP and was obtained by using ANSYS FEM software. In computation, finite element SOLID185 was used and the number of finite elements was about 3000.
In our computation, we accept the von Mises strength criterion for the fibers and the binder. In this case, , where is the deviator of the local stress tensor .
In general cases, the failure surface {} and the safety zone {} are 6D objects and cannot be visualized in an appropriate way. In the case of macroscopic in-plane deformation, the problem depends on the three macroscopic strains and the corresponding objects may be visualized. For this reason, we consider the macroscopic in-plane deformation in the plane. In this case, ().
The periodicity cell
is formed of fibers
F and matrix
M:
. The strength conditions for the fibers and the matrix separately) take the following form:
where
,
, and
The equation introduces the failure surface for the fibers and introduces the failure surface for the matrix. The destruction of the fiber starts at the point(s) at which the first maximum in (8) reaches values 1. The destruction of the matrix starts at the point(s) , at which the second maximum in (8) reaches values 1.
The failure surface for the periodicity cell (i.e., for composite as a whole) is determined as the {max(,) = 1}. This definition reflects the fact that destruction may start both in the fibers and in the matrix.
The failure surface may be determined in another equivalent way. , , and are the safety zones for the fibers, matrix, and the composite as whole. The domain is determined by the inequality max ( and ) 1 is the intersection of the domains and determined by the inequalities (7). The failure surface is the boundary of the domain .
Figure 2 displays the failure surfaces
and
. The safety surfaces are colored in red/blue and pink-yellow for the best visibility of their 3D geometry. In
Figure 2a,b, the half of the failure surfaces
and
are displayed. The remaining parts of the failure surfaces are the mirror reflections of the displayed fragments with respect to the
-plane. The mark 0.01 in the figures corresponds to 1% tension/shift strains. In
Figure 2, the failure surface
and a fragment of the failure surface
are displayed. The surface
is long in the
-direction. This means that the fibers are strong against the macroscopic shift in the
-plane. This is in agreement with the “scissor-like” deformation of the fibers [
1].
The interactions of the failure surfaces with the coordinate planes provide us with additional information concerning the geometry of the failure surfaces. There are three coordinate planes:
,
, and
.
Figure 3 displays the failure surfaces
and
for zero shift (
) and arbitrary axial tensions
and
. The planar failure surface
is colored in yellow for the best visibility. The fracture surface
is entirely contained within the fracture surface
(it is also seen in
Figure 2). This means that the weakest element of the composite is matrix. At the same time,
and
are placed at a small distance one from another. This means that carbon fibers and epoxy matrix are nearly equal in strength as the components of the composite. Note that the strengths of carbon fibers and epoxy individually are very different, see
Table 1.
Figure 4a displays the failure surfaces for zero axial stress
and arbitrary axial stress
and shift
.
Figure 4c displays the failure surfaces for the zero axial stress
and the arbitrary axial stress
and shift
. The central fragments of
Figure 4a,c are enlarged. In both cases, the matrix is the weakest element of the composite. In
Figure 4, we see that the strengths of the fibers and the matrix are similar against the macroscopic axial tension and rather different against the macroscopic shift.
We demonstrate the way in which the change of characteristics of components of the composite changes the HSC. The strength limit of epoxy may change from
to
[
50].
Figure 5 displays the interactions of the failure surfaces with the coordinate planes when the epoxy strength limit is changed to
=
Pa (the other characteristics are the same). Then, the failure surface
is placed inside the failure surface
(see
Figure 5a for the axial tension mode). This means that the fibers become the weakest element of the composite. In the tension-shift case, the failure surface
goes beyond the failure surfaces
,
Figure 5b. The failure surface of the composite as a whole consists of the fragments of the failure surfaces of the fibers and the matrix, see
Figure 5c. When strains
belong to
or
in
Figure 5c, the fibers are destroyed. When the strains
belong to
or
, the matrix is destroyed.
Figure 6 displays the 3D failure surfaces for composite as a whole when
=
Pa. It consists of the fragments of the failure surface for the fibers and matrix.
Although we use the von Mises strength criterion for the components of the composite,
Figure 2,
Figure 3,
Figure 4,
Figure 5 and
Figure 6 demonstrate no similarity with the von Mises criterium at the macrolevel. The proportions of the right angle sides in
Figure 6 are 1:4 for the matrix and approximately 1:6 for the fibers. One can treat the long length of the failure surfaces along the line
as a residual form of the hydrostatic axis.
The failure surface of the composite is the intersection of the domains and (8). Usually, intersection has a non-smooth boundary, even if and are smooth surfaces. This kind of non-smoothness results from the change of the failure mode from the “matrix failure” mode to the “fiber failure” mode (or vice versa) when the macroscopic SSS changes.
Figure 3,
Figure 4 and
Figure 5 show that the fracture surface
of the fibers can also be non-smooth. Thus, the fracture surfaces of fibers as an element of a composite are not directly related to the strength criteria of fiber materials. In the example above, we used the von Mises strength criteria for the material of the fibers. The well-known von Mises fracture surface is smooth and shows no similarity with
.