230 likes | 262 Views
Explore simulations of fragmenting cylinders using the Material Point Method (MPM) for fluid-structure interaction without mesh entanglement. Learn about stress update, plasticity modeling, damage/failure modeling, fracture simulation, and validation processes with different materials and cell counts.
E N D
Material Point Method Simulations of Fragmenting Cylinders Biswajit Banerjee Department of Mechanical Engineering University of Utah 17th ASCE Engineering Mechanics Conference, 2004
Outline • Scenario • Material Point Method (MPM) • Approach • Validation • Simulations of fragmentation
Simulation Requirements • Fire-container interaction • Large deformations • Strain-rate/temperature dependence • Failure due to void growth/shear bands
Tightly-coupled fluid-structure interaction. No mesh entanglement. Convenient contact framework. Mesh generation trivial. Easily parallelized. No tensile instabilities. First-order accuracy. High particle density for tension dominated problems. Computationally more expensive than FEM. Why MPM ? Advantages Disadvantages
Stress update • Hypoelastic-plastic material • Corotational formulation (Maudlin & Schiferl,1996) • Semi-implicit (Nemat-Nasser & Chung, 1992) • Stress tensor split into isotropic/deviatoric • Radial return plasticity • State dependent elastic moduli, melting temperature
Plasticity modeling • Isotropic stress using Mie-Gruneisen Equation of State. • Deviatoric stress : • Flow stress : Johnson-Cook, Mechanical Threshold Stress, Steinberg-Cochran-Guinan • Yield function : von Mises, Gurson-Tvergaard-Needleman, Rousselier • Temperature rise due to plastic dissipation • Associated flow rule
Damage/Failure modeling • Damage models: • Void nucleation/growth (strain-based) • Porosity evolution (strain-based) • Scalar damage evolution: Johnson-Cook/Hancock-MacKenzie • Failure • Melt temperature exceeded • Modified TEPLA model (Addessio and Johnson, 1988) • Drucker stability postulate • Loss of hyperbolicity (Acoustic tensor)
Fracture Simulation • Particle mass is removed. • Particle stress is set to zero. • Particle converted into a new material that interacts with the rest of the body via contact.
Validation: Plasticity Models 635 K 194 m/s 718 K 188 m/s JC MTS SCG JC MTS SCG 655 K 354 m/s 727 K 211 m/s 6061-T6 Aluminum EFC Copper
Validation: Mesh dependence 18,900 cells 151,000 cells 1,200,000 cells OFHC Copper 298 K 177 m/s MTS 11,500 cells 91,800 cells 735,000 cells 6061-T6 Al 655 K 354 m/s JC
Validation: Penetration/Failure 160,000 cells 1,280,000 cells
Validation: 2D Fragmentation JC (steel), ViscoScram (PBX 9501) MTS (steel), ViscoScram (PBX 9501) Gurson-Tvergaard-Needleman yield, Drucker stability, Acoustic tensor, Gaussian porosity, fragments match Grady equation, gases with ICE-CFD code.