*heading
 LOW DENSITY FOAM
 Uniaxial test with stepped velocity conditions
*preprint, model=yes
*node,nset=all
*nset,nset=face11
*nset,nset=face12
*nset,nset=face13
*nset,nset=face14
*nset,nset=face15
*nset,nset=face16
*nset, nset=ndisp1
*element,type=c3d8r,elset=c3d1
**
*node,nset=all
*nset,nset=face111
*nset,nset=face112
*nset,nset=face113
*nset,nset=face114
*nset,nset=face115
*nset,nset=face116
*nset, nset=ndisp101
*element,type=c3d8r,elset=c3d2
**
*node,nset=all
*nset,nset=face21
*nset,nset=face22
*nset,nset=face23
*nset,nset=face24
*nset, nset=ndisp2
*element,type=cpe4r,elset=cpe1
**
*node,nset=all
*nset,nset=face121
*nset,nset=face122
*nset,nset=face123
*nset,nset=face124
*nset, nset=ndisp102
*element,type=cpe4r,elset=cpe2
**
*node,nset=ntruss1
*nset, nset=ndisp4
*nset, nset=nfix1
*element,type=t3d2,elset=truss1
**
*node,nset=ntruss2
*nset, nset=ndisp104
*nset, nset=nfix2
*element,type=t3d2,elset=truss2
**
*elset, elset=solid1
 c3d1, cpe1, truss1
*elset, elset=solid2
 c3d2, cpe2, truss2
*elset, elset=eall
 solid1, solid2
*elset, elset = QA_TEST
eall,
*nset, nset=QA_TEST
**
*solid section,elset=c3d1,material=lowdensfoam_vol
 1.0,
*solid section,elset=cpe1,material=lowdensfoam_vol
 1.0,
*solid section,elset=truss1,material=lowdensfoam_vol
 1.0,
*solid section,elset=c3d2,material=lowdensfoam_prin
 1.0,
*solid section,elset=cpe2,material=lowdensfoam_prin
 1.0,
*solid section,elset=truss2,material=lowdensfoam_prin
 1.0,
**
*material,name=lowdensfoam_vol
*DENSITY
1.824E-10
*LOW DENSITY FOAM
**, STRAIN RATE=VOLUMETRIC (default)
** Use default values of mu0, mu1, alpha
**0.0001, 0.005, 2.0
*UNIAXIAL TEST DATA, DIRECTION=TENSION
*UNIAXIAL TEST DATA, DIRECTION=COMPRESSION
*material,name=lowdensfoam_prin
*DENSITY
1.824E-10
*LOW DENSITY FOAM, STRAIN RATE=PRINCIPAL
** Use default of mu0, mu1, alpha
**0.0001, 0.005, 2.0
*UNIAXIAL TEST DATA, DIRECTION=TENSION
*UNIAXIAL TEST DATA, DIRECTION=COMPRESSION
*equation
 face14,1,1,ndisp1,1,-1
 face114,1,1,ndisp101,1,-1
 face22,1,1,ndisp2,1,-1
 face122,1,1,ndisp102,1,-1
*nset, nset=ndisp
 ndisp1, ndisp2, ndisp4
 ndisp101, ndisp102, ndisp104
*boundary,type=velocity
 face11,3
 face13,2
 face16,1
**
 face111,3
 face113,2
 face116,1
**
 face21,2
 face24,1
**
 face121,2
 face124,1
**
 ntruss1,2,3
 nfix1,1
**
 ntruss2,2,3
 nfix2,1
**
*amplitude, name=vel_step
0.0, 10, 0.03, 10, 0.031, 50, 0.036, 50,
0.037, 100, 0.038, 100, 0.039, -10, 0.1, -10
*step
 uniaxial loading
*dynamic, explicit
 ,0.1
*boundary,type=velocity,amplitude=vel_step
 ndisp,1,1,-1
*output,field, variable=preselect
*element output, elset=eall
 s, le
*node out, nset=all
 u, rf
*output,field, number=2
*element output, elset=QA_TEST
 s,
*node output, nset=QA_TEST
 u,
*output, history, variable=preselect
*element output, elset=eall
 s11, le11, le22, ne11, ne22, ener
*node out, nset=ndisp
 rf1
*end step