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. 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: foam.k
The used Keywords in alphabetical order.
*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
Reduced Input
Reduced Input
This example has been taken from the SPH class of Jean Luc Lacome. Example name: foam.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 = 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
Figures
Figures
This example has been taken from the SPH class of Jean Luc Lacome. Example name: foam.k
Animated Result
Animated Result
This example has been taken from the SPH class of Jean Luc Lacome. Example name: foam.k
Download
Download
This example has been taken from the SPH class of Jean Luc Lacome. Example name: foam.k
