*HEADING
*NODE
1, 0.5, 0.5, 0.
2, 0., 0.5, 0.5
3, 0.5, 0.5, 0.5
4, 0., -0.5, 0.5
5, 0.5, -0.5, 0.
6, 0.5, -0.5, 0.5
7, -0.5, -0.5, 1.
8, -0.5, -0.5, 0.5
9, -0.5, 0.5, 1.
10, -0.5, 0.5, 0.5
11, 0.5, 0.5, 1.
12, 0., 0.5, 1.
13, 0.5, -0.5, 1.
14, 0., -0.5, 1.
15, -0.5, 0.5, 0.
16, 0., 0.5, 0.
17, -0.5, -0.5, 0.
18, 0., -0.5, 0.
*ELEMENT, TYPE=C3D6, ELSET=CUBE
1, 3, 2, 1, 6, 4, 5
2, 8, 7, 4, 10, 9, 2
3, 12, 2, 11, 14, 4, 13
4, 16, 2, 15, 18, 4, 17
5, 14, 4, 7, 12, 2, 9
6, 10, 15, 2, 8, 17, 4
7, 18, 4, 5, 16, 2, 1
8, 6, 4, 13, 3, 2, 11
*ELSET,ELSET=EONE
2
*NSET,NSET=NONE
2
*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
3, DISTANGLE
*SOLID SECTION, ELSET=CUBE, MATERIAL=MATERIAL-1, ORIENTATION=ORI
1.,
*NSET, NSET=FIX
4, 5, 6, 7, 8, 11, 13, 14, 17, 18
*NSET, NSET=LOAD
1, 2, 3, 9, 10, 12, 15, 16
*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, PERTURBATION, NAME=STEP-0
DEFINE RESIDUAL MODES
*STATIC
*LOAD CASE, NAME=LC1
*CLOAD
LOAD, 2, 1.
*END LOAD CASE
*LOAD CASE, NAME=LC2
*CLOAD
LOAD, 3, 1.
*END LOAD CASE
*OUTPUT,FIELD,FREQUENCY=0
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-1
FREQUENCY ANALYSIS
*FREQUENCY,RESIDUAL MODES
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,120,4
*LOAD CASE,NAME=LC1B
*CLOAD
LOAD, 2, 1.E5
LOAD, 3, 1.E5
*END LOAD CASE
*LOAD CASE,NAME=LC2B
*CLOAD
LOAD, 2, 2.E5
LOAD, 3, 2.E5
*OUTPUT, FIELD
*NODE OUTPUT
U
*OUTPUT, HISTORY,FREQ=1
*ELEMENT OUTPUT,ELSET=EONE
S11,NFORC2
*END LOAD CASE
*LOAD CASE,NAME=LC3B
*CLOAD
LOAD, 2, 2.2E5
LOAD, 3, 2.2E5
*OUTPUT, FIELD
*ELEMENT OUTPUT
S,
*OUTPUT, HISTORY
*NODE OUTPUT,NSET=NONE
U1
*OUTPUT, HISTORY,FREQ=2
*ELEMENT OUTPUT,ELSET=EONE
NFORC2
*END LOAD CASE
*EL PRINT, FREQ=999
S,
*NODE PRINT, FREQ=0
*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
*EL PRINT, FREQ=999
S,
*NODE PRINT, FREQ=0
*OUTPUT, FIELD, OP=NEW
*ELEMENT OUTPUT
S,
*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