*HEADING
 TRANSIENT LOADING OF A VISCOELASTIC BUSHING -- 
 CPE4R
*NODE
1,0.5,0.
7,1.0,0.
999,0.,0.
*NGEN,NSET=RADIAL
1,7
*NCOPY,CHANGE NUMBER=10,OLD SET=RADIAL,
NEW SET=ALL,SHIFT,MULTIPLE=55

0.,0.,0.,0.,0.,1.,6.4285714
*NSET,NSET=INSIDE,GENERATE
1,551,10
*NSET,NSET=OUTSIDE,GENERATE
7,557,10
*ELEMENT,TYPE=CPE4R,ELSET=ONE
1, 1,2,12,11
*ELGEN,ELSET=ALL
1,6,1,1,55,10,10
*ELSET,ELSET=TWO
1,2
*ELSET,ELSET=RIM,GENERATE
1,21,10
2,22,10
401,421,10
402,422,10
*ELEMENT,TYPE=CPE4R
551,551,552,2,1
*ELGEN,ELSET=ALL
551,6,1,1                                  
*KINEMATIC COUPLING, REF NODE=999
INSIDE,1,2
*BOUNDARY
OUTSIDE,1,2
*MATERIAL,NAME=RUBBER
*HYPERELASTIC,N=1,MODULI=INSTANTANEOUS
4000.,1000.,0.00002
*VISCOELASTIC,TIME=PRONY
0.2,0.5,.1
0.1,0.2,.2
*SOLID SECTION,ELSET=ALL,MATERIAL=RUBBER
1., 
*NSET,NSET=CENTER
999, 
*STEP,NLGEOM
*STATIC
0.001,0.001,0.0001
*CLOAD
999,1,50000.
*EL PRINT,FREQUENCY=0
*NODE PRINT,NSET=CENTER
*NODE FILE,FREQUENCY=150,NSET=RADIAL
U,RF
*EL FILE,FREQUENCY=150,ELSET=TWO
S,E,CE,CEP
*OUTPUT,VAR=PRESELECT,FIELD,FREQ=999
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT,NSET=CENTER
U, 
*OUTPUT,FIELD,Frequency=150
*NODE OUTPUT,NSET=RADIAL
U,RF
*OUTPUT,HISTORY,Frequency=150
*NODE OUTPUT,NSET=RADIAL
U,RF
*OUTPUT,FIELD,Frequency=150
*ELEMENT OUTPUT,ELSET=TWO
S,E,CE,CEP
*OUTPUT,HISTORY,Frequency=150
*ELEMENT OUTPUT,ELSET=TWO
S,E,CE,CEP
*END STEP
*STEP,NLGEOM,INC=150
*VISCO,CETOL=5.E-4
0.04,1.
*END STEP
*STEP,NLGEOM
*STATIC
0.001,0.001
*CLOAD
999,6,10000.
*NODE FILE,FREQUENCY=150,NSET=RADIAL
U,RF
*EL FILE,FREQUENCY=150,ELSET=RIM
S,E,CE,CEP
*OUTPUT,FIELD,Frequency=150
*NODE OUTPUT,NSET=RADIAL
U,RF
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT,NSET=CENTER
U, 
*OUTPUT,HISTORY,Frequency=150
*NODE OUTPUT,NSET=RADIAL
U,RF
*OUTPUT,VAR=PRESELECT,FIELD,FREQ=999
*OUTPUT,FIELD,Frequency=150
*ELEMENT OUTPUT,ELSET=RIM
S,E,CE,CEP
*OUTPUT,HISTORY,Frequency=150
*ELEMENT OUTPUT,ELSET=RIM
S,E,CE,CEP
*END STEP
*STEP,NLGEOM,INC=150
*VISCO,CETOL=5.E-4
0.04,1.
*END STEP