*HEADING
 :SPRING-MASS SYSTEM WITH MASS PROPORTIONAL DAMPING AND ANISOTROPIC MASS
****        testing *mass, alpha=
****        The constants are such set that the results for U1 of NODE 3
****        should be the same as running 2.8.4-example A (exsd3glx.inp)
****        
*NODE,NSET=NALL
3,10.,0.
**
** Spring Definitions
**
*ELEMENT,TYPE=SPRING1,ELSET=SPR1
1,3
*ELEMENT,TYPE=SPRING1,ELSET=SPR2
2,3
*ELEMENT,TYPE=SPRING1,ELSET=SPR3
3,3
**
** Point Mass
**
*ELEMENT,TYPE=MASS,ELSET=MA
9,3
*NSET,NSET=MASSES
3,12,13,102,103
*SPRING,ELSET=SPR1
1, 
30.0, 
*SPRING,ELSET=SPR2
2, 
30.0,
*SPRING,ELSET=SPR3
3, 
30.0, 
*MASS,ELSET=MA,ALPHA=4.6367852,TYPE=ANISOTROPIC
0.02588,0.05176,0.07764
*RESTART,WRITE
**
**
**
**
**      ----------- STEP 1 ----------- 
**
*STEP,NLGEOM
  STRETCH OUT
*STATIC
1.E2,1.E2
*BOUNDARY
3,1,1,1.
3,2,3,0.
*NODE PRINT,FREQ=2,SUMMARY=NO,GLOBAL=NO
U
*END STEP
**
**
**      ------------ STEP 2 ------------
**
*STEP,INC=200,NLGEOM
 RELEASE
*DYNAMIC,HAFTOL=0.1
5.E-4,0.7,,0.7
*BOUNDARY,OP=NEW
**
** Removed previous BC's; let node 3 free in x dir
**
3,2,3,0.
*MONITOR,NODE=3,DOF=1
*NODE PRINT,FREQ=2,SUMMARY=NO,GLOBAL=NO
U,V,A
*EL PRINT, FREQ=0
*ENERGY PRINT
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*EL PRINT, SUMMARY=NO
 ELKE,ELVD,ELSE
*END STEP
**
**      ----------- STEP 3 ----------- 
**
*STEP,NLGEOM
  STRETCH OUT
*STATIC
1.E2,1.E2
*BOUNDARY,OP=NEW
3,2,2,1.
3,1,1,0.
3,3,3,0.
*NODE PRINT,FREQ=2,SUMMARY=NO,GLOBAL=NO
U
*END STEP
**
**
**      ------------ STEP 4 ------------
**
*STEP,INC=200,NLGEOM
 RELEASE
*DYNAMIC,HAFTOL=0.1
5.E-4,0.7,,0.7
*BOUNDARY,OP=NEW
**
** Removed previous BC's; let node 3 free in y dir
**
3,1,1,0.
3,3,3,0.
*MONITOR,NODE=3,DOF=1
*NODE PRINT,FREQ=2,SUMMARY=NO,GLOBAL=NO
U,V,A
*EL PRINT, FREQ=0
*ENERGY PRINT
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*EL PRINT, SUMMARY=NO
 ELKE,ELVD,ELSE
*END STEP
**
**      ----------- STEP 5 ----------- 
**
*STEP,NLGEOM
  STRETCH OUT
*STATIC
1.E2,1.E2
*BOUNDARY,OP=NEW
3,3,3,1.
3,1,2,0.
*NODE PRINT,FREQ=2,SUMMARY=NO,GLOBAL=NO
U
*END STEP
**
**
**      ------------ STEP 6 ------------
**
*STEP,INC=200,NLGEOM
 RELEASE
*DYNAMIC,HAFTOL=0.1
5.E-4,0.7,,0.7
*BOUNDARY,OP=NEW
**
** Removed previous BC's; let node 3 free in z dir
**
3,1,2,0.
*MONITOR,NODE=3,DOF=1
*NODE PRINT,FREQ=2,SUMMARY=NO,GLOBAL=NO
U,V,A
*EL PRINT, FREQ=0
*ENERGY PRINT
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*EL PRINT, SUMMARY=NO
 ELKE,ELVD,ELSE
*END STEP