*HEADING MULTISTEP VERIFICATION TEST FOR S4 ELEMENTS WITH NODAL THICKNESS *NODE 1, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00 2, 0.3333333E+00, 0.0000000E+00, 0.0000000E+00 3, 0.6666667E+00, 0.0000000E+00, 0.0000000E+00 4, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00 5, 0.0000000E+00, 0.3333333E+00, 0.0000000E+00 6, 0.3333333E+00, 0.3333333E+00, 0.0000000E+00 7, 0.6666667E+00, 0.3333333E+00, 0.0000000E+00 8, 0.1000000E+01, 0.3333333E+00, 0.0000000E+00 9, 0.0000000E+00, 0.6666667E+00, 0.0000000E+00 10, 0.3333333E+00, 0.6666667E+00, 0.0000000E+00 11, 0.6666667E+00, 0.6666667E+00, 0.0000000E+00 12, 0.1000000E+01, 0.6666667E+00, 0.0000000E+00 13, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00 14, 0.3333333E+00, 0.1000000E+01, 0.0000000E+00 15, 0.6666667E+00, 0.1000000E+01, 0.0000000E+00 16, 0.1000000E+01, 0.1000000E+01, 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, 7, 3, *ELSET, ELSET=SHELLSECL2, GENERATE 2, 8, 3, *ELSET, ELSET=SHELLSECL3, GENERATE 3, 9, 3, *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, **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=EALL 1, 1, 2, 6, 5, 2, 2, 3, 7, 6, 3, 3, 4, 8, 7, 4, 5, 6, 10, 9, 5, 6, 7, 11, 10, 6, 7, 8, 12, 11, 7, 9, 10, 14, 13, 8, 10, 11, 15, 14, 9, 11, 12, 16, 15, *DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE COORD3D, COORD3D *DISTRIBUTION, NAME=DISTORIENT, LOCATION=ELEMENT, TABLE=ORIENTATIONTABLE , 1., 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 *ORIENTATION, NAME=ORI DISTORIENT 3, 0.0 *DISTRIBUTION TABLE, NAME=THICKNESSTABLE LENGTH *DISTRIBUTION, NAME=DISTTHICK, LOCATION=ELEMENT, TABLE=THICKNESSTABLE , 1.0 SHELLSECL2, 0.01 *SHELL SECTION, COMPOSITE, ELSET=MIDEL, OFFSET=-0.50000E+00, ORIENT=ORI, SHELL THICKNESS=DISTTHICK 1.0, 3, MAT1, 90, PLY-1 1.0, 3, MAT1, 0, PLY-2 *SHELL SECTION, COMPOSITE, ELSET=SIDEEL, NODAL THICKNESS, OFFSET=-0.5, ORIENTATION=ORI 1.0, 3, MAT1, 90, PLY-1 1.0, 3, MAT1, 0, PLY-2 *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 3, TSHR, *NODE PRINT U, *OUTPUT, FIELD *NODE OUTPUT U, *ELEMENT OUTPUT S, SF, SSAVG, SE, STH, *ELEMENT OUTPUT 3, 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 3, TSHR, *NODE PRINT, FREQ=999 U, *OUTPUT, FIELD *NODE OUTPUT U, *ELEMENT OUTPUT S, SF, SSAVG, SE, STH, *ELEMENT OUTPUT 3, 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 3, TSHR, *NODE PRINT, FREQ=999 U, *OUTPUT, FIELD *NODE OUTPUT U, *ELEMENT OUTPUT S, SF, SSAVG, SE, STH, *ELEMENT OUTPUT 3, TSHR, *OUTPUT, HISTORY *ELEMENT OUTPUT,ELSET=SHELLSECL1 S11, SF1, SSAVG1, SE1, SK1, STH, *ELEMENT OUTPUT,ELSET=SHELLSECL1 3, TSHR13, *END STEP ** ** -------------------------------------------------- ** *STEP, PERTURBATION, NAME=STEP-9 STATIC *STATIC *EL PRINT, FREQ=999 S, SF, SSAVG, SE, STH, *EL PRINT, FREQ=999 3, TSHR, *NODE PRINT, FREQ=999 U, *OUTPUT, FIELD *NODE OUTPUT U, *ELEMENT OUTPUT S, SF, SSAVG, SE, STH, *ELEMENT OUTPUT 3, 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