3 Crashbox
3 Crashbox
LS-DYNA / LS-PrePost
Ex. 3. Crash test
Contents
1 Introduction ............................................................................. Error! Bookmark not defined.
1.1 Purpose ..........................................................................................................................2
1.2 Prerequisites...................................................................................................................2
1.3 Problem Description........................................................................................................2
1.4 Data files ........................................................................................................................2
2 Create the model ....................................................................................................................3
2.1 Mesh ..............................................................................................................................3
2.2 Boundary conditions .......................................................................................................6
2.3 Rigid plate ......................................................................................................................7
2.4 Material properties ..........................................................................................................8
2.5 Element properties..........................................................................................................9
2.6 Motion .......................................................................................................................... 10
2.7 Contact ......................................................................................................................... 11
2.8 Tied contact .................................................................................................................. 12
2.9 Termination time ........................................................................................................... 13
2.10 Output .......................................................................................................................... 13
2.11 Control cards ................................................................................................................ 13
2.12 Save ............................................................................................................................. 13
3 Run the simulation................................................................................................................ 14
4 Post processing .................................................................................................................... 14
5 Summary and comments ...................................................................................................... 15
5.1 Contact definition .......................................................................................................... 15
5.2 Contact parameters ...................................................................................................... 15
5.3 Contact recommendations ............................................................................................ 16
5.4 Time step...................................................................................................................... 16
6 Optional exercises ................................................................................................................ 17
7 Record of revisions............................................................................................................... 18
8 Copyright and Trademark Notice .......................................................................................... 18
9 Disclaimer ............................................................................................................................ 18
1 Introduction
1.1 Purpose
• Learn to use mesh tools.
• Get familiar with contact definitions.
• Understand how to affect the time step
1.2 Prerequisites
• Basic knowledge in the finite element method.
• Understand the steps in tutorial 1 - Getting Started.
Symmetry
50 mm 200 mm
**
*
• Select the nodes in ***, change Sym Plane to All Fix and click Apply
• The constraints to the left should now have been created.
• Click Done
***
Even though the part cannot be deformed, the Young’s modulus E and Poison’s ratio PR should be
set. This is because they are needed to compute the contact stiffness, if a contact definition is used.
The parameters CMO, CON1 and CON2 are used to constrain the center of mass of the rigid body
in translational and rotational degrees of freedom:
2.6 Motion
The rigid plate is to crush the crash box. Apply the motion of the plates as follows. First create a
curve defining the motion:
Now use the defined curve to apply a motion to the rigid plate:
Note: The options above imply that a constant velocity motion will be applied in the negative (negative
scale factor) x-direction for the rigid Plate.
2.7 Contact
By using contact definitions, the structural domains will interact instead of just passing through each
other. For self-contact, meaning that nodes and elements from the same part are able to contact
each other, a single surface type of contact in LS-DYNA has to be chosen:
Note: The slave nodes in SSID (crash box) will be tied to a master surface in MSID (for shell
elements) if the distance between them are smaller than 𝛿.
In this case:
𝛿 = 2.16
Thus, for the tying to occur, the distance between the slave node and the master surface must
therefore be smaller than 2.16 mm. Check this as follows:
2.10 Output
Set the output of d3plot data from the simulation:
• *CONTROL_ACCURACY: Set OSU = 1 and INN = 4 (2 would also work in this case)
• *CONTROL_SHELL: Set ESORT = 1.
2.12 Save
The model is now ready to be saved, use File > Save As > Save Keyword As. Chose a folder path
and name your file crashbox.k for example. Note that the folder path cannot contain any spaces.
Close LS-PrePost.
4 Post processing
Open the d3plot. Click Forward in the Animation Toolbar to see what happens.
Plot the contact energy, i.e. frictionally dissipated energy and penalty spring energy:
Also plot all energies from 1-Kinetic Energy to 9-External Work to get an idea of the size of the
energies and work involved.
The contact force applied from the interface springs are calculated as:
𝑓𝑠 = −𝑙𝑘𝑐 𝑛𝑐
𝑓𝑠 = 𝑎𝑝𝑝𝑙𝑖𝑒𝑑 𝑓𝑜𝑟𝑐𝑒
𝑙 = 𝑙𝑒𝑛𝑔𝑡ℎ 𝑡ℎ𝑎𝑡 𝑠𝑙𝑎𝑣𝑒 𝑛𝑜𝑑𝑒 ℎ𝑎𝑠 𝑝𝑒𝑛𝑒𝑡𝑟𝑎𝑡𝑒𝑑
𝑘𝑐 = 𝑐𝑜𝑛𝑡𝑎𝑐𝑡 𝑖𝑛𝑡𝑒𝑟𝑓𝑎𝑐𝑒 𝑠𝑡𝑖𝑓𝑓𝑛𝑒𝑠𝑠
𝑛𝑐 = 𝑛𝑜𝑟𝑚𝑎𝑙 𝑑𝑖𝑟𝑒𝑐𝑡𝑖𝑜𝑛 𝑖𝑛 𝑤ℎ𝑖𝑐ℎ 𝑡ℎ𝑒 𝑓𝑜𝑟𝑐𝑒 𝑖𝑠 𝑎𝑝𝑝𝑙𝑖𝑒𝑑
This contact interface stiffness 𝑘𝑐 will be calculated in different ways, depending on the used
formulation. Which formulation to be used are set by the SOFT parameter in the CONTACT keyword
• IGNORE: Ignore initial penetrations. “Initial” in this context refers to the first timestep that a
penetration is encountered. By setting IGNORE = 2, LS-DYNA allow initial penetrations to
exist by tracking the initial penetrations. However, penetration warning messages are printed
with the original coordinates and the recommended coordinates of each slave node given.
For explicit simulations, the time step is determined for all deformable elements in the model. The
expression varies depending on if which type of element that are used. The general idea is:
𝐿𝑒 √𝜌
∆𝑡𝑒 ≈
√𝐸
𝐿𝑒 = 𝐶ℎ𝑎𝑟𝑡𝑒𝑟𝑖𝑠𝑡𝑖𝑐 𝑙𝑒𝑛𝑔𝑡ℎ 𝑜𝑓 𝑎𝑛 𝑒𝑙𝑒𝑚𝑒𝑛𝑡
𝜌 = 𝐷𝑒𝑛𝑠𝑖𝑡𝑦
𝐸 = 𝑌𝑜𝑢𝑛𝑔 ′ 𝑠 𝑚𝑜𝑑𝑢𝑙𝑢𝑠
The critical time step is determined as the minimum of all calculated time steps:
∆𝑡𝑐 = min(∆𝑡𝑒 )
The actual time step ∆𝒕 is calculated by scaling the critical time step, default in LS-DYNA is
∆𝒕 = 𝟎. 𝟗∆𝒕𝒄 . The scale factor TSSFAC (0.9 as default) can be changed in CONTROL_TIMESTEP.
To get more detailed information about the critical time step, see section “Time Step Control” in the
Theory Manual.
By observing the time step calculation, one can see that by changing certain material properties, the
time step can be increased or decreased. A regular approach is mass scaling i.e. modify the mass
of the parts in the model. One simple way is to change the density of the material, which is not
recommended since it in some cases require a large increase in mass in order to get a noticeable
change in the time step. A more stable and more commonly used mass scaling method is applied in
Optional exercise no. 2.
6 Optional exercises
Optional exercises:
1. At which thickness of the plate will the crashbox NOT be tied to the plate? Use the
following equations.
2. Add the keyword CONTROL_TIMESTEP and set DT2MS = -5e-7. If the time step for an
element becomes smaller than |DT2MS|, mass will be added to those elements so the time
step becomes equal to |DT2MS|. The minimum time step allowed is then TSSFAC x
|DT2MS|. In this case, 0.9 x |-5e-7| = 4.5e-7.
3. While running the simulation, notice that the time step never becomes lower than
4.5e-7. Use glstat to see how much mass that have been added to the model, both
in mass unit and percentage. To obtain information about the added mass for the
different parts, add matsum in DATABASE_ASCII_option.
7 Record of revisions
Rev. no Release Author Comment
date
1.0 Robin Thuresson Release of document.
2.0 2019-03-20 Anders Jernberg Formatting, copyright, disclaimer
9 Disclaimer
By using this Tutorial, you hereby consent to this disclaimer and agree to its terms.
All the information in this Tutorial, comprised of this document and the accompanying simulation
models, is published in good faith and for general information purpose only. Neither DYNAmore
Nordic AB nor the authors make any warranties about the completeness, reliability and accuracy of
the information in this Tutorial. Any action you take upon the information you find in this Tutorial is
strictly at your own risk. Neither DYNAmore Nordic AB nor the authors will be liable for any losses
and/or damages in connection with the use of the Tutorial.