*heading
 ISOTROPIC HYPERELASTIC, VISCOELASTICITY and MULLINS EFFECT
 Uniaxial test, displacement control
*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=face41
*nset,nset=face42
*nset,nset=face43
*nset,nset=face44
*nset,nset=face45
*nset,nset=face46
*nset, nset=ndisp1
*nset, nset=ndisp4
*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
*nset, nset=QA_TEST
*elset, elset=solid
c3d,cpe
*elset, elset = QA_TEST
c3d,cpe,cps
*solid section,elset=solid,material=holzapfel
 1.0,
*parameter
 density=5.0e-4
 C10=7.64
 D=1.e-6
*material,name=holzapfel
*density
 <density>
*hyperelastic,neoh
 <C10>, <D>
*mullins effect
*viscoelastic, time=prony
 0.2,,0.05
*solid section,elset=cps,material=holzinc
 1.0,
*material,name=holzinc
*density
 <density>
*hyperelastic,neoh
 <C10>,
*mullins effect
*viscoelastic, time=prony
 0.2,,0.05
*equation
 face14,1,1,ndisp1,1,-1
 face22,1,1,ndisp2,1,-1
 face32,1,1,ndisp3,1,-1
 face44,1,1,ndisp4,1,-1
*nset, nset=ndisp
 ndisp1,ndisp2,ndisp3,ndisp4
*boundary
 face11,3
 face13,2
 face16,1
**
 face41,3
 face42,3
 face43,2
 face46,1
**
 face21,2
 face24,1
**
 face31,2
 face34,1
*amplitude,name=bc
0.0,0.0,0.25,0.75,0.5,0.0,0.75,1.0,
1.0,0.0
**
*step,nlgeom,inc=1000
 uniaxial loading
*visco,cetol=0.01
0.01,1.0,,0.1
*boundary,amplitude=bc
 ndisp,1,1,0.3
*output,field, variable=preselect
*element output, elset=qa_test
 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=qa_test
 s11, le11, le22, le33, ne11, ne22, ener
*node out, nset=ndisp
 rf1, u1
*end step