*HEADING MULTISTEP REFERENCE SOLUTION TEST FOR S4 ELEMENTS *PREPRINT,MODEL=YES *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, *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=S4, ELSET=SHELLSECL1 1, 1, 2, 6, 5, *SHELL GENERAL SECTION, ORIENT=ORIENTL1, COMPOSITE, ELSET=SHELLSECL1, OFFSET=-0.50000E+00, LAYUP=LAYUP-1 0.005, , MAT1, 90, PLY-1 0.005, , 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.5 **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL2 2, 2, 3, 7, 6, *SHELL GENERAL SECTION, ORIENT=ORIENTL2, COMPOSITE, ELSET=SHELLSECL2, OFFSET=-0.47500E+00, LAYUP=LAYUP-1 0.01, , MAT1, 90, PLY-1 0.01, , MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL2 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00 3,1.0 **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL3 3, 3, 4, 8, 7, *SHELL GENERAL SECTION, ORIENT=ORIENTL3, COMPOSITE, ELSET=SHELLSECL3, OFFSET=-0.45000E+00, LAYUP=LAYUP-1 0.015, , MAT1, 90, PLY-1 0.015, , MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL3 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00 3,2.0 **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL4 4, 5, 6, 10, 9, *SHELL GENERAL SECTION, ORIENT=ORIENTL4, COMPOSITE, ELSET=SHELLSECL4, OFFSET=-0.42500E+00, LAYUP=LAYUP-1 0.02, , MAT1, 90, PLY-1 0.02, , MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL4 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00 3,4.0 **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL5 5, 6, 7, 11, 10, *SHELL GENERAL SECTION, ORIENT=ORIENTL5, COMPOSITE, ELSET=SHELLSECL5, OFFSET=-0.40000E+00, LAYUP=LAYUP-1 0.025, , MAT1, 90, PLY-1 0.025, , MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL5 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00 3,8.0 **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL6 6, 7, 8, 12, 11, *SHELL GENERAL SECTION, ORIENT=ORIENTL6, COMPOSITE, ELSET=SHELLSECL6, OFFSET=-0.37500E+00, LAYUP=LAYUP-1 0.03, , MAT1, 90, PLY-1 0.03, , MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL6 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00 3,9.0 **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL7 7, 9, 10, 14, 13, *SHELL GENERAL SECTION, ORIENT=ORIENTL7, COMPOSITE, ELSET=SHELLSECL7, OFFSET=-0.35000E+00, LAYUP=LAYUP-1 0.035, , MAT1, 90, PLY-1 0.035, , MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL7 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00 3,10.0 **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL8 8, 10, 11, 15, 14, *SHELL GENERAL SECTION, ORIENT=ORIENTL8, COMPOSITE, ELSET=SHELLSECL8, OFFSET=-0.32500E+00, LAYUP=LAYUP-1 0.04, , MAT1, 90, PLY-1 0.04, , MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL8 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00 3,-8.0 **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S4, ELSET=SHELLSECL9 9, 11, 12, 16, 15, *SHELL GENERAL SECTION, ORIENT=ORIENTL9, COMPOSITE, ELSET=SHELLSECL9, OFFSET=-0.30000E+00, LAYUP=LAYUP-1 0.045, , MAT1, 90, PLY-1 0.045, , MAT1, 0, PLY-2 *ORIENTATION, NAME=ORIENTL9 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00 3,-10.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 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, 2, 2, 0.025 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, *OUTPUT, HISTORY *ELEMENT OUTPUT,ELSET=SHELLSECL1 S11, SF1, SSAVG1, SE1, SK1, STH, *ELEMENT OUTPUT,ELSET=SHELLSECL1 2, TSHR13, *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