*HEADING
:NAFEMS FORCED VIBRATION; TEST 21H; [DEGENERATE S4]
(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=S4, ELSET=ALLE
1, 11, 10, 1, 11
2, 1, 2, 11, 1
3, 2, 3, 11, 2
4, 12, 11, 3, 12
5, 13, 12, 3, 13
6, 3, 4, 13, 3
7, 4, 5, 13, 4
8, 14, 13, 5, 14
9, 15, 14, 5, 15
10, 5, 6, 15, 5
11, 6, 7, 15, 6
12, 16, 15, 7, 16
13, 17, 16, 7, 17
14, 7, 8, 17, 7
15, 8, 9, 17, 8
16, 18, 17, 9, 18
17, 10, 11, 19, 10
18, 20, 19, 11, 20
19, 21, 20, 11, 21
20, 11, 12, 21, 11
21, 12, 13, 21, 12
22, 22, 21, 13, 22
23, 23, 22, 13, 23
24, 13, 14, 23, 13
25, 14, 15, 23, 14
26, 24, 23, 15, 24
27, 25, 24, 15, 25
28, 15, 16, 25, 15
29, 16, 17, 25, 16
30, 26, 25, 17, 26
31, 27, 26, 17, 27
32, 17, 18, 27, 17
33, 29, 28, 19, 29
34, 19, 20, 29, 19
35, 20, 21, 29, 20
36, 30, 29, 21, 30
37, 31, 30, 21, 31
38, 21, 22, 31, 21
39, 22, 23, 31, 22
40, 32, 31, 23, 32
41, 33, 32, 23, 33
42, 23, 24, 33, 23
43, 24, 25, 33, 24
44, 34, 33, 25, 34
45, 35, 34, 25, 35
46, 25, 26, 35, 25
47, 26, 27, 35, 26
48, 36, 35, 27, 36
49, 28, 29, 37, 28
50, 38, 37, 29, 38
51, 39, 38, 29, 39
52, 29, 30, 39, 29
53, 30, 31, 39, 30
54, 40, 39, 31, 40
55, 41, 40, 31, 41
56, 31, 32, 41, 31
57, 32, 33, 41, 32
58, 42, 41, 33, 42
59, 43, 42, 33, 43
60, 33, 34, 43, 33
61, 34, 35, 43, 34
62, 44, 43, 35, 44
63, 45, 44, 35, 45
64, 35, 36, 45, 35
65, 47, 46, 37, 47
66, 37, 38, 47, 37
67, 38, 39, 47, 38
68, 48, 47, 39, 48
69, 49, 48, 39, 49
70, 39, 40, 49, 39
71, 40, 41, 49, 40
72, 50, 49, 41, 50
73, 51, 50, 41, 51
74, 41, 42, 51, 41
75, 42, 43, 51, 42
76, 52, 51, 43, 52
77, 53, 52, 43, 53
78, 43, 44, 53, 43
79, 44, 45, 53, 44
80, 54, 53, 45, 54
81, 46, 47, 55, 46
82, 56, 55, 47, 56
83, 57, 56, 47, 57
84, 47, 48, 57, 47
85, 48, 49, 57, 48
86, 58, 57, 49, 58
87, 59, 58, 49, 59
88, 49, 50, 59, 49
89, 50, 51, 59, 50
90, 60, 59, 51, 60
91, 61, 60, 51, 61
92, 51, 52, 61, 51
93, 52, 53, 61, 52
94, 62, 61, 53, 62
95, 63, 62, 53, 63
96, 53, 54, 63, 53
97, 65, 64, 55, 65
98, 55, 56, 65, 55
99, 56, 57, 65, 56
100, 66, 65, 57, 66
101, 67, 66, 57, 67
102, 57, 58, 67, 57
103, 58, 59, 67, 58
104, 68, 67, 59, 68
105, 69, 68, 59, 69
106, 59, 60, 69, 59
107, 60, 61, 69, 60
108, 70, 69, 61, 70
109, 71, 70, 61, 71
110, 61, 62, 71, 61
111, 62, 63, 71, 62
112, 72, 71, 63, 72
113, 64, 65, 73, 64
114, 74, 73, 65, 74
115, 75, 74, 65, 75
116, 65, 66, 75, 65
117, 66, 67, 75, 66
118, 76, 75, 67, 76
119, 77, 76, 67, 77
120, 67, 68, 77, 67
121, 68, 69, 77, 68
122, 78, 77, 69, 78
123, 79, 78, 69, 79
124, 69, 70, 79, 69
125, 70, 71, 79, 70
126, 80, 79, 71, 80
127, 81, 80, 71, 81
128, 71, 72, 81, 71
**
** 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
1.,
*MATERIAL,NAME=MAT
*DENSITY
8000.,
*ELASTIC
200.E9,0.3
*DAMPING,ALPHA=5.772, BETA=6.929E-5
*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
1.,78.17,200,,
*DLOAD,LOAD CASE=1
ALLE,P,1.E6
*NODE PRINT,NSET=NMID
U3,
*EL FILE,FRE=122,POSITION=AVERAGED AT NODES,ELSET=EMID
S,
*NODE FILE,NSET=NMID
U,
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID,freq=122
S,
*ENDSTEP