*HEADING
**
** Creep analysis of a rod. -- B31 elements
** Verification of *HYPERELASTIC, MARLOW and *VISCOELASTIC
** compressible
** 
*PREPRINT, MODEL=YES
*NODE,NSET=ALL
 1,
 2,1.
*NSET, NSET=NDISP
 2
*ELEMENT,TYPE=B21,ELSET=EL
1,1,2
*BEAM SECTION, MATERIAL=MARLOW,SECTION=CIRC,ELSET=EL,POISSON=0.475
 0.564189588
*MATERIAL,NAME=MARLOW
*HYPERELASTIC,MARLOW, MODULI = INSTANTANEOUS,POISSON=0.475
*UNIAXIAL TEST DATA
  0.0, 0.0
  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
*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 
*STEP, NLGEOM
 UNIAXIAL LOADING
*STATIC
 0.01,1.0,,0.05
*BOUNDARY
 1,1,2
 2,2
*CLOAD
 2,1,65.0E5
*OUTPUT,FIELD, VARIABLE=PRESELECT
*OUTPUT, HISTORY, VARIABLE=PRESELECT
*ELEMENT OUTPUT, ELSET=EL
 S, LE11, NE11, EE11, CE11
*NODE OUTPUT, NSET=NDISP
 RF1,U,CF1
*END STEP
*STEP, NLGEOM
 SUSTAIN LOAD
*VISCO
 5.,100.,2.,10.,
*END STEP