*HEADING     
 HYPERELASTIC WITH VAN DER WAALS STRAIN ENERGY POTENTIAL
 TEST DATA INPUT, UNIAXIAL TENSION (SI UNIT)
 USING UNIAXIAL TENSION DATA ONLY TO DO THE CURVE FIT 
**
*RESTART,TIMEMARKS=NO,WRITE,NUM=18
*NODE,NSET=ALL
 1,
 2, 1.
 3, 1., 1.,
 4, 0., 1.,
 5, 0., 0., 1.
 6, 1., 0., 1.
 7, 1., 1., 1.
 8, 0., 1., 1.
*NSET,NSET=FACE1
 1,2,3,4
*NSET,NSET=FACE2
 5,6,7,8
*NSET,NSET=FACE3
 1,2,5,6
*NSET,NSET=FACE4
 2,3,6,7
*NSET,NSET=FACE5
 3,4,7,8
*NSET,NSET=FACE6
 4,1,8,5
*ELEMENT,TYPE=C3D8R,ELSET=ONE
 1,1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=ONE,MATERIAL=TRELOAR
*MATERIAL,NAME=TRELOAR
*DENSITY
 1000.,
*HYPERELASTIC,VAN DER WAALS,TEST DATA INPUT,BETA=0.
*UNIAXIAL TEST DATA
  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
*AMPLITUDE,NAME=RAMP
 0., 0., 2., 6.
**
*STEP
 UNIAXIAL TENSION
*DYNAMIC,EXPLICIT
 ,2.0
*BOUNDARY,OP=NEW,AMP=RAMP,TYPE=VELOCITY
 FACE1,3,3,0.0
 FACE3,2,2,0.0
 FACE6,1,1,0.0
 FACE4,1,1,1.0
**
*NSET,NSET=ONE
 2,
*OUTPUT,HISTORY,TIME INTERVAL=0.0
*ELEMENT OUTPUT,ELSET=ONE,VAR=PRESELECT
 NE11,
*NODE OUTPUT,NSET=ONE,VAR=PRESELECT
 RF1,
*OUTPUT,FIELD,VAR=PRESELECT,NUMBER INTERVAL=5
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST
ALL,
*OUTPUT,FIELD,TIMEMARKS=YES,NUMBERINTERVAL=4
*NODE OUTPUT,NSET=QA_TEST
U,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.5
*FILE OUTPUT,TIMEMARKS=YES,NUMBER INTERVAL=4
*NODE FILE
 U,
**
*END STEP