550 likes | 648 Views
Composites Forming Analysis. Remko Akkerman www.utwente.nl/ctw/pt r.akkerman@utwente.nl 26 th September 2013. Introduction. Scope Mechanisms Constitutive Models Implementation. Freedom of Design. The sky is the limit? Limits in FORMABILITY Which, why, where & how?.
E N D
Composites Forming Analysis RemkoAkkerman www.utwente.nl/ctw/pt r.akkerman@utwente.nl 26th September 2013
Introduction • Scope • Mechanisms • Constitutive Models • Implementation
Freedom of Design • The sky is the limit? • Limits in FORMABILITY • Which, why, where & how?
Composite Life line • What is a material, what is a structure? • What is a Forming Process? • Micro is close to Meso is close to Macro...
Impregnation& consolidation quality residualstresses productdistortions joining,welding & bonding environmental loading mechanically induced stresses recycling crack initiation & crack growth After life Composite life line
product geometry processsettings fibreorientation product properties composite properties fibre/matrixproperties Interrelations:Processing, Properties & Performance
FormingProcesses • Consolidation • Drape (pre-forming) • Press Forming • Compression Molding • ....
Deformation Limits • “Form ability” • Low resistance to shear & bending • High anisotropy • Negligible fibre extension • Low compressive “strength” (fibre buckling)
Formability Analysis... • From deformation mechanisms • ... to material characterisation • ... to constitutive modelling • ... to process modelling • ... and formability prediction
MaterialCharacterisation • Intra-ply shear (a) Picture frame. (b) Bias extension.
Material Characterisation • Bi-axial response Crimp leads to non-linear behaviour depending on the warp/weft strain ratio
Material Characterisation • Ply/tool and Ply/ply Friction Tool/ply friction (glass/PP) Shear stress vs pressure.
ContinuumMechanics RECAP: ContinuumMechanics = Balance equations + Material ‘Laws’ + Formalism
ContinuumMechanics Balance Equations • Conservation of mass • Conservation of energy • Conservation of momentum Material ‘Laws’ • Constitutiveequations,relatingforces & fluxes Formalism • Scalars, vectors, tensors • Deformationtheories
Balance Equations • Conservation of mass • Conservation of momentum • Conservationof energy (1stLaw)
ConstitutiveEquations • Relations betweenFluxes(transport of anextensivequantity) e.g. • andForces(gradient of an intensive quantity) e.g. • or, indeed, between stressesandstrains / strainrates e.g.
Formalism • Scalars: e.g. • Vectors: e.g. • Tensors: e.g.
Formalism • Single contraction, • Dyadicproduct,
CompositesFormingProcessesbalanceequations • Viscous & elasticforces dominant (low Reynoldsnumber)Neglectinertia: • Neglectalso body forces: Stress equilibrium • Neglectcoolingduringforming(at leastinitially)
CompositesFormingProcessesconstitutiveequations • Matrix response:Viscousvisco-elasticelasticlow modulus, O(1 MPa) • Fibre response: Elastichigh modulus, O(100 GPa) • Prepreg/laminate response:Elastic/high modulus - in fibredirectionVisco-elastic/low modulus - transverse dir.
CompositesFormingProcessesconstitutiveequations Concluding: • Very high anisotropy • Large rotations & deformationspossibleexcept in the fibredirection udply woven fabric
Reinforcement structures… some terminology • Unidirectional • Biaxial (weft & warp) • Triaxial • ….
warp fill 1 2 Textiles: Woven Fabrics plain 3x1 twill 2x2 twill 5H satin
b a Fibre Directions • unit vectors a, b • deformation gradient F • rate of deformation D
b a a' b' Fibre Directions deformation
ConstitutiveEquationsdefinition of strain Strain definition: Frame of reference: Which“l”? Total Lagrange or Updated Lagrange? Differential calculus:
ConstitutiveEquationsdefinition of strain 3D Straindefinition: Goodforlinearelasticity But does itworkforCompositesForming?
ConstitutiveEquationsdefinition of strain Rigidrotation: Often non-zero axialstrain Exceptfor the “averageconfiguration”
ConstitutiveEquationsdefinition of strain Averageconfiguration: But in whichdirection does the stress act? Shouldbe in the FinalConfiguration! (considering the high anisotropy) INCONSISTENCY
ConstitutiveEquationsdefinition of strain Result(tensile test simulation, E1/E2=105): Exact straindefinitionrequired
ConstitutiveEquationsdefinition of strain Large deformationtheory Deformationgradient: andalso:
ConstitutiveEquationsdefinition of strain The usualpolardecomposition: (Rorthogonal, V&Usymmetric) maintainsanorthogonal basis which is usually wrong!
ConstitutiveEquationsdefinition of strain Solution: multiplicative split (Rorthogonal, G non-symmetric), knowing suchthat andhence leadingto as the scalarfibrestrainϵ in directiona with
Continuum model • Recall incompressible isotropic viscous fluids: • Now directional properties f (a,b)
Continuum model • Inextensibility: • or introduce • leads to
Continuum model • Incompressibility: • Combine with • leads to
Continuum model • extra stress t • Form-invariance under rigid rotations:isotropic function of its arguments • Assume linearity, leads to: with
Continuum model • Fabric Reinforced Fluid (FRF) model • Can be simplified by symmetry considerations (sense of a, b, fabric symmetry)
ConstitutiveModelling • Continuummechanics • Alternative: Discrete approach (resin + fibre + structure) forinstanceusingmesoscopicmodelling
Mesoscopic modelling • Composite property prediction from mesostructure Shear response from FE model
Mesoscopic modelling • Composite property prediction from mesostructure 3D Biaxial 2D Triaxial 2D Knit Multiaxial 2D (NCF) TexGen, WiseTex, etc
Implementation issues • Accuracy especially concerning fibre directions • Consistent tangent (as above) • Shear locking (due to large stiffness differences)
Shear Locking Linear triangle (N1, N2, N3) Linear strains & rotations
Shear Locking Fibres in x and y direction (inextensibility) Eliminate rigid body displacements
N3 y N1 x N2 Shear Locking N1 in the origin (0,0) Remaining d.o.f.s
Shear Locking Suppress a single node Ni(i=2,3) Shear locking ! Unless: xi=0or yi=0(i=2,3) Edge coincides with fibre direction!
Shear Locking • Result of locking: • Far too high stiffness • Spuriouswrinkles • Incorrect deformations • Example: bias extension
Shear Locking • Alignedvsunalignedmesh (quads)