*heading
 ANISOTROPIC HYPERELASTIC, FUNG form
 Uniaxial test, displacement control
 Base problem for Explicit->Explicit Import
*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=c3d
**
*node,nset=all
*nset,nset=face21
*nset,nset=face22
*nset,nset=face23
*nset,nset=face24
*nset, nset=ndisp2
*element,type=cpe4r,elset=cpe
**
*node,nset=all
*nset,nset=face31
*nset,nset=face32
*nset,nset=face33
*nset,nset=face34
*nset, nset=ndisp3
*element,type=cps4r,elset=cps
**
*node,nset=all
*nset,nset=face41
*nset,nset=face42
*nset,nset=face43
*nset,nset=face44
*nset, nset=ndisp4
*element,type=m3d4r,elset=m3d
**
*node,nset=all
*nset,nset=face51
*nset,nset=face52
*nset,nset=face53
*nset,nset=face54
*nset, nset=ndisp5
*nset, nset=QA_TEST
*element,type=s4r,elset=s4r
**
*elset, elset=solid
 c3d, cpe, cps
*elset, elset=eall
 solid, m3d, s4r
*elset, elset = QA_TEST
eall,
**
*orientation,name=ori-1
**
*solid section,elset=solid,material=fung,orientation=ori-1
 1.0,
*membrane section,elset=m3d,material=fung,orientation=ori-1
 1.0,
*shell section,elset=s4r,material=fung,orientation=ori-1
 1.0,
**
*parameter
 c = 26.95e3
 d = 1.5e-7
 b1 = 0.9925
 b2 = 0.4180
 b3 = 0.0089
 b4 = 0.0749
 b5 = 0.0295
 b6 = 0.0193
 b7 = 5.0000
 b8 = b7
 b9 = b7
 b1111 = b1
 b1122 = b4
 b2222 = b2
 b1133 = b5
 b2233 = b6
 b3333 = b3
 b1212 = b7
 b1313 = b8
 b2323 = b9
 density=1.0e-2
*material,name=fung
*density
 <density>
*anisotropic hyperelastic,fung-ortho
 <b1111>,<b1122>,<b2222>,<b1133>,<b2233>,<b3333>,<b1212>,<b1313>,
 <b2323>,<c>,<d>
*equation
 face14,1,1,ndisp1,1,-1
 face22,1,1,ndisp2,1,-1
 face32,1,1,ndisp3,1,-1
 face42,1,1,ndisp4,1,-1
 face52,1,1,ndisp5,1,-1
*nset, nset=ndisp
 ndisp1, ndisp2, ndisp3, ndisp4, ndisp5
*boundary,type=velocity
 face11,3
 face13,2
 face16,1
**
 face21,2
 face24,1
**
 face31,2
 face34,1
**
 face41,2
 face44,1
 face41,3
 face43,3
**
 face51,2
 face54,1
 face51,3
 face53,3
**
*amplitude,name=bc
0.0,0.0,0.25,0.75,0.5,0.0,0.75,1.0,
**
*restart, write
**
*step
 uniaxial loading
*dynamic, explicit
 ,0.75
*boundary,amplitude=bc
 ndisp,1,1,1.0
*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, le33, ne11, ne22, ener
*node out, nset=ndisp
 rf1, u1
*end step