*heading
 Test initial conditions TYPE=STRESS with 
 HYPERELASTIC, MARLOW + MULLINS EFFECT   
** all displacements are constrained during the analysis:
** 1) the stress field should remain constant
** 2) the damage energy should also be constant
*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.
*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.,
*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.,
*element,type=cps4r,elset=cps
 3,31,32,33,34
**
*node,nset=all
 41,0.,6.,
 42,1.,6.,
*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
*nset, nset=QA_TEST
17,22,32,42,52,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=solid, 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
*initial conditions,type=stress
  eall,100000.
*step
*dynamic, explicit
 ,1.0
*boundary,type=velocity
 all,1,6
*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
 rf,
*output, history, variable=preselect
*element output, elset=eall
 s11, le11, ne11, ener
*end step