*heading
 HYPERELASTIC, MARLOW  with VISCOELASTICITY
 Creep test
*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=nload1
 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=nload2
 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=nload3
 32
*element,type=cps4r,elset=cps
 3,31,32,33,34
**
*node,nset=ntruss
 41,0.,6.,
 42,1.,6.,
*nset, nset=nload4
 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=nload5
 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=nload6
 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
 solid,m3d,62,63,64
**
*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
 10., 
*hyperelastic,marlow,moduli=instantaneous
*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
*volumetric test data
   0.000, 1.0000
  60.0e5, 0.9703
 118.2e5, 0.9412
 175.2e5, 0.9127
 231.1e5, 0.8847
*viscoelastic,time=relaxation test data
*combined test data,shrinf=0.50,volinf=0.50
 0.99256 , 0.99256 , 0.10000    
 0.98525 , 0.98525 , 0.20000    
 0.97805 , 0.97805 , 0.30000    
 0.96402 , 0.96402 , 0.50000    
 0.93730 , 0.93730 , 0.90000    
 0.88886 , 0.88886 , 1.7000     
 0.80894 , 0.80894 , 3.3000     
 0.69864 , 0.69864 , 6.5000     
 0.58776 , 0.58776 , 12.900     
 0.52060 , 0.52060 , 25.700     
 0.50149 , 0.50149 , 51.300     
 0.50011 , 0.50011 , 76.900     
 0.50001 , 0.50001 , 100.00 
*equation
 2,
 face14,1,1,nload1,1,-1
 2,
 face22,1,1,nload2,1,-1
 2,
 face32,1,1,nload3,1,-1
 2,
 face52,1,1,nload5,1,-1
 2,
 face62,1,1,nload6,1,-1
*nset, nset=nload
 nload1, nload2, nload3, nload4, nload5, nload6
*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
**
*amplitude, definition=smooth step, name=creep
 0.0, 0.0, 0.1, 1.0, 10.0, 1.0
**
*step
 uniaxial loading
*dynamic, explicit
 ,10.0
*cload, amplitude=creep
 nload,1,6.0e5
*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
*end step