Next Article in Journal
The Mechanism of Dendrite Formation in a Solid-State Transformation of High Aluminum Fe-Al Alloys
Next Article in Special Issue
Feasibility of Stress Wave-Based Debond Defect Detection for RCFSTs Considering the Influence of Randomly Distributed Circular Aggregates with Mesoscale Homogenization Methodology
Previous Article in Journal
Fracturing Behaviors and Mechanism of Serial Coal Pillar Specimens with Different Strength
Previous Article in Special Issue
Molecular Dynamics Simulation on Nanoindentation of M50 Bearing Steel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiscale Study of the Effect of Fiber Twist Angle and Interface on the Viscoelasticity of 2D Woven Composites

1
School of Mechanical and Automotive Engineering, Shanghai University of Engineering Science, Shanghai 201620, China
2
Department of Civil Engineering, Zhejiang College of Construction, Hangzhou 311231, China
3
School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(7), 2689; https://doi.org/10.3390/ma16072689
Submission received: 5 February 2023 / Revised: 11 March 2023 / Accepted: 26 March 2023 / Published: 28 March 2023

Abstract

:
Time and temperature affect the viscoelasticity of woven composites, and thus affect their long-term mechanical properties. We develop a multiscale method considering fiber twist angle and interfaces to predict viscoelasticity. The multiscale approach is based on homogenization theory and the time–temperature superposition principle (TTSP). It is carried out in two steps. Firstly, the effective viscoelasticity properties of yarn are calculated using microscale homogenization; yarn comprises elastic fibers, interface, and a viscoelastic matrix. Subsequently, the effective viscoelasticity properties of woven composites are computed by mesoscale homogenization; it consists of homogenized viscoelastic yarns and matrix. Moreover, the multiscale method is verified using the Mechanics of Structure genome (MSG) consequence. Finally, the effect of temperature, fiber twist angle, fiber array, and coating on either the yarn’s effective relaxation stiffness or the relaxation moduli of the woven composite is investigated. The results show that increased temperature shortens the relaxation time of viscoelastic woven composites, and fiber twist angle affects tensors in the relaxation stiffness matrix of the yarn; the coating affects the overall mechanical properties of woven composites as well.

1. Introduction

Woven composites have the advantage of high strength, high stiffness, good oxidation resistance, and excellent thermal stability. They have been widely utilized in aerospace, national defense, biomedical, and other industrial fields. Woven composites are often subjected to long-term load in engineering practice, which seriously affects their viscoelastic behavior leading to the failure of the structure, especially at high temperature. Stress relaxation occurs when viscoelastic materials are subjected to static or variable strain continuously and this phenomenon is called creep [1]. Researchers have now developed a number of approaches to investigate the viscoelasticity of composites. Hashin [1,2] proposes the elastic viscoelastic correspondence principle. Via this principle, the parametric three-dimensional finite-volume direct averaging micromechanics (FVDAM) by Chen et al. [3] and the elastic-based locally exact homogenization theory (LEHT) by Wang and Pindera [4] are developed to accommodate linearly viscoelastic phase response. As an extension of the Eshelby-based Mori–Tanaka (MT) model, Weng et al. [5] investigated the overall viscoelastic behavior of composites with different shapes according to different aspect ratios. Katouzian et al. [6] employ the MT method to study the response of viscoelastic composites vs. the time compared with the experiment data. Moreover, Yang et al. [7,8,9] experimentally studied the long-term creep behavior of fiber-reinforced composite tubes subjected to flexural loading. Martynenko et al. [10,11] conduct numerical simulations and experiments to investigate the effect of the temperature and the lasting time on the effective viscoelasticity of fiber-reinforced composites. Kwok et al. [12] establish a viscoelastic model of single-layer plain woven carbon fiber-reinforced epoxy composites after being folded for a period of time compared with experimental measurements to study the deployment of tape springs at different times and temperatures. However, the classical models, elasticity-based homogenization approaches, or mechanics experiments among techniques are presently limited. For instance, the MT model may provide a reasonable estimate of elastic homogenized moduli but not the local stress distributions [3]; the LEHT is limited to composites reinforced by long cylindrical fibers, and experiments cost much money, time, and effort.
In recent years, homogenization theory (HT), the finite element method (FEM), and multiscale analysis have been widely used to predict woven composites’ stiffness and strength. For example, Zhao et al. [13] used multiphysics locally exact homogenization theory to investigate the multiscale homogenized thermal conductivity and thermomechanical properties of Advanced European superconductors filament groups (EAS). Pathan et al. [14,15] studied the sensitivity of the viscoelastic response of fiber-reinforced polymers (FRPs) to fiber shapes, fiber volume fraction, interphase volume fraction, and interphase properties via FEM. Also, they compare the Monte Carlo simulations’ results on random RVEs to those obtained by analyzing periodic square and hexagonal unit cells. Deviredy et al. [16] investigate the square array RVE (S-RVE) and hexagonal array RVE (H-RVE) of FPRs with the circular and square cross section of fiber and calculate the elastic modulus and thermal conductivity. Moreover, Liu et al. [17] used mechanics of structure genome (MSG) solid and plate model to capture the long-term viscoelastic behaviors of textile composites. Rique et al. [18] extend MSG to construct a linear thermo-viscoelastic model to analyze three-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. Seifert et al. [19] proposed a finite element-based micromechanical model to obtain the viscoelastic properties of glass fiber composites at high temperatures. Cai et al. [20] propose a multiscale model of three-dimensional four-way woven composites to analyze the influence of braiding angles and fiber volume fractions on the viscoelastic properties by FEM and creep experiment.
Further, the microscopic parameters of woven composites relating to viscoelasticity should be studied more comprehensively. As shown in Figure 1a,b, a traditional twist yarn in a fabric consists of fibers at an angle with the direction of the yarn axis. Twisted fiber bundles characterized by the twist angel, distinguishing significantly from untwisted ones in the mechanical property, contribute to the mechanical properties of the composite through fiber–yarn–fabric sequence. In Figure 1d, a region exists at the boundary between fibers and matrix, which possesses mechanical properties different from those of the fibers and the matrix due to physical and chemical reactions between the two main phases; this region is usually modeled as a physical coating of the continuum [21], so that the thickness, the material property, and the bonding condition of the coating significantly affect the overall mechanical properties of the woven composite. Some relevant studies are as follows. Miao et al. [22,23] attempt to optimize the yarn structure of fiber-reinforced polymer composites and establish the relationship between the fiber twist angle and mechanical properties of unidirectional fiber composites. Xiong et al. [24] construct a multiscale mechanical model and study the influence of fiber twist angle on the mechanical properties of plain woven composites. Fisher and Brinson [25] analyze the mechanical response of fiber-reinforced polymer matrix composites with viscoelastic interface regions by Mori–Tanaka micromechanical model and study the physical aging of viscoelastic composites. To obtain the viscoelastic behavior of polymer-based heterogeneous materials, Huang et al. [26] developed 3D viscoelastic calculated grains (CGs) containing spherical inclusions, interfacial phases/coatings, and non-interfacial phases/coatings. Yang et al. [27] theoretically study the frequency- and temperature-dependent viscoelastic behavior of the short fiber-reinforced polymers (SFRPs) and consider interface/interface conditions. Although these applications provide various aspects of viscoelastic behavior of composites, they are limited to composites with infinitesimally small microstructures.
Usually, multiscale methods are classified as hierarchical and concurrent methods by Belystchko and Song [31]. Based on the structure of woven composite having hierarchy and feature, this paper develops a hierarchical multiscale method to predict the viscoelastic of woven composite. The technique can study the influence of microscopic parameters on the long-term viscoelastic behavior of two-dimensional woven composites. Considering the non-isotropic viscoelastic materials cannot be calculated directly by commercial finite element software, the mesoscale estimate in the multiscale method introduces the discretization theory. Therefore, the long-term viscoelastic properties of yarns are divided into multiple transient properties with equal intervals. The long-term anisotropic viscoelastic simulation of RVE2 is equivalent to the instantaneous elastic simulation. This method is easy to calculate the viscoelastic properties of 2D woven composites. We were, moreover, using the multiscale approach to investigate the effect of more factors, such as twist yarn, coating interface, fiber array, and ambient temperature.
The remainder part of the paper is organized as follows: Section 2 establishes a multiscale RVE structure of a 2D woven composite and defines the fiber twist angle. Section 3 describes a viscoelastic multiscale model considering fiber twist angle and coating interface. Section 4 verifies the accuracy of the multiscale method via MSG methods. In Section 5 we investigate the combined effects of temperature, fiber twist angle, coating interface, and array type on the homogenized viscoelasticity moduli of woven composites, reporting new results. Conclusions are presented in Section 6.

2. The Microstructure, Multiscale Framework and RVEs of Woven Composites

2.1. Multiscale Framework and RVEs

The 2D and 3D woven composites have complex hierarchical structures composed with matrix and interwoven yarns and the fill yarn and warp yarn have periodical characteristics, respectively. The mechanical properties of woven composites mainly depend on the intrinsic properties of the fiber and matrix, the yarn, and the woven composite structure. The yarn structure, such as its geometry, internal fiber orientation, and fiber volume fraction, determines the main properties of the woven composite, especially mechanical properties [32]. As shown in Figure 2, it can be seen that the hierarchical structures of 2D woven composites consist of three scales: microscale, mesoscale, and macroscale. We start from the modeling of the yarn as RVE1 at microscale, which consists of twisted or untwisted fibers distributed in the matrix and there is the interface between the fiber and matrix; in this paper, we account the fiber with hexagonal or square array in matrix. The mesoscale as RVE2 consists of warp and fill yarns embedding in the matrix. The macroscale is 2D woven composites which is composed with RVE2s. The numerical examples of the latter two scales used in this paper were generated and analyzed by Texgen software (v3.12.0, University of Nottingham, Nottingham, UK) [33] and finite element analysis software ABAQUS 2020. Based on the “bottom-up” analysis process, we predict the overall performance of macroscopic 2D woven composites from the inherent properties of microscopic fibers and matrices. Fibers are assumed to be linear elastic and transversely isotropic, while the matrix is viscoelastic and isotropic.

2.2. Twist Angle and Coating of Fiber in the Yarn

The viscoelastic properties of every scale model are subject to the components and structure. Yarns play an important role in the woven composites. Effective properties of yarn are computed based on fibers and matrix. The twisted yarn axis has an angle θ with the fiber axis, and the angle is called a torsion angle or twist angle as shown in Figure 1a,b. During the injection molding and extrusion process [34], the twist angle of a fiber, related to the fiber orientation, depends on its radical location in the yarn’s cross section and changes gradually and continuously from the interior of the yarn to its surface, as shown in Figure 3, so the fiber orientation distribution function can be defined as seen in the following. Xiong et al. [24] used a spatial coordinate system to define the orientation of fibers in the yarn in order to describe the fiber orientation distribution (FOD). As shown in Figure 3b [34], 1-axis is the principal axis of the fiber, and the 2-axis and the 3-axisareradial-axes perpendicular to the 1-axis. x 1 , y 1 and z 1 are the coordinate system of the RVE1. θ is defined as the angle between 1-axis and the z 1 -axis of the spatial coordinate system; the value range is ( θ 0 , θ 0 ) . θ 0 is the twisting of the fibers on the surface of the yarn. φ between the projection of the fiber on the x 1 y 1 plane and the y 1 -axis; the value range is ( 0 , 2 π ) .
The twist angle of the fiber in the yarn depends on the radial position of the yarn cross section, and the twist angle of the fiber near the surface of the yarn is greater than the twist angle inside the yarn, as shown in Figure 3b. The fiber orientation distribution function can be defined as [24]:
f ( θ , φ ) = g ( θ ) g ( φ ) ( θ 0 θ θ 0 , 0 φ 2 π )
g ( θ ) = | tan θ | sec 2 θ tan 2 θ 0 ( θ 0 θ θ 0 )
g ( φ ) = 1 2 π
θ 0 θ 0 0 2 π f ( θ , φ ) d θ d φ = 1
Although fibers in the yarn are randomly arranged, the microstructure of the yarn is usually idealized in the literature as a square or hexagonal pack, and representative volume element is S-RVE or H-RVE, as shown in Figure 4. To be considered are the interfaces between the fiber and matrix (shown in Figure 1d) and it is assumed the interfaces or coating are perfect enough to connect the fiber and matrix together. a is the diameter of the fiber and b is the exterior diameter of the coating. l 1 , l 2 and l 3 are the lengths of RVEs. In Figure 5, the mesoscale RVE2 model was created by Texgen software. w and h are the width and thickness of the yarn embedding in the matrix of L z thickness and L x length in the periodical interval c .

3. Multiscale Viscoelastic Model of Woven Composites

3.1. Viscoelastic Constitutive Model

The constitutive relationship of viscoelastic materials depends on the time and the temperature, and the uniaxial stress–strain relationship of isotropic viscoelastic materials can be represented by the Boltzmann’s heredity integral:
σ ( t ) = 0 t E ( t τ ) d ε ( τ ) d τ d τ
ε ( t ) = 0 t D ( t τ ) d σ ( τ ) d τ d τ
0 t E ( t τ ) d D ( t τ ) d τ d τ = 1
t is a time variable, τ is a passed time, σ ( t ) is the time-dependent stress, ε ( t ) is the time-dependent strain, D ( t ) is the creep modulus, and the elastic modulus E ( t ) can be defined by the Prony series over a wide range of timescale [35]:
E ( t ) = E + i = 1 n E i e t / ρ i
E is the long-term modulus, E i are the Prony coefficients, and ρ i are the relaxation times, whose value is the larger, the slower the stress relaxation decays.
Temperature has a significant influence on viscoelastic behavior, and TTSP (the time–temperature superposition principle) provides a theoretical foundation to study the long-term behavior of viscoelastic materials by short-term characterization experiments. This principle relates the temperature-dependent relaxation modulus to the time-dependent one, just as E ( t , T ) = E ( t , T 0 ) , which indicates that the modulus at temperature T and time t is equal to that at the reference temperature T 0 and shortened time t . Therefore, we can cut down the observation time to study the long-term viscoelastic behavior of the polymer by increasing the temperature of the experimental environment, and this equivalence can be achieved by means of a conversion factor a T , which is the ratio of the relaxation time at two different temperatures:
a T = ρ ( T ) ρ ( T 0 )
where T is the general temperature, and T 0 is the reference temperature that is often taken as the glass transition temperature. In order to illustrate the dependence of viscoelasticity on temperature, the viscoelastic material studied in this paper is a thermos-rheological simplicity material and the conversion factor a T of the material is applicable over the whole relaxation time. When the temperature changes with time, the shortened time is obtained by integration:
t ( t ) = 0 t d τ a T ( T ( τ ) )
According to TTSP, the temperature shift function can be expressed by the Williams–Landel–Ferry (WLF) equation [36]:
log a T = c 1 ( T T 0 ) c 2 + ( T T 0 )
here c 1 and c 2 are material-dependent constants and the logarithm base is 10. Thus, the viscoelastic constitutive relation Equation (5) at temperature T can be rewritten as:
σ ( t ) = 0 t E ( t τ , T ) d ε ( τ ) d τ d τ
For the anisotropic material, the stress–strain relationship can be expressed in the time domain as:
σ ( t ) = 0 t C ( t τ , T ) d ε ( τ ) d τ d τ
Then, the relaxation modulus tensor can be expressed as
C i j k l = C i j k l , + m = 1 n C i j k l , m e t / ρ m

3.2. Multiscale Homogenization

The homogenization theory (HT) is employed by [3,4,13]. The theory can harvest the equivalent material properties of the composite material having periodic feature. In the twisted yarn, the helical fibers in different radial positions have different twist angles (shown in Figure 3). Thus, the stiffness matrix of each fiber and coating at the twist angle can be calculated as [24]:
[ C ] = [ M ] [ C 0 ] [ M ] T
[ C 0 ] is the stiffness matrix of a single fiber and coating at θ = 0 , [ M ] is the stiffness transformation matrix between the offset coordinate system and the global coordinate system, and [ M ] T is the transpose matrix of [ M ] . The definitions of [ C 0 ] and [ M ] are described in Appendix A Equations (A1)–(A3).
Since the fiber is an elastic material, its stiffness matrix does not change with time. A series of helical fibers located at the point of the same radius in the yarn cross section can be assumed as a layer, whose orientation changes continuously from the center to the surface of the yarn. Based on the directional average method and the iso-strain hypothesis [34,37], the stiffness matrix of each helix layer is integrated by helix angle to obtain the average stiffness matrix after fiber and coating torsion [24]:
[ C f ] = 0 2 π θ 0 θ 0 [   C ] f ( θ , φ ) d θ d φ
The yarn is composed of fibers, coatings, and the matrix. The first two of these are elastic and the matrix is made from the polymer of viscoelastic elasticity, so the yarn is transversely isotropic. According to the model of the RVE1 established in the previous section, the overall average stress component of the twisted yarn can be expressed as:
σ ˜ i R V E 1 ( t ) = 1 V R V E 1 0 t V m C i j m ( t τ , T ) d ε j ( t ) d τ d τ d V m + 1 V R V E 1 V f σ i j f d V f + 1 V R V E 1 V i n t σ i j i n t d V i n t
C ˜ i y ( t ) = σ ˜ i R V E 1 ( t )   with   ε j ( t ) = 1
where f , i n t and m represent the fiber, interfacial phase and matrix, and the superscript “~” in the formulation indicates the homogenized sign of linear viscoelasticity. σ ˜ i R V E 1 ( t ) is the average time-varying stress field in RVE1, C ˜ i y ( t ) is the time-varying relaxation stiffness tensor in RVE1, and V R V E 1 is the volume.
The in-plane structure of the 2D woven composite material is symmetrical and periodical, and the mesoscale model is in-plane isotropic. The relaxation stiffness matrix of warp and fill yarns in the global coordinate system can be described as:
[ C y ( t ) ] = { [ T w a r p ] [ C ˜ y ( t ) ] [ T w a r p ] T ( w a r p ) [ T f i l l ] [ C ˜ y ( t ) ] [ T f i l l ] T ( f i l l )
Here [ C ˜ y ( t ) ] is the equivalent stiffness matrix of the yarn in RVE2. [ T w a r p ] and [ T f i l l ] are the transformation matrix between the microscale coordinate system of the warp yarn (fill yarn) and the mesoscale coordinate system; [ T w a r p ] T and [ T f i l l ] T are the transpose matrices of [ T w a r p ] and [ T f i l l ] , respectively, which are described in Equation (A4) of Appendix A.
Based on the RVE2 built in the previous section, the average stress component of the 2D woven composite can be expressed as:
σ ˜ i j R V E 2 , λ ( t ) = 1 V 0 t V C ˜ i j R V E 2 , λ ( t τ , T ) d ε j ( t ) d τ d τ d V , ( λ = w a r p , f i l l , m )
σ ˜ i j R V E 2 = v w a r p R V E 2 σ ˜ i j ( R V E 2 , w a r p ) + v f i l l R V E 2 σ ˜ i j ( R V E 2 , f i l l ) + v m R V E 2 σ ˜ i j ( R V E 2 , m )
In the formulation, the subscript λ means that the variable belongs to the corresponding phase, such as the warp, the fill, and the matrix.
We can apply the corresponding normal strain or shear strain to the single-cell model of different scales, and then solve the viscoelastic partial differential equation according to the above analytical models to restore the stress distribution in the numerical microelement, and finally calculate the macroscopic overall average stress component of the 2D woven composite material in the form of integration:
σ ˜ i j ( t ) = 1 V 0 t V C ˜ i j R V E 2 ( t τ , T ) d ε j ( t ) d τ d τ d V
The constitutive equations of 2D woven composites are represented in the form of relaxation stiffness matrix engineering:
{ σ ˜ x x σ ˜ y y σ ˜ z z σ ˜ y z σ ˜ x z σ ˜ x y } = [   C ˜ x x   C ˜ x y   C ˜ x z 0 0 0   C ˜ x y   C ˜ y y   C ˜ x z 0 0 0   C ˜ x z   C ˜ x z   C ˜ z z 0 0 0 0 0 0   C ˜ y z 0 0 0 0 0 0   C ˜ x z 0 0 0 0 0 0   C ˜ x y ] { ε x x ε y y ε z z ε y z ε x z ε x y }

3.3. Periodic Boundary Conditions of RVEs

Periodic boundary conditions must be applied to RVE of different scales to ensure uniform deformation and continuity of RVE under various loading conditions. In this paper, the periodic boundary condition proposed by [38] is adopted, and the displacement equation for each pair of nodes on the parallel boundary surface of RVE is as follows:
u i j + = ε ¯ i k x k j + + u i *
u i j = ε ¯ i k x k j + u i *
where the superscripts j + and j refer to the two corresponding nodes on the boundary of RVE and ε ¯ i k is the average strain of RVE. x k is the coordinate value of k -axis and u i * is the periodic displacement of i -axis. The relative displacement between two corresponding nodes on a parallel boundary surface can be expressed as:
u i j + u i j = ε ¯ i k ( x k j + x k j ) = ε ¯ i k Δ x k j . ( i , k = 1 , 2 , 3 )
Equation (26) can be easily implemented by setting the displacement constraint on each pair of coupled nodes on the boundary of RVE.

4. Validation

To verify the result based on the multiscale model and homogenization viscoelasticity to every scale RVE, the results at the microscale predicted in this paper are compared with that of the literature [17]. In this case, the yarn has no coating and twist angel, and the fiber volume fraction v f is fixed at 0.64 in the yarn. The elastic properties of the fiber are shown in Table 1, the viscoelasticity of the matrix is expressed by the relaxation time and Prony coefficient, as shown in Table 2, and the Poisson’s ratio is assumed to be constant at 0.33.
FE models of S-RVE1 and H-RVE1 are established in ABAQUS 2020, as shown in Figure 4a. Each analysis consists of two steps and the relaxation stiffness is defined with the option ∗VISCO in ABAQUS. The first step lasting a short time period (0.1 s) is to apply a unit strain on RVE1. In the second step, the strain remains constant and lasts 10 10 s.
In Figure 6, the results of microscale S-RVE1 and H-RVE1 in this paper are compared with that of Liu-MSG in [17], whose work only focuses on S-RVE1 and both results agree well. Additionally, results of two RVE1s have a small difference in axial relaxation stiffness tensor C ˜ 11 R V E 1 , the maximum difference of which is between 0.05~2.2%, while the maximum difference in radial relaxation stiffness tensor C ˜ 22 R V E 1 is about 2.1~8.4%. The effect of time for tensor C ˜ 11 R V E 1 is almost negligible, because the behavior in this direction is dominated by the fibers, whose elastic behavior is independent of time. Other tensors of the relaxation stiffness decrease with the time increased, exhibiting a trend of variation over time. More importantly, the microstructure has more significant effect on the radial relaxation tensor than that on the axial relaxation one, which can be explained by the dominated contribution of fibers to the axial direction of the yarn and the behavioral gap between the matrix and the fiber which is not obvious in the radial direction.
We employed the FEM to get the homogenization viscoelasticity prediction for the woven composite. The width and thickness L z of the warp or fill in Figure 5 are 0.9 mm and 0.06 mm, respectively, while the thickness of the fabric is 0.12 mm, and the interval of neighboring warps or fills is 1.75 mm. Progressive meshes with C3D8R and C3D10 solid elements are adopted in ABAQUS to guarantee the convergence, and co-node grids are set up at the interface between the yarn and the matrix. It is well established that yarn and RVE2 are not isotropic, and finite element software cannot directly define the viscoelasticity of the anisotropic material [17]. However, the long-term relaxation viscoelasticity of RVE2 has consisted of many transient responses. Based on the discretization theory, the long-term relaxation response of RVE2 can be decomposed into multiple transient reactions with the same interval, simplifying the finite element calculation. The yarn’s transient properties are summarized in Table A1 and Table A2 of Appendix A. Then, periodic boundary conditions can be applied to RVE2, and the woven composite’s overall mechanical properties, E i j R V E 2 ( G i j R V E 2 ) , are obtained via Equations (20)–(23) and (A5) of the Appendix A. In Figure 7, the results of the paper are compared with that of Rique–MSG in [18], and both of results are in good agreement.
Figure 7 shows the relaxation moduli of plain woven composites by viscoelastic multiscale homogenization. S-RVE2 and H-RVE2 represent the relaxation moduli of the mesoscale corresponding to S-RVE1 or H-RVE1, respectively. Rique–MSG represents the relaxation moduli of plain woven composites calculated by Rique et al. [18] using the Mechanics of Structure genome MSG under periodic boundary conditions. In Rique’s research, the arrangement of fibers in the yarn is a square array. Therefore, in Figure 7, the results of S-RVE2 are closer to those of Rique–MSG, and the consequences of S-RVE2 and H-RVE2 have noticeable differences only in the in-plane tensile moduli. In contrast, the differences in other moduli are slight.

5. Numerical Investigations

Based on the multiscale homogenization, this section investigates the effect of the temperature and several other microstructural parameters, such as the twist angle of fibers, interfacial parameters, and fiber arrangements, on the relaxation modulus of two-dimensional woven composites.

5.1. The Effect of Temperature

The relaxation modulus curve of PMT-F4 is obtained by creep experiment [20]. The material-dependent constants c 1 , c 2 are found by fitting data of the conversion factor vs. temperature data to the WLF equation (Equation (11)) and then we can infer that c 1 = 28.3816 , c 2 = 93.291 at the reference temperature T 0 = 40 °C. This section assumes that the elastic properties of fibers are independent of the temperature and the fiber volume fraction v f and the twist angle θ in each yarn are chosen to be 0.64 and 0 ° , taking no account of the coating’s influence, in order to investigate the effect of the temperature on the viscoelastic properties of 2D woven composites. In Abaqus 6.14, the relationship between the viscoelasticity of the matrix and the temperature is defined in the form of the WLF equation and the temperature field is applied to the matrix.
Figure 8 and Figure 9 illustrate the relaxation response of yarns and 2D woven composites to different temperatures. With the temperature increasing, the relaxation time of yarn and woven composites gradually shortens, so the effect of temperature on the viscoelasticity of the yarn and the woven composite conforms to TTSP. Moreover, the RVE1’s (yarn) relaxation tensor C ˜ i j R V E 1 and the RVE2′s (woven composite) relaxation modulus E i j R V E 2 ( G i j R V E 2 ) decrease sharply in the time domain with the temperature increasing, and the higher the temperature leads to the faster relaxation. The phenomenon is due to the molecular movement stimulated by the temperature. When the temperature rises, the viscous frictional resistance between molecules becomes smaller and the relaxation reaction of molecules is easier. In addition, there is no significant distinction between H-RVE1 and S-RVE1 in the trend of stiffness vs. time curve and the relaxation time. Figure 8a–c,f show that tensors C ˜ 11 R V E 1 , C ˜ 12 R V E 1 , C ˜ 22 R V E 1 , and C ˜ 55 R V E 1 in the H-RVE1 are significantly lower than that of S-RVE1, while the tensors C ˜ 23 R V E 1 and C ˜ 44 R V E 1 of H-RVE1 are significantly higher than that of S-RVE1 in Figure 8d,e. At the mesoscale, as illustrated in Figure 9, the impact of RVE1’s structure on RVE2’s homogenized viscoelastic moduli E i j R V E 2 ( G i j R V E 2 ) is insignificant.

5.2. The Effect of the Fiber Twist Angle

After studying the effect of the temperature on the viscoelasticity of RVEs of different scales, the influence of the twist angle is investigated considering θ = 0 , 30 , and 60 at the temperature fixed at 60 °C in this section. The effective stiffness matrix of the twist fiber calculated by Equation (15) in Section 3.2 is shown in Table 3. In spite of the twist angle, the transverse elastic moduli E f 2 equals to E f 3 as well as the shear moduli G f 12 equals to G f 13 , so the twist fiber can be considered transverse isotropic. Table 3 shows that the twist angle of the fiber is increasing, and the effective elastic properties of the fiber change dramatically where these radial properties are effectively improved.
Figure 10 shows the variation of the relaxation stiffness of the yarn affected by a wide range of the twist angle. The relaxation stiffness tensor   C ˜ 11 R V E 1 decreases sharply with the increase of the twist angle, but the other relaxation stiffness components increase. In Figure 10a, the effective normal stiffness tensor   C ˜ 11 R V E 1 of the yarn with the twist angle of 30° is only about 84% of the twist less yarn, while the ratio decrease to 30% when the twist angle is 60°. In Figure 10b,f, the twist angle promotes a large improvement on   C ˜ 12 R V E 1 and   C ˜ 55 R V E 1 . Surprisingly, the value of   C ˜ 12 R V E 1 on the twist angle 60° gradually decreases with the time increasing and finally is lower than that on the twist angle 30°. The effect of the high twist angle is more prominent on   C ˜ 22 R V E 1 ,   C ˜ 23 R V E 1 , and   C ˜ 44 R V E 1 , but this effect will diminish after a period of time. Additionally, curve trends of the two microstructures H-RVE1 and S-RVE1 are almost the same, but their values are different.
Figure 11 shows the effect of the fiber twist angle on the relaxation moduli of a two-dimensional woven composite. It can be seen that with the increase of the fiber twist angle, the in-plane relaxation moduli E x x R V E 2 and E y y R V E 2 of the woven composite decrease sharply, but the variation of the out-plane relaxation modulus E z z R V E 2 and the shear moduli G x y R V E 2 , G x z R V E 2 , and G y z R V E 2 induced by the fiber twist angle is not obvious. Therefore, the fiber twist angle can be realized as a non-main parameter in the optimal design of 2D woven composites. In addition, there is no significant difference in E x x R V E 2 , E y y R V E 2 , G x z R V E 2 , and G y z R V E 2 whenever fibers are distributed in square array or hexagonal array, but the value of E z z R V E 2 and G x y R V E 2 of S-RVE2 are higher than that of H-RVE2.

5.3. The Effect of the Coating Thickness

The physical coating between the fiber and matrix is accepted as a kind of component existing in the yarn [21]. To demonstrate the effect of the coating thickness which is represented by the ratio b/a, b is the coating diameter and a is fiber diameter, respectively, shown in Figure 4, b/a = 1 indicates that there is no coating between the fiber and the matrix, while b/a > 1 means the coating is surrounding the fiber. The field temperature is fixed at 60 °C and the volume fraction v f of the untwist fiber in each yarn is 0.64. The mechanical properties of Pyrolytic carbon used for the coating are given in Table 4.
The coating thickness significantly improves the viscoelastic properties of the yarn and affects the axial and radial relaxation components in varying degrees, just as shown in Figure 12. The normal stiffness tensor   C ˜ 11 R V E 1 of b/a = 1.05 is greatly improved compared with that of b/a = 1,and the microstructure of the hexagonal array has greater impact on   C ˜ 11 R V E 1 than that of the square array. However, the effect of the coating on   C ˜ 12 R V E 1 is very insignificant. At the initial moment t = 1 s,   C ˜ 12 R V E 1 at the ratio b/a = 1.05 is lower than that at b/a = 1; with the increase of the time,   C ˜ 12 R V E 1 gradually decreases and becomes higher than the value of b/a = 1. The coating thickness also improves the radial relaxation stiffness tensor   C ˜ 22 R V E 1 and   C ˜ 23 R V E 1 to an extent as illustrated in Figure 12c,d. In Figure 12e,f, the coating has a greater effect on the axial shear tensor   C ˜ 55 R V E 1 than the radial shear tensor   C ˜ 44 R V E 1 , which depends on the property of the coating, whose axial mechanical stiffness is much higher than the radial one.
In Figure 13, the coating thickness significantly improves the in-plane relaxation moduli E x x R V E 2 and E y y R V E 2 of the 2D woven composite, but the out-plane relaxation modulus E z z R V E 2 and the shear moduli G x y R V E 2 , G y z R V E 2 are slightly promoted. In general, excellent coating properties are necessary to improve the overall mechanical properties of the 2D woven composite. In Figure 13a–c, S-RVE2’s curves are slightly higher than H-RVE2, while both curves ( G x z R V E 2 and G y z R V E 2 ) of the two arrays have no apparent distinction in Figure 13d.

6. Conclusions

This paper investigates the viscoelastic properties of 2D woven composite employing elastic-based homogenization theory to establish a multiscale model considering the coating and the fiber twist angle. According to TTSP, the influence of temperature on the relaxation time of the 2D woven composites is studied. The effects of the fiber twist angle and the coating on the viscoelasticity are discussed in detail. The conclusions are summarized as below:
(1)
The results using the multiscale method show that the fiber array considerably affects stiffness relaxation of the yarn. C ˜ 11 R V E 1 , C ˜ 12 R V E 1 , and C ˜ 22 R V E 1 of the square array are higher than H-RVE1, while the H-RVE1’s tensors C ˜ 23 R V E 1 and C ˜ 44 R V E 1 are higher than that of S-RVE1. At the same temperature, the relaxation time and variation trend of S-RVE1 and H-RVE1 are almost identical.
(2)
The multiscale solutions show that the yarn surface twist angle has significant affection to the viscoelastic properties of the composites. The negative effect of high twist angle on moduli E x x R V E 2 and E y y R V E 2 are more important, while the improvement on other modulus is minor, and gradually disappears with the time. In addition, the lower twist angle has a more significant effect on the axial stiffness of the yarn, while the radial stiffness is more sensitive to the higher angle.
(3)
The coating, the material property, and the thickness can effectively improve the overall viscoelasticity of 2D woven composites, especially the in-plane relaxation moduli. When the stiffness of the coating is higher than that of the matrix, the coating will effectively improve the overall mechanical properties of the composite. Designing the coating is significant in exploiting the potentiality of 2D woven composites.
(4)
The multiscale method facilitates the calculation of the viscoelasticity of woven composites. Combining with the discrete theory and FEM, this paper provides an appropriate approach for analyzing non-isotropic composites. In addition, the effect of more microscopic parameters on the mesoscopic properties is considered and calculated.

Author Contributions

Conceptualization, B.L., C.L. and X.Z.; methodology, B.L and X.Z.; software, B.L.; validation, B.L., X.Z. and F.G.; formal analysis, B.L. and F.G.; investigation, B.L.; resources, B.L.; writing—original draft preparation, B.L.; writing—review and editing, B.L., C.L., X.Z. and J.Y.; visualization, B.L., J.Y. and F.G.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Research and Development Program of China (No. 2020YFA0711700) and National Natural Science Foundation of China (No. 12002303) and ZJU-ZCCC Institute of Collaborative Innovation (No. ZDJG2021002). The authors have benefited from the communication with Guannan Wang of Zhejiang University.

Data Availability Statement

My research data are presented in the table in Appendix A.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Each fiber is assumed to be transverse isotropic and their stiffness matrix is given as follows [39]:
[ C 0 ] = [ C 11 0 C 12 0 C 12 0 0 0 0 C 12 0 C 22 0 C 23 0 0 0 0 C 12 0 C 23 0 C 22 0 0 0 0 0 0 0 C 44 0 0 0 0 0 0 0 C 55 0 0 0 0 0 0 0 C 55 0 ]
{ C 11 0 = 1 υ 23 υ 32 E 2 E 3 A C 12 0 = C 13 0 = C 31 0 = υ 12 + υ 32 υ 13 E 1 E 2 A C 22 0 = C 33 0 = 1 υ 13 υ 31 E 1 E 3 A C 23 0 = C 32 0 = υ 23 + υ 21 υ 13 E 1 E 2 A C 44 0 = C 22 0 C 23 0 2 = G 23 C 55 0 = C 66 0 = G 12 A = 1 υ 12 υ 21 υ 23 υ 32 υ 13 υ 31 2 υ 12 υ 23 υ 31 E 1 E 2 E 3
A coordinate transformation matrix is defined as [40]:
[ M ] = [ l 1 2 l 2 2 l 3 2 2 l 2 l 3 2 l 3 l 1 2 l 1 l 2 m 1 2 m 2 2 m 3 2 2 m 2 m 3 2 m 3 m 1 2 m 1 m 2 n 1 2 n 2 2 n 3 2 2 n 2 n 3 2 n 3 n 1 2 n 1 n 2 m 1 n 1 m 2 n 2 m 3 n 3 m 2 n 3 + m 3 n 2 m 1 n 3 + m 3 n 1 m 1 n 2 + m 2 n 1 n 1 l 1 n 2 l 2 n 3 l 3 n 2 l 3 + n 3 l 2 n 1 l 3 + n 3 l 1 n 1 l 2 + n 2 l 1 l 1 m 1 l 2 m 2 l 3 m 3 l 2 m 3 + l 3 m 2 l 1 m 3 + l 3 m 1 l 1 m 2 + l 2 m 1 ]
where l 1 = cos θ , l 2 = 0 , l 3 = sin θ , m 1 = sin θ sin φ , m 2 = cos φ , m 1 = cos θ cos φ ; n 1 = sin θ cos φ , n 2 = sin φ , n 3 = cos θ cos φ .
The coordinate transformation matrix [ T w a r p ] of the warp yarn and [ T f i l l ] of the fill yarn are defined as:
[ T w a r p ] = [ m 2 0 n 2 0 2 m n 0 0 1 0 0 0 0 n 2 0 m 2 0 2 m n 0 0 0 0 m 0 n m n 0 m n 0 m 2 n 2 0 0 0 0 n 0 m ] [ T f i l l ] = [ 0 1 0 0 0 0 u 2 0 v 2 0 2 u v 0 v 2 0 u 2 0 2 u v 0 u v 0 u v 0 u 2 v 2 0 0 0 0 u 0 v 0 0 0 v 0 u ]
where m = cos α 0 , n = sin α 0 , u = cos α 0 , v = sin α 0 and α 0 is undulation angle.
The relaxation stiffness [ C ˜ R V E 2 ] of RVE2 to calculate relaxation moduli is as follows:
[   C ˜ x x   C ˜ x y   C ˜ x z 0 0 0   C ˜ x y   C ˜ y y   C ˜ x z 0 0 0   C ˜ x z   C ˜ x z   C ˜ z z 0 0 0 0 0 0   C ˜ y z 0 0 0 0 0 0   C ˜ x z 0 0 0 0 0 0   C ˜ x y ] 1 = [ 1 E x x υ y x E y y υ z x E z z 0 0 0 υ x y E x x 1 E y y υ z y E z z 0 0 0 υ x z E x x υ y z E y y 1 E z z 0 0 0 0 0 0 1 G y z 0 0 0 0 0 0 1 G z x 0 0 0 0 0 0 1 G x y ]
Table A1. The effective viscoelasticity of S-RVE1 at microscale.
Table A1. The effective viscoelasticity of S-RVE1 at microscale.
s 1234567
λ s ( s ) -102103104105106107108
C ˜ 11 R V E 1 (MPa)150,185152,861152,790.3152,719.9152,558.5152,310.8152,068.5151,699.5
C ˜ 12 R V E 1 (MPa)14774049.664000.23951.223835.213648.913457.233149.23
C ˜ 22 R V E 1 (MPa)4839.212,127.811,991.11,849.8111,525.2910,999.110,451.49559.33
C ˜ 23 R V E 1 (MPa)1392.645824519.54454.744305.734063.473811.253401.7
C ˜ 44 R V E 1 (MPa)10464634.24571.994507.854360.734123.693879.073485
C ˜ 55 R V E 1 (MPa)15813335.453285.563234.273116.832928.932736.782430.8
Table A2. The effective viscoelasticity of H-RVE1 at microscale.
Table A2. The effective viscoelasticity of H-RVE1 at microscale.
s 1234567
λ s ( s ) -102103104105106107108
C ˜ 11 R V E 1 (MPa) 152,780.2152,708.9152,636.8152,473.6152,222.5151,975.7151,598
C ˜ 12 R V E 1 (MPa) 4035.733984.313931.253809.393612.23407.643076.14
C ˜ 22 R V E 1 (MPa) 11,860.111,713.2211,561.4511,212.6510,646.6610,057.749100.42
C ˜ 23 R V E 1 (MPa) 4772.604714.714654.694516.454290.54053.693666.18
C ˜ 44 R V E 1 (MPa) 4531.394466.614399.914247.024001.313748.723343.95
C ˜ 55 R V E 1 (MPa) 3615.553573.313529.743429.713267.893099.822826.58

References

  1. Hashin, Z. Viscoelastic fiber reinforced materials. AIAA J. 1966, 4, 1411–1417. [Google Scholar]
  2. Hashin, Z. Viscoelastic behavior of heterogeneous media. J. Appl. Mech. 1964, 32, 630–636. [Google Scholar]
  3. Chen, Q.; Wang, G.N.; Chen, X.F.; Geng, J. Finite-volume homogenization of elastic/viscoelastic periodic materials. Compos. Struct. 2017, 182, 457–470. [Google Scholar]
  4. Wang, G.N.; Pindera, M.J. Locally-exact homogenization of viscoelastic unidirectional composites. Mech. Mater. 2016, 103, 95–109. [Google Scholar]
  5. Weng, Y.M.; Wang, G.J. The Influence of Inclusion Shape on the Overall Viscoelastic Behavior of Compoisites. J. Appl. Mech. 1992, 59, 510–518. [Google Scholar] [CrossRef]
  6. Katouzian, M.; Vlase, S. Mori–Tanaka Formalism-Based Method Used to Estimate the Viscoelastic Parameters of Laminated Composites. Polymers 2020, 12, 2481. [Google Scholar]
  7. Mayookh Lal, H.; Xian, G.J.; Thomas, S.; Zhang, L.; Zhang, Z.H.; Wang, H.L. Experimental Study on the Flexural Creep Behaviors of Pultruded Unidirectional Carbon/Glass Fiber-Reinforced Hybrid Bars. Materials 2020, 13, 976. [Google Scholar]
  8. Doan, H.G.; Mertiny, P. Creep Testing of Thermoplastic Fiber-Reinforced Polymer Composite Tubular Coupons. Materials 2020, 13, 4637. [Google Scholar] [CrossRef]
  9. Yang, Z.; Wang, H.; Ma, X.; Shang, F.; Ma, Y.; Shao, Z.; Hou, D. Flexural creep tests and long-term mechanical behavior of fi ber-reinforced polymeric composite tubes. Compos. Struct. 2018, 193, 154–164. [Google Scholar] [CrossRef]
  10. Martynenko, V.G.; Lvov, G.I. Numerical prediction of temperature-dependent anisotropic viscoelastic properties of fiber reinforced composite. J. Reinf. Plast. Compos. 2017, 36, 1790–1801. [Google Scholar]
  11. Martynenko, V.G.; Lvov, G.I.; Ulianov, Y.N. Experimental investigation of anisotropic viscoelastic properties of glass fi ber-reinforced polymeric composite material. Polym. Polym. Compos. 2019, 27, 323–336. [Google Scholar]
  12. Kwok, K.; Pellegrino, S. Micromechanics models for viscoelastic plain-weave composite tape springs. AIAA J. 2016, 55, 309–321. [Google Scholar]
  13. Zhao, X.; Wang, G.; Chen, Q.; Duan, L.; Tu, W. An effective thermal conductivity and thermomechanical homogenization scheme for a multiscale Nb3Sn filaments. Nanotechnol. Rev. 2021, 10, 187–200. [Google Scholar]
  14. Pathan, M.V.; Tagarielli, V.L.; Patsias, S. Effect of fibre shape and interphase on the anisotropic viscoelastic response of fiber composites. Compos. Struct. 2017, 162, 156–163. [Google Scholar] [CrossRef]
  15. Pathan, M.V.; Tagarielli, V.L.; Patsias, S. Numerical predictions of the anisotropic viscoelastic response of uni-directional fiber composites. Compos. Part A Appl. Sci. Manuf. 2017, 93, 18–32. [Google Scholar]
  16. Devireddy, S.B.R.; Biswas, S. Effect of Fiber Geometry and Representative Volume Element on Elastic and Thermal Properties of Unidirectional Fiber-Reinforced Composites. J. Compos. 2014, 2014, 629175. [Google Scholar]
  17. Liu, X.; Tang, T.; Yu, W.; Pipes, R.B. Multiscale modeling of viscoelastic behaviors of textile composites. Int. J. Eng. Sci. 2018, 130, 175–186. [Google Scholar]
  18. Rique, O.; Liu, X.; Yu, W.; Pipes, R. Constitutive modeling for time- and temperature-dependent behavior of composites. Compos. Part B 2020, 184, 107726. [Google Scholar]
  19. Seifert Ole, E.; Schumacher, S.; Hansen, A. Viscoelastic properties of a glass fabric composite at elevated temperatures: Experimental and numerical results. Compos. Part B-Eng. 2003, 34, 571–586. [Google Scholar]
  20. Cai, Y.; Sun, H. Prediction on viscoelastic properties of three-dimensionally braided composites by multi-scale model. J. Mater. Sci. 2013, 48, 6499–6508. [Google Scholar] [CrossRef]
  21. Wang, G.N.; Pindera, M.J. Locally-exact homogenization of unidirectional composites with coated or hollow reinforcement. Mater. Des. 2016, 93, 514–528. [Google Scholar] [CrossRef]
  22. Gu, H.B.; Miao, M.H. Optimising fiber alignment in twisted yarns for natural fiber composites. J. Compos. Mater. 2013, 48, 2993–3002. [Google Scholar] [CrossRef]
  23. Zaidi, B.M.; Zhang, J.; Magniez, K.; Gu, H.; Miao, M.H. Optimizing twisted yarn structure for natural fiber-reinforced polymeric composites. J. Compos. Mater. 2018, 52, 373–381. [Google Scholar] [CrossRef]
  24. Xiong, X.S.; Hua, L.; Miao, M.H.; Shen, S.Z.; Li, X.; Wan, X.J.; Guo, W. Multi-scale constitutive modeling of natural fiber fabric reinforced composites. Compos. Part A Appl. Sci. Manuf. 2018, 115, 383–396. [Google Scholar] [CrossRef]
  25. Fisher, F.T.; Brinson, L.C. Viscoelastic interphases in polymer–matrix composites: Theoretical models and finite-element analysis. Compos. Sci. Technol. 2001, 61, 731–748. [Google Scholar] [CrossRef]
  26. Huang, Y.; Wang, G.; Dong, L.; Atluri, S.N. 3D Viscoelastic Computational Grains with Spherical Inclusions with or without Interphases/Coatings for Micromechanical Modeling of Heterogeneous Materials. Int. J. Numer. Methods Eng. 2021, 122, 4966–4986. [Google Scholar] [CrossRef]
  27. Yang, Y.; He, Q.; Rao, Y.N.; Dai, H.L. Estimation of dynamic thermos-viscoelastic moduli of short fiber-reinforced polymers based on a micromechanical model considering interphases/interfaces conditions. Polym. Compos. 2020, 41, 788–803. [Google Scholar] [CrossRef]
  28. Ma, H.; Li, Y.; Shen, Y.; Xie, L.; Wang, D. Effect of linear density and yarn structure on the mechanical properties of ramie fiber yarn reinforced composites. Compos. Part A Appl. Sci. Manuf. 2016, 87, 98–108. [Google Scholar] [CrossRef]
  29. Wang, C.G.; Tang, M.; Liu, W.K.; Zhu, T. Study on Microstructure Characteristics of Axially Braided Carbon/Carbon Composites Based on SEM and Micro-CT. Materials 2020, 13, 1414. [Google Scholar] [CrossRef] [Green Version]
  30. Xu, Y.; Zhang, W.; Bassir, D. Stress analysis of multi-phase and multi-layer plain weave composite structure using glob al/local approach. Compos. Struct. 2010, 92, 1143–1154. [Google Scholar] [CrossRef]
  31. Belytschko, T.; Song, J.H. Coarse-graining of multiscale crack propagation. Int. J. Numer. Methods Eng. 2010, 81, 537–563. [Google Scholar] [CrossRef]
  32. Matthews, F.L.; Rawlings, R.D. Composite Materials: Engineering and Science; Chapman & Hall: London, UK, 1994. [Google Scholar]
  33. Lin, H.; Brown, L.P.; Long, A.C. Modelling and simulating textile structures using TexGen. Adv. Mater. Res. 2011, 331, 44–47. [Google Scholar]
  34. Fu, S.Y.; Lauke, B. The elastic modulus of misaligned short-fiber-reinforced polymers. Compos. Sci. Technol. 1998, 58, 389–400. [Google Scholar] [CrossRef]
  35. Barbero, E.J. Finite Element Analysis of Composite Materials Using Abaqus; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  36. Williams, M.L.; Landel, R.F.; Ferry, J.D. Temperature Dependence of Relaxation Mechanisms the Temperature Dependence of Relaxation Mechanisms in Amorphous Polymers and Other Glass-Forming Liquids. J. Am. Chem. Soc. 1955, 77, 3701–3707. [Google Scholar] [CrossRef]
  37. Modniks, J.; Andersons, J. Modeling the non-linear deformation of a short-flax-fiber-reinforced polymer composite by orientation averaging. Compos. B Eng. 2013, 54, 188–193. [Google Scholar] [CrossRef]
  38. Xia, Z.; Zhou, C.W.; Yong, Q.; Wang, X.W. On selection of repeated unit cell model and application of unified periodic boundary conditions in micro-mechanical analysis of composites. Int. J. Solids Struct. 2006, 43, 266–278. [Google Scholar] [CrossRef] [Green Version]
  39. Udhayaraman, R.; Mulay, S.S. Multi-scale approach based constitutive modelling of plain woven textile composites. Mech. Mater. 2017, 112, 172–192. [Google Scholar] [CrossRef]
  40. Huang, Z.M. The mechanical properties of composites reinforced with woven and braided fabrics. Compos. Sci. Technol. 2000, 60, 479–498. [Google Scholar] [CrossRef]
Figure 1. SEM photographs: (a) Natural fiber fabric [24]; (b) The orientation of sisal fibers in the yarn [28]; (c) The fiber rod end surface [29]; (d) The fiber and coating in the yarn [30]. (Reproduced with permission ref. [24]. Copyright 2018 Elsevier; ref. [28]. Copyright 2016 Elsevier; ref. [29]. Copyright 2020 MDPI; ref. [30]. Copyright 2010 Elsevier).
Figure 1. SEM photographs: (a) Natural fiber fabric [24]; (b) The orientation of sisal fibers in the yarn [28]; (c) The fiber rod end surface [29]; (d) The fiber and coating in the yarn [30]. (Reproduced with permission ref. [24]. Copyright 2018 Elsevier; ref. [28]. Copyright 2016 Elsevier; ref. [29]. Copyright 2020 MDPI; ref. [30]. Copyright 2010 Elsevier).
Materials 16 02689 g001
Figure 2. Multiscale framework of the 2D woven composites considering the twisted angle and interface in RVE. The H-RVE is fiber arrayed in hexagonal and the S-RVE is fiber arrayed in square.
Figure 2. Multiscale framework of the 2D woven composites considering the twisted angle and interface in RVE. The H-RVE is fiber arrayed in hexagonal and the S-RVE is fiber arrayed in square.
Materials 16 02689 g002
Figure 3. Definitions of the fiber orientation angles. (a) The twisting angle of fiber from inner periphery to outer periphery in the yarn; (b) The relation between the fiber orientation with the RVE1 coordinate system.
Figure 3. Definitions of the fiber orientation angles. (a) The twisting angle of fiber from inner periphery to outer periphery in the yarn; (b) The relation between the fiber orientation with the RVE1 coordinate system.
Materials 16 02689 g003
Figure 4. Microscale RVE of the yarn consisting of fibers and the matrix: (a) H-RVE; (b) S-RVE.
Figure 4. Microscale RVE of the yarn consisting of fibers and the matrix: (a) H-RVE; (b) S-RVE.
Materials 16 02689 g004
Figure 5. FE model of mesoscale RVE: (a) the woven yarn; (b) the matrix.
Figure 5. FE model of mesoscale RVE: (a) the woven yarn; (b) the matrix.
Materials 16 02689 g005
Figure 6. The effective relaxation stiffness of the yarn with no coating and twist angle.
Figure 6. The effective relaxation stiffness of the yarn with no coating and twist angle.
Materials 16 02689 g006
Figure 7. Equivalent viscoelasticity of 2D woven composite.
Figure 7. Equivalent viscoelasticity of 2D woven composite.
Materials 16 02689 g007
Figure 8. Effect of temperature on the effective relaxation stiffness of the yarn.
Figure 8. Effect of temperature on the effective relaxation stiffness of the yarn.
Materials 16 02689 g008
Figure 9. Effect of temperature on the effective relaxation modulus of 2D woven composites.
Figure 9. Effect of temperature on the effective relaxation modulus of 2D woven composites.
Materials 16 02689 g009
Figure 10. Effect of the twist angle on the effective relaxation stiffness tensors of the yarn.
Figure 10. Effect of the twist angle on the effective relaxation stiffness tensors of the yarn.
Materials 16 02689 g010
Figure 11. Effect of the twist angle on the effective relaxation moduli of 2D woven composite.
Figure 11. Effect of the twist angle on the effective relaxation moduli of 2D woven composite.
Materials 16 02689 g011
Figure 12. Effect of the coating thickness on effective relaxation stiffness of the yarn.
Figure 12. Effect of the coating thickness on effective relaxation stiffness of the yarn.
Materials 16 02689 g012
Figure 13. Effect of the coating thinkness on effective relaxation moduli of 2D woven composites.
Figure 13. Effect of the coating thinkness on effective relaxation moduli of 2D woven composites.
Materials 16 02689 g013
Table 1. Material parameters of fiber.
Table 1. Material parameters of fiber.
ParameterYoung’s Modulus (GPa)Shear Modulus (GPa)Poisson’s Ratio
E f 1 E f 2 = E f 3 G f 12 = G f 13 G f 23 υ f 12 = υ f 13 υ f 23
Value233158.9635.6390.20.33
Table 2. Relaxation times and Prony coefficients for PMT-F4 epoxy.
Table 2. Relaxation times and Prony coefficients for PMT-F4 epoxy.
i 1234567
E i ( MPa ) 1000224.1450.8406.1392.7810.4203.71486.0
ρ i ( s ) -1.0 × 1031.0 × 1051.0 × 1061.0 × 1071.0 × 1081.0 × 1091.0 × 1010
Table 3. Effective elastic properties of fiber calculated from Equation (15).
Table 3. Effective elastic properties of fiber calculated from Equation (15).
ParameterTwist Angle
θ ( ° )
Young’s Modulus (GPa)Shear Modulus (GPa)Poisson’s Ratio
E f 1 E f 2 = E f 3 G f 12 = G f 13 G f 23 υ f 12 = υ f 13 υ f 23
Value233158.9635.6390.20.33
30°179.958317.6730722.1166.80.650.299
60°53.0767437.77733.20217.0560.470.107
Table 4. Mechanical properties of the coating [30] (Reprinted/adapted with permission from Ref. [30]. 2009, Elsevier).
Table 4. Mechanical properties of the coating [30] (Reprinted/adapted with permission from Ref. [30]. 2009, Elsevier).
ParameterYoung’s Modulus (GPa)Shear Modulus (GPa)Poisson’s Ratio
E 1 int E 2 int = E 3 int G y 12 int = G y 13 int G y 23 int υ 12 int = υ 13 int υ 23 int
Pyrolytic carbon301224.30.120.4
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

Li, B.; Liu, C.; Zhao, X.; Ye, J.; Guo, F. Multiscale Study of the Effect of Fiber Twist Angle and Interface on the Viscoelasticity of 2D Woven Composites. Materials 2023, 16, 2689. https://doi.org/10.3390/ma16072689

AMA Style

Li B, Liu C, Zhao X, Ye J, Guo F. Multiscale Study of the Effect of Fiber Twist Angle and Interface on the Viscoelasticity of 2D Woven Composites. Materials. 2023; 16(7):2689. https://doi.org/10.3390/ma16072689

Chicago/Turabian Style

Li, Beibei, Cheng Liu, Xiaoyu Zhao, Jinrui Ye, and Fei Guo. 2023. "Multiscale Study of the Effect of Fiber Twist Angle and Interface on the Viscoelasticity of 2D Woven Composites" Materials 16, no. 7: 2689. https://doi.org/10.3390/ma16072689

APA Style

Li, B., Liu, C., Zhao, X., Ye, J., & Guo, F. (2023). Multiscale Study of the Effect of Fiber Twist Angle and Interface on the Viscoelasticity of 2D Woven Composites. Materials, 16(7), 2689. https://doi.org/10.3390/ma16072689

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