# Geometrical and mechanical draping of composite fabric

Abdelhakim Cherouat, Houman Borouchaki, Jean-Louis Billoët

► To cite this version:

Abdelhakim Cherouat, Houman Borouchaki, Jean-Louis Billoët. Geometrical and mechanical draping of composite fabric. Revue Européenne des Éléments Finis, 2005, 14 (6-7), pp.693-707. 10.3166/reef.14.693-707. hal-00020293

HAL Id: hal-00020293

https://hal.science/hal-00020293v1

Submitted on 27 Mar 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Geometrical and mechanical draping of composite fabric

Abdelhakim Cherouat* — Houman Borouchaki*
Jean-Louis Billoët**

* LASMIS - CNRS FRE 2719, Université de Technologie de Troyes
BP 2060, F-10010 Troyes
abdelhakim.cherouat, houman.borouchaki@utt.fr

** LMSP – CNRS UMR 8106, ENSAM-ESEM et Ecole polytechnique
de l'Université d'Orléans, 151 Bd de l'hôpital, F-75013 Paris
jean-louis.billoet@univ-orleans.fr

**ABSTRACT.** This paper presents an optimization based method for simulation of forming processes of woven fabric reinforced composites. Two approaches are proposed for the simulation of the forming of woven fabric: geometrical and mechanical. The geometrical approach is based on a fishnet model. It is well adapted to pre-dimensioning fabrics and to give a suitable quantification of the resulting flat patterns. The mechanical approach is based on a mesostructural model. It allows us to take into account the mechanical properties of pre-preg fabric and the various dominating mode of deformation of woven fabrics during the forming process. Some numerical simulations of forming process are proposed and compared with the experimental results in order to demonstrate the efficiency of our approaches.

**RÉSUMÉ.** Ce papier traite le problème d'optimisation de la mise en forme des tissus composites pré-imprégnés. Deux approches complémentaires sont proposées pour modéliser le drapage de ces matériaux : une approche géométrique et une approche mécanique. La première est basée sur le modèle du filet et s'inscrit dans une démarche de pré-dimensionnement, la seconde est basée sur une approche mesostructurale et s'inscrit dans une perspective d'optimisation de la conception et la fabrication de pièces composites.

**KEYWORDS:** Pre-preg woven fabric, fibres and resin, forming process, MOSAIC algorithm, meso-structural formulation, truss and membrane finite elements.

**MOTS-CLÉS :** pré-imprégnés, fibres et résine, mise en forme, algorithme MOSAIC, formulation mésostructurale, éléments finis de barre et de membrane.
# 1. Introduction

Composite materials with glass, carbon or aramid fibres and polymeric matrices are known to have high specific strengths and stiffness and, in combination with automatic manufacturing processes, make it possible to fabricate composite structures with high level of weight and cost efficiency. As known, the substitution of metal alloys by composite materials, in general, reduces structures mass by 20-30%. This increasing is due also to the numerous variety of semi-products (roving, fabrics, knitted fabrics, braids pre-impregnated or not) permitting the elaboration of structures. Fabrication processes, also, have undergone substantial evolution in recent years. Although the traditional hand lay-up process will remain the process of choice for some applications, new developments in resin transfer moulding, sheet moulding compound, low temperature curing prepregs and low pressure moulding compounds are reached new heights of sophistication, and are now being exploited in high technology areas such as aerospace industry (Binétruy, 2004).

The choice of manufacturing process depends on the type of matrix and fibres, the temperature required to form the part and the cost effectiveness of the process. Often, the manufacturing process is the initial consideration in the design of composite parts. Each manufacturing process imposes particular limitations on the structural design. It is necessary to develop specific algorithms taking into account not only the appropriated mechanics of the semi-products but also the associated process. Different levels of modelling intervene in the design process: architecture design level, pre-dimensioning level, mechanical level by computational software and optimisation level (Kawabata et al., 1973, Cherouat et al., 1995, Lim et al., 1999, Rozant et al., 2000, Boisse, 2004) and (Liu et al., 2005).

Most of these levels are integrated and take into account specificities of manufacturing processes. In the case of prepreg composite fabrics, these materials give flexibility in the design manufacture and play an important role in increasing the use of advanced composites. But for the manufacturing of non-developable pre-impregnated fabric composite parts, a new problem intervenes in the design chain resulting of the number of parameters influencing the global behaviour of composite forming process. The ability to define, in advance, the ply shapes and material orientation allowed the engineers to optimise the composite structural properties of the composite products for maximum strength, maximum material utilisation and maximum lay-up efficiency.

Hand lay-up is the most common method of producing thin composite parts (Rozant et al., 2000) and (Binétruy, 2004). The primary methods of automation in hand lay-up relate to computer software. Software is used to generate flat patterns from the layer surface and the ply boundary and to find the most efficient nest of cut plies to minimize the scrap. The formulation of the more efficient numerical models for the simulation of the forming composite processes must delays in manufacturing of complex parts and an optimisation of costs in an integrated design approach. The
composite manufacturing process involves large displacements and rotations and large shears of weft and warp fibres, which can have a significant effect on the processing and structural properties of the finished product. We propose two different numerical approaches to simulate the deformation of prepreg woven fabrics by shaping process: geometrical approach and mechanical approach (Van Der Ween, 1999, Cherouat et al., 2001, Borouchaki et al., 2003, Boisse, 2004) and Harrison et al., 2004). The numerical simulation is an efficient means of evaluating factors related to manufacturing processes and an efficient help to design pre-forming sequence for the manufacturing of fabric reinforced composites. It is possible to detect main problem occurring during the shaping deformation and to obtain good quantitative information on the forming process.

The geometrical approach is well adapted to pre-dimensioning level. It is based on geometrical aspects of the warping. Our method is based on a modified "MOSAIC" algorithm, which is suitable to generate a regular quad mesh representing the lay-up of the curved surfaces (giving the exact fibre orientations). The method is implemented in the GeomDrap software (Borouchaki et al., 1999) which is now integrated in the ESI-Pam software. This software provides a fibre quality chart (showing the fibre distortions, the rate of falling and the rate of draped surface) to predict difficult impregnated regions. It can be used to optimise the draping process (with respect to the above quality measure) by improving the lay-up directions or the marker data location. The lay-up of complex curved surfaces can be made in a few seconds.

Our mechanical approach is based on a mesostructural level. It permits to take into account the various dominating mechanics in the physics of the mechanical transformation of prepreg fabrics during the shaping process, namely large angular variations of yarns, viscoelasticity of resin and evolution of possible damages in yarns. The different advantages of this modelling is first to give the ability to obtain a very good material orientation of yarns and consequently to introduce good data in the pre-processing of the calculation of the final piece after polymerisation, secondly to give the mechanical limits of the fabric during the forming process and finally to limit the falls of the fabric by on better definition of the flat form (Cherouat et al., 2001). Different numerical simulations and experiments of shaping process have been carried out in order to validate the proposed computational formulation. The various forming parameters examined have included the initial shape of fabric, the fibre orientations and the viscosity of resin.

## 2. Discrete geometrical approach

First, we recall the mathematical formulation of the geometrical draping and then we propose an algorithm to solve the draping problem. Let denote by $\Sigma$ the surface to drape. We assume that a geometrical mesh $T_\Sigma$ of surface $T_\Sigma$ is known. Let F be the woven fabric modelled by two families of mutually orthogonal and inextensible
fibre described by the coordinates x = (ξ, η). These families constitute regular quad
mesh T<sub>F</sub> of the fabric F (see Figure 1). The problem of draping of woven fabric F
onto the surface Σ consists of associating each node of T<sub>F</sub> with a point of the
surface mesh T<sub>Σ</sub> such that the lengths of the edge of the corresponding mesh α on
the surface are preserved. This problem presents an infinity of solutions depending
on the first node of T<sub>F</sub><sup>Σ</sup> (called the point of impact) associated with a node of T<sub>F</sub>
and on the direction of the fabric co-ordinate axes at this node. Thus, to ensure a
unique solution, we suppose that the points of impact on the surface as well as the
fabric orientation at this point are known. The draping process is given by the
following scheme (Borouchaki *et al.*, 2003):

1. associate a node x₀ = (ξ₀, η₀) of T_F with the specified point of impact on the surface,

2. compute step by step the nodes of T<sub>F</sub><sup>Σ</sup>, classified as α-nodes, from the point of impact, associated with nodes (ξ, η<sub>0</sub>) of T<sub>F</sub>,

3. compute step by step the nodes of T<sub>F</sub><sup>Σ</sup>, classified also as α-nodes, from the point of impact, associated with nodes (ξ<sub>0</sub>, η) of T<sub>F</sub>,

4. compute cell by cell all the other nodes of T<sub>F</sub><sup>Σ</sup>, classified as β-nodes, from
x<sub>0</sub> and the nodes associated with nodes (ξ, η<sub>0</sub>) and (ξ<sub>0</sub>, η) of T<sub>F</sub>.

The nodes of T<sub>F</sub><sup>Σ</sup> associated with nodes (ξ, η<sub>0</sub>) and (ξ<sub>0</sub>, η) of T<sub>F</sub> (the α-nodes) are putted on the surface along the geodesic lines emanating from the point of impact. Regarding the β-nodes, various algorithms are proposed (Bergsma *et al.*, 1988) and (Van Der Ween, 1999). Most of them use an analytical expression of the surface and formulate the draping problem in terms of non-linear partial differential equations. Other algorithms are also proposed to simplify these equations by using a discrete approximation of the surface by flat triangular face (i.e. a mesh of the surface).

Based on this latter approach we propose a new algorithm. The geodesic lines on
the surface are approximated by the polylines plotted on the surface using linear
orthogonal transformation in $\mathbb{R}^3$ to setting flat the surface locally (these polylines
become a straight line segment after these transformations). This allows us to
determine the α-nodes. Let P be an α-node of T<sub>F</sub><sup>Σ</sup>, and K be the element of T<sub>Σ</sub>
containing P. To determine the geodesic line from P in a direction $\vec{u}$, we compute
the projection $\vec{u}_K$ of $\vec{u}$ on element K and we proceed as follow:

- go through K in the direction $\vec{u}_K$: we determine element K' adjacent to K intersecting direction $\vec{u}_K$,

- go from K to K': we compute the new direction $\vec{u}_{K'}$ following $\vec{u}_K$ such that
the angles of directions $\bar{u}_K$ and $\bar{u}_{K'}$ with the common edge of K and K' are the same (cf. Figure 2),

- replace K by K', $\bar{u}_K$ by $\bar{u}_{K'}$, and iterate.

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Geometrical draping of pre-preg woven fabric</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Computing Geodesic lines</figcaption></figure>

The $\beta$-nodes are computed by solving an optimisation problem corresponding to determine a vertex of an equilateral quad plotted on the surface from the data of the three other vertices. This optimisation problem formulates the direction of the geodesic lines emanating from the searched vertex. Suppose that the specified length of the edges of mesh $T_F$ is 1. Let P and Q be two nodes of $T_F^\Sigma$. To determine the geodesic segments emanating from P and Q, intersecting mutually with a length 1, the following algorithm is applied iteratively (directions $\bar{u}$ and $\bar{v}$ are appropriately initialized to obtain the desired solution between two possible solutions):
- compute the extremity $P'$ of the geodesic segment of length 1 emanating from $P$ in direction $\bar{u}$,

- compute the extremity $Q'$ of the geodesic segment of length 1 emanating from $Q$ in direction $\bar{v}$,

- if $P'$ is different from $Q'$, change $\bar{u}$ and $\bar{v}$ towards the mid point of segment $P'Q'$ and iterate (cf. Figure 3).

<figure><img src="image_4.png" /><figcaption>**Figure 3.** *Computing *β*-nodes*</figcaption></figure>

Nodes of $T_F^\Sigma$ associated with nodes $(\xi, \eta)$ where $\xi \neq \xi_0$ and $\eta \neq \eta_0$ complete the mesh elements for which three nodes associated to $(\xi^-, \eta^-)$, $(\xi^-, \eta)$ and $(\xi, \eta^-)$ are already placed on the surface mesh $T_\Sigma$ (node of $T_F^\Sigma$ associated with $(\xi^-, \eta)$ (resp. $(\xi, \eta^−)$) represents the node which precedes the node associated with $(\xi, \eta)$ in direction $\xi$ (resp. $\eta$) and node of $T_F^\Sigma$ associated with $(\xi^-, \eta^−)$ is the node which precedes the node associated with $(\xi, \eta)$ in the two directions $\xi$ and $\eta$). In the case where the mesh element touches the surface boundary, one of the nodes associated with nodes $(\xi^-, \eta)$ and $(\xi, \eta^−)$ can be undefined and the proposed algorithm to determine node associated with $(\xi, \eta)$ can not be applied. In this case, we identify the node of $T_F^\Sigma$ associated with one of the node $(\xi^-, \eta)$ or $(\xi, \eta^−)$ which is well defined (i.e. placed on the surface mesh) and we determine the node associated with $(\xi, \eta)$ on the geodesic line emanating from the well defined node associated with $(\xi^-, \eta)$ or $(\xi, \eta^−)$ in direction $\xi$ or $\eta$. Let us suppose that the node associated with $(\xi^-, \eta)$ is well defined and that the node associated with $(\xi, \eta^−)$ couldn't be placed on the surface mesh. Node associated with $(\xi, \eta)$ is then placed on the geodesic line emanating from node associated with $(\xi^-, \eta)$ in direction $\xi^{--}\xi^-$ where $\xi^{--}$ is such that node associated with $(\xi^{--}, \eta)$ precedes node associated with $(\xi^-, \eta)$ in direction $\xi$ (cf. Figure 4).
<figure><img src="image_5.png" /><figcaption>**Figure 4.** *Computing boundary $\beta$-nodes*</figcaption></figure>

**3. Continuum mechanical approach**

The deformation modes of prepreg woven fabrics during the forming process are different than the classical continuous sheet metal (Cherouat *et al.*, 1994), (Boisse, 2004) and (Liu *et al.*, 2005). A number of mechanisms are available including shear deformation between warp and weft yarns, yarn straightening, relative slip deformation and yarn buckling. In the deep-drawing and the draping processes of prepreg woven fabrics, the evolution of two straight line grille draw alternatively on warp and weft yarns, shows that the fabric can be considered as a continuous 3D surface domain. During the deformation, the draw lines become curved but remain continuous (see Figure 5). The absence of inter-yarn sliding (ensured by the weaving fabric, the friction fiber/fiber and the friction fiber/resin) can be observed over the main areas of the fabric (i.e. far enough from the free edges of the fabric). The buckling mechanism will appear when certain zones are locally subjected to compressive forces along yarns axis. For the prepreg fabrics based on high modulus yarns, the compressive as well as bending stiffness are negligible compared to the in-plane membrane stiffness. The function of reinforcement is to cover complex shapes by angular variations between warp and weft fibers without excessive extension of yarns. The function of resin is to avoid inter-fiber sliding and transmit the effort to fibers by shearing strain. The consequence is that the shearing strain energy is small for dry woven fabric but not for prepreg woven fabric.

<figure><img src="image_6.png" /><figcaption>**Figure 5.** *Deformation of woven fabric: straight lines before and after deformation*</figcaption></figure>
The warp and weft yarns are assimilated as a truss which connecting points are hinged and the resin, considered as a membrane, is coupled kinematically to the prepreg fabric at these connecting points. For each connecting point $\vec{X}^f$ of warp and weft yarns, in the reference configuration $C_0$, is associated a material position space of a resin $\vec{X}^m$ (see Figure 6). At the connecting points (I,J,K,L) we have $\vec{X}^f = \vec{X}^m = \vec{X}$ before deformation. The current space position of these points after deformation (i,j,k,l) is obtained by

$$ \left\{ \begin{array}{ll} d\vec{x}^f = \mathbf{F}^f(\vec{X}, t)d\vec{X}^f & \text{fibres} \\ d\vec{x}^m = \mathbf{F}^m(\vec{X}, t)d\vec{X}^m & \text{resin} \end{array} \right. \qquad [1] $$

where $\mathbf{F}^f$ and $\mathbf{F}^m$ are the deformation gradient tensor of fibre and resin respectively. The relationship of the no sliding inter-fibre can write at each connecting point as $\vec{x} = \vec{x}(\vec{X}, t)$ with $x^f = x^m = x$. For these points we have $\mathbf{F}^f = \mathbf{F}^m$ but for all other point $\mathbf{F}^f \neq \mathbf{F}^m$. It is possible to decompose the deformation gradient tensor in terms of the rigid rotation tensor followed by a stretch. By the spectral theorem, the stretch component (scalar for the fiber $\lambda_L^f$ and tensor for the resin $\mathbf{U}^m$) is defined in the reference configuration $C_0$ as :

$$ \left\{ \begin{array}{ll} \lambda_L^f = \sqrt{N_L^f \mathbf{F}^{fT} \mathbf{F}^f N_L^f} & \text{fibres} \\ \mathbf{U}^m = \sqrt{\mathbf{F}^{mT} \mathbf{F}^m} & \text{resin} \end{array} \right. \qquad [2] $$

Using the above assumptions, the mechanical behaviour of prepreg woven fabric depends on the deformation of yarns and the deformation of resin. The problem of the integration of strain rate tensors is a central one in large deformations. The rate equations for finite strains use objective derivatives (Gilormini et al., 1993). The frame associated with Green-Naghdi's derivative is defined, at the material point considered, by the rigid body rotation of the triplet orthogonal material directions ($N_i^f \otimes N_i^f$ for the fibre and $e_{0i} \otimes e_{0i}$ for the resin). The stretching tensors are written in the material frames.

$$ \left\{ \begin{array}{ll} \overline{\mathbf{D}}_L^f = \left( \frac{\lambda_L^f}{\lambda_L^m} \right) (N_L^f \otimes N_L^m) & \text{fibres} \\ \overline{\mathbf{D}}^m = \mathbf{U}^m \mathbf{U}^{m-1} (e_{0i} \otimes e_{0i}) & \text{resin} \end{array} \right. \qquad [3] $$

The longitudinal fibres stress $\overline{\sigma}_L^f$, depending on the fibre stretching elongation $\overline{D}_L^f$, and the resin membrane stress tensor $\overline{\sigma}^m$ can be written at each time as:
$$ \left\{ \begin{array}{l} \dot{\bar{\sigma}}_L^f = E_L^f (\lambda_L^f) \bar{D}_L^f \quad \text{fibres} \\ \dot{\bar{\sigma}}^m = C^m(\tau) \bar{D}^m \quad \text{resin} \end{array} \right. \qquad [4] $$

Because of the weaving, the behaviour in tension of the fibres is elastic non-linear and is related to elastic modulus $\bar{E}_f$ and yarns undulation factor $\varepsilon_{sh}$. In compression, the behaviour of fibres is supposed negligible. The viscoelastic behaviour of resin is formulates in the time domain by the hereditary integral and using the relaxation time $\tau_k$ and the fourth order relaxation tensor, which are material parameters $C_{ij}^{m\kappa}$:

$$ \left\{ \begin{array}{l} E_L^f (\lambda_L^f) = \bar{E}_f \left( 1 + \exp\left(\frac{-\dot{\lambda}_L^f}{\varepsilon_{sh}}\right) \right) \quad \text{fibres} \\ C_{ij}^m (t) = C_{ij}^{m\infty} + \sum_{l=1}^k C_{ij}^{m\kappa} e^{-t/\tau_k} \quad \text{resin} \end{array} \right. \qquad [5] $$

The governing equilibrium equations are solved as a dynamic problem using explicit integration. This approach has proven to be, in particular, suitable to highly non-linear geometric and material problems, particularly where a large amount of contact between different structural parts occurs. The analyses are carried out using the FEA computer code Abaqus/Explicit (Abaqus, 2002). The detailed formulations are described elsewhere (Cherouat et al., 1999).

<figure><img src="image_7.png" /><figcaption>Figure 6. Mechanical deformation of prepreg woven fabric</figcaption></figure>

## 4. Numerical examples

Four draping simulation examples are given. These simulations are performed using the geometrical analysis computer code GeomDrap (Borouchaki et al., 1999) and the FEA computer code Abaqus/explicit. For each example, we assume that a
mesh of the piece to drape is given. The first and fourth examples show the influence of the fibre orientations draping in the draping process. The first one uses the pure geometrical approach while the fourth concerns both the mechanical approach and experimental results. The last example concerns the 3D deep-drawing of complex part using mechanical approach and experimental results.

The first example concerns the discrete geometrical draping of a base plate piece. The centroid of this piece is chosen as the point of impact for which two different fibre orientations: (0°/90°) and (-45°/45°) are specified. Figure 7 shows the resulting 3D surface lay-up for the (0°/90°) fibre orientation and the 2D corresponding flat patterns. Likewise, Figure 8 shows the draping results for the (-45°/45°) fibre orientation. One can notice that, in the considered cases, the surface of the piece is draped globally. However, in the second case, a smaller area of the flat fabric is used. This result shows the importance of the fabric orientation in the draping process.

<figure><img src="image_8.png" /><figcaption>**Figure 7.** Geometrical draping of mechanical piece</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 8.** 2D Flat patterns of the composite parts</figcaption></figure>

The second example is also a discrete geometrical draping. The considered part is a car hood. The centroid of the part is chosen as the point of impact from which the (0°/90°) and (-45°/45°) fibre orientations are specified. Figure 9 shows the
resulting 3D draping without applying the particular procedure for boundary β -node. We can note that the surface part is not totally covered. Figure 10 shows the 3D draping of the part using the boundary β -node treatment. All of the surface part is now draped. Figure 11 presents shaded contours interpolated from the map of the fiber distorsions of (0°/90°) and (-45°/45°) fiber orientations. The fiber distorsions for both (0°/90°) and (-45°/45°) draping are very small.

<figure><img src="image_10.png" /><figcaption>**Figure 9.** Geometrical draping at (0°/90°) and (-45°/45°)</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 10.** Geometrical whole draping at (0°/90°) and (-45°/45°)</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 11.** Fibre distorsions map at (0°/90°) and (-45°/45°)</figcaption></figure>
The third example is the mechanical draping of a prepreg fabric (carbon fibre) on a hemispherical shape. The resin is modelled with linear membrane finite elements (triangular 3 nodes M3D3 and quadrilateral 4 nodes M3D4 of Abaqus elements), the warp and weft fibres are initially modelled with two orthogonal linear truss finite elements (T2D2) and the hemisphere tool is modelled by rigid surfaces (triangular 3 nodes R3D3 and quadrilateral 4 nodes R3D4). The behaviour of the resin is isotropic viscoelastic and the behaviour of the fibres is elastic non-linear. As for the first example, we consider the two fibre orientations (0°/90°) and (-45°/45°). Figure 12 shows the experimental draping with respect to these orientations. Notice that the initial directions of the fibre influence the final resulting shape. Figure 13 illustrates the numerical final shapes (for each orientation) obtained by mechanical approach. The numerical results are in good agreement with the experimental one.

<figure><img src="image_13.png" /><figcaption>**Figure 12.** *Experimental draping at (0°/90°) and (-45°/45°)*</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 13.** *Numerical draping at (0°/90°) and (-45°/45°)*</figcaption></figure>

The last example shows the 3D deep-drawing of aramid prepreg fabric with conical tools. Figures 14a, 14b, 14c and 14d report the experimental obtained shapes with respect to (-45°/45°) fibre orientation for different punch displacements (u = 0 mm, 20 mm, 60 mm and 100 mm, respectively). Figures 14e and 14f show the experimental shapes with respect to (0°/90°) and (-45°/45°) fiber orientations at 100 mm of punch displacement. Notice that, the initial directions of the fiber influence considerably the final resulting shape. Figure 15 presents shaded contours interpolated from the map of the fiber distortions of (0°/90°) and (-45°/45°) fiber
orientations. The fiber distorsion is represented by the orientation of the deformed truss elements (warp and weft). The angular fibre distortion exceeds 38° along the diagonal axis for (0°/90°) fibre orientations (Figure 16a) and along the diagonal axis for (-45°/45°) fibre orientations (Figure 16b).

<figure><img src="image_15.png" /><figcaption>**Figure 14.** *Experimental deep-drawing at (0°/90°) and (-45°/45°)*</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 15.** *Fibre distortions map for 100mm of punch displacement*</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 16.** Angular distorsions along the fabric lines (a) 0°/90° and (b) (-45°/45°)</figcaption></figure>

**5. Conclusions**

An efficient numerical approach (geometrical and mechanical) has been presented to simulate accurately the draping of composite fabric. This numerical methodology is shown to be very helpful for composite forming process.

**6. References**

Abaqus theory, User's Manual, 2002.

Bergsma O.K., Huisman J., “Deep Drawing of fabric reinforced thermoplastic”, 2nd *Inter. Conf. Comp. Aided Design in Composite Material Technology*, 1988, pp. 323-333.

Binétruy C., “Physique du moulage des composites avancés: aspects théoriques”, *Techniques de l'ingénieur*, AM 3718, 2004, pp. 1-16.

Boisse P., “Mise en forme des renforts fibreux de composites”, *Techniques de l'ingénieur*, AM 3734, 2004, pp. 1-10.

Borouchaki H., Cherouat A., Billoët J.L., *GeomDrap New Computer Aided Design and Manufacturing for Advanced Textile Composites*, Version 1, 1999.

Borouchaki H., Cherouat A., “Drapage géométrique des composites”, C.R. Acad. Sci. Paris, Série II B, *Mécanique des solides et des structures*, Vol. 331, 2003, pp. 437-442.

Cherouat A., Gelin J.C., Boisse P., Sabhi H., “Modélisation de l'emboutissage des tissus de fibres de verre par la méthode des éléments finis”, *Revue Européenne des Eléments Finis*, Vol. 4, 1995, pp. 159-182.

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. Mat. Proc. Technology*, Vol. 118, 2001, pp. 460-471.

Gilormini P., Roudier P., Abaqus and Finite Strain, Rapport interne No. 140, 1993, Cachan.
Harrison P., Clifford M.J., Long A.C., Rudd C.D., "A constituent-based predictive approach to modelling the rheology of viscous textile", *Composites Part A*, Vol. 35, 2004, pp. 915-931.

Kawabata S., Niwa M., Kawai H., "The Finite deformation theory of plain-weave fabrics", *J. Text. Inst*, Parts I, II and III, Vol. 64, 1973, pp. 21-83.

Lim T.C., Ramakrishna S., Shang H.M., "Optimization of the formability of knitted fabric composite sheet by means of combined deep drawing and stretch forming", *J. of Materials Processing Technology*, Vol. 89-90, 1999, pp. 99-103.

Liu L., Chen J., Li X., Sherwood J., "Two-dimensional macro-mechanics shear models of woven fabrics", *Composites Part A*, Vol. 36, 2005, pp. 105-114.

Rozant O., Bourban P.E., Manson J.A.E., "Drapability of dry textile fabrics for stampable thermoplastic performs", *Composites Part A*, Vol. 31, 2000, pp. 1167-1177.

Van Der Ween F., "Algorithms for draping fabrics on doubly curved surfaces", *Int. J. Num. Meth. Engng*, 31, 1999, pp. 1415-1426.