Next Article in Journal
Microstructural Assessment of a Multiple-Intermetallic-Strengthened Aluminum Alloy Produced from Gas-Atomized Powder by Hot Extrusion and Friction Extrusion
Next Article in Special Issue
Towards the Prediction of Sandwich Composites Durability in Severe Condition of Temperature: A New Numerical Model Describing the Influence of Material Water Content during a Fire Scenario
Previous Article in Journal
Research of Wood Waste as a Potential Filler for Loose-Fill Building Insulation: Appropriate Selection and Incorporation into Polyurethane Biocomposite Foams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Textile Composite Damage Analysis Taking into Account the Forming Process

1
INSA-Lyon, LaMCoS CNRS, University de Lyon, 69621 Lyon, France
2
LMA CNRS, Aix Marseille University, 13007 Marseille, France
*
Author to whom correspondence should be addressed.
Materials 2020, 13(23), 5337; https://doi.org/10.3390/ma13235337
Submission received: 17 October 2020 / Revised: 19 November 2020 / Accepted: 23 November 2020 / Published: 25 November 2020
(This article belongs to the Special Issue Organic Matrix Composites and Multifunctional Materials)

Abstract

:
The internal structure of composite materials is modified during manufacturing. The formation of woven prepregs or dry preforms changes the angle between the warp and weft yarns. The damage behaviour of the consolidated composite is modified by these changes of angle. It is important when designing a composite part to consider this modification when calculating the damage in order to achieve a correct dimensioning. In this paper, a damage calculation approach of the consolidated textile composite that takes into account the change in orientation of the yarns due to forming is proposed. The angles after forming are determined by a simulation of the draping based on a hypoelastic behaviour of the woven fabric reinforcement. Two orthogonal frames based on the warp and weft directions of the textile reinforcement are used for the objective integration of stresses. Damage analysis of the cured woven composite with non-perpendicular warp and weft directions is achieved by replacing it with two equivalent Unidirectional (UD) plies representing the yarn directions. For each ply, a model based on Continuum Damage Mechanics (CDM) describes the progressive damage. Two examples are presented, a bias extension specimen and the hemispherical forming coupon. In both cases, the angles between the warp and weft yarns are changed. It is shown that the damage calculated by taking into account these angle changes is greatly modified.

1. Introduction

When the ratio between stiffness and weight is important in the design of a structure, composite materials are a preferred solution. This is particularly the case for many aeronautical structures [1,2] as well as in the automotive industry [3,4,5]. The reinforcements and matrix can be adapted to a given application. However, the manufacturing process often changes the final state of the composite, especially when forming is performed. The directions of the reinforcement fibres are modified, both in the case of Liquid Composite Moulding processes (LCM) [6,7] as well as in the thermoforming of thermoset or thermoplastic prepregs [4,8,9]. The simulation of the forming process of composite parts has developed significantly in recent years [10,11]. These advancements make it possible to predict the state of the composite after manufacture in particular the orientation of the fibres and certain defects such as wrinkles [12,13,14]. Some studies have been conducted on the influence of the manufacturing processes on the mechanical characteristics of composites [15,16,17]. The influence of defects and in particular of wrinkles and waviness on the mechanical characteristics of composites were analysed in [18,19,20,21].
Damage is one of the major degradations of composite parts during service. Micro-cracking of the matrix can lead to failure of the composite. Damage analysis and modelling is an important aspect of composite science, especially since damage modes can be complex due to the fibre-matrix composition of composites. This paper analyses the influence of angle changes between warp and weft yarns created by the forming of the textile reinforcement on the development of damage in the consolidated composite. The simulation of the forming of a woven reinforcement is based on a hypoelastic approach. One of the main results of this simulation is the prediction of the fibre orientations after the forming process. A damage model for reinforced woven-composites is presented for the analysis of the consolidated composite. This model is based on Continuum Damage Mechanics (CDM) and can describe the progressive damage evolution of the material which consists in the propagation of small matrix cracks parallel to the fibre directions. The model, initially proposed for unidirectional plies [22], has been extended to woven plies [23,24]. The woven ply is modelled by two equivalent UD plies that correspond to the warp and the weft directions, directions that may be non-orthogonal after the forming process. After this simplification, it is possible to directly apply the classical unidirectional model [22] that describes the linear elastic behaviour in the fibre direction and the nonlinear damaged behaviour in the transverse and shear directions [23].
To our knowledge, there is a lack of simulation techniques suitable for the estimation of matrix damage in a loaded composite structure accounting for the effect of the forming process. This represents the objective of the present article in the case of textile composites. It will also be shown that this consideration of the forming process has important consequences on the calculated damage. To demonstrate the potential of coupling the two approaches of forming simulation and damage calculation, first a bias extension test (in-plane shear test obtained by traction at 45° [25,26,27]) at high temperature on a thermoplastic prepreg with a woven reinforcement is performed and simulated. Afterwards, a tensile test is performed on the consolidated composite with modified warp and weft angles due to the forming process. The damage generated in the coupon for this configuration is then simulated and analysed. In the last part, a tension test on a coupon obtained by hemispherical forming is simulated. It is shown that, in this case, the damage level in zones where the angle variations between warp and weft are large is strongly affected by the forming process.

2. Forming Process Simulation

The simulation of composite forming has been the subject of much work in recent decades. Reviews in this domain are presented in [10,11]. In LCM processes, the textile reinforcements are dry (resin-free) during the forming process, which takes place before the resin is injected. The thermoforming of prepregs is carried out at a temperature such that the matrix is melted to allow the forming. The presence of the resin in this case conditions the mechanical behaviour of the composite during forming by making it dependent on temperature and giving it a viscous character [28,29,30,31]. In all cases, the mechanical behaviour during forming is mainly driven by the reinforcement and the position of the fibres, which are quasi inextensible. Concerning simulation techniques, one can distinguish between kinematic methods [32,33,34,35] which are numerically fast but do not take into account the nature and the forces acting on the reinforcement such as those generated by a blank-holder, and approaches based on the equations of mechanics, the behaviour of the reinforcement and the frictional contact between the blank and the tools. Because the textile reinforcements are thin, some approaches consider them as membrane where the bending stiffness is considered low enough to be neglected [36,37,38,39,40,41]. It has nevertheless been shown that the bending stiffness of the textile reinforcement plays an important role in particular on the development of wrinkles [14,42,43]. In this article, the objective of the forming simulation is mainly the determination of warp and weft yarns angle variation, which is related to the membrane deformation of the textile reinforcement. Therefore, a membrane approach will be used. The textile reinforcement is assumed to be continuous and a hypoelastic behaviour model is used. Hyperelastic approaches [44,45,46] and elastic approaches [11,47,48,49] have also been developed.

2.1. An Hypoelastic Approach Based on the Warp and Weft Yarn Rotations

Hypoelastic behaviour models (also known as ‘rate constitutive equations’ [50]) have been developed to model the mechanical behaviour of materials under large deformations. These behaviour laws are frequently used and are often proposed by default in finite element codes. The constitutive law presented below is implemented in the explicit Abaqus code via a Vumat user subroutine. From the strain rate D _ _ , a constitutive tensor C _ _ _ _ gives the objective derivative of the Cauchy stress tensor. This is the strain rate for an observer fixed to the material.
σ _ _ = C _ _ _ _ : D _ _
with
σ _ _ = σ _ _ + σ _ _ . Ω _ _ Ω _ _ . σ _ _
Here, σ _ _ = d σ d t is the time derivative of the Cauchy stress. The spin Ω _ _ = Q _ _ . Q _ _ T is associated to the rotation Q _ _ of the frame that rotates with the material. The use of the objective derivative of stresses is necessary to ensure that movements of rigid bodies do not create stresses. There are several objective derivatives. The most classical are Green Naghdi and Jaumann [51,52]. It has been shown that these derivatives are not suitable for fibrous materials. In this case, the objective derivative must be that defined by the rotation of the fibre [53,54]. The directions f _ 1 and f _ 2 of warp and weft yarns during deformation are given by the initial positions of the yarns f _ 10 and f _ 20 and the deformation gradient F _ _ .
f _ 1 = F _ _ . f _ 10 F _ _ . f _ 10 f _ 2 = F _ _ . f _ 20 F _ _ . f _ 20
Two orthogonal frames based on these two yarn directions are defined: g (g1 = f1, g2) and h (h1, h2 = f2) (Figure 1).

2.2. Stress Update in the Frames of the Yarns

The strain increment between times tn and tn+1 is denoted d ε _ _ . Its components in the frames g (g1, g2) and h (h1, h2) are considered (α and β are indexes taking the value 1 or 2).
d ε _ _ = d ε α β g g _ α g _ β = d ε α β h h _ α h _ β
The incremental elongation in the warp yarn direction g _ 1 gives the incremental axial stress component in g _ 1 direction from tensile stiffness Eg.
d ε 11 g = g _ 1 . d ε _ _ . g _ 1 d σ 11 g = E g d ε 11 g
In the same way, the stress increment in the weft direction g _ 2 is obtained by:
d ε 22 h = h _ 2 . d ε _ _ . h _ 2 d σ 22 h = E h d ε 22 h
The in-plain shear stress increment is obtained from the in- plain shear strain increment and the shear modulus G.
d ε 12 g = g _ 1 . d ε _ _ . g _ 2 d ε 12 h = h _ 1 . d ε _ _ . h _ 2 d σ 12 g = G d ε 12 g d σ 12 h = G d ε 12 h
The shear modulus G is not constant and varies significantly with the shear strain [25,26,27]. The stress components are integrated between tn and tn+1 following the scheme of Hughes and Winget [49] in g (g1, g2) and h (h1, h2) frames, respectively.
( σ 11 g ) n + 1 = ( σ 11 g ) n + d σ 11 g n + 1 / 2 ( σ 12 g ) n + 1 = ( σ 11 g ) n + d σ 12 g n + 1 / 2 ( σ 22 h ) n + 1 = ( σ 22 h ) n + d σ 22 h n + 1 / 2 ( σ 12 h ) n + 1 = ( σ 11 h ) n + d σ 12 h n + 1 / 2
The Cauchy stress tensor at time tn+1 is then determined by:
σ _ _ n + 1 = ( σ _ _ g ) n + 1 + ( σ _ _ h ) n + 1
The components of this stress tensor can be computed in any frame. In particular, the Abaqus explicit solver requires that the Vumat user routine gives the components in the Green–Naghdi’s frame that rotates with the polar rotation.
In addition to the stresses, the approach described above provides the directions f _ 1 and f _ 2 of warp and weft yarns (Equation (3)). These directions, and thus the angle change between warp and weft yarns, are considered in the damage calculation of the consolidated composite presented in the next sections. The following two tests will be considered to analyse the influence of the forming and more particularly of the warp-weft angle change in the calculation of the damage of the consolidated composite. The Young modulus Eh and Eg of Equations (5) and (6) (in warp and weft directions) and the shear modulus G of Equation (7) are given in Table 1 for the woven reinforcement considered in the simulations in Section 2.2 and Section 2.3.

2.3. Bias Extension Test

The Bias Extension Test is an experiment fairly simple to perform that imposes in-plane shearing on a specimen with a textile reinforcement [25,26,27]. Warp and weft yarns must be initially oriented at ±45°. The specimen is extended creating constant shear zones (A,B,C) in Figure 2. For the test to be carried out, the length of the specimen must be greater than or equal to twice its width. For the specimen utilised here, the dimensions are 210 × 70 mm2. The prepreg utilised is a PA66 thermoplastic matrix reinforced with an 8-Harness Satin glass woven. The bias extension test can be carried out either on dry textile reinforcements (without resin) or on prepregs at a sufficiently high temperature for the resin to be molten. This last technique is applied in our case. The test is performed at 270 °C, i.e., slightly above the melting temperature of the PA66 matrix. The deformation of the specimen is obtained by a 50 mm displacement of the moving jaw of the traction machine (d in Figure 2a) which corresponds to a warp and weft angle of 38°. The specimen is left cooling obtaining the final configuration where the warp and weft yarn directions have been modified by the bias extension.

2.4. Hemispherical Forming

The simulation of a hemispherical forming is then performed (Figure 3).
The hemispherical shape is widely used to perform draping analysis for all materials [55,56] but especially for textile reinforcements of composites [12,44,57]. It is simple in shape but clearly double curved and requires large plane shear deformations for its shape to be achieved [57,58,59,60,61,62]. This forming process leads to angles between warp and weft yarns, which evolve with the position between 0° and 45° (Figure 4). This will allow the influence of a variable angle between warp and weft to be analysed below.
The two tests presented in this section, bias extension test and hemispherical forming, both strongly modify the angles between warp and weft yarns. The objective of the following sections is to assess the influence of these modifications on the damage behaviour of the consolidated composite.

3. Damage Model for Woven Ply Laminates

This section concerns the cured composite after the manufacturing process. The failure of laminate composite structures is due to many mechanisms acting at different scales and depending on the type of ply, woven or unidirectional (UD), concern the matrix or the fibres and are a function of many parameters (thickness, orientations, loads, etc.) [63,64,65,66,67,68,69,70]. Matrix damage usually starts with fibre/matrix decohesions that propagate in to form transverse cracks over the entire thickness of the ply. In the case of UD plies, these transverse cracks can spread over great lengths and can initiate delamination [71,72,73,74]. In the case of woven plies, these transverse cracks develop over the whole thickness of the yarns, but their propagation is blocked by the crossing of the warp and weft fibres. The models based on Continuum Damage Mechanics (CDM) [22] can described the progressive damage of the material. These CDM models are suitable for woven composites because matrix cracks are distributed more regularly than for UD plies [23]. Also, laminates made of woven plies exhibit better resistance to delamination. It is therefore possible to remain within the framework of the theory of laminates in plane stresses and CDM to describe the evolution of matrix damage for static and fatigue loads [75,76].
A simplified model has been proposed to describe the behaviour of woven plies within the framework of 2D laminate theory. In order to describe the different evolutions of matrix damage in the warp and weft, the woven ply is replaced by two Equivalent UD (EqUD) plies whose thicknesses respect the ratios between the number of warp and weft yarns. In the case of a balanced woven ply, the thicknesses of the EqUD plies are the same. In the case of unbalanced woven the thicknesses of the EqUD plies will be different. This modelling strategy allows for the simulation of woven reinforcement where the warp and weft yarns are non-perpendicular, this way the orientations of the yarns can correspond to the ones obtained after forming; it is this original point that will be used in this paper. After this simplification, it is possible to directly apply classical UD models, [22] for example, describing the linear elastic behaviour in the fibre direction and the damaged non-linear matrix behaviour in the transverse and shear directions.

3.1. Equivalent Unidirectional Plies Model of the Woven Ply

The woven ply is replaced by two Equivalent Unidirectional (EqUD) plies with different thicknesses to take into account the different ratios of fibres in the warp and the weft directions (Figure 5). This model applies to both balanced and unbalanced woven fabrics because it takes into account a parameter δ which represents the percentage of fibres in the warp direction, 1-δ therefore being the percentage in the weft direction. In this study, we consider only balanced woven plies, so δ will be always equal to 0.5.
The elastic parameters and the damage law coefficients of the EqUD plies were evaluated from the mechanical response of the woven ply laminates. The identification of the material properties required by the model is detailed in [23,77] for an unbalanced woven glass/epoxy plies and static loads. The elastic and damage laws have been tested and validated for balanced woven carbon/epoxy plies in static [76,78] and extended for fatigue loads [75,76]. The proposed simplified model of using two equivalent UD plies to model the behaviour of a woven ply has been used recently by other authors [79,80,81].

3.2. Matrix Damage Model at the Unidirectional Ply Scale

3.2.1. Assumptions

A model based on Continuum Damage Mechanics (CDM) has been developed to describe the damage evolution in composite materials at the ply scale [22].
Material damage can be described by its effects on elastic coefficients [82,83]. In the case of a ply consisting of continuous fibres bound by a brittle epoxy matrix, matrix cracks parallel to the direction of the fibres will appear very early, spreading rapidly throughout the thickness of the ply and also in the plane (Figure 5a). The density of these cracks will gradually increase with loading and lead to a decrease in transverse and shear moduli at the ply scale. While the development of these matrix cracks is gradual and generally not catastrophic for the laminate, the evolution of fibre cracks is much more abrupt. In the context of CDM, it is conventional to introduce internal damage variables di to describe these decreases:
E 1 = E 1 0 ( 1 d 1 )   ;   E 2 = E 2 0 ( 1 d 2 )   ;   G 12 = G 12 0 ( 1 d 12 )
where E 1 0 ,   E 2 0   and   G 12 0 are the initial moduli and the damage variables d1, d2, d12 vary from 0 for the initial material to 1 for the fully degraded material:
-
d1 whose evolution represents the linear elastic behaviour and the brittle fracture of the fibres observed during a tensile test in the longitudinal direction of the fibres;
-
d2 representing the effect of the matrix damage on the stiffness in the transverse direction;
-
d12 representing the effect of the matrix damage on the shear stiffness.
The damage is assumed to be uniform in the thickness of the ply.
The progressive development of the damage variables d2 and d12 depends on the tensile and shear stresses, which generate cracks in the matrix. Under the assumption of plane stress and small perturbations, the strain energy in the ply can be written as follows [22,23]:
E D p s = 1 2 [ σ 1 + 2 E 1 0 ( 1 d 1 ) + σ 1 2 E 1 0 + σ 2 + 2 E 2 0 ( 1 d 2 ) + σ 2 2 E 2 0 2 ν 12 0 E 1 0 σ 1 σ 2 + σ 12 2 G 12 0 ( 1 d 12 ) ]
where < . >+ and < . >- are the MacAuley brackets which select respectively the positive or the negative part of the quantity that is within. This way, the tensile energy and the compressive energy are separated in order to describe the unilateral nature of the damage progression given by the opening of the cracks in tension and their closing in compression. The thermodynamic forces associated with the internal tensile and shear variables d1, d2 and d12 are defined as follows:
{ Y d i = E D p s d i = σ i + 2 2 E i 0 ( 1 d i ) 2 w i t h i = 1 , 2 Y d 12 = E D p s d 12 = ( σ 12 ) 2 2 G 12 0 ( 1 d 12 ) 2

3.2.2. Damage Evolution Laws

The development of the internal damage variables depends on these thermodynamic forces and, more specifically, their maximum value during the loading history.
Under tensile load, d1 increases suddenly to model the brittleness behaviour in the direction of the fibres.
By noting Y 1 m a x the resistance in the fibre direction, we then have:
d 1 = 0   if   . Y d 1 < Y 1 m a x   ;   d 1 = 1   otherwise .
In the case of static loading, the tension/shear coupling during the progressive development of d2 is accounted for by the following equivalent thermodynamic force:
Y e q = a ( Y d 2 ) n + b ( Y d 12 ) m
where a, b, m and n are material parameters specifying the tension/shear coupling. The evolution law for the damage is written as:
d 2 = 1 e ( Y e q Y 0 s ) +
where the constant   Y 0 s corresponds to the threshold value for the development of d2.
As proposed in [23,77], the shear damage d12 is taken proportional to the transverse damage d2. This choice is based on the fact that the cracks are parallel to the fibres and their effect on the transverse and shear modulus is the same.
d 12 = c d 2
By considering that the moduli are proportional to the effective surfaces (initial surfaces minus the cracked areas) [83], one can take as a first approximation the coefficient c = 1. This coefficient probably varies as a function of the crack density and numerical micromechanical analyses [81] showed that it is less than 1 but increases for high damage levels.

3.2.3. Inelastic Strain in the Shear Direction

After loading has been applied to a [45,−45]s laminate for example, inelastic strains are observed. These strains may result from the slipping/friction process occurring between fibre and matrix as the result of the damage. Although inelastic strain can also be observed during tensile test on 0° and 90° laminates, only the inelastic strain in the shear direction is relevant. A kinematic hardening model was used to describe the inelastic shear strain evolution.
The coupling between the damage and the plasticity is accounted for by the effective stress and the effective strain [22,23] which are respectively written as:
σ ˜ 12 = σ 12 ( 1 d 12 ) and ε ˜ 12 p = ε 12 p ( 1 d 12 )
It is assumed that stresses σ 11 and σ 22 do not influence the elastic field domain defined by:
f = | σ ˜ 12 h ε ˜ 12 p | R 0
Corresponding to a kinematic linear hardening law where R0 is the initial inelastic strain threshold and h is the linear law coefficient.
The identification of parameters of the damage model presented above are given in details in [23]. The identification of the behaviour typically requires three loading-unloading-reloading tests: a test at 0° for the behaviour in the fibre direction (Figure 5b) [84], a test at 90° for the transverse behaviour and a test on a [45,−45] s laminate to identify the behaviour in shear (Figure 5b). The numerical values of the parameters adopted in the simulations are given in the Table 2. The papers [23,84] present several validation tests on different laminates for static loads. The model has been extended to fatigue loads [75,76]. In this case, the matrix damage evolves during the cycles and the model has been identified and validated. It should be noted that for each type of fibre, resin or weaving, it is necessary to repeat these tests to identify the damage model coefficients of the equivalent UD ply.
The part of the model presented here relates mainly to the nonlinear evolution of the matrix damage for static loadings. Laminate failure is not discussed here. This more complex part is presented in the cited papers. Catastrophic laminate failure is rarely due to matrix damage localization mechanisms because the laminate usually consists of plies oriented in the loading direction. On the other hand, the fibre failure is catastrophic for the laminate and the structure. The approach used in this case is described in [85] and takes into account the resistance decrease in the fibre direction as a function of the level of matrix damage and also takes into account the increase in fibre strength in the presence of stress concentrations using non-local criteria [84,86]. It should be noted that these non-local criteria make it possible to regularize the problem and thus avoids the mesh-size dependence for the numerical approaches.

4. Damage Analysis Taking into Account the Yarns Angle Variations Due to Forming

4.1. Bias Extension Test

A tension test simulation was performed on the consolidated specimen obtained after a bias extension test. As described in Section 2.3, after forming the specimens presents different yarn orientations, the original [45,−45] s layup is conserved only at the extremities and a [19,−19] s orientation is found in the central zone. To perform the simulation the left end is fixed (Figure 6). The right end is subjected to a horizontal force of 7000 N and vertical displacement are blocked.

4.1.1. Transverse Damage d2

Figure 6 presents the matrix transverse damage map. It can be observed that the damage is very low in the central zone (dark blue), where the orientation of the fibres is [19,−19] s with respect to the loading direction. In this area the orientation of the yarns is close to the loading direction, therefore, the transverse and shear stresses responsible for matrix damage are the lowest. In the zones at the ends of the test specimen, the yarn angle remains unchanged after forming and therefore close to the initial orientation of [45,−45] s. In these zones, where the loading direction is far from the orientation of the fibres, the damage reaches its maximum values in the order of 0.5, which corresponds to a loss of rigidity of 50% (in green) because the shear stresses at the scale of the EqUD plies are maximum.
In this example, the damage shows a strong variation depending on the angles between warp and weft fibres and the loading direction. It should be noted that with the described damage model, two damage maps are obtained: one for the EqUD ply corresponding to the warp direction and another for the EqUD ply to the weft direction. As the applied force follows the bisector between the original warp and weft directions, the results obtained for the two EqUD plies are very similar for this example.

4.1.2. Analysis of Inelastic Strains

For the studied case inelastic strains are present only in the zones of strong yarn angle variation especially at the meeting point between the A and B zones where the yarn orientation is [45,−45] s (Figure 7).

4.2. Damage Analysis on a Specimen from Hemispherical Forming

A rectangular specimen is extracted from the flat zone of the part shaped by hemispherical stamping. This specimen is subjected to tension and the damage is analysed. In this coupon, the angle between the yarns of the first ply varies gradually but strongly (Zone A, in Figure 8). The boundary conditions are the same as in the case of the bias extension test. The left end of the specimen is fixed, and the right end is subjected to a tension force while transverse displacements are blocked.
Figure 9 clearly shows that the direction of the fibres plays a major role in the development of damage when the specimen is subjected to tension. The left part where the warp and weft yarns have a direction quite close to that of the loading direction is little or not damaged. The damage starts and then becomes significant in the right side of the specimen where the weft and warp directions are close to [45,−45] s. This example shows that taking into account changes in orientation due to forming is essential in the damage analysis of composites structures obtained by forming techniques. We can see that the strongly nonlinear evolution of the damage which increases rapidly in the right part of the specimen when the force varies only from 4000 to 5000 N.
The inelastic strains associated with the damage also develop in the right part where the warp and wefts yarns are in directions close to ±45° (Figure 10) and, like the damage, this evolution is nonlinear since the strains double for an increase in force of 25%.
Figure 11 displays the non-linear nature of the force-displacement curve due to damage and inelastic strain, especially for imposed forces in excess of 2000 N.

5. Discussion and Conclusions

In this study, a simulation model for the prediction of yarns angle after forming process was coupled with a damage model for woven composite materials. The association of the two models demonstrated the importance of taking into account the change in angle between warp and weft yarns due to the manufacturing process when estimating the damage behaviour of composite structures. The forming simulation is based on a hypoelastic model for the textile reinforcement and it provides the yarn directions and therefore the warp and weft angle changes after the process.
The analysis of the damage behaviour of the cured textile composite is carried out considering two UD equivalent plies in the yarns’ directions. It has thus been shown that this approach to the estimation of matrix damage makes it possible to take into account the reorientation of the fibres consequence of the forming process. The influence of this reorientation is very important and must be considered for a realistic analysis.
This study is part of the work carried out with the aim of taking into account the manufacturing processes when calculating composite structures. One of the first future developments concerns the extension of damage analysis to curved areas (single or double curvature) which are numerous in composite parts. In addition, experimental investigations are necessary to confirm and validate the influence of the manufacturing process on the damage.

Author Contributions

Conceptualization, P.B. and C.H.; methodology, M.J., A.C. and N.N.; software N.N., A.C., M.J., validation M.J.; supervision. P.B. and C.H.; original draft preparation, P.B., C.H. and A.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Middendorf, P.; Metzner, C. Aerospace applications of non-crimp fabric composites. In Non-Crimp Fabric Composites-Manufacturing, Properties and Applications; Lomov, S.V., Ed.; Woodhead Publishing: Cambridge, UK, 2011; pp. 441–449. [Google Scholar]
  2. Irving, P.E.; Soutis, C. Polymer Composites in the Aerospace Industry; Woodhead Publishing: Cambridge, UK, 2019. [Google Scholar]
  3. Friedrich, K.; Almajid, A.A. Manufacturing aspects of advanced polymer composites for automotive applications. Appl. Compos. Mater. 2013, 20, 107–128. [Google Scholar] [CrossRef]
  4. Henning, F.; Kärger, L.; Dörr, D.; Schirmaier, F.J.; Seuffert, J.; Bernath, A. Fast processing and continuous simulation of automotive structural composite components. Compos. Sci. Technol. 2019, 171, 261–279. [Google Scholar] [CrossRef]
  5. Martulli, L.M.; Creemers, T.; Schöberl, E.; Hale, N.; Kerschbaum, M.; Lomov, S.V.; Swolfs, Y. A thick-walled sheet moulding compound automotive component: Manufacturing and performance. Compos. Part A Appl. Sci. Manuf. 2020, 128, 105688. [Google Scholar] [CrossRef]
  6. Sozer, E.M.; Simacek, P.; Advani, S.G. Resin transfer molding (RTM) in polymer matrix composites. In Manufacturing Techniques for Polymer Matrix Composites (PMCs); Woodhead Publishing: Cambridge, UK, 2012; pp. 245–309. [Google Scholar]
  7. Bickerton, S.; Govignon, Q.; Kelly, P. Resin infusion/liquid composite moulding (LCM) of advanced fibre-reinforced polymer (FRP). In Advanced Fibre-Reinforced Polymer (FRP) Composites for Structural Applications; Woodhead Publishing Series in Civil and Structural Engineering: Cambridge, UK, 2013; pp. 155–186. [Google Scholar]
  8. Lukaszewicz, D.J.; Potter, K.D. The internal structure and conformation of prepreg with respect to reliable automated processing. Compos. Part A Appl. Sci. Manuf. 2011, 42, 283–292. [Google Scholar] [CrossRef]
  9. Dirk, H.J.L.; Ward, C.; Potter, K.D. The engineering aspects of automated prepreg layup: History, present and future. Compos. Part B Eng. 2012, 43, 997–1009. [Google Scholar]
  10. Gereke, T.; Döbrich, O.; Hübner, M.; Cherif, C. Experimental and computational composite textile reinforcement forming: A review. Compos. Part A Appl. Sci. Manuf. 2013, 46, 1–10. [Google Scholar] [CrossRef]
  11. Bussetta, P.; Correia, N. Numerical forming of continuous fibre reinforced composite material: A review. Compos. Part A Appl. Sci. Manuf. 2018, 113, 12–31. [Google Scholar] [CrossRef]
  12. Skordos, A.A.; Aceves, C.M.; Sutcliffe, M.P. A simplified rate dependent model of forming and wrinkling of pre-impregnated woven composites. Compos. Part A Appl. Sci. Manuf. 2007, 38, 1318–1330. [Google Scholar] [CrossRef] [Green Version]
  13. Boisse, P.; Hamila, N.; Vidal-Sallé, E.; Dumont, F. Simulation of wrinkling during textile composite reinforcement forming. Influence of tensile, in-plane shear and bending stiffnesses. Compos. Sci. Technol. 2011, 71, 683–692. [Google Scholar] [CrossRef] [Green Version]
  14. Lightfoot, J.S.; Wisnom, M.R.; Potter, K. A new mechanism for the formation of ply wrinkles due to shear between plies. Compos. Part A Appl. Sci. Manuf. 2013, 49, 139–147. [Google Scholar] [CrossRef]
  15. Vieille, B.; Albouy, W.; Chevalier, L.; Taleb, L. About the influence of stamping on thermoplastic-based composites for aeronautical applications. Compos. Part B Eng. 2013, 45, 821–834. [Google Scholar] [CrossRef]
  16. Vieille, B.; Albouy, W.; Taleb, L. Influence of stamping on the compressive behavior and the damage mechanisms of C/PEEK laminates bolted joints under severe conditions. Compos. Part B Eng. 2015, 79, 631–638. [Google Scholar] [CrossRef]
  17. Aridhi, A.; Arfaoui, M.; Mabrouki, T.; Naouar, N.; Denis, Y.; Zarroug, M.; Boisse, P. Textile composite structural analysis taking into account the forming process. Compos. Part B Eng. 2019, 166, 773–784. [Google Scholar] [CrossRef] [Green Version]
  18. Potter, K.; Khan, B.; Wisnom, M.; Bell, T.; Stevens, J. Variability, fibre waviness and misalignment in the determination of the properties of composite materials and structures. Compos. Part A Appl. Sci. Manuf. 2008, 39, 1343–1354. [Google Scholar] [CrossRef]
  19. Mahadik, Y.; Hallett, S.R. Effect of fabric compaction and yarn waviness on 3D woven composite compressive properties. Compos. Part A Appl. Sci. Manuf. 2011, 42, 1592–1600. [Google Scholar] [CrossRef]
  20. Bloom, L.D.; Wang, J.; Potter, K.D. Damage progression and defect sensitivity: An experimental study of representative wrinkles in tension. Compos. Part B Eng. 2013, 45, 449–458. [Google Scholar] [CrossRef]
  21. Bender, J.J.; Hallett, S.R.; Lindgaard, E. Investigation of the effect of wrinkle features on wind turbine blade sub-structure strength. Compos. Struct. 2019, 218, 39–49. [Google Scholar] [CrossRef] [Green Version]
  22. Ladeveze, P.; Ledantec, E. Damage modelling of the elementary ply for laminated composites. Compos. Sci. Technol. 1992, 43, 257–267. [Google Scholar] [CrossRef]
  23. Thollon, Y.; Hochard, C. A general damage model for woven fabric composite laminates up to first failure. Mech. Mater. 2009, 41, 820–827. [Google Scholar] [CrossRef]
  24. Payan, J.; Hochard, C. Damage modelling of laminated carbon/epoxy composites under static and fatigue loadings. Int. J. Fatigue 2002, 24, 299–306. [Google Scholar] [CrossRef]
  25. Lebrun, G.; Bureau, M.N.; Denault, J. Evaluation of bias-extension and picture-frame test methods for the measurement of intraply shear properties of PP/glass commingled fabrics. Compos. Struct. 2003, 61, 341–352. [Google Scholar] [CrossRef]
  26. Harrison, P.; Clifford, M.J.; Long, A.C. Shear characterisation of viscous woven textile composites: A comparison between picture frame and bias extension experiments. Compos. Sci. Technol. 2004, 64, 1453–1465. [Google Scholar] [CrossRef]
  27. Boisse, P.; Hamila, N.; Guzman-Maldonado, E.; Madeo, A.; Hivet, G.; Dell’Isola, F. The bias-extension test for the analysis of in-plane shear properties of textile composite reinforcements and prepregs: A review. Int. J. Mater. Form. 2017, 10, 473–492. [Google Scholar] [CrossRef]
  28. McGuinness, G.B.; ÓBrádaigh, C.M. Development of rheological models for forming flows and picture-frame shear testing of fabric reinforced thermoplastic sheets. J. Non-Newton. Fluid Mech. 1997, 73, 1–28. [Google Scholar] [CrossRef]
  29. Haanappel, S.P.; Ten Thije, R.H.W.; Sachs, U.; Rietman, B.; Akkerman, R. Formability analyses of uni-directional and textile reinforced thermoplastics. Compos. Part A Appl. Sci. Manuf. 2014, 56, 80–92. [Google Scholar] [CrossRef]
  30. Guzman-Maldonado, E.; Hamila, N.; Boisse, P.; Bikard, J. Thermomechanical analysis, modelling and simulation of the forming of pre-impregnated thermoplastics composites. Compos. Part A Appl. Sci. Manuf. 2015, 78, 211–222. [Google Scholar] [CrossRef]
  31. Dörr, D.; Joppich, T.; Kugele, D.; Henning, F.; Kärger, L. A coupled thermomechanical approach for finite element forming simulation of continuously fiber-reinforced semi-crystalline thermoplastics. Compos. Part A Appl. Sci. Manuf. 2019, 125, 105508. [Google Scholar] [CrossRef]
  32. Long, A.C.; Rudd, C.D. A simulation of reinforcement deformation during the production of preforms for liquid moulding processes. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 1994, 208, 269–278. [Google Scholar] [CrossRef]
  33. Wang, J.; Paton, R.; Page, J.R. The draping of woven fabric preforms and prepregs for production of polymer composite components. Compos. Part A Appl. Sci. Manuf. 1999, 30, 757–765. [Google Scholar] [CrossRef]
  34. Hancock, S.G.; Potter, K.D. The use of kinematic drape modelling to inform the hand lay-up of complex composite components using woven reinforcements. Compos. Part A Appl. Sci. Manuf. 2006, 37, 413–422. [Google Scholar] [CrossRef]
  35. Cherouat, A.; Borouchaki, H.; Giraud-Moreau, L. Mechanical and geometrical approaches applied to composite fabric forming. Int. J. Mater. Form. 2010, 3, 1189–1204. [Google Scholar] [CrossRef]
  36. Gelin, J.C.; Cherouat, A.; Boisse, P.; Sabhi, H. Manufacture of thin composite structures by the RTM process: Numerical simulation of the shaping operation. Compos. Sci. Technol. 1999, 56, 711–718. [Google Scholar] [CrossRef]
  37. Xue, P.; Peng, X.; Cao, J. A non-orthogonal constitutive model for characterizing woven composites. Compos. Part A Appl. Sci. Manuf. 2003, 34, 183–193. [Google Scholar] [CrossRef]
  38. Yu, W.R.; Zampaloni, M.; Pourboghrat, F.; Chung, K.; Kang, T.J. Analysis of flexible bending behavior of woven preform using non-orthogonal constitutive equation. Compos. Part A Appl. Sci. Manuf. 2005, 36, 839–850. [Google Scholar] [CrossRef]
  39. Ten Thije, R.H.W.; Akkerman, R.; Huétink, J. Large deformation simulation of anisotropic material using an updated Lagrangian finite element method. Comput. Methods Appl. Mech. Eng. 2007, 196, 3141–3150. [Google Scholar] [CrossRef]
  40. Chen, S.; Harper, L.T.; Endruweit, A.; Warrior, N.A. Formability optimisation of fabric preforms by controlling material draw-in through in-plane constraints. Compos. Part A Appl. Sci. Manuf. 2015, 76, 10–19. [Google Scholar] [CrossRef]
  41. Chen, S.; McGregor, O.P.L.; Harper, L.T.; Endruweit, A.; Warrior, N.A. Defect formation during preforming of a bi-axial non-crimp fabric with a pillar stitch pattern. Compos. Part A Appl. Sci. Manuf. 2016, 91, 156–167. [Google Scholar] [CrossRef]
  42. Boisse, P.; Colmars, J.; Hamila, N.; Naouar, N.; Steer, Q. Bending and wrinkling of composite fiber preforms and prepregs. A review and new developments in the draping simulations. Compos. Part B Eng. 2018, 141, 234–249. [Google Scholar] [CrossRef]
  43. Huang, J.; Boisse, P.; Hamila, N.; Zhu, Y. Simulation of Wrinkling during Bending of Composite Reinforcement Laminates. Materials 2020, 13, 2374. [Google Scholar] [CrossRef]
  44. Aimène, Y.; Vidal-Sallé, E.; Hagège, B.; Sidoroff, F.; Boisse, P. A hyperelastic approach for composite reinforcement large deformation analysis. J. Compos. Mater. 2010, 44, 5–26. [Google Scholar] [CrossRef]
  45. Peng, X.; Guo, Z.; Du, T.; Yu, W.R. A simple anisotropic hyperelastic constitutive model for textile fabrics with application to forming simulation. Compos. Part B Eng. 2013, 52, 275–281. [Google Scholar] [CrossRef]
  46. Yao, Y.; Huang, X.; Peng, X.; Liu, P.; Youkun, G. An anisotropic hyperelastic constitutive model for plain weave fabric considering biaxial tension coupling. Text. Res. J. 2019, 89, 434–444. [Google Scholar] [CrossRef]
  47. Boisse, P.; Cherouat, A.; Gelin, J.C.; Sabhi, H. Experimental study and finite element simulation of a glass fiber fabric shaping process. Polym. Compos. 1995, 16, 83–95. [Google Scholar] [CrossRef]
  48. Hamila, N.; Boisse, P. A meso–macro three node finite element for draping of textile composite preforms. Appl. Compos. Mater. 2007, 14, 235–250. [Google Scholar] [CrossRef]
  49. Boisse, P.; Hamila, N.; Helenon, F.; Hagège, B.; Cao, J. Different approaches for woven composite reinforcement forming simulation. Int. J. Mater. Form. 2008, 1, 21–29. [Google Scholar] [CrossRef]
  50. Hughes, T.J.; Winget, J. Finite rotation effects in numerical integration of rate constitutive equations arising in large-deformation analysis. Int. J. Numer. Methods Eng. 1980, 15, 1862–1867. [Google Scholar] [CrossRef]
  51. Dafalias, Y.F. Corotational rates for kinematic hardening at large plastic deformations. Trans. Asme J. Appl. Mech. 1983, 50, 561–565. [Google Scholar] [CrossRef]
  52. Dienes, J.K. On the analysis of rotation and stress rate in deforming bodies. Acta Mech. 1979, 32, 217–232. [Google Scholar] [CrossRef]
  53. Badel, P.; Vidal-Sallé, E.; Boisse, P. Large deformation analysis of fibrous materials using rate constitutive equations. Comput. Struct. 2008, 86, 1164–1175. [Google Scholar] [CrossRef]
  54. Khan, M.A.; Mabrouki, T.; Vidal-Sallé, E.; Boisse, P. Numerical and experimental analyses of woven composite reinforcement forming using a hypoelastic behaviour. Application to the double dome benchmark. J. Mater. Process. Technol. 2010, 210, 378–388. [Google Scholar] [CrossRef]
  55. Efrati, E.; Sharon, E.; Kupferman, R. Elastic theory of unconstrained non-Euclidean plates. J. Mech. Phys. Solids 2009, 57, 762–775. [Google Scholar] [CrossRef] [Green Version]
  56. Cogun, F.; Darendeliler, H. Comparison of different yield criteria in various deep drawn cups. Int. J. Mater. Form. 2017, 10, 85–98. [Google Scholar] [CrossRef]
  57. Cherouat, A.; Billoët, J.L. Mechanical and numerical modelling of composite manufacturing processes deep-drawing and laying-up of thin pre-impregnated woven fabrics. J. Mater. Process. Technol. 2001, 118, 460–471. [Google Scholar] [CrossRef]
  58. Cao, J.; Xue, P.; Peng, X.; Krishnan, N. An approach in modeling the temperature effect in thermo-stamping of woven composites. Compos. Struct. 2003, 61, 413–420. [Google Scholar] [CrossRef]
  59. Ten Thije, R.H.W.; Akkerman, R. A multi-layer triangular membrane finite element for the forming simulation of laminated composites. Compos. Part A Appl. Sci. Manuf. 2009, 40, 739–753. [Google Scholar] [CrossRef]
  60. Duhovic, M.; Mitschang, P.; Bhattacharyya, D. Modelling approach for the prediction of stitch influence during woven fabric draping. Compos. Part A Appl. Sci. Manuf. 2011, 42, 968–978. [Google Scholar] [CrossRef]
  61. Bardl, G.; Nocke, A.; Hübner, M.; Gereke, T.; Pooch, M.; Schulze, M.; Cherif, C. Analysis of the 3D draping behavior of carbon fiber non-crimp fabrics with eddy current technique. Compos. Part B Eng. 2018, 132, 49–60. [Google Scholar] [CrossRef]
  62. Liu, K.; Zhang, B.; Xu, X.; Ye, J. Experimental characterization and analysis of fiber orientations in hemispherical thermostamping for unidirectional thermoplastic composites. Int. J. Mater. Form. 2019, 12, 97–111. [Google Scholar] [CrossRef]
  63. Laurin, F.; Carrère, N.; Maire, J.-F. A multiscale progressive failure approach for composite laminates based on thermodynamical viscoelastic and damage models. Compos. Part A Appl. Sci. Manuf. 2007, 38, 198–209. [Google Scholar] [CrossRef]
  64. Wisnom, M.R.; Atkinson, J.W. Reduction in tensile and flexural strength of unidirectional glass fibre-epoxy with increasing specimen size. Compos. Struct. 1997, 38, 405–411. [Google Scholar] [CrossRef]
  65. Melro, A.R.; Camanho, P.P.; Andrade Pires, F.M.; Pinho, S.T. Micromechanical analysis of polymer composites reinforced by unidirectional fibres: Part II—Micromechanical analyses. Int. J. Solids Struct. 2013, 50, 1906–1915. [Google Scholar] [CrossRef] [Green Version]
  66. Reifsnider, K. Durability and damage tolerance of fibrous composite systems. In Handbook of Composites; Springer: Berlin/Heidelberg, Germany, 1998; pp. 794–809. [Google Scholar]
  67. Talreja, R. Fatigue of composite materials: Damage mechanisms and fatigue-life diagrams. Proc. R. Soc. Lond. A 1981, 378, 461–475. [Google Scholar]
  68. Jen, M. Strength and life in thermoplastic composite laminates under static and fatigue loads. Part I: Experimental. Int. J. Fatigue 1998, 20, 605–615. [Google Scholar] [CrossRef]
  69. Shokrieh, M.M.; Lessard, L.B. Multiaxial fatigue behaviour of unidirectional plies based on uniaxial fatigue experiments—I. Modelling. Int. J. Fatigue 1997, 19, 201–207. [Google Scholar] [CrossRef]
  70. Gao, Z. A Cumulative Damage Model for Fatigue Life of Composite Laminates. J. Reinf. Plast. Compos. 1994, 13, 128–141. [Google Scholar]
  71. Tohgo, K.; Wang, A.S.D.; Chou, T.-W. A Criterion for Splitting Crack Initiation in Unidirectional Fiber-Reinforced Composites. J. Compos. Mater. 1993, 27, 1054–1076. [Google Scholar] [CrossRef]
  72. Ladevèze, P.; Lubineau, G.; Violeau, D. A Computational Damage Micromodel of Laminated Composites. Int. J. Fract. 2006, 137, 139–150. [Google Scholar] [CrossRef]
  73. Hochard, C.; Aubourg, P.-A.; Charles, J.-P. Modelling of the mechanical behaviour of woven-fabric CFRP laminates up to failure. Compos. Sci. Technol. 2001, 61, 221–230. [Google Scholar] [CrossRef]
  74. Camanho, P.P.; Davila, C.G.; de Moura, M.F. Numerical Simulation of Mixed-Mode Progressive Delamination in Composite Materials. J. Compos. Mater. 2003, 37, 1415–1438. [Google Scholar] [CrossRef]
  75. Hochard, C.; Thollon, Y. A generalized damage model for woven ply laminates under static and fatigue loading conditions. Int. J. Fatigue 2010, 32, 158–165. [Google Scholar] [CrossRef]
  76. Hochard, C.; Payan, J.; Bordreuil, C. A progressive first ply failure model for woven ply CFRP laminates under static and fatigue loads. Int. J. Fatigue 2006, 28, 1270–1276. [Google Scholar] [CrossRef]
  77. Thollon, Y. Analyse du Comportement à Rupture de Composites Stratifiés Constitués de Plis Tissés sous Chargements Statique et de Fatigue. Ph.D. Thesis, Aix-Marseille Université, Marseille, France, July 2009. [Google Scholar]
  78. Miot, S. Rupture de Structures Composites Stratifiées sous Chargements Statique et de Fatigue. Ph.D. Thesis, Aix-Marseille Université, Marseille, France, July 2009. [Google Scholar]
  79. Pannier, Y.; Foti, F.; Gigliotti, M. High temperature fatigue of carbon/polyimide 8-harness satin woven composites. Part I: Digital Image Correlation and Micro-Computed Tomography damage characterization. Compos. Struct. 2020, 244. [Google Scholar] [CrossRef]
  80. Sevenois, R.D.B.; Garoz, D.; Gilabert, F.A.; Hochard, C.; Van Paepegem, W. Influence of tab debonding on measured stiffness evolution in Compression-Compression and Tension-Compression fatigue testing of short gauge length coupons. Compos. Sci. Technol. 2019, 180, 1–13. [Google Scholar] [CrossRef]
  81. 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]
  82. Ladeveze, P. A damage computational approach for composites: Basic aspects and micromechanical relations. Comput. Mech. 1995, 17, 142–150. [Google Scholar] [CrossRef]
  83. Lemaitre, J.; Chaboche, J.L. Mechanics of Solid Materials; Cambridge University Press: Cambridge, UK, 1994. [Google Scholar]
  84. Miot, S.; Hochard, C.; Lahellec, N. A non-local criterion for modelling unbalanced woven ply laminates with stress concentrations. Compos. Struct. 2010, 92, 1574–1580. [Google Scholar] [CrossRef]
  85. Hochard, C.; Miot, S.; Thollon, Y. Fatigue of laminated composite structures with stress concentrations, Part B. Engineering 2014, 65, 11–16. [Google Scholar]
  86. Hochard, C.; Miot, S.; Lahellec, N.; Mazerolle, F.; Herman, M. Behaviour up to rupture of woven ply laminate structures under static loading conditions. Compos. Part A 2009, 40, 1017–1102. [Google Scholar] [CrossRef]
Figure 1. Orthogonal frames g (g1 = f1, g2) and h (h1, h2 = f2) defined from the direction of the warp and weft yarns.
Figure 1. Orthogonal frames g (g1 = f1, g2) and h (h1, h2 = f2) defined from the direction of the warp and weft yarns.
Materials 13 05337 g001
Figure 2. Bias extension test. (a) initial position; (b) development of sheared zone (C), half-sheared zone (B) and non-sheared zone (A); (c) finite element simulation and experimental thermoforming; (d) experiments.
Figure 2. Bias extension test. (a) initial position; (b) development of sheared zone (C), half-sheared zone (B) and non-sheared zone (A); (c) finite element simulation and experimental thermoforming; (d) experiments.
Materials 13 05337 g002aMaterials 13 05337 g002b
Figure 3. Hemispherical forming: (a) initial textile blank and (b) tools.
Figure 3. Hemispherical forming: (a) initial textile blank and (b) tools.
Materials 13 05337 g003
Figure 4. Shear angles after a hemispherical forming.
Figure 4. Shear angles after a hemispherical forming.
Materials 13 05337 g004
Figure 5. (a) Damage in a composite with woven reinforcement [23]. (b) Stresses versus strains in tension in different directions in experiments and model. (c) Assumption for woven plies.
Figure 5. (a) Damage in a composite with woven reinforcement [23]. (b) Stresses versus strains in tension in different directions in experiments and model. (c) Assumption for woven plies.
Materials 13 05337 g005
Figure 6. Boundaries condition for the simulation of tension test and simulated damage map.
Figure 6. Boundaries condition for the simulation of tension test and simulated damage map.
Materials 13 05337 g006
Figure 7. Inelastic shear strain of a specimen obtained via bias extension test after loading in tension.
Figure 7. Inelastic shear strain of a specimen obtained via bias extension test after loading in tension.
Materials 13 05337 g007
Figure 8. A specimen from the hemispherical forming (Zone A).
Figure 8. A specimen from the hemispherical forming (Zone A).
Materials 13 05337 g008
Figure 9. Matrix damage d2 for for 4000 and 5000 N.
Figure 9. Matrix damage d2 for for 4000 and 5000 N.
Materials 13 05337 g009
Figure 10. Inelastic shear strain for 4000 and 5000 N.
Figure 10. Inelastic shear strain for 4000 and 5000 N.
Materials 13 05337 g010
Figure 11. Non-linear force-displacement curve.
Figure 11. Non-linear force-displacement curve.
Materials 13 05337 g011
Table 1. Mechanical properties of the textile reinforcement during forming.
Table 1. Mechanical properties of the textile reinforcement during forming.
Young’s modulus in warp direction Eg:35,400 MPa
Young’s modulus in weft direction Eh:35,400 MPa
Shear modulus G (the shear angle γ is in radians) G = 6.7135 γ 4 9.8228 γ 3 + 6.3822 γ 2 1.5928 γ + 0.1948 GPa
Table 2. Damage properties of the cured textile composite [80].
Table 2. Damage properties of the cured textile composite [80].
Initial Young’s modulus in the fibre direction in a UD virtual ply: E 1 0 = 40,000 MPa
Initial Young’s modulus in the transverse direction in a UD virtual ply: E 2 0 = 13,000 MPa
Initial Shear modulus in a UD virtual ply: G 12 0 = 4000 MPa
Poisson ratio of UD virtual ply: ν 12 = 0.25
Parameters of the damage law:a = 0.9, b = 0.32, c = 1, m = 0.75, n = 0.75, Y0 = 0 MPa
Inelastic strain in the shear directionR0 = 60 MPa, h = 4500 MPa.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jauffret, M.; Cocchi, A.; Naouar, N.; Hochard, C.; Boisse, P. Textile Composite Damage Analysis Taking into Account the Forming Process. Materials 2020, 13, 5337. https://doi.org/10.3390/ma13235337

AMA Style

Jauffret M, Cocchi A, Naouar N, Hochard C, Boisse P. Textile Composite Damage Analysis Taking into Account the Forming Process. Materials. 2020; 13(23):5337. https://doi.org/10.3390/ma13235337

Chicago/Turabian Style

Jauffret, Marjorie, Aldo Cocchi, Naim Naouar, Christian Hochard, and Philippe Boisse. 2020. "Textile Composite Damage Analysis Taking into Account the Forming Process" Materials 13, no. 23: 5337. https://doi.org/10.3390/ma13235337

APA Style

Jauffret, M., Cocchi, A., Naouar, N., Hochard, C., & Boisse, P. (2020). Textile Composite Damage Analysis Taking into Account the Forming Process. Materials, 13(23), 5337. https://doi.org/10.3390/ma13235337

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