Next Article in Journal
Reflecting the Quality Degradation of Engine Oil by the Thermal Diffusivity: Radiative and Nonradiative Analyses
Previous Article in Journal
Design and Optimization of Thermal Field for PVT Method 8-Inch SiC Crystal Growth
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

True Stress Theory of Matrix in A Composite: A Topical Review

School of Aerospace Engineering and Applied Mechanics, Key Laboratory of the Ministry of Education for Advanced Civil Engineering Materials, Tongji University, 1239 Sipping Road, Shanghai 200092, China
Materials 2023, 16(2), 774; https://doi.org/10.3390/ma16020774
Submission received: 29 November 2022 / Revised: 3 January 2023 / Accepted: 9 January 2023 / Published: 12 January 2023
(This article belongs to the Section Advanced Composites)

Abstract

:

Highlights

Non-uniform matrix stresses in composites cannot be assessed against matrix strengths obtained in uniform stress states.
The non-uniform matrix stresses must be converted into equivalent uniform ones called true stresses before estimating any matrix-dominated composite failure.
Conversion is achieved through matrix true stress theory summarized in this paper.
Selected applications for resolving several challenging problems in composite failures are highlighted.

Abstract

Whereas mechanics theories for isotropic materials are almost matured, only linear elastic theories for composites were essentially established. This is because only homogenized or approximated stresses are obtainable for a composite. Its mechanical properties must be estimated on a true stress level. According to Eshelby, the true stresses of the fiber are the same as its homogenized counterparts. The true stress theory for the matrix was systematically established by the author, and is reviewed and summarized in the paper. An Excel table-based program for calculating all of the possible true stress components is provided as a supplement for the reader to download. As most composite failures are caused by matrix failures, the true stress theory plays a predominant role in estimating the composite properties outside a linear elastic range. Some challenging composite failures were resolved upon the matrix true stresses, and are highlighted in the paper.

Graphical Abstract

1. Introduction

The objectives investigated by solid mechanics can be classified into isotropic and anisotropic. The mechanics theories of isotropic materials are almost matured. Namely, almost every mechanics problem of an isotropic material is able to be well dealt with using existing theories. In contrast, the matured mechanics theories for anisotropic or composite materials are limited to linear elasticity only. All of the other mechanics problems of the composites, such as nonlinear constitutive relations or deformations, and failure and strength predictions, are overall not able to be resolved satisfactorily.
A fundamental reason is that only homogenized stresses are obtainable for a composite by the current theories, and the resulting internal stresses of the fiber and matrix are the homogenized values as well. They must be converted into true stresses before a composite property can be analyzed, as is the case for an elastic property of the composite as well. Only because the homogenized and the true stresses are in a linear elastic range are the resulting composite elastic properties predicted on both of them the same, giving rise to a false impression that an elastic property of the composite can be independent of its true stresses. A schematic diagram to illustrate this concept is shown in Figure 1.
The stress field of the fiber (standing for the reinforcement subsequently) in a continuous fiber composite is uniform [1,2,3]. Its true and homogenized stresses are the same. The true stresses of the fiber in a short fiber/particle composite are immaterial, since the fiber/particle is generally linear elastic up to rupture and less possible to fail before a matrix failure. Furthermore, a short fiber can be sufficiently well considered as an ellipsoid of a finite length. According to Eshelby [1], the stress field in the ellipsoid is uniform no matter what kind of load is applied to the composite. Thus, the bottleneck for advancing mechanics of composites is to evaluate the true stresses of the matrices in the composites.
We found a true stress of the matrix is determined by multiplying its homogenized counterpart with a stress concentration factor (SCF) of the matrix in the composite [4,5]. A plate with a hole generates a stress concentration at the hole edge. When the hole is filled with a fiber of different properties, a stress concentration occurs as well. However, such an SCF cannot be defined through a classical method any more. Otherwise, the resulting SCF would be infinite if an interface crack between the fiber and matrix occurred, as the matrix stress field at the crack tip is singular. Even though the interface is in perfect bonding up to failure, the SCF given by the classical method is still less pertinent since a predicted composite failure incorporated with it is still far away from a measured counterpart. We succeeded in achieving the new definition [4,5], and almost all of the possible SCF formulae were derived [4,5,6,7,8,9,10,11]. The matrix true stresses were then determined accordingly.
It is evident that the load share of the fiber attains the highest when a continuous fiber composite is subjected to a longitudinal load, and thus a fiber fracture can occur most possibly. Even in such a case, however, a fiber splitting failure is often seen at a longitudinal tension (Figure 2a), and a fiber kinking failure occurs most frequently at a longitudinal compression (Figure 2b [12]), both of which are caused by a matrix shear failure [13,14]. Due to an inevitable fiber misalignment in a composite fabrication, a longitudinal load generates a shear stress component in the misaligned coordinate system (Figure 2c). The axial stress component will cause the fiber to possibly fail, whereas the shear component can bring a failure to the matrix. Apparently, both the splitting and kinking failure modes indicate that the matrix shear failures occur earlier than a fiber fracture. Otherwise, the composite failure surface would be perpendicular to the loading direction.
The above evidences demonstrate that almost all of the composite failures are essentially caused by matrix failures. Without the true stress theory, none of the matrix failures can be analyzed properly. An Excel table-based program to calculate almost all of the possible matrix SCFs derived by the author et al. [4,5,6,7,8,9,10,11] was worked out, and is attached as a supplement with this paper for a reader to download (Supplementary Materials). One can automatically obtain the desired results by running the program with input of the monolithic fiber and matrix property data and the fiber geometric parameters. For a continuous fiber composite, only the fiber volume fraction is required. For a discontinuous fiber composite, two other geometric parameters, i.e., the fiber aspect ratio and the fiber length ratio, should be provided as well. Using the true stress concept, a number of long-standing and challenging problems in composite failures are resolved satisfactorily and easily. They include a necessary and sufficient condition for a fiber and matrix interface debonding to occur given a composite subjected to an arbitrary load [8], a fiber splitting [13] or kinking [14] failure versus a fiber fracture, and a large shear deformation of the composite induced from relative slippage displacements between debonded fiber and matrix interfaces [15]. All of these failures are analyzed with no iteration. All of these failures are analyzed using fiber and matrix properties plus limited composite data measured independently and following existing standards. The matrix SCF formulae, together with their applications to the analyses of the challenging composite failures, are briefly reviewed and summarized in the paper.

2. Homoginized Stresses

In mechanics of continuum media, a stress at a media point is defined as the average of those on an infinitesimally small unit element containing the point through
σ i = 1 V V σ ˜ i d V , V 0     ( i = 11 , 22 , 33 , 23 , 13 , 12 )
where σ ˜ i is called a point-wise stress. For a composite, however, its unit element (called a representative volume element or RVE) cannot be infinitesimal, since both the fiber and matrix must be contained in it. Equation (1) becomes
σ i = 1 V V σ ˜ i d V = V f V ( 1 V f V f σ ˜ i d V ) + V m V ( 1 V m V m σ ˜ i d V ) = V f σ i f + V m σ i m
where V′ = V f + V m . V f and V m are the fiber and matrix volumes in the RVE, respectively, and V is a volume fraction. Hence, a stress or strain in a composite is a homogenized or approximated quantity by nature. No exact value of it is available, although, for a constituent, e.g., matrix, both the point-wise and the homogenized stresses, σ ˜ i m and σ i m , are able to be determined.
Mechanics is a rigorous science. Only an infinitesimally small element is mathematically rigorous enough to denote a material point, regardless of any geometric shape for the element. Although an RVE for a composite cannot be infinitesimal, its purpose is to define a stress for the composite. It is equivalent to an infinitesimally small element for defining a stress in an isotropic solid. Therefore, an RVE cannot be literally considered as a representative of the entire composite, since at different points the composite microstructure and the stresses can be different.
From the very fundamental definition, it is seen that the smaller the chosen RVE is, the more accurate the homogenized stress will be in general. Indeed, this was verified through a comparison study with experiments in [16]. One of the smallest RVEs for a continuous fiber composite is shown in Figure 3a, whereas that for a short fiber composite is indicated in Figure 3b, both of which are transversely isotropic and resemble a unidirectional (UD) and a uniaxially aligned short fiber composite, respectively. It is noticed that when the fiber aspect ratio, l/a, is close to 1, Figure 3b becomes an RVE of a particle composite.
Using a bridging equation, { σ i m } = [ a i j ] { σ j f } , to correlate the homogenized stresses of the matrix with those of the fiber, one obtains [17]
{ σ i f } = ( V f [ I ] + V m [ a i j ] ) 1 { σ j }
{ σ i m } = [ a i j ] ( V f [ I ] + V m [ a i j ] ) 1 { σ j }
[ S i j ] = ( V f [ S i j f ] + V m [ S i j m ] [ a i j ] ) ( V f [ I ] + V m [ a i j ] ) 1
Equation (5) implies that the compliance tensor [Sij] in a constitutive relation, {εi} = [Sij]{σj}, of the composite is equivalent to its bridging tensor [aij], in which [ S i j f ] and [ S i j m ] are, respectively, the compliance tensors of the fiber and matrix known in advance.
There are a number of micromechanics models, either analytical or numerical, which can be used to determine an elastic bridging tensor of the composite. The bridging model established by this author [17,18] has several unique features. First, it is the only available constitutive theory that is consistent in calculating the internal stresses [16]. Let [Sij]2D, [aij]2D, [ S i j f ] 2 D , and [ S i j m ] 2 D denote the 2D (two-dimensional) quantities deteriorated from the corresponding 3D ones, [Sij], [aij], [ S i j f ] , and [ S i j m ] . Let the composite be subjected to a planar stress state, {σj}2D = {σ11,σ22,σ12}T. The internal stresses of the fiber and matrix can be calculated from Equations (3) and (4) in either the 2D or the 3D formulae, with a 3D stress vector given by, e.g., {σj} = {σ11,σ22,0,0,0,σ12}T. If the resulting internal stresses from the two sets of formulae are exactly the same, the constitutive theory is consistent in the internal stress calculation. A necessary and sufficient condition for a constitutive theory to be consistent is that its bridging tensor is upper-triangular. Only the bridging tensor of the bridging model is always upper-triangular, whereas all of the others known yet are not [16]. For a non-consistent theory, its 3D formulae have to be used to calculate the internal stresses of the fiber and matrix more accurately. The calculation amount is at least one order higher than that by using the 2D ones. Second, it is a unified analytical constitutive model for a composite reinforced with either continuous fibers, or short fibers, or particles, no matter the constituent, e.g., matrix, undergoes an elastic or plastic deformation [17,18,19,20]. Third, it is overall more accurate than most other micromechanics including numerical models [21,22,23].
When both the fiber and matrix are in elastic deformations, the homogenized internal stresses by the bridging model are expressed as [7]
σ 11 f = σ 11 V f + V m a 11 V m a 12 ( σ 22 + σ 33 ) ( V f + V m a 11 ) ( V f + V m a 22 )
σ 11 m = a 11 σ 11 V f + V m a 11 + V f a 12 ( σ 22 + σ 33 ) ( V f + V m a 11 ) ( V f + V m a 22 )
σ i j f = σ i j V f + V m a 22 ,   i j = 22 , 33 , 23
σ i j m = a 22 σ i j V f + V m a 22 ,   i j = 22 , 33 , 23
σ i j f = σ i j V f + V m a 66 ,   i j = 12 , 13
σ i j m = a 66 σ i j V f + V m a 66 ,   i j = 12 , 13
In Equations (6)–(11), {σ11,σ22,σ33,σ23,σ13,σ12} are any loads applied on the RVE or a UD/uniaxially aligned short fiber composite. Explicit expressions for the nonzero bridging tensor elements are as follows [17,18]:
a 11 = E m / E 11 f , f o r   a   c o n t i n u o u s   f i b e r   c o m p o s i t e V f E m V m E 11 f σ 11 L E 11 f [ ε L 1 l + ε L 2 ( L l ) ] E m [ ε L 1 l + ε L 2 ( L l ) ] σ 11 L , f o r   a   s h o r t   f i b e r / p a r t i c l e   c o m p o s i t e
a 22 = a 33 = a 44 = 0.3 + 0.7 E m / E 22 f
a 55 = a 66 = 0.3 + 0.7 G m / G 12 f
a 12 = a 13 = S 12 f S 12 m S 11 f S 11 m ( a 11 a 22 ) = E 11 f ν m E m ν 12 f E 11 f E m ( a 22 a 11 )
E 11 f , E 22 f , G 12 f and ν 12 f are, respectively, the longitudinal, transverse, in-plane shear moduli and longitudinal Poisson’s ratio of the fiber. Em, Gm, and νm are the matrix tensile and shear moduli and Poisson’s ratio; ε L 1 and ε L 2 are the homogenized uniaxial strains of the regions Ω1 and Ω2 in Figure 3b, which are determined as follows [18]:
0 0 l L [ g E 11 f + ( 1 g ) E m ] ( 1 l / L ) E m 2 ν m E m ( 1 + ν m ) ( 1 2 ν m ) 2 g ν m E m ( 1 + ν m ) ( 1 2 ν m ) 0 ( 1 ν m ) E m ( 1 + ν m ) ( 1 2 ν m ) E m [ g ν 12 f + ν m ( 1 g ) ] 0.5 ( 1 + ν m ) ( 1 2 ν m ) 2 g ν 12 f E m 1 + ν m k 1 0 k 2 k 3 k 4 0 C 1 C 2 a 2 ε L 1 ε L 2 = σ 11 σ 11 σ 11 0
k 1 = 2 g ν m ν 12 f + ( 1 ν m ) ( 1 g ) ( 1 + ν m ) ( 1 2 ν m ) E m + g E 11 f , g = ( a / b ) 2
k 2 = E 11 f E m ν 12 f ( 1 + ν 23 f ) + ν 12 f E 22 f ( 1 + 2 ν m ) ν m ν 12 f E 22 f ( 1 + ν m ) ( 1 2 ν m ) [ E 11 f E 22 f ( ν 12 f ) 2 ] 2 E m ν 12 f ( 1 + ν m ) ( 1 2 ν m )
k 3 = E 11 f ν 12 f ( 1 + ν m ) E 22 f ( 1 + ν 23 f ) E m ν 12 f ( 1 + ν m ) [ E 11 f E 22 f ( ν 12 f ) 2 ] + 2 E m ν 12 f 1 + ν m
k 4 = E 11 f E 22 f ( ν 12 f ) 2 ( 1 + ν m ) ( 1 2 ν m ) + E m ν m ν 12 f ( 1 + ν 23 f ) ( 1 + ν m ) ( 1 2 ν m ) [ E 11 f E 22 f ( ν 12 f ) 2 ] 2 ν m E m ν 12 f ( 1 + ν m ) ( 1 2 ν m )
ν 23 f is the fiber transverse Poisson’s ratio. Without loss of any generality, one can set σ11 = 1 on the right hand side of Equations (12) and (16). When the fiber aspect ratio l/a → ∞ or its length ratio l/L → 1 (see Figure 3b for definition), the expression of a11 for a short fiber composite in Equation (12) is near to that for a continuous one. The five elastic moduli of the UD or uniaxially aligned short fiber composite are given by [7]
E 11 = ( V f + V m a 11 ) E 11 f E m V f E m + V m a 11 E 11 f
ν 12 = E 11 ( V f ν 12 f E m + V m a 11 ν m E 11 f ) E 11 f E m ( V f + V m a 11 )
E 22 = ( V f + V m a 11 ) ( V f + V m a 22 ) ( V f + V m a 11 ) ( V f S 22 f + a 22 V m S 22 m ) + V f V m ( S 21 m S 21 f ) a 12
G 12 = G 12 f G m ( V f + V m a 66 ) V f G m + V m G 12 f a 66
G 23 = G 23 f G m ( V f + V m a 22 ) V f G m + V m G 23 f a 22 .
It is noted that among the three geometric parameters used to define the RVE geometry of a short fiber composite (Figure 3b), Vf, l/a and l/L, only Vf and l/a are measurable in advance. The length ratio l/L is not directly measurable. An empirical expression for it is proposed in [19], which reads
l L = V f l a tan arctan ( a / l ) c ( l a ) arctan ( a / l ) arctan ( a / ( l V f ) ) c ( l a ) + a 2 / 3
where c = 0.03 is an empirical parameter [19].
For any structure that is not a UD, nor a uniaxially aligned short fiber composite, as schematically shown in Figure 4, a subdivision of it into finite elements is necessary. With a further application of a lamination theory [7], the load share by each lamina layer in the laminated element is obtained. Then, the lamina layer is cut in an arbitrary way into a series of sheets, each of which contains at most one straight fiber yarn segment. Such a sheet is considered as an RVE of either Figure 3a or Figure 3b in its local coordinate system, with possibly a different Vf. Some sheets may even contain no fiber. Namely, they are pure matrix unit elements. The bridging model formulae summarized above are then applied to determine the homogenized internal stresses in the fiber and matrix, as well as the compliance tensor of the RVE. The homogenized internal stresses are then converted into true values as per the methods summarized in the subsequent section. An assemblage in terms of, e.g., an iso-strain or iso-stress scheme [24,25], from all of the RVEs gives rise to the overall internal true stresses in the fiber and matrix and the compliance tensor of the lamina layer.

3. True Stress

3.1. Background

Although the elastic properties of a composite are calculated by the bridging model formulae, Equations (21)–(25), correlate well with the experiments [21,22,23], and a predicted composite failure based on the internal stresses by Equations (6)–(11) can be very much different from reality. For instance, the transverse tensile strength of a UD composite thus obtained can be more than five times bigger than the measured counterpart [16,20]. A similar conclusion on the poor agreement is applicable to any other micromechanics theory.
This demonstrates that the homogenized internal stresses must be converted into true values before a failure detection can be made. The fiber stress field in Figure 3a is uniform, and the fiber true and homogenized stresses are the same; the fiber in Figure 3b generally does not fail, and its true stresses are immaterial. Further, a short fiber can be well regarded as an ellipsoid of a finite length, and its stresses are uniform [1]. We only need to evaluate the true stress of the matrix, which is equal to its homogenized counterpart multiplied by a factor. This factor was still called an SCF by the author [4,5], as it has all the physical characteristics of a classical SCF. Namely, it is dimensionless, defined by a stress ratio, relevant to a material failure, and dependent on a loading manner.

3.2. Definition

As an SCF of the matrix in a composite can no longer be defined by a classical method, which is given as a point-wise stress divided by an overall applied one, an averaging of the stresses in both the numerator and denominator must be made. The “point-wise stress” in the numerator has a geometrical characteristic of zero dimensions (0D), whereas the “overall applied stress” in the denominator is in fact a surface-averaged quantity (2D) with respect to the surface where the load is being applied. By similarity, the new definition must be given by 1D over 3D geometries, i.e., a line-averaged stress divided by a volume-averaged one. Further, the line averaging in the numerator should be along the outward normal of the failure surface of the composite under the given load, from the fiber cylinder to the matrix one in the RVE [4,5].

3.3. Matrix SCFs with Perfect Interface Bonding [4,5,6,7,8]

The matrix SCFs under transverse tension, transverse compression, transverse shear, and longitudinal shear, K 22 t , K 22 c , K23, and K12, are derived as
K 22 t = 1 + V f 2 A + V f 2 ( 3 V f V f ) B ( V f + 0.3 V m ) E 22 f + 0.7 V m E m 0.3 E 22 f + 0.7 E m
K 22 c = 1 V f 2 A σ u , c m σ u , t m 2 σ u , c m + B 2 ( 1 V f ) V f 2 1 2 σ u , c m σ u , t m 2 σ u , c m 2 +
( σ u , c m + σ u , t m ) V f σ u , c m 1 + σ u , c m σ u , t m σ u , c m V f σ u , c m σ u , t m σ u , c m + 1 2 σ u , c m σ u , t m 2 σ u , c m 2 ×
( V f + 0.3 V m ) E 22 f + 0.7 V m E m 0.3 E 22 f + 0.7 E m
K 23 = 2 σ u , s m K 22 t K 22 c σ u , t m σ u , c m
K 12 = 1 V f G 12 f G m G 12 f + G m { W ( V f ) 1 3 } ( V f + 0.3 V m ) G 12 f + 0.7 V m G m 0.3 G 12 f + 0.7 G m
A = 2 E 22 f E m ( ν 12 f ) 2 + E 11 f { E m ( ν 23 f 1 ) E 22 f [ 2 ( ν m ) 2 + ν m 1 ] } E 11 f [ E 22 f + E m ( 1 ν 23 f ) + E 22 f ν m ] 2 E 22 f E m ( ν 12 f ) 2
B = E m ( 1 + ν 23 f ) E 22 f ( 1 + ν m ) E 22 f [ ν m + 4 ( ν m ) 2 3 ] E m ( 1 + ν 23 f )
W ( V f ) = π V f [ 1 4 V f 1 32 1 256 V f 5 4096 V f 2 ]
σ u , c m and σ u , s m are the matrix compressive and shear strengths, respectively.

3.4. Matrix SCFs after Interface Debonding [8,9]

The in-plane shear and transverse tensile SCFs after the interface debonding are
K ^ 12 = G m ( 1 + V f ) G m + G 12 f + G 12 f G m + G 12 f π V f 2 + V f 1 V f Im 1 V f 1 2 i 1 V f 1
+ ln 1 V f 2 i 1 V f V f i 1 1 V f ln ( V f 1 V f )
+ ln 1 V f 1 2 i 1 V f 1 1 V f 1 + i 1 + V f 1 2 V f 2 i V f ( V f ) 2
ln 1 1 V f + i V f V f G 12 f + V m ( 0.3 G 12 f + 0.7 G m ) 0.3 G 12 f + 0.7 G m
K ^ 22 t = Re e 2 i ψ M ( b e i ψ ) ( a 2 / b b ) e i ψ N 2 N 1 ( a 2 b e i ψ )
+ e i ψ ( 2 + e 2 i ψ ) [ N ( b e i ψ ) N 3 ] ( V f + 0.3 V m ) E 22 f + 0.7 V m E m 2 ( b a ) ( 0.3 E 22 f + 0.7 E m )
N ( z ) = F z + a 2 k z ( z a e i ψ ) 0.5 + i λ ( z a e i ψ ) 0.5 i λ [ ( F 0.5 ) D a 2 z ]
N 1 ( z ) = F z + a 2 k z + 1 ξ ( z a e i ψ ) 0.5 + i λ ( z a e i ψ ) 0.5 i λ [ ( F 0.5 ) D a 2 z ]
N 2 = a F e i ψ + a k e i ψ ,   N 3 = F a e i ψ + e i ψ a k
M ( z ) = F a 2 k z 2 [ ( F 0.5 ) z + H + C z + D z 2 ] χ ( z )
F = 1 ( cos ψ + 2 λ sin ψ ) exp [ 2 λ ( π ψ ) ] + ( 1 k ) ( 1 + 4 λ 2 ) sin 2 ψ 4 k 2 2 ( cos ψ + 2 λ sin ψ ) exp [ 2 λ ( π ψ ) ]
H = a ( cos ψ + 2 λ sin ψ ) ( 0.5 F )
C = ( k 1 ) ( cos ψ 2 λ sin ψ ) a 2 exp [ 2 λ ( ψ π ) ]
D = ( 1 k ) a 3 exp [ 2 λ ( ψ π ) ]
χ ( z ) = ( z a e i ψ ) 0.5 + i λ ( z a e i ψ ) 0.5 i λ
k = G m 1 + κ 2 1 + ξ ( G m + κ 1 G 23 f ) ,   λ = - ( ln ξ ) / ( 2 π ) ,   ξ = ( G 23 f + κ 2 G m ) / ( G m + κ 1 G 23 f )
κ 1 = 3 4 v m , κ 2 = 3 ν 23 f 4 ν 12 f ν 21 f 1 + ν 23 f ,
In the above, z = x2 + ix3 and i = 1 ; ψ is the half central angle corresponding to a steady state interface crack, which is determined through the following equations [8]:
Re G 0 1 k 2 ( 1 k ) k exp ( i ϕ ) exp [ 2 λ ( ψ π ) ] R ( e i ϕ ) ϕ = ψ γ = 0
R ( exp ( i ϕ ) ) = [ exp ( i ( ϕ ) ) e i ψ ] 0.5 + i λ [ exp ( i ( ϕ ) ) e i ψ ] 0.5 i λ exp ( i ( ϕ ) )
G 0 = 1 ( cos ψ + 2 λ sin ψ ) exp [ 2 λ ( π ψ ) ] + ( 1 k ) ( 1 + 4 λ 2 ) sin 2 ψ 2 k k ( cos ψ + 2 λ sin ψ ) exp [ 2 λ ( π ψ ) ]
γ = 2 λ ( J 1 2 + J 2 2 ) J 1 2 + J 2 2 2 J 2 J 3 , i f ξ < 1 2 λ ( J 1 2 + J 2 2 ) J 1 2 + J 2 2 2 J 2 J 3 , i f ξ > 1
J 1 = k G 0 1 2 ( 1 k ) ξ exp ( 2 λ ψ ) cos ( ψ )
J 2 = 2 ( 1 k ) ξ exp ( 2 λ ψ ) sin ( ψ )
J 3 = 2 ( 1 k ) ξ exp ( 2 λ ψ ) [ J 1 cos ( ψ ) J 2 sin ( ψ ) ] / J 2
When ξ = 1, the crack angle ψ is indeterminate and the crack is named a singular crack. However, adjusting a fiber or matrix property parameter slightly, one can always attain ξ ≠ 1, since a deviation exists in measuring the parameter.

3.5. Matrix Longitudinal SCFs [10,11]

Only in a short fiber or particle composite (Figure 3b), can a longitudinal SCF of the matrix exist. No such SCF occurs in a continuous fiber composite, since the resulting matrix point-wise stress field is uniform [3]. Different from the derivations for all of the other matrix SCFs which are essentially based on the RVE geometry of Figure 3a, the RVE of Figure 3b used to derive the longitudinal SCFs of the matrix in a short fiber composite must be separated into three segments schematically shown in Figure 5, so that a longitudinal SCF of the matrix in a short fiber composite can become that in a continuous one when the fiber aspect ratio tends to infinity or when the fiber length ratio equals to 1. According to the condition that the matrix true stress resultant in the longitudinal direction in Figure 5a should equal to that in Figure 5b plus those in Figure 5c, one obtains (noting that the cross-sectional area in each figure is the same)
2 L K 22 t ( σ 11 m ) B M L = 2 ( l - h ) ( σ 11 m ) B M c + 2 ( L - l + h ) K 22 t , I ( σ 11 m ) B M h ,
where
( σ 11 m ) B M L = a 11 L σ 11 V f + ( 1 V f ) a 11 L ,
( σ 11 m ) B M c = a 11 c σ 11 V f I + ( 1 V f I ) a 11 c ,
( σ 11 m ) B M h = a 11 h σ 11 V f I I + ( 1 V f I I ) a 11 h ,
V f = ( a / b ) 2 ( l / L ) ,   V f I = ( a / b ) 2 = V f ( L / l ) ,   V f I I = ( a / b ) 2 ( h / ( L l + h ) )
a 11 c = E m / E 11 f ,
a 11 L = V f E m ( 1 V f ) E 11 f σ 11 L E 11 f [ ε L 1 l + ε L 2 ( L l ) ] E m [ ε L 1 l + ε L 2 ( L l ) ] σ 11 L ,
a 11 h = V f I I E m ( 1 V f I I ) E 11 f σ 11 L E 11 f [ ε L 1 h + ε L 2 ( L h ) ] E m [ ε L 1 h + ε L 2 ( L h ) ] σ 11 L ,
L = L l + h .
It is noted that Equation (61) is obtained from the second part of Equation (12), i.e., Equation (60), with L and l in Figure 3b and in Equations (16) replaced by L’ and h, respectively, in which L’ is given by Equation (62). The longitudinal tensile and compressive SCFs of the matrix in the short fiber composite are expressed below.
K 11 t = K 11 t , I , i f l h ( L l + h ) ( σ 11 m ) B M h L ( σ 11 m ) B M L K 11 t , I + l h L ( σ 11 m ) B M c ( σ 11 m ) B M L , i f l > h
K 11 c = K 11 c , I , i f l h ( L l + h ) ( σ 11 m ) B M h L ( σ 11 m ) B M L K 11 c , I + l h L ( σ 11 m ) B M c ( σ 11 m ) B M L , i f l > h
K 11 t , I = ( 1 V f ) E m σ 11 ε L 1 E 11 f V f σ 11 E m D I I sinh [ n I I ( L l ) ] [ J 0 ( n I I a ) J 1 ( n I I b ) Y 1 ( n I I b ) Y 0 ( n I I a ) ] ( 1 + ν m ) ( L l )
K 11 c , I = 1 V f σ 11 ε L 1 E 11 f V f σ 11 E m 2 ( 1 + ν m ) D I I a k = 1 N λ k f b + a 2 + b a 2 t k
f r = cosh n I I ( r a ) cot ϕ J 0 n I I r J 1 n I I b Y 1 n I I b Y 0 n I I r
φ = π 4 + 1 2 arcsin σ u , c m σ u , t m 2 σ u , c m
D I I a = σ F D σ z z , I I f , 0 n I I a cosh [ n I I ( l L ) ] ω f , I I
σ F D = σ z z , I f , 0 ω I tanh ( n I l ) + σ z z , I I f , 0 ω I I tanh [ n I I ( L l ) ] ω I tanh ( n I l ) + ω I I tanh [ n I I ( L l ) ]
σ z z , I f , 0 = E 11 f ε L 1 + 2 ν 12 f E m ( 1 + ν m ) ( 1 2 ν m ) ε L 1 ν m + C 1 ( 1 2 ν m ) C 2 a 2
σ z z , I I f , 0 = E m ε L 2 + 2 ν m E m ( 1 + ν m ) ( 1 2 ν m ) ε L 2 ν m + ( 1 + ν m ) ( 1 2 ν m ) σ 11 ( 1 ν m ) ε L 2 E m 2 ν m E m
ω I = ω r z , I ω f , I ,   ω I I = ω r z , I I ω f , I I ,   ω f , I = 2 n I a ω r z , I ,   ω f , I I = 2 n I I a ω r z , I I
ω r z , I = E m 1 + ν m n I a J 0 ( n I a ) J 1 ( n I b ) Y 1 ( n I b ) Y 0 ( n I a ) +   B 53 J 1 ( n I a ) + B 63 Y 1 ( n I a )
+ 2 ( 1 ν m ) J 1 ( n I a ) J 1 ( n I b ) Y 1 ( n I b ) Y 1 ( n I a )
ω r z , I I = E m 1 + ν m J 1 ( n I I a ) J 1 ( n I I b ) Y 1 ( n I I b ) Y 1 ( n I I a )
D 11 D 12 D 21 D 22 B 5 3 B 6 3 = M 1 M 2
D 11 = 2 ν 12 f E m 1 + ν m E 11 f n I a J 0 ( n I a ) + 2 ( 1 ν 12 f ) E m 1 + ν m J 1 ( n I a )
D 12 = 2 ν 12 f E m 1 + ν m E 11 f n I a Y 0 ( n I a ) + 2 ( 1 ν 12 f ) E m 1 + ν m Y 1 ( n I a )
D 21 = E 11 f n I a J 0 ( n I a ) ( 1 + ν m ) [ E 11 f E 22 f ( ν 12 f ) 2 ] [ ν 12 f E m ( 1 + ν 23 f ) ( 1 + ν m ) E 11 f ]
+ E 11 f J 1 ( n I a ) ν 12 f ( 1 + ν m ) [ E 11 f E 22 f ( ν 12 f ) 2 ] [ ( 1 + ν m ) E 22 f E m ( 1 + ν 23 f ) ] + 2 E m J 1 ( n I a ) 1 + ν m
D 22 = E 11 f n I a Y 0 ( n I a ) ( 1 + ν m ) [ E 11 f E 22 f ( ν 12 f ) 2 ] [ ν 12 f E m ( 1 + ν 23 f ) ( 1 + ν m ) E 11 f ]
+ E 11 f Y 1 ( n I a ) ν 12 f ( 1 + ν m ) [ E 11 f E 22 f ( ν 12 f ) 2 ] [ ( 1 + ν m ) E 22 f E m ( 1 + ν 23 f ) ] + 2 E m Y 1 ( n I a ) 1 + ν m
M 1 = 2 n I a E m [ ( 1 2 ν m ) ν 12 f 1 ] 1 + ν m + 2 ( 1 ν m ) E 11 f J 0 ( n I a ) J 1 ( n I b ) Y 1 ( n I b ) Y 0 ( n I a )
2 E m [ 2 ( 1 ν m ) + ( n I a ) 2 ν 12 f ] 1 + ν m + ( n I a ) 2 E 11 f J 1 ( n I a ) J 1 ( n I b ) Y 1 ( n I b ) Y 1 ( n I a )
M 2 = n I a 1 + ν m m 1 J 0 ( n I a ) J 1 ( n I b ) Y 1 ( n I b ) Y 0 ( n I a )
+ m 2 J 1 ( n I a ) J 1 ( n I b ) Y 1 ( n I b ) Y 1 ( n I a )
m 1 = E m ν 12 f ( 1 + ν 23 f ) ( 1 2 ν m ) 4 E 11 f [ 1 ( ν m ) 2 ] + ν 12 f E 22 f ( 1 + ν m ) 1 E 22 f ( ν 12 f ) 2 / E 11 f 2 E m
m 2 = n I a [ E m ν 12 f ( 1 + ν 23 f ) E 11 f ( 1 + ν m ) ] 1 E 22 f ( ν 12 f ) 2 / E 11 f 4 E m ( 1 ν m ) n I a
B 53 J 1 ( n I b ) Y 1 ( n I b ) + B 63 Y 1 2 ( n I b ) n I b [ Y 0 ( n I b ) J 1 ( n I b ) J 0 ( n I b ) Y 1 ( n I b ) ] = 0
n I I a J 0 ( n I I a ) 2 J 1 ( n I I a ) Y 1 ( n I I b ) n I I a Y 0 ( n I I a ) 2 Y 1 ( n I I a ) J 1 ( n I I b ) = 0
[ σ F D σ z z , I f , 0 ] cosh [ n I ( l h ) ] cosh ( n I l ) + σ z z , I f , 0 = 0.99 σ ¯ z z f ( 0 ) , if   l / a 100
[ σ F D σ z z , I f , 0 ] cosh [ n I a ( 100 h / a ) ] cosh ( 100 n I a ) + σ z z , I f , 0 = 0.99 σ ¯ z z f ( 0 ) , if   l / a > 100
σ ¯ z z f ( 0 ) = [ σ F D σ z z , I f , 0 ] / cosh ( n I l ) + σ z z , I f , 0 , i f l / a 100 σ z z , I f , 0 , i f l / a > 100
In the above, J0, J1, Y0, and Y1 are, respectively, the zero- and one-order Bessel functions of the first and second classes, and a and b are the radii of the fiber and matrix cylinders in the RVE (Figure 3b), respectively. C1, C2, ε L 1 , and ε L 2 are the coefficients solved from Equations (16) to (20), and g is defined in Equation (17). In Equation (66), N is the number of a Gaussian integration, and tk and λk are the integration points and weighing coefficients, respectively. In general, taking N = 5 can attain an enough high accuracy.

3.6. Numerical Examples

The matrix SCF formulae under a transverse tension after interface debonding and under a longitudinal load in a short fiber composite are extremely complicated, as the elasticity solutions to the matrix stress fields themselves are very complicated [18,26,27]. After line integration, the complication level is elevated one order higher. Fortunately, this author implemented all of the SCF formulae, Equations (27)–(53) and Equations (63)–(89), into an Excel table-based program to calculate them. For any question in using the program, the reader can write to the author. The SCFs of the matrices in the nine independent material systems used in three world-wide failure exercises (WWFEs) [28,29,30] are summarized in Table 1, whereas the input data for the calculations taken from [28,29,30] are listed in Table 2.
Simply speaking, a relative difference between the predicted strengths of a composite made from the nine material systems based on the homogenized and the true stresses will be at least 34% and at most 669%, as indicated in Table 1.

3.7. True Stresses

Let { d σ i m } = { d σ 11 m , d σ 22 m , d σ 33 m , d σ 23 m , d σ 13 m , d σ 12 m } T be the homogenized stress increments of the matrix due to an arbitrary load increments {j} applied on the RVE of either Figure 3a or Figure 3b. It is evaluated from Equation (4) by replacing {σj} with {j} on its right hand side. The overall true stresses of the matrix are updated through
{ σ ¯ i m } n = { σ ¯ 11 m , σ ¯ 22 m , σ ¯ 33 m , σ ¯ 23 m , σ ¯ 13 m , σ ¯ 12 m } n = { σ ¯ i m } n 1 + { d σ ¯ i m } , n = 1 , 2 ,
{ d σ ¯ i m } = { K 11 d σ 11 m , K 22 d σ 22 m , K 33 d σ 33 m , K 23 d σ 23 m , K ¯ 12 d σ 13 m , K ¯ 12 d σ 12 m }  
where
K 11 = 1 , f o r   c o n t i n u o u s   f i b e r   c o m p o s i t e K 11 t , f o r   o t h e r   c o m p o s i t e   w i t h   d σ 11 m 0 K 11 c , f o r   o t h e r   c o m p o s i t e   w i t h   d σ 11 m < 0
K 22 = K 22 t , i f   d σ 22 m > 0   w i t h   p e r f e c t   int e r f a c e K ^ 22 t , i f   d σ 22 m > 0   w i t h   d e b o n d e d   int e r f a c e K 22 c , i f d σ 22 m < 0
K 33 = K 22 t , i f   d σ 33 m > 0   w i t h   p e r f e c t   int e r f a c e K ^ 22 t , i f   d σ 33 m > 0   w i t h   d e b o n d e d   int e r f a c e K 22 c , i f   d σ 33 m < 0
K ¯ 12 = K 12 , with   p e r f e c t   int e r f a c e K ^ 12 , with   d e b o n d e d   int e r f a c e
The fiber true stresses are updated from:
{ σ ¯ i f } n = { σ ¯ i f } n 1 + { d σ j f } , n = 1 , 2 ,
It is noted that the internal stresses occurred in Figure 4 must be the internal true stresses.

3.8. Roles of the True Stresses

Eshelby [1] found that if a solid embedded in an unbounded matrix is an ellipsoid, the stress field in it is uniform no matter what kind of load is applied to the matrix. A continuous fiber in a composite can be regarded as an infinitely long ellipsoid, whereas a short fiber or particle can be considered as an ellipsoid of a finite length. Thus, the Eshelby’s finding implies that one can estimate any failure of the fiber in a composite by directly comparing its stresses against the strengths of the monolithic fiber determined independently, since the stress fields in obtaining the fiber strengths are uniform. This also explains why a fiber controlled composite failure can be estimated reasonably by anybody in general.
However, the stress field in the matrix outside the fiber is not uniform. Any failure of the matrix in a composite cannot be assessed by directly comparing its stresses against the strengths of the monolithic matrix measured independently. The reason is in that the matrix strengths are measured in uniform stress fields, whereas the matrix stresses in the composite are not uniform. The two kinds of quantities are not comparable. To estimate any matrix-induced composite failure, one must firstly convert the non-uniform stresses of the matrix into equivalent, uniform ones. Otherwise, one cannot obtain the critical stress or deformation data of the matrix to detect any of the matrix-dominated composite properties outside a linear elastic range, such as a plastic behavior, failure, strength, etc. Simply speaking, the true stresses of the matrix determined through Equation (90) stand for the uniform values.

4. Selected Applications

4.1. Uniaxial Strengths of UD Composites

First, substitute the material property data of Table 2 into Equations (6)–(11) to calculate the homogenized internal stresses of the composite subjected to uniaxial load one by one, and let the major internal stress attain either the fiber or matrix strength of Table 2. The resulting external load is taken as the composite strength. Except for the longitudinal strength, which is assumed to be caused by a fiber failure, all of the other uniaxial strengths are resulted from the matrix failures. They are named homogenized stresses-based predictions. The averaged errors between the predictions and the experiments are summarized in Table 3. Second, after converting the homogenized stresses into the true values and assuming that the interface bonding for each of the composites is perfect up to failure, similar predictions can be made. Again, the averaged errors for the true stresses-based predictions are listed in the table as well.
From Table 3, the following conclusions can be made: (1) Only a longitudinal strength is independent of a true stress effect. (2) The overall averaged error of the predictions based on the homogenized stresses is 3.92 times bigger than that upon the true stresses. If the longitudinal strengths are excluded, the difference becomes 5.12 times. (3) Using the true stresses, the averaged prediction errors bigger than 20% are for the composite strengths under longitudinal compression, transverse tension, and transverse compression. The error for the longitudinal compressive strength comes from an initial fiber misalignment; the error for the transverse compressive strength is due to the fact that the outward normal of the resulting failure surface is always perpendicular to the fiber axis; and the error for the transverse tensile strength is attributed to the interface debonding. How to reduce these errors will be highlighted subsequently. (4) The predictions for the in-plane and transverse shear strengths are high enough, indicating that an interface debonding has a negligible effect on the shear load sustaining ability of the composite.
In fact, if an interface debonding is taken into account, the in-plane shear strength formula is changed to
σ 12 u = ( V f + V m a 66 ) [ σ ^ e m / 3 + ( σ u , s m K 12 σ ^ e m / 3 ) / K ^ 12 ] / a 66
where σ ^ e m is determined by a subsequent Equation (101). Having done so, the overall error for the in-plane shear strengths is changed from 13.1% (Table 3) to 11.1%.
Any experiment has a deviation. The measurement errors exist for not only the constituent properties and the fiber volume content, but also the composite property data. All of them gathered together generate a cumulative effect. Thus, an expected overall error between measured and predicted elastic moduli of the composites using independent constituent properties and fiber geometric parameters can only be 10%. Indeed, this error is 10.4% [16] if the data of Table 2 are substituted into Equations (21)–(25). An expectation for an even higher accuracy is unrealistic in general. Similarly, an expected overall error for strength prediction of a composite can be as small as 20% if only the independent constituent properties and the fiber content are available, since the prediction needs more input data and a larger cumulative error is involved. A higher prediction accuracy depends on the input of some composite ultimate data. Table 3 indicates that a failure and strength prediction for a composite in terms of the bridging model and the matrix true stress theory is near to the expectation.

4.2. Interface Debonding

It is known in the composite community that most composite failures are initiated from fiber and matrix interface debonding. However, given an arbitrary load, when does the interface debonding occur? Both the existing studies [31,32,33] and Table 3 indicate that the interface debonding only has an effect on the transverse tensile load sustaining ability of the composite.
Let the composite under a transverse tension σ22 attain an interface debonding at σ ^ 22 . Due to a Poisson’s retardant force in the direction perpendicular to the loading, as shown in Figure 6, the interface debonding soon reaches stable, and a load variation from the initial to the stable debonding stages is negligible. One can assume that the interface before the critical load σ ^ 22 is in perfect bonding, and after it, in stable debonding. The transverse failure condition reads
K 22 t σ ^ 22 m + K ^ 22 t ( σ 22 m , Y σ ^ 22 m ) = σ u , t m
In Equation (98), σ ^ 22 m and σ 22 m , Y are the matrix homogenized stresses by replacing σ22 in Equation (9) with σ ^ 22 and Y, respectively, where Y is the transverse tensile strength of the composite. From Equation (98), one obtains the critical transverse load as
σ ^ 22 = K ^ 22 t Y K ^ 22 t K 22 t ( V f + 0.3 V m ) E 22 f + 0.7 V m E m ( 0.3 E 22 f + 0.7 E m ) ( K ^ 22 t K 22 t ) σ u , t m
A necessary and sufficient condition for an interface debonding in any composite subjected to an arbitrary load is given by
σ ¯ e m σ ^ e m and σ ¯ m 1 > 0
σ ^ e m = ( σ ^ 11 m ) 2 + ( K 22 t σ ^ 22 m ) 2 K 22 t σ ^ 11 m σ ^ 22 m
In the above equation, σ ^ 11 m is calculated from Equation (7) with σ22 replaced by σ ^ 22 and σ11 = σ33 = 0; σ ¯ e m and σ ¯ m 1 are the Mises and the first principal true stresses calculated from Equation (90), respectively. Therefore, in addition to the fiber and matrix properties, one only needs to know the transverse tensile strength Y of the composite.

4.3. Interface Slippage

An advanced, e.g., carbon fiber-reinforced composite is generally considered to be brittle, since its load-deflection curves in both the longitudinal and transverse directions are essentially linear up to failure. However, significant non-linear deformation can occur in a composite, especially when it is subjected to an in-plane shear or an off-axial load [34]. The in-plane shear nonlinear deformation of the composite can be even bigger than the elastic–plastic deformation of the pure matrix. Shear deformation of the composite is one of the most difficult deformations to be analyzed [35]. We found [15] that a relative slippage displacement between debonded fiber and matrix interfaces is the key reason for the composite under an in-plane shear to generate a large deformation, as indicated in Figure 7. Before the interface debonding, a longitudinal plane denoted by the line FD in a transverse cross-section (Figure 7b) deforms as schematically shown in Figure 7a. After the debonding, the point D further moves to D’ (Figure 7a).
Under a transverse tension, an interface debonding will soon reach stable due to a Poisson’s internal retardant (Figure 6), resulting in hardly any effect on a transverse tensile stress–strain curve. In the contrast, no such retardant exists in the case of an in-plane shear loading. An interface debonding will propagate until there is a complete separation between the fiber and the matrix. Thus, as long as a fiber pulling out is seen, there must be a shear load component or a bi-axial transverse tension. Explicit formulae for an interface debonding-induced shear slippage displacement were derived in [15,20]. Its application to the prediction of shear deformation for the E-glass/MY750 UD composite [28] is shown in Figure 8. More details can refer to [20].
Therefore, a nonlinear constitutive relation of the composite is attributed to three contributions: (1) the matrix plasticity or nonlinear deformation, which is well known in the composite community, (2) the matrix true stress effect, which occurs in an instantaneous compliance tensor of the matrix (by any classical theory, a matrix instantaneous compliance tensor outside linear elastic range contains the current matrix stresses, which must be true values in the composite), and (3) a relative slippage displacement between debonded fiber and matrix interfaces.

4.4. Fiber Kinking or Splitting Failure

A particular phenomenon was already recognized in engineering. Namely, the fiber strength from, e.g., T300 to T1100 is elevated significantly, but the longitudinal compressive strength of the resulting composite almost keeps unchanged. Perhaps the most well known explanation for this was attributed to a matrix shear nonlinearity [36]. However, even for a very brittle matrix, the phenomenon still exists. Using the matrix true stress theory, the phenomenon can be displayed accurately.
Let θ c f , 0 be an initial fiber misalignment angle in the composite, which is almost inevitable to occur during a composite fabrication. A longitudinal compression in the misaligned coordinate system (Figure 3c) has a shear stress component, which will cause a misalignment angle increment θ c f , 1 . Referring to Figure 9, external loads {σ11,σ22, σ12} generate three stress components in the local coordinate system, ( x 1 I , x 2 I ), as
σ 11 I = σ 11 + σ 22 2 + σ 11 σ 22 2 cos ( 2 θ c f ) + σ 12 sin ( 2 θ c f )
σ 22 I = σ 11 + σ 22 2 σ 11 σ 22 2 cos ( 2 θ c f ) σ 12 sin ( 2 θ c f )
σ 12 I = σ 11 σ 22 2 sin ( 2 θ c f ) + σ 12 cos ( 2 θ c f )
where θ c f = θ c f , 0 + θ c f , 1 . Equations (102)–(104) are substituted into Equations (6)–(11) to calculate the homogenized stresses, and further into Equations (90) and (96) to obtain the true stresses. Letting the matrix attain a failure status, the misalignment angle increment θ c f , 1 is determined as [20]
θ c f , 1 = σ 12 e q ( 1 V f ) π 2 G m [ V f 1 V f + sin 1 ( V f ) ]
σ 12 e q = V f + V m a 66 K 12 a 66 ( σ ¯ 11 m + σ ¯ 22 m ) c 1 1 2 c 1 e 1 2 c 1 2
c 1 = 2 σ u , c m 4 ( σ u , s m ) 2 ( σ u , c m ) 2 ,
e 1 = [ 4 ( σ u , s m ) 2 + ( σ u , c m ) 2 ] 2 8 ( σ u , c m ) [ 4 ( σ u , s m ) 2 ( σ u , c m ) 2 ] .
It is noted that the matrix true stresses in Equation (106), σ ¯ 11 m and σ ¯ 22 m , should be expressed in the local system ( x 1 I , x 2 I ). If the corresponding load at the matrix failure is smaller than the fiber fracture load, the composite is said to assume a fiber kinking failure. Otherwise, the fiber fracture is attained. In the former case, a further increase in the fiber strength has essentially no contribution to the composite longitudinal strength. In the latter one, the longitudinal compressive strength of the composite reaches its maximum.
The initial fiber misalignment angle, θ c f , 0 , is not easy to measure. Yurgartis reported that more than 80% of the initial misalignment angles were around 10 [37]. Through a numerical test, we found that θ c f , 0 =1.50 resulted in the best prediction for the longitudinal compressive strengths of the nine composites in Table 1, with an overall correlation error from 25.1% in Table 3 (where θ c f , 0 =00) reduced to 16.1% [20].
A fiber splitting failure can be analyzed similarly [13]. The only difference is in that the misalignment increment θ c f , 1 is positive under a longitudinal compression, but negative under a longitudinal tension [13].

4.5. Matrix Compression-Induced Composite Failure

Except for a longitudinal compression on a composite, a failure of which can be caused by either a fiber fracture or a matrix shear failure, a matrix compression-dominated composite failure has two distinct characteristics (Figure 10): (1) the failure surface has an inclined angle with the loading direction [38,39]; and (2) the outward normal to the failure surface is generally perpendicular to the fiber axis.
In accordance with the above failure characteristics, a physics-based failure criterion was developed by the author et al. to detect a matrix compression-dominated composite failure [40,41]:
f ( σ n m , τ n m ) = τ n m + μ σ n m = ( τ n s m ) 2 + ( τ n t m ) 2 + μ σ n m σ u , s m / K n
K n = ( K 12 τ n s m ) 2 + ( K 23 τ n t m ) 2 ( τ n s m ) 2 + ( τ n t m ) 2
μ = 0.5 , i f   U D   c o m p o s i t e   i s   u n d e r   2 D   c o m p r e s s i o n 0.2 , i f   U D   c o m p o s i t e   i s   u n d e r   3 D   c o m p r e s s i o n
σ n m ( θ ) = σ 22 m cos 2 θ + σ 33 m sin 2 θ + 2 σ 23 m sin θ cos θ τ n t m ( θ ) = ( σ 33 m σ 22 m ) sin θ cos θ + σ 23 m ( cos 2 θ sin 2 θ ) τ n s m ( θ ) = σ 13 m sin θ + σ 12 m cos θ
in which θ is the inclined angle between the failure surface (parallel to the fiber axis) and the x3 coordinate (in thickness direction of the UD composite) shown in Figure 11.
Using Equations (109)–(112), the predicted transverse compressive strengths of the nine composites in Table 1 have a much better correlation with the experiments, and the relative error of 23.2% in Table 3 is reduced to 13.8% [20,40].

4.6. Applications to Other Areas

The matrix true stress theory was also applied by many other researchers in their analysis for the failure and strength behaviors of composites under different load conditions. An incomplete list of the published papers by the other people using the matrix true stress theory is given in [42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63].

5. Conclusions

Eshelby’s work explained why a fiber-dominated composite failure can be estimated reasonably, since the uniform fiber stresses in any composite under arbitrary loads are comparable with the fiber strengths determined in uniform stress fields. The matrix stresses in the composite, however, are not uniform, and hence any matrix-dominated composite failure, as well as nonlinear behavior, cannot be assessed by comparing the non-uniform matrix stresses with its monolithic strengths or other critical parameters measured in uniform stress states. The two kinds of quantities of the matrix are not comparable. A conversion of the non-uniform matrix stresses into equivalent uniform ones must be carried out. The matrix true stresses stand for its equivalent uniform values, and the conversion is achieved through the author’s matrix true stress theory summarized in this paper. As most composite failures resulted from matrix failures, none of them can be properly estimated without the matrix true stresses. On the other hand, almost all of the mechanics of composite problems can be efficiently resolved based on the matrix true stress theory. Only some typical applications of the theory are highlighted in the paper. Analytical formulae for almost all of the true stress components were derived by the author, and are summarized in the paper. They were implemented into an Excel table-based program. If the reader has any question relevant to the Excel table-based program for calculating all of the matrix SCFs in a composite, he/she is welcome to write to the author at [email protected].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ma16020774/s1.

Funding

Financial support from the National Natural Science Foundation of China with a grant No. 11832014 is acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Available upon request.

Acknowledgments

Great helps and contributions from my postgraduate students including, among others, Lin Liu and Yi Zhou, Li-Min Xin, Chun-Chun Zhang and Wei-Jing Guo and Feng Jiang, Li-Sheng Wang and Hong-Bo Huang are highly appreciated.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Eshelby, J.D. The determination of the elastic field of an ellipsoidal inclusion and related Problems. Proc. Royal Soc. 1957, A240, 367–396. [Google Scholar]
  2. Cheng, S.; Chen, D. On the Stress Distribution in Laminae. J. Reinf. Plast. Compos. 1988, 7, 136–144. [Google Scholar] [CrossRef]
  3. Chen, T.; Dvorak, G.J.; Benveniste, Y. Stress fields in composites reinforced by coated cylindrically orthotropic fibers. Mech. Mater. 1990, 9, 17–32. [Google Scholar] [CrossRef]
  4. Huang, Z.-M.; Liu, L. Predicting strength of fibrous laminates under triaxial loads only upon independently measured constituent properties. Int. J. Mech. Sci. 2014, 79, 105–129. [Google Scholar] [CrossRef]
  5. Huang, Z.-M.; Xin, L.M. In situ strengths of matrix in a composite. Acta Mech. Sin. 2017, 33, 120–131. [Google Scholar] [CrossRef]
  6. Huang, Z.-M.; Xin, L.M. Stress Concentration Factor in Matrix of a Composite Subjected to Transverse Compression. Int. J. Appl. Mech. 2016, 8, 1650034. [Google Scholar] [CrossRef]
  7. Huang, Z.-M. Failure and Strength of Composite Materials; Science Press: Beijing, China, 2018. (In Chinese) [Google Scholar]
  8. Zhou, Y.; Huang, Z.-M.; Liu, L. Prediction of Interfacial Debonding in Fiber-Reinforced Composite Laminates. Polym. Compos. 2019, 40, 1828–1841. [Google Scholar] [CrossRef]
  9. Zhou, Y.; Huang, Z.-M. Prediction of in-plane shear properties of a composite with debonded interface. Appl. Compos. Mater. 2022, 29, 901–935. [Google Scholar] [CrossRef]
  10. Huang, Z.-M.; Guo, W.J.; Huang, H.B.; Zhang, C.C. Tensile Strength Prediction of Short Fiber Reinforced Composites. Materials 2021, 14, 2708. [Google Scholar] [CrossRef]
  11. Guo, W.J.; Huang, Z.-M. Compressive strength of a short fiber reinforced composite. Chinese J. Solid Mech. 2022, 43, 1–15. (In Chinese) [Google Scholar]
  12. Sun, Q.; Guo, H.; Zhou, G.; Meng, Z.; Chen, Z.; Kang, H.; Keten, S.; Su, X. Experimental and computational analysis of failure mechanisms in unidirectional carbon fiber reinforced polymer laminates under longitudinal compression loading. Comp. Struct. 2018, 203, 335–348. [Google Scholar] [CrossRef]
  13. Zhou, Y.; Huang, Z.-M.; Wang, J.X. Splitting failure of a composite under longitudinal tension. Poly. Comp. 2022. [Google Scholar] [CrossRef]
  14. Zhou, Y.; Huang, Z.-M. Failure of fiber-reinforced composite laminates under longitudinal compression. J. Comp. Mater. 2019, 53, 3395–3411. [Google Scholar] [CrossRef]
  15. Zhou, Y.; Huang, Z.-M. Shear deformation of a composite until failure with a debonded interface. Comp. Struct. 2020, 254, 112797. [Google Scholar] [CrossRef]
  16. Huang, Z.-M. On micromechanics approach to stiffness and strength of unidirectional composites. J. Reinf. Plastics Comp. 2019, 38, 167–196. [Google Scholar] [CrossRef]
  17. Huang, Z.-M. Simulation of the Mechanical Properties of Fibrous Composites by the Bridging Micromechanics Model. Comp. Part A 2001, 32, 143–172. [Google Scholar] [CrossRef]
  18. Huang, Z.-M.; Zhang, C.C.; Xue, Y.D. Stiffness prediction of short fiber reinforced composites. Int. J. Mech. Sci. 2019, 161–162, 105068. [Google Scholar] [CrossRef]
  19. Huang, H.B.; Huang, Z.-M. Micromechanical prediction of elastic-plastic behavior of a short fiber or particle reinforced composite. Comp. Part A 2020, 134, 105889. [Google Scholar] [CrossRef]
  20. Huang, Z.-M. Constitutive Relation, Deformation, Failure and Strength of Composites Reinforced with Continuous/Short Fibers or Particles. Comp. Struct. 2021, 262, 113279. [Google Scholar] [CrossRef]
  21. Ryan, S.; Wicklein, M.; Mouritz, A.; Riedel, W.; Schäfer, F.; Thoma, K. Theoretical prediction of dynamic composite material properties for hypervelocity impact simulations. Int. J. Impact Eng. 2009, 36, 899–912. [Google Scholar] [CrossRef] [Green Version]
  22. Younes, R.; Hallal, A.; Fardoun, F.; Chehade, F.H. Comparative review study on elastic properties modeling for unidirectional composite materials. In Composites and Their Properties; Hu, N., Ed.; IntechOpen: London, UK, 2012; pp. 391–408. [Google Scholar]
  23. Vignoli, L.L.; Savi, M.A.; Pacheco, P.M.C.L.; Kalamkarov, A.L. Comparative analysis of micromechanical models for the elastic composite laminae. Comp. Part B 2019, 174, 106961. [Google Scholar] [CrossRef]
  24. Huang, Z.-M.; Ramakrishna, S. Micromechanical modelling approaches for the stiffness and strength of knitted fabric composites: A review & comparative study. Comp. Part A 2000, 31, 479–501. [Google Scholar]
  25. Huang, Z.-M.; Fujihara, K.; Ramakrishna, S. Flexural Failure Behavior of Laminated Composites Reinforced with Braided Fabrics. AIAA J. 2002, 40, 1415–1420. [Google Scholar] [CrossRef]
  26. Toya, M. A crack along the interface of a circular inclusion embedded in an infinite solid. J. Mech. Phys. Solids. 1974, 22, 325–348. [Google Scholar] [CrossRef]
  27. Jiang, Z.; Liu, X.; Li, G.; Lian, J. A new analytical model for three-dimensional elastic stress field distribution in short fibre composite. Mater. Sci. Engin. A. 2004, 366, 381–396. [Google Scholar]
  28. Soden, P.D.; Hinton, M.J.; Kaddour, A.S. Lamina properties, lay-up configurations and loading conditions for a range of fibre-reinforced composite laminates. Comp. Sci. Technol. 1998, 58, 1011–1022. [Google Scholar] [CrossRef]
  29. Kaddour, A.S.; Hinton, J.M. Input data for test cases used in benchmarking triaxial failure theories of composites. J. Comp. Mater. 2012, 46, 2295–2312. [Google Scholar] [CrossRef]
  30. Kaddour, A.S.; Hinton, M.J.; Smith, P.A.; Li, S. Mechanical properties and details of composite laminates for test cases used in the third world-wide failure exercise. J. Comp. Mater. 2013, 47, 2427–2442. [Google Scholar] [CrossRef]
  31. Hobbiebrunken, T.; Hojo, M.; Adachi, T.; Jong, C.D.; Fiedler, B. Evaluation of interfacial strength in CF/epoxies using FEM and in-situ experiments. Comp. Part A 2006, 37, 2248–2256. [Google Scholar] [CrossRef]
  32. Mortell, D.J.; Tanner, D.A.; McCarthy, C.T. In-situ SEM study of transverse cracking and delamination in laminated composite materials. Comp. Sci. Technol. 2014, 105, 118–126. [Google Scholar] [CrossRef]
  33. Koyanagi, J.; Ogihara, S.; Nakatani, H.; Okabe, T.; Yoneyama, S. Mechanical properties of fiber/ matrix interface in polymer matrix composites. Adv. Comp. Mater. 2014, 23, 551–570. [Google Scholar] [CrossRef]
  34. Fallahi, H.; Taheri-Behrooz, F.; Asadi, A. Nonlinear Mechanical Response of Polymer Matrix Composites: A Review. Polym. Rev. 2020, 60, 42–85. [Google Scholar] [CrossRef]
  35. Totry, E.; Molina-Aldareguía, J.M.; González, C.; LLorca, J. Effect of fiber, matrix and interface properties on the in-plane shear deformation of carbon-fiber reinforced composites. Comp. Sci. Technol. 2010, 70, 970–980. [Google Scholar] [CrossRef]
  36. Budiansky, B.; Fleck, N.A. Compressive failure of fibre composites. J. Mech. Phy. Solids. 1993, 41, 183–211. [Google Scholar] [CrossRef]
  37. Yurgartis, S. Measurement of small angle fiber misalignments in continuous fiber composites. Comp. Sci. Technol. 1987, 30, 279–293. [Google Scholar] [CrossRef]
  38. Puck, A.; Shurmann, H. Failure analysis of FRP laminates by means of physically based phenomenological models. Comp. Sci. Technol. 1998, 58, 1045–1067. [Google Scholar] [CrossRef]
  39. Pinho, S.T.; Davila, C.G.; Camnaho, P.P.; Lannucci, L.; Robinson, P. Failure Models and Criteria for FRP under In-Plane or Three-Dimensional Stress States Including Shear Non-Linearity; NASA/TM-2005-213530; NASA Langley Research Center: Hampton, VA, USA, 2005. [Google Scholar]
  40. Huang, Z.-M.; Wang, L.S.; Jiang, F.; Xue, Y.D. Detection on matrix induced composite failures. Comp. Sci. Technol. 2021, 205, 108670. [Google Scholar] [CrossRef]
  41. Wang, L.S.; Huang, Z.-M. On Strength Prediction of Laminated Composites. Comp. Sci. Technol. 2022, 219, 109206. [Google Scholar] [CrossRef]
  42. Hafiychuk, V. Modeling of Microstructure for Uncertainty Assessment of Carbon Fiber Reinforced Polymer Composites. In Proceedings of the 2016 IEEE Aerospace Conference, Big Sky, MT, USA, 5–12 March 2016; pp. 1–9. [Google Scholar]
  43. Xin, H.; Liu, Y.; Mosallam, A.; He, J.; Du, A. Evaluation on material behaviors of pultruded glass fiber reinforced polymer (GFRP) laminates. Comp. Struct. 2017, 182, 283–300. [Google Scholar] [CrossRef] [Green Version]
  44. Zhu, X.; Chen, X.; Zhai, Z.; Yang, Z.; Chen, Q. The effects of thermal residual stresses and interfacial properties on the transverse behaviors of fiber composites with different microstructures. Sci. Eng. Compos. Mater. 2017, 24, 41–51. [Google Scholar] [CrossRef]
  45. Toh, W.; Tan, L.B.; Tse, K.M.; Giam, A.; Raju, K.; Lee, H.P.; Tan, V.B.C. Material characterization of filament-wound composite pipes. Comp. Struct. 2018, 206, 474–483. [Google Scholar] [CrossRef]
  46. Vignoli, L.L.; Savi, M.A. Multiscale Failure Analysis of Cylindrical Composite Pressure Vessel: A Parametric Study. Latin Amer. J. Solids Struct. 2018, 15, e63. [Google Scholar] [CrossRef] [Green Version]
  47. Xin, H.; Mosallam, A.S.; Liu, Y.; Veljkovic, M.; He, J. Mechanical characterization of a unidirectional pultruded composite lamina using micromechanics and numerical homogenization. Constr. Build. Mater. 2019, 216, 101–118. [Google Scholar] [CrossRef] [Green Version]
  48. Ren, M.-F.; Zhang, X.-W.; Huang, C.; Wang, B.; Li, T. An integrated macro/micro-scale approach for in situ evaluation of matrix cracking in the polymer matrix of cryogenic composite tanks. Comp. Struct. 2019, 216, 201–212. [Google Scholar] [CrossRef]
  49. Du, Z.; Ni, X.; Liu, X.; Yu, J.; Wu, Y. Effect of nano-interface stress concentration on fracture stress of composite grains. J. Harbin Inst. Technol. 2019, 51, 118–124. (In Chinese) [Google Scholar]
  50. Vignoli, L.L.; Savi, M.A.; Pacheco, P.M.C.L.; Kalamkarov, A.L. Micromechanical analysis of transversal strength of composite laminae. Comp. Struct. 2020, 250, 112546. [Google Scholar] [CrossRef]
  51. Vignoli, L.L.; Savi, M.A.; Pacheco, P.M.C.L.; Kalamkarov, A.L. Multiscale approach to predict strength of notched composite plates. Comp. Struct. 2020, 253, 112827. [Google Scholar] [CrossRef]
  52. Vignoli, L.L.; Savi, M.A.; Pacheco, P.M.C.L.; Kalamkarov, A.L. Micromechanical analysis of longitudinal and shear strength of composite laminae. J. Comp. Mater. 2020, 54, 4853–4873. [Google Scholar] [CrossRef]
  53. Jiang, H.; Ren, Y.; Jin, Q. A novel synergistic multi-scale modeling framework to predict micro- and meso-scale damage behaviors of 2D triaxially braided composite. Int. J. Damage Mech. 2022, 31, 108–141. [Google Scholar] [CrossRef]
  54. Liu, W.L.; Chen, P.H. Theoretical analysis and experimental investigation of the occurrence of fiber bridging in unidirectional laminates under Mode I loading. Comp. Struct. 2021, 257, 113383. [Google Scholar] [CrossRef]
  55. Sun, Y.; Liu, Y.; Wang, C.; Zuo, Y.; Xin, H. Web buckling mechanism of pultruded GFRP bridge deck profiles subjected to concentrated load. Structures 2021, 34, 3789–3805. [Google Scholar] [CrossRef]
  56. Zhang, X.W.; Xing, H.Z.; Zhang, J.K.; Ren, M.F. Multi-scale analysis of cryogenic composite structures considering interfacial phase. J. Mech. Strength 2021, 6, 1417–1421. (In Chinese) [Google Scholar]
  57. Lv, Q.; Zhao, Z.; Li, C.; Zhang, C. Progressive damage and failure simulation of 2.5D woven composites. Acta Mater. Compos. Sin. 2021, 38, 2747–2757. (In Chinese) [Google Scholar]
  58. Wang, Y.C.; Zhou, Y.X.; Jin, C.H. An analytical nonlinear model for plain-woven composites under off-axis loads. Comp. Struct. 2022, 296, 115905. [Google Scholar] [CrossRef]
  59. Liu, L.B.; Zhang, X.H.; Wang, Z.B.; Wang, Y.N.; Guo, J.Z. Micromechanics Modeling of Transverse Tensile Strength for Unidirectional CFRP Composite. Materials 2022, 15, 8577. [Google Scholar] [CrossRef]
  60. Sha, Y.; Huang, J.; Luo, L.; Zhou, Y. Establishment and Modification of Structural Strength Calculation Model for Fiber Reinforced Titanium Matrix Composites. J. Propuls. Technol. 2023, 44, 201–210. (In Chinese) [Google Scholar]
  61. Wang, M. Damage and failure analysis of carbon fiber reinforced composites across macro-meso-micro scales. Ph.D. Thesis, Southeast University, Nanjing, China, 2020. (In Chinese). [Google Scholar]
  62. Zhang, X.W. Multi-scale analysis method of composite tanks at cryogenic temperatures considering fiber/matrix interface. Ph.D. Thesis, Dalian University of Technology, Dalian, China, 2020. (In Chinese). [Google Scholar]
  63. Qi, S.J. Investigation on the VAT reinforcing technique for open-hole CFRP laminates. Ph.D. Thesis, University of Chinese Academy of Sciences, Beijing, China, 2020. (In Chinese). [Google Scholar]
Figure 1. Schematic composite properties at different scales of stresses.
Figure 1. Schematic composite properties at different scales of stresses.
Materials 16 00774 g001
Figure 2. (a) Longitudinal tension-induced fiber splitting failure, (b) longitudinal compression-induced fiber kinking failure [12], (c) load decomposition in fiber misligned coordinate system.
Figure 2. (a) Longitudinal tension-induced fiber splitting failure, (b) longitudinal compression-induced fiber kinking failure [12], (c) load decomposition in fiber misligned coordinate system.
Materials 16 00774 g002
Figure 3. RVEs of (a) continuous fiber, and (b) short fiber reinforced composites.
Figure 3. RVEs of (a) continuous fiber, and (b) short fiber reinforced composites.
Materials 16 00774 g003
Figure 4. Analysing steps for any composite structure: 1. the structure is discretized into finite elements, and a laminate theory is applied to determine load share by each layer of the laminated element; 2. the layer is cut into a series of sheets each containing at most one straight fiber yarn segment, the bridging model is applied to determine the homogenized internal stresses and compliance tensor of each sheet, and the homogenized stresses are converted into true values as per the methods described subsequently; and 3. an assemblage from all of the sheets give the internal true stresses and compliance tensor of the layer.
Figure 4. Analysing steps for any composite structure: 1. the structure is discretized into finite elements, and a laminate theory is applied to determine load share by each layer of the laminated element; 2. the layer is cut into a series of sheets each containing at most one straight fiber yarn segment, the bridging model is applied to determine the homogenized internal stresses and compliance tensor of each sheet, and the homogenized stresses are converted into true values as per the methods described subsequently; and 3. an assemblage from all of the sheets give the internal true stresses and compliance tensor of the layer.
Materials 16 00774 g004
Figure 5. RVE of a short fiber composite (a) is separated into the RVE of a continuous fiber composite (b) plus two end segments of a short fiber composite with a small fiber aspect ratio (c).
Figure 5. RVE of a short fiber composite (a) is separated into the RVE of a continuous fiber composite (b) plus two end segments of a short fiber composite with a small fiber aspect ratio (c).
Materials 16 00774 g005
Figure 6. Schematic of (a) UD composite subjected to a transverse tension, and (b) stable interface debonding due to Poisson’s retardant.
Figure 6. Schematic of (a) UD composite subjected to a transverse tension, and (b) stable interface debonding due to Poisson’s retardant.
Materials 16 00774 g006
Figure 7. Schematic diagram for relative slippage of a RVE subjected to an in-plane shear shown in (a) a longitudinal cross-section, (b) a transverse cross-section.
Figure 7. Schematic diagram for relative slippage of a RVE subjected to an in-plane shear shown in (a) a longitudinal cross-section, (b) a transverse cross-section.
Materials 16 00774 g007
Figure 8. In-plane shear stress–strain curve of E-glass/MY750 up to failure.
Figure 8. In-plane shear stress–strain curve of E-glass/MY750 up to failure.
Materials 16 00774 g008
Figure 9. (a) Longitudinal compression-induced kink-band failure, (b) load analysis in kink band coordinate system.
Figure 9. (a) Longitudinal compression-induced kink-band failure, (b) load analysis in kink band coordinate system.
Materials 16 00774 g009
Figure 10. Failure surface of a UD composite under a compression.
Figure 10. Failure surface of a UD composite under a compression.
Materials 16 00774 g010
Figure 11. Schematic of matrix compression-induced composite failure surface.
Figure 11. Schematic of matrix compression-induced composite failure surface.
Materials 16 00774 g011
Table 1. Matrix SCFs with other relevant parameters of the nine composite systems used in WWFEs.
Table 1. Matrix SCFs with other relevant parameters of the nine composite systems used in WWFEs.
E-Glass
LY556
E-Glass
MY750
AS4
3501-6
T300
BSL914C
IM7
8511-7
T300
PR319
AS
Epoxy
S2-Glass
Epoxy
G400-800
5260
K121.5201.4911.4241.4301.4751.5101.4491.5001.483
K233.0202.9361.3372.4212.0342.1671.9992.9822.469
K 22 t 3.3393.2532.0982.1432.3273.1232.3393.3172.464
K 22 c 2.2492.1811.4691.5701.7612.0351.7432.1721.732
K ^ 12 1.8701.8401.7601.7601.8301.7601.7601.8501.830
K ^ 22 t 7.6907.2204.9505.0405.4106.9705.4307.3405.680
ψ71.8071.9073.9073.9073.4072.0073.3071.8072.80
h/a0.6090.6660.5630.5740.6480.6730.5910.6800.627
K 11 t , I *3.6603.5743.7333.8774.4864.1883.9533.6884.024
K 11 t *1.5981.6751.5161.5881.9651.8441.6361.7261.720
K 11 c , I 3.5103.4193.5503.6854.2713.9953.7593.5283.833
K 11 c *1.5041.5771.4271.4961.8601.7571.5451.6311.630
* l/a = 2, l/L is determined from Equation (26) in which c = 0.03.
Table 2. Constituent properties and other parameters of the nine composite systems used in WWFEs [28,29,30].
Table 2. Constituent properties and other parameters of the nine composite systems used in WWFEs [28,29,30].
E-Glass
LY556
E-Glass
MY750
AS4
3501-6
T300
BSL914C
IM7
8511-7
T300
PR319
AS
Epoxy
S2-Glass
Epoxy
G400-800
5260
E 11 f (GPa)807422523027623023187290
E 22 f (GPa)807415151915158719
ν 12 f 0.20.20.20.20.20.20.20.20.2
G 12 f (GPa)33.3330.8151527151536.327
ν 23 f 0.20.20.070.070.360.070.070.20.357
σ u , t f (MPa)215021503350250051802500350028505860
σ u , c f (MPa)145014502500200032002000300024503200
Em (GPa)3.353.354.244.080.953.23.23.45
νm0.350.350.340.350.380.350.350.350.35
σ u , t m (MPa)808069759970857370
σ u , c m (MPa)120120250150130130120120130
σ u , s m (MPa)545450705741505257
Y (MPa)354048277340386375
Vf0.620.60.60.60.60.60.60.60.6
Table 3. Averaged errors between predictions and measurements for uniaxial strengths of the nine UD composites used in WWFEs.
Table 3. Averaged errors between predictions and measurements for uniaxial strengths of the nine UD composites used in WWFEs.
Longitudinal Tensile strengths Longitudinal Compressive StrengthsTransverse Tensile StrengthsTransverse Compressive StrengthsIn-Plane Shear StrengthsTransverse Shear StrengthsOverall Error
Prediction on homo-stresses11%25.1%241%65.5%48.1%105.4%82.7%
Prediction on true-stresses11%25.1%39.2%23.2%13.1%14.3%21.1%
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

Huang, Z.-M. True Stress Theory of Matrix in A Composite: A Topical Review. Materials 2023, 16, 774. https://doi.org/10.3390/ma16020774

AMA Style

Huang Z-M. True Stress Theory of Matrix in A Composite: A Topical Review. Materials. 2023; 16(2):774. https://doi.org/10.3390/ma16020774

Chicago/Turabian Style

Huang, Zheng-Ming. 2023. "True Stress Theory of Matrix in A Composite: A Topical Review" Materials 16, no. 2: 774. https://doi.org/10.3390/ma16020774

APA Style

Huang, Z. -M. (2023). True Stress Theory of Matrix in A Composite: A Topical Review. Materials, 16(2), 774. https://doi.org/10.3390/ma16020774

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop