*HEADING
 - HYPERFOAM-VISCOELASTICITY, VISCO, LINEAR PERT. and Multiple load case STEPS,  CPE4
 ELSET TWO ADDED AS A TEST TO GUARD AGAINST v55_104
*NODE, NSET=NALL
1,
2,10.
3,10.,5.,
4,0.,5.,
*NODE, NSET=NAL
5,20.0, 0.0
6,30.0, 0.0
7, 30.0, 5.0
8, 20.0, 5.0
*ELEMENT,TYPE=CPE4,ELSET=ONE
1,1,2,3,4
*ELEMENT, TYPE=CPE4, ELSET=TWO
2, 5, 6, 7, 8
*SOLID SECTION,ELSET=ONE,MATERIAL=SIMPLE
1., 
*SOLID SECTION, ELSET=TWO, MATERIAL=S2
*MATERIAL,NAME=SIMPLE
*HYPERFOAM,N=3,MODULI=INSTANTANEOUS
-17.4,-1.22,548.2,17.3,10.47,-1.775

*VISCOELASTIC,TIME=PRONY
0.5,,3.
*DENSITY
4.0E-13, 
*MATERIAL,NAME=S2
*HYPERFOAM,N=3,MODULI=LONG TERM
-17.4,-1.22,548.2,17.3,10.47,-1.775

*VISCOELASTIC,FREQUENCY=TABULAR
0.,0.,0.,0.,0.
0.,0.,0.,0.,25.
*DENSITY
4.0E-13, 
*AMPLITUDE,NAME=STEP
0.,1.,10.,1.
*BOUNDARY
1,1,2
2,2,2
4,1,1
5,1,2
6,2,2
8,1,1 
*STEP,PERTURBATION
** Step 1: Initial Multiple Load Cases-1
*STATIC
*LOAD CASE,name=initial-1
*BOUNDARY
1,1,2
2,2,2
4,1,1
3,1,1,.1
2,1,1,.1
5,1, 2
6,2, 2
8,1, 1
7, 1, 1, 0.1
6, 1, 1, 0.1
*END LOAD CASE
*LOAD CASE,name=initial-2
*BOUNDARY
1,1,2
2,2,2
4,1,1
3,1,1,.1
2,1,1,.1
5,1, 2
6,2, 2
8,1, 1
7, 1, 1, 0.2
6, 1, 1, 0.2
*END LOAD CASE
*EL PRINT
S,E
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,E
*NODE OUTPUT
U,RF
*END STEP
*STEP,NLGEOM
** Step 2: Apply instantaneous uniaxial tensile stretch and allow to relax
*VISCO
2.,6.
*BOUNDARY,AMP=STEP
3,1,1,2.
2,1,1,2.
7,1,1,2.
6,1,1,2.
*EL PRINT,F=10
S,E
EE, 
CE, 
*EL FILE,F=10,ELSET=ONE
S,E,EE,CE
*NODE FILE,F=10,NSET=NALL
U,RF
*END STEP
*STEP,PERT
** Step 3: Intermediate Multiple Load Case
*STATIC
*LOAD CASE,name=intermediate-1
*BOUNDARY
1,1,2
2,2,2
4,1,1
3,1,1,.1
2,1,1,.1
5, 1, 2
6, 2, 2
8, 1, 1
7,1, 1,0.1
6,1, 1,0.1
*END  LOAD CASE
*LOAD CASE,name=intermediate-2
*BOUNDARY
1,1,2
2,2,2
4,1,1
3,1,1,.1
2,1,1,.1
5, 1, 2
6, 2, 2
8, 1, 1
7,1, 1,0.2
6,1, 1,0.2
*END  LOAD CASE
*EL PRINT
S,E
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,E
*NODE OUTPUT
U,RF
*END STEP
*STEP,NLGEOM
** Step 4: Continue stress relaxation
*VISCO 
2.,4.
*END STEP
*STEP,PERT
*STATIC
** Step 5: Final Multiple load case-1
*LOAD CASE,name=Final-1
*BOUNDARY
1,1,2
2,2,2
4,1,1
3,1,1,.1
2,1,1,.1
5, 1, 2
6, 2, 2
8, 1, 1
6, 1,1,0.1
7,1,1,0.1
*END LOAD CASE
*LOAD CASE,name=Final-2
*BOUNDARY
1,1,2
2,2,2
4,1,1
3,1,1,.1
2,1,1,.1
5, 1, 2
6, 2, 2
8, 1, 1
6, 1,1,0.2
7,1,1,0.2
*END LOAD CASE
*EL PRINT
S,E

*END STEP
*STEP,PERT
*STATIC
** Step 6: Final linear perturbation-2
*BOUNDARY
1,1,2
2,2,2
4,1,1
3,1,1,.1
2,1,1,.1
5, 1, 2
6, 2, 2
8, 1, 1
6, 1,1,0.2
7,1,1,0.2
*EL PRINT
S,E
*OUTPUT,FIELD
*ELEMENT OUTPUT,VARIABLE=PRESELECT
*NODE OUTPUT,VARIABLE=PRESELECT
*END STEP
*STEP,NLGEOM,INC=1000
*STEADY STATE DYNAMICS,DIRECT,FREQUENCY SCALE=LINEAR
 0.,5.,5,1
*BOUNDARY
1,1,2
2,2,2
4,1,1
3,1,1,.1
2,1,1,.1
5, 1, 2
6, 2, 2
8, 1, 1
6, 1, 1, 0.1
7, 1, 1, 0.1
*NODE FILE, FREQ=10,NSET=NAL
 U, 
*EL FILE,ELSET=TWO,FREQ=10
S, 
*END STEP