*HEADING
 NAFEMS FORCED VIBRATION; TEST 21T; [SC8R]
 (DIRECT INTEGRATION)
*NODE
   1,  0.,  0., -0.5
   9, 10.,  0., -0.5
  81,  0., 10., -0.5
  89, 10., 10., -0.5
 101,  0.,  0., +0.5
 109, 10.,  0., +0.5
 181,  0., 10., +0.5
 189, 10., 10., +0.5
*NGEN,NSET=BOT
   1,  9
 101,109
*NGEN,NSET=TOP
  81, 89
 181,189
*NFILL,NSET=NALL
 BOT,TOP,8,10
*NSET,NSET=LEFT,GEN
 1,81,10
 101,181,10
*NSET,NSET=RIGHT,GEN
 9,89,10
 109,189,10
*NSET,NSET=EDGE
 TOP,BOT,RIGHT,LEFT
*NSET,NSET=NMID
  45, 
** 145
*ELEMENT,TYPE=SC8R
 1, 1,2,12,11, 101,102,112,111
*ELGEN,ELSET=ALLE
 1, 8,1,1, 8,10,10
*ELSET,ELSET=EMID
 34,35,44,45
*SHELL SECTION,MATERIAL=MAT,ELSET=ALLE, POISSON=0.0, THICK=200.E9
 1.0, 
*MATERIAL,NAME=MAT
*DENSITY
 8000.,
*ELASTIC
 200.E9,0.3
*DAMPING,ALPHA=5.772,BETA=6.929E-5
*NSET, NSET=BOT_FACE, GEN
 1, 9, 1
 11,19,1
 21,29,1
 31,39,1
 41,49,1
 51,59,1
 61,69,1
 71,79,1
 81,89,1
*NSET, NSET=TOP_FACE, GEN
 101,109,1
 111,119,1
 121,129,1
 131,139,1
 141,149,1
 151,159,1
 161,169,1
 171,179,1
 181,189,1
*NSET, NSET=BOT_EDGE
 1,2,3,4,5,6,7,8,9
 11,19,21,29,31,39,41,49,51,59,
 61,69,71,79,
 81,82,83,84,85,86,87,88,89
*NSET, NSET=TOP_EDGE
 101,102,103,104,105,106,107,108,109
 111,119,121,129,131,139,141,149,151,159,
 161,169,171,179,
 181,182,183,184,185,186,187,188,189
*NSET, NSET=BOT_INT_FACE, GEN
 12,18,1
 22,28,1
 32,38,1
 42,48,1
 52,58,1
 62,68,1
 72,78,1
*NSET, NSET=TOP_INT_FACE, GEN
 112,118,1
 122,128,1
 132,138,1
 142,148,1
 152,158,1
 162,168,1
 172,178,1
*NSET, NSET=BOT_EDGE_A, GEN
 1, 9, 1
*NSET, NSET=TOP_EDGE_A, GEN
 101, 109, 1
*NSET, NSET=BOT_EDGE_C, GEN
 81, 89, 1
*NSET, NSET=TOP_EDGE_C, GEN
 181, 189, 1
*NSET, NSET=BOT_EDGE_B, GEN
 9, 89, 10
*NSET, NSET=TOP_EDGE_B, GEN
 109, 189, 10
*NSET, NSET=BOT_EDGE_D, GEN
 1, 81, 10
*NSET, NSET=TOP_EDGE_D, GEN
 101, 181, 10
*NSET, NSET=EDGE_A
 BOT_EDGE_A, TOP_EDGE_A
*NSET, NSET=EDGE_B
 BOT_EDGE_B, TOP_EDGE_B
*NSET, NSET=EDGE_C
 BOT_EDGE_C, TOP_EDGE_C
*NSET, NSET=EDGE_D
 BOT_EDGE_D, TOP_EDGE_D
**
** EQUIVALENT BOUNDARY CONDITIONS
**
** Fix interior nodes in the 1- and 2-directions
*EQUATION
 2,
 TOP_INT_FACE, 1, 1.0, BOT_INT_FACE, 1, 1.0
 2,
 TOP_INT_FACE, 2, 1.0, BOT_INT_FACE, 2, 1.0
** Fix edge nodes
*BOUNDARY
 EDGE_A, 1, 1
 EDGE_A, 3, 3
 EDGE_C, 1, 1
 EDGE_C, 3, 3
 EDGE_B, 2, 2
 EDGE_B, 3, 3
 EDGE_D, 2, 2
 EDGE_D, 3, 3
*EQUATION
 2,
 TOP_EDGE_A, 2, 1.0, BOT_EDGE_A, 2, 1.0
 2,
 TOP_EDGE_C, 2, 1.0, BOT_EDGE_C, 2, 1.0
 2,
 TOP_EDGE_B, 1, 1.0, BOT_EDGE_B, 1, 1.0
 2,
 TOP_EDGE_D, 1, 1.0, BOT_EDGE_D, 1, 1.0
** Remove thickness constraints on interior nodes
*EQUATION
 2,
 TOP_INT_FACE, 3, 1.0, BOT_INT_FACE, 3, -1.0
** **
** ** Frequency
** **
** *STEP,PERT
**  Frequency extraction
** *FREQUENCY, EIGENSOLVER=LANCZOS
**  12, , 
** *OUTPUT,FIELD
** *NODE OUTPUT
**  U,
** *END STEP
**
** Transient Load
** --------------
**
*STEP,AMP=STEP,INC=1000
*DYNAMIC,NOHAF
 0.0001,0.05
*DLOAD
 ALLE,P2,-1.E6
*NODE FILE,NSET=NMID
 U, 
*EL FILE,FRE=108,POSITION=AVERAGED AT NODES,ELSET=EMID
 S, 
*NODE PRINT,NSET=NMID
 U3, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID
 S11, 
*OUTPUT, HISTORY
*NODE OUTPUT, NSET=NMID
 U3,
*ELEMENT OUTPUT, ELSET=EMID
 S11,
*ENDSTEP
**
** Increase the damping and timestep to find steady state
**
*STEP,AMP=STEP,INC=1000
*DYNAMIC,NOHAF
 0.01,1.0
*DLOAD
 ALLE,P2,-1.E6
*NODE PRINT,NSET=NMID,F=999
 U3, 
*NODE FILE,NSET=NMID,F=999
 U, 
*OUTPUT, HISTORY
*NODE OUTPUT, NSET=NMID
 U3,
*ELEMENT OUTPUT, ELSET=EMID
 S11,
*ENDSTEP