*heading
 HYPERELASTIC, MARLOW with MULLINS EFFECT   
*preprint, model=yes
*node,nset=all
 11,
 12,1.
 13,1.,1.,
 14,0.,1.,
 15,0.,0.,1.
 16,1.,0.,1.
 17,1.,1.,1.
 18,0.,1.,1.
*nset,nset=face11
 11,12,13,14
*nset,nset=face12
 15,16,17,18
*nset,nset=face13
 11,12,15,16
*nset,nset=face14
 13,16,17
*nset,nset=face15
 13,14,17,18
*nset,nset=face16
 14,11,18,15
*nset, nset=ndisp1
 12
*element,type=c3d8r,elset=c3d
 1,11,12,13,14,15,16,17,18
**
*node,nset=all
 21,0.,2.,
 22,1.,2.,
 23,1.,3.,
 24,0.,3.,
*nset,nset=face21
 21,22
*nset,nset=face22
 23
*nset,nset=face23
 23,24
*nset,nset=face24
 21,24
*nset, nset=ndisp2
 22
*element,type=cpe4r,elset=cpe
 2,21,22,23,24
**
*node,nset=all
 31,0.,4.,
 32,1.,4.,
 33,1.,5.,
 34,0.,5.,
*nset,nset=face31
 31,32
*nset,nset=face32
 33
*nset,nset=face33
 33,34
*nset,nset=face34
 31,34
*nset, nset=ndisp3
 32
*element,type=cps4r,elset=cps
 3,31,32,33,34
**
*node,nset=ntruss
 41,0.,6.,
 42,1.,6.,
*nset, nset=ndisp4
 42
*nset, nset=nfix
 41
*element,type=t3d2,elset=truss
 4,41,42
**
*node,nset=all
 51,0.,7.,
 52,1.,7.,
 53,1.,8.,
 54,0.,8.,
*nset,nset=face51
 51,52
*nset,nset=face52
 53
*nset,nset=face53
 53,54
*nset,nset=face54
 51,54
*nset, nset=ndisp5
 52
*element,type=m3d4r,elset=m3d
 5,51,52,53,54
**
*node,nset=all
 61,0.,9.,
 62,1.,9.,
 63,1.,10.,
 64,0.,10.,
*nset,nset=face61
 61,62
*nset,nset=face62
 63
*nset,nset=face63
 63,64
*nset,nset=face64
 61,64
*nset, nset=ndisp6
 62
*element,type=s4r,elset=s4r
 6,61,62,63,64
**
*elset, elset=solid
 c3d, cpe, cps, truss
*elset, elset=eall
 solid, m3d, s4r
*elset, elset=QA_TEST
 eall,
**
*solid section,elset=c3d,material=hyper
 1.0,
*solid section,elset=cpe,material=hyper
 1.0,
*solid section,elset=cps,material=hyper
 1.0,
*solid section,elset=truss,material=hyper
 1.0,
*membrane section, elset=m3d, material=hyper
 1.0,
*shell section, elset=s4r, material=hyper
 1.0,
**
*material,name=hyper
*density
 100., 
*hyperelastic,marlow,poisson=0.49
*uniaxial test data, smooth
  0.0000  ,  0.0000
  1.5506e5,  0.1338
  2.4367e5,  0.2675
  3.1013e5,  0.3567
  4.2089e5,  0.6242
  5.3165e5,  0.8917
  5.9810e5,  1.1592
  6.8671e5,  1.4268
  8.8608e5,  2.0510
 10.6329e5,  2.5860
 12.4051e5,  3.0318
 16.1709e5,  3.7898
 19.9367e5,  4.3694
 23.4810e5,  4.8153
 27.4684e5,  5.1720
 31.0127e5,  5.4395
 34.5570e5,  5.7070
 38.3228e5,  5.9299
 42.0886e5,  6.0637
 45.6329e5,  6.1975
 49.3987e5,  6.3312
 53.1646e5,  6.4650
 56.9304e5,  6.5541
 64.2405e5,  6.6433
*mullins effect
 1.6854, 10000.0, 0.1
*equation
 2,
 face14,1,1,ndisp1,1,-1
 2,
 face22,1,1,ndisp2,1,-1
 2,
 face32,1,1,ndisp3,1,-1
 2,
 face52,1,1,ndisp5,1,-1
 2,
 face62,1,1,ndisp6,1,-1
*nset, nset=ndisp
 ndisp1, ndisp2, ndisp3, ndisp4, ndisp5, ndisp6
*nset, nset=QA_TEST
 ndisp,
*boundary,type=velocity
 face11,3
 face13,2
 face16,1
**
 face21,2
 face24,1
**
 face31,2
 face34,1
**
 ntruss,2,3
 nfix,1
**
 face51,2
 face54,1
**
 face61,2
 face64,1
**
*step
 uniaxial loading
*dynamic, explicit
 ,1.0
*boundary,type=velocity
 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, ne11, ener
*node out, nset=ndisp
 rf1
*energy output
 alldmd
*end step
*step
 unloading
*dynamic, explicit
 ,1.0
*boundary,type=velocity
 ndisp,1,1,-1.
*end step