Bar III
This models the impact of two equivalent bars against a shell structure. The left bar is discretized using Lagrangian solid elements, the right one is discretized using SPH method.
https://www.dynaexamples.com/sph/bar-iii
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Bar III
This models the impact of two equivalent bars against a shell structure. The left bar is discretized using Lagrangian solid elements, the right one is discretized using SPH method.
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_TITLE *CONTROL_CONTACT *CONTROL_CPU *CONTROL_DAMPING *CONTROL_ENERGY *CONTROL_OUTPUT *CONTROL_SHELL *CONTROL_SPH *CONTROL_STRUCTURED *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_EXTENT_BINARY *DATABASE_GLSTAT *DATABASE_MATSUM *DATABASE_SLEOUT *DEFINE_BOX *DEFINE_CURVE *ELEMENT_SHELL_THICKNESS *ELEMENT_SOLID *ELEMENT_SPH *END *EOS_GRUNEISEN *INITIAL_VELOCITY *KEYWORD *MAT_NULL *MAT_PIECEWISE_LINEAR_PLASTICITY *NODE *PART *SECTION_SHELL *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 = 10848,Largest node ID = 1002500
$* Total Number of Shells = 5000,Largest Shell ID = 5000
$* Total Number of Solids = 2500,Largest Solid ID = 2500
$* Total Number of Parts = 4,Largest Part ID = 101
$* Total Number of Mats = 2 Largest Mat ID = 2
*TITLE
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_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
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
20.00000 2 -2 1 11 2 1
*CONTROL_SPH
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *CONTROL_SPH $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* ncbs boxid dt dim memory form start maxv
0 0 0.0 0 0 0 0.0 0.0
*CONTROL_STRUCTURED
*CONTROL_TERMINATION
$* endtim endcyc dtmin endeng endmas
0.100000 0 0.0 0.0 0.0
*CONTROL_TIMESTEP
$* dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.0 0 0.0 0.0
*DATABASE_GLSTAT
1.0000e-04
*DATABASE_MATSUM
1.0000e-04
*DATABASE_SLEOUT
1.0000e-04
*DATABASE_BINARY_D3PLOT
0.001000
*DATABASE_EXTENT_BINARY
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_TITLE
$* cidtitle
1
$* ssid msid sstyp mstyp sboxid mboxid spr mpr
2 1 3 3
$* fs fd dc vc vdc penchk bt dt
0.0 0.0 0.0 0.0 0.0 0 0.0 0.0
$* sfs sfm sst mst sfst sfmt fsf vsf
0.0 0.0 0.400000 0.0 0.0 0.0 0.0 0.0
$* soft sofscl lcidab maxpar edge depth bsort frcfrq
1 0.100000 0 0.0 0.0
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_TITLE
$* cidtitle
2
$* ssid msid sstyp mstyp sboxid mboxid spr mpr
101 3 4 3
$* fs fd dc vc vdc penchk bt dt
0.0 0.0 0.0 0.0 0.0 0 0.0 0.0
$* sfs sfm sst mst sfst sfmt fsf vsf
0.0 0.0 0.400000 0.0 0.0 0.0 0.0 0.0
$* soft sofscl lcidab maxpar edge depth bsort frcfrq
1 0.100000 0 0.0 0.0
*SET_NODE_LIST_GENERATE
101 0.0 0.0 0.0 0.0
1000001 1002500
*DEFINE_BOX
1 -4.000000 4.000000 -50.00000 50.00000 0.001000 21.00000
*DEFINE_CURVE
2 0 1.000000 1.000000 0.0 0.0
0.0 1100.000000
0.01000000 1170.000000
0.04000000 1240.000000
0.06000000 1300.000000
0.09000000 1400.000000
0.11000000 1470.000000
0.14000000 1570.000000
0.15000001 1600.000000
*ELEMENT_SOLID
$
$ ELEMENT CARDS FOR SOLID ELEMENTS
$
*ELEMENT_SHELL_THICKNESS
$
$ ELEMENT CARDS FOR SHELL ELEMENTS
$
*ELEMENT_SPH
$
*NODE
$
*INITIAL_VELOCITY
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ INITIAL CONDITION
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
0 0 1
0.0 0.0 -50.00000 0.0 0.0 0.0
*PART
$* title
plaque sous SPH
$* pid secid mid eosid hgid grav adpopt tmid
1 1 1
*SECTION_SHELL
$* secid elform shrf nip propt qr/irid icomp setyp
1 11 0.0 5 0 0.0
$* t1 t2 t3 t4 nloc marea
0.002000 0.002000 0.002000 0.002000 0 0.0
*MAT_PIECEWISE_LINEAR_PLASTICITY
$* mid ro e pr sigy etan fail tdel
1 0.0040002.0000e+05 0.330000 1.000000 0.0 0.0 0.0
$* c p lcss lcsr vp
0.0 0.0 2 0 0.0
$* eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$* es1 es2 es3 es4 es5 es6 es7 es8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
*PART
$* title
barre lagrangienne
$* pid secid mid eosid hgid grav adpopt tmid
2 2 2 2
*SECTION_SOLID
$* secid elform aet
2
*MAT_NULL
$* mid ro pc mu terod cerod ym pr
2 0.001000 -0.100000 0.001000 0.0 0.0 0.0 0.0
*EOS_GRUNEISEN
$* eosid c s1 s2 s3 gamao a e0
2 1647.000 2.480000 0.0 0.0 0.0 0.0 0.0
$* v0
0.0
*PART
$* title
plaq pour brique
$* pid secid mid eosid hgid grav adpopt tmid
3 1 1
*PART
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ SPH PART $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* title
barre sph
$* pid secid mid eosid hgid grav adpopt tmid
101 101 2 2
*SECTION_SPH
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *SECTION_SPH $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* secid cslh hmin hmax sphini death start
101 1.200000 0.0 0.0 0.0 0.0 0.0
*END