*HEADING
:NAFEMS FORCED VIBRATION; TEST 13H;  [S3R]     
(DIRECT INTEGRATION) 
*RESTART,WRITE,FREQ=999
*NODE
       1,          0.,          0.
       2,        1.25,          0.
       3,         2.5,          0.
       4,        3.75,          0.
       5,          5.,          0.
       6,        6.25,          0.
       7,         7.5,          0.
       8,        8.75,          0.
       9,         10.,          0.
      10,          0.,        1.25
      11,        1.25,        1.25
      12,         2.5,        1.25
      13,        3.75,        1.25
      14,          5.,        1.25
      15,        6.25,        1.25
      16,         7.5,        1.25
      17,        8.75,        1.25
      18,         10.,        1.25
      19,          0.,         2.5
      20,        1.25,         2.5
      21,         2.5,         2.5
      22,        3.75,         2.5
      23,          5.,         2.5
      24,        6.25,         2.5
      25,         7.5,         2.5
      26,        8.75,         2.5
      27,         10.,         2.5
      28,          0.,        3.75
      29,        1.25,        3.75
      30,         2.5,        3.75
      31,        3.75,        3.75
      32,          5.,        3.75
      33,        6.25,        3.75
      34,         7.5,        3.75
      35,        8.75,        3.75
      36,         10.,        3.75
      37,          0.,          5.
      38,        1.25,          5.
      39,         2.5,          5.
      40,        3.75,          5.
      41,          5.,          5.
      42,        6.25,          5.
      43,         7.5,          5.
      44,        8.75,          5.
      45,         10.,          5.
      46,          0.,        6.25
      47,        1.25,        6.25
      48,         2.5,        6.25
      49,        3.75,        6.25
      50,          5.,        6.25
      51,        6.25,        6.25
      52,         7.5,        6.25
      53,        8.75,        6.25
      54,         10.,        6.25
      55,          0.,         7.5
      56,        1.25,         7.5
      57,         2.5,         7.5
      58,        3.75,         7.5
      59,          5.,         7.5
      60,        6.25,         7.5
      61,         7.5,         7.5
      62,        8.75,         7.5
      63,         10.,         7.5
      64,          0.,        8.75
      65,        1.25,        8.75
      66,         2.5,        8.75
      67,        3.75,        8.75
      68,          5.,        8.75
      69,        6.25,        8.75
      70,         7.5,        8.75
      71,        8.75,        8.75
      72,         10.,        8.75
      73,          0.,         10.
      74,        1.25,         10.
      75,         2.5,         10.
      76,        3.75,         10.
      77,          5.,         10.
      78,        6.25,         10.
      79,         7.5,         10.
      80,        8.75,         10.
      81,         10.,         10.
**
**
*ELEMENT, TYPE=S3R, ELSET=ALLE
       1,      11,      10,       1
       2,       1,       2,      11
       3,       2,       3,      11
       4,      12,      11,       3
       5,      13,      12,       3
       6,       3,       4,      13
       7,       4,       5,      13
       8,      14,      13,       5
       9,      15,      14,       5
      10,       5,       6,      15
      11,       6,       7,      15
      12,      16,      15,       7
      13,      17,      16,       7
      14,       7,       8,      17
      15,       8,       9,      17
      16,      18,      17,       9
      17,      10,      11,      19
      18,      20,      19,      11
      19,      21,      20,      11
      20,      11,      12,      21
      21,      12,      13,      21
      22,      22,      21,      13
      23,      23,      22,      13
      24,      13,      14,      23
      25,      14,      15,      23
      26,      24,      23,      15
      27,      25,      24,      15
      28,      15,      16,      25
      29,      16,      17,      25
      30,      26,      25,      17
      31,      27,      26,      17
      32,      17,      18,      27
      33,      29,      28,      19
      34,      19,      20,      29
      35,      20,      21,      29
      36,      30,      29,      21
      37,      31,      30,      21
      38,      21,      22,      31
      39,      22,      23,      31
      40,      32,      31,      23
      41,      33,      32,      23
      42,      23,      24,      33
      43,      24,      25,      33
      44,      34,      33,      25
      45,      35,      34,      25
      46,      25,      26,      35
      47,      26,      27,      35
      48,      36,      35,      27
      49,      28,      29,      37
      50,      38,      37,      29
      51,      39,      38,      29
      52,      29,      30,      39
      53,      30,      31,      39
      54,      40,      39,      31
      55,      41,      40,      31
      56,      31,      32,      41
      57,      32,      33,      41
      58,      42,      41,      33
      59,      43,      42,      33
      60,      33,      34,      43
      61,      34,      35,      43
      62,      44,      43,      35
      63,      45,      44,      35
      64,      35,      36,      45
      65,      47,      46,      37
      66,      37,      38,      47
      67,      38,      39,      47
      68,      48,      47,      39
      69,      49,      48,      39
      70,      39,      40,      49
      71,      40,      41,      49
      72,      50,      49,      41
      73,      51,      50,      41
      74,      41,      42,      51
      75,      42,      43,      51
      76,      52,      51,      43
      77,      53,      52,      43
      78,      43,      44,      53
      79,      44,      45,      53
      80,      54,      53,      45
      81,      46,      47,      55
      82,      56,      55,      47
      83,      57,      56,      47
      84,      47,      48,      57
      85,      48,      49,      57
      86,      58,      57,      49
      87,      59,      58,      49
      88,      49,      50,      59
      89,      50,      51,      59
      90,      60,      59,      51
      91,      61,      60,      51
      92,      51,      52,      61
      93,      52,      53,      61
      94,      62,      61,      53
      95,      63,      62,      53
      96,      53,      54,      63
      97,      65,      64,      55
      98,      55,      56,      65
      99,      56,      57,      65
     100,      66,      65,      57
     101,      67,      66,      57
     102,      57,      58,      67
     103,      58,      59,      67
     104,      68,      67,      59
     105,      69,      68,      59
     106,      59,      60,      69
     107,      60,      61,      69
     108,      70,      69,      61
     109,      71,      70,      61
     110,      61,      62,      71
     111,      62,      63,      71
     112,      72,      71,      63
     113,      64,      65,      73
     114,      74,      73,      65
     115,      75,      74,      65
     116,      65,      66,      75
     117,      66,      67,      75
     118,      76,      75,      67
     119,      77,      76,      67
     120,      67,      68,      77
     121,      68,      69,      77
     122,      78,      77,      69
     123,      79,      78,      69
     124,      69,      70,      79
     125,      70,      71,      79
     126,      80,      79,      71
     127,      81,      80,      71
     128,      71,      72,      81
**
** emid
**
*ELSET, ELSET=EMID
      55,      56,      57,      58,      71,      72,      73,       74
**
** nmid
**
*NSET, NSET=NMID
     ,   41
**
** right
**
*NSET, NSET=RIGHT
       9,      18,      27,      36,      45,      54,      63,      72,
     ,   81
**
** left
**
*NSET, NSET=LEFT
       1,      10,      19,      28,      37,      46,      55,      64,
     ,   73
**
** top
**
*NSET, NSET=TOP
      73,      74,      75,      76,      77,      78,      79,      80,
     ,   81
**
** bot
**
*NSET, NSET=BOT
       1,       2,       3,       4,       5,       6,       7,       8,
     ,    9
**
** edge
**
*NSET, NSET=EDGE
       1,       2,       3,       4,       5,       6,       7,       8,
       9,      10,      18,      19,      27,      28,      36,      37,
      45,      46,      54,      55,      63,      64,      72,      73,
      74,      75,      76,      77,      78,      79,      80,       81
**
** nall
**
*NSET, NSET=NALL
       1,       2,       3,       4,       5,       6,       7,       8,
       9,      10,      11,      12,      13,      14,      15,      16,
      17,      18,      19,      20,      21,      22,      23,      24,
      25,      26,      27,      28,      29,      30,      31,      32,
      33,      34,      35,      36,      37,      38,      39,      40,
      41,      42,      43,      44,      45,      46,      47,      48,
      49,      50,      51,      52,      53,      54,      55,      56,
      57,      58,      59,      60,      61,      62,      63,      64,
      65,      66,      67,      68,      69,      70,      71,      72,
      73,      74,      75,      76,      77,      78,      79,      80,
     ,   81
**
** all
**
*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
NALL,1,2
NALL,6
EDGE,1,3
BOT,5
TOP,5
LEFT,4
RIGHT,4
**
** Harmonic Load (direct)
** ----------------------
**
*STEP
*STEADY STATE DYNAMICS,DIRECT,FREQUENCY SCALE=LINEAR
0.1,4.16,200,,
*DLOAD,LOAD CASE=1
ALLE,P,100.
*NODE PRINT,NSET=NMID
U3, 
*EL PRINT,FREQUENCY=118,POSITION=AVERAGED AT NODES,ELSET=EMID
S11, 
*NODE FILE,NSET=NMID
U, 
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NMID
U,
*EL FILE,FREQUENCY=118,POSITION=AVERAGED AT NODES,ELSET=EMID
S, 
*ENDSTEP