*HEADING
 NAFEMS FORCED VIBRATION; TEST 13H;  [SC8R]     
 (DIRECT INTEGRATION) 
*PREPRINT, MODEL=YES, HISTORY=YES
*NODE
 90000001, 0.
 90000009, 10.
 90000081, 0., 10.
 90000089, 10.,10.
 90001001, 0. ,0.0, 0.05
 90001009, 10.,0.0, 0.05
 90001081, 0., 10., 0.05
 90001089, 10.,10., 0.05
*NGEN,NSET=EDGE_1_BOT
 90000001,90000009, 1
*NGEN,NSET=EDGE_1_TOP
 90001001,90001009, 1
*NGEN,NSET=EDGE_3_BOT
 90000081,90000089, 1
*NGEN,NSET=EDGE_3_TOP
 90001081,90001089, 1
*NFILL,NSET=NTOP
 EDGE_1_TOP,EDGE_3_TOP,8,10
*NFILL,NSET=NBOT
 EDGE_1_BOT,EDGE_3_BOT,8,10
*NSET,NSET=NALL
 NTOP,NBOT
*NSET,NSET=EDGE_4_BOT,GEN
 90000001, 90000081, 10
*NSET,NSET=EDGE_4_TOP,GEN
 90001001, 90001081, 10
*NSET,NSET=EDGE_2_BOT,GEN
 90000009, 90000089, 10
*NSET,NSET=EDGE_2_TOP,GEN
 90001009, 90001089, 10
*NSET,NSET=NEDGE_BOT
 EDGE_1_BOT,EDGE_2_BOT,EDGE_3_BOT,EDGE_4_BOT
*NSET,NSET=NEDGE_TOP
 EDGE_1_TOP,EDGE_2_TOP,EDGE_3_TOP,EDGE_4_TOP
*NSET,NSET=NMID
 90000045, 90001045
*NSET,NSET=NINTERIOR_BOT, GEN
 90000012, 90000018, 1
 90000022, 90000028, 1
 90000032, 90000038, 1
 90000042, 90000048, 1
 90000052, 90000058, 1
 90000062, 90000068, 1
 90000072, 90000078, 1
*NSET,NSET=NINTERIOR_TOP, GEN
 90001012, 90001018, 1
 90001022, 90001028, 1
 90001032, 90001038, 1
 90001042, 90001048, 1
 90001052, 90001058, 1
 90001062, 90001068, 1
 90001072, 90001078, 1
*ELEMENT,TYPE=SC8R
 100000000, 90000001, 90000002, 90000012, 90000011, 
     90001001, 90001002, 90001012, 90001011
*ELGEN,ELSET=ALLE
 100000000,8,1,1, 8,10,10
*ELSET,ELSET=EMID
 100000044,100000045,100000054,100000055
*SHELL SECTION,MATERIAL=MAT,ELSET=ALLE
  0.05, 
*MATERIAL,NAME=MAT
*DENSITY
 8000.,
*ELASTIC
 200.E9,0.3
*DAMPING,ALPHA=0.299,BETA=1.339E-3
*BOUNDARY
 EDGE_1_TOP, 1, 1
 EDGE_1_BOT, 1, 1
 EDGE_3_TOP, 1, 1
 EDGE_3_BOT, 1, 1
 EDGE_2_TOP, 2, 2
 EDGE_2_BOT, 2, 2
 EDGE_4_TOP, 2, 2
 EDGE_4_BOT, 2, 2
*EQUATION
 2,
 EDGE_1_TOP, 2, 1.0, EDGE_1_BOT, 2, 1.0
 2,
 EDGE_3_TOP, 2, 1.0, EDGE_3_BOT, 2, 1.0
 2,
 EDGE_2_TOP, 1, 1.0, EDGE_2_BOT, 1, 1.0
 2,
 EDGE_4_TOP, 1, 1.0, EDGE_4_BOT, 1, 1.0
 2,
 NINTERIOR_TOP, 1, 1.0, NINTERIOR_BOT, 1, 1.0
 2,
 NINTERIOR_TOP, 2, 1.0, NINTERIOR_BOT, 2, 1.0
 2,
 NINTERIOR_TOP, 3, 1.0, NINTERIOR_BOT, 3, -1.0
 2,
 NEDGE_TOP, 3,  1.0, NEDGE_BOT, 3, 1.0
**
** Harmonic Load (direct)
** ----------------------
**
*STEP
*STEADY STATE DYNAMICS,DIRECT,FREQUENCY SCALE=LINEAR
 0.1,4.16,200,,
*DLOAD,LOAD CASE=1
 ALLE, P2, -100.
*NODE PRINT,NSET=NMID
 U3, 
*EL PRINT,FREQUENCY=114,POSITION=AVERAGED AT NODES,ELSET=EMID
 S11, 
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=NMID
 U3,
*ELEMENT OUTPUT,ELSET=EMID
 S11,
*ENDSTEP