13e rencontre smai math industrie
play

13e rencontre SMAI math industrie Physically correct and timely - PowerPoint PPT Presentation

13e rencontre SMAI math industrie Physically correct and timely simulations: the challenge towards Virtual Prototyping Argiris Kamoulakos, ESI Group 22/03/2011 SMAI HPC 1 Introduction Today's industry trend is towards reduction of


  1. 13e rencontre SMAI math ‐ industrie Physically correct and timely simulations: the challenge towards Virtual Prototyping Argiris Kamoulakos, ESI Group 22/03/2011 SMAI ‐ HPC 1

  2. Introduction • Today's industry trend is towards reduction of physical prototyping through simulation driven design in order to reduce drastically the time ‐ to ‐ market of new products. This is called Virtual Prototyping. However, it relies heavily on the capacity to reproduce through simulation the correct physics and in a timely way to respect the industrial development process deadlines. 22/03/2011 SMAI ‐ HPC 2

  3. ESI Group Offering Software & Services for End ‐ to ‐ End Virtual Prototyping 22/03/2011 SMAI ‐ HPC 3

  4. Enjeux / problématique industrielle • Crash Failure Analysis as an industrially affordable solution • Must be performed “over-night” 22/03/2011 SMAI ‐ HPC 4

  5. A typical problem: automotive side ‐ impact B ‐ Pillar model 22/03/2011 SMAI ‐ HPC

  6. The Wilkins (EWK) metal rupture model - Mathematical details - A cumulative strain damage model Non proportional loading Plasticity based Triaxiality dependence Lode angle dependence 22/03/2011 SMAI ‐ HPC

  7. Importance of "Size effect" A severe condition imposed on model element size in order to be predictive Size effect: • critical volume to be saturated with damage in order to start a “crack” • Typical values • Aluminum: (0.15mm)**3 • ~6 solid elements per mm • High Strength Steel: (0.05mm)**3 • 15 ‐ 20 solid elements per mm 22/03/2011 SMAI ‐ HPC 7

  8. High fidelity modelling of the B ‐ Pillar ~1,000,000 solids ~4,000 shells 22/03/2011 SMAI ‐ HPC

  9. The industrial challenge The industrial challenge  Crash Failure Analysis as an industrially affordable solution  Example  Must be performed “over-night” 22/03/2011 SMAI ‐ HPC

  10. Towards extreme local refinement Towards extreme local refinement Car Model: Subframe (520k-Model): 520,000 TET10-Elements 879,000 Elements Time Step: 1 μ s Time Step: 0.05 μ s + Time step ratio : R  t = 20 Local model element fraction : f L = 59 % Problem: It is not just a matter of “size”, but algorithmic, too. (Disproportionate distribution of solution effort)

  11. Casted parts pose a special issue Gas Porosity, Shrink Porosity, Micro Porosity no porosity with porosity SMAI ‐ HPC 22/03/2011

  12. Very detailed models are needed for the prediction of the grain structure Uniform Solidification in Ni Turbocharger CAFE (Cellular Automata FE) was developed in collaboration with: Howmet, Snecma (France), ABB, EPFL (Switzerland), AETC, Rolls Royce (UK), PCC (USA) SMAI ‐ HPC 22/03/2011 12

  13. In Aeronautics the problem of Birdstrike is particularly accute SPH Courtesy of EC project CRAHVI 22/03/2011 SMAI ‐ HPC

  14. Incorporation of details is critical for performance prediction You cannot predict what you do not model for ! Courtesy of EC project CRAHVI The rivets and their effect in rupture have to be sufficiently represented by the model. SMAI ‐ HPC 14 22/03/2011

  15. Les points de blocage • Mesh size – huge meshes needed • Computational algorithm – Dissimilar meshes in size and type (FE / meshless) • Initial state of structure – Manufacture (filamentary composites) – Connectors (spotwelds, rivets etc.) 22/03/2011 SMAI ‐ HPC 15

  16. • Manufacture of Fillamentary Composites 22/03/2011 SMAI ‐ HPC 16

  17. New Manufacturing methods for advanced composites 17

  18. New Manufacturing methods for advanced composites 18

  19. Virtual Properties prediction (through virtual coupon testing) Particle methods for the resin SMAI ‐ HPC 19 22/03/2011

  20. Virtual Properties prediction (through virtual coupon testing) Particle methods for the resin SMAI ‐ HPC 20 22/03/2011

  21. • Modelling connectors (Spotwelds, Rivets etc.) 22/03/2011 SMAI ‐ HPC 21

  22. Effect of connectors in automotive crash Spotewelded beam modelling Deformed Deformed spotwelds beam 22 22/03/2011 SMAI ‐ HPC

  23. Effect of connectors in automotive crash Spotewelded beam modelling 23 22/03/2011 SMAI ‐ HPC

  24. The challenge of Self Piercing Rivets: knowing their shape and state 22/03/2011 SMAI ‐ HPC

  25. The challenge of Self Piercing Rivets: knowing their shape and state 22/03/2011 SMAI ‐ HPC

  26. Les points de blocage • Mesh size – huge meshes needed • Computational algorithm – Dissimilar meshes in size and type (FE / meshless) 22/03/2011 SMAI ‐ HPC 26

  27. • Modelling substructures and systems 22/03/2011 SMAI ‐ HPC 27

  28. Application to Battery Modelling Pr Prob oblema lematics: ma tics: material issues terial issues Mic Microscale oscale structur cture ( e (20 0 μ m) m) Discontin Discontinuous uous Nano-par Nano-particles les, nano-f nano-fibers Modeling Modeling Continuous mec Contin uous mechanics anics Issues Iss Contin Continuous mec uous mechanics anics pr proper operti ties es Interpreta Inter tation tion Calibr Calibration tion Ca Cath thodes Separat rators rs Anod Anodes SMAI ‐ HPC 22/03/2011 28

  29. Application to Battery Modelling Structur ctural la al layout: ut: shaping: bending, stretching (?), drawing (?) 200 layers 4mm 20 µm plastic covers clamping: pressure Requirements: Re Str Stresses sses and str and strains in eac ins in each h la layer Result : 200 sult : 200 elements/mm elements/mm A 20cm x 20cm ensemb A 20cm x 20cm ensemble le giv gives s 32 Million 32 Million elements elements SMAI ‐ HPC 22/03/2011 29

  30. Ultimate goal: be predictive at the system level EC project HUMOS demonstrator of real car crashes HUMOS Bullet car speed = 48 km/h Target car speed = 24 km/h HUMOS 50% male model: Driver of target car Hands attached on steering wheel Right foot on breaking pedal Left foot on the floor Restraint devices: belts Injury assessment SMAI ‐ HPC 22/03/2011 30

  31. Marine applications are the next frontier Motion in Waves: Regular head waves of 8 meter height and a wave length of 294 meter 22/03/2011 SMAI ‐ HPC

  32. • What are the answers ? • What needs to be assessed / improved ? 22/03/2011 SMAI ‐ HPC 32

  33. Multiscale ??? 22/03/2011 SMAI ‐ HPC

  34. Parallelisation ??? (DMP etc.) 4 3.5 3 2.5 Time (Hours) 2 Series1 1.5 1 0.5 0 0 20 40 60 80 100 120 140 Number of CPUs 8 processors : a bit 128 processors : less more than 3.5 hours than half hour 22/03/2011 SMAI ‐ HPC

  35. Infrastructure ??? 1 ms 22/03/2011 SMAI ‐ HPC

  36. GPUs ??? For FE and mesheless !!! 22/03/2011 SMAI ‐ HPC 36

  37. Methodology ??? FE versus meshless for fracture or XFEM etc. 22/03/2011 SMAI ‐ HPC 37

  38. Ultimate goal: be predictive at the system level EC project HUMOS demonstrator of real car crashes HUMOS Mid ‐ sagittal section SMAI ‐ HPC 22/03/2011 38

  39. Towards virtual prototyping at large scale … Courtesy of EADS ‐ CASA 39

  40. …and "full system" predictability Model 1 Model 2 32 processes 128 processes 22/03/2011 SMAI ‐ HPC 40

  41. 22/03/2011 SMAI ‐ HPC 41

Recommend


More recommend