Bar II
This example shows the impact of two equivalent bars made of plastic material against a planar rigid wall. The upper bar is calculated using Lagrangian discretization, the lower one is calculated using SPH method.
https://www.dynaexamples.com/sph/bar-ii
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Bar II
This example shows the impact of two equivalent bars made of plastic material against a planar rigid wall. The upper bar is calculated using Lagrangian discretization, the lower one is calculated using SPH method.
*CONTROL_CONTACT *CONTROL_CPU *CONTROL_DAMPING *CONTROL_ENERGY *CONTROL_OUTPUT *CONTROL_SHELL *CONTROL_SPH *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_EXTENT_BINARY *DATABASE_GLSTAT *DATABASE_MATSUM *DATABASE_RCFORC *DATABASE_RWFORC *DATABASE_SLEOUT *ELEMENT_SHELL *ELEMENT_SPH *END *INITIAL_VELOCITY *KEYWORD *MAT_PLASTIC_KINEMATIC *NODE *PART *RIGIDWALL_PLANAR *SECTION_SHELL *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 = 6081,Largest node ID = 1004800
$* Total Number of Shells = 1200,Largest Shell ID = 1200
$* Total Number of Parts = 2,Largest Part ID = 101
$* Total Number of Mats = 2 Largest Mat ID = 2
*TITLE
SPH / TP1 bis Impact barre 2D contre SW
*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_CPU
$* cputim
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
2 2 2 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
*CONTROL_SPH
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *CONTROL_SPH $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* ncbs boxid dt dim memory form start maxv
0 0 0.0 2 0 0 0.0 0.0
*CONTROL_TERMINATION
$* endtim endcyc dtmin endeng endmas
0.035000 0 0.0 0.0 0.0
*CONTROL_TIMESTEP
$* dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.900000 0 0.0 0.0
*DATABASE_GLSTAT
4.0000e-04
*DATABASE_MATSUM
4.0000e-04
*DATABASE_RCFORC
4.0000e-04
*DATABASE_RWFORC
4.0000e-04
*DATABASE_SLEOUT
4.0000e-04
*DATABASE_BINARY_D3PLOT
0.004000
*DATABASE_EXTENT_BINARY
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0
*ELEMENT_SHELL
$
*ELEMENT_SPH
$
*NODE
$
*INITIAL_VELOCITY
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ INITIAL CONDITION
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
0 0 0
-50.00000 0.0 0.0 0.0 0.0 0.0
*PART
$* title
$* pid secid mid eosid hgid grav adpopt tmid
1 1 1
*SECTION_SHELL
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *LAGRANGIAN PART $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* secid elform shrf nip propt qr/irid icomp setyp
1 13 0.0 0 0 0.0
$* t1 t2 t3 t4 nloc marea
1.000000 1.000000 1.000000 1.000000 0 0.0
*MAT_PLASTIC_KINEMATIC
$* mid ro e pr sigy etan beta
1 0.0027106.8900e+04 0.330000 117.0000 6.800000 1.000000
$* src srp fs vp
0.0 0.0 0.0 0.0
*PART
$* title
$* pid secid mid eosid hgid grav adpopt tmid
101 101 1
*SECTION_SPH
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *SPH PART $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* secid cslh hmin hmax sphini death start
101 1.300000 0.0 0.0 0.0 0.0 0.0
*MAT_PLASTIC_KINEMATIC
$* mid ro e pr sigy etan beta
2 0.0027106.8900e+04 0.330000 117.0000 6.800000 1.000000
$* src srp fs vp
0.0 0.0 0.0 0.0
*RIGIDWALL_PLANAR
0 0 0 0.0
0.0 0.0 0.0 1.000000 0.0 0.0 0.0 0.0
*SET_NODE_LIST_GENERATE
101 0.0 0.0 0.0 0.0
1000001 1004800
*END