*HEADING                                                                                            
 MULTISTEP VERIFICATION TEST FOR S8R5 ELEMENTS USING DISTRIBUTIONS                                  
*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.6666667E+00,  0.0000000E+00,  0.0000000E+00
        6,  0.8333333E+00,  0.0000000E+00,  0.0000000E+00
        7,  0.1000000E+01,  0.0000000E+00,  0.0000000E+00
        8,  0.0000000E+00,  0.1666667E+00,  0.0000000E+00
        9,  0.1666667E+00,  0.1666667E+00,  0.0000000E+00
       10,  0.3333333E+00,  0.1666667E+00,  0.0000000E+00
       11,  0.5000000E+00,  0.1666667E+00,  0.0000000E+00
       12,  0.6666667E+00,  0.1666667E+00,  0.0000000E+00
       13,  0.8333333E+00,  0.1666667E+00,  0.0000000E+00
       14,  0.1000000E+01,  0.1666667E+00,  0.0000000E+00
       15,  0.0000000E+00,  0.3333333E+00,  0.0000000E+00
       16,  0.1666667E+00,  0.3333333E+00,  0.0000000E+00
       17,  0.3333333E+00,  0.3333333E+00,  0.0000000E+00
       18,  0.5000000E+00,  0.3333333E+00,  0.0000000E+00
       19,  0.6666667E+00,  0.3333333E+00,  0.0000000E+00
       20,  0.8333333E+00,  0.3333333E+00,  0.0000000E+00
       21,  0.1000000E+01,  0.3333333E+00,  0.0000000E+00
       22,  0.0000000E+00,  0.5000000E+00,  0.0000000E+00
       23,  0.1666667E+00,  0.5000000E+00,  0.0000000E+00
       24,  0.3333333E+00,  0.5000000E+00,  0.0000000E+00
       25,  0.5000000E+00,  0.5000000E+00,  0.0000000E+00
       26,  0.6666667E+00,  0.5000000E+00,  0.0000000E+00
       27,  0.8333333E+00,  0.5000000E+00,  0.0000000E+00
       28,  0.1000000E+01,  0.5000000E+00,  0.0000000E+00
       29,  0.0000000E+00,  0.6666667E+00,  0.0000000E+00
       30,  0.1666667E+00,  0.6666667E+00,  0.0000000E+00
       31,  0.3333333E+00,  0.6666667E+00,  0.0000000E+00
       32,  0.5000000E+00,  0.6666667E+00,  0.0000000E+00
       33,  0.6666667E+00,  0.6666667E+00,  0.0000000E+00
       34,  0.8333333E+00,  0.6666667E+00,  0.0000000E+00
       35,  0.1000000E+01,  0.6666667E+00,  0.0000000E+00
       36,  0.0000000E+00,  0.8333333E+00,  0.0000000E+00
       37,  0.1666667E+00,  0.8333333E+00,  0.0000000E+00
       38,  0.3333333E+00,  0.8333333E+00,  0.0000000E+00
       39,  0.5000000E+00,  0.8333333E+00,  0.0000000E+00
       40,  0.6666667E+00,  0.8333333E+00,  0.0000000E+00
       41,  0.8333333E+00,  0.8333333E+00,  0.0000000E+00
       42,  0.1000000E+01,  0.8333333E+00,  0.0000000E+00
       43,  0.0000000E+00,  0.1000000E+01,  0.0000000E+00
       44,  0.1666667E+00,  0.1000000E+01,  0.0000000E+00
       45,  0.3333333E+00,  0.1000000E+01,  0.0000000E+00
       46,  0.5000000E+00,  0.1000000E+01,  0.0000000E+00
       47,  0.6666667E+00,  0.1000000E+01,  0.0000000E+00
       48,  0.8333333E+00,  0.1000000E+01,  0.0000000E+00
       49,  0.1000000E+01,  0.1000000E+01,  0.0000000E+00
*NSET, NSET=SIDE1, GENERATE                                                                         
        1,       7,       1,
*NSET, NSET=SIDE3, GENERATE                                                                         
       43,      49,       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=S8R5, ELSET=EALL                                                                     
        1,       1,       3,      17,      15,       2,      10,      16,       8,
        2,       3,       5,      19,      17,       4,      12,      18,      10,
        3,       5,       7,      21,      19,       6,      14,      20,      12,
        4,      15,      17,      31,      29,      16,      24,      30,      22,
        5,      17,      19,      33,      31,      18,      26,      32,      24,
        6,      19,      21,      35,      33,      20,      28,      34,      26,
        7,      29,      31,      45,      43,      30,      38,      44,      36,
        8,      31,      33,      47,      45,      32,      40,      46,      38,
        9,      33,      35,      49,      47,      34,      42,      48,      40,
*DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE
 COORD3D, COORD3D
*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,
        1,   0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00,
        2,   0.9998477E+00,   0.1745241E-01,   0.0000000E+00,  -0.1745241E-01,   0.9998477E+00,   0.0000000E+00,
        3,   0.9993908E+00,   0.3489950E-01,   0.0000000E+00,  -0.3489950E-01,   0.9993908E+00,   0.0000000E+00,
        4,   0.9986295E+00,   0.5233596E-01,   0.0000000E+00,  -0.5233596E-01,   0.9986295E+00,   0.0000000E+00,
        5,   0.9975641E+00,   0.6975647E-01,   0.0000000E+00,  -0.6975647E-01,   0.9975641E+00,   0.0000000E+00,
        6,   0.9961947E+00,   0.8715574E-01,   0.0000000E+00,  -0.8715574E-01,   0.9961947E+00,   0.0000000E+00,
        7,   0.9945219E+00,   0.1045285E+00,   0.0000000E+00,  -0.1045285E+00,   0.9945219E+00,   0.0000000E+00,
        8,   0.9925462E+00,   0.1218693E+00,   0.0000000E+00,  -0.1218693E+00,   0.9925462E+00,   0.0000000E+00,
        9,   0.9902681E+00,   0.1391731E+00,   0.0000000E+00,  -0.1391731E+00,   0.9902681E+00,   0.0000000E+00,
*DISTRIBUTION TABLE, NAME=THICKNESSTABLE
 LENGTH
*DISTRIBUTION, NAME=DISTTHICK, LOCATION=ELEMENT, TABLE=THICKNESSTABLE                                        
         ,   0.1000000E-01,
        1,   0.1000000E-01,
        2,   0.2000000E-01,
        3,   0.3000000E-01,
        4,   0.4000000E-01,
        5,   0.5000000E-01,
        6,   0.6000000E-01,
        7,   0.7000000E-01,
        8,   0.8000000E-01,
        9,   0.9000000E-01,
*DISTRIBUTION TABLE, NAME=OFFSETTABLE
 RATIO
*DISTRIBUTION, NAME=DISTOFFSET, LOCATION=ELEMENT, TABLE=OFFSETTABLE                                       
         ,  -0.5000000E+00,
        1,  -0.5000000E+00,
        2,  -0.4750000E+00,
        3,  -0.4500000E+00,
        4,  -0.4250000E+00,
        5,  -0.4000000E+00,
        6,  -0.3750000E+00,
        7,  -0.3500000E+00,
        8,  -0.3250000E+00,
        9,  -0.3000000E+00,                                                      
*DISTRIBUTION TABLE, NAME=ORIANGLETABLE
 ANGLE
*DISTRIBUTION, NAME=DISTORIANGLE, LOCATION=ELEMENT, TABLE=ORIANGLETABLE
         ,   0.0,                                     
        1,  -0.5
        2,  -0.5
        3,  -1.0
        4,  -1.0
        5,  -1.5
        6,  -1.5,
        7,  -2.0,
        8,  -2.0,
        9,  -2.5,
*ORIENTATION, NAME=ORI
 DISTORIENT
 3, DISTORIANGLE
*SHELL SECTION, COMPOSITE, ELSET=EALL, ORIENTATION=ORI, OFFSET=DISTOFFSET, SHELL THICKNESS=DISTTHICK                                                           
 DISTTHICKPLY1, 3, MAT1,           -45, PLY-1
            1., 3, MAT1, DISTANGLEPLY2, PLY-2                 
*DISTRIBUTION, NAME=DISTTHICKPLY1, LOCATION=ELEMENT, TABLE=THICKNESSTABLE                                        
     , 0.1000000E-01
 EALL, 1.0
*DISTRIBUTION, NAME=DISTANGLEPLY2, LOCATION=ELEMENT, TABLE=ORIANGLETABLE
     , 0.1000000E-01
 EALL, 45
*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,  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                                                                                 
 3,                                                                                             
 TSHR,                                                                                              
*NODE PRINT, FREQ=999                                                                               
 U,                                                                                                 
*OUTPUT, FIELD                                                                                      
*NODE OUTPUT                                                                                        
 U,                                                                                                 
*ELEMENT OUTPUT                                                                                     
 S,                                                                                                 
 SF,
 SSAVG,
 SE,
 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,
 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