*HEADING
MULTISTEP VERIFICATION TEST FOR CPE8R ELEMENTS USING DISTRIBUTIONS
*NODE
1, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00
2, 0.1666667E+00, 0.0000000E+00, 0.0000000E+00
3, 0.3333333E+00, 0.0000000E+00, 0.0000000E+00
4, 0.5000000E+00, 0.0000000E+00, 0.0000000E+00
5, 0.6666667E+00, 0.0000000E+00, 0.0000000E+00
6, 0.8333333E+00, 0.0000000E+00, 0.0000000E+00
7, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00
8, 0.0000000E+00, 0.1666667E+00, 0.0000000E+00
9, 0.1666667E+00, 0.1666667E+00, 0.0000000E+00
10, 0.3333333E+00, 0.1666667E+00, 0.0000000E+00
11, 0.5000000E+00, 0.1666667E+00, 0.0000000E+00
12, 0.6666667E+00, 0.1666667E+00, 0.0000000E+00
13, 0.8333333E+00, 0.1666667E+00, 0.0000000E+00
14, 0.1000000E+01, 0.1666667E+00, 0.0000000E+00
15, 0.0000000E+00, 0.3333333E+00, 0.0000000E+00
16, 0.1666667E+00, 0.3333333E+00, 0.0000000E+00
17, 0.3333333E+00, 0.3333333E+00, 0.0000000E+00
18, 0.5000000E+00, 0.3333333E+00, 0.0000000E+00
19, 0.6666667E+00, 0.3333333E+00, 0.0000000E+00
20, 0.8333333E+00, 0.3333333E+00, 0.0000000E+00
21, 0.1000000E+01, 0.3333333E+00, 0.0000000E+00
22, 0.0000000E+00, 0.5000000E+00, 0.0000000E+00
23, 0.1666667E+00, 0.5000000E+00, 0.0000000E+00
24, 0.3333333E+00, 0.5000000E+00, 0.0000000E+00
25, 0.5000000E+00, 0.5000000E+00, 0.0000000E+00
26, 0.6666667E+00, 0.5000000E+00, 0.0000000E+00
27, 0.8333333E+00, 0.5000000E+00, 0.0000000E+00
28, 0.1000000E+01, 0.5000000E+00, 0.0000000E+00
29, 0.0000000E+00, 0.6666667E+00, 0.0000000E+00
30, 0.1666667E+00, 0.6666667E+00, 0.0000000E+00
31, 0.3333333E+00, 0.6666667E+00, 0.0000000E+00
32, 0.5000000E+00, 0.6666667E+00, 0.0000000E+00
33, 0.6666667E+00, 0.6666667E+00, 0.0000000E+00
34, 0.8333333E+00, 0.6666667E+00, 0.0000000E+00
35, 0.1000000E+01, 0.6666667E+00, 0.0000000E+00
36, 0.0000000E+00, 0.8333333E+00, 0.0000000E+00
37, 0.1666667E+00, 0.8333333E+00, 0.0000000E+00
38, 0.3333333E+00, 0.8333333E+00, 0.0000000E+00
39, 0.5000000E+00, 0.8333333E+00, 0.0000000E+00
40, 0.6666667E+00, 0.8333333E+00, 0.0000000E+00
41, 0.8333333E+00, 0.8333333E+00, 0.0000000E+00
42, 0.1000000E+01, 0.8333333E+00, 0.0000000E+00
43, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00
44, 0.1666667E+00, 0.1000000E+01, 0.0000000E+00
45, 0.3333333E+00, 0.1000000E+01, 0.0000000E+00
46, 0.5000000E+00, 0.1000000E+01, 0.0000000E+00
47, 0.6666667E+00, 0.1000000E+01, 0.0000000E+00
48, 0.8333333E+00, 0.1000000E+01, 0.0000000E+00
49, 0.1000000E+01, 0.1000000E+01, 0.0000000E+00
*NSET, NSET=SIDE1, GENERATE
1, 7, 1,
*NSET, NSET=SIDE3, GENERATE
43, 49, 1,
*ELSET, ELSET=EALL, GENERATE
1, 9, 1,
**------------------------------------------------------------------------------------------
*MATERIAL, NAME=MAT1
*ELASTIC, TYPE=ENGINEERING CONSTANTS
0.300E+07, 0.400E+07, 0.500E+07, 0.100E+00, 0.200E+00, 0.300E+00, 0.300E+06, 0.400E+06
0.500E+06,
*DENSITY
0.750E-03,
**------------------------------------------------------------------------------------------
*ELEMENT, TYPE=CPE8R, ELSET=EALL
1, 1, 3, 17, 15, 2, 10, 16, 8,
2, 3, 5, 19, 17, 4, 12, 18, 10,
3, 5, 7, 21, 19, 6, 14, 20, 12,
4, 15, 17, 31, 29, 16, 24, 30, 22,
5, 17, 19, 33, 31, 18, 26, 32, 24,
6, 19, 21, 35, 33, 20, 28, 34, 26,
7, 29, 31, 45, 43, 30, 38, 44, 36,
8, 31, 33, 47, 45, 32, 40, 46, 38,
9, 33, 35, 49, 47, 34, 42, 48, 40,
*DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE
COORD3D, COORD3D
*DISTRIBUTION, NAME=DISTORIENT, LOCATION=ELEMENT, TABLE=ORIENTATIONTABLE
, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00,
1, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00,
2, 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00,
3, 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00,
4, 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00,
5, 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00,
6, 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00,
7, 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00,
8, 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00,
9, 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00,
*DISTRIBUTION TABLE, NAME=ANGLETABLE
ANGLE
*DISTRIBUTION, NAME=DISTANGLE, TABLE=ANGLETABLE, LOCATION=ELEMENT
, 0.0
1, 1.0
2, -1.0
3, 2.0
4, -2.0
5, 3.0
6, -3.0
7, 4.0
8, -4.0
9, 5.0
*ORIENTATION, NAME=ORI
DISTORIENT
3, DISTANGLE
*SOLID SECTION, MATERIAL=MAT1, ELSET=EALL, ORIENTATION=ORI
1.,
*BOUNDARY
SIDE1, 1, 2
*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.
**
** --------------------------------------------------
**
*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
SIDE3, 2, 100.
*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
SIDE3, 2, 100.
*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
SIDE3, 2, 100.
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-5
SSD SUBSPACE
*STEADY STATE DYNAMICS, SUBSPACE PROJECTION, INTERVAL=EIGENFREQUENCY
60.,
*CLOAD
SIDE3, 2, 100.
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-6
RANDOM RESPONSE
*RANDOM RESPONSE
60.,
*SELECT EIGENMODES, GENERATE
1, 5
*CLOAD, LOAD CASE=1
SIDE3, 2, 100.
*CLOAD, LOAD CASE=2
SIDE3, 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
SIDE3, 1, 1, 0.025
SIDE3, 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