Pam-Composites Iacmi November2015

ESI’s Composites Simulation Solution
Integrated solution to simulate the manufacturing of

structural composites components
Dr. Xiaoshi Jin

November 2015
www.esi-group.com 1
Copyright © ESI
Copyright
Group, ©
2015.
ESI Group,
All rights
2015.
reserved.
All rights reserved.
What is it?
FEA Solution to define and Optimize process parameters

• Dedicated to composites structural components made of continuous fibers
Dedicated to automated processes

• For mass production
• For full control on costs, delays and reproducibility
Address a wide range of materials
• Carbon, Glass and natural fibers
• Thermoplastic and thermoset resins
• Including core, inserts
• Woven fabrics, NCF, UD…
• Prepregs or dry textiles
www.esi-group.com 2
Copyright © ESI Group, 2015. All rights reserved.
Thermoforming of organosheets
Post-cure distortion of a fuselage panel PAM-COMPOSITES
Performance
Analysis
Infusion of a wind blade

CPD-to-ESI in CATIA V5 XML CATIA-RSO export
Fibersim XML
CAD / DESIGN
www.esi-group.com 3
Compensated mold (CATIA RSO)
Temperature history
Degree of cure history Fiber orientations
www.esi-group.com 4
PAM-FORM Simulation
What for?
To Simulate
• Preforming
• Draping of thermoset prepregs
• Thermoforming of Organosheet
To Determine and Optimize process parameters such as:

• Kinematic of the tools
• Temperature cycle
• Pressure cycle
• Clamping conditions
• Clamping forces
• Initial flat pattern
Through the prediction of

• Wrinkles
• Thicknesses
• Bridging
• Strains (shearing and in fibers)
• Stresses (Shearing and in fibers)
• Fibers orientation
www.esi-group.com 5
PAM-FORM Simulation
Mechanical deformations
modeling
Non-linear elasto-plastic material with damage
• Possible Temperature, Strain rate and Curvature dependency
• Behavior defined through following input:
• Tension and Compression deformation in fiber directions
• Tension is elastic with damage or plastic (transverse
direction of UD)
• Compression is elastic
• In-plane shear deformation
• Elastic with damage or plastic
• Bending deformation in fiber directions
• Thickness deformation through normal pressure
www.esi-group.com 6
PAM-FORM Simulation
modeling
Material compressibility
• Whether the material is defined as incompressible
• The volume of each element remains constant during the simulation:
epsilon1+epsilon2+epsilon3=0
• Whether the material is defined as compressible
• Thickness is then dependant of the shear deformation and/or pressure
• Possibility to visualize compaction ratio (=initial volume of
element/final volume of element) as an output
www.esi-group.com 7
PAM-FORM Simulation
modeling
Elastic damage behavior
• Damage behavior is automatically activated after each necking point
detected in the stress-strain curves
• Computed damage amount express the broken fiber ratio
• It is approximated as the ratio between the damaged stress amount
and the assumed ideal tensile stress when no damage
www.esi-group.com 8
PAM-FORM Simulation
Thermal modeling
Thermal phenomena coupled to mechanical analysis

• Interesting when thermal dependency in material properties
(prepreg materials)
• The following phenomena can be taken into account:
• Conductivity in the composites material
• Conductivity in the tools
• Heat transfer between composites layers
• Heat transfer between composites and tools
• Convection exchanges
• Radiation exchanges
www.esi-group.com 9
Thermoplastic forming Application
Flap rib thermoforming
Rubber pad forming simulation
Bridging effect Fiber Shearing

www.esi-group.com 10
Wrinkling prediction
Wrinkle
Smooth
surface
Wrinkling prediction
20 plies carbon UD / APC2-AS4
Flat pattern optimization
Upper tool
4 plies
Lower tool
Initial flat Optimized flat
pattern pattern Simulation setup
Poor part
quality Improved part
quality
Flat pattern optimisation

4 plies / thermoplastic matrix
Wing box thermoforming
Laminate thickness
Ultrasonic measurement versus simulation
Top section
Thickness per ply
Bottom section
X=0 480 X=480
Process animation
1,2 web thicknes (top section)
1
0,8
ick nes
(m m)
0,6
0,4 thicknes (trial)
0,2 web thickness (bottom section)
1.2
0 1
th (mm)
0.8
s
thickness
0.6
0.4 thickness (trial)
0.2 thickness (PAM-FORM)
0
0 100 200 300 400
chordwise location
Preforming Application
Fiber shearing
PAM-RTM Simulation
What for?
To Simulate
• RTM
• VARI
• Light RTM
• C-RTM
• Curing
To Determine and Optimize process parameters such as:

• Location of injection gates
• Location of vents
• Position and size of flow media
• Heating of the mold
• Cure cycle

• Air traps
• Micro porosity
• Injection time
• Curing time
• Temperature evolution
• Degree of cure evolution
• Pressure in the mold
PAM-RTM Simulation
Flow in porous media
Darcy’s law
• Darcy’s law: • Resin mass conservation:
Boundary conditions required

to solve the equation:
- imposed pressure
- or imposed flow rate at the inlet
PAM-RTM Simulation
Thermal phenomena
Temperature of the resin

• Governs the reactivity of the polymerization reaction
• Influences the filling since the viscosity of the resin is temperature
dependent
• Temperature field is governed by following equation:
PAM-RTM Simulation
Material data
Permeability of the reinforcements

• Permeability is fiber volume content dependent (degree of compaction)
• It also depends on the draping of the reinforcements (fiber “shearing”)
• Fiber volume content and draping results can be computed with PAM-
FORM and imported in PAM-RTM
Viscosity of the resin

• It is temperature and degree of conversion dependent
Kinetics of the resin

• Defined through pre-defined Kamal-Sourour models or User defined
models
PAM-RTM Simulation
Summary
Coupling of physical phenomena
PAM-RTM Simulation
To be released in 2016
Fluid-Mechanics coupling solver

• Includes Navier-Stokes law
• Introduced to treat in a more accurate way processes as
Compression-RTM and Infusion RESIN
TOOL VELOCITY
INJECTION
TOOL VELOCITY
Moving Tool Resin gate TOOL
Moving Tool
Moving Tool PRESSURE
Void
Void
Gap
GAP
Preform
Resin + Preform
Resin + Preform
Resin + Preform
Vent Dry Preform
Dry Preform
Fixed tool
Mold
Infusion Application
Large dimensions parts
Filling time Curing time
30 meters wind spar cap infusion Fan blade of an aircraft engine RTM
Filling degree
Fuselage panel infusion

Light RTM Application
Bus body side
• 13*6.5 feet
• fiberglass Chopped Strand Mat (CSM)
• integrated flow media
• Polyester resin
"In this project, the RTM simulation helped us to

secure and to optimize the process. Today, we are
using ESI’s PAM- RTM not only to assess process
parameters, including injection time and pressure
in mold, but also to fine-tune mold design.”
Optimum resin flow pattern computed
Jérôme RAYNAL with PAM-RTM
Sales and Export Director
Pôle de Plasturgie de l’Est
Inner liner for hull reinforcement
Very complex part including high shapes (1.2m)

Injection analysis allows determination of injection strategy
Initial injection points
(injection points/channels and vents location as well as
open/closing sequence)
to minimize:
Dry spots
Filling and curing times
Fiber washing
Pressure in the mold
Secondary injection points/channels
Flow front during injection

Vents location
TECABS project
Renault – Mines de Douai
RTM Application
Automotive floor panel
Floor pan
injection
Wind blade
Flow media
Zoom-in on
flow media
influence on
resin flow front
Aeronautic part
Simulation helped to divide

infusion time by almost 5
while producing better quality
parts with less scraps
(better reproducibility)
RTM Application
Automatic reduction of porosities
Porosity prediction and reduction:

Principle: Critical impregnation velocity:  
max v front  vcrit
RTM Application

PAM-RTM input data:
14%
Ca= viscosity*velocity/surface tension/contact angle
Chomarat - Roviply
12%
12.824-1573.7*(V) 20 psi
10% 30 psi
45 psi
Optimum velocity
Voids content
2,5 ml/s
8%
Based on experiment
6%
4% Macro-void 1.26+100.55*(V)
2%
Micro-void
0%
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Flow velocity (m/sec)
Experimental data PAM-RTM input curve
RTM Application

PAM-RTM output: injection flow rate curve
(a) Constant injection pressure

7 1E-02
Capillary number at flow front

injection flow rate 1e-7 [m3/s]
6
4
1E-03
3
(b) Constant injection flow rate

2
0 1E-04
0 50 100 150
injection time [sec]
(c) Optimized injection flow rate
PAM-DISTORTION Simulation
What for?
To simulate
• Process induced distortions such as spring-in and warping
To Define and Optimize parameters such as:

• Stacking definition
• Curing cycle
• Mold material and design
• Mold geometrical compensation

• Internal stresses during curing
• Residual stresses after de-molding
• Deformation during curing
• Deformation after de-molding
Takes into account

• Material history during curing process (temperature and degree of cure)
• Thermal and Mechanical interaction with the mold
Physical phenomena
PAM-DISTORTION computes
• Thermal strains
• Chemical strains
It takes into account resin phase transformation during curing

• Initially the resin is liquid: no stress, no strain
• When the resin reaches gelation (alpha=alpha gel) the resin becomes rubbery.
From this point the resin can sustain stresses.
• When the resin reaches the glass transition temperature (Temperature of the
resin = glass transition temperature) the resin becomes glassy.
• Glass transition temperature evolution is defined using Di Benedetto function.
Physical phenomena
Phases of the resin for 1 element during curing
PAM-DISTORTION simulation
Fuselage Panel
Deformed shape / springback

Temperature evolution during curing
Initial and deformed geometries
Degree of cure evolution during curing

33
KTM Race Car Application
Front Splitter
Front Splitter
Laminate definition
Zoom-in on structure
Fiber: Tenax–J STS 40 F13 24 K 1800 tex; Binding type: plain
Resin+ Hardener: EpikoteResin04695/1 and EpikureCuringAgent 05357
Core: Büfadur67 –15, PUR Foam
Preforming of upper plies
With PAM-FORM
Upper tool’s
displacement
Upper
tool
Lower tool
(lock)
Preforming Simulation
To predict and optimize
Thickness distribution
Thinning contour
wrinkles
Shearing
Shear angle
Fiber orientations
Fiber orientations
Export your fiber

orientation for injection
simulation or structural
analysis
Injection Simulation
Dry spots & Porosities
3D Injection with PAM-RTM

• Takes into account fiber orientations (impact on permeability) from preforming simulation
Bottom view Top view
Injection Simulation
Filling time
Filling time
Bottom view Top view
Selected Industrial References
T H A N K YO U
[email protected]
Copyright © ESI Group, 2014.
2015. All rights reserved.

Pam-Composites Iacmi November2015

Uploaded by

Copyright:

Available Formats

Pam-Composites Iacmi November2015

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Pam-Composites Iacmi November2015

Uploaded by

Copyright:

Available Formats

What is ESI's Composites Simulation Solution and what does it do?

What is ESI's Composites Simulation Solution and what does it do?

What types of materials and processes does the solution address?

What types of materials and processes does the solution address?

ESI’s Composites Simulation Solution

Integrated solution to simulate the manufacturing of

Dr. Xiaoshi Jin

FEA Solution to define and Optimize process parameters

Dedicated to automated processes

Post-cure distortion of a fuselage panel PAM-COMPOSITES

Infusion of a wind blade

Compensated mold (CATIA RSO)

To Determine and Optimize process parameters such as:

Through the prediction of

Thermal phenomena coupled to mechanical analysis

Rubber pad forming simulation

Bridging effect Fiber Shearing

Flat pattern optimisation

To Determine and Optimize process parameters such as:

Through the prediction of

• Darcy’s law: • Resin mass conservation:

Boundary conditions required

Temperature of the resin

• Temperature field is governed by following equation:

Permeability of the reinforcements

Viscosity of the resin

Kinetics of the resin

Coupling of physical phenomena

Fluid-Mechanics coupling solver

Fuselage panel infusion

"In this project, the RTM simulation helped us to

Very complex part including high shapes (1.2m)

Secondary injection points/channels

Flow front during injection

Simulation helped to divide

Porosity prediction and reduction:

Porosity prediction and reduction:

Flow velocity (m/sec)

Experimental data PAM-RTM input curve

Porosity prediction and reduction:

(a) Constant injection pressure

Capillary number at flow front

(b) Constant injection flow rate

(c) Optimized injection flow rate

To Define and Optimize parameters such as:

Through the prediction of

Takes into account

It takes into account resin phase transformation during curing

Phases of the resin for 1 element during curing

Deformed shape / springback

Initial and deformed geometries

Degree of cure evolution during curing

Export your fiber

3D Injection with PAM-RTM

Bottom view Top view

Bottom view Top view

You might also like