Hindawi Publishing Corporation Advances in Materials Science and Engineering Volume 2013, Article ID 709495, 18 pages http://dx.doi.org/10.1155/2013/709495

Research Article

Numerical Tools for Composite Woven Fabric Preforming

Abel Cherouat and Houman Bourouchaki

GAMMA3-INRIA Project Team, Charles Delaunay Institute, University of Technology of Troyes, 12 Rue Marie-Curie, BP 2060, 10010 Troyes, France

Correspondence should be addressed to Abel Cherouat; abel.cherouat@utt.fr Received 13 September 2012; Accepted 12 November 2012 Academic Editor: Philip Harrison

Copyright © 2013 A. Cherouat and H. Bourouchaki. His is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

An important step in the manufacturing processes of thin composite components is the layingup of the reinforcement onto the mould surface. "Hie prediction of the angular distortion of the reinforcement during draping and the changes in fibre orientation are essential for the understanding of the manufacture process and the evaluation of the mechanical properties of the composite structures. His paper presents an optimization-based method for the simulation of the forming processes of woven fabric reinforced composites. Two different approaches are proposed for the simulation of the draping of woven fabric onto complex geometries: geometrical and mechanical approaches. He geometrical approach is based on a fishnet model. It is well adapted to predimensioning fabrics and to give a suitable quantification of the resulting flat patterns. He mechanical approach is based on a mesostructural model. It allows us to take into account the mechanical properties of fibres and resin and the various dominating mode of deformation of woven fabrics during the forming process. Some numerical simulations of the forming process are proposed and compared with the experimental results in order to demonstrate the efficiency of our approaches.

1. Introduction

Composite reinforced by woven fabric is known to have high specific stiffness and, in combination with automatic manufacturing processes, makes it possible to fabricate complex components (aircraft, boat, automotive, and military) with high level of weight and cost efficiency due to their excellent mechanical properties. One of the most attractive properties of the fabric reinforcements is that they are easily handled and automatically processed, which can lower the production cost. As known, the substitution of metal alloys by composite materials, in general, reduces structural mass by 20-30%. He mass increase is also due to the numerous variety of semiproducts (roving, fabrics, knitted fabrics, and braids preimpregnated or not) permitting the development of new structures. Fabrication processes, also, have undergone substantial evolution in recent years. Although the traditional layup process will remain the process of choice for some applications, new developments in resin transfer molding (RTM) or sheet molding compound (SMC), low temperature curing prepregs, and low pressure molding compounds

have matured significantly are reached, and are now being exploited in high technology areas such as aerospace industry. For example, by using such composites, the automotive industry can realize improved fuel economy through vehicle-weight reduction by replacing the currently used steel and aluminum parts with thermoplastic composites with the added benefit of a corrosion-resistant material. He choice of the manufacturing process depends on the type of matrix and fabric, the temperature required to form the part, and the cost effectiveness of the process. In particular, thermoforming is a promising manufacturing process for producing highvolume low-cost composite parts using commingled fibre glass/polypropylene woven fabrics [1].

He simulation of the manufacturing of a textile-reinforced composite part with a liquid composite moldinglike process, which involves draping (or deep drawing) and impregnation of the preform, includes several stages (see Figure 1). First, a mould is designed with CAD software, and the CAD model is meshed. Hen, a draping (or deep drawing) simulation tool is used to compute the deformations of the textile layers inside the mould. As a result, for every element

CAD design

Mould meshing

Geometrical or mechanical draping

Flat pattern of the fabric

Mechanical simulation after polymerization Fibre orientations

Figure 1: Illustration of the manufacturing of a textile reinforced composite part.

of the mesh, textile parameters like the shear angle and the thickness of the layer are available. With these parameters given, the local (mesoscale) properties for every element are determined (preprocessing) and given as input for the macroscopic structure after resin polymerization simulation. tte result of the macrosimulation is then postprocessed to optimize the mechanical properties of the composite structure. tte numerical simulation of composite forming is an efficient means of evaluating the factors related to the manufacturing processes and an efficient help to design preforming sequence for the manufacturing of fabric-reinforced composites. It is possible to detect the main problem occurring during the shaping deformation and to obtain good quantitative information on the forming process [2-5].

Different levels of modeling intervene in the simulation of woven reinforcement forming: (1) architecture design level, (2) preliminary design level, (3) mechanical level by computational software, and (4) optimization level [4-6]. Most of these levels are integrated and take into account specific constraints of the manufacturing processes. tte particular form of composite fabrication (prepreg) begins with the preimpregnation of reinforcement materials with a resin. tte combining of these two materials occurs prior to the molding process and therefore enables a very accurate reinforcement to resin ratio to be achieved. Prepreg materials are used extensively in the aerospace or automotive industry due to their ability to maximize strength to weight ratios. Prepregs are pliable and therefore are able to be cut into various shapes or patterns prior to processing into the molded products. But for the manufacturing of nondevelopable composite part (part that cannot be made flat and unstretched), a new problem intervenes in the design chain resulting from the number of parameters influencing the global behavior of composite forming process. tte ability to define, in advance, the ply shapes and material orientation allowed the engineers to optimize the composite structural properties of

the composite products for maximum strength, maximum material utilization, and maximum layup efficiency [7-15].

tte composite manufacturing process involves large displacements and rotations and large shear of weft and warp fibres, which can have a significant effect on the processing and structural properties of the finished product. tte effective modeling of the forming properties of woven fabric polymer composite materials has been the focus of research for more than a decade. tte formulation of new and more efficient numerical models for the simulation of the shaping composite processes must allow for reduction in the delay in manufacturing of complex parts and an optimization of costs in an integrated design approach [15-20]. Several modeling approaches have been developed to account for the evolution of the orthotropic directions during high shearing, and these approaches include the geometrical and the finite element approaches.

In this context, simulation methods are needed to anticipate the performance of the final part but also to predict the reinforcement preforming and the resin injection. Early methods, based on a geometric approach of the fishnet type algorithm, emerged a few decades ago. tte geometrical approach based on the so-called fishnet algorithms is used to determine the deformed shape of draped fabrics. tte main assumptions are that the fibers are inextensible, the intersection points between warp and weft yarns are fixed during preforming, and the rotations between warp and weft yarns are free. ttis method, where the fabric is placed progressively from an initial line, provides a close enough resemblance to handmade draping [18-21]. ttese methods were originally developed for prepreg hand draping. ttey are very fast and fairly efficient in many prepreg draping cases. Nevertheless, this method has major drawbacks. ttey account neither for the mechanical behavior of the fabric nor for the static boundary conditions. ttis last point is very important in the case of forming with punch and die

(such as in the preforming of the RTM process). He loads on the tools, especially on the blank holder, influence the quality of the shaping operation and therefore need to be considered in simulations [21,22]. More recently finite element methods have been used for the draping simulation of composite reinforcements.

He alternative to the geometrical approach consists of a mechanical analysis of the fabric deformation under the boundary conditions prescribed by the forming process. His requires a specific model of the woven reinforcement and its mechanical behavior. He mechanical behavior of woven fabrics is complex due to the intricate interactions of the fibres. It is a multiscale problem. He macroscopic behavior is very much dependent on the interactions of fibres at the mesoscale (scale of the woven unit cell) and at the microscale (level of the fibres constituting yarns). Despite of the great amount of work in the field, there is no widely accepted model that accurately describes all the main aspects of fabric mechanical behavior. He main model families come from the multiscale nature of the textile. A first family of models is obtained by homogenizing the mechanical behavior of the underlying mesostructure and considering the fabric as an anisotropic continuum [9-15]. If these models can easily be integrated in standard finite element using conventional shell or membrane elements, then the identification of homogenized material parameters is difficult, especially because these parameters change when the fabric is strained and when, consequently, the directions and the geometry (crimp, transverse sections, etc.) of the fibres change. Some of these approaches will be described, especially a nonorthogonal constitutive model [23-31] and an anisotropic hypoelastic continuous behavior for fibrous material based on an objective derivative using the rotation of the fibre [24]. Conversely, some authors present finite element models to model the warp and the weft fibres behavior [28]. Each yarn or each fibre is modeled and is assumed to be a straight or a curved beam or truss. Sometimes they are modeled as 3D domains [28]. Springs are often used to model warp and weft yarn interactions. In the objective of fabric forming simulations, some authors extend the finite element modeling to the whole textile structure that is represented by a network of interwoven trusses or beams with different tensional and rotational springs. Accounting for the simplicity of each component, the whole textile structure deformation can be computed.

Nevertheless, the computational effort needed is relatively significant. At present, this method is restricted to simple geometry of the local yarn and relatively simple mechanical behavior. When a fine model of the fibrous yarns is used, the analysis can only consider a small part of the textile reinforcement such as a few woven or knitted cells. He semidiscrete approach is a compromise between the above continuous and discrete approaches [14]. A finite element method is associated to a mesoscopic analysis of the woven unit cell. Specific finite elements are defined and are made of a discrete number of woven unit cells. He mechanical behaviour of these woven cells is obtained by experimental analyses or from 3D FE computations of the woven cell. He nodal interior loads are deduced from this local

behavior and the corresponding strain energy in the element deformation.

For experimental technique, Kenkare and May-Plumlee [32] distinguished draping from paperiness (paper bending behaviour) by accounting for 3D deformation behavior using a drapemeter that is capable of distorting fabric samples in all three dimensions over a circular disc. Providing an accurate vertical projection of the draped sample by optical means allowed the determination of the drape coefficient (DC), which they defined as "the fraction of the area of the annular ring covered by the projection of the draped sample." Hence, a high drape coefficient reflects poor or little deformation.

Over several years, it was the main concern of drape researchers to accurately record the draped contour. From complex projection techniques [33] to the use of photovoltaic cells [34] and "cut-and-weigh" techniques, image analysis techniques [35] have definitely simplified investigations, providing higher accuracies. Most measuring methods produced similar results. However, image analysis techniques greatly reduced result dispersion as well as measuring time.

In this study, geometrical and finite element approaches are used to simulate the deformation of preforming of woven fabric. Hese approaches, while giving good results and being efficient in terms of computing time, are generally somewhat complex and sometimes very challenging to implement into commercially available FEA packages.

(i) He geometrical approach is well adapted to preliminary design level. It is based on geometrical aspects of the warping. It is based on a modified MOSAIC algorithm, which is suitable to generate a regular quad mesh representing the layup of the curved surfaces (giving the exact fibre orientations). He method is implemented in the GeomDrap software which is now integrated in the ESI-Pam software. His software provides a fibre quality chart (showing the fibre distortions, the rate of falling, and the rate of draped surface) to predict local folding due to that the overlapping of fibers in the shear exceeds limit value (up to 60° in some cases). It can be used to optimize the draping process (with respect to the above quality measure) by improving the layup directions or the marker data location. He layup of complex curved surfaces can be made in a few seconds [18, 22].

(ii) He mechanical approach is based on a mesostruc-tural description for finite deformations and geometrical nonlinearity. He unpolymerized resin has a viscous behavior, and the reinforced fibers are treated as either unidirectional or nonlinear elastic behavior. He unit cell of the mesoscopic model used here for a plain weave fabric consists of bicomponent finite elements. He tensile load is carried by the 1D elements that will capture the changes in the orthotropic directions during the shearing. He 2D element accounts only for the shearing resistance of the fabric and hence has no tensile stiffness. He bicomponent finite elements for modeling composite fabric behavior are based on 3D membrane finite elements representative

of resin behavior and truss finite elements representative of warp and weft fibers behavior. tte efficiency of the proposed model resides in the simplicity of its finite element discretization and the performance of its mechanical background [30, 31].

2. Geometrical Approach

tte draping of composite fabric using a mechanical approach requires the resolution of equilibrium PDE's problems by the finite element method (see Section 3). In general, in the case of complex surfaces, the boundary conditions are not well defined, and the contact between the surface and the fabric is difficult to manage. Furthermore, the resolution of such a problem can be too long in CPU time and is detrimental to the optimization stage of draping regarding the initial fibre directions. All of these facts lead us to consider rather a geometrical approach which is very fast and more robust allowing simultaneously to define the stratification sequences and the flat pattern for different plies and to predict difficult impregnated areas which involve manual operations like dart insertion or, on the contrary, the shortage of fabric. Based on technical criteria (mold surface covering, fabric drape covering, and fibre angular distortion), this approach can constitute the predimensioning or the preoptimization stage for the manufacturing of complex composite parts. tte geometrical approach is based, in general, on the fishnet method for which a fabric mesh element is subjected only to shear deformations. tte difficulty of such a method is the mapping of the fabric mesh element onto any surface. Within this context, several algorithms approximating the geometry of a fabric mesh element plotted onto the surface are proposed. In particular, the edges of the fabric mesh element are approximated by line segments representing a pure estimate in the highly curved area (see [19-22] for a synthesis).

In this paper, we propose a new geometrical algorithm which takes into account the true geometry of the fabric mesh element plotted onto the surface. Such a fabric mesh element is then defined by a curved quadrilateral whose edges are geodesic lines with the same length plotted onto the surface to drape. Given three vertices of the fabric mesh element on the surface, we propose an optimization algorithm to define the fourth vertex of the fabric mesh element. ttis algorithm allows us to drape the complex surface using an advancing front approach from the data of an initial start point between the fabric and the surface and the initial fibre directions at this point. In this section, geometrical approach based on the mould mesh discretization is presented. For the second formulation, we propose an algorithm of composite fabric draping without any approximation on the geometry of surface to be draped.

2.1. Geometrical Formulation. First, we present the mathematical formulation of the geometrical draping, and then we propose an algorithm scheme to solve the draping problem. Let Z denote the surface of the part to drape, and we assume that a geometrical mesh TZ of surface is known. Let O be the woven composite fabric modeled by two families (warp and

weft) of mutually orthogonal and inextensible fibre described by the local coordinates x = ttese families constitute

regular quadrilateral fabric mesh TF of the fabric O (Figure 2 gives example of draping steps of complex part). tte problem of geometrical draping of O onto the surface Z consists of calculating each node displacement of fabric mesh TF with a point of the surface mesh TZ such that the lengths of the edge of the corresponding mesh TF on the surface are preserved (not extensible). ttis problem presents infinity of solutions depending on

(i) starting point associated with a node of fabric TF,

(ii) initial warp and weft orientation a.

ttus, to ensure a unique solution, we suppose that the points of impact on the part surface as well as the fabric orientation are given. tte draping scheme is given by the following step [36]:

(1) associate a starting point (corresponding to the point of impact of the machine to drape) on the surface on geometrical part mesh xZ = (£0, (Figure 2(a)),

(2) compute numerically step by step the warp nodes of TF, classified as a-nodes, from the starting point, associated with nodes of TF (Figure 2(b)),

(3) compute numerically step by step the weft nodes of TF, classified also as a-nodes, from the starting point, associated with nodes of TF (Figure 2(c)),

(4) compute numerically cell by cell all the other nodes of TF, classified as p-nodes, from x0 and the nodes associated with nodes and (£0, of TF (Figure 2(c)).

tte nodes of TF associated with nodes and (£0, of TF and the a-nodes are located on the surface along the geodesic lines emanating from the point of impact. Regarding the nodes, various algorithms are proposed [20]. Most of them use an analytical expression of the surface and formulate the draping problem in terms of nonlinear partial differential equations. Other algorithms are also proposed to simplify these equations by using a finite element discretization of the surface by flat triangular face (i.e., a mesh of the surface). Based on this latter approach, we propose a new algorithm. tte p-nodes are computed by solving an optimization problem corresponding to determine a vertex of an equilateral quad plotted on the surface from the data of the three other vertices. ttis optimization problem formulates the direction of the geodesic lines emanating from the searched vertex.

Consequently, two problems arise:

(1) problem 1: determine the geodetic exit of a given point of surface according to a given orientation,

(2) problem 2: determine the geodetic exits of these points intersecting themselves mutually according to given two points of surface and lengths (this geodetic is given according to their orientations).

tte developed algorithm is implemented in GeomDrape tool [18]. ttis software provides a fibre quality chart

Geodesic lines

(a) Starting point and draping direction «-Nodes fi-Nodes

(b) Geodesic lines and fibres directions

(c) Computing the four nodes of fabric cell

(d) Flat pattern of the initial fabric

N^^^WffrVT

(e) Shear angles between fibres

(f) Stretching of fibres

Figure 2: Geometrical draping steps of woven composite fabric.

(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 layup directions or the marker data location. ^e layup of complex curved surfaces can be made in a few seconds. ^e use of draping tool should not only allow for a more detailed assessment of the draped fabric (Figure 2(c)), but also for an optimization in terms of total fabric shear deformation (Figure 2(e)), flat pattern of the initial fabric (Figure 2(d)) and fibre stretching

(Figure 2(f)), and the optimization of the best layup start point (Figure 2(a)).

2.2. Application of Geometrical Approach. Two draping simulation examples are given. ^ese simulations are performed using the geometrical analysis computer code GeomDrap [18]. For each example, we assume that a mesh of the piece to drape is given. ^e first example shows the influence of the fibre orientations draping in the draping process. ^e second example shows the efficiency of the proposed

method to simulate geometrically the draping composite part.

^e first example concerns the geometrical draping of a complex piece composed of two half hemispheres with a radius of 38.8 mm that are connected by a half cylinder with a length of 170 mm [23]. Figure 3 shows the piece as well as a mesh of the mechanical part. ^e centroid of the part is chosen as the impact point for which two different fibre orientations are specified. Figure 4 shows the resulting 3D surface layup for the 0°/90° fibre orientation (Figure 4(a)) and the corresponding 2D flat patterns (Figure 4(b)). Likewise, Figure 5 shows the draping results for the ±45° fibre orientation (Figure 5(a) 3D surface layup and Figure 5(b) the corresponding 2D flat pattern). ^e shears angles between fibres are given in Figure 6. One can notice that, in the considered cases, the surface of the piece is draped globally. However, to drape completely the proposed surface without excessive fibre distortion <50°, it is necessary to make cuts chisel along the direction 45° for 0°/90° fibre orientation (Figure 6(a)) and along the direction 0° ± 45° fibre orientation (Figure 6(b)). In case we can allow excessive distortions of the fibers without cutting chisel, the results given in Figure 7 show the distribution of fibers in both directions 0°/90° (Figure 7(a)) and ±45° (Figure 7(a)). We can see that the shear limit reached 52° in the case of 0°/90° fibre orientation and 89° in the case of ±45° fibre orientation.

To compare the experimentally obtained fibre angle result [23] with the geometrical predicted ones, two cross sections where most shearing occurs are examined. ^e amount of shear along the edge of these cross sections is used to compare with the draping approach. Figure 3 gives an overview of the two cross sections along the symmetrical L1 and along the median line L2 for different ply orientations. ^e line L2 corresponds with the examined cross section for a 0°/90° fabric draping orientation and the line L1 for a 45 ± 45° fabric draping orientation. From Figure 8(a), it can be concluded that for the 0°/90° orientation, the agreement between the experimental shear angles and the geometrical draping results is good. On the other side, the drape orientation ±45°, the results from the geometrical model do not agree at all with the experimental results [23] (see Figure 8(b)). ^e oversimplification of the fabric deformation in this geometrical model gives shear angles up to 89° in case of the 45° fabric orientation, which is impossible for fabric composite. ^e geometrical model is used with a cutoff shear angle based either on an experimentally determined locking angle or the maximum orientation that the designer is prepared to tolerate, when defining a cutoff angle of 38°, which equals the experimentally determined locking angle. From the above discussed results, it can be concluded that the geometrical draping method is not appropriate for predicting the fiber directions of a fabric sheet after stamp forming.

^e second example is the draping of a complex shape (car hood). ^e centroid of the part is chosen as the starting point from which the 0°/90° and ±45° fibre orientations are specified. Figure 9 shows the resulting 3D draping for the two orientations (Figure 9(a) for 0°/90° fabric orientation and Figure 9(b) for 45°/45° fabric orientation). We can note that all the part surface is completely draped. Figure 10

Figure 3: CAD of complex part.

presents shaded contours interpolated from the map of the fiber distortions of 0°/90° (Figure 10(a)) and ±45° fiber orientations (Figure 10(b)). ^e fiber distortions for both 0°/90° and ±45° draping are very small but the maximum shear angle localization is different.

3. Mechanical Finite Element Approach

^e mechanical behaviour of composite fabric forming is complex, and it is a multiscale problem due to the fibres behavior and the interactions of the fibres or yarns [3741]. ^e macroscopic behaviour is very much dependent on the interactions of fibres at the mesoscale (scale of the woven unit cell) and at the microscale (level of the fibres constituting yarns). Despite of the great amount of work in the field, there is no widely accepted model that accurately describes all the main aspects of fabric mechanical behaviour. ^e main model families come from the multiscale nature of the textile. During the forming process of woven fabric, the two main modes of deformation at the mesoscopic scale are the stretching of the fibres due fibres undulation and the in-plane shearing of the fabric resulting in a change of the angle between the warp and the left yarns. In the deep drawing or the draping of woven fabrics, the in-plane shear of fibres is the principal mode of deformation and is very different than the sheet metal [37, 40]. Figure 11 shows the evolution of two straight lines draw alternatively on warp and weft fibre directions during the forming deformation. ^ese lines become curved but remain continuous. ^e absence of interyarn sliding ensured by the fabric weaving, viscoelastic behavior of resin, and friction fiber/fiber and fiber/resin can be observed over the main areas of the fabric (i.e., far enough from the free edges of the fabric). Also, for the composite fabrics based on high modulus, the compressive as well as bending stiffness is negligible compared to the inplane membrane stiffness. ^e assumption is that each cross connection of straight warp and weft fibres before deformation remains cross connected during the deformation. ^e basic assumptions for the mechanical forming are that the woven fabric is considered as a continuous 3D material. ^e warp and weft fibres are represented by a truss to which connecting points are hinged, and the membrane resin

V -Q 20

Distance along L2 (mm)

♦ Experimental [23] - Geometrical draping

90 80 70 60 50 40 30 20 10 0

Distance along L1 (mm)

♦ Experimental [23] - Geometrical draping

Figure 8: Shear angles of optimal geometrical draping: (a) along the line L2 of 0°/90° and (b) along the line L1 of-45°/45°.

(a) (0°/90°)

(b) (-45°/45°)

Figure 9: 3D geometrical draping of hood car part: (a) 0°/90° fabric orientation and (b) -45°/45° fabric orientation

(a) (0°/90°) (b) (-45745°)

Figure 10: Isovalues of shear angles between fibres: (a) 0°/90° fabric orientation and (b) -45°/45° fabric orientation.

Figure 11: Fibres deformation during the forming process.

is coupled kinematically to the fabric at these connecting points.

3.1. Model Description. In the unit cell of the mesoscopic model used here for a plain weave fabric, each fibre and resin is modelled and is assumed to be a straight or a curved beam or truss. ^e tensile load is carried by the 3D truss or beam elements that will capture the changes in the orthotropic directions during the shearing. ^e 3D membrane or shell element accounts only for the shearing resistance of the fabric and hence has no tensile stiffness. ^e interaction warp-weft yarn and resin fibres is negligible. In this approach, the stress and strain of a continuous material are related to fibrous reinforcement using the constitutive relation in a nonorthogonal frame directed by the fibre directions. In this study, we consider two yarn directions, and we use them to define the nonorthogonal frame. ^is approach uses the Green-Naghdi frame. It is an orthonormal frame which is rotated by R, the rotation of the polar decomposition in which the local stress increment computations at finite strain are made [42-44].

Two reference frames have to be considered. ^e ej unit vectors define the local orthogonal reference frame that rotates with the continuum material, and the gj basis vectors form a nonorthogonal frame that follows the fibre direction. Here, g1, g2 correspond to warp and weft directions, respectively. In this case, for each connecting point X of warp and weft yarns is associated a material position space of a resin Xm. At the connecting points, we have Xf = Xm = X before deformation. ^e current position of these points is obtained by

dxf = Ff(X,t)dXf fibres,

dxm = Fm (X,t)dXm resin,

where Ff and Fm are the deformation gradient tensor of fibre and resin, respectively. ^e relationship of the no sliding interfibre can be written at each connecting point as x = x(X, t)/xf = xm = x. ^e gradient of transformation of the

preimpregnated woven fabric F and the pseudo gradient of

transformation of the fibre F are defined by the function

W/ = ® 9oj fibres « dxm (2)

FTj = ® % resin

where A£ is the longitudinal elongation of each fibre and gQi, and gi are, respectively, the fibre orientations in the initial C0 and the current Ct configurations and are the resin. ^e relative rotation of fibre can be associated to the rotation of the rigid body of the median line of the fibre R. We can note this assumption by the following kinematic relation:

R = &®0of. (3)

Using the previous assumptions, the mechanical deformation of composite fabric depends on the relative movement of fibres and the deformation of unpolymerized resin. Consider

V = fib^

Um = VFmTFm resin.

^e shaping problem imposes the use of incremental formulation in finite deformations. In finite deformation analysis, a careful distinction has to be made between the coordinate systems that can be chosen to describe the behaviour of the body. ^e rate constitutive equations for finite strain use objective derivatives [43]. ^e problem of the integration of strain rate tensors is a central one in large deformations. ^e rate of deformation tensor of woven fabric is obtained by

D/ = fibres,

Dm = 2 (-F-F-1 + Co,-) resin.

^e rate equations for finite strains use objective derivatives [44]. ^e approaches traditionally developed in finite element codes for anisotropic metal at large strains are based on Jaumann corotational formulation or the Green-Naghdi approach. In these models, a rotation is used both to define an objective derivative for the hypoelastic model law and to update the orthotropic frame. ^e rotations used in Green-Naghdi and Jaumann derivatives are average rotations of the material (polar rotation and corotational rotations, respectively). ^e frame associated with Green-Naghdi's derivative is defined, at the material point considered. ^e rotation R^ is used to update the initial constitutive axes of warps or weft fibres gQi ® to the current constitutive axes, and the rotation R is used to update the initial constitutive axes of resin % ® %. ^e stretching tensors in (5) are written in the rigid body rotation frames. ^e longitudinal component

—/ ■ f f —/ —/ is Dl = A^/A^, and the transversal components (DT,D3)

are obtained by the unidirectional behaviour of fibre as

= Df = -vltd{R ^ A^ = A{ = (A{)-Vir, and vLT is the Poison's ratio of fibre. Using Green-Naghdi's objective

tensor stress, the stress rate depending on the stretching

deformation d{, and the stress rate tensor of the membrane resin CTm, depending on the tensor deformation rate Dm and the elastic properties Cm, can be written at each time as

follows. He stress and strain relationship of the woven fabric composites are defined in the nonorthogonal material coordinate frame. Herefore, coordinate transformations of stress and strain into the orthogonal coordinate system should be considered in the nonorthogonal constitutive model as o = ct^ + äm as follows:

/ I£î(AowE{(A{)I

weft V ' A{ wef

E —m Ai r^m vmE

1- v 2 -L Vm 1

r^m vmE —m Du rm E

1 - v 2 ^ vm 1 -vm:

warP v ' \{ weft y A,

He constitutive law of fibres is nonlinear and is written in terms of longitudinal modulus of stretching e£(A{), and

the compressive stiffness (A{ < 0) of fiber is supposed negligible E^ = 0. Later is the function of elongation of warp and weft fibre (A{, A{), effective elastic modulus of fibre Ey, and undulation factor £sh; Em, vm are the membrane elastic properties. To determine the stress state in a membrane material at a given time, the deformation history must be considered. For linear viscoelastic materials, a superposition of hereditary integrals describes the time-dependent response. Let Gm(i) be the shear stress relaxation modulus

of the unpolymerized resin and Gm = Gm (i = œ) the limit value. He viscoelastic behaviour of unpolymerized resin is formulated in the time domain by the hereditary integral

and using the relaxation time and the shear modulus

relaxation, which are material parameters Gm . Hereditary integrals with Prony series kernels can be applied to model the shear behavior of the unpolymerized resin. He behavior of fibre and resin can be written as

^M1-^!: ))fibrs

Gm (i) = Gm (œ) + ^G™ exp ( —

3.2. Finite Element Formulation. Each material point is moving as in a continuum, ensured by the nonsliding of

fibers due to fabric weaving, and resin behavior. Herefore, a nodal approximation for the displacement can be used. He deformation of composite fabric is described within the frame of membrane assumptions. He energy of deformation n(u) is obtained by a summation of membrane strain energy of unpolymerized resin and elastic tensile strain energy of fibers

(ù) = fy, J

Sn(ù) = fio| äm:<5D ds +

fibres JL

: c5D{dZ

f i-ôùidr,

where denote the initial thickness of fabric, l/ the length of fiber, and Sq the initial effective cross section of fibre, and t is the external surface load applied along rCT of woven fabric. He effective cross section of the fibre Sq and the effective surface S™ of the membrane resin, that assumed no void between the fibre and the resin, were used and are calculated by using the fibre volume fraction Vy of the

woven fabric (Vfabric = VyS^ l/ + (1 - Vy)fo0S™). He global equilibrium of the fabric is obtained by minimizing the total potential energy n(u). He effect of spatial equilibrium of composite material on the actual configuration is established in terms of nonlinear equations: kinematic nonlinearity, material nonlinearity, and contact with friction nonlinearity. It is linearized for each load increment by an iterative Newton method. It should be emphasized that during the motion, nodes and elements are permanently attached to the material points with which they were initially associated.

Consequently, the subsequent motion is fully described in terms of the current nodal positions as

x=^Nk (£ rç) Xk + £Nk (Ç, ri) uk,

where u are the nodal displacements of each connecting point, are the standard shape functions (of mem-

brane of truss element), and (m) denotes the number of nodes.

tte discretization of rate deformation tensor can be obtained by introducing (5) into the definition of d{ and D™ and given in (8) to give

—f A{ 1 (dxT du\ rn , . ,,

= d^ J = [Bfibres] fibres'

—m lT( du df 9uT d^', r D + )R = [BresinK resin>

dÇ dx dÇ dx

_ [dX/T dX/ 2dX/T du du^ du (21)

""" di di + di di + d£ d£'

and [Bresin] and [Bfibres] are the geometric or strain-displacement matrix of membrane resin and truss fibres.

In a finite element approximation, the only independent variables in the equations of linearized virtual work are the displacements of the material points. Substituting the element coordinate and displacement interpolations into the equilibrium equations, for a given set of elements, we obtain that the state that forces acting on a fabric equals the mass times the acceleration of the body

+ X ({RU-ICt}) = {0}, (12)

where [Me] is the consistent composite mass matrix, and ({Rfnt} - {lO) is the so-called quasistatic equilibrium residual

m=prh0 \m[Ne]T[Ne] ds+YpfSf0 f [ivf[Ne] dl,

JS fibres JL

{Rfnt}-{R:xt} = ho[m[B?eSin]T{^mK} d5

- XS^[fË{[BfibreSf-K*} dl

fibres

[ [Afffl ds.

[Ne] and [Ne] are the matrix of the nodal interpolation functions of both associated elements, and h0, Sq are the initial thickness and surface of resin and fibres, respectively. tte index e refers to the eth element.

According to the different modes of deformation occurring in the prepreg fabric during the shaping process, bicomponent finite elements are developed to characterize the mechanical behavior of thin composite structures. tte bicomponent element is based on an association of 3D linear membrane finite elements (T3 and Q4) combined with a complementary truss linear finite elements. tte global stiffness of composite fabric is obtained by the summation of elementary stiffness matrix of warp fiber, elementary stiffness matrix of weft fiber, and elementary stiffness matrix of resin. ttese finite elements are complementary in the finite element discretization, isoparametric and use three dof per node, and use the same mechanical formulation in finite deformations (Green-Naghdi's approach). tte nonlinear constitutive equation of fibre behavior is implemented in the Abaqus/Explicit using VUMAT user's subroutine [1315]. tte governing equilibrium equation (12) is solved as a dynamic problem using explicit integration. ttis is achieved by using the central difference method to approximate the velocity and the acceleration in the next time step, using only information from the previous step, where all state variables are known. ttis approach has proven to be, in particular, suitable to highly nonlinear geometric and material problems, particularly where a large amount of contact between different structural parts occurs. In the present work, the dynamic explicit (DE) resolution procedure is used within the general purpose FE code Abaqus/Explicit. tte DE algorithm available in Abaqus/Explicit for solving the algebraic system works by using the lumped form of the mass matrix [26]. tte major disadvantage of the explicit scheme is that it is only stable for time steps small enough. ttis can cause, for instance, the energy balance to be changed. One simple and conservative criterion for a stable time increment is

1 stable _f '

where Le is the minimal length of membrane element, pr is

the resin density, and E^ is Young's modulus of fibre. tte interpretation of this requirement is that the time increment must be shorter than the time it takes for a propagating wave in the material to cross the shortest side of an element. ttus, the element size and the critical time increment are connected. tte smaller the elements are, the shorter the time increment must be. Because of the difference in wave speed for different materials, the critical element size is larger for a stiff material (e.g., steel) than for a softer material (e.g., polymer).

3.3. Application of Mechanical Finite Element

3.3.1. Uniaxial Test. Due to the importance of the composite fabric behaviour on material formability, tensile test of preim-pregnated fabric is proposed in order to study the influence of fibre orientation, fibre undulations due to fabric weaving and resin behaviour. At processing, the impregnated fabric is idealized as a viscous material subject to the kinematic constraints of incompressibility and inextensibility in the fibre direction. tte preimpregnated fabric tested in this study

(a) (b)

Figure 12: Tensile test of preimpregnated woven fabric.

Table 1: Mechanical properties of the pre-impregnated composite fabric.

0 10 20 30 40 50 60 70 80 90 100 Displacement (mm)

Figure 13: Load force versus displacement for different fibre orientations.

was a satin 5 with aramid woven fabric. He fabric was impregnated with epoxy resin using a hot-melt prepregging process (the mechanical properties are given in Table 1). A lower loading velocity will generate lower viscous forces at the intraply shearing of fabric. In this experiment, a displacement in the vertical direction is imposed at the moving extremity of the rectangular composite specimen with three layers of fibres (length = 150 mm, width = 30 mm, and thickness = 2 mm). He uniaxial tensile test (see Figure 12) is carried out for different orientations of fibres with the loading direction (0, 15, 30, and 45°). He experimental effort imposed by the tensile machine is compared to the numerical values for different fibre orientations in Figure 13. In this figure, we can show the good correlation between the model and the experimental results. He agreement between predicted and

Ef (MPa) £sh p (g/cm3)

130000 0.005 1,45

Time (s) k Shear modulus of resin Gm

0.01 0.02332

0.1 0.023332

1 0.083509

10 0.11723

100 0.14423

1000 0.178

experimental values is good and proves the validity of the proposed model of preimpregnated woven fabric behaviour. He numerical model described previously clearly shows the strong nonlinearity of this behaviour law. It takes into account the mechanical characteristics of a viscoelastic resin, the anisotropic behaviour of fabric, and the geometrical nonlin-earities due to the high deformability of fibres (straightening and relative rotation).

3.3.2. Deep Drawing of Dry Woven Fabric. He second example is the 3D deep drawing of dry glass woven fabric by hemispherical punch by experimental and mechanical approaches. He initial shape of the glass fibre fabric is a square (360 x 360 mm). Its edges are free, but a pressure equal to 2 MPa is applied on the binder and the friction coefficient between the glass fabric and the steel tools is 0.27. He forming simulation has been performed in Abaqus/Explicit

\V v >[ <sf№ J

■v>< :<• > 1 i i v . \ f >o

A VVft b- - ' y< ' > X X ' jSv • ' X y V K-Bs" fcjfS, V/v'Y'V 4 7X -"TWPWVA 1 \X>iTn}vV \ « » \m\'V V x v >04

c^v VI w . f

vVXVWy JOw' vvjffl^ / v \

Figure 14: Experimental shape: (a) 0°/90° fibre orientations and (b) ±45° fibre orientations.

(a) (b)

Figure 15: Predicted shape and shear angles: (a) 0°/90° and (b) ±45° fibre orientations.

*-250-„

Experiment

Figure 16: Geometry of deep drawing tools.

(a) (0790°) (b) (-457+45°)

Figure 17: Experimental results: (a) 0°/90° fibre orientation and (b) -45745° fibre orientation.

(a) FEM (0°/90°) (b) FEM (-45°/45°)

(c) GEOMETRICAL (0°/90°) (d) GEOMETRICAL (-45745°)

Figure 18: Geometrical and FE predicted results 0°/90° and -45°/45° fibre orientations.

Distance to the center (mm) Distance to the center (mm)

^ Experimental

---Geometrical

♦ Experimental

---Geometrical

using rigid tools. As mentioned before, the shear resistance model is assigned to reduced membrane elements, and truss elements parallel to the membrane edges represent the high tensile stiffness in the yarn direction. Figures 14 and 15 report, respectively, the experimental obtained shapes with respect to 0°/90° (Figure 14(a)) and ±45° fibre orientations (Figure 14(a)) and the corresponding final shape using the numerical approach (Figure 15(a) for 0°/90° and Figure 15(b) for ±45° fibre orientations). Likewise, the numerical simulation agrees with the experimental results. ^e angular fibre distortion exceeds 38° along the diagonal axis for 0°/90° fibre orientations and 50° along the diagonal axis for ±45° fibre orientations. Furthermore, we can notice that the final shape obtained with 0°/90° fibre direction is very different from the -45°/+45°. Another interesting result of the numerical calculation is the angular distortion variations between warp and weft fibres.

3.3.3. Deep Drawing ofPrepreg Woven Fabric. ^e numerical analysis of prepreg composite fabric deformation by the deep drawing process is performed by utilizing the commercial FEM package Abaqus/Explicit. ^e resin is modeled by using 1600 membrane finite elements (linear triangular element

M3D3), and warp and weft fibres are modeled by 3200 truss finite elements (linear element T3D2). ^e rigid surface is modeled by 1600 Bezier patches (three nodes R3D3 and four nodes R3D4). ^e behavior of the resin is assumed to be isotropic viscoelastic, and the behaviour of the fibre is supposed to be elastic.

^e first forming example concerns the 3D deep drawing of aramid preimpregnated fabric with conical tools (see Figure 16). ^e mechanical properties of the used material are given in Table 1. Figure 17 reports the experimentally obtained shapes with respect to 0°/90° (Figure 17(a)) and ±45° fibre orientation (Figure 17(b)) for different punch displacements (initial, and final for 100 mm of punch displacement). ^e corresponding final computed shapes are shown in Figure 18(a) for 0°/90° and in Figure 18(b) for ±45° fibre orientations. ^e corresponding final geometrical shapes are shown in Figure 18(c) for 0°/90° and in Figure 18(d) for ±45° fibre orientations. ^ere is good agreement between the experimental and the predicted (mechanical and geometrical) results.

We notice that these shear angle values are very large >38° for mechanical approach along the median line for ±45° fabric (Figure 19(b)) and along the diagonal line for 0°/90° fabric

OE+OO l.i

Figure 22: ] 40 -| ?^0786E+00

0£+00 1.8157.

35 - m _

_2./25bh+[)L

57E+01 ^ ^lAF+m ^-35°

sovalues of shear angles.

2.72365+01

[2] M. Hou, L. Ye, and Y. W. Mai, "Manufacturing process and mechanical properties of thermoplastic composite components," Journal of Materials Processing Technology, vol. 63, pp. 334-338, 1997.

[3] F. Trochu, E. Ruiz, V Achim, and S. Soukane, "Advanced numerical simulation of liquid composite molding for process analysis and optimization," Composites A, vol. 37, no. 6, pp. 890-902, 2006.

[4] G. Marsh, "Composites get in deep with new-generation engine," Reinforced Plastics, vol. 50, no. 11, pp. 26-29, 2006.

[5] P. Boisse, A. Cherouat, J. C. Gelin, and H. Sabhi, "Experimental study and finite element simulation of a glass fiber fabric shaping process," Polymer Composites, vol. 16, no. 1, pp. 83-95, 1995.

[6] A. G. Prodromou, S. V. Lomov, and I. Verpoest, "He method of cells and the mechanical properties of textile composites," Composite Structures, vol. 93, no. 4, pp. 1290-1299, 2011.

[7] S. Kawabata, M. Niwa, and H. Kawai, "He finite-deformation theory of plain-weave fabrics. III. He shear-deformation theory," Journal of the Textile Institute, vol. 64, pp. 21-83, 1973.

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

[9] L. Liu, J. Chen, X. Li, and J. Sherwood, "Two-dimensional macro-mechanics shear models of woven fabrics," Composites A, vol. 36, no. 1, pp. 105-114, 2005.

[10] O. Rozant, P. E. Bourban, and J. A. E. Manson, "Drapability of dry textile fabrics for stampable thermoplastic preforms," Composites A, vol. 31, no. 11, pp. 1167-1177, 2000.

[11] P. Potluri, I. Parlak, R. Ramgulam, and T. V. Sagar, "Analysis of tow deformations in textile preforms subjected to forming forces," Composites Science and Technology, vol. 66, no. 2, pp. 297-305, 2006.

[12] A. Cherouat, J. C. Gelin, P. Boisse, and H. Sabhi, "Numerical modeling of glass woven fabric deep-drawing using finite element method," European Journal of Computational Mechanics, vol. 4, pp. 159-182, 1995.

[13] A. Cherouat and J. L. Billoet, "Mechanical and numerical modelling of composite manufacturing processes deep-drawing and laying-up of thin pre-impregnated woven fabrics," Journal of Materials Processing Technology, vol. 118, no. 1-3, pp. 460-471, 2001.

[14] P. Boisse, K. Buet, A. Gasser, and J. Launay, "Meso/macro-mechanical behaviour of textile reinforcements for thin composites," Composites Science and Technology, vol. 61, no. 3, pp. 395-401,2001.

[15] A. El Hami, B. Radi, and A. Cherouat, "Treatment of the composite fabric's shaping using a Lagrangian formulation," Mathematical and Computer Modelling, vol. 49, no. 7-8, pp. 1337-1349,2009.

[16] M. K. Warby, J. R. Whiteman, W. G. Jiang, P. Warwick, and T. Wright, "Finite element simulation of thermoforming processes for polymer sheets," Mathematics and Computers in Simulation, vol. 61, no. 3-6, pp. 209-218, 2003.

[17] T. C. Lim and S. Ramakrishna, "Modelling of composite sheet forming: a review," Composites A, vol. 33, no. 4, pp. 515-537, 2002.

[18] A. Cherouat, H. Borouchaki, and J. L. Billoet, "Geometrical and mechanical draping of composite fabric," European Journal of Computational Mechanics, vol. 14, pp. 693-707, 2005.

[19] C. Mark and H. M. Taylor, "He fitting of woven cloth to surfaces," Journal of the Textile Institute, vol. 47, pp. 477-488, 1956.

[20] F. Van Der Ween, "Algorithms for draping fabrics on doubly-curved surfaces," International Journal for Numerical Methods in Engineering, vol. 31, no. 7, pp. 1415-1426, 1991.

[21] A. C. Long and C. D. Rudd, "Simulation of reinforcement deformation during the production of preforms for liquid moulding processes," Proceedings of the Institution of Mechanical Engineers B, vol. 208, no. 4, pp. 269-278, 1994.

[22] S. G. Hancock and K. D. Potter, "Inverse drape modelling—an investigation of the set of shapes that can be formed from continuous aligned woven fibre reinforcements," Composites A, vol. 36, no. 7, pp. 947-953, 2005.

[23] K. Vanclooster, S. V. Lomov, and I. Verpoest, "Experimental validation of forming simulations of fabric reinforced polymers using an unsymmetrical mould configuration," Composites A, vol. 40, no. 4, pp. 530-539, 2009.

[24] P. Potluri, S. Sharma, and R. Ramgulam, "Comprehensive drape modelling for moulding 3D textile preforms," Composites A, vol. 32, no. 10, pp. 1415-1424, 2001.

[25] B. Zhu, T. X. Yu, H. Zhang, and X. M. Tao, "Experimental investigation of formability of commingled woven composite preform in stamping operation," Composites B, vol. 42, no. 2, pp. 289-295,2011.

[26] P. Xue, X. Peng, and J. Cao, "A non-orthogonal constitutive model for characterizing woven composites," Composites A, vol. 34, no. 2, pp. 183-193,2003.

[27] R. H. W. Ten Hije, R. Akkerman, and J. Huetink, "Large deformation simulation of anisotropic material using an updated Langrangian finite element method," Computer Methods in Applied Mechanics and Engineering, vol. 196, no. 33-34, pp. 3141-3150,2007.

[28] S. B. Sharma and M. P. F. Sutcliffe, "A simplified finite element model for draping of woven material," Composites A, vol. 35, no. 6, pp. 637-643, 2004.

[29] V. Frishfelds, T. S. Lundstrom, and A. Jakovics, "Bubble motion through non-crimp fabrics during composites manufacturing," Composites A, vol. 39, no. 2, pp. 243-251, 2008.

[30] J. P. Fan, C. Y. Tang, C. P. Tsui, L. C. Chan, and T. C. Lee, "3D finite element simulation of deep drawing with damage development," International Journal of Machine Tools and Manufacture, vol. 46, no. 9, pp. 1035-1044, 2006.

[31] D. Jauffres, C. D. Morris, J. A. Sherwood, and J. Chen, "Simulation of the thermostamping of woven composites: mesoscopic modelling using explicit FEA codes," International Journal of Material Forming, vol. 2, no. 1, pp. 173-176, 2009.

[32] N. Kenkare and T. May-Plumlee, "Evaluation of drape characteristics in fabrics," International Journal of Clothing Science and Technology, vol. 17, no. 2, pp. 109-123, 2005.

[33] N. Kenkare and T. May-Plumlee, "Fabric drape measurement: a modified method using digital image processing," Journal of Textile and Apparel, Technology and Management, vol. 4, no. 3, 2005.

[34] L. Vangheluwe and P. Kiekens, "Time dependence of the drape coefficient of fabrics," International Journal of Clothing Science and Technology, vol. 5, no. 5, pp. 5-8, 1993.

[35] W. R. Yu, M. Zampaloni, F. Pourboghrat, K. Chung, and T. J. Kang, "Analysis of flexible bending behavior of woven preform using non-orthogonal constitutive equation," Composites A, vol. 36, no. 6, pp. 839-850, 2005.

[36] A. Cherouat, H. Borouchaki, and L. Giraud-Moreau, "Mechanical and geometrical approaches applied to composite fabric forming," International Journal of Material Forming, vol. 3, no. 2, pp. 1189-1204, 2010.

[37] O. K. Bergsma and J. Huisman, "Deep Drawing of fabric reinforced thermoplastic," in Proceedings of the 2nd International Conference on Computer Aided Design in Composite Material Technology, pp. 323-333, 1988.

[38] B. Gommers, I. Verpoest, and P. Van Houtte, "Modelling the elastic properties of knitted-fabric-reinforced composites," Composites Science and Technology, vol. 56, no. 6, pp. 685-694, 1996.

[39] Y. Luo and I. Verpoest, "Biaxial tension and ultimate deformation of knitted fabric reinforcements," Composites A, vol. 33, no. 2, pp. 197-203, 2002.

[40] K. A. Padmanabhan, "Metal forming at very low strain rates," in Encyclopedia of Materials: Science and Technology, pp. 5384-5389, 2008.

[41] A. Cherouat and H. Borouchaki, "Present state of the art of composite fabric forming," Geometrical and Mechanical Approaches, Materials, vol. 2, no. 4, pp. 1835-1857, 2009.

[42] R. H. W. ten ^hije, R. Akkerman, and J. Huetink, "Large deformation simulation of anisotropic material using an updated Lagrangian finite element method," Computer Methods in Applied Mechanics and Engineering, vol. 196, no. 33-34, pp. 3141-3150,2007.

[43] Y. F. Dafalias, "Corotational rates for kinematic hardening at large plastic deformations," Journal of Applied Mechanics, vol. 50, pp. 561-565, 1983.

[44] J. K. Dienes, "On the analysis of rotation and stress rate in deforming bodies," Acta Mechanica, vol. 32, no. 4, pp. 217-232, 1979.

Copyright of Advances in Materials Science & Engineering is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.