*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
*DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE
COORD3D, COORD3D
*DISTRIBUTION, NAME=DISTORI, TABLE=ORIENTATIONTABLE, 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.
*DISTRIBUTION TABLE, NAME=ANGLETABLE
ANGLE
*DISTRIBUTION, NAME=DISTANGLE, TABLE=ANGLETABLE, LOCATION=ELEMENT
,1.0
1, 1.0
2, 2.0
3, 3.0
4, 4.0
5, 5.0
6, 6.0
7, 7.0
8, 8.0
*ORIENTATION, NAME=ORI
DISTORI
3, DISTANGLE
*SOLID SECTION, ELSET=CUBE, MATERIAL=MATERIAL-1, ORIENTATION=ORI
1.,
*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, 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
*PSD-DEFINITION, NAME=WHITE, TYPE=FORCE
1.0,0.1,0.0
1.0,0.1,60.0
*SPECTRUM, NAME=SPEC,TYPE=ACCELERATION
0., 0., 0.
1.E5, 10., 0.
1.E5, 100., 0.
*BOUNDARY
FIX, 1,3
**
** --------------------------------------------------
**
*STEP, NAME=STEP-1
FREQUENCY ANALYSIS
*FREQUENCY,SIM,normalization=mass
5,
*EL PRINT
S,
*NODE PRINT
U,
*OUTPUT, FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-2
SSD MODE-BASED
*STEADY STATE DYNAMICS, INTERVAL=EIGENFREQUENCY
60,
*EL PRINT, FREQ=999
S,
*NODE PRINT, FREQ=0
*OUTPUT, FIELD, OP=NEW
*ELEMENT OUTPUT
S,
*CLOAD
LOAD, 2, 1.E5
LOAD, 3, 1.E5
*SELECT EIGENMODES, GENERATE
1, 5
*MODAL DAMPING, RAYLEIGH
1, 5, 0.5
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-3
MODAL DYNAMIC
*MODAL DYNAMIC
0.1,1.
*CLOAD
LOAD, 2, 1.E5
LOAD, 3, 1.E5
*SELECT EIGENMODES, GENERATE
1, 5
*MODAL DAMPING, RAYLEIGH
1, 5, 0.5
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-4
SSD DIRECT
*STEADY STATE DYNAMICS, DIRECT, INTERVAL=EIGENFREQUENCY
60.,
*CLOAD
LOAD, 2, 100.
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-5
SSD SUBSPACE
*STEADY STATE DYNAMICS, SUBSPACE PROJECTION, INTERVAL=EIGENFREQUENCY
60.,
*CLOAD
LOAD, 2, 100.
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-6
RANDOM RESPONSE
*RANDOM RESPONSE
60.,
*SELECT EIGENMODES, GENERATE
1, 5
*CLOAD, LOAD CASE=1
LOAD, 2, 100.
*CLOAD, LOAD CASE=2
LOAD, 2, 100.
*CORRELATION, PSD=WHITE, TYPE=UNCORRELATED, COMPLEX=NO
1,1
2,1
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-7
RESPONSE SPECTRUM
*RESPONSE SPECTRUM, SUM=ABS, COMP=ALG
SPEC,0.,1.,0.,1.
*SELECT EIGENMODES, GENERATE
1, 5
*END STEP
**
** --------------------------------------------------
**
*STEP, NLGEOM=YES, NAME=STEP-8
STATIC
*STATIC
1.0, 1.0
*BOUNDARY
LOAD, 1, 1, 0.025
LOAD, 2, 2, 0.050
*EL PRINT, FREQ=999
S,
*NODE PRINT, FREQ=999
U,
*OUTPUT, FIELD, OP=NEW
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
*END STEP
**
** --------------------------------------------------
**
*STEP, PERTURBATION, NAME=STEP-9
STATIC
*STATIC
*EL PRINT, FREQ=999
S,
*NODE PRINT, FREQ=0
*OUTPUT, FIELD, OP=NEW
*ELEMENT OUTPUT
S,
*LOAD CASE, NAME=LOAD1
*DLOAD
EALL, BX, 1.E3
*END LOAD CASE
*LOAD CASE, NAME=LOAD2
*DLOAD
EALL, GRAV, 1.E6, 0.,1.,0.
*END LOAD CASE
*LOAD CASE, NAME=LOAD3
*DLOAD
EALL, CENT, 100.,0.,0.,0.,1.,0.
*END LOAD CASE
*LOAD CASE, NAME=LOAD4
*DLOAD
EALL, CENTRIF, 1.E6,0.,0.,0.,1.,0.
*END LOAD CASE
*END STEP