Macro-meso scale simulations of 3D woven composite reinforcements during the forming process

Online since 02 April 2021

Article

Abstract



During the forming stage in the RTM process, deformations and orientations of yarns at the mesoscopic scale are essential to evaluate mechanical behaviors of final composite products and calculate the permeability of the reinforcement. However, due to the high computational cost, it is very difficult to carry out a mesoscopic draping simulation for the entire reinforcement. In this paper, a macro-meso scale simulation of composite reinforcements is presented in order to predict mesoscopic deformations of the fabric in a reasonable calculation time. The proposed multi-scale method allows linking the macroscopic simulation of the reinforcement with the mesoscopic modelling of the RVE through a macromeso embedded analysis. On the base of macroscopic simulations using a hyperelastic constitutive law of the reinforcement, an embedded mesoscopic geometry is first deduced from the macroscopic simulation of the draping. To overcome the inconvenience of the macro-meso embedded solution which leads to unreal excessive yarn extensions, local mesoscopic simulations based on the embedded analysis are carried out on a single RVE by defining specific boundary conditions. Finally, the multi-scale forming simulations are investigated in comparison with the experimental results, illustrating the efficiency of the proposed approach, in terms of accuracy and CPU time.

Keywords

Macro-meso Analyses, Forming Simulation, 3D Textile Reinforcement, RVE, Hyperelasticity

Table of contents

Text

1 Introduction

In the manufacturing process of textile composites, especially the RTM process, the forming step is important because it strongly influences the mechanical behavior of the composite in service. In order to better predict the geometrical and mechanical characteristics of composite materials, numerical methods are increasingly developed in view of the duration and cost of testing. Different types of simulations can be distinguished at three different scales: macroscopic, mesoscopic and microscopic. 

For the shaping simulations of woven reinforcements, the approaches are generally carried out at the macroscopic scale [1, 2], in which the reinforcement is considered as a continuous medium. At this scale, it is able to predict the appearance of wrinkles, which is the major defect appearing during this phase [3, 4]. However, some important phenomena may appear only on a smaller scale during a forming process of reinforcements [5-8]. Mesoscopic models, which consider reinforcements as a set of yarns in contact with their neighbors, are able to predict possible defects, such as the rupture and the local buckling of yarns. In addition, the knowledge of the deformations and orientations of the yarns allows to determine the mechanical behaviour of the final part and to calculate the permeability tensor of the deformed reinforcement which can be used to simulate the resin flow in the injection stage [9-11]. Nevertheless, it is very difficult to simulate the shaping of woven reinforcements at this mesoscopic scale for the whole part by taking into account the huge number of yarns and their complex interactions, especially for thick textile composites.

This is why researchers are currently turning to macro-meso simulations of reinforcements. Gatouillat [7] has proposed a simplified model which consists in a mesoscopic scale description of the reinforcements by shell finite elements in contact, allowing the simulation of the fabric on a macroscopic scale, which can predict the wrinkling and sliding of yarns related to the decohesion of the yarns. In Iwata's study [12], visualizations of yarns at the local level of the reinforcement is achieved by coupling a macroscopic simulation of draping and a local simulation at the mesoscopic scale in the positions where the most probable defects are found in the macroscopic simulation. In these models, one of the most critical problems is to compensate for the increase in computing time. In this context, the objective of this paper is therefore to propose an efficient multi-scale method that allows to obtain the deformations and orientations of the yarns for the whole part during the shaping of the woven reinforcements, with a reasonable computation time.

The multi-scale method developed consists in creating a passage of information between the two different scales [13]. This passage makes it possible to transform the deformed information between the macroscopic and mesoscopic scales. In other words, the deformations and orientations of the yarns, as well as the variations of the yarns sections for the whole part can be calculated from macroscopic forming simulations of woven reinforcements. The method is mainly divided into three parts: macroscopic simulation of the woven reinforcements, macro-meso embedded analysis, and local mesoscopic simulation of RVE. Three types of woven reinforcements with different architectures, natures and weaves are studied, and only the data related to orthogonal 3D glass woven reinforcements (Fig. 1) will be presented here.

Image 10000000000003030000012A476F184B.pngFig. 1. Structure of the 3D non-crimp orthogonal woven fabric and its 3D architecture by X-ray tomography.

2 Macroscopic simulation of woven reinforcements

The first part consists in implementing numerical simulations of the woven reinforcements on a macroscopic scale using the PlastFib [14] numerical simulation tool developed by the LaMCoS laboratory. This finite element code uses a temporal scheme in explicit dynamics. Because of the possible slippage between fibers and yarns, the mechanical behaviors of textile reinforcements are very specific compared to traditional materials such as metals. It leads to a high tensile stiffness in the direction of the fibers compared to other stiffnesses. In macroscopic approaches, woven reinforcements are considered as a continuous medium with a constitutive law of homogeneous behavior. In order to describe the mechanical characteristics of thick textile reinforcements during draping, a 3D hyperelastic behaviour law [15, 16] is applied in macroscopic simulations.

In this law of hyperelastic behaviour, six modes of deformation are considered: elongation of the reinforcement in the warp and weft direction, transverse compression of the reinforcement, shear of the reinforcement in the plane, and transverse shear in the warp and weft direction. For each deformation mode, the corresponding physical invariant based on physical observations of the six above deformation modes can be defined by combinations of the invariants defined in Equation (1-4). I1, I2, I3 are the classical invariants of the right Cauchy-green tensor, and I41, I42 are the mixed invariants concerning to the structural tensors.

Image 1000000000000366000000BD0E4D4244.png

The stress-strain relationship is determined through the derivative of a strain energy density. Based on the decoupling hypothesis, the volume energy of deformation can then be expressed as a function of the physical invariants relative to the different modes of deformation. In addition, the contribution of curvature in the hexahedral finite elements is taken into account by adding a local flexural stiffness.

Image 10000000000004A00000002BE99A61C3.png

The characterization of the material parameters used in this law are carried out by different classical tests: Tensile test, Transverse Compression, In-Plane Shear, Transverse Shear and Bending test. The identifications are carried out by the direct method or the inverse method with the Levenberg-Marquardt algorithm. Then, the second Piola-Kirchhoff stress tensor S͇ and ̳σ Cauchy stress tensor can be calculated by the differentiation Әwk/ ӘIk:

Image 100000000000051F0000006D08662C59.png

Once the material parameters have been identified, macroscopic finite element simulations can be applied for the woven reinforcements in two cases: three-point bending and hemispherical forming (Fig. 2). Finally, in comparison with the corresponding experiments, the results obtained by simulation are both in good agreement with the experimental ones.

Image 1000000000000456000001872447C05A.pngFig. 2. Comparisons between the numerical and experimental results: a) three-point bending and b) hemispherical forming.

3 Macro-meso embedded analysis

3.1 Macro-meso embedded approach

The mesoscopic information of the deformed part is very important for the simulations of the resin flow in the injection stage. Mesoscopic analyses of the entire reinforcement have been performed [5, 12], but the size of the numerical model and the computational time are very important. The objective of the macro-meso embedded approach is to visualize the deformation of the yarns of the whole part during draping with a reasonable computational time. This developed approach proposes to perform mesoscopic analyses using the results obtained at macroscopic scale. If the macroscopic results are relevant, they should be able to provide some basic elements at the mesoscopic scale. In short, the embedded mesoscopic geometry is thus determined from the macroscopic simulation.

The macro-meso embedded approach allows to link macroscopic simulations of interlock reinforcements and mesoscopic modeling of RVEs (Representative Volume Element). For periodic textile reinforcements initially, the geometrical modeling of the RVE is first reconstructed according to X-ray tomography images based on the real geometry of the textile armor (Fig. 3) [17]. As the smallest volume of material, a RVE allows to reconstruct the whole reinforcement by repetitions and translations. The modeling of the entire reinforcement at mesoscopic scale can then be generated in the initial configuration. Then, each mesoscopic node of yarns of the woven reinforcement can be located or integrated into a macroscopic element (Fig. 4). During deformation, the integrated mesoscopic element deforms together with the macroscopic element where this mesoscopic element is located and its mesoscopic nodes have constant coordinates in the reference macroscopic element.

Image 10000000000002F3000001182FCA52EF.pngFig. 3. a) Geometry obtained by X-ray tomography of the 3D non-crimp woven reinforcement, b) Corresponding prism mesh.

Image 10000000000002BB00000121749408E2.pngFig. 4. Meso model embedded in a macro element. Real and reference frames.

Image 10000000000005120000002B312D1FD0.png

where (ξ, η ,ζ) are the naturel coordinates of the mesoscopic node in the macroscopic reference element. Xmeso and Xmacro are the spatial coordinates of the mesoscopic node and macroscopic element nodes.

During the transformation, the natural coordinates (ξ, η, ζ) of a mesoscopic node in the corresponding macroscopic element remain identical. In other words, the position of an integrated mesoscopic node remains constant in the macroscopic element. Therefore, the deformation obtained by the macroscopic simulation can provide a first deformed geometry of RVE at the mesoscopic scale in the deformed configuration.

3.2 Applications

Fig. 5 and Figure-6 show the deformed geometries of the reinforcement at the mesoscopic scale obtained by the macro-meso embedded approach in the two cases: three points bending and hemispherical stamping. In addition, the post-treatments associated with each deformation mode: elongation in the warp and weft directions, transverse compression, in-plane shear and transverse shear in the warp and weft directions, are implemented in order to visualize the deformations of yarns.

Image 1000000000000388000001704406FFAA.pngFig. 5. Mesoscopic analysis of three points bending. Physical invariants of yarns: (a) elongation, (b) compaction and (c) elongation in the binder yarn.

The above mesoscopic results on the basis of macroscopic simulations allow to determine the deformations of the yarns, their final orientations and also the voids in an elementary cell. These results can be used for the calculation of the permeability during the next step of resin injection. The calculations concern only geometrical calculations, the calculation time is very fast, from a few minutes to one hour depending on the reinforcement size. This may be sufficient for some mesoscopic scale analyses. Nevertheless, as shown in the figure, it generally has drawbacks, in particular excessive elongation of yarns. A local mesoscopic analysis below can avoid these problems.

Image 1000000000000235000001A83A6932C3.pngFig. 6. Mesoscopic analysis of hemispherical stamping. Physical invariants of yarns: compaction.

4 Local mesoscopic simulation

As mentioned above, mesoscopic analyses of the woven reinforcement obtained by the macro-meso embedded approach have shown their efficiency in providing a mesoscopic scale solution during shaping. However, the equations of mechanics at the mesoscopic scale are not strictly verified. The mesoscopic recessed deformed mesoscopic configuration is not obtained by solving a mechanical problem but is directly deduced from the macroscopic analysis, due to the assumptions on the continuous medium and the finite element approach. From the point of view of the constitutive law, the stresses do not check the local equilibrium. This approach can lead to excessive elongations of the yarns because it does not take into account local slippage between the yarns. To overcome this difficulty, a local mesoscopic simulation based on macroscopic analysis and embedded mesoscopic analysis is proposed.

The mesoscopic simulation is performed on a local area (e.g. a RVE) based on the results of the macro-meso embedded approach. The deformed configuration obtained by the embedded mesoscopic simulation constitutes the initial state of the local mesoscopic simulation via the boundary conditions in displacement. Then, the displacements imposed on all nodes are released by embedding the edges of the RVE. A relatively simple transverse isotropic hyperelastic behaviour law (neo-hookean) is used to describe the mechanical behaviour of the yarns [18, 19]. Consequently, stress equilibrium is generally achieved within the yarns. The results obtained by the local mesoscopic simulation are compared with those of the macro-meso embedded approach for the two cases: three points bending and hemispherical forming (Fig. 7 and Fig. 8).

Image 1000000000000241000001B7A80B3D18.pngFig. 7. Three points bending: mesoscopic simulations of a RVE. Elongation in a) macro-meso embedded analysis b) local mesoscopic simulation analysis.

Image 100000000000035C0000016C3DDA3879.pngFig. 8. Hemispherical forming: mesoscopic simulations of a RVE. Elongation in a) macro-meso embedded analysis b) local mesoscopic simulation analysis.

The elongation of ten elements located in different areas of this RVE is shown in Fig. 9 for three points bending and Fig. 10 for hemispherical forming. This shows that parasitic elongations in the macro-meso embedded approach are globally suppressed by the local mesoscopic simulation.

Image 100000000000039200000133B95A4856.png

Fig. 9. Elongation change of different locations of a RVE in local mesoscopic simulation for three points bending.

Image 100000000000036F00000147EB29D8E1.pngFig. 10. Elongation change of different locations of a RVE in local mesoscopic simulation for hemispherical forming.

Finally, mesoscopic numerical results are compared with experimental results. The deformations of the longitudinal and transverse sections of the yarns can be observed and analyzed by a microscope. In comparison with the experimental results (Fig. 11 b), it can be seen that the numerical results (Fig. 11 a) based on the multi-scale method are in good agreement with the experimental ones from regarding the deformed geometry.

Image 100000000000035C00000189759AB381.png

Fig. 11. Woven unit cells in the same position obtained by: a) Local mesoscopic simulation b) experiment.

5 Conclusions

The deformation information of woven reinforcements on the mesoscopic scale is essential to evaluate the mechanical properties of the final composite part, as well as to simulate resin flow during the next injection step. The first macro-meso embedded method proposed allows to obtain mesoscopic scale deformations of the woven reinforcements during forming. However, it has drawbacks, especially excessive elongation of the yarns. To overcome this problem, local mesoscopic simulation based on the results obtained by macro-meso embedded method is carried out by defining specific boundary conditions. This step takes a reasonable computational cost. The results are satisfactory by comparing with the experiment. In addition, some other work is carrying out to optimize and verifier this multi-scale method, such as the study of influences of the number and the positions of RVEs, comparison with the full meso simulation at large scale.

Bibliography

[1] Gereke, T. Döbrich, O. Hübner, M. Cherif, C. Experimental and computational composite textile reinforcement forming: a review. Compos A Appl Sci Manuf 2013, 46, 1–10.

[2] Bussetta, P. Correia, N. Numerical forming of continuous fibre reinforced composite material: a review. Compos A Appl Sci Manuf 2018, 113, 12–31.

[3] 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–92.

[4] Dangora, LM. Mitchell, CJ. Sherwood, JA. Predictive model for the detection of out- of-plane defects formed during textile-composite manufacture. Compos A Appl Sci Manuf 2015, 78, 102–12.

[5] Gatouillat, S. Bareggi, A. Vidal-Sallé, E. Boisse, P. Meso modelling for composite preform shaping–simulation of the loss of cohesion of the woven fibre network. Compos A Appl Sci Manuf. 2013, 54, 135–44.

[6] Allaoui, S. Hivet, G. Soulat, D. Wendling, A. Ouagne, P. Chatel, S. Experimental preforming of highly double curved shapes with a case corner using an interlock reinforcement. Int J Mater Form, 2014, 7(2), 155–65.

[7] Schirmaier, FJ. Dörr, D. Henning, F. Kärger, L. A macroscopic approach to simulate the forming behaviour of stitched unidirectional non-crimp fabrics (UD-NCF). Compos A Appl Sci Manuf, 2017, 102, 322–35.

[8] Schirmaier, FJ. Weidenmann, KA. Kärger, L. Henning, F. Characterisation of the draping behaviour of unidirectional non-crimp fabrics (UD-NCF). Compos A Appl Sci Manuf, 2016, 80, 28–38.

[9] Bréard, J. Henzel, Y. Trochu, F. Gauvin, R. Analysis of dynamic flows through porous media. Part I: comparison between saturated and unsaturated flows in fibrous reinforcements. Polym Compos 2003, 24(3), 391–408.

[10] Loix, F. Badel, P. Orgéas, L. Geindreau, C. Boisse, P. Woven fabric permeability: From textile deformation to fluid flow mesoscale simulations. Compos Sci Technol, 2008, 68(7–8), 1624–30.

[11] Tran, T. Comas-Cardona, S. Abriak, NE. Binetruy, C. Unified microporomechanical approach for mechanical behavior and permeability of misaligned unidirectional fiber reinforcement. Compos Sci Technol 2010, 70(9), 1410–8.

[12] Iwata, A. Inoue, T. Naouar, N. Boisse, P. Lomov, SV. Coupled meso-macro simulation of woven fabric local deformation during draping. Compos A Appl Sci Manuf,  2019, 118, 267–80.

[13]Wang, J. Wang, P. Hamila, N. Boisse, P. Mesoscopic analyses of the draping of 3D woven composite reinforcements based on macroscopic simulations. Compos Struct, 2020, 250, 112602.

[14] Software PLASFIB, Paris, Inter Deposit Certification, agence pour la protection des programmes. 2011 and 2015.

[15] Charmetant, A. Orliac, JG. Vidal-Sallé, E. Boisse, P. Hyperelastic model for large deformation analyses of 3D interlock composite preforms. Compos Sci Technol, 2012, 72(12), 1352–60.

[16] Mathieu, S. Hamila, N. Bouillon, F. Boisse, P. Enhanced modeling of 3D composite preform deformations taking into account local fiber bending stiffness. Compos Sci Technol, 2015, 117, 322–33.

[17] Naouar, N. Vidal-Salle, E. Schneider, J. Maire, E. Boisse, P. 3D composite reinforcement meso FE analyses based on X-ray computed tomography. Compos Struct, 2015, 132, 1094–104.

[18] Bonet, J. Burton, AJ. A simple orthotropic, transversely isotropic hyperelastic constitutive equation for large strain computations. Comput Methods Appl Mech Eng, 1998, 162(1–4), 151–64.

[19] Florimond, C. Contribution à la modélisation mécanique du comportement de mèches de renforts tissés à l’aide d’un schéma éléments finis implicite PhD Thesis. Lyon, France: INSA; 2015.

Illustrations

References

Electronic reference

Jie Wang, Peng Wang, Nahiène Hamila and Philippe Boisse, « Macro-meso scale simulations of 3D woven composite reinforcements during the forming process », ESAFORM 2021 [Online], Online since 02 April 2021, connection on 29 March 2024. URL : https://popups.uliege.be/esaform21/index.php?id=496

Authors

Jie Wang

Univ Lyon, INSA-Lyon, CNRS UMR5259, LaMCoS, F-69621, France.

Corresponding author: jie.wang@insa-lyon.fr

Peng Wang

University of Haut-Alsace, ENSISA, LPMT, Mulhouse, F-68000, France.

By this author

Nahiène Hamila

Univ Lyon, INSA-Lyon, CNRS UMR5259, LaMCoS, F-69621, France.

IRDL - UMR CNRS 6027 – ENIB, Brest, France.

By this author

Philippe Boisse

Univ Lyon, INSA-Lyon, CNRS UMR5259, LaMCoS, F-69621, France.

By this author

Details

Title
Macro-meso scale simulations of 3D woven composite reinforcements during the forming process
Language
en
Author, co-author
Jie Wang, Peng Wang, Nahiène Hamila and Philippe Boisse,
Publication date
02 April 2021
Journal title
ESAFORM 2021
Copyright
CC-BY
DOI
10.25518/esaform21.496
Permanent URL
https://popups.uliege.be/esaform21/index.php?id=496

Statistics

Views : 628 (1 ULiège)
Downloads : 0 (0 ULiège)

Cite

Wang, J., Wang, P., Hamila, N., & Boisse, P. (2021). Macro-meso scale simulations of 3D woven composite reinforcements during the forming process. Paper presented at ESAFORM 2021. 24th International Conference on Material Forming, Liège, Belgique. doi: 10.25518/esaform21.496