*Heading                                                                                            
 Multistep verification test for cpe4r elements using distributions                                 
*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
       13,      14,       15, 16
*Nset, nset=ams_out, generate
        9,      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=cpe4r, 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,
*Solid Section, material=mat1, elset=eall, orientation=ori                                                         
1.,                                                         
*Orientation,name=ori
 distOri
 3,0.                                  
*Distribution Table, name=distTableOri
 coord3d,coord3d      
*Distribution, name=distOri, table=distTableOri, location=element 
        ,  1., 0., 0.,  0., 1., 0.                             
        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,
*Boundary                                                                                           
 side1, 1, 2                                                                                        
**
** --------------------------------------------------
**
*step, name=freq
 frequency analysis
*frequency, eigensolver=ams, residual modes, nset=ams_out
 , ,14000.
*output, field
*node output, nset=ams_out
U,
*end step
**
** --------------------------------------------------
**
*step, name=step-3
 SSD mode-based
*steady state dynamics, interval=eigenfrequency
 13000,
*cload
 side3, 2, 100.
** 16,    2, 100.
*output, field
*node output, nset=ams_out
U,
*end step
**
** --------------------------------------------------
**