Preload
Example showing the use of dynamic relaxation to preload a system with a force.
https://www.dynaexamples.com/show-cases/preload
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Preload
Example showing the use of dynamic relaxation to preload a system with a force.
*BOUNDARY_SPC_SET *CONTACT_AUTOMATIC_SINGLE_SURFACE_ID *CONTACT_FORCE_TRANSDUCER_PENALTY_ID *CONTROL_MPP_IO_NODUMP *CONTROL_TERMINATION *DATABASE_GLSTAT *DATABASE_MATSUM *DATABASE_RCFORC *DATABASE_SLEOUT *DATABASE_BINARY_D3PLOT *DATABASE_HISTORY_NODE *DEFINE_CURVE_TITLE *ELEMENT_SOLID *ELEMENT_SHELL *END *KEYWORD_ID *LOAD_RIGID_BODY *MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE *MAT_RIGID_TITLE *NODE *PART *SECTION_SHELL *SECTION_SOLID_TITLE *SET_NODE_LIST *SET_PART_LIST_TITLE *TITLE
*KEYWORD_ID PRELOAD *TITLE SHOW CASE: PRELOAD *CONTROL_MPP_IO_NODUMP *CONTROL_TERMINATION 25.0 0 0.0 0.0 0.0 0 *DATABASE_GLSTAT 3.0E-3 1 0 1 *DATABASE_MATSUM 3.0E-3 1 0 1 *DATABASE_RCFORC 3.0E-3 1 0 1 *DATABASE_SLEOUT 3.0E-3 1 0 1 *DATABASE_BINARY_D3PLOT 0.5 0 0 0 0 *DATABASE_HISTORY_NODE 3016 3124 4144 4198 4256 4310 4368 4422 *MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE steel 1 7.86E-6 210.0 0.3 0.25 0.0 0.0 0.0 0.0 0.0 0 0 0.0 0.0 4.0E-2 9.0E-2 0.17 0.26 0.37 0.4 0.0 0.25 0.3 0.34 0.38 0.43 0.52 0.6 0.0 *MAT_RIGID_TITLE rigid fixed 2 7.86E-6 210.0 0.3 0.0 0.0 0.0 1.0 7.0 7.0 0 0.0 0.0 0.0 0.0 0.0 *MAT_RIGID_TITLE rigid x free 3 7.86E-6 210.0 0.3 0.0 0.0 0.0 1.0 5.0 7.0 0 0.0 0.0 0.0 0.0 0.0 *MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE plastic 4 1.0E-6 1.0 0.3 0.1 0.25 0.0 0.0 0.0 0.0 0 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 *SECTION_SHELL 1 16 0.833 5 1.0 0.0 0 1 0.4 0.4 0.4 0.4 0.0 0.0 0.0 0 *SECTION_SOLID_TITLE solid 2 2 0 *PART D rigid fix 1 1 2 0 0 0 0 0 *PART Y rigid fix 2 1 2 0 0 0 0 0 *PART A rigid fix 3 1 2 0 0 0 0 0 *PART N rigid fix 4 1 2 0 0 0 0 0 *PART m deformable 5 1 1 0 0 0 0 0 *PART r rigid x free 7 1 3 0 0 0 0 0 *PART e rigid fix 8 1 2 0 0 0 0 0 *PART cylinder rigid fix 9 1 2 0 0 0 0 0 *PART o deformable solid 10 2 4 0 0 0 0 0 *NODE *ELEMENT_SHELL *ELEMENT_SOLID *DEFINE_CURVE_TITLE load during initialization 6 1 0.0 0.0 0.0 0.0 0 0.0 0.1 5.0000000 0.10000000 50.000000 0.10000000 *DEFINE_CURVE_TITLE load during transient simulation 7 0 0.0 0.0 0.0 0.0 0 0.0 0.10000000 5.0000000 0.10000000 20.000000 1.0000000 50.000000 1.0000000 *BOUNDARY_SPC_SET 1 0 1 1 1 0 0 0 *SET_NODE_LIST *LOAD_RIGID_BODY 7 1 6 -1.0 0 7 1 7 -1.0 0 *CONTACT_AUTOMATIC_SINGLE_SURFACE_ID 1contact m o r 2 0 2 0 0 0 1 1 0.2 0.0 0.0 0.0 20.0 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1 0.0 0 0.0 0.0 0 0 0 *SET_PART_LIST_TITLE parts in contact 2 0.0 0.0 0.0 0.0 5 7 10 *CONTACT_FORCE_TRANSDUCER_PENALTY_ID 2force between r o 7 10 3 3 0 0 1 1 0.0 0.0 0.0 0.0 0.0 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 *CONTACT_FORCE_TRANSDUCER_PENALTY_ID 3force between o m 10 5 3 3 0 0 1 1 0.0 0.0 0.0 0.0 0.0 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 *END
contact forces versus time:
vmises stress at time 0ms:
vmises stress at time 23ms:
contact forces versus time