Short Introduction to LS-DYNA and LS-PrePost

Jimmy Forsberg
Content

■ DYNAmore Nordic presentation

■ Introduction to LS-DYNA
■ General work with different solvers.
■ LS-DYNA capabilities
■ Keywordfile structure

■ Introduction to LS-PrePost
■ Layout
■ Pre-processing
■ Post-processing
■ Special features
■ Composite tool

Test 2013-09-09 4
DYNAmore Group

■ CAE Software
■ Engineering services
■ Distributor for LSTC
■ Personnel: 70
■ LSTC code developer: 10
■ Head office in
Stuttgart, Germany

5
DYNAmore Group

■ Sweden
■ 17 Employees
■ 37 years in average
■ 9 Ph.D.
■ 8 M.Sc.
■ 1 Economics/Adm
■ Office in Linköping
■ Office in Göteborg

6
DYNAmore GmbH

Germany
■ ~60 Employees
■ Headquarters in Stuttgart-Vaihingen
■ Offices
■ Ingolstadt
■ Dresden
■ Wolfsburg
■ Fürstenwalde (Berlin)
■ On-site Offices
■ Sindelfingen
■ Untertürkheim
■ Weissach
■ Ingolstadt

Stuttgart [Headquarters]

7
Business model

Technical Software Software development Consultancy work

• Sales • Development • Non-linear analysis

• Support • Research • Linear analysis
• Training • Implementation • Dynamic analysis
• Improvement • Static analysis
• Support • Optimization

DYNAmore Nordic AB
• LS-DYNA • Vehicle safety
• LS-OPT • Material modeling • Explosion analysis
• Ansa • Contacts • Metal forming
• Crash dummies • Element technology • Offshore
• Crash barriers • Training • Energy
• Oasys Primer • GUI development • Roadside safety
• DynaForm • HPC Cluster • Accident reconstruction
• FormingSuite • Vibration and NVH
• Femzip • Thermo-mechanical
• On-site

8
DYNAmore Nordic - Selected customers

9
DYNAmore Group – Selected customers

10
Contact
■ Software Products
■ Dr. Marcus Redhe
■ E-mail: [email protected]
■ Mobile: +46 – (0)70 55 131 42
■ Engineering Service and Support
■ Dr. Daniel Hilding
■ E-mail: [email protected]
■ Mobile: +46 – (0)70 65 366 85
■ Address:
DYNAmore Nordic
Brigadgatan 14
587 58 Linköping
Sweden
■ Web: http://www.dynamore.se
■ Phone: +46 – (0)13 23 66 80

11
Introduction to LS-DYNA

Test 2013-09-09 12
 LS-DYNA One code strategy
“Combine the multi-physics capabilities into one scalable code for solving
highly nonlinear transient problems to enable the solution of coupled multi-
physics and multi-stage problems”

Explicit/Implicit

Heat Transfer Incompressible Fluids R7

Mesh Free CESE Compressible Fluid R7

EFG,SPH,Airbag Particle Solver

User Interface Electromagnetism R7

Elements, Materials, Loads

Acoustics Frequency
Response, Modal Methods

Discrete Element Method

SBD – Simulation Based Design
■ Instead of a physical prototype, a virtual
model is created. The purpose of the
model is to resemble the behaviour of the
physical product.
■ All development/testing is made in the
virtual product. Thus, you treat the
model as you would if it was a physical
product.
■ The benefits are several:
■ Shorter time to market
■ Reduce number of costly prototypes
■ Increased innovation
■ Lower development costs
■ Higher quality
■ … but also the challenges
■ Rethink development process
■ Trust the results
■ Educate personnel, new partners..
Volvo XC60

15
What do you need?

History?
PRE-PROCESSOR
Generates the FE-model Geometry
Applies boundary conditions etc Material
Process
LS-PrePost
SOLVER
Solves the numerical model

LS-DYNA

POST-PROCESSOR
View the results

LS-PrePost Dependence on analysis

17
Simulation process
Build FE-model
-Parts
-Material LS-DYNA
-Element

LS-PrePost/ANSA
Modify
-Process
- Initial powder volume Pre – simulation?
-Geometry -Initial stress/stress?
- Bolts etc.
LS-PrePost
LS-PrePost
Evaluate results

LS-PrePost
LS-DYNA

Test 2013-09-09 18
Introduction to LS-DYNA

Test 2013-09-09 19
Keywords and Elements
 Keywords - Define Geometry
Input file (.k)

Newton’s second law, F=ma, requires consistent units

S1 S2 S3
length meter millimeter millimeter
time second second millisecond
mass kilogram tonne kilogram
force Newton Newton kiloNewton
Young’s modulus of steel 210.0E+09 210.0E+03 210.0
density of steel 7.85E+03 7.85E-09 7.85E-06
gravitation 9.81 9.81E+03 9.81E-03

Keyword and Elements 2013-09-09 21

Keyword User’s manual

Keyword and Elements 2013-09-09 22

 Input file - Keywords
Input file (.k) Element PART
*KEYWORD Mandatory
*TITLE
Test example
$ Control cards govern entire model / simulation Section Material
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
$ Define output of results
*DATABASE_BINARY_D3PLOT
*DATABASE_GLSTAT
Comment card $ Define section and material
*PART
begins with $ $ Define element types and integration
*SECTION_SHELL
$ Define material properties
*MAT_ELASTIC
*MAT_FIBER
$ Define nodes and elements
*NODE
*ELEMENT_SHELL
$ Define loads and BC
*LOAD_NODE
*END
Mandatory

Keyword and Elements 2013-09-09 23

 Keyword Format
Input file (.k)
■ Similar functions are grouped together under
the same keyword
■ A data block begins with a keyword and ends
with the next keyword
■ Keywords are left justified
■ No distinction between lower and upper case
letters
■ Variables are right justified in their fields
■ A ‘0’ or blank means that the variable will get
the default value
■ The decimal point is always written out for
floating point variables

Keyword and Elements 2013-09-09 24

 Keyword Format
Input file (.k)
■ Comments rows are written after a dollar sign in
the first position
■ *COMMENT keyword exist
■ Do not use ‘tabs’ when editing or creating your
file
■ Line feed signs may cause problems when
transferring files from Dos to Unix

Keyword and Elements 2013-09-09 25

 Keywords - Define Geometry
Input file (.k) Element PART

Section Material
*NODE
y $ NID * x * y * z *

Fixed format
1 0.00 0.00 0.00
0.02 2 1.0E-2 0 0
3 0.02 0 0
n4 n5 n6 4 2.0E-2 1.0E-2 0
0.01
e1 e2 5, 0.01, 0.01, 0.0
0 6, 0.02, 0.01, 0 Free format
n1 n2 n3 $
$
x *ELEMENT_SHELL
$ID, PID, n1, n2, n3, n4
0 0.01 0.02 1, 1, 1, 2, 5, 4
2, 1, 2, 3, 6, 5
4 5 Local coordinate system:
ye1 Xe1: from n1 to n2
e1
ye1: perpendicular to xe1
xe1
directed towards n3
1 2

Keyword and Elements 2013-09-09 26

Elements

■ Some element formulations are more

costly than others

■ Stresses and strains are calculated at the

integration points
■ Accelerations, velocities and
displacements are evaluated at the nodes

Keyword and Elements 2013-09-09 27

Under Integrated and Fully Integrated Elements

■ Most element formulations in LS-DYNA are under-

integrated, i.e. the stresses and strains are only
calculated in the mid-point of each element.
■ Advantage: Computational efficiency. The material model
is called once per integration point and time step.
■ Disadvantage: The element formulation contains zero-
energy modes (hourglass modes)

Integration
point (s)

Keyword and Elements 2013-09-09 28

Under Integrated and Fully Integrated Elements

■ The following element deformation does not yield any

strains in the integration point, and thus no stress
■ There is deformation, but no associated internal energy,
hence the name zero-energy modes.
■ These modes have to be suppressed using ”hourglass
control”

 

“No strain”
x  x 

Keyword and Elements 2013-09-09 29

Hourglass Control

■ Zero energy modes = Hourglass modes

■ Hourglass controlled by *CONTROL_HOURGLASS and
*HOURGLASS
■ Hourglass modes for 1 point integration Q4 shell
elements:

■ Hourglass modes for 1 point integration solid elements:

+ 8 more!

Keyword and Elements 2013-09-09 30

SECTION_SHELL

■ Element formulation
■ Belytschko-Tsay
■ Belytschko-Wong-Chiang
■ Hughes-Liu
■ Belytschko-Leviathan
■ Fully integrated shells
■ Higher order shells 6/8 noded tria/quad
■ ……
■ Element thickness
■ Number of integration points through shell thickness

Keyword and Elements 2013-09-09 31

Elements (shell) - NIP

■ 1 point integration through the thickness gives a

membrane element
■ 2 point integration through the thickness is the default
(sufficient for a linearly elastic material)
■ For plastic bending behaviour, at least 3 points are
needed through the thickness
■ 5 points recommended for sheet metal stamping.
7 points for springback
■ Use odd numbers to include the neutral axis

Keyword and Elements 2013-09-09 32

 0
5
10
15
20
25

1
Belytschko-Tsay

Belytschko-Tsay +

1.07
warping stiffness
Element Performance

Belytschko- 1.25
Leviathan
Belytschko-Wong-
1.28

Chiang

Hughes-Liu
1.49

Hughes-Liu
2.45

(corotational)

Keyword and Elements

Fully integrated BT
2.8

(type 16)

Fully integrated HL
8.8

Fully integrated HL
20

(corotational)
2013-09-09
33
*CONTROL_ACCURACY

■ Invariant node numbering

■ particularly important when large shear forces are present in an
element
■ 2nd order stress update
■ spinning bodies such as turbine blades, rotating tires
■ sometimes for stiffness hourglass control
■ implicit solutions with large strains in each step

Keyword and Elements 2013-09-09 34

Material Models
Material Models

■ Over 200 models for various applications exists in

LS-DYNA.
■ Determine the stress based on strain, strain-rate,
temp etc.
■ Not materials, but models subject to restrictions:
■ Load magnitude
■ Deformation speed (strain rate)
■ Temperature
■ The models are defined by material parameters
■ E, , , etc.

Material Models 2013-09-09 36

Hypoelasticity

Hypoelasticity relates a strain rate to a corresponding stress

rate
σ  C : D Hooke’s law:
  E
Stress is incrementally updated from the strain rate with
aid of the constitutive tensor C
  t
  Dt

Most of the materials in LS-DYNA are based on this

formulation for the elastic response.

Material Models 2013-09-09 37

Merits and drawbacks (theoretical)

+
■ It is fairly straightforward to use and easy to implement
in a finite element code
-
■ The response is path-dependent, the stress for a closed
strain cycle can be nonzero, it should be used when the
elastic deformation is relatively small
■ It is difficult to deal with anisotropic constitutive models
because the constitutive tensor C is restricted to be
isotropic for nonlinear analysis. This is however solved in
LS-DYNA with a co-rotational update.

Material Models 2013-09-09 38

Hyperelasticity - definition

■ A material is hyperelastic if the internal work is

independent of the deformation path.
■ It is characterized by the existence of a strain energy
function that is a potential for the stress.

 (C) W (E)
S2 
C E
S Second Piola Kirchhoff stress tensor
E Green strain tensor
C Right Cauchy - Green tensor
■ Typically used when elastic deformation is substantial,
e.g. rubber.

Material Models 2013-09-09 39

Stress and strain – uni-axial deformation

Tensile test:
Engineering stress E  F / A 0 L
L0
Engineering strain  E  (L  L0 ) / L0
F F

A0 A
In LS-DYNA:
True stress   F/ A
True strain   ln( L / L0 ) σ

Elastic response:
Hooks law: σ  Eε
E
Area reduction: A  A0 (1  2)

ε
Material Models 2013-09-09 40
Elasto-plasticity in 3-D – multi-axial deformation

Deviatoric stress
Volumetric stress
   
Stress decomposition……….: ij  Sij  ij kk / 3
Von Mises yield criterion….: f  3
2 sij sij   y ( p )
Plastic strain…………………: f
ijp  
Sij

1 1 1

ε ijp ε ijp ε ijp

2 3 2 3 2 3

Perfect plasticity Isotropic hardening Kinematic hardening

Material Models 2013-09-09 41
Elastic-visco-plastic material

*MAT_PIECEWISE_LINEAR_PLASTIC
MID RO E PR SIGY ETAN FAIL TDEL
C P LCSS LCSR VP
EPS1 EPS2 .... ..... .
EP1 ES2 .... ..... .
For: Metals, loading exceeding yielding stress, rate effects
In: All element types
Theory: Isotropic plasticity model with visco-plasticity option
E Young's Modulus C,P Strain-rate parameters
RO Density LCSS Load curve for
PR Poisson's Ratio LCSR Load curve for strain-rate scaling
SIGY Yield stress VP Visco-plastic flag
ETAN Tangent modulus EPS1… Piecewise linear def.

Material Models 2013-09-09 42

Elastic-visco-plastic material

Activating visco-plasticity:
(  y =static yield stress )
s

C, P  0  No visco-plastic effects

 s 1
C, P  0, VP  1  Scale y by:  e  P
1    , eij  ij  ij kk
 C
1

C, P  0, VP  0  Scale 
s
   P 
y by: 1    ,   ijij
C 
   p 
1

 
P

C, P  0, VP  1  Yield stress is given by:  y   ys 1   eff 

 
  C  

VP=1 is recommended as it uses a consistent

visco-plastic theory
Material Models 2013-09-09 43
Elastic-plastic material with Bauschinger efftect

*MAT_PLASTIC_KINEMATIC
MID RO E PR SIGY ETAN BETA
SRC SRP FS VP
For: Metals under large loading
In: All element types
Theory: Isotropic and kinematic hardening plasticity, viscoplastic
E Young's Modulus
RO Density
PR Poisson's Ratio
SIGY Yield stress
ETAN Tangent modulus
BETA Hardening parameter FS Failure strain
SRC Strain rate parameter C VP Rate formulation flag
SRP Strain rate parameter P

Material Models 2013-09-09 44

Elastic-plastic material with Bauschinger efftect

Definition of material hardening:


 y0
ETAN

2y1
E 2y0
2y1

Kinematic hardening 0

Isotropic hardening  1

Other models with kinematic hardening:

*MAT_PLASTIC_GREEN-NAGHDI_RATE
*MAT_ANISOTROPIC_VISCOPLASTIC

Material Models 2013-09-09 45

*EOS

■ Certain material models only solve for the deviatoric part

of the stress tensor
■ An Equation of State (EOS) is required to find the
pressure part of the stress tensor
■ Mostly used in conjunction with fluid-like behaviour (high
explosives, airbag inflation …)
■ Solid elements only

Material Models 2013-09-09 46

Boundary/Initial Conditions
Initial and Boundary Conditions

■ Variation in time using load curves

■ Variation in space Traction( t )
■ Arbitrary directions using
■ Local coordinate systems
■ Vectors
But limited to cartesian coordinates
b( t )

u(t)

Initial/Boundary Conditions 2013-09-09 48

*LOAD

*LOAD_NODE[_SET|_POINT]
NODE/NSID DOF LCID SF CID M1 M2 M3

Nodal loads for one node or a set of nodes

DOF Direction of load in current coordinate system
LCID Load curve ID for variation in time
SF Scale load curve amplitude
CID Define a local coordinate system
F
M1-M3 Follower force definition
Singularities at point loads may be a problem. F
M3
Multiple load cards are accumulated.

M1

M2

Initial/Boundary Conditions 2013-09-09 49

*INITIAL
*INITIAL_STRESS[_BEAM|_SHELL|_SOLID]
*INITIAL_STRAIN[_SHELL|_SOLID]

Initialise the state of stress and strain in elements

Normally used to carry results obtained in one simulation to

another.
- Multistage forming
- Forming -> Crash
Keyword data normally generated automatically by
preprocessors.

Initial/Boundary Conditions 2013-09-09 50

Kinematic Conditions

■ Prescribe motion in the model

■ *BOUNDARY: w.r.t cartesian coordinates

■ Fixed supports
■ Symmetric boundaries
■ *CONSTRAINED: internal definitions
■ Mechanical Joints
■ Merging shell-brick elements
■ Define rigid bodies

Initial/Boundary Conditions 2013-09-09 51

*BOUNDARY_PRESCRIBED_MOTION_**

*BOUNDARY_PRESCRIBED_MOTION[_NODE|_SET|_RIGID]
ID DOF VAD LCID SF VID DEATH BIRTH

Apply nodal displacement, velocity, or acceleration to the model;

translations or rotations

DOF Direction of load, global or local direction, see manual!

VAD Type of load
LCID Load curve ID for variation in time
SF Scale amplitude of the loadcurve
VID Vector ID for vector to be used if DOF=4 or 8
DEATH/BIRTH Active range of time for this boundary condition

Use the _RIGID option for rigid bodies.

For local directions with rigid bodies see the MAT_RIGID keyword.

Initial/Boundary Conditions 2013-09-09 52

*CONSTRAINED

*CONSTRAINED_NODAL_RIGID_BODY
PID CID NSID PNODE IPRT

Create a new rigid body using existing nodes

PID Part id req. is a unique one

CID Coordinate system for output
NSID Node set
PNODE Optional centre node
IPRT Print flag RB

For spot-welds and other types of rigid connections.

Initial/Boundary Conditions 2013-09-09 53

*CONSTRAINED

*CONSTRAINED_JOINT_”JOINTTYPE”
N1 N2 N3 N4 N5 N6 RPS DAMP

Define mechanical joints between rigid bodies

N1-N6 Nodes in the rigid bodies
RPS Scale the penalty stiffness
DAMP Dynamic damping

N1,N3,N6 in RB1. N2,N4,N6 in RB2.

Place the nodes in one RB far apart.
(N1,N2) etc. initially coincident, except universal joint, read the
manual!
Motor and gear joints are available for advanced mechanisms.

Initial/Boundary Conditions 2013-09-09 54

*CONSTRAINED

Cylindrical Revolute Planar

Spherical Translational Locking Universal

Initial/Boundary Conditions 2013-09-09 55

Contacts
Some of the available contacts *CONTACT_option_option_…

AIRBAG_SINGLE_SURFACE NODES_TO_SURFACE_INTERFERENCE
AUTOMATIC_GENERAL ONE_WAY_SURFACE_TO_SURFACE
AUTOMATIC_GENERAL_INTERIOR RIGID_NODES_TO_RIGID_BODY
AUTOMATIC_NODES_TO_SURFACE RIGID_BODY_ ONE_WAY_TO_RIGID_BODY
AUTOMATIC_NODES_TO_SURFACE_TIEBREAK RIGID_BODY_TWO_WAY_TO_RIGID_BODY
AUTOMATIC_ONE_WAY_SURFACE_TO_SURFACE SINGLE_EDGE
AUTOMATIC_SINGLE_SURFACE SINGLE_SURFACE
AUTOMATIC_SURFACE_TO_SURFACE SLIDING_ONLY
AUTOMATIC_SURFACE_TO_SURFACE_TIEBREAK SLIDING_ONLY_PENALTY
CONSTRAINT_NODES_TO_SURFACE SURFACE_TO_SURFACE
CONSTRAINT_SURFACE_TO_SURFACE SURFACE _TO_SURFACE_INTERFERENCE
DRAWBEAD TIEBREAK_NODES_TO_SURFACE
ERODING_NODES_TO_SURFACE TIEBREAK_ SURFACE _TO_SURFACE
ERODING_SURFACE_TO_SURFACE TIED_NODES_TO_SURFACE
FORCE_TRANSDUCER_CONSTRAINT TIED_NODES_TO_SURFACE_OFFSET
FORCE_TRANSDUCER_PENALTY TIED_SHELL_EDGE_TO_SURFACE
FORMING_NODES_TO_SURFACE_TIEBREAK SPOTWEALD
FORMING _ONE_WAY_SURFACE_TO_SURFACE SPOTWEALD_WITH_TORSION
FORMING _SURFACE_TO_SURFACE TIED_ SURFACE _TO_SURFACE
NODES_TO_SURFACE TIED_ SURFACE _TO_SURFACE_OFFSET

Contacts 2013-09-09 57
Contact
■ A way of treating interaction between different parts
■ Contacts are defined by sets (node/part/segments) or a
box
■ Generally there is a master side and a slave side of the
contact
■ The master side can be a mathematically described with
a geometrical surface (rigid)
■ The thickness of shells are normally taken into account
■ Most recommended contacts are based on the penalty
method
■ Several contacts treating special applications exists
■ Old contact types kept for
compatibility reasons Motion

SLAVE
MASTER

Contacts 2013-09-09 58
Interesting Keywords for Contacts
■ Contacts in LS-DYNA is affected by many different
keywords
■ *SECTION_SHELL (Shell thicknesses,
middle/top/bottom surface meshed)
■ *MAT_xxx (Penalty stiffness, E, pr, dens)
■ *DEFINE_FRICTION (Friction behavior between parts)
■ *PART_CONTACT (Contact behavior for parts)
■ *CONTROL_CONTACT (Overall contact behavior)
■ *CONTACT_xxx (Contact definition)

■ The different parameters on different keywords slave

might be used, depending on contact type.
■ The parameter might have a different meaning
depending on contact type use.
■ Makes contact definition tricky in LS-DYNA!
Some of the most interesting parameters found on master
different cards will be examined in this
presentation.

Contacts 2013-09-09 59
Important contact parameters: penalty method (default)

Motion Contact force

Fi= δi k
k= interface spring stiffness
Solid elements Shell elements

cKA cKA
k k
V diagonal
A
K= bulk modulus
c = penalty factor
δi=
penetration

The time step of the analysis is determined by LS-DYNA from the

elements of the FE-mesh without considering the contact interfaces!

Contacts 2013-09-09 60
Contact Thickness and Initial Penetrations

d1

d2
Initial
penetration

Change of shell thickness

only for contact treatment

d1‘

d2‘

Contacts 2013-09-09 61
Important contact parameters: friction
■ Sliding friction – FS, FD, DC and VC
■ Defined in keyword *CONTACT
■ Based on Coulomb friction
■ Default values gives no friction
 DC Vrel
c  FD  ( FS  FD)e
■ FS and FD are static respectively dynamic friction coefficient
■ DC - decay coefficient
■ If FD and FS not are equal, then FD should be less than FS
and DC nonzero
■ VC is the coefficient for viscous friction and limits the friction
force (typically 3-½ of yield stress)
■ Viscous damping VDC improves stability. For metal contacts use
20% and for soft material 40-60%

Contacts 2013-09-09 62
Automatic contacts without self contact
■ *CONTACT_AUTOMATIC_NODES_TO_SURFACE
■ *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE

■ Thickness taken into account

■ Contact surface is offset by half thickness from mid-plane
■ Orientation of segments not needed
■ Contact from both sides
■ Handles disjoint meshes
■ Applies a smooth surface based on a
radii at the edges (including free edges)
■ Initial penetrations are detected
■ Possible to change or scale contact thickness
■ Friction and damping available

Contacts 2013-09-09 63
Single-surface contacts (self contact)
■ *CONTACT_AUTOMATIC_SINGLE_SURFACE
■ *CONTACT_AUTOMATIC_GENERAL

■ Same features as the

automatic contacts
■ Only require definition of
the slave surface
■ Include self contact
■ Sensitive to initial penetrations
■ Possible to use only one contact
definition forthe complete model
■ Beam and edge to edge contacts
are included
*CONTACT_AUTOMATIC_GENERAL

Contacts 2013-09-09 64
Edge/Beam Contacts
■ *CONTACT_AUTOMATIC_GENERAL (26)
■ exclude interior edges
■ entire length of each exterior edge is checked for contact
■ OBS, the edge cylinder is not affected by OPTT or TH when d
using part_contact.
■ *CONTACT_AUTOMATIC_GENERAL_INTERIOR (i26)
■ like *CONTACT_AUTOMATIC_GENERAL,
■ but interior edges are treated like exterior edges d/2
■ Alternative way to treat edge contact:
■ creating null beam elements (*ELEMENT_BEAM,*MAT_NULL)
approximately 1mm in diameter along every edge wished to be
considered for edge-to-edge contact and including these null
beams in a separate AUTOMATIC_GENERAL contact
■ *CONTACT_SINGLE_EDGE (22)
■ Treats only edge-to-edge contact
■ no thickness offset at the contact edge
■ *CONTACT_xxx_MORTAR ()
■ edge-to-edge contact
■ no thickness offset at the contact edge

Contacts 2013-09-09 65
Tied contacts
■ CONTACT_TIED_NODES_TO_SURFACE
■ *CONTACT_TIED_SURFACE_TO_SURFACE
■ *CONTACT_TIED_SHELL_EDGE_TO_SURFACE
…._OFFSET

■ Possibility to “tie” nodes to a surface (segment)

■ NODES_... and SURFACE_... ties translational d.o.f
■ SHELL_EDGE_.. ties translational and rotational d.o.f
■ Constraint based. Thus, will not work with rigid bodies.
■ …_OFFSET allows for a segment thickness and is penalty based.
■ …_TIEBREAK_... has failure options.
■ Can be used to model glue, spotwelds etc.

Contacts 2013-09-09 66
Control Cards & Execution
Control Cards

■ The purpose is:

■ Activate solution options;
implicit solution, adaptive remeshing, mass scaling …
■ Change default values on options and parameters

■ Remember that:
■ Ordering between them and position are arbitrary
Good practise is to put them first in your input file
■ Do not use more then one control card of each type
■ All control cards are optional except
*CONTROL_TERMINATION

Control Cards & Execution 2013-09-09 68

Control Card Default Values

■ Default values exist for all options and most parameters

■ Control cards change default values globally
■ Default values are defined hierarchically
The order between them are:
■ LS-DYNA defaults
■ Control card input
■ Individual Keyword input
■ Set your defaults with the control cards and change the
keyword input where default values not should be used
■ Input of ‘0’ will normally give the default value which is
shown in the manual

Control Cards & Execution 2013-09-09 69

Most Important Control Cards

■ Always consider the following control cards since

they can strongly affect your results or output

■ *CONTROL_ACCURACY
■ *CONTROL_CONTACT
■ *CONTROL_ENERGY
■ *CONTROL_HOURGLASS
■ *CONTROL_SHELL
■ *CONTROL_SOLID
■ *CONTROL_TERMINATION
■ *CONTROL_TIMESTEP

Control Cards & Execution 2013-09-09 70

Implicit Solution Types

■ Linear Analysis
● static or dynamic
● single, multi-step
■ Eigenvalue Analysis
● frequencies and mode shapes
● linear buckling loads and modes
● modal analysis: extraction and superposition
● Dynamic analysis by modal superposition (971)
■ Nonlinear Analysis
● Newton, Quasi-Newton, Arclength solution
● static or dynamic

■ default LS-DYNA: static and nonlinear!

Control Cards & Execution 2013-09-09 71

Output Files

■ Binary files (can be viewed in LS-PrePost)

*DATABASE_BINARY_Option
■ ASCII files for more detailed output
(graphs can be shown in LS-PrePost)
*DATABASE_Option
■ Data in the binary and ASCII files is controlled by
*DATABASE_EXTENT_Option
*DATABASE_HISTORY_Option
■ Control files (d3hsp)
■ Message files (messag)

Control Cards & Execution 2013-09-09 72

Output Files

■ D3PLOT (database for complete output states)

■ D3DUMP (complete database for restart)
■ RUNRSF (running restart file, overwritten)
■ D3PART (as D3PLOT but includes just specified parts)
■ D3THDT (database for time history data of element
subsets)
■ D3DRLF (dynamic relaxation database)
■ D3MEAN (CFD database)
■ INTFOR (database for output of contact interface data)
■ XTFILE (extra time history data)
■ D3EIGV (modal data from eigenvalue analysis)
■ D3CRCK (crack data from Winfrith concrete model)

Control Cards & Execution 2013-09-09 73

ASCII Output Files
■ ABSTAT (airbag statistics)
■ GLSTAT (global data)
■ BNDOUT (boundary condition
■ MATSUM (material energies) force/ energy)
■ RCFORC (resultant interface ■ RBDOUT (rigid body data)
forces)
■ GCEOUT (geometric contact
■ SLEOUT (sliding interface energy) entities)
■ NODOUT (nodal point data) ■ JNTFORC (joint force)
■ ELOUT (element data) ■ SBTOUT (seat belt output)
■ SECFORC (cross section forces) ■ AVSFLT (AVS database)
■ RWFORC (rigid wall forces) ■ SWFORC (nodal constraint
■ SSSTAT (subsystem data) reaction forces)
■ DEFORC (discrete elements) ■ MOVIE
■ NCFORC (nodal interface forces) ■ MPGS
■ DEFGEO (deformed geometry) ■ TRHIST (trace particle
■ SPCFORC (SPC reaction forces) history)
■ NODFOR (nodal force groups) ■ TPRINT (thermal output)
■ SPHOUT (SPH data)

Control Cards & Execution 2013-09-09 74

Test 2013-09-09 75
Demonstrate LS-PrePost

■ PreProcessing
■ PostProcessing

Test 2013-09-09 76
Thank you!

Test 2013-09-09 77