Bar I
This example shows the deformation of two equivalent bars made of plastic material. Initial conditions are specified by *INITIAL_VELOCITY keyword. The upper bar is calculated using Lagrangian discretization, the lower one is calculated using SPH method.
https://www.dynaexamples.com/sph/bar-i
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Bar I
This example shows the deformation of two equivalent bars made of plastic material. Initial conditions are specified by *INITIAL_VELOCITY keyword. The upper bar is calculated using Lagrangian discretization, the lower one is calculated using SPH method.
*BOUNDARY_SPC_SET *CONTROL_CONTACT *CONTROL_DAMPING *CONTROL_ENERGY *CONTROL_OUTPUT *CONTROL_SHELL *CONTROL_SPH *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_EXTENT_BINARY *DATABASE_GLSTAT *DATABASE_MATSUM *ELEMENT_SOLID *ELEMENT_SPH *END *INITIAL_VELOCITY *INITIAL_VELOCITY_NODE *KEYWORD *MAT_PLASTIC_KINEMATIC *NODE *PART *SECTION_SOLID *SECTION_SPH *SET_NODE_LIST_GENERATE *TITLE
*KEYWORD
$* LS-DYNA Keyword file: LSTC LS-PRE/POST Ver 1.0(Beta)- 31OCT2002(1553)
$* Total Number of Nodes = 10652,Largest node ID = 2002700
$* Total Number of Solids = 2500,Largest Solid ID = 2500
$* Total Number of Parts = 4,Largest Part ID = 102
$* Total Number of Mats = 1 Largest Mat ID = 1
*TITLE
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ NODES AND ELEMENTS
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
Untitled
*CONTROL_CONTACT
$* slsfac rwpnal islchk shlthk penopt thkchg orien enmass
0.0 0.0
$* usrstr usrfrc nsbcs interm xpene ssthk ecdt tiedprj
0 0 0 0 0.0
*CONTROL_DAMPING
$* nrcyck drtol drfctr drterm tssfdr irelal edttl idrflg
0 0.0 0.0 0.0 0.0 0 0.0
*CONTROL_ENERGY
$* hgen rwen slnten rylen
1 2 1 1
*CONTROL_OUTPUT
$* npopt neecho nrefup iaccop opifs ipnint ikedit iflush
0 0 0 0 0.0
*CONTROL_SHELL
$* wrpang esort irnxx istupd theory bwc miter proj
0.0 0 0 0 9
*CONTROL_SPH
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *CONTROL_SPH $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* ncbs boxid dt dim memory form start maxv
1 0 0.0 2 0 0 0.0 0.0
*CONTROL_TERMINATION
$* endtim endcyc dtmin endeng endmas
20.00000 0 0.0 0.0 0.0
*CONTROL_TIMESTEP
$* dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.500000 0 0.0 0.0
*DATABASE_GLSTAT
0.001000
*DATABASE_MATSUM
0.001000
*DATABASE_BINARY_D3PLOT
0.400000
*DATABASE_EXTENT_BINARY
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0
*ELEMENT_SOLID
$
$ ELEMENT CARDS FOR SOLID ELEMENTS
$
*ELEMENT_SPH
$
*NODE
$
*INITIAL_VELOCITY_NODE
$
*INITIAL_VELOCITY
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ INITIAL CONDITION
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
1
0.010000 0.0 0.0 0.0 0.0 0.0
2
-0.010000 0.0 0.0 0.0 0.0 0.0
101
0.010000 0.0 0.0 0.0 0.0 0.0
102
-0.010000 0.0 0.0 0.0 0.0 0.0
*BOUNDARY_SPC_SET
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
4 0 0 0 1 0 0 1
*PART
$* title
$* pid secid mid eosid hgid grav adpopt tmid
1 1 1
*SECTION_SOLID
$* secid elform aet
1 1
*MAT_PLASTIC_KINEMATIC
$* mid ro e pr sigy etan beta
1 2.700000 0.165600 0.330000 0.002900 0.0 0.0
$* src srp fs vp
0.0 0.0 0.0 0.0
*PART
$* title
$* pid secid mid eosid hgid grav adpopt tmid
2 1 1
*PART
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ SPH PART $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* title
material type #101
$* pid secid mid eosid hgid grav adpopt tmid
101 101 1
*SECTION_SPH
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *SECTION_SPH $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* secid cslh hmin hmax sphini death start
101 1.100000 0.0 0.0 0.0 0.0 0.0
*PART
$* title
material type #102
$* pid secid mid eosid hgid grav adpopt tmid
102 102 1
*SECTION_SPH
$* secid cslh hmin hmax sphini death start
102 1.100000 0.0 0.0 0.0 0.0 0.0
*SET_NODE_LIST_GENERATE
1 0.0 0.0 0.0 0.0
1 2652
*SET_NODE_LIST_GENERATE
2 0.0 0.0 0.0 0.0
2653 5304
*SET_NODE_LIST_GENERATE
4 0.0 0.0 0.0 0.0
1 5304
*SET_NODE_LIST_GENERATE
101 0.0 0.0 0.0 0.0
1000001 1002700
*SET_NODE_LIST_GENERATE
102 0.0 0.0 0.0 0.0
2000001 2002700
*END