*HEADING
: CIRCULAR FRAME SECTION, 
                         WIND LOADING (*TRANSFORM)  [FRAME2D]
*NODE
 1,  0.0, 60.0
11, 64.95, 97.5
*NGEN,NSET=ALL
 1,11,2
*AQUA
 -100., 50., 10.0, 8000.0E-6
 15., -25.98, 0., -100.
 15., -25.98, 0.,  50.
*WIND
8.0E-3, 10., 11.209982, 0.000001, 1.0, 0., 0.000000001
*TRANSFORM,TYPE=R,NSET=ALL
  64.95, 37.5, 0.0, -37.5, 64.95, 0.0
*ELEMENT,TYPE=FRAME2D, ELSET=EALL
1,1,3
2,3,5
3,5,7
4,7,9
5,9,11
*ELSET, ELSET=FOUR
4, 
*FRAME SECTION,SECTION=CIRC, ELSET=EALL, DENSITY=8000.0E-4
2., 
0., 0., -1.
3.E6, 1.1538E6, .001
*BOUNDARY
1, 1, 2
1, 6
*STEP,PERTURBATION
*STATIC
*CLOAD
11,2,25.E-4
*EL PRINT
SF,
SEE,
*NODE PRINT, GLOBAL=YES
 U,
 RF, 
*NODE PRINT, GLOBAL=NO
 U,
 RF, 
*OUTPUT, FIELD, FREQ=99, VAR=ALL 
*OUTPUT, HISTORY, FREQ=99, VAR=PRESELECT 
*ENERGY OUTPUT, VAR=PRESELECT 
*END STEP
*********************** STEP 2
*STEP
*STATIC
*DLOAD
EALL, GRAV,1.0E-3, 37.5, -64.95, 0.
*EL PRINT
SF,
SEE,
*NODE PRINT, GLOBAL=YES
 U,
 RF, 
*NODE PRINT, GLOBAL=NO
 U,
 RF, 
*OUTPUT, FIELD, FREQ=99, VAR=PRESELECT
*OUTPUT, HISTORY, FREQ=99, VAR=ALL
*END STEP
*********************** STEP 3
*STEP
same as 2nd step 
*STATIC
*DLOAD, OP=NEW 
EALL, PX,   5.026548E-3 
EALL, PY,  -8.705982E-3 
*END STEP
*********************** STEP 4
*STEP,PERTURBATION
*STATIC
*DLOAD
EALL, WDD,  1.0E-2, 4., 2., , 
*EL PRINT
SF,
SEE,
*NODE PRINT, GLOBAL=YES
 U,
 RF, 
*NODE PRINT, GLOBAL=NO
 U,
 RF, 
*OUTPUT, FIELD, FREQ=99, OP=ADD
*ELEMENT OUTPUT, ELSET=EALL
SF, SEE, THE
*NODE OUTPUT, NSET=ALL
U, RF, 
*OUTPUT, HISTORY, FREQ=99, OP=REPLACE
*ELEMENT OUTPUT, ELSET=EALL
SF, SEE, 
*NODE OUTPUT, NSET=ALL
U, RF, 
*ENERGY OUTPUT, VAR=ALL  
*END STEP
*********************** STEP 5
*STEP 
 axial tensile loading: 
 The RF at the support node should match the 
 RF magnitude in 4th step at the support.
*STATIC
*DLOAD, OP=NEW 
EALL, WD1, 1.0E-2, 20., 2., ,
*END STEP