140 likes | 251 Views
Adaptive Multiscale Modeling and Simulation for Munitions Simulations* Progress Report. PIs: Jacob Fish and Mark S. Shephard Post-docs: Gal Davidi, Caglar Oskay Students: Zheng Yuan, Rong Fan. *AFRL support leveraged by support from NSF, ONR and General Motors. Roadmap of Developments.
E N D
Adaptive Multiscale Modeling and Simulation for Munitions Simulations*Progress Report PIs: Jacob Fish and Mark S. Shephard Post-docs: Gal Davidi, Caglar Oskay Students: Zheng Yuan, Rong Fan *AFRL support leveraged by support from NSF, ONR and General Motors
Roadmap of Developments • Assessment of commercial code capabilities • Mesh sensitivity studies (Gal Davidi) • Validation studies (Rong Fan) • Fragmentation capabilities for metals • Homogenization based approach (Gal Davidi) • Integration of homogenization in ABAQUS (Zhen Yuan) • PUM based (Zhen Yuan and Rong Fan)
Roadmap of Developments (cont) • Fragmentation capabilities for composites • Reduced order methodology (Oskay) • Validation studies (Oskay) • Integration in ABAQUS • Multiscale Enrichment based PUM • Applications
Target Assembly Impactor in Sabot DH36 Steel Plate Fragmentation in MetalsExperimental setup The experimental parameters considered: • Steel target plate: DH36 steel; 3/16 inch thick; 6 inch diameter; • Impact velocity: In the range between 920 ft/sec. • Backing material: Polyurea: 0.215 inch • Impactor: non-deformable
Experiment vs ABAQUS simulation (without backing) Mises stress (without backing) Equivalent plastic strain (without backing)
Shell 21 layers 3D model (21 layers) Drawbacks of commercial software • Cost of 3D simulations (4 days for 21 layer-model, r-adaptivity) • Mesh dependency of both 3D and shell models 3.00E+05 3D models (4-8 layers) 2.50E+05 2.00E+05 1.50E+05 Fine (160) 1.00E+05 Coarse (80) Very Coarse (40) 5.00E+04 0.00E+00 0 0.00005 0.0001 0.00015 0.0002 0.00025 0.0003 0.00035 0.0004 0.00045
Remedy: Multiscale Enrichment • Global (structure) Enrichment • Enrich the kinematics of the global mesh with failure characteristic (delamination, shear banding, fragmentation) characteristic computed on the local patch • For computational efficiency • Local (material) Enrichment • Embed discontinuities (strong or weak) into material (micromechanical) model • For regularization of failure models
Global deformation modes Failure deformation mode-shapes Cell problems on delamination fracture Global Enrichment (MEPU) Better (Superposition) (Domain decomposition)
DH36 & ERC (Shell) Global Enrichment (metals)3D simulations 300 DH36 & ERC (3D-21 layers) 250 MEPU 200 Velocity of Impactor (m/s) 150 100 50 0 0 0.0001 0.0002 0.0003 0.0004 Time (s)
4 7 8 3 1 6 5 2 Discontinuity plane RVE Local Enrichment (metals)(in progress) • Calculate discontinuity direction at each Gauss point • Align the RVE local coordinate system with one of the axis normal to the localization plane • Develop a 3-point RVE model as follows: Constrained RGB Gauss point Constrained periodicity Shell master
Phenomenological Advantages Fast Disadvantages Reliability Experiments architecture dependent Eigendeformation-based Reduced Order Homogenization Material Point Component Matrix point (s) • Engineering Accuracy • Fast • Architecture independent Experiments Interface point (s) Fiber point (s) Impact Fragmentation of composites Direct Homogenization • Advantages • Reliability • Architecture independent Exp. • Disadvantages • Computationally formidable
Validation: Tube Crush Experiment • Experiments by Oak Ridge (Starbuck et al.) • Impact Velocity: 4000 mm/sec • Microstructure: Woven composite