Numerical analysis and optimal design of composite thermoforming process

Abstract In this paper, we present a numerical methodology to simulate and design the composite thermoforming process. First, in the process simulation, the FEM modeling including flow, heat transfer and residual stress is developed to analyze the thermoforming process of composite laminates with woven fabric microstructures. The governing equations and material properties of thermoplastic composites at the forming temperature are obtained by the homogenization method based on the assumption of instantaneously rigid solid fibers suspended in a viscous non-Newtonian polymer melt. The processing theology of the composites is characterized by a power-law constitutive model for this non-isothermal and shear thinning fluid. To simulate the shaping and cooling stages of the entire forming process, the three-dimensional finite element method incorporating a fiber orientation model of woven-fabric microstructures is developed. Based on kinematic considerations, the fiber orientation model is used to keep track of the microstructure change due to large deformation. This global—local numerical methodology is capable of predicting macroscopic and microscopic deformation mechanics during processing. Second, in the process optimization, the design objective is to achieve a uniform final thickness distribution by designing the preheated composite temperature field. The finite difference approximation is utilized to evaluate the design sensitivities. Using the FEM developed in the first part as the analysis tool, the optimization scheme can allow process designers to optimally design this process. A number of numerical examples and experimental results are presented to demonstrate this combined analysis and optimization algorithm.

[1]  A. Ogale,et al.  Influence of microstructure on elastic and viscoelastic properties of polyether ether ketone , 1987 .

[2]  T. Osswald,et al.  Prediction of Shrinkage and Warpage of Fiber Reinforced Thermoset Composite Parts , 1994 .

[3]  G. Jeronimidis,et al.  Residual Stresses in Carbon Fibre-Thermoplastic Matrix Laminates , 1988 .

[4]  First-Order Approximations for the Effective Shearing Viscosities of Continuous-Fiber Suspensions , 1995 .

[5]  K. Potter,et al.  The influence of accurate stretch data for reinforcements on the production of complex structural mouldings , 1979 .

[6]  Garrett B. McGuinness,et al.  Numerical analysis of stresses and deformations in composite materials sheet forming: central indentation of a circular sheet , 1993 .

[7]  N. Pagano,et al.  Curing Stresses in Composite Laminates , 1975 .

[8]  Jan-Anders E. Månson,et al.  Prediction of Process-Induced Residual Stresses in Thermoplastic Composites , 1990 .

[9]  Jean-Loup Chenot,et al.  OPTIMAL DESIGN FOR NON-STEADY-STATE METAL FORMING PROCESSES. II: APPLICATION OF SHAPE OPTIMIZATION IN FORGING , 1996 .

[10]  van der Erik Giessen,et al.  SIMULATION OF MATERIALS PROCESSING: THEORY, METHODS AND APPLICATIONS , 1998 .

[11]  D. Allman A compatible triangular element including vertex rotations for plane elasticity analysis , 1984 .

[12]  R. E. Robertson,et al.  Continuous fiber rearrangements during the molding of fiber composites. II. Flat cloth to a rounded cone , 1984 .

[13]  F. N. Cogswell,et al.  The Processing Science of Thermoplastic Structural Composites , 1987 .

[14]  O. C. Zienkiewicz,et al.  Flow of solids during forming and extrusion: Some aspects of numerical solutions , 1978 .

[15]  Nicholas Zabaras,et al.  A sensitivity analysis for the optimal design of metal-forming processes , 1996 .

[16]  Richard J. Day,et al.  Fibre deformation and residual thermal stresses in carbon fibre reinforced PEEK , 1989 .

[17]  T. Lévy,et al.  Application of the homogenization method to a suspension of fibres , 1994 .

[18]  Scott J. Hollister,et al.  Digital-image-based finite element analysis for bone microstructure using conjugate gradient and Gaussian filter techniques , 1993, Optics & Photonics.

[19]  E. Sanchez-Palencia Non-Homogeneous Media and Vibration Theory , 1980 .

[20]  A. Spencer,et al.  Deformations of fibre-reinforced materials, , 1972 .

[21]  T. G. Rogers Rheological characterization of anisotropic materials , 1989 .

[22]  J. Seferis,et al.  Crystallization kinetics of polyetheretherketone (PEEK) matrices , 1986 .

[23]  J. Kirkhope,et al.  An alternative explicit formulation for the DKT plate-bending element , 1985 .

[24]  A. Dasgupta,et al.  Orthotropic Thermal Conductivity of Plain-Weave Fabric Composites Using a Homogenization Technique , 1992 .

[25]  D. Groves A characterization of shear flow in continuous fibre thermoplastic laminates , 1989 .

[26]  Y. Weitsman,et al.  Residual Thermal Stresses Due to Cool-Down of Epoxy-Resin Composites , 1979 .

[27]  The use of friction in the shaping of a flat sheet into a hemisphere , 1996 .

[28]  J. A. Barnes,et al.  The formation of residual stresses in laminated thermoplastic composites , 1994 .

[29]  E. Sanchez-Palencia,et al.  Suspension of solid particles in a newtonian fluid , 1983 .

[30]  Frederic Neil Cogswell,et al.  Thermoplastic aromatic polymer composites : a study of the structure, processing, and properties of carbon fibre reinforced polyetheretherketone and related materials , 1992 .

[31]  R. Byron Pipes,et al.  Finite element analysis of composite sheet-forming process , 1991 .

[32]  Eugenio Oñate,et al.  A viscous shell formulation for the analysis of thin sheet metal forming , 1983 .