*HEADING
MULTISTEP VERIFICATION TEST FOR S3R ELEMENTS WITH NODAL THICKNESS
*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.0000000E+00, 0.1666667E+00, 0.0000000E+00
6, 0.1666667E+00, 0.1666667E+00, 0.0000000E+00
7, 0.3333333E+00, 0.1666667E+00, 0.0000000E+00
8, 0.5000000E+00, 0.1666667E+00, 0.0000000E+00
9, 0.0000000E+00, 0.3333333E+00, 0.0000000E+00
10, 0.1666667E+00, 0.3333333E+00, 0.0000000E+00
11, 0.3333333E+00, 0.3333333E+00, 0.0000000E+00
12, 0.5000000E+00, 0.3333333E+00, 0.0000000E+00
13, 0.0000000E+00, 0.5000000E+00, 0.0000000E+00
14, 0.1666667E+00, 0.5000000E+00, 0.0000000E+00
15, 0.3333333E+00, 0.5000000E+00, 0.0000000E+00
16, 0.5000000E+00, 0.5000000E+00, 0.0000000E+00
*NSET, NSET=SIDE1, GENERATE
1, 4, 1,
*NSET, NSET=SIDE3, GENERATE
13, 16, 1,
*NSET, NSET=NODETHICK1, GENERATE
1, 13, 4,
4, 16, 4,
*NSET, NSET=NODETHICK2, GENERATE
2, 14, 4,
3, 15, 4,
*ELSET, ELSET=SHELLSECL1, GENERATE
1, 13, 6,
2, 14, 6,
*ELSET, ELSET=SHELLSECL2, GENERATE
3, 15, 6,
4, 16, 6,
*ELSET, ELSET=SHELLSECL3, GENERATE
5, 17, 6,
6, 18, 6,
*ELSET, ELSET=SIDEEL
SHELLSECL1, SHELLSECL3
*ELSET, ELSET=MIDEL
SHELLSECL2,
*NODAL THICKNESS
NODETHICK1, 0.01
NODETHICK2, 0.02
**------------------------------------------------------------------------------------------
*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,
*MATERIAL, NAME=MAT2
*ELASTIC, TYPE=ENGINEERING CONSTANTS
0.600E+07, 0.200E+07, 0.200E+07, 0.100E+00, 0.200E+00, 0.300E+00, 0.300E+06, 0.400E+06
0.200E+06,
*DENSITY
0.700E-03,
**------------------------------------------------------------------------------------------
*ELEMENT, TYPE=S3R, ELSET=EALL
1, 1, 2, 6,
2, 6, 5, 1,
3, 2, 3, 7,
4, 7, 6, 2,
5, 3, 4, 8,
6, 8, 7, 3,
7, 5, 6, 10,
8, 10, 9, 5,
9, 6, 7, 11,
10, 11, 10, 6,
11, 7, 8, 12,
12, 12, 11, 7,
13, 9, 10, 14,
14, 14, 13, 9,
15, 10, 11, 15,
16, 15, 14, 10,
17, 11, 12, 16,
18, 16, 15, 11,
*DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE
COORD3D, COORD3D
*DISTRIBUTION, NAME=DISTORIENT, LOCATION=ELEMENT, TABLE=ORIENTATIONTABLE
, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
SHELLSECL1, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00
SHELLSECL2, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, -0.1000000E+01, 0.0000000E+00, 0.0000000E+00
SHELLSECL3, -0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, -0.1000000E+01, 0.0000000E+00
*DISTRIBUTION TABLE, NAME=ORIANGLETABLE
ANGLE
*DISTRIBUTION, NAME=DISTORIANGLE, LOCATION=ELEMENT, TABLE=ORIANGLETABLE
, 0
SHELLSECL1, 10
SHELLSECL2, 15
SHELLSECL3, 20
*DISTRIBUTION TABLE, NAME=THICKNESSTABLE
LENGTH
*DISTRIBUTION, NAME=DISTTHICK, LOCATION=ELEMENT, TABLE=THICKNESSTABLE
, 1.0
SHELLSECL2, 0.01
*ORIENTATION, NAME=ORI
DISTORIENT
3, DISTORIANGLE
*SHELL GENERAL SECTION, COMPOSITE, ELSET=MIDEL, OFFSET=-0.50000E+00, SHELL THICKNESS=DISTTHICK, ORIENTATION=ORI
1.0, , MAT1, -30, PLY-1
1.0, , MAT1, +45, PLY-2
1.0, , MAT2, -45, PLY-3
1.0, , MAT2, +45, PLY-4
*SHELL GENERAL SECTION, COMPOSITE, ELSET=SIDEEL, NODAL THICKNESS, OFFSET=-0.5, ORIENTATION=ORI
1.0, , MAT1, -30, PLY-1
1.0, , MAT1, +45, PLY-2
1.0, , MAT2, -45, PLY-3
1.0, , MAT2, +45, PLY-4
*BOUNDARY
SIDE1, 1, 6
*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,
SF,
SSAVG,
SE,
*EL PRINT
2,
TSHR,
*NODE PRINT
U,
*OUTPUT, FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
SF,
SSAVG,
SE,
STH,
*ELEMENT OUTPUT
2,
TSHR,
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-2
SSD MODE-BASED
*STEADY STATE DYNAMICS, INTERVAL=EIGENFREQUENCY
60,
*EL PRINT, FREQ=999
S,
SF,
SSAVG,
SE,
STH,
*EL PRINT, FREQ=999
2,
TSHR,
*NODE PRINT, FREQ=999
U,
*OUTPUT, FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
SF,
SSAVG,
SE,
STH,
*ELEMENT OUTPUT
2,
TSHR,
*CLOAD
SIDE3, 4, 0.1
SIDE3, 3, 2.0
SIDE3, 2, 0.5
*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, 4, 0.1
SIDE3, 3, 0.2
SIDE3, 2, 0.5
*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, 4, 0.1
SIDE3, 3, 0.2
SIDE3, 2, 0.5
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-5
SSD SUBSPACE
*STEADY STATE DYNAMICS, SUBSPACE PROJECTION, INTERVAL=EIGENFREQUENCY
60.,
*CLOAD
SIDE3, 4, 0.1
SIDE3, 3, 0.2
SIDE3, 2, 0.5
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-6
RANDOM RESPONSE
*RANDOM RESPONSE
60.,
*SELECT EIGENMODES, GENERATE
1, 5
*CLOAD, LOAD CASE=1
SIDE3, 4, 0.1
SIDE3, 3, 0.2
SIDE3, 2, 0.5
*CLOAD, LOAD CASE=2
SIDE3, 4, 0.1
SIDE3, 3, 0.2
SIDE3, 2, 0.5
*CORRELATION, PSD=WHITE, TYPE=UNCORRELATED, COMPLEX=NO
1,1
2,1
*EL PRINT, FREQ=999
SF,
SE,
*NODE PRINT, FREQ=999
U,
*OUTPUT, FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
SF,
SE,
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-7
RESPONSE SPECTRUM
*RESPONSE SPECTRUM, SUM=ABS, COMP=ALG
SPEC,0.,0.,1.,1.
*SELECT EIGENMODES, GENERATE
1, 5
*END STEP
**
** --------------------------------------------------
**
*STEP, NLGEOM=YES, NAME=STEP-8
STATIC
*STATIC
1.0, 1.0
*BOUNDARY
SIDE3, 3, 3, 0.050
*EL PRINT, FREQ=999
S,
SF,
SSAVG,
SE,
STH,
*EL PRINT, FREQ=999
2,
TSHR,
*NODE PRINT, FREQ=999
U,
*OUTPUT, FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
SF,
SSAVG,
SE,
STH,
*ELEMENT OUTPUT
2,
TSHR,
*END STEP
**
** --------------------------------------------------
**
*STEP, PERTURBATION, NAME=STEP-9
STATIC
*STATIC
*EL PRINT, FREQ=999
S,
SF,
SSAVG,
SE,
STH,
*EL PRINT, FREQ=999
2,
TSHR,
*NODE PRINT, FREQ=999
U,
*OUTPUT, FIELD
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
SF,
SSAVG,
SE,
STH,
*ELEMENT OUTPUT
2,
TSHR,
*LOAD CASE, NAME=LOAD1
*DLOAD
EALL, BZ, 1.E3
*END LOAD CASE
*LOAD CASE, NAME=LOAD2
*DLOAD
EALL, GRAV, 1.E6, 0.,0.,1.
*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