*Heading
*Node
1, 0.5, 0., 1.
2, 0.5, 0., 0.5
3, 0.5, 0.5, 0.5
4, 0.5, 0.5, 1.
5, 0., 0., 1.
6, 0., 0., 0.5
7, 0., 0.5, 0.5
8, 0., 0.5, 1.
9, 0.5, 0., 0.
10, 0.5, 0.5, 0.
11, 0., 0., 0.
12, 0., 0.5, 0.
13, -0.5, 0., 0.5
14, -0.5, 0.5, 0.5
15, -0.5, 0., 0.
16, -0.5, 0.5, 0.
17, 0., -0.5, 0.5
18, -0.5, -0.5, 0.5
19, 0., -0.5, 1.
20, -0.5, -0.5, 1.
21, -0.5, 0., 1.
22, -0.5, 0.5, 1.
23, 0., -0.5, 0.
24, -0.5, -0.5, 0.
25, 0.5, -0.5, 0.5
26, 0.5, -0.5, 0.
27, 0.5, -0.5, 1.
28, 0., 0.25, 1.
29, 0., 0.5, 0.75
30, 0., 0.25, 0.5
31, 0., 0., 0.75
32, 0.5, 0., 0.75
33, 0.5, 0.25, 0.5
34, 0.5, 0.5, 0.75
35, 0.5, 0.25, 1.
36, 0.25, 0., 0.5
37, 0.25, 0., 1.
38, 0.25, 0.5, 0.5
39, 0.25, 0.5, 1.
40, 0., 0.5, 0.25
41, 0., 0.25, 0.
42, 0., 0., 0.25
43, 0.5, 0., 0.25
44, 0.5, 0.25, 0.
45, 0.5, 0.5, 0.25
46, 0.25, 0., 0.
47, 0.25, 0.5, 0.
48, -0.25, 0.5, 0.5
49, -0.5, 0.25, 0.5
50, -0.25, 0., 0.5
51, -0.25, 0., 0.
52, -0.5, 0.25, 0.
53, -0.25, 0.5, 0.
54, -0.5, 0., 0.25
55, -0.5, 0.5, 0.25
56, -0.5, -0.25, 1.
57, -0.25, -0.5, 1.
58, 0., -0.25, 1.
59, -0.25, 0., 1.
60, 0., -0.25, 0.5
61, -0.25, -0.5, 0.5
62, -0.5, -0.25, 0.5
63, -0.5, 0., 0.75
64, 0., -0.5, 0.75
65, -0.5, -0.5, 0.75
66, -0.5, 0.5, 0.75
67, -0.25, 0.5, 1.
68, -0.5, 0.25, 1.
69, 0., -0.5, 0.25
70, -0.5, -0.5, 0.25
71, -0.25, -0.5, 0.
72, 0., -0.25, 0.
73, -0.5, -0.25, 0.
74, 0.5, -0.25, 0.5
75, 0.5, -0.5, 0.25
76, 0.5, -0.25, 0.
77, 0.25, -0.5, 0.5
78, 0.25, -0.5, 0.
79, 0.5, -0.25, 1.
80, 0.5, -0.5, 0.75
81, 0.25, -0.5, 1.
*Element, type=C3D20R, elset=cube
1, 5, 6, 7, 8, 1, 2, 3, 4, 31, 30, 29, 28, 32, 33, 34,
35, 37, 36, 38, 39
2, 6, 11, 12, 7, 2, 9, 10, 3, 42, 41, 40, 30, 43, 44, 45,
33, 36, 46, 47, 38
3, 6, 13, 14, 7, 11, 15, 16, 12, 50, 49, 48, 30, 51, 52, 53,
41, 42, 54, 55, 40
4, 21, 5, 19, 20, 13, 6, 17, 18, 59, 58, 57, 56, 50, 60, 61,
62, 63, 31, 64, 65
5, 5, 6, 13, 21, 8, 7, 14, 22, 31, 50, 63, 59, 29, 48, 66,
67, 28, 30, 49, 68
6, 11, 6, 13, 15, 23, 17, 18, 24, 42, 50, 54, 51, 69, 61, 70,
71, 72, 60, 62, 73
7, 11, 6, 17, 23, 9, 2, 25, 26, 42, 60, 69, 72, 43, 74, 75,
76, 46, 36, 77, 78
8, 6, 5, 19, 17, 2, 1, 27, 25, 31, 58, 64, 60, 32, 79, 80,
74, 36, 37, 81, 77
*Solid Section, elset=cube, material=Material-1, orientation=ori
1.,
*Distribution Table, name=distTableOri
coord3d, coord3d
*Distribution, name=distOri, table=distTableOri, location=element
, 1.,0.,0., 0.,1.,0.
1, 1., 1., 1., -1., 1.,-1.
2, 1., 1.,-1., -1., 1., 1.
3, -1., 1.,-1., 1., 1., 1.
4, -1.,-1., 1., 1.,-1.,-1.
5, -1., 1., 1., 1., 1.,-1.
6, -1.,-1.,-1., 1.,-1., 1.
7, 1.,-1.,-1., -1.,-1., 1.
8, 1.,-1., 1., -1.,-1.,-1.
*orientation,name=ori
distOri
3, 0.
*Nset, nset=Fix
4, 17, 18, 19, 20, 23, 24, 25, 26, 27, 57, 61, 64, 65, 69, 70
71, 75, 77, 78, 80, 81
*Nset, nset=Load
3,
** 3, 7, 8, 10, 12, 14, 16, 22, 29, 34, 38, 39, 40, 45, 47
** 48, 53, 55, 66, 67
*Material, name=Material-1
*Elastic, type=ENGINEERING CONSTANTS
3e+06, 4e+06, 5e+06, 0.1, 0.2, 0.3,300000.,400000.
500000.,
*Density
0.750E-3
*Elset, elset=eall, generate
1, 8, 1
*Boundary
Fix, 1,3
**
**
**
*step, perturbation, name=step-1
perturbation step for residual modes
*static
*load case, name=step3
*cload
load, 1, 1.e5
*end load case
*load case, name=step4
*cload
load, 2, 1.e5
*end load case
*load case, name=step5to7
*cload
load, 3, 1.e5
*end load case
*end step
**
**
*step, name=freq
frequency analysis with residual modes
*frequency, eigensolver=lanczos, normalization=mass, residual modes
, , 14000.
*output, field
*node output
U,
*end step
**
** --------------------------------------------------
**
*step, name=step-3
SSD mode-based
*steady state dynamics, interval=eigenfrequency
13000,
*cload
load, 2, 1.e5
load, 3, 1.e5
*output, field
*node output
U,
*end step