*HEADING
**
** Creep analysis of a rod. -- B31 elements
** Verification of *HYPERELASTIC, Neo Hooke and *VISCOELASTIC
** incompressible
** 
*PREPRINT, MODEL=YES
*NODE,NSET=nsolid
1,0.,0.
2,1.,0.
3,1.,1
4,0.,1
5,0,0,1
6,1,0,1
7,1,1,1
8,0,1,1
*NODE,NSET=nbeam
11,3.
12,4.
*nset,nset=alln
nbeam,nsolid
*ELEMENT,TYPE=B31,ELSET=ebeam
11,11,12
*element,type=c3d8RH,elset=esolid
1,1,2,3,4,5,6,7,8
*BEAM SECTION, MATERIAL=TEST,SECTION=CIRC,ELSET=ebeam, poisson=0.5
 0.564189588
*solid section,elset=esolid,material=test
*elset,elset=alle
ebeam,esolid
*MATERIAL, NAME=TEST
*HYPERELASTIC, NEO HOOKE, MODULI=INSTANTANEOUS
1, 0
*VISCOELASTIC,TIME=PRONY, FORMULATION=SOLID
0.5,,1
**
**
*BOUNDARY
 1,1,3
 2,2,3
 3,3,3
 4,1,1
 4,3,3
 5,1,2
 6,2,2
 8,1,1
11,1,6
12,2,6
**
*AMPLITUDE, NAME=PULL
0,0,1,1,2,0,3,1
4,0,5,1,6,0,7,1
8,0,9,1,10,0
**
*STEP, NLGEOM, inc=100000
 SUSTAIN LOAD
*VISCO, cetol=0.01
 0.001,10
*CLOAD, AMPLITUDE=PULL
5,3,0.25
6,3,0.25
7,3,0.25
8,3,0.25
12,1,1
*OUTPUT,FIELD, VARIABLE=PRESELECT
*OUTPUT, HISTORY, VARIABLE=PRESELECT
*ELEMENT OUTPUT, ELSET=ebeam
 S11, LE11, NE11,elcd
*ELEMENT OUTPUT, ELSET=esolid
 S33, LE33, NE33,elcd
*NODE OUTPUT, NSET=ALLN
U
*NODE OUTPUT, NSET=nbeam
 RF1,CF1
*NODE OUTPUT, NSET=nsolid
 RF3,CF3
*END STEP