*HEADING MULTISTEP REFERENCE SOLUTION 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=EALL, GENERATE 1, 9, 1, *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=SHELLSECL1 1, 1, 2, 6, 5, 4, 5, 6, 10, 9, 7, 9, 10, 14, 13, *SHELL SECTION, ORIENT=ORIENTL1, COMPOSITE, ELSET=SHELLSECL1, OFFSET=-0.50000E+00, NODAL THICKNESS 1.0, 3, MAT1, 90, PLY-1 1.0, 3, MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL1 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00 3,0. **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL2 2, 2, 3, 7, 6, 5, 6, 7, 11, 10, 8, 10, 11, 15, 14, *SHELL SECTION, ORIENT=ORIENTL2, COMPOSITE, ELSET=SHELLSECL2, OFFSET=-0.50000E+00 0.005, 3, MAT1, 90, PLY-1 0.005, 3, MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL2 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, -0.1000000E+01, 0.0000000E+00, 0.0000000E+00 3,0. **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL3 3, 3, 4, 8, 7, 6, 7, 8, 12, 11, 9, 11, 12, 16, 15, *SHELL SECTION, ORIENT=ORIENTL3, COMPOSITE, ELSET=SHELLSECL3, OFFSET=-0.50000E+00, NODAL THICKNESS 1.0, 3, MAT1, 90, PLY-1 1.0, 3, MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL3 -0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, -0.1000000E+01, 0.0000000E+00 3,0. *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