LS-DYNA General Update - Oasys Software

Status and Development Plan Oasys LS-­‐DYNA 6th Annual Update Mee

LS-­‐TaSC v2Was LS-­‐OPT/Topology for V1; renamed as LS-­‐TaSC, Topology and Shape Computa

LS-­‐DYNA Applica

LS-DYNAOne code strategy “Combine the multi-physics capabilities into one scalable code for solving highlynonlinear transient problems to enable the solution of coupled multi-physics andmulti-stage problems”Explicit/Implicit Heat Transfer ALE & Mesh Free i.e., EFG, SPH, Airbag Par

Accommodates coupled simula

One Code Strategy One Model LS-­‐DYNA Mul

One code strategy Specialized codes for each problem Mul

One Model for All Applica

Mul

Mul

One code strategy: mul

One code strategy: mul

One code strategy: mul

Adap

New single user license • Node locked SMP Window license for single user O/S to allow usage of 16 processor cores • 40% price reduc=on compared to 16 cores with network license – 16 one core simultaneous jobs, – 8 two core simultaneous jobs, – 4 four core simultaneous jobs, – etc. • Extension of security sojware to single user Window’s O/S to license MPP version is underway to take advantage of beler scaling – The SMP version does not scale well ajer 6 to 8 cores.

Development Speakers LS-­‐PrePost Dummies Incompressible CFD Electromagne

LS-DYNA ® R7.0 and R611Release Highlights

*AIRBAG_PARTICLE • You can now use – heat convec=on coefficient HCONV and – fabric thermal conduc

*AIRBAG_PARTICLE • Nega

*AIRBAG_PARTICLE • *DATABASE_CPM_SENSOR • Spherical, rectangular and cylindrical shapes to capture pressure and mass flow • When using IAIR=2 mass of ini

AIRBAG Par

CONTACT • *CONTACT_CONTACT/Op

CONTACT • Smother responses in TIEBREAK contacts with OFFSET – Interface displacements are now incrementally computed rather than total displacements elimina

Improved PentaHedron Element • Implemented Element type 115 single point penta element with HG control • Robustness enhanced when pentahedron elements (depicted in brown) are run with element type 115 compared to element type 15 • CONTROL_SOLID automa

Cohesive Elements • Allow degenerated hexahedrons (pentas) for cohesive solid elements (ELFORM=19, 20) that evolve from an extrusion of triangular shells. The input of nodes on the element cards for such a pentahedron is given by: N1,N2,N3,N3,N4,N5,N6,N6 • Fix for bug due to mass scaling due to ROFLG flag in MAT_138

Fix for MAT_SPOTWELD • Modify the behavior of spot weld assemblies failure used with the TRUE_T parameter on *MAT_SPOTWELD. • There was a failure to measure and reduce only the moment due to shear loading. • Scales the Moment value in Lap Shear and Coach Peel Spot Weld coupon simula

HAZ proper

EFG • New EFG solid type 43 (called Meshfree-­‐Enriched FEM, MEFEM) for both implicit and explicit. • This element formula

One-­‐Step Forming Simula

Framework for One-­‐Step Forming Crash Model Selected Parts LS-­‐DYNA One-­‐Step Forming Simula

Comparison TWB parts

One Step Forming • Two input variables for *CONTROL_FORMING_ONESTEP simula

DEFINE_TABLE_3DSimplified Table Definitionsε*DEFINE_TABLE_2D– Unlike the *DEFINE_TABLE keyword, a curve ID isspecified for each abscissa value defined in thetable.– The same curve ID to be referenced by multipletables, and the curves may be defined anywhere inthe input file.*DEFINE_TABLE_3D– A table ID is specified for each abscissa valuedefined for the 3d table

Consider a thermal material model. For each temperature, T, wehave a table of hardening curves of stress versus strain at 3 strainrates, i.e, σ = f ( , , T)ε ε*Define_table_3dσεT=1000, TID=300T=500, TID=200T=20, TID=100εT=10000.0 0.1 0.4 1.00.1 98 190 244 268T=5001.0 109 212 273 30010. 115 224 287 316εεεε0.0 0.1 0.4 1.00.1 130 253 325 3571.0T=20146 283 364 40010. 153 298 383 421εε0.0 0.1 0.4 1.00.1 σ = 162 σ = 316 σ = 406 σ = 4461.0 182 354 455 50010. 192 373 479 527

*DEFINE_TABLE_{2,3}D Exampleεε*DEFINE_TABLE_3D$ tbid2000$ temperature tbid20. 100500. 2001000. 300*DEFINE_TABLE_2D$ tbid100$ strain_rate lcid0.1 1011.0 10210.0 103*DEFINE_CURVE$ lcid101$ strain stress0.0 1621.0 446For eachtemperature, wespecify a table with3 strain ratesFor each strainrate, we specifya curve of σ vs ε

*ELEMENT_SHELL_COMPOSITE A way to define elements for a general composite shell part where the shellswithin the part can have an arbitrary number of layersThe material ID, thickness, and material angle are specified for thethickness integration points for each shell in the part– The number of composite layers and overall shell thickness canchange from shell to shell– The thickness of each shell is the summation of the integrationpoint thicknesses. The total number of integration points isarbitrary– Implementation works with all standard shell formulationsOnly one part ID is needed as apposed to multiple parts for each38different layup

Contact ID’s are added to INTFOR39

Part titles are added to D3PLOT40

*DATABASE_BINARY_D3PLOTCard 1 2 3 4 5 6 7 8Variable DT/CYCL LCDT/NR BEAM NPLTC PSETIDType F I I I IDefault - - - - -RemarksPSETID: All parts defined in the set will be removed from D3PLOT41

MAT 224 • Tabulated Johnson-­‐Cook Model • Yield is a func

Plas

LS-­‐PrePost

New Graphics Rendering in version 4.0 • Taken from a visualiza

New Rendering Performance • 5.65million elements (4.29m Shells, 1.36m solids, some beams, 1680 parts), 59 states • On HP Z800 8-­‐core, with Nvidia Quadro 6000,

New Rendering Performance • 10.65million elements (8.44m Shells, 2.21m solids, 5223 beams, 816 parts), 49 states • Spot weld beam was drawn as circle • On HP Z800 8-­‐core, with Nvidia Quadro 6000,

User group and Online Documenta

Other new features and improvements inLS-PrePost3.2/4.0

Fast Rendering in ver 4.0/4.1• Fast rendering is the default renderingmode for versions 4.0/4.1• If graphics hardware is not capable,“Normal Rendering” will be usedautomatically• To switch between “Fast rendering” and“Normal rendering” mode, enter cntl-Ltwice before loading the data• Rendering mode will be memorized andrecorded in the configuration file50

Batch Mode Operations• Use -nographics on the command line:– LsPrePost c=cmd.cfile -nographics• -nographics means no graphics windows popup, however, still requires graphics hardwareand software to perform operations• -nographics works well on local desktopmachine. On remote machine, it requirescompatible graphics hardware between localmachine and remote machine• Use runc command for truly no graphicsoperations (same as LS-PrePost 2.4 or older):– LsPrePost runc=cmd.cfile51

Cutting Plane for CPM (Particle)• A special new cuttingplane interfaces hasbeen developed forSPH, CPM (particle),DES, and CFDanalyses• Multiple planes can bedefined and visualized• Must set env variablelspp_developemode toyes or onMultiple planes definitionsClick this icon to activate theplane interface52

Cutting Plane for CPM (Particle)Fringe particle data oncutting planes with gridSmoothed fringe particle dataon multiple planes with grid53

Scripting Command Language• Scripting Command Language – is a C-like programminglanguage to be executed within LS-PrePost• Executes LS-PrePost commands• Allows “if then else”, for, and while loop operations• Provides API (Application Programming Interface) toextract model and result data from LS-PrePost Data base• Operations can be done on extracted data to form newdata. New data can be output to file or fringed on screen• Most suitable to perform same operations over differentpart of the model• Documentation and tutorial for Scripting CommandLanguage will be available on LSTC ftp site soon54

LS-DYNA long format Support• A special double precision version of LS-PrePost is built tosupport LS-DYNA long format• Long format uses 64 bits for all floating point and integerdata• Long format allows user ID to be as large as 15 digits• Floating point number will be saved in 15 digits• This special double precision version will double thememory requirement• This version will process regular keyword file as well as32bit d3plot files• Not recommend for normal use and post-processing dueto the large memory requirement55

Other Miscellaneous Features• Print picture in PDF format• Solid element splitting options with all 3 directions (2x2x2)or only one direction• General selection for element can choose sub types:• Support 3D Connexion device• Support LS-DYNA Random Vibration Analysis, ShockSpectrum Analsys.56

Batch mode Opera

LS-­‐PREPOST Features for LS-­‐980 • Support for Mul

CESE with stochas

Fuel Tank Fluid Surface shown by Levelset part. Levelset can be fringed with CFD variables, and with velocity vectors on the surface

New XYPLOT layout • New XY plot interface allows xy plot to be drawn to main graphics windows, or to a separate page with mul

Fringing by Script • In the fringe expression interface, use script (a programming code) instead of expression • Assign components to variables • User write the script (code) to perform whatever data manipula

DynFold Applica

DynFold Setup Process • Define Parameters: Define Project Step Name, Termina

• Spc_Birth_Death, BPMF(Box), S

Define Part Mo

ISO-Geometry Element• To Create iso-geometry element, go to Mesh->Nurbs->Create,• Current development on capability to modify theisogemetry element within LS-PrePost. Allows user torefine the no. of patches and modify the control points68

Helix finite element model Creation• Both Shell element orSolid element modelcan be generated• Go to Mesh->EleGen->Shell or Solid, SelectHelix and type inparamters• Warning message willbe given if it is notpossible to create thehelix model69

Helix finite element model Creation• Pitches for both the Tdirection and R directioncan be used• Size is the element sizein the spiral direction• No. of loops has to bean integer number70

*Airbag_shell_reference_geometry • *Airbag_shell_reference_geometry is the required data for airbag deployment in LS-­‐DYNA • LS-­‐Prepost creates this data by asking user to pick the parts that make up the airbag in 3D final configura

*Airbag_shell_reference_geometry Pick this part to beunrolled*Airbag_shell_ref_geometry

Part Replace • Model-­‐>PartD-­‐>Replace • To replace a part with another part • The 2 parts do not need to be the same in no. of elements/nodes. • Connec

Part Replacement Old partNew partBeams are connected properlyautomatically

User wrilen script • C-­‐like programming language script to execute LS-­‐PrePost commands • Allows “if then else”, for, and while loop opera

Solid Meshing with Hex Element • Solid meshing by blocks -­‐ using cut and dice method and then sweeping

Metal Forming -­‐ Die System Module Complete metal forming Die design system

Metal Forming -­‐ Die System Module • Provides a user friendly interface to design the complete tooling system – Star

THUMS Posi

Summary • New GUI provides beler look and feel also yields maximum windows space for graphics, at the same

LS-­‐PrePost Recap LSTC is commiled to con

Other Miscellaneous Improvements • Many bugs have been fixed in geometry engine • Improved mid-­‐surface genera

Dummies Update

Released LSTC Dummy Models Detailed ModelsHYBRID III 5 thHYBRID III 50 thHYBRID III 95 th (scaled)SID IIs DEuroSID 2EuroSID 2reUSSIDFAST ModelsHYBRID III 5 thHYBRID III 50 thHYBRID III 95 thSID IIs DHYBRID III 5 th Lower BodyHYBRID III 50 th Lower BodyHYBRID III 50 th standingHYBRID III 6-year-oldFree Motion HeadformPedestrian LegformsBioRID II (ALPHA)84

LSTC Dummy Models in Development ModelHYBRID III 3-year-oldStatusMaterial OptimizationHYBRID III 95 thHYBRID III 95 th FASTBioRID IIWorldSID 50 thModel Improvements and MaterialOptimizationModel Calibration and SledVerificationModel Improvements and MaterialOptimizationModel Build-upTHOR NTMeshingEjection Mitigation HeadformMaterial OptimizationHYBRID IIMeshing85

Planned LSTC Dummy Models • Pedestrian Headforms • FAST versions of EuroSID 2 and EuroSID 2re • Q-­‐series child dummies • Flex PLI • WorldSID 5 th percen

LSTC Dummy Models Recap We commiled to the con

LSTC Barrier Models UpdateLSTC Barrier Models Update88

LSTC family of barriersEEVC WG13 /Japan214 sideIIHS sideEuro. sideEuro NCAPAE-­‐MDB V3.1 Shell V1.0 214 Shell V2.0 Solid V3.0 ODB Shell V2.0 Hybrid V2.0 IIHS Shell V1.0 Solid V3.0 MDB, ECE-­‐R95 Shell V2.0 89

LSTC barrier modelsODB 2.0Version Shell Solid214 3.0IIHS 3.0MDB, ECER95 2.0AE-MDB, EEVC WG13 3.1RCARBetaContact dilip@lstc.com for more information90

Thank You !

Incompressive CFD

Introduc

Coupling with other LS-­‐DYNA solvers • Scope of the new 980 solvers to be coupled with LS-­‐DYNA solvers in order to solve complex fluid-­‐structure or thermal problems, • Strong coupling is available for implicit mechanics. More robust but more costly, • Loose coupling for explicit mechanics. Less robust and less costly. Suitable for simpler couplings. • E.g. aeroelas

Applica

Applica

Sloshing Water Tank example : Moving Water Tank coming to a brutal halt, Sloshing occurring, Study of pendulum oscilla

Strong FSI Coupling High viscosity Liquid coming out of bolle due to finger pressure Simula

Conjugate Heat Transfer Monolithic strong coupling between the solid and fluid thermal solvers providing good stabilityThermal pipe flow exampleSolid Mesh Fluid veloci=es and Conjugate heat analysis (cut plane of fluid vel.) Steady state fluid temperature (cut plane through the sec=on)

Range of Applica

Automa

Error Control and Adap

MPP Scalability: Real Car Model The results show a speedup of 40 for 128 cpus in the CFD only case (2.1 elements) and a speedup of 55 for 128 cpus in the case of FSI (3.6 elements). For the next development cycle further improvements will be implemented

Incompressible CFD Roadmap Validation/Benchmarking process under wayproblems will include FSI, Conjugate Heat transfer cases as well as moreAerodynamics and Free surface analyzesAdditional post treatments and Tools with LSPP 3.2.See Website for additional documentationhttp://www.lstc.com/applications/icfd

Thank You !

Electromagne

Presenta

Coupling with other LS-­‐DYNA Solvers Scope of the new 980 solvers : to be coupled with LS-DYNA solvers in order to solvecomplex multi-physics problems

Electromagne

Electromagne

Electromagne

Electromagne

Electromagne

Electromagne

Other Possible Applica

Advancement Status • All EM solvers work on solid elements (hexahedral, tetrahedral, wedges) for conductors. • Shells can be used for insulator materials. • Serial and MPP versions available. • 2D axi-­‐symmetric available. • The EM fields as well as EM force and Joule hea

Plan for Future Introduction of Magnetic materials.Further optimization of the FEM / BEM calculations.Continue the validation process (T.E.A.M. problems).Wishes from users. Please let us know !

Thank you for your Alen

SPH, ALE, DEM, Airbag Par

SPH Thermal Solver • An explicit thermal conduc

Metal Cu|ng with Heat Ini

Metal Cu|ng with Heat Heat source: *BOUNDARY_FLUX *MAT_JOHNSON_COOK (stress flow depends on the temperature) Temperature Plas

Fric

Fric

ALE and Thermal Coupling ALE *MAT_GAS_MIXTURE coupled with shell structure using *CONSTRAINED_LAGRANGE_IN_SOLID X=132cmX=0cmGas mixture*DATABASE_PROFILE Energy is removed from gas and deposited to shell via heat convec

ALE Dynamic Adap

Dynamic Adap

Dynamic Adap

Par

Adiaba

Adiaba

Discrete Element Sphere (DES) Dry Wet

LSTC DES Bond Model • All par