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. Download is available in the download section of this document.
Keywords
Keywords
This example has been taken from the SPH class of Jean Luc Lacome. Example name: impact.k
The used Keywords in alphabetical order.
*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
Reduced Input
Reduced Input
This example has been taken from the SPH class of Jean Luc Lacome. Example name: impact.k
*Node and *Element Cards are removed.
*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
Figures
Figures
This example has been taken from the SPH class of Jean Luc Lacome. Example name: impact.k
Animated Result
Animated Result
This example has been taken from the SPH class of Jean Luc Lacome. Example name: impact.k
Download
Download
This example has been taken from the SPH class of Jean Luc Lacome. Example name: impact.k
