*HEADING MULTISTEP VERIFICATION TEST FOR S3 ELEMENTS USING DISTRIBUTIONS *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, 18, 1, *ELSET, ELSET=SHELLSECL1 1, **------------------------------------------------------------------------------------------ *MATERIAL, NAME=MAT1 *ELASTIC, TYPE=ENGINEERING CONSTANTS 0.310E+07, 0.410E+07, 0.510E+07, 0.110E+00, 0.210E+00, 0.310E+00, 0.310E+06, 0.410E+06 0.510E+06, *DENSITY 0.751E-03, *MATERIAL, NAME=MAT2 *ELASTIC, TYPE=ENGINEERING CONSTANTS 0.320E+07, 0.420E+07, 0.520E+07, 0.120E+00, 0.220E+00, 0.320E+00, 0.320E+06, 0.420E+06 0.520E+06, *DENSITY 0.752E-03, *MATERIAL, NAME=MAT3 *ELASTIC, TYPE=ENGINEERING CONSTANTS 0.330E+07, 0.430E+07, 0.530E+07, 0.130E+00, 0.230E+00, 0.330E+00, 0.330E+06, 0.430E+06 0.530E+06, *DENSITY 0.753E-03, *MATERIAL, NAME=MAT4 *ELASTIC, TYPE=ENGINEERING CONSTANTS 0.340E+07, 0.440E+07, 0.540E+07, 0.140E+00, 0.240E+00, 0.340E+00, 0.340E+06, 0.440E+06 0.540E+06, *DENSITY 0.754E-03, **------------------------------------------------------------------------------------------ *ELEMENT, TYPE=S3, ELSET=E1 1, 1, 2, 6, 2, 6, 5, 1, *ELEMENT, TYPE=S3, ELSET=E2 3, 2, 3, 7, 4, 7, 6, 2, *ELEMENT, TYPE=S3, ELSET=E3 5, 3, 4, 8, 6, 8, 7, 3, *ELEMENT, TYPE=S3, ELSET=E4 7, 5, 6, 10, 8, 10, 9, 5, *ELEMENT, TYPE=S3, ELSET=E5 9, 6, 7, 11, 10, 11, 10, 6, *ELEMENT, TYPE=S3, ELSET=E6 11, 7, 8, 12, 12, 12, 11, 7, *ELEMENT, TYPE=S3, ELSET=E7 13, 9, 10, 14, 14, 14, 13, 9, *ELEMENT, TYPE=S3, ELSET=E8 15, 10, 11, 15, 16, 15, 14, 10, *ELEMENT, TYPE=S3, ELSET=E9 17, 11, 12, 16, 18, 16, 15, 11, *SHELL SECTION, COMPOSITE, ELSET=EALL, OFFSET=DISTOFFSET, SHELL THICKNESS=DISTTHICKSECT, ORIENTATION=ORI,LAYUP=LAYUP-1 DISTTHICKPLY1, 3, MAT1, DISTANGLEPLY1, PLY-1 DISTTHICKPLY2, 3, MAT2, -60, PLY-2 .01253, 3, MAT3, DISTANGLEPLY3, PLY-3 DISTTHICKPLY4, 3, MAT4, DISTANGLEPLY4, PLY-4 *ORIENTATION, NAME=ORI DISTORIENT 3, DISTORIENTANGLE *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, E1, 0.1000000E+01, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.1000000E+01, 0.0000000E+00, E2, 0.9998477E+00, 0.1745241E-01, 0.0000000E+00, -0.1745241E-01, 0.9998477E+00, 0.0000000E+00, E3, 0.9993908E+00, 0.3489950E-01, 0.0000000E+00, -0.3489950E-01, 0.9993908E+00, 0.0000000E+00, E4, 0.9986295E+00, 0.5233596E-01, 0.0000000E+00, -0.5233596E-01, 0.9986295E+00, 0.0000000E+00, E5, 0.9975641E+00, 0.6975647E-01, 0.0000000E+00, -0.6975647E-01, 0.9975641E+00, 0.0000000E+00, E6, 0.9961947E+00, 0.8715574E-01, 0.0000000E+00, -0.8715574E-01, 0.9961947E+00, 0.0000000E+00, E7, 0.9945219E+00, 0.1045285E+00, 0.0000000E+00, -0.1045285E+00, 0.9945219E+00, 0.0000000E+00, E8, 0.9925462E+00, 0.1218693E+00, 0.0000000E+00, -0.1218693E+00, 0.9925462E+00, 0.0000000E+00, E9, 0.9902681E+00, 0.1391731E+00, 0.0000000E+00, -0.1391731E+00, 0.9902681E+00, 0.0000000E+00, *DISTRIBUTION, NAME=DISTORIENTANGLE, LOCATION=ELEMENT, TABLE=ORIANGLETABLE , 0., E1, 1., E2, 2., E3, 3., E4, 4., E5, 5., E6, 6., E7, 7., E8, 8., E9, 9., *DISTRIBUTION, NAME=DISTTHICKSECT, LOCATION=ELEMENT, TABLE=THICKNESSTABLE , 1. E1, 2.2110E-02 E2, 3.1053E-02 E3, 3.8610E-02 E4, 5.4483E-02 E5, 5.5110E-02 E6, 6.4053E-02 E7, 7.1610E-02 E8, 8.7483E-02 E9, 8.8110E-02 *DISTRIBUTION, NAME=DISTOFFSET, LOCATION=ELEMENT, TABLE=OFFSETTABLE , 0. E1, -0.5000000E+00, E2, -0.4750000E+00, E3, -0.4500000E+00, E4, -0.4250000E+00, E5, -0.4000000E+00, E6, -0.3750000E+00, E7, -0.3500000E+00, E8, -0.3250000E+00, E9, -0.3000000E+00, *DISTRIBUTION, NAME=DISTTHICKPLY1, LOCATION=ELEMENT, TABLE=THICKNESSTABLE , 0.00251 E2, 0.0051 E3, 0.00751 E4, 0.011 E5, 0.01251 E6, 0.0151 E7, 0.01751 E8, 0.021 E9, 0.02251 *DISTRIBUTION, NAME=DISTTHICKPLY2, LOCATION=ELEMENT, TABLE=THICKNESSTABLE , .02252 E1, 0.00252 E2, 0.0052 E3, 0.00752 E4, 0.012 E5, 0.01252 E6, 0.0152 E7, 0.01752 E8, 0.022 *DISTRIBUTION, NAME=DISTTHICKPLY4, LOCATION=ELEMENT, TABLE=THICKNESSTABLE , .02254 E1, 0.00254 E2, 0.0054 E3, 0.00754 E4, 0.014 E5, 0.01254 E6, 0.0154 E7, 0.01754 E8, 0.024 E9, 0.02254 *DISTRIBUTION, NAME=DISTANGLEPLY1, LOCATION=ELEMENT, TABLE=ORIANGLETABLE , 60. E2, 60.2 E3, 60.3 E4, 60.4 E5, 60.5 E6, 60.6 E7, 60.7 E8, 60.8 E9, 0. *DISTRIBUTION, NAME=DISTANGLEPLY3, LOCATION=ELEMENT, TABLE=ORIANGLETABLE , 0. *DISTRIBUTION, NAME=DISTANGLEPLY4, LOCATION=ELEMENT, TABLE=ORIANGLETABLE , 0. E1, 90. E2, 90.2 E3, 90.3 E4, 90.4 E5, 90.5 E6, 90.6 E7, 90.7 E8, 90.8 E9, 0. *DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE COORD3D, COORD3D *DISTRIBUTION TABLE, NAME=ORIANGLETABLE ANGLE *DISTRIBUTION TABLE, NAME=THICKNESSTABLE LENGTH *DISTRIBUTION TABLE, NAME=OFFSETTABLE RATIO *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, 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, 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, *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 3, TSHR, *NODE PRINT, FREQ=999 U, *OUTPUT, FIELD *NODE OUTPUT U, *ELEMENT OUTPUT S, SF, SSAVG, STH, *ELEMENT OUTPUT 3, TSHR, *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, 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