Material Modelling of Fabric Deformation in Forming Simulation of Fiber-Metal Laminates – A Review on Modelling Fabric Coupling Mechanisms

Online since 12 April 2021



During forming of complex fiber-metal laminates (FML), compressive stress zones occur. In pure textile forming, these compressive stresses typically lead to extensive wrinkling. In FML forming, however, wrinkling is partly hindered by the metal layers. Thus, combined stress states occur, where compression influences the deformation. In forming simulation, these compressive stresses can lead to erroneous formation of shear bands within the fabric layer, if the deformation behavior is not modelled correctly. Simple fabric models neither consider interactions between roving directions nor model interactions between membrane strains and shear strains. More advanced invariant-based hyperelastic material models are able to capture these interactions, but only consider tension and shear, while disregarding compression. A common assumption is to set the fabric compression stiffness close to zero. Experimentally, the in-plane fabric compression stiffness has not been determined so far. However, in FML forming, the compression stiffness and the combined compression-tension-shear behavior becomes relevant. In this article, the authors summarize and analyze the capacity of state-of-the-art fabric material models to predict the deformation behavior of fabrics under combined loading. Based on these findings, conclusions are drawn for a new macroscopic modeling approach for woven fabrics, including coupling of tension, compression and shear.


Fabric Forming, Fabric Behavior, Fabric Material Modelling, FE-Forming Simulation, Fiber-metal Laminates, FML, Hybrid, Deep Drawing

Table of contents


1 Introduction

In recent years, sustainability and climate protection have become key issues. This is also apparent in the design process of lightweight structures. Modern lightweight constructions are realized in a multi-material design. Therefore, joining or combining materials of different material classes has become a significant part of the development process. A special type of multi material design are fiber-metal laminates (FML), which combine the advantageous properties of metals and fiber reinforced composites in a layer-wise material. However, this makes the manufacturing of parts from this material an enormous challenge. Classical FMLs, like GLARE, were developed for the aviation industry and are expensive in their production [1]. For mass production of FML parts, the costs of material and process have to be decreased and the producible complexity of the parts has to be increased.

One possible process could be the combination of deep drawing and thermoplastic resin transfer molding (T-RTM) [2]. The deep drawing of the metal sheet, the draping and the infiltration of the fabric are performed simultaneously in this one-step process (Fig. 1). In previous work, the authors presented a finite element (FE) model and a parametric study to simulate this process [3, 4]. They have shown that the common known compressive stresses in deep drawing processes require an in-plane compression sensitive membrane modelling of the fabric. Without accounting this in-plane compression sensitive fabric behavior, unrealistic deformation modes occur, like shear bands. At first, this paper summarizes the membrane behavior of fabrics and then it reviews mechanical macroscopic models for fabrics. Finally, conclusions are drawn for a new fabric material model, which can account for in-plane compressive stresses and full membrane coupling.

Image 10000201000006400000028B3BE11AA9309271EF.png

Fig. 1. Process flow of combined deep drawing and resin transfer molding

2 Mechanical behavior of woven fabrics

A woven fabric is a plane textile structure produced by the technique of weaving. The consideration and modelling of woven fabrics can be done on three scales. On the microscale, several thousand individual fibers form a fiber bundle. These bundles are called threads or rovings and are assumed to be homogeneous at the meso-level. The most common form of fabrics are two-dimensional fabrics consisting of interwoven fiber bundles in two directions. At the macroscale, individual rovings are not considered, but the fabric is assumed to be a homogeneous solid instead [5]. Rovings in the direction of production are referred to as warp yarns (here 0° or a), while those running perpendicular to them are referred to as weft yarns (here 90° or b). The waviness of the rovings that results from the weaving is called crimp. Despite the resulting three-dimensional structure, this type of weave provides only in-plane reinforcement and is therefore referred to as two-dimensional [6]. The proportions of straight and curved roving sections change with the weave type, thereby influencing the achievable fiber volume fraction and the macroscopic material behavior significantly [7, 9].

Due to the structure of fabrics, their macroscopic material behavior differs from that of continuous solids. The fiber directions cause a strong anisotropy of the material. Boisse et al. [10] have shown that macroscopic nonlinear behavior is essentially attributable to geometry changes at the meso-level and depends only to a small extent on nonlinearities in the contact points. The crimped and interacting rovings also cause the material behavior under tensile loading to differ from that under compressive loading. This tension-compression anisotropy has hardly been explored so far.

2.1 Tension behavior

Some working groups are dedicated to the computational determination of the effective stiffness taking into account the fabric architecture. Starting from a geometric model at the meso-level, material properties depending on the local fiber orientation and fiber volume fraction are assigned to the constituents (fiber, matrix and interface). Then, using classical laminate theory (CLT), the effective macro-level material properties are determined [6, 11]. Instead of using CLT, the properties can also be determined using the "theory of effective moduli" [9]. Other studies differ mainly in their assumptions on the cross-sectional geometry of the fibers. Dixit and Mali [7] give a detailed overview of the different approaches to predict the mechanical properties, especially the tensile stiffness. However, none of these analytical approaches takes into account the actual processes during deformation. For a better understanding of the fabric material behavior under tensile loading, it is important to analyze the processes at the meso-level, since here the behavior is mainly determined by geometric changes. In their work, Kashani et al. [12] describe the macroscale behavior due to mesoscale processes in plain weave fabrics under different loads. Based on experimental data, different coupling mechanisms and influences can be inferred. Due to the interwoven structure of the rovings, interactions occur between them. Strain in one fiber direction consequently induces not only longitudinal normal stresses but also stresses in the transverse direction. This biaxial behavior leads to an effect corresponding to transverse contraction in homogeneous solids. This is evidenced by experimental results of biaxial tensile tests on fabrics with different strain ratios between the two roving directions [10]. The uniaxial stress-strain curve without transverse pre-tensioning can be divided into two regions, as shown in Fig. 2. Both regions are attributed to different physical phenomena.

Image 1000020100000633000002576D259804771D89A1.png

Fig. 2. Characteristic tensile stress-strain diagram of a fabric according to Kashani et al. [12]

First, under tensile loading, the crimping in the rovings is reduced along the tensile direction. The reduction of crimping in one direction causes increased waviness of the transverse yarns. This effect is often referred to as "crimp interchange". Under uniaxial tensile loading without hindrance to transverse contraction, this change in crimping can easily occur. The slightly ascending uniaxial tension curve at small strains describes this smoothing mechanism. After the reduction of crimping in the rovings along the tensile direction, they are stretched as the load increases. Due to the high tensile stiffness of the fibers, there is a sharp increase in tensile stresses with increasing strain. Transverse pre-strain makes crimp interchange more difficult and thus increases the energy required to achieve a given longitudinal strain. In the stress-strain diagram, this is noticeable by a steeper progression in the smoothing zone. In the region of large longitudinal strains (εi > 3 %), experimental studies show a more compliant behavior for fabrics with increasing transverse pre-tensioning. One explanation for this could be that with increasing transverse pre-stretch, the crimping in the rovings is reduced to a lesser extent along the tensile direction, resulting in a more compliant strain behavior [12]. Since large displacements in fabric forming are realized mainly by large shear deformations, only the range of small longitudinal strains is often considered [5, 13].

In further experimental investigations, the influence of the current shear deformation on the tensile behavior was analyzed. Compared to the influence of transverse strains, the effects of shear deformation on normal stresses are very small. Accordingly, the biaxial coupling has the largest influence on the tensile behavior and should therefore be considered in a material model for draping simulations [12].

2.2 Compression behavior

There have been no published studies on the material behavior of fabrics under in-plane compressive loading. Various publications report that the compressive stiffness is very low compared to the tensile stiffness, which is why it is assumed to be almost zero in many numerical calculations [14, 15, 16, 17, 18].

2.3 Shear behavior

Shear deformation is the central deformation mechanism in fabrics, especially in doubly curved structures [19]. Idealized, the contact points of warp and weft rovings act as hinge points allowing only in-plane rotations, justifying the assumption of pure shear. The shear angle γ describes the angular change between the two roving directions. Fig. 3 shows the characteristic shear stress curve without the influence of normal strains along the fiber directions. It can be divided into three zones.

Image 10000201000006300000026FAA872F1DA52023CF.png

Fig. 3. Characteristic shear stress – shear angle diagram according to Launay et al. [20]

In the first zone (I), pure rigid-body rotations of the rovings take place. The low shear force can be attributed solely to the friction in the contact points. As soon as the rovings come into lateral contact with each other, an increase in force follows, which is visible in zone II. The higher force can be explained by the lateral compaction of the rovings. The transition between zones II and III is marked by the critical shear angle γL, which is often also referred to as the shear-locking angle. The state of maximum compaction is called shear-locking and is characteristic for the third region (III). Further shear deformation is possible only under large force application. In addition, high compressive stresses occur across the fibers, which often leads to wrinkling due to the low bending stiffness of fabrics [20, 21]. Characterizing the shear behavior of fabrics using analytical models is difficult, so complex meso-level simulations [22] or experimental studies [19, 20, 23] are mainly used. The paper published by Cao et al. [23] contains benchmark tests of various research groups characterizing the shear behavior using picture frame (PF) and bias extension test (BET). They explain the necessity to use normalization methods to compensate the influence of different test setups on the results. Among other things, preconditioning is suggested, because the fabric shows a stiffer behavior during the first deformation cycle compared to further cycles. For the description of industrial applications, the behavior in the first cycle is decisive, since the forming processes occur without preconditioning [23].

No information is currently available on the shear behavior under compressive loading along the fiber directions, so further investigations are necessary for this area. Therefore, only the influences of positive normal strains are described in the following. Similar to the tensile behavior, the shear behavior shows strain-dependence. In their work, Nosrat-Nezami et al. [19] demonstrate that an increase in membrane strains increases the shear stiffness of carbon fiber fabrics. In the shear stress diagram, this effect is shown by a steeper curve as well as a sharp initial increase of shear stress. This initial jump can be attributed to the fact that for strains in the fiber direction, the rovings interact more strongly with the transversely extending rovings. The resulting larger contact forces at the crossing points of the rovings increase the friction between warp and weft yarns, thus complicating the rotations, which dominate the behavior in the range of small shear angles. Biaxial pre-stretching of the fabric maximizes this effect. For very large shear angles in the range of the shear-locking angle, the influence of the normal strains decreases [12, 20]. Thus, the coupling is nonlinear, and the actual coupling behavior depends on the fabric structure. For draping simulations, the range of small and medium shear angles is of particular interest. Furthermore, globally applied strains along the rovings increase the resistance to wrinkling, whereas local strains facilitate wrinkling [19]. The influence of normal strains on the shear behavior is particularly important in manufacturing processes that operate with high forming and blank holder forces [20]. The dependence of the shear behavior on the normal strains is called the tension-shear coupling and is the strongest coupling in the deformation behavior of fabrics. Its influence even exceeds that of biaxial coupling, so that consideration of tension-shear coupling in material models is necessary [12].

2.4 Bending behavior

The bending behavior of fabrics influences the formation and shape of wrinkles during forming. If fiber structures, such as woven fabrics, are subjected to a bending load, relative displacements occur between the individual, tensile-stiff fibers. Accordingly, the bending behavior of a fabric differs significantly from that of a continuous solid. Fiber sliding leads to a comparatively lower bending stiffness, so that severe bending is possible without fiber damage [24]. Boisse et al. [25] have shown that classical plate or shell theories, which establish a relationship between tensile and bending stiffness, cannot realistically describe the bending behavior of fabrics. In order to provide an acceptable representation of bending behavior in numerical models, bending and membrane behavior are usually modeled independently [7, 8]. For this reason, the bending behavior is not taken into account in the further consideration of this work.

3 Review on mechanical macro-models for woven fabrics

In macro-models, individual rovings are not represented and the entire fabric is considered as a continuum. This requires a lower computational effort compared to mesoscopic models and thus allows the numerical calculation of larger models. The geometric effects on the mesoscopic level described in Section 2 are not reproduced on the macro-level. For this reason, the nonlinear macroscopic material behavior must be embedded in the constitutive equations. For this purpose, two approaches are mainly followed in the literature.

3.1 Cauchy-elastic material models

The first approach is based on the definition of Cauchy-elastic material properties along the material directions in the covariant reference frame. The Cauchy-stresses are coupled directly with the strains. Since fabrics can undergo large shear deformations during forming, it follows that a non-orthogonal constitutive model is needed [16]. The core of the work by Xue et al. [16] is the decoupled description of tensile and shear properties in the covariant coordinate system along the fiber directions. The material properties are first represented as parameters in the constitutive equations. It is also assumed that the material is incompressible. Using transformation matrices, the material equations are transformed into the initial orthogonal reference frame. FE simulations are then performed with the model, and the material parameters are determined via comparison with experimental studies. With this model, the representation of nonlinear behavior is possible, but coupling mechanisms are not considered. Peng and Cao [14] use an alternative representation of stresses and strains under unchanged assumptions. Furthermore, they implement a method to track fiber directions to guarantee alignment of material behavior along them. In further work by Lee et al. [22], the model is extended individually to include membrane strain-dependent shear behavior and biaxial tensile behavior, respectively. The parameters describing these couplings are determined from simulation results of mesoscopic models. Subsequently, the material model is applied in FE simulations at the macro-level and compared with experimental data [27, 28].

The FABRIC model available in ABAQUS also uses the Cauchy-elastic approach. Experimental data on uniaxial tensile behavior and shear behavior must be passed as input variables. Alternatively, the material properties can be provided via a subroutine. With the help of an algorithm that tracks the fiber directions, the material properties are always aligned along the rovings. It is possible to specify different data for loading and unloading. A consideration of couplings is not provided, but could be implemented via the subroutine VFABRIC [29].

3.2 Hyperelastic material models

The second frequently used approach is based on hyperelastic material modeling. The most important component of hyperelastic material models is the formulation of a potential function from which a stress measure can be calculated by derivation with respect to a kinematic quantity. In some cases, this type of modeling allows to set up the material equations with respect to the initial configuration, so that no complex coordinate transformations have to be performed. All here presented models use the strain energy density W as potential function and formulate it according to

Image 10000201000006270000006009F142D7D2422509.png

as a function of the invariants of the right Cauchy-Green tensor C (RCG) in combination with the initial fiber directions a0, b0. The invariants are described by

Image 1000020100000617000000AE7BE767073D01316D.png

Here, the invariants I4 and I8 represent the quadratic stretch 𝜆i2 (i = a, b) of one roving direction and I10 is the shear angle between the roving directions in accordance to the initial orientation. It is also assumed that the total strain energy density W can be decomposed additively into a tension-compression dependent component for each fiber direction WTi and a shear component WS according to

Image 100002010000061600000039C0840CB241D8EDF6.png

as shown in Fig. 4. Forming the derivative of the strain energy density according to the RCG tensor yields the second Piola-Kirchhoff stresses. It is also assumed that there is no energy dissipation due to friction. The first hyperelastic material models consider neither the in-plane compression strains, nor coupling mechanisms. In each case, the tensile-dependent components of the distortion energy density depend only on the elongation in the fiber direction, and the shear-dependent component is determined by an invariant associated with the shear angle. The parameters of the individual components are determined by means of uniaxial tensile tests and picture frame tests [17, 30, 31]. Based on these publications, further work has extended the model to include the influence of viscous matrix [15, 32], the biaxial coupling of the tensile-dependent components [33, 34, 35], and the tension-shear coupling [18, 32, 36]. Depending on the considered coupling, uni- and biaxial tensile tests as well as picture frame tests with and without membrane strains are necessary to determine the parameters. The models already show good agreement with experimental data, even for complex geometries. Model quality increases with the representation of coupling mechanisms.

Image 10000201000006330000020C36D3FE8816D3E422.png

Fig. 4. Energy density functions WT and WS for biaxial coupling with positive elongation according to Yao et al. [33 - 36] calculated with data from Komeili et al. [5]

3.3 Other approaches

Galliot and Luchsinger [37] propose another macroscale material model. It is based on an orthotropic modeling approach fitted via experimentally obtained data. The material model focuses on the biaxial coupling under tension loading while the complex shear behavior is not represented. Boisse et al. [38] use a newly developed shell element for FE simulations in various publications [10, 21, 39, 40]. The element takes into account the biaxial tensile behavior, but the influence of the membrane strains on the shear behavior is not shown [21, 40]. Another work follows a semi discrete approach, where a finite element consists of a certain number of unit cells. This approach ranks between the meso- and macro-levels [39]. In addition, Iwata et al. [41] propose a combined approach of a meso- and a macro-model in their work. They refine a macro-level numerical model with local numerical calculations at the meso-level.

A publication that represents both, the transverse strain-dependent tension behavior and the membrane strain-dependent shear behavior, is not yet known. Moreover, the deformation behavior under combined in-plane compression has not been considered yet in macroscopic material modeling.

4 Conclusions and research demands for macroscopic modeling of fabrics with in-plane compression and full membrane coupling

The combined deep drawing and T-RTM process requires high blank holder forces for forming of the stiff metal sheets. Without these, the metal sheet would buckle. This induces combined stress states with in-plane compression stresses in the intermediate fabric layers. To model this complex multiaxial deformation behavior, a material model is needed, which considers the tension-compression anisotropy, the biaxial coupling resulting from crimp interchange of the rovings and the shear stiffening or softening due to membrane tensile or compressive strains. A conceivable approach would be a hyperelastic, invariant-based modelling approach, e.g., based on the work of Yao et al. [33, 34, 35, 36]. However, the equations contributing to the strain energy density have to be extended to negative strains. Furthermore, the equations have to be modified to fit the more complex behavior, where care has to be taken on the mathematical continuity in case of a piecewise description of the material behavior. Depending on the level of the material equations (stiffness level, stress-strain level or energy level), multiple continuously differentiable functions are necessary. Based on these findings, the authors intend to develop such a model in future work and implement a corresponding material subroutine in Abaqus/Explicit.


The German Research Foundation (DFG) kindly supports this research project (HE 6154/4-1, HE 6154/4-2). Moreover, the authors thank their colleagues for support, especially Constantin Krauß and Christian Poppe.


[1] Sinke, J. Manufacturing of GLARE Parts and Structures. Applied Composite Materials 2003;10(4/5):293–305.

[2] Mennecart, T. Werner, H. Khalifa, N.B. Weidenmann, K.A. Developments and Analyses of Alternative Processes for the Manufacturing of Fiber Metal Laminates. Proceedings of the 13th Manufacturing Science and Engineering Conference (MSEC2018), June 18-22, 2018, College Station, Texas, USA.

[3] Werner, H.O. Dörr, D. Henning, F. Kärger, L. Numerical modeling of a hybrid forming process for three-dimensionally curved fiber-metal laminates. AIP Conference Proceedings 2019;2113(1):20019.

[4] Werner, H.O. Poppe, C. Henning, F. Kärger, L. Material Modeling in Forming Simulation of Three-Dimensional Fiber-Metal-Laminates – A Parametric Study. Procedia Manufacturing, Volume 47, 2020, Pages 154-161.

[5] Komeili, M. Multi-Scale Characteization and Modeling of Shear-Tension Interaction in Woven Fabrics for Composite Forming and Structural Applications. Dissertation, The University of British Columbia, Okanagan, 2014.

[6] Adumitroaie, A. Barbero, E.J. Beyond plain weave fabrics – I. geometrical model. Composite Structures, 93(5):1424–1432, 2011. ISSN 02638223. doi: 10.1016/j.compstruct.2010.11.014.

[7] Dixit, A. Mali, H.S. Modeling techniques for predicting the mechanical properties of woven-fabric textile composites: a review. Mechanics of Composite Materials, 49(1):1–20, 2013. ISSN 1573-8922. doi: 10.1007/s11029-013-9316-8.

[8] Dörr, D.; Schirmaier, F.; Henning, F.; Kärger, L.: A viscoelastic approach for modeling bending behavior in Finite Element forming simulation of continuously fiber reinforced composites. Composites Part A 94: 113-123, 2017.

[9] Aitharaju, V.R. Averill, R.C. Three-dimensional properties of woven-fabric composites. Composites Science and Technology, (59):1901–1911, 1999. ISSN 02663538.

[10] Boisse, P. Buet, K. Gasser, A. Launay, J. Meso/macro-mechanical behaviour of textile reinforcements for thin composites. Composites Science and Technology, 61(3):395–401, 2001. ISSN 02663538. doi: 10.1016/s0266-3538(00)00096-8.

[11] Adumitroaie, A. Barbero, E.J. Beyond plain weave fabrics – II. Mechanical properties. Composite Structures, 93(5):1449–1462, 2011. ISSN 02638223. doi: 10.1016/j.compstruct.2010.11.016.

[12] Haghi Kashani, M. Hosseini, A. Sassani, F. Ko, F.K. Milani, A.S. Understanding different types of coupling in mechanical behavior of woven fabric reinforcements:A critical review and analysis. Composite Structures, 179:558–567, 2017. ISSN02638223. doi: 10.1016/j.compstruct.2017.06.069.

[13] Hivet, G. Boisse, P. Consistent mesoscopic mechanical behaviour model for woven composite reinforcements in biaxial tension. Composites Part B: Engineering, 39(2):345–361, 2008. ISSN 13598368. doi: 10.1016/j.compositesb.2007.01.011.

[14] Peng, X.Q. Cao, J. A continuum mechanics-based non-orthogonal constitutive model for woven composite fabrics. Composites Part A: Applied Science and Manufacturing, 36(6):859–874, 2005. ISSN 1359835X. doi: 10.1016/j.compositesa.2004.08.008.

[15] Gong, Y. Peng, X. Yao, Y. Guo, Z. An anisotropic hyperelastic constitutive model for thermoplastic woven composite prepregs. Composites Science and Technology, 128:17–24, 2016. ISSN 02663538. doi: 10.1016/j.compscitech.2016.03.005.

[16] Xue, P. Peng, X. Cao, J. A non-orthogonal constitutive model for characterizing woven composites. Composites Part A: Applied Science and Manufacturing, (34):183–193, 2003. ISSN 1359835X.

[17] Peng, X.Q. Guo, Z.Y. Zia-Ur-Rehman, Harrison, P. A simple anisotropic fiber reinforced hyperelastic constitutive model for woven composite fabrics. International Journal of Material Forming, 3(S1):723–726, 2010. ISSN 1960-6206. doi:10.1007/s12289-010-0872-3.

[18] Gong, Y. Yan, D. Yao, Y. Wei, R. Hu, H. Xu, P. Peng, X. An anisotropic hyperelastic constitutive model with tension–shear coupling for woven composite reinforcements. International Journal of Applied Mechanics, 09(06):1750083, 2017. ISSN 1758-8251. doi: 10.1142/S1758825117500831.

[19] Nosrat-Nezami, F. Gereke, T. Eberdt, C. Cherif, C. Characterization of the shear–tension coupling of carbon-fibre fabric under controlled membrane tensions for precise simulative predictions of industrial preforming processes. Composites Part A: Applied Science and Manufacturing, 67:131–139, 2014. ISSN 1359835X. doi:10.1016/j.compositesa.2014.08.030.

[20] Launay, J. Hivet, G. Duong, A.V. Boisse, P. Experimental analysis of the influence of tensions on in plane shear behaviour of woven composite reinforcements. Composites Science and Technology, 68(2):506–515, 2008. ISSN 02663538. doi:10.1016/j.compscitech.2007.06.021.

[21] Boisse, P. Zouari, B. Daniel, J.-L. Importance of in-plane shear rigidity in finite element analyses of woven fabric composite preforming. Composites Part A: Applied Science and Manufacturing, 37(12):2201–2212, 2006. ISSN 1359835X. doi:10.1016/j.compositesa.2005.09.018.

[22] Lee, W. Cao, J. Badel, P. Boisse, P.. Non-orthogonal constitutive model for woven composites incorporating tensile effect on shear behavior. International Journal of Material Forming, 1(S1):891–894, 2008. ISSN 1960-6206. doi: 10.1007/s12289-008-0239-1.

[23] Cao, J. Akkerman, R. Boisse, P. Chen, J. Cheng, H.S. de Graaf, E.F. Gorczyca, J.L. Harrison, P. Hivet, G. Launay, J. Lee, W. Liu, L. Lomov, S.V. Long, A. de Luycker, E. Morestin, F. Padvoiskis, J. Peng, X.Q. Sherwood, J. Stoilova, Tz. Tao, X.M. Verpoest, I. Willems, A. Wiggers, J. Yu, T.X. Zhu, B. Characterization of mechanical behavior of woven fabrics: Experimental methods and benchmark results. Composites Part A: Applied Science and Manufacturing, 39(6):1037–1053, 2008. ISSN1359835X. doi: 10.1016/j.compositesa.2008.02.016.

[24] 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. Composites Part B: Engineering, 141:234–249, 2018. ISSN13598368. doi:10.1016/j.compositesb.2017.12.061.

[25] Boisse, P. Bai, R. Colmars, J. Hamila, N. Liang, B. Madeo, A. The need to use generalized continuum mechanics to model 3D textile composite forming. Applied Composite Materials, 25(4):761–771, 2018. ISSN 0929189X. doi: 10.1007/s10443-018-9719-8.

[26] Kärger, L. Bernath, A. Fritz, F. Galkin, S. Magagnato, D. Oeckerath, A. Schön, A. Henning, F. Development and validation of a CAE chain for unidirectional fibre reinforced composite components. Composite Structures, 132:350–358, 2015. ISSN 02638223. doi:10.1016/j.compstruct.2015.05.047.

[27] Lee, W. Byun, J.-H. Um, M.-K. Cao, J. Boisse, P. Coupled non-orthogonal constitutive model for woven fabric composites. ICCM-17 - 17th International Conference on Composite Materials, 2009.

[28] Lee, W. Cao, J. Numerical simulations on double-dome forming of woven composites using the coupled non-orthogonal constitutive model. International Journal of Material Forming, 2(S1):145–148, 2009. ISSN 1960-6206. doi: 10.1007/s12289-009-0499-4.

[29] Dassault Systèmes. SIMULIA User Assistance 2018: Fabric material behavior, 23/12/2020.

[30] Aimène, Y. Vidal-Sallé, E. Hagège, B. Sidoroff, F. Boisse, P. A hyperelastic approach for composite reinforcement large deformation analysis. Journal of Composite Materials, 44(1):5–26, 2010. doi: 10.1177/0021998309345348.

[31] Peng, X. Guo, Z. Du, T. Yu, W.-R. A simple anisotropic hyperelastic constitutive model for textile fabrics with application to forming simulation. Composites Part B: Engineering, 52:275–281, 2013. ISSN 13598368. doi: 10.1016/j.compositesb.2013.04.014.

[32] Gong, Y. Xu, P. Peng, X. Wei, R. Yao, Y. Zhao, K. A lamination model for forming simulation of woven fabric reinforced thermoplastic prepregs. Composite Structures, 196:89–95, 2018. ISSN 02638223. doi: 10.1016/j.compstruct.2018.05.004.

[33] Yao, Y. Huang, X. Peng, X. Gong, Y. An anisotropic constitutive model with biaxial-tension coupling for woven composite reinforcements. In AIP Conference Proceedings 1769. American Institute of Physics. doi: 10.1063/1.4963562.

[34] Yao, Y. Youkun, G. Peng, X. Huang, X. Modeling biaxial-tension coupling in woven fabric. In Proceedings of the 21st International Conference on Composite Materials, 2017. URL

[35] Yao, Y. Huang, X. Peng, X. Liu, P.I. Youkun, G. An anisotropic hyperelastic constitutive model for plain weave fabric considering biaxial tension coupling. Textile Research Journal, 89(3):434–444, 2019. ISSN 0040-5175. doi:10.1177/0040517517748495.

[36] Yao, Y. Peng, X. Gong, Y. Influence of tension–shear coupling on draping of plain weave fabrics. Journal of Materials Science, 54(8):6310–6322, 2019. ISSN 0022-2461. doi: 10.1007/s10853-019-03334-w.

[37] Galliot, C. Luchsinger, R.H. A simple model describing the non-linear biaxial tensile behaviour of PVC-coated polyester fabrics for use in finite element analysis. Composite Structures, 90(4):438–447, 2009. ISSN 02638223. doi: 10.1016/j.compstruct.2009.04.016.

[38] Boisse, P. Gasser, A. Hivet, G. Analyses of fabric tensile behaviour: Determination of the biaxial tension-strain surfaces and their use in forming simulations. Composites Part A: Applied Science and Manufacturing, 32(10):1395–1414, 2001. ISSN 1359835X.

[39] Boisse, P. Hamila, N. Helenon, F. Hagege, B. Cao, J. Different approaches for woven composite reinforcement forming simulation. International Journal of Material Forming, 1(1):21–29, 2008. ISSN 1960-6206. doi: 10.1007/s12289-008-0002-7.

[40] 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. Composites Science and Technology, 71(5):683–692, 2011. ISSN 02663538. doi: 10.1016/j.compscitech.2011.01.011.

[41] Iwata, A. Inoue, T. Naouar, N. Boisse, P. Lomov, S.V. Coupled meso-macro simulation of woven fabric local deformation during draping. Composites Part A: Applied Science and Manufacturing, 118:267–280, 2019. ISSN 1359835X. doi:10.1016/j.compositesa.2019.01.004.



Electronic reference

Henrik O. Werner, Florian Schäfer, Frank Henning and Luise Kärger, « Material Modelling of Fabric Deformation in Forming Simulation of Fiber-Metal Laminates – A Review on Modelling Fabric Coupling Mechanisms », ESAFORM 2021 [Online], Online since 12 April 2021, connection on 28 May 2024. URL :


Henrik O. Werner

Karlsruhe Institute of Technology (KIT), Institute of Vehicle System Technology, Karlsruhe, Germany

Karlsruhe Institute of Technology (KIT), Institute of Applied Materials, Karlsruhe, Germany

Corresponding author:

Florian Schäfer

Karlsruhe Institute of Technology (KIT), Institute of Vehicle System Technology, Karlsruhe, Germany

Corresponding author: florian.schä

Frank Henning

Karlsruhe Institute of Technology (KIT), Institute of Vehicle System Technology, Karlsruhe, Germany

Fraunhofer Institute for Chemical Technology (ICT), Pfinztal, Germany

By this author

Luise Kärger

Karlsruhe Institute of Technology (KIT), Institute of Vehicle System Technology, Karlsruhe, Germany

By this author


Material Modelling of Fabric Deformation in Forming Simulation of Fiber-Metal Laminates – A Review on Modelling Fabric Coupling Mechanisms
Author, co-author
Henrik O. Werner, Florian Schäfer, Frank Henning and Luise Kärger,
Publication date
14 April 2021
Journal title
Permanent URL


Views : 818 (3 ULiège)
Downloads : 0 (0 ULiège)


Werner, H., Schäfer, F., Henning, F., & Kärger, L. (2021). Material Modelling of Fabric Deformation in Forming Simulation of Fiber-Metal Laminates – A Review on Modelling Fabric Coupling Mechanisms. Paper presented at ESAFORM 2021. 24th International Conference on Material Forming, Liège, Belgique. doi: 10.25518/esaform21.2056