Foam
This models compression of a foam between a planar rigid wall and a solid block. The foam material is specified with *MAT_LOW_DENSITY_FOAM keyword and discretized using SPH method.
https://www.dynaexamples.com/sph/foam
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Foam
This models compression of a foam between a planar rigid wall and a solid block. The foam material is specified with *MAT_LOW_DENSITY_FOAM keyword and discretized using SPH method.
*CONSTRAINED_RIGID_BODIES *CONTACT_AUTOMATIC_NODES_TO_SURFACE_TITLE *CONTROL_BULK_VISCOSITY *CONTROL_CONTACT *CONTROL_CPU *CONTROL_DAMPING *CONTROL_ENERGY *CONTROL_HOURGLASS *CONTROL_OUTPUT *CONTROL_SHELL *CONTROL_SPH *CONTROL_STRUCTURED *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_BINARY_D3THDT *DATABASE_BINARY_RUNRSF *DATABASE_DEFORC *DATABASE_ELOUT *DATABASE_EXTENT_BINARY *DATABASE_GLSTAT *DATABASE_HISTORY_NODE *DATABASE_HISTORY_SOLID *DATABASE_MATSUM *DATABASE_NODOUT *DATABASE_RBDOUT *DATABASE_RCFORC *DEFINE_CURVE *ELEMENT_MASS *ELEMENT_SHELL *ELEMENT_SOLID *ELEMENT_SPH *END *INITIAL_VELOCITY_NODE *KEYWORD *MAT_LOW_DENSITY_FOAM *MAT_RIGID *NODE *PART *RIGIDWALL_PLANAR *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 = 1725,Largest node ID = 1001000
$* Total Number of Shells = 180,Largest Shell ID = 1743
$* Total Number of Solids = 432,Largest Solid ID = 1359
$* Total Number of Parts = 4,Largest Part ID = 101
$* Total Number of Mats = 4 Largest Mat ID = 6
*TITLE
$ curve # 1 for foam property
$ 1. when y-scale is 100. - contact OK
$ 2. when y-scale is 1.0 - negative volume (element compressed to pancake)
test of material 57 - foam
*CONTROL_BULK_VISCOSITY
$* q1 q2 type
1.500000 0.060000
*CONTROL_CONTACT
$* slsfac rwpnal islchk shlthk penopt thkchg orien enmass
0.100000 1.000000 2 0 1 0 1
$* 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 2
*CONTROL_HOURGLASS
$* ihq qh
5 0.100000
*CONTROL_OUTPUT
$* npopt neecho nrefup iaccop opifs ipnint ikedit iflush
0 0 1 1 0.0 1
*CONTROL_SHELL
$* wrpang esort irnxx istupd theory bwc miter proj
0.0 1 0 0 2
*CONTROL_SPH
$* ncbs boxid dt dim memory form start maxv
1 0 0.0 0 0 0 0.0 0.0
*CONTROL_STRUCTURED
*CONTROL_TERMINATION
$* endtim endcyc dtmin endeng endmas
10.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_DEFORC
0.200000
*DATABASE_ELOUT
0.100000
*DATABASE_GLSTAT
0.200000
*DATABASE_MATSUM
0.200000
*DATABASE_NODOUT
0.100000
*DATABASE_RBDOUT
0.100000
*DATABASE_RCFORC
0.100000
*DATABASE_BINARY_D3PLOT
1.000000
*DATABASE_BINARY_D3THDT
0.100000
*DATABASE_BINARY_RUNRSF
10000
*DATABASE_EXTENT_BINARY
0 0 0 1
0 0 0 0 0 0 0
*DATABASE_HISTORY_SOLID
928
*DATABASE_HISTORY_NODE
1
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_TITLE
$ INTERFACE NAME: 2 $$$
$ 101 5 4 3 0 0 0 0
$ .000E+00 .000E+00 .000E+00 .000E+00 .000E+00 0 .000E+00 .000E+00
$ 10.000000 0.0100000 .000E+00 .000E+00 .000E+00 .000E+00
$ 1 0.1000000 3
$ 1
$ 4 5
$* cidtitle
1
$* ssid msid sstyp mstyp sboxid mboxid spr mpr
101 4 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
10.00000 0.010000 0.0 0.0 0.0 0.0 0.0 0.0
*SET_NODE_LIST_GENERATE
101 0.0 0.0 0.0 0.0
1000001 1001000
*DEFINE_CURVE
$
*ELEMENT_SOLID
$
*ELEMENT_SHELL
$
*ELEMENT_SPH
$
*ELEMENT_MASS
1 25 20.00000000
2 621 20.00000000
*NODE
$
*INITIAL_VELOCITY_NODE
$
*CONSTRAINED_RIGID_BODIES
$ INTERFACE NAME: 2 $$$
$ 101 5 4 3 0 0 0 0
$ .000E+00 .000E+00 .000E+00 .000E+00 .000E+00 0 .000E+00 .000E+00
$ 10.000000 0.0100000 .000E+00 .000E+00 .000E+00 .000E+00
$ 1 0.1000000 3
$ 1
$ 4 5
4 6
*PART
$* title
PSOLID : 1 CHEXA :large.fm
$* pid secid mid eosid hgid grav adpopt tmid
4 4 4
*SECTION_SOLID
$* secid elform aet
4
*MAT_RIGID
$ 1 2.400E-07 2.500E-01 0.0000000 1.000E+00 4.000E-02 .000E+00
$* mid ro e pr n couple malias
41.0000e-07 207.0000 0.280000 0.0 0.0 0.0
$* cmo con1 con2
0.0 0.0 0.0
$lco or a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
*PART
$* title
PSHELL : 61 CQUAD4:small.rg
$* pid secid mid eosid hgid grav adpopt tmid
5 5 5
*SECTION_SHELL
$* secid elform shrf nip propt qr/irid icomp setyp
5 2 0.0 0 0 0.0
$* t1 t2 t3 t4 nloc marea
1.000000 1.000000 1.000000 1.000000 0 0.0
*MAT_RIGID
$* mid ro e pr n couple malias
51.0000e-07 207.0000 0.280000 0.0 0.0 0.0
$* cmo con1 con2
0.0 0.0 0.0
$lco or a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
*PART
$* title
PSHELL : 61 CQUAD4:large.rg
$* pid secid mid eosid hgid grav adpopt tmid
6 6 6
*SECTION_SHELL
$1.000E-04 1.000E-04 1.000E-04 1.000E-04 .000E+00
$* secid elform shrf nip propt qr/irid icomp setyp
6 2 0.0 0 0 0.0
$* t1 t2 t3 t4 nloc marea
1.0000e-041.0000e-041.0000e-041.0000e-04 0 0.0
*MAT_RIGID
$* mid ro e pr n couple malias
61.0000e-07 207.0000 0.280000 0.0 0.0 0.0
$* cmo con1 con2
0.0 0.0 0.0
$lco or a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
*PART
$* title
$* pid secid mid eosid hgid grav adpopt tmid
101 101 1
*SECTION_SPH
$* secid cslh hmin hmax sphini death start
101 1.200000 0.0 0.0 0.0 0.0 0.0
*MAT_LOW_DENSITY_FOAM
$* mid ro e lcid tc hu beta damp
12.4000e-07 0.250000 1 0.040000 0.500000 0.0 0.0
$* shape fail bvflag ed beta1 kcon ref
0.0 0.0 0.0 0.0 0.0 0.0 0.0
*RIGIDWALL_PLANAR
101 0 0 0.0
0.0 0.0 0.0 0.0 0.0 1.000000 0.0 0.0
*END