$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$  This models a steel buoy with a Rod Element connecting it to a
$  subsea manifold on the seafloor.
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$**************************************          program settings
$
&title Buoy on Rod Sample
&device -oecho yes -g_default device
&dimen -dimen feet kips
&default -depth 6000
&pb_parm -R_TOLLE  2e-5  2e-5
$
$**************************************          set up
$
inmodel
&set draft  = 130
&set bottom = 6000-20
&describe body buoy -ignore rz
&instate -condition buoy %draft% 0.0 0.0
$
$**************************************          define connectors
$
medit
   &set props = -wtpl 1.789 -displ 1.769 -sect 147.7 41865 41865 2*41865 \
                -buoy 0
   &set l1 = 400
   &set l3 = 200
   &set l2 = &number(real %bottom%-%draft%-%l1%-%l3%)
   ~tubebp rod 60 1 -refine 20 -len %l1 %props -depan %bottom%
   ~tubebp rod 60 1 -refine 50 -len %l2 %props
   ~tubebp rod 60 1 -refine 10 -len %l3 %props
   connector tube ~tubebp -anchor 0.0 0.0 *pt1
end_medit
&connector tube   -st_addition 2 1572.5 1572.5 0 \
                  -a_stiff 600 60 60  \
                  -top_moment yes
&equi -omega .25
$
$**************************************        Define  environment
$
&data env
    profile loopcur     0.0  6.10 \
                      300.0  6.10 \
                      500.0  5.30 \
                     1000.0  2.00 \
                     2000.0  0.50 \
                     6000.0  0.50
   env gomloop -sea issc 0 39.0 15.2/1.289 -current loopcur 0 \
               -time 100 .25
end_&data
&env gomloop
$
$**************************************         find equi
$
&status force
&instate -locate buoy  -173.09 .00   -135   0 0 0
&equi -omega .25
&status force
&status g_rod
&status s_rod
&status f_rod
$
$**************************************          frequency domain
$
freq_res
   rao  -heading 0 -period 7 8 9 10 11 12 13 14 15 16 17 18 19 20 22.4 \
                           25
   fr_point 0.0 0.0 50.0
   report
   end
   fr_cforce
   report
   end
   st_point gomloop
   report
   end
   st_cforce tube gomloop
   report
   end
   st_rforce tube gomloop
   report
   end
end
$
$**************************************          time domain
$
tdom -newmark
prcpost
     trajectory  -events 0 900 .5
       report location
        end
     r_envelope tube
        report
        end
     r_view  tube -events 0 900 .5
        report
        end
end
$
$**************************************         all done
$
&fini
$