Next Article in Journal
Static and Fatigue Tensile Properties of Cross-Ply Carbon-Fiber-Reinforced Epoxy-Matrix-Composite Laminates with Thin Plies
Previous Article in Journal
Physical and Chemical Foam Injection Moulding of Natural-Fibre-Reinforced Polypropylene—Assessment of Weight-Reduction Potential and Mechanical Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Construction of the Homogenized Strength Criterion for Fiber-Reinforced Composite

by
Alexander G. Kolpakov
1,2,* and
Sergei I. Rakin
2
1
SysAn, A. Nevskogo Str., 12a, 34, 630075 Novosibirsk, Russia
2
Mathematics Department, Siberian Transport University, Dusi Kovalchuk Str., 191, 630049 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
J. Compos. Sci. 2023, 7(4), 145; https://doi.org/10.3390/jcs7040145
Submission received: 24 February 2023 / Revised: 16 March 2023 / Accepted: 3 April 2023 / Published: 6 April 2023
(This article belongs to the Section Fiber Composites)

Abstract

:
In this paper, we investigate whether the strength characteristics of composite materials can be described through the predictions of the homogenization theory concerning local stresses. We establish the homogenized strength criterion (HSC) of composite materials, following the general scheme developed in the homogenization theory. Since the homogenization theory involves solving the so-called periodicity cell problem (PCP), HSC can be constructed in the form of a computer procedure only. We developed the HSC computer program and carried out numerical calculations for fiber-reinforced material. We conclude that HSC can be used to calculate safety zones and the first failure strength criteria (see detailed definitions below). We present numerically calculated safety zones and fracture surfaces for several cases.

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 O x 1 x 2 -plane and the fibers in the layers are parallel to O x 1 or O x 2 -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: ε < < 1 .
To describe this two-scale material, the “fast” (microscopic) variables y = x / ε are introduced in addition to the “slow” (macroscopic) variables x [15,16,17]. Figure 1b displays the periodicity cell (PC) P = [ 0 , h 1 ] × [ 0 , h 2 ] × [ 0 , h 3 ] of the composite in the “fast” variables. Denote Γ i = [ 0 , h j ] × [ 0 , h 3 ] and Γ i + h i e i ( i = 1 , 2 and j = 1 , 2 ) 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 σ p q l o c ( y ) in the composite, subjected to the macroscopic strains ε m n , are computed by the formula:
σ p q l o c ( y ) = ε m n ( x ) a p q k l ( y ) Z k , l m n ( y )
where elastic constants a i j k l ( y ) = a i j k l F in the fibers, a i j k l ( y ) = a i j k l M in the matrix, and the functions Z α β ν ( y ) ( α , β = 1   or   3 , ν = 0 ,   1 ) are solutions to the boundary-value problem (4) from [1] (often referred to PCP).
By using the local stresses σ p q m n ( y ) = a p q k l ( y ) Z k , l m n ( y ) corresponding to PCP (4) from [1], we can write (1) as the following:
σ p q l o c ( y ) = ε m n ( x ) σ p q m n ( y )
After the functions Z m n ( y ) have been computed, the local stresses σ p q l o c ( y ) in the composite are computed following (1) or (2). The effective (homogenized) constants A p q m n of the composite are given by the formulas [15,16,17]
A p q m n = 1 m e s P P a p q k l ( y ) Z k , l m n ( y ) d y
Hereafter, m e s P means volume in 3D cases and square in 2D cases.
The homogenized stress and strains are connected by the homogenized Hook’s low σ k l = A k l m n ε m n .

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 ε m n or homogenized stresses σ k l . We assume that the strength criterion of the material of the fibers and the binder may be written in the form f ( y , σ p q l o c ) σ * ( y ) , where
σ * ( y ) = σ F *   in   fiber σ M *   in   binder   ,   f ( y , σ p q l o c ) = f F ( σ p q l o c )   in   fiber f M ( σ p q l o c )   in   binder  
In (4), σ F * is the strength limit of the fibers and σ M * is the strength limit of the matrix.
The following condition
F ( ε m n ) = max y P f ( y , ε m n a p q k l ( y ) Z k , l m n ( y ) ) σ * ( y ) < 1
ensures no damage at all the points of PC P . For this reason, we call the strains { ε m n : F ( ε m n ) < 1 } the safety zone V .
The destruction of composite starts when the condition F ( ε m n ) = 1 is satisfied and occurs at point(s) y 0 P , at which the maximum in (5) reaches values 1:
f ( y 0 , ε m n a p q k l ( y 0 ) Z k , l m n ( y 0 ) ) σ * ( y 0 ) = 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) F ( ε m n ) = 1 is the “first crack” criterion and y 0 indicates the “weakest element” of the composite. The fulfillment of Equation (6) F ( ε m n ) = 1 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 F ( ε m n ) = 1 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 { ε m n : F ( ε m n ) = 1 } form a fracture surface S .
The condition (5) F ( ε m n ) < 1 guarantees that the composite does not fail when strains ε m n are applied to the composite. For this reason, we say that deformations { ε m n : F ( ε m n ) < 1 } form a safety zone V .
The safety zone V and the destruction surface S belong to R 3 . The destruction surface S is the boundary of the safety zone V .
Formulas (5) and (6) can be written in terms of homogenized stresses σ m n if the homogenized strains ε m n are expressed through σ m n by using the homogenized Hook’s low ε m n = A m n k l 1 σ k l 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 σ p q l o c ( y ) 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 ε m n = δ m n , compute the local stresses σ p q m n ( y ) = a p q k l ( y ) Z k , l m n ( y ) in all finite elements, and save σ p q m n ( y ) into files.
 Jcs 07 00145 i001
Block HSC: compute the local stresses σ p q l o c ( y ) in all finite elements for the given macroscopic strains ε m n 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 ε m n that satisfy the equation F ( ε m n ) = 1 (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 ε m n . The second block program—HSC itself—computes the strength criterion value for prescribed macroscopic ε m n . 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 σ p q m n ( y ) -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): h 1 = 1.1 , h 2 = 3 , h 3 = 1.1 , h = 0.1, δ = 0.1 , and R = 0.45, see Figure 1. These values are included in the non-dimensional “fast” variables y . The corresponding dimensional values are computed by multiplying by ε . For carbon fibers, ε varies from 5 to 20 μm ( 5 20 10 6 m).
The solution to PCP and σ p q m n 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, f ( σ m n ) = 3 2 s i j s i j , where s i j is the deviator of the local stress tensor σ m n .
In general cases, the failure surface { ε m n   :   F ( ε m n ) = 1 } and the safety zone { ε m n   :   F ( ε m n ) < 1 } 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 O x 1 x 3 plane. In this case, ε i 3 = 0 ( i = 1 , 2 , 3 ).
The periodicity cell P is formed of fibers F and matrix M: P = F M . The strength conditions for the fibers and the matrix separately) take the following form:
B F ( ε α β ) 1 , B M ( ε α β ) 1  
where ( ε α β ) = ( ε 11 , ε 22 , ε 12 ) , α , β = 1 , 3 , and
B F ( ε α β ) = max y F f F ( ε α β a p q k l F Z k , l α β ( y ) ) σ F * , B M ( ε α β ) = max y M f M ( ε α β a p q k l M Z k , l α β ( y ) ) σ M *
The equation B F ( ε α β ) = 1   introduces the failure surface S F for the fibers and B M ( ε α β ) = 1 introduces the failure surface S M for the matrix. The destruction of the fiber starts at the point(s) y 0 F F at which the first maximum in (8) reaches values 1. The destruction of the matrix starts at the point(s) y 0 M M , at which the second maximum in (8) reaches values 1.
The failure surface S for the periodicity cell (i.e., for composite as a whole) is determined as the { ε m n   :   max( B F ( ε α β ) , B M ( ε α β ) ) = 1}. This definition reflects the fact that destruction may start both in the fibers and in the matrix.
The failure surface S may be determined in another equivalent way. V F , V M , and V are the safety zones for the fibers, matrix, and the composite as whole. The domain V is determined by the inequality max ( B F ( ε α β ) and B M ( ε α β ) ) 1 is the intersection of the domains V F and V M determined by the inequalities (7). The failure surface S is the boundary of the domain V .
Figure 2 displays the failure surfaces S M and S F . 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 S M and S F are displayed. The remaining parts of the failure surfaces are the mirror reflections of the displayed fragments with respect to the O ε x x ε z z -plane. The mark 0.01 in the figures corresponds to 1% tension/shift strains. In Figure 2, the failure surface S M and a fragment of the failure surface S F are displayed. The surface S F is long in the O ε 13 -direction. This means that the fibers are strong against the macroscopic shift in the O x z -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: ε 11 = 0 , ε 33 = 0 , and ε 13 = 0 . Figure 3 displays the failure surfaces S M and S F for zero shift ( ε 13 = 0 ) and arbitrary axial tensions ε 33 and ε 13 . The planar failure surface S M is colored in yellow for the best visibility. The fracture surface S M is entirely contained within the fracture surface S F (it is also seen in Figure 2). This means that the weakest element of the composite is matrix. At the same time, S M and S F 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 ε 33 and arbitrary axial stress ε 11 and shift ε 13 . Figure 4c displays the failure surfaces for the zero axial stress ε 11 and the arbitrary axial stress ε 33 and shift ε 13 . 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 30 10 6 to 90 10 6 [50]. Figure 5 displays the interactions of the failure surfaces with the coordinate planes when the epoxy strength limit is changed to σ M * = 85 10 6 Pa (the other characteristics are the same). Then, the failure surface S F is placed inside the failure surface S M (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 S M goes beyond the failure surfaces S F , 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 ( ε 11 , ε 33 ) belong to A B or D A in Figure 5c, the fibers are destroyed. When the strains ( ε 11 , ε 33 ) belong to B C or C D , the matrix is destroyed.
Figure 6 displays the 3D failure surfaces for composite as a whole when σ M * = 85 10 6 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 ε 11 + ε 33 = 0 as a residual form of the hydrostatic axis.
The failure surface of the composite is the intersection of the domains V F and V M (8). Usually, intersection V F V M has a non-smooth boundary, even if S F and S M 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 S F 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 S F .

6. Conclusions

We present an implementation of HSC for fiber-reinforced material based on the homogenization theory. In general, an HSC for a fiber-reinforced composite can be developed as a computer program. The original HSC-authoring program was developed in C/C++. The program uses auxiliary data calculated using the ANSYS FEM program. The developed HSC program causes it to determine the “safety zone” and the “the first crack failure surface”, to identify the “weakest” component of the composite (fiber or matrix), as well as the weakest point in the composite.
HSC (at least in its current form) is not designed to analyze progressive damage or the complete failure of a composite. Considering the reliability of ANSYS FEM (or similar FEM software) for solving problems in the theory of elasticity, the authors consider it possible to use HSC to predict the “safety zone” for composites.

Author Contributions

A.G.K. methodology; S.I.R. computations, programming. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All the data are included into paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kolpakov, A.G.; Rakin, S.I. Homogenized strength criterion for composite reinforced with orthogonal systems of fibers. Mech. Mater. 2020, 20, 103489. [Google Scholar] [CrossRef]
  2. Hashin, Z.; Rotem, A. A fatigue failure criterion for fiber reinforced materials. J. Compos. Mater. 1973, 7, 448–464. [Google Scholar] [CrossRef] [Green Version]
  3. Narayanaswami, R. Evaluation of the tensor polynomial and Hoffman strength theories for composite materials. J. Compos. Mater. 1977, 11, 366–377. [Google Scholar] [CrossRef]
  4. Parry, T.V.; Wroski, A.S. Kinking and tensile, compressive and interlaminar shear failure in carbon-fiber-reinforced plastic beams tested in flexure. J. Mater. Sci. 1981, 16, 439–450. [Google Scholar] [CrossRef]
  5. Soden, P.D.; Leadbetter, D.; Griggs, P.R.; Eckold, G.C. The strength of a filament wound composites under biaxial loading. Composites 1978, 9, 247–250. [Google Scholar] [CrossRef]
  6. Azzi, V.D.; Tsai, S.W. Anisotropic Strength of Composites. Exp. Mech. 1965, 5, 283–288. [Google Scholar] [CrossRef]
  7. Puck, A.; Schneider, W. On failure mechanisms and failure criteria of filament-wound glass-fiber/resin composites. Plast. Polym. Technol. 1969, 37, 33–43. [Google Scholar]
  8. Hashin, Z. Failure Criteria for unidirectional fiber composites. J. Appl. Mech. 1980, 47, 329–334. [Google Scholar] [CrossRef]
  9. Kalamkarov, A.L.; Kolpakov, A.G. Analysis, Design and Optimization of Composite Structures; Wiley: Chichester, UK, 1997. [Google Scholar]
  10. Bakhvalov, N.S.; Panasenko, G.P. Homogenisation: Averaging Processes in Periodic Media; Kluwer: Dordrecht, The Netherlands, 1989. [Google Scholar]
  11. Agarwal, B.D.; Broutman, L.J.; Chandrashekhara, K. Analysis and Performance of Fiber Composites, 4th ed.; Wiley: Hoboken, NJ, USA, 2017. [Google Scholar]
  12. Dorworth, L.C.; Gardiner, G.L.; Mellema, G.M. Essentials of Advanced Composite Fabrication &Repair, 2nd ed.; Aviation Supplies & Academics Inc.: Newcastle, WA, Australia, 2019. [Google Scholar]
  13. Lomov, S.V.; Bogdanovich, A.E.; Ivanov, D.S.; Mungalov, D. A comparative study of tensile properties of non-crimp 3D orthogonal weave and multi-layer plain weave E-glass composites. Part 1: Materials, methods and principal results. Compos. Part A 2009, 40, 1134–1143. [Google Scholar] [CrossRef]
  14. Lopes, R.S.; Moreira, C.S.; Nunes, L.C.S. Modeling of an elastic matrix reinforced with two families of fibers under simple shear: A mimic of annulus fibrosus. J. Braz. Soc. Mech. Sci. Eng. 2019, 41, 385–391. [Google Scholar] [CrossRef]
  15. El Hage, C.H.; Younes, R.; Aboura, Z.; Benzeggagh, M.L.; Zoaeter, M. Analytical and numerical modeling of mechanical properties of orthogonal 3D CFRP. Compos. Sci. Technol. 2009, 69, 111–116. [Google Scholar] [CrossRef] [Green Version]
  16. Molker, H.; Wilhelmsson, D.; Gutkin, R.; Asp, L.E. Orthotropic criteria for transverse failure of non-crimp fabric-reinforced composites. J. Compos. Mater. 2015, 50, 2445–2458. [Google Scholar] [CrossRef] [Green Version]
  17. Berlyand, L.; Kolpakov, A.G.; Novikov, A. Introduction to the Network Approximation Method for Materials Modeling; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  18. Kolpakov, A.A.; Kolpakov, A.G. Capacity and Transport in Contrast Composite Structures: Asymptotic Analysis and Applications; CRC Press: Boca Raton, FL, USA, 2009. [Google Scholar]
  19. Keller, J.B.; Flaherty, J.E. Elastic behavior of composite media. Commun. Pure Appl. Math. 1973, 26, 565–580. [Google Scholar]
  20. Kolpakov, A.A. Numerical verification of existence of the energy-concentration effect in a high-contrast high-filled composite material. J. Eng. Phys. Thermophys. 2007, 80, 812–819. [Google Scholar] [CrossRef]
  21. Rakin, S.I. Numerical verification of the existence of the elastic energy localization effect for closely spaced rigid disks. J. Eng. Phys. Thermophys. 2014, 87, 246–252. [Google Scholar] [CrossRef]
  22. Sanchez-Palencia, E. Non-Homogeneous Media and Vibration Theory; Springer: Berlin/Heidelberg, Germany, 1980. [Google Scholar]
  23. Bensoussan, A.; Lions, J.-L.; Papanicolaou, G. Asymptotic Analysis for Periodic Structures; North-Holland: Amsterdam, The Netherlands, 1978. [Google Scholar]
  24. Jikov, V.V.; Kozlov, S.M.; Oleinik, O.A. Homogenization of Differential Operators and Integral Functionals; Springer: Berlin/Heidelberg, Germany, 1994. [Google Scholar]
  25. Chechkin, G.A.; Piatnitski, A.L.; Shamaev, A.S. Homogenization; AMS: Providence, RI, USA, 2007. [Google Scholar]
  26. Braun, M.; Iváñez, I.; Ariza, M. A numerical study of progressive damage in unidirectional composite materials using a 2D lattice model. Eng. Fract. Mech. 2021, 249, 107767. [Google Scholar] [CrossRef]
  27. Lei, Y.-P.; Wang, H.; Qin, Q.-H. Micromechanical properties of unidirectional composites filled with single and clustered shaped fibers. Sci. Eng. Compos. Mater. 2018, 25, 143–152. [Google Scholar] [CrossRef]
  28. Orlik, J.; Panasenko, G.; Shiryaev, V. Optimization of textile-like materials via homogenization and beam approximations SIAM. J. Multiscale Model. Simul. 2016, 14, 637–667. [Google Scholar] [CrossRef]
  29. Rana, S.; Fangueiro, R. (Eds.) Fibrous and Textile Materials for Composite Applications; Springer: Singapore, 2016. [Google Scholar]
  30. Carvelli, V.; Poggi, C. A homogenization procedure for the numerical analysis of woven fabric composites. Compos. Part A Appl. Sci. Manuf. 2001, 32, 1425–1432. [Google Scholar] [CrossRef]
  31. Huang, Z.M. The mechanical properties of composites reinforced with woven and braided fabrics. Compos. Sci. Technol. 2000, 60, 479–498. [Google Scholar] [CrossRef]
  32. Mahmoud, B.; Manseri, L.; Rogani, A.; Navarro, P.; Marguet, S.; Ferrero, J.-F.; Tawk, I. Experimental and numerical study of the damage mechanisms in hybrid unidirectional/woven composites under impact loading. Compos. Struct. 2019, 209, 606–615. [Google Scholar] [CrossRef] [Green Version]
  33. Angioni, S.; Meo, M.; Foreman, A. A comparison of homogenization methods for 2-D woven composites. Compos. Part B Eng. 2011, 42, 181–189. [Google Scholar] [CrossRef]
  34. Füssl, J.; Lackner, R. Homogenization of Strength—A Numerical Limit Analysis Approach. In ECCOMAS Multidisciplinary Jubilee Symposium. Computational Methods in Applied Sciences; Eberhardsteiner, J., Hellmich, C., Mang, H.A., Périaux, J., Eds.; Springer: Dordrecht, The Netherlands, 2009; Volume 14. [Google Scholar]
  35. Bignonnet, F.; Hassen, G.; Dormieux, L. Fourier-based strength homogenization of porous media. Comput. Mech. 2016, 58, 833–859. [Google Scholar] [CrossRef] [Green Version]
  36. Dormieux, L.; Lemarchand, E.; Kondo, D.; Brach, S. Strength criterion of porous media: Application of homogenization techniques. J. Rock Mech. Geotech. Eng. 2017, 9, 62–73. [Google Scholar] [CrossRef]
  37. Ke, L.; van der Meer, F.P. A computational homogenization framework with enhanced localization criterion for macroscopic cohesive failure in heterogeneous materials. J. Theor. Comput. Appl. Mech. March 2022. [Google Scholar] [CrossRef]
  38. Bendsoe, M.P.; Sigmund, O. Topology Optimization; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  39. Bendsoe, M.P.; Kikuchi, N. Generating optimal topologies in structural design using a homogenization method. Comput. Methods Appl. Mech. Eng. 1988, 71, 197–224. [Google Scholar] [CrossRef]
  40. Molnár, G.; Blal, N. Topology Optimization of Periodic Beam Lattices Using Cosserat Elasticity. hal-03831937. 2022. Available online: https://hal.science/hal-03831937/ (accessed on 2 March 2022).
  41. Kolpakov, A.G. On determination of the average characteristics of elastic frameworks. J. Appl. Math. Mech. 1985, 49, 739–745. [Google Scholar] [CrossRef]
  42. Almgren, R.F. An isotropic three-dimensional structure with Poisson’s ratio = −1. J. Elast. 1985, 15, 427–430. [Google Scholar]
  43. Luo, C.H.; Han, C.H.Z.; Zhang, X.; Zhang, X.G.; Ren, X.; Xie, Y.M. Design, manufacturing and applications of auxetic tubular structures: A review. Thin-Walled Struct. 2021, 163, 107682. [Google Scholar] [CrossRef]
  44. Athul, J.; Vinyas, M.; Dineshkumar, H. On the application of additive manufacturing methods for auxetic structures: A review. Adv. Manuf. 2021, 9, 342–368. [Google Scholar]
  45. Ren, X.; Das, R.; Tran, J.P.; Ngo, T.D.; Xie, Y.M. Auxetic metamaterials and structures: A review. Smart Mater. Struct. 2018, 27, 023001. [Google Scholar] [CrossRef]
  46. Annin, B.D.; Kolpakov, A.G.; Kalamkarov, A.L. Analysis of local stresses in high modulus fiber composites. In Localized Damage Computer-Aided Assessment and Control. V.2. 1990; Comput. Mechanics Publ.: Southampton, UK, 1990; pp. 131–144. [Google Scholar]
  47. Beaumont, P.W.R. The structural integrity of composite materials and long-life implementation of composite structures. Appl. Comp. Mater. 2020, 27, 449–478. [Google Scholar] [CrossRef]
  48. Thomas, D.J. Complexity of understanding the failure of aerospace composite structures. J. Fail. Anal. Prev. 2016, 16, 513–514. [Google Scholar] [CrossRef] [Green Version]
  49. Rakin, S.I. Strength Analysis of Fibers Composites with ANSYS. AIP Conf. Proc. 2022, 2647, 060040. [Google Scholar]
  50. Lubin, G. (Ed.) Handbook of Composites; Van Nostrand: New York, NY, USA, 1982. [Google Scholar]
Figure 1. Layers of fibers in “slow” variables x and PC P of composite in the “fast” variables y .
Figure 1. Layers of fibers in “slow” variables x and PC P of composite in the “fast” variables y .
Jcs 07 00145 g001
Figure 2. Failure surfaces: (a)—for the matrix, (b)—for the fibers.
Figure 2. Failure surfaces: (a)—for the matrix, (b)—for the fibers.
Jcs 07 00145 g002
Figure 3. Failure surface for matrix and for fiber for the axial tension modes.
Figure 3. Failure surface for matrix and for fiber for the axial tension modes.
Jcs 07 00145 g003
Figure 4. Failure surface for the tension-shift modes: (a)—in plane O ε 11 ε 13 , (b)—enlarged, and (c)—in plane O ε 13 ε 33 , (d)—enlarged.
Figure 4. Failure surface for the tension-shift modes: (a)—in plane O ε 11 ε 13 , (b)—enlarged, and (c)—in plane O ε 13 ε 33 , (d)—enlarged.
Jcs 07 00145 g004
Figure 5. Failure surface: (a)—for the tension modes; (b)—for tension-shift mode; (c)—for composite as a whole (the central part of Figure (b), enlarged).
Figure 5. Failure surface: (a)—for the tension modes; (b)—for tension-shift mode; (c)—for composite as a whole (the central part of Figure (b), enlarged).
Jcs 07 00145 g005
Figure 6. The failure surfaces when the epoxy strength limit σ M * = 85 10 6 Pa.
Figure 6. The failure surfaces when the epoxy strength limit σ M * = 85 10 6 Pa.
Jcs 07 00145 g006
Table 1. The material parameters of the fibers and the matrix used in our computations.
Table 1. The material parameters of the fibers and the matrix used in our computations.
Young’s Modulus GPaPoisson’s RatioStrength Limit Pa
Fibers E F = 170 GPa ν F = 0.3 σ F * = 1.5 10 9
Matrix E M = 2 GPa ν M = 0.36 σ M * = 60 10 6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kolpakov, A.G.; Rakin, S.I. Numerical Construction of the Homogenized Strength Criterion for Fiber-Reinforced Composite. J. Compos. Sci. 2023, 7, 145. https://doi.org/10.3390/jcs7040145

AMA Style

Kolpakov AG, Rakin SI. Numerical Construction of the Homogenized Strength Criterion for Fiber-Reinforced Composite. Journal of Composites Science. 2023; 7(4):145. https://doi.org/10.3390/jcs7040145

Chicago/Turabian Style

Kolpakov, Alexander G., and Sergei I. Rakin. 2023. "Numerical Construction of the Homogenized Strength Criterion for Fiber-Reinforced Composite" Journal of Composites Science 7, no. 4: 145. https://doi.org/10.3390/jcs7040145

APA Style

Kolpakov, A. G., & Rakin, S. I. (2023). Numerical Construction of the Homogenized Strength Criterion for Fiber-Reinforced Composite. Journal of Composites Science, 7(4), 145. https://doi.org/10.3390/jcs7040145

Article Metrics

Back to TopTop