Next Article in Journal
Guidelines for the LTO Noise Assessment of Future Civil Supersonic Aircraft in Conceptual Design
Previous Article in Journal
Real-Time Precise Orbit Determination for LEO between Kinematic and Reduced-Dynamic with Ambiguity Resolution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analytical Model for the Uniaxial Tensile Modulus of Plain-Woven Fabric Composites and Its Application and Experimental Validation

Department of Astronautics Science and Mechanics, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(1), 26; https://doi.org/10.3390/aerospace9010026
Submission received: 3 December 2021 / Revised: 29 December 2021 / Accepted: 30 December 2021 / Published: 4 January 2022

Abstract

:
With advantages in efficiency and convenience, analytical models using experimental inputs to predict the mechanical properties of plain-woven fabric (PWF) composites are reliable in guaranteeing the composites’ engineering applications. Considering the importance of the aspect above, a new analytical model for predicting the uniaxial tensile modulus of PWF is proposed in this article. The composite yarns are first simplified as the lenticular-shaped cross-sections undulate along arc-composed paths. Force analyses of the yarn segments are then carried out with the internal interactions simplified, and the analytical model is subsequently deduced from the principle of minimum potential energy and Castigliano’s second theorem. The PWF of T300/Cycom970 is chosen as the study object to which the proposed analytical model is applied. Microscopic observations and thermal ablation experiments are conducted on the specimens to obtain the necessary inputs. The uniaxial tensile modulus is calculated and tensile experiments on the laminates are performed to validate the analytical prediction. The small deviation between the experimental and analytical results indicates the feasibility of the proposed analytical model, which has good prospects in validating the effectiveness of the experimentally obtained modeling parameters and guaranteeing the accuracy of mesoscale modeling for the PWF.

1. Introduction

With the characteristics of lighter weight in addition to higher stiffness and strength, fiber-reinforced composites have been widely used in the aerospace industry. In order to reduce production costs and simplify the machining processes plain-woven fabric composites are becoming increasingly important in the design and manufacturing of lightweight aerospace structures.
With the aim of improving the level of structural design and analysis, it is necessary to understand the loading response and failure mechanism of plain-woven composites. Research has been widely carried out and various experimental techniques as well as methods have been utilized in the study of plain-woven composites. X-ray tomography and microscale measurement systems are important in measuring the 3D coordinates of plain-woven composites’ mesoscale structures and establishing mesoscale geometric models. Experimental observations based on the above methods [1,2] have contributed to a direct understanding of the composites’ woven morphology and offered enlightenment in the accurate modeling of the interlaced yarns. Wang et al. [3] used micro-CT equipment with a high-resolution lens-coupled detector to reconstruct 3D images of C/epoxy plain-woven composites. Based on the obtained clear 3D image of the T300/epoxy material, the statistical properties of the yarn feature parameters that were then analyzed. Scanning electron microscopes (SEMs) and digital imaging systems are also frequently used in the microscopic observations of composites. With the assistance of microscopic observations, Carvalho et al. [4] investigated the compressive failure of orthogonal 2D woven composites. The regional failure of the load-aligned tows was studied and the damage propagation process was analyzed in detail. Specialized fixtures have also been designed for different test targets. With a picture-frame apparatus Nguyen et al. [5] studied the shearing properties of plain-woven fabrics and prepregs. Experimental investigations on the shear properties of the fabric composites were conducted and the relationships between the shearing angle and the applied load for different weave architectures were analyzed.
Experimental results provide the basis on which plain-woven composites could be modeled. With the numerical models, the finite-element (FE) method has been widely applied in the analysis of plain-woven fabrics due to its advantages in calculation efficiency and low cost. In the FE modeling of plain-woven fabrics the industrial application of advanced meshing methods has attracted a growing interest due to its operational convenience [6]. The applicability of voxel meshes in the numerical modeling of plain-woven composites was studied by Doitrand et al. [7], and the calculation accuracy of the voxel method was evaluated. Grail et al. [8] developed an automated method to generate consistent FE meshes for plain-woven composites. Based on the proposed method, a representative unit cell (RUC) of a compacted and nested multilayer woven composite was generated.
Obtaining the mechanical properties of the plain-woven composites and evaluating the properties’ effects are necessary preconditions for the following structural analysis. With the assistance of RUC models and FE analyses, Goyal et al. [9] investigated the effects of geometric parameters and basic material properties on the effective engineering properties of plain-woven braids. Matveev et al. [10] numerically studied the effects of fiber strength on the macroscale strength of textile composites and found that a wide fiber strength distribution led to a narrow distribution of composite strength. In evaluating the mechanical properties’ effects and predicting the mechanical response of braided composites, Richter et al. [11] proposed a simulation framework with a meshing methodology that considered yarn-to-yarn delamination. The validation of the framework underlined the great potential for further damage modeling of the fabric composites.
The simulation of failure is also an important concern in the FE analysis of fabrics. In order to understand the elementary damage mechanisms and investigate the development of damage in detail, advanced techniques such as improved two-dimensional generalized plane strain FE models [12], the cracking kinetics simulation [13] and shear nonlinearity incorporated failure process modeling [14] have been utilized in investigations. Important damage details have been captured and the relationships between ply damage behaviors and elementary damage mechanisms have been linked.
Macro fabric composite structures could be further studied at the mesoscale of the fiber tows and the microscale of the constituents. In the multiscale failure analyses of the fabric composites micromechanics-based failure theories [15,16,17,18] are becoming increasingly appealing for researchers due to their advantage in explaining the failure of composites at the constituent level [19]. By combining macroscale FE models with RUCs of meso- and microscales the multiscale analysis method takes the advantages of both the efficiency of the macro analysis and the accuracy of the micromechanics-based failure criteria [20]. The effectiveness of the multiscale analysis in predicting the fabric composites’ damage initiation and failure process [21,22,23] indicated the huge potential of the micromechanics-based methods’ further engineering applications.
The mesoscale RUC is an important stress transformation medium of the plain-woven fabric composites’ multiscale analysis, and its accurate establishment is necessary for further failure analysis. Microscopic observations and basic tests could be conducted, making it convenient to obtain the geometric and mechanical parameters of the fabric composites. However, the feasibility of the experimentally obtained parameters in modeling the RUC needs further validation considering the dispersibility and errors in the experiments. In predicting the properties of textile composites, analytical solutions suitable for engineering applications have been proposed [24,25,26,27,28]. Although the theoretical foundations are solid the geometric parameters and internal force analyses have been oversimplified and idealized in previous research [24,28], and the prediction accuracy could be further raised with the improvement of the previous models. Using the fiber tows’ geometric and mechanical parameters as inputs and having an advantage in calculation efficiency, reliable analytical models have broad prospects in engineering applications in validating the effectiveness of the experimentally obtained modeling parameters and guaranteeing the accuracy of the established RUC.
In this article attempts have been made to develop a new reliable analytical model with experimental inputs in calculating the uniaxial tensile modulus of plain-woven fabric composites.
In the Section 2 the cross-section and undulation path of the fiber tows are reasonably simplified and geometrically expressed with parameters and formulations. Within the mesoscale RUC the studied yarn segments are assigned.
In the Section 3 the interactions between the warp and weft yarns are simplified and the internal forces of the studied yarn segments are analyzed in detail.
In the Section 4 the analytical model is developed. The analytical solution for the uniaxial tensile modulus is derived from the elastic deformation energy analysis with the fiber tows’ internal forces more comprehensively considered. The necessary parameters of inputs are referred and the matched experiments to obtain the parameters are designed.
Finally, in the Section 5 the proposed analytical model is applied to predict the uniaxial tensile modulus of T300/Cycom970, and the application process is introduced in detail. Tensile experiments on the T300/Cycom970 laminates are conducted and the effectiveness of the analytical model is validated.

2. Simplification of the Cross-Section and Undulation Path of the Fiber Tows

Due to the randomness of the woven prepregs and the influence of the hot-pressing manufacturing process there are difficulties in accurately describing the irregular-shaped cross-section of the fiber tows with certain parameters. Therefore, the simplification of the cross-section is needed as an important basis for the analytical solution of the fabric composites. Compared to the observations from computed tomography (CT), as illustrated in Figure 1a, previous research results [29,30] indicate that the lenticular shape best approaches the cross-section of the plain-woven composites’ fiber tows. Besides, the lenticular-shaped cross-section exceeds the other frequently used sections in ensuring the fiber tows of the warp and weft directions realistically contact without interpenetration. Considering the above advantages, the lenticular shape is chosen to simplify the cross-section of the fiber tows in this article.
The warp yarn is taken as the example and the cross-section is projected onto the Y–Z plane, as presented in Figure 1b. The controlling parameters in defining the lenticular-shaped cross-section are a, b and rT, which represent the width, thickness and controlling radius of the cross-section, respectively:
r T r T 2 1 4 a 2 = 1 2 b
In Figure 1b, z stands for the ordinate value of the Y–Z plane and L(z) stands for the z-corresponded horizontal intercept of the cross-section. The area, A, and the inertia moment of the cross-section, I, for the fiber tows in the warp and the weft directions are then expressed as follows:
A = A d A = 2 0 b / 2 L ( z ) d z
I = I x = I y = A z 2 d A = 2 0 b / 2 z 2 L ( z ) d z
Referring to Equation (1), the area and the inertia moment are further expressed as:
A = 4 0 b / 2 ( a 2 + b 2 4 b ) 2 [ ( a 2 + b 2 4 b ) 2 ( 1 2 a ) 2 + z ] 2 d z
I = 4 0 b / 2 z 2 ( a 2 + b 2 4 b ) 2 [ ( a 2 + b 2 4 b ) 2 ( 1 2 a ) 2 + z ] 2 d z
Corresponding to the lenticular-shaped cross-section, the undulation path composed of arcs with the same radius is used in this article. With the periodic characteristics of the plain-woven fabric composites, the analysis is performed on the mesoscale RUC, inside of which the segments of the warp and weft yarns are taken as the study objects, as illustrated in Figure 2.
The undulation path of the studied warp yarn segment is presented in Figure 3a. rU represents the radius of the arcs and LRUC represents the length of the mesoscale RUC. H stands for the distance between the horizontal center line of the RUC and the center of the arcs in the undulation path:
H = r U 2 ( 1 4 L R U C ) 2
For the studied weft yarn segment in Figure 2 the geometric parameters of the undulation path in Figure 3b remain the same, and the width of the mesoscale RUC (WRUC) is equal to LRUC, considering the plain-woven structure of the fabric composites.
In the x range of ( 0 , 1 2 L R U C ) the undulation path function z1(x) of the warp yarn segment corresponding to Figure 3a is expressed as:
z 1 ( x ) = { r U r U 2 x 2 , x = 0 ~ 1 4 L R U C r U 2 H + r U 2 ( x 1 2 L R U C ) 2 , x = 1 4 L R U C ~ 1 2 L R U C
Similarly, for the weft yarn segment in Figure 3b the undulation path function z2(y) is expressed as:
z 2 ( y ) = { r U 2 H + r U 2 y 2 , y = 0 ~ 1 4 W R U C r U r U 2 ( y 1 2 W R U C ) 2 , y = 1 4 W R U C ~ 1 2 W R U C

3. Force Analysis of the Fiber Tows under Uniaxial Tension Load

The monolayer of the plain-woven fabric composites under the uniaxial tension load along the warp direction is presented in Figure 4a. As the warp yarns are straightened the distributed interaction force of q between the weft and warp yarns exists in hindering the straightening of the warp yarns and causing the crimping of the weft yarns. For the convenience of the following analyses the distributed interaction force of q is simplified as the concentrate force of Q, as illustrated in Figure 4b.
For the studied mesoscale RUC of the plain-woven fabric composites in Figure 2, the force analyses of the warp and weft yarn segments are illustrated in Figure 5.
In Figure 5a the axial load of N and the constraint moment of M1 are loaded at the boundary cross-sections of the warp yarn segment. The interaction loads of Q and ½Q from the weft yarns concentrate at the middle location and the left and right endpoints, respectively. dl1 stands for the differential length of the warp yarn segment and θ stands for the angle between dl1 and the X-axis. At the x-corresponded cross-section in the warp yarn segment the horizontal internal load of N and the vertical internal load of ½Q, which satisfy the force balance condition, are further converted into the internal tensile load N1(x) and the internal shearing load Q1(x), which are perpendicular and parallel to the cross-section, respectively.
For the cross-sections of the studied warp yarn segment, within the range of (0, ½LRUC) along the X-axis, the internal bending moment M1(x) is expressed as:
M 1 ( x ) = M 1 + N z 1 ( x ) 1 2 Q x
The internal tensile load N 1 ( x ) is expressed as:
N 1 ( x ) = N cos θ + 1 2 Q sin θ = N 1 + tan 2 θ + Q tan θ 2 1 + tan 2 θ
Additionally, the internal shearing load Q 1 ( x ) is expressed as:
Q 1 ( x ) = N sin θ 1 2 Q cos θ = N tan θ 1 + tan 2 θ Q 2 1 + tan 2 θ
For the weft yarn segment in Figure 5b the constraint moments of M2 are loaded at the boundary cross-sections and the interaction loads of Q and ½Q from the warp yarns concentrate at the middle location and the left and right endpoints, respectively. The angle between the differential length of the weft yarn segment dl2 and the Y-axis is presented as φ. Similar to the warp yarn segment, the vertical internal load of 1 2 Q at the y-corresponded cross-section in the weft yarn segment is converted into the internal tensile load N2(y) and the internal shearing load Q2(y).
For the cross-sections of the studied weft yarn segment, within the range of (0, ½WRUC) along the Y-axis, the internal bending moment M2(y) is expressed as:
M 2 ( y ) = M 2 + 1 2 Q y
The internal tensile load N2(y) is expressed as:
N 2 ( y ) = 1 2 Q sin φ = Q tan φ 2 1 + tan 2 φ
Additionally, the internal shearing load Q2(y) is expressed as:
Q 2 ( y ) = 1 2 Q cos φ = Q 2 1 + tan 2 φ
Considering the arc-shaped undulation paths of the warp and weft yarns, the tangents of θ and φ are, respectively, expressed as:
tan θ = d z 1 d x = { x r U 2 x 2 , x = 0 ~ 1 4 L R U C 1 2 L R U C x r U 2 ( x 1 2 L R U C ) 2 , x = 1 4 L R U C ~ 1 2 L R U C
tan φ = d z 2 d y = { y r U 2 y 2 , y = 0 ~ 1 4 W R U C y 1 2 W R U C r U 2 ( y 1 2 W R U C ) 2 , y = 1 4 W R U C ~ 1 2 W R U C

4. Analytical Solution in Predicting the Uniaxial Tensile Modulus of the Plain-Woven Fabric Composites

The parameters of E11(tow), G12(tow) and G13(tow) stand for the longitudinal modulus, the in-plane shear modulus and the out-of-plane shear modulus of the plain-woven fabric composites’ fiber tows, respectively. Under the hypothesis of transverse isotropy G12(tow) is equal to G13(tow). L1 and L2 stand for the lengths of the warp and weft yarn segments illustrated in Figure 2, respectively. Under the hypothesis that the shearing stress is uniformly distributed on the fiber tow’s cross-section the strain energy of the studied warp yarn segment is presented as U1:
U 1 = 1 2 E 11 ( t o w ) I 0 L 1 M 1 ( x ) 2 d l 1 + 1 2 E 11 ( t o w ) A 0 L 1 N 1 ( x ) 2 d l 1 + 1 2 G 13 ( t o w ) A 0 L 1 Q 1 ( x ) 2 d l 1 = 1 2 E 11 ( t o w ) I 0 L R U C M 1 ( x ) 2 cos θ d x + 1 2 E 11 ( t o w ) A 0 L R U C N 1 ( x ) 2 cos θ d x + 1 2 G 13 ( t o w ) A 0 L R U C Q 1 ( x ) 2 cos θ d x
For the studied weft yarn segment the strain energy U2 is expressed as:
U 2 = 1 2 E 11 ( t o w ) I 0 L 2 M 2 ( y ) 2 d l 2 + 1 2 E 11 ( t o w ) A 0 L 2 N 2 ( y ) 2 d l 2 + 1 2 G 13 ( t o w ) A 0 L 2 Q 2 ( y ) 2 d l 2 = 1 2 E 11 ( t o w ) I 0 W R U C M 2 ( y ) 2 cos φ d y + 1 2 E 11 ( t o w ) A 0 W R U C N 2 ( y ) 2 cos φ d y + 1 2 G 13 ( t o w ) A 0 W R U C Q 2 ( y ) 2 cos φ d y
In the above equations the cosines of θ and φ are, respectively, expressed as:
cos θ = 1 1 + tan 2 θ = { r U 2 x 2 r U , x = 0 ~ 1 4 L R U C r U 2 ( x 1 2 L R U C ) 2 r U , x = 1 4 L R U C ~ 1 2 L R U C
cos φ = 1 1 + tan 2 φ = { r U 2 y 2 r U , y = 0 ~ 1 4 W R U C r U 2 ( y 1 2 W R U C ) 2 r U , y = 1 4 W R U C ~ 1 2 W R U C
The simplification functions of C1(x), C2(y), S1(x) and S2(y) are defined as:
C 1 ( x ) = r U 2 x 2 , C 2 ( y ) = r U 2 y 2
S 1 ( x ) = r U 2 ( x 1 2 L R U C ) 2 , S 2 ( y ) = r U 2 ( y 1 2 W R U C ) 2
Putting the internal loads of Equations (9)–(11) into the strain energy formulation in Equation (17), U1 could be further expressed with the simplified intermediate coefficients of X1, X2, X31, X5, X6 and X7 as follows, considering the symmetrical characteristics of the studied warp yarn segment:
U 1 = 1 E 11 ( t o w ) I { 0 L R U C 4 { M 1 + N [ r U C 1 ( x ) ] 1 2 Q x } 2 r U C 1 ( x ) d x + L R U C 4 L R U C 2 { M 1 + N [ r U 2 H + S 1 ( x ) ] 1 2 Q x } 2 r U S 1 ( x ) d x } + 1 E 11 ( t o w ) A { 0 L R U C 4 [ 2 N C 1 ( x ) + Q x ] 2 4 r U C 1 ( x ) d x + L R U C 4 L R U C 2 [ 2 N S 1 ( x ) Q ( x L R U C 2 ) ] 2 4 r U S 1 ( x ) d x } + 1 G 13 ( t o w ) A { 0 L R U C 4 [ 2 N x Q C 1 ( x ) ] 2 4 r U C 1 ( x ) d x + L R U C 4 L R U C 2 [ 2 N ( x L R U C 2 ) + Q S 1 ( x ) ] 2 4 r U S 1 ( x ) d x } = M 1 2 X 1 + N 2 X 2 + Q 2 X 31 + M 1 N X 5 M 1 Q X 6 N Q X 7
Similarly, for the studied weft yarn segment U2 could be further expressed with the simplified intermediate coefficients of X32, X4 and X8 as:
U 2 = 1 E 11 ( t o w ) I [ 0 W R U C 4 ( M 2 + 1 2 Q y ) 2 r U C 2 ( y ) d y + W R U C 2 W R U C 4 ( M 2 + 1 2 Q y ) 2 r U S 2 ( y ) d y ] + 1 E 11 ( t o w ) A [ 0 W R U C 4 Q 2 y 2 4 r U C 2 ( y ) d y + W R U C 4 W R U C 2 Q 2 ( y 1 2 W R U C ) 2 4 r U S 2 ( y ) d y ] + 1 G 13 ( t o w ) A [ 0 W R U C 4 Q 2 C 2 ( y ) 4 r U d y + W R U C 4 W R U C 2 Q 2 S 2 ( y ) 4 r U d y ] = Q 2 X 32 + M 2 2 X 4 + M 2 Q X 8
The simplified intermediate coefficients in the above equations of (23) and (24) are expressed as follows:
X 1 = 0 L R U C 4 r U E 11 ( t o w ) I C 1 ( x ) d x + L R U C 4 L R U C 2 r U E 11 ( t o w ) I S 1 ( x ) d x
X 2 = 0 L R U C 4 { r U [ r U C 1 ( x ) ] 2 E 11 ( t o w ) I C 1 ( x ) + C 1 ( x ) E 11 ( t o w ) A r U + x 2 G 13 ( t o w ) A r U C 1 ( x ) } d x + L R U C 4 L R U C 2 { r U [ r U 2 H + S 1 ( x ) ] 2 E 11 ( t o w ) I S 1 ( x ) + S 1 ( x ) E 11 ( t o w ) A r U + ( x 1 2 L R U C ) 2 G 13 ( t o w ) A r U S 1 ( x ) } d x
X 31 = 0 L R U C 4 [ r U x 2 4 E 11 ( t o w ) I C 1 ( x ) + x 2 4 E 11 ( t o w ) A r U C 1 ( x ) + C 1 ( x ) 4 G 13 ( t o w ) A r U ] d x + L R U C 4 L R U C 2 [ r U x 2 4 E 11 ( t o w ) I S 1 ( x ) + ( x 1 2 L R U C ) 2 4 E 11 ( t o w ) A r U S 1 ( x ) + S 1 ( x ) 4 G 13 ( t o w ) A r U ] d x
X 32 = 0 W R U C 4 [ r U y 2 4 E 11 ( t o w ) I C 2 ( y ) + y 2 4 E 11 ( t o w ) A r U C 2 ( y ) + C 2 ( y ) 4 G 13 ( t o w ) A r U ] d y + W R U C 4 W R U C 2 [ r U y 2 4 E 11 ( t o w ) I S 2 ( y ) + ( y 1 2 W R U C ) 2 4 E 11 ( t o w ) A r U S 2 ( y ) + S 2 ( y ) 4 G 13 ( t o w ) A r U ] d y
X 4 = 0 W R U C 4 r U E 11 ( t o w ) I C 2 ( y ) d y + W R U C 4 W R U C 2 r U E 11 ( t o w ) I S 2 ( y ) d y
X 5 = 0 L R U C 4 2 r U [ r U C 1 ( x ) ] E 11 ( t o w ) I C 1 ( x ) d x + L R U C 4 L R U C 2 2 r U [ r U 2 H + S 1 ( x ) ] E 11 ( t o w ) I S 1 ( x ) d x
X 6 = 0 L R U C 4 r U x E 11 ( t o w ) I C 1 ( x ) d x + L R U C 4 L R U C 2 r U x E 11 ( t o w ) I S 1 ( x ) d x
X 7 = 0 L R U C 4 { r U x [ r U C 1 ( x ) ] E 11 ( t o w ) I C 1 ( x ) x E 11 ( t o w ) A r U + x G 13 ( t o w ) A r U } d x + L R U C 4 L R U C 2 { r U x [ r U 2 H + S 1 ( x ) ] E 11 ( t o w ) I S 1 ( x ) + x 1 2 L R U C E 11 ( t o w ) A r U x 1 2 L R U C G 13 ( t o w ) A r U } d x
X 8 = 0 W R U C 4 r U y E 11 ( t o w ) I C 2 ( y ) d y + W R U C 4 W R U C 2 r U y E 11 ( t o w ) I S 2 ( y ) d y
By defining the intermediate coefficient of X3 as:
X 3 = X 31 + X 32
the total strain energy, U, for the part of the fiber tows in the mesoscale RUC of Figure 2 could then be obtained as:
U = 2 ( U 1 + U 2 ) = 2 ( M 1 2 X 1 + N 2 X 2 + Q 2 X 3 + M 2 2 X 4 + M 1 N X 5 M 1 Q X 6 N Q X 7 + M 2 Q X 8 )
Constraint at the left end, the right end’s displacement of the deformed warp yarn segment along the X-axis corresponding to the axial tensile load of N is presented as u1 in Figure 6.
For the studied warp yarn segment in Figure 2 U1* stands for its complementary energy. According to Castigliano’s second theorem:
u 1 = U 1 * N
Within the framework of the ideal linear elasticity U1* is equal to the strain energy U1, and u1 in Figure 6 could be expressed as:
u 1 = U 1 N = 0 L R U C M 1 ( x ) E 11 ( t o w ) I cos θ M 1 ( x ) N d x + 0 L R U C N 1 ( x ) E 11 ( t o w ) A cos θ N 1 ( x ) N d x + 0 L R U C Q 1 ( x ) G 13 ( t o w ) A cos θ Q 1 ( x ) N d x = 2 X 2 N + X 5 M 1 X 7 Q
The constraint moments of M1 and M2 and the interaction load of Q are all introduced by N , which is the only external load applied to warp yarn in the monolayer of the plain-woven fabric composites. Under the linear hypothesis within the small deformation extents the constraint moments and the interaction load are supposed linearly correlated with the external load, N, through the coefficients of k1, k2 and k3:
M 1 = k 1 N
M 2 = k 2 N
Q = k 3 N
Putting Equations (38)–(40) into Equation (35), the total strain energy of U then turns into:
U = 2 ( k 1 2 X 1 + X 2 + k 3 2 X 3 + k 2 2 X 4 + k 1 X 5 k 1 k 3 X 6 k 3 X 7 + k 2 k 3 X 8 ) N 2
For the part of fiber tows in the mesoscale RUC of Figure 2, the total potential energy is defined as Π and the external work is defined as W. Considering the total strain energy of U, Π could be expressed as:
Π = U W
According to the principle of minimum potential energy, among all the displacement solutions which satisfy the geometric boundary conditions the real displacement solution minimizes the structural total potential energy. Corresponding to the minimum total potential energy, Πmin, under real structural displacement the reasonable choice of the coefficients k1, k2 and k3 minimizes the total strain energy, U, which leads to:
U k i = 0 , i = 1 ~ 3
The coefficients of k1, k2 and k3 are then determined as:
k 1 = X 7 X 6 2 X 1 X 5 X 6 2 4 X 1 2 2 X 3 X 6 2 2 X 1 X 8 2 2 X 4 X 5 2 X 1
k 2 = ( X 5 X 6 2 X 1 X 7 ) X 8 ( 4 X 3 X 6 2 X 1 X 8 2 X 4 ) X 4
k 3 = X 7 X 5 X 6 2 X 1 2 X 3 X 6 2 2 X 1 X 8 2 2 X 4
The relationships in Equations (38)–(40) with k1, k2 and k3 in Equations (44)–(46) could also be deduced from the principle of the minimum complementary energy as presented in Equations (47)–(49), indicating the rationality of the linear hypothesis:
U 1 * M 1 = 0
U 1 * M 2 = 0
U 1 * Q = 0
Considering Equations (44)–(46) and putting Equations (38)–(40) into Equation (37), the relationship between u1 and N in Figure 6 could be further determined with the intermediate coefficient of X9:
u 1 = X 9 N
where:
X 9 = 2 X 2 X 5 2 2 X 1 + 2 X 5 X 6 X 7 2 X 1 X 7 2 X 5 2 X 6 2 2 X 1 4 X 1 X 3 X 6 2 X 1 X 8 2 X 4
The part composed of fiber tows in the mesoscale RUC is isolated and studied in Figure 7. With the left-end constraint the displacement of the loading end reaches u1 under the total axial tension of 2N loaded at the cross-sections of the warp yarns with the area of 2A, corresponding to the result of the single warp yarn segment in Figure 6.
For the studied part in Figure 7, define the equivalent modulus along the X-axis as Ett and u1 could be expressed as:
u 1 = N L R U C E t t A
Put Equation (50) into Equation (52) and the modulus of Ett is determined as:
E t t = L R U C X 9 A
Em stands for the elastic modulus of the matrix and Vt stands for the volume fraction of the tows in the plain-woven fabric composites. Based on the rule of mixtures, the uniaxial tensile modulus, E11, of the plain-woven fabric composites could be expressed as:
E 11 = E t t V t + E m ( 1 V t ) = L R U C X 9 A V t + E m ( 1 V t )
With the above deductions, the proposed analytical model for the uniaxial tensile modulus of the plain-woven fabric composites is illustrated in Figure 8. Based on the experimental observations, the necessary geometric parameters of the fiber tows and the RUC, including a, b, rU and LRUC, are first obtained as inputs. The corresponding cross-section area of A and the inertia moment of I are subsequently calculated. Series of experiments including the thermal ablation of the laminates and microscopic observations are needed in order to achieve the fiber volume fraction in tows, Vf, and the volume fraction of tows, Vt. Besides, parameters of the elastic modulus of the matrix, Em, the shear modulus of the matrix, Gm, the in-plane shear modulus of the fiber, Gf12, and the longitudinal modulus of the fiber, Ef11, need to be obtained, with which the longitudinal modulus of E11(tow) and the out-of-plane shear modulus of G13(tow) for the fiber tows could be calculated through the Chamis equations:
E 11 ( t o w ) = V f E f 11 + ( 1 V f ) E m
G 13 ( t o w ) = G m 1 V f ( 1 G m G f 12 )
Subsequently, the intermediate coefficients of X 1 X 9 could be achieved and the uniaxial tensile modulus of E 11 is analytically solved. Using the platform of an Intel Core i7-9750H CPU, the single run time of the analytical model in calculating E 11 is 2 s with the necessary inputs prepared.

5. Application of the Analytical Model and the Experimental Validation

The plain-woven fabric composite of T300/Cycom970 is chosen as the study object in this article and the uniaxial tensile modulus of E11 is predicted with the method proposed above. The mechanical properties of the T300 carbon fiber and the Cycom970 epoxy matrix [31] are presented in Table 1.
Specimens of T300/Cycom970 laminates with a stacking sequence of [0]16 and dimensions of 25 mm × 25 mm × 3.456 mm are manufactured and scanned first with micro-CT equipment to achieve the necessary geometric parameters. As a carbon fiber/epoxy composite, the similar atomic numbers of the fiber and the matrix in the T300/Cycom970 laminates lead to a small difference in the components’ X-ray absorption rates. In order to solve the blurry imaging problem micro-CT equipment with a high-resolution lens-coupled detector is utilized in the current scanning. A typical Y–Z slice image of the fabric laminate is illustrated in Figure 9 and the statistical average values of a, b, rU and LRUC from different samples of sections are presented in Table 2, together with the corresponding A and I.
Thermal ablation experiments have been conducted on the specimens under the ASTM Standard D3171 [32]-based procedure in order to separate the component of the fiber from the matrix and obtain the total fiber volume fraction of VfT. The specimens are first weighed with an electronic balance and the initial weight of m0 for each specimen is achieved. After heating with a chamber furnace under a temperature of 600 °C for four hours, as shown in Figure 10, the epoxy matrix component in the specimens is totally ablated and the residual weight of fibers, mf, for each specimen is recorded. For each of the specimens, with the densities of the matrix (ρm) and the carbon fiber (ρf) presented in Table 1, the total fiber volume fraction, VfT, of the plain-woven fabric composite could be obtained as:
V f T = m f ρ m m f ρ m + ( m 0 m f ) ρ f
The results of the thermal ablation experiments on the specimens, numbered one–five, are presented in Table 3. For the studied composite of T300/Cycom970, VfT is determined as the statistical average value of the experimental results.
In obtaining the fiber volume fraction in tows of Vf, microscopic observations have been conducted on the plain-woven fabric composites. In order to achieve ideal observation results focused ion beam (FIB) polishing has been conducted on the specimen, and the regions of the fiber tows’ cross-sections are then observed with a scanning electron microscope (SEM), as shown in Figure 11. Within a fixed observation area the regions of the fiber sections are distinguished from the adjacent matrix regions through further digital image processing of the observation result. A f and A m , which represent the area of the fibers and the area of the matrix, respectively, are achieved, and the regional fiber volume fraction, which is equal to the area ratio of A f / ( A f + A m ) , is then calculated. A typical microscopic observation on the cross-section of the fiber tows and the corresponding digital image processing result are presented in Figure 11, and the value of V f for the plain-woven fabric composite is determined as the statistical average value of the obtained different regional fiber volume fractions.
According to the Equation (58) the volume fraction of tows, Vt, of the T300/Cycom970 fabric composite could be obtained with VfT and Vf:
V t = V f T V f
The finally determined values of V f T , V f and V t are listed in Table 4. Based on the fiber volume fraction in tows of V f and the constituents’ properties in Table 1, the calculations of the fiber tows’ mechanical properties based on the Chamis equations are presented in Table 5.
With the above inputs the intermediate coefficients of X1X9 are calculated according to the flowchart illustrated in Figure 8, and the result of E11 for the T300/Cycom970 plain-woven fabric composite is predicted as 54.39 MPa.
Tensile experiments have been performed on the laminates made of T300/Cycom970 plain-woven fabrics. The experimental result is compared with the above prediction of E11 in order to validate the proposed analytical model.
Considering the purpose of validation, the manufacturing process of the composite laminates has been intentionally controlled to enhance the specimens’ consistency and reduce the effects of uncertainty. The specimens with a stacking sequence of [0]16 and dimensions of 250 mm × 25 mm × 3.456 mm are shown in Figure 12, which are tested by an electronic universal testing machine with a maximum loading capacity of 100 kN. Six specimens numbered CL-1-01 to CL-1-06 have been tested under the ASTM Standard D3039M-17 [33]-based procedure. Subjected to tensile loading, the stress–strain curves of the specimens have been achieved, with longitudinal strain results obtained from the strain gages. For two different points on the stress–strain curve, the difference in the applied tensile stress is defined as Δσ and the difference between the strains is defined as Δε. Stress–strain data within the lower half of the stress–strain curves are utilized ito calculate the uniaxial tensile modulus of E11:
E 11 = Δ σ Δ ε
In Table 6 the experimental E11 of the specimens is listed, along with the average value, the standard deviation and the coefficient of dispersion. The small values of the standard deviation and the coefficient of dispersion indicate the good repeatability of the experimental results. With neglectable effects of uncertainty, the specimens behave well in the aspect of validation. Compared to the average value of the experimental E11 the calculation result from the analytical model of 54.39 MPa has a deviation of −4.33%. The small deviation validates the effectiveness of the analytical model in predicting the unidirectional tensile modulus of the plain-woven fabric composites and further indicates the rationality of the inputs of the experimentally obtained geometric and mechanical parameters, which could be utilized in the establishment of the composites’ mesoscale RUC.

6. Conclusions

(1) With a solid theoretical foundation and validated by tensile experiments on the laminates, the analytical model based on the energy principles has a small margin of error and a high calculation efficiency. The model provides a new and effective method to predict the uniaxial tensile modulus of plain-woven fabric composites.
(2) The matched experiments were designed to obtain inputs for the analytical model and were conducted to predict the uniaxial tensile modulus of the T300/Cycom970 composite. The experimental procedure provides a reference for the further study of PWF composites’ geometric and mechanical properties.
(3) The proposed analytical model relies heavily on the inputs obtained from experimental methods; the costly experiments limit the wide application of the analytical model to some extent. Besides, the cross-section of fiber tows and the interaction between the warp and weft yarns are simplified in the analytical model for the convenience of the analyses. Although the current margin of error is experimentally validated a more comprehensive and detailed consideration of the fiber tows’ geometric characteristics and internal interactions is still needed in further research.
(4) With the advantages in efficiently validating the choosing of the fiber tows’ geometric and mechanical parameters, the analytical model has good application prospects in guiding the mesoscale modeling of plain-woven fabric composites. In future work the mesoscale RUC will be built on the current basis and multiscale failure analyses could be further carried out on the T300/Cycom970 PWF composite laminates and structures.

Author Contributions

Conceptualization, R.Z. and W.G.; methodology, W.L. and J.X.; software, R.Z.; validation, R.Z., W.L. and W.G.; formal analysis, R.Z.; investigation, R.Z.; resources, W.G. and J.X.; data curation, W.L.; writing—original draft preparation, R.Z.; writing—review and editing, W.G.; visualization, W.L.; supervision, R.Z.; project administration, R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The first author would like to thank Xing Li for the thermal ablation experiments in this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Naouar, N.; Vidal-Sallé, E.; Schneider, J.; Maire, E.; Boisse, P. Meso-scale FE analyses of textile composite reinforcement deformation based on X-ray computed tomography. Compos. Struct. 2014, 116, 165–176. [Google Scholar] [CrossRef]
  2. Zheng, T.; Zhang, X.; Zhao, Z.; Wu, Z.; Li, T. Geometric structure model of plain woven fabric based on progressive spring-slide mechanics. Text. Res. J. 2014, 84, 1803–1819. [Google Scholar] [CrossRef]
  3. Wang, H.; Wang, Z.-W. Statistical Analysis of Yarn Feature Parameters in C/Epoxy Plain-Weave Composite Using Micro CT with High-Resolution Lens-Coupled Detector. Appl. Compos. Mater. 2016, 23, 601–622. [Google Scholar] [CrossRef]
  4. De Carvalho, N.; Pinho, S.; Robinson, P. An experimental study of failure initiation and propagation in 2D woven composites under compression. Compos. Sci. Technol. 2011, 71, 1316–1325. [Google Scholar] [CrossRef] [Green Version]
  5. Nguyen, M.; Herszberg, I.; Paton, R. The shear properties of woven carbon fabric. Compos. Struct. 1999, 47, 767–779. [Google Scholar] [CrossRef]
  6. Kim, H.J.; Swan, C.C. Voxel-based meshing and unit-cell analysis of textile composites. Int. J. Numer. Methods Eng. 2003, 56, 977–1006. [Google Scholar] [CrossRef]
  7. Doitrand, A.; Fagiano, C.; Irisarri, F.-X.; Hirsekorn, M. Comparison between voxel and consistent meso-scale models of woven composites. Compos. Part A Appl. Sci. Manuf. 2015, 73, 143–154. [Google Scholar] [CrossRef]
  8. Grail, G.; Hirsekorn, M.; Wendling, A.; Hivet, G.; Hambli, R. Consistent Finite Element mesh generation for meso-scale modeling of textile composites with preformed and compacted reinforcements. Compos. Part A Appl. Sci. Manuf. 2013, 55, 143–151. [Google Scholar] [CrossRef]
  9. Goyal, D.; Tang, X.; Whitcomb, J.D.; Kelkar, A.D. Effect of Various Parameters on Effective Engineering Properties of 2 × 2 Braided Composites. Mech. Adv. Mater. Struct. 2005, 12, 113–128. [Google Scholar] [CrossRef]
  10. Matveev, M.; Long, A.; Jones, I. Modelling of textile composites with fibre strength variability. Compos. Sci. Technol. 2014, 105, 44–50. [Google Scholar] [CrossRef]
  11. Wehrkamp-Richter, T.; De Carvalho, N.V.; Pinho, S. A meso-scale simulation framework for predicting the mechanical response of triaxial braided composites. Compos. Part A Appl. Sci. Manuf. 2018, 107, 489–506. [Google Scholar] [CrossRef]
  12. Le Page, B.; Guild, F.; Ogin, S.; Smith, P. Finite element simulation of woven fabric composites. Compos. Part A Appl. Sci. Manuf. 2004, 35, 861–872. [Google Scholar] [CrossRef]
  13. Obert, E.; Daghia, F.; Ladevèze, P.; Ballere, L. Micro and meso modeling of woven composites: Transverse cracking kinetics and homogenization. Compos. Struct. 2014, 117, 212–221. [Google Scholar] [CrossRef]
  14. Zhou, Y.; Lu, Z.; Yang, Z. Progressive damage analysis and strength prediction of 2D plain weave composites. Compos. Part B Eng. 2013, 47, 220–229. [Google Scholar] [CrossRef]
  15. Gosse, J.; Christensen, S. Strain invariant failure criteria for polymers in composite materials. In Proceedings of the 19th AIAA Applied Aerodynamics Conference, Fluid Dynamics and Co-located Conferences, Anaheim, CA, USA, 11–14 June 2001. [Google Scholar] [CrossRef]
  16. Mayes, J.; Hansen, A.C. Composite laminate failure analysis using multicontinuum theory. Compos. Sci. Technol. 2004, 64, 379–394. [Google Scholar] [CrossRef]
  17. Ha, S.K.; Jin, K.K.; Huang, Y. Micro-Mechanics of Failure (MMF) for Continuous Fiber Reinforced Composites. J. Compos. Mater. 2008, 42, 1873–1895. [Google Scholar] [CrossRef]
  18. Ha, S.K.; Huang, Y.; Han, H.H.; Jin, K.K. Micromechanics of Failure for Ultimate Strength Predictions of Composite Laminates. J. Compos. Mater. 2010, 44, 2347–2361. [Google Scholar] [CrossRef]
  19. Li, X.; Guan, Z.; Li, Z.; Liu, L. A new stress-based multi-scale failure criterion of composites and its validation in open hole tension tests. Chin. J. Aeronaut. 2014, 27, 1430–1441. [Google Scholar] [CrossRef] [Green Version]
  20. Lou, X.; Han, X.; Cai, H. A high-efficient model for interface debonding analysis of carbon fiber-reinforced polymer composite. Compos. Interfaces 2019, 27, 87–109. [Google Scholar] [CrossRef]
  21. Zhao, C.; Huang, Y.; Chen, Z.; Ha, S.K. Progressive failure prediction of a landing gear structure of braided composites. Compos. Struct. 2017, 161, 407–418. [Google Scholar] [CrossRef]
  22. Xu, L.; Jin, C.Z.; Ha, S.K. Ultimate strength prediction of braided textile composites using a multi-scale approach. J. Compos. Mater. 2014, 49, 477–494. [Google Scholar] [CrossRef]
  23. Wang, M.; Zhang, P.; Fei, Q.; Guo, F. Modified micro-mechanics based multiscale model for progressive failure prediction of 2D twill woven composites. Chin. J. Aeronaut. 2020, 33, 2070–2087. [Google Scholar] [CrossRef]
  24. Xiong, J.; Shenoi, R.; Cheng, X. A modified micromechanical curved beam analytical model to predict the tension modulus of 2D plain weave fabric composites. Compos. Part B Eng. 2009, 40, 776–783. [Google Scholar] [CrossRef]
  25. Sheng, S.Z.; Van Hoa, S. Three Dimensional Micro-Mechanical Modeling of Woven Fabric Composites. J. Compos. Mater. 2001, 35, 1701–1729. [Google Scholar] [CrossRef]
  26. Sun, H.; Pan, N. Shear deformation analysis for woven fabrics. Compos. Struct. 2005, 67, 317–322. [Google Scholar] [CrossRef] [Green Version]
  27. Basit, M.M.; Luo, S.-Y. A simplified model of plain weave fabric reinforcements for the pure shear loading. Int. J. Mater. Form. 2018, 11, 445–453. [Google Scholar] [CrossRef]
  28. Cheng, X.; Xiong, J.; Bai, J. Analytical Solution for Predicting In-plane Elastic Shear Properties of 2D Orthogonal PWF Composites. Chin. J. Aeronaut. 2012, 25, 575–583. [Google Scholar] [CrossRef] [Green Version]
  29. Nguyen, Q.; Vidal-Sallé, E.; Boisse, P.; Park, C.H.; Saouab, A.; Bréard, J.; Hivet, G. Mesoscopic scale analyses of textile composite reinforcement compaction. Compos. Part B Eng. 2013, 44, 231–241. [Google Scholar] [CrossRef]
  30. Xu, L.; Huang, Y.; Zhao, C.; Ha, S.K. Progressive failure prediction of woven fabric composites using a multi-scale approach. Int. J. Damage Mech. 2016, 27, 97–119. [Google Scholar] [CrossRef]
  31. Tran, T.; Kelly, D.; Prusty, B.G.; Pearce, G.; Gosse, J. A micromechanical sub-modelling technique for implementing Onset Theory. Compos. Struct. 2013, 103, 1–8. [Google Scholar] [CrossRef]
  32. ASTM Standard D3171-15. Standard Test Methods for Constituent Content of Composite Materials; ASTM International: West Conshohocken, PA, USA, 2015. [Google Scholar] [CrossRef]
  33. ASTM Standard D3039/D3039M—17. Standard Test Method for Tensile Properties of Polymer Matrix Composite Materials; ASTM International: West Conshohocken, PA, USA, 2017. [Google Scholar] [CrossRef]
Figure 1. The lenticular-shaped cross-section of the fiber tows and the corresponding geometric parameters: (a) experimental CT image of the fiber tow [29]; (b) geometric parameters of the lenticular-shaped cross-section.
Figure 1. The lenticular-shaped cross-section of the fiber tows and the corresponding geometric parameters: (a) experimental CT image of the fiber tow [29]; (b) geometric parameters of the lenticular-shaped cross-section.
Aerospace 09 00026 g001
Figure 2. RUC of the plain-woven fabric composites and the studied warp and weft yarn segments.
Figure 2. RUC of the plain-woven fabric composites and the studied warp and weft yarn segments.
Aerospace 09 00026 g002
Figure 3. The undulation paths of the studied yarn segments: (a) the undulation path of the warp yarn segment; (b) the undulation path of the weft yarn segment.
Figure 3. The undulation paths of the studied yarn segments: (a) the undulation path of the warp yarn segment; (b) the undulation path of the weft yarn segment.
Aerospace 09 00026 g003
Figure 4. Interactions between the fiber tows: (a) the monolayer of the plain-woven fabric composites under uniaxial tension load; (b) simplification of the interactions between the warp and weft tows.
Figure 4. Interactions between the fiber tows: (a) the monolayer of the plain-woven fabric composites under uniaxial tension load; (b) simplification of the interactions between the warp and weft tows.
Aerospace 09 00026 g004
Figure 5. Force analyses of the studied yarn segments: (a) force analysis of the warp yarn segment; (b) force analysis of the weft yarn segment.
Figure 5. Force analyses of the studied yarn segments: (a) force analysis of the warp yarn segment; (b) force analysis of the weft yarn segment.
Aerospace 09 00026 g005
Figure 6. X-directional deformation of the deformed warp yarn segment.
Figure 6. X-directional deformation of the deformed warp yarn segment.
Aerospace 09 00026 g006
Figure 7. The part composed of fiber tows in the mesoscale RUC of the plain-woven fabric composites.
Figure 7. The part composed of fiber tows in the mesoscale RUC of the plain-woven fabric composites.
Aerospace 09 00026 g007
Figure 8. Flowchart of the analytical model for the uniaxial tensile modulus of the plain-woven fabric composites.
Figure 8. Flowchart of the analytical model for the uniaxial tensile modulus of the plain-woven fabric composites.
Aerospace 09 00026 g008
Figure 9. A typical Y–Z slice image of the T300/Cycom970 plain-woven fabric composite.
Figure 9. A typical Y–Z slice image of the T300/Cycom970 plain-woven fabric composite.
Aerospace 09 00026 g009
Figure 10. Thermal ablation experiments on the specimens of T300/Cycom970 fabric laminates.
Figure 10. Thermal ablation experiments on the specimens of T300/Cycom970 fabric laminates.
Aerospace 09 00026 g010
Figure 11. Microscopic observation of the specimen of plain-woven fabric composites and digital image processing for the SEM observation result.
Figure 11. Microscopic observation of the specimen of plain-woven fabric composites and digital image processing for the SEM observation result.
Aerospace 09 00026 g011
Figure 12. Tensile specimens of the T300/Cycom970 fabric laminates.
Figure 12. Tensile specimens of the T300/Cycom970 fabric laminates.
Aerospace 09 00026 g012
Table 1. Mechanical properties of the fiber and matrix constituents in the T300/Cycom970 composite.
Table 1. Mechanical properties of the fiber and matrix constituents in the T300/Cycom970 composite.
ConstituentMaterial PropertyValue
Fiber (T300)Longitudinal modulus, Ef11 (GPa)230.00
Transverse modulus, Ef22 = Ef33 (GPa)13.80
In-plane shear modulus, Gf12 (GPa)8.97
Out-of-plane shear modulus, Gf13 (GPa)8.97
Out-of-plane shear modulus, Gf23 (GPa)4.38
In-plane Poisson’s ratio, νf120.20
Out-of-plane Poisson’s ratio, νf130.20
Out-of-plane Poisson’s ratio, νf230.25
Density, ρf (g/cm3)1.76
Matrix (Cycom970)Elastic modulus, Em (GPa)3.35
Elastic Poisson’s ratio, νm0.35
Density, ρm (g/cm3)1.24
Table 2. Geometric parameters of fiber tows in the T300/Cycom970 plain-woven fabric composite.
Table 2. Geometric parameters of fiber tows in the T300/Cycom970 plain-woven fabric composite.
Geometric ParameterValue
Width of the fiber tows’ lenticular cross-section, a (mm)1.82
Thickness of the fiber tows’ lenticular cross-section, b (mm)0.13
Radius of the arcs in fiber tows’ undulation path, rU (mm)7.18
Length of the mesoscale RUC, LRUC (mm)4.00
Area of the fiber tows’ lenticular cross-section, A (mm2)0.16
Inertia of the fiber tows’ lenticular cross-section, I (10−4 mm4)1.67
Table 3. Results of the thermal ablation experiments.
Table 3. Results of the thermal ablation experiments.
SpecimenInitial Weight m 0 (g)Residual Weight m f (g)Total Fiber Volume Fraction V f T
No. 13.212.070.56
No. 23.232.090.56
No. 33.222.040.55
No. 43.202.050.56
No. 53.202.050.56
Table 4. Volume fraction parameters of the T300/Cycom970 plain-woven fabric composite.
Table 4. Volume fraction parameters of the T300/Cycom970 plain-woven fabric composite.
ParameterValue
Total fiber volume fraction, VfT0.56
Fiber volume fraction in tows, Vf0.78
Volume fraction of tows, Vt0.72
Table 5. Mechanical properties of the fiber tows in the T300/Cycom970 plain-woven fabric composite.
Table 5. Mechanical properties of the fiber tows in the T300/Cycom970 plain-woven fabric composite.
PropertyValue
Longitudinal modulus, E11(tow) (GPa)180.24
Transverse modulus, E22(tow) = E33(tow) (GPa)10.27
In-plane shear modulus, G12(tow) (GPa)9.66
Out-of-plane shear modulus, G13(tow) (GPa)9.66
Out-of-plane shear modulus, G23(tow) (GPa)3.78
Table 6. Experimental uniaxial tensile modulus of the T300/Cycom970 fabric laminates and comparison with the analytical solution.
Table 6. Experimental uniaxial tensile modulus of the T300/Cycom970 fabric laminates and comparison with the analytical solution.
SpecimenExperimental E11 (GPa)Average Value (GPa)Standard Deviation (GPa)Coefficient of Dispersion (%)Analytical Solution of E11 (GPa)Relative Deviation (%)
CL-1-0157.3556.850.721.2754.39−4.33
CL-1-0257.21
CL-1-0356.52
CL-1-0456.21
CL-1-0555.85
CL-1-0657.96
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, R.; Gao, W.; Liu, W.; Xu, J. An Analytical Model for the Uniaxial Tensile Modulus of Plain-Woven Fabric Composites and Its Application and Experimental Validation. Aerospace 2022, 9, 26. https://doi.org/10.3390/aerospace9010026

AMA Style

Zhou R, Gao W, Liu W, Xu J. An Analytical Model for the Uniaxial Tensile Modulus of Plain-Woven Fabric Composites and Its Application and Experimental Validation. Aerospace. 2022; 9(1):26. https://doi.org/10.3390/aerospace9010026

Chicago/Turabian Style

Zhou, Rui, Weicheng Gao, Wei Liu, and Jianxun Xu. 2022. "An Analytical Model for the Uniaxial Tensile Modulus of Plain-Woven Fabric Composites and Its Application and Experimental Validation" Aerospace 9, no. 1: 26. https://doi.org/10.3390/aerospace9010026

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