*HEADING
*PREPRINT,MODEL=YES
*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.
*ELEMENT, TYPE=C3D8, ELSET=CUBE
1, 5, 6, 7, 8, 1, 2, 3, 4
2, 6, 11, 12, 7, 2, 9, 10, 3
3, 6, 13, 14, 7, 11, 15, 16, 12
4, 21, 5, 19, 20, 13, 6, 17, 18
5, 5, 6, 13, 21, 8, 7, 14, 22
6, 11, 6, 13, 15, 23, 17, 18, 24
7, 11, 6, 17, 23, 9, 2, 25, 26
8, 6, 5, 19, 17, 2, 1, 27, 25
*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
,0
1, 0.0
2, 1.0
3, 2.0
4, 3.0
5, 4.0
6, -1.0
7, -2.0
8, -3.0
*ORIENTATION, NAME=ORI
DISTORI
1, DISTANGLE
*SOLID SECTION, ELSET=CUBE, MATERIAL=MATERIAL-1, ORIENTATION=ORI
1.,
*NSET, NSET=FIX
4, 17, 18, 19, 20, 23, 24, 25, 26, 27
*NSET, NSET=LOAD
3, 7, 8, 10, 12, 14, 16, 22
*ELSET, ELSET=EALL, GENERATE
1, 8, 1
**
** MATERIALS
**
*DISTRIBUTION TABLE,NAME=DISTTABLEORTHO
MODULUS,MODULUS,MODULUS,MODULUS,MODULUS,MODULUS,MODULUS,
MODULUS,MODULUS
*DISTRIBUTION, NAME=DISTORTHO, LOCATION=ELEMENT, TABLE=DISTTABLEORTHO
,4.03846153846154E6,1.73076923076923E6,4.03846153846154E6,1.73076923076923E6,1.73076923076923E6,4.03846153846154E6,1.15384615384615E6
1.15384615384615E6,1.15384615384615E6
CUBE,3.38095238095238E6,0.88888888888889E6,4.74074074074074E6,1.46031746031746E6,2.07407407407407E6,6.26455026455027E6,3.E5,
4.E5,5.E5
*MATERIAL, NAME=MATERIAL-1
*ELASTIC, TYPE=ORTHOTROPIC
DISTORTHO
*DENSITY
0.750E-03,
*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