$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$     ----------     Flag Pole With Generalized DOF    ---------
$                                                file 3
$
$ This is the second in a series of samples which deals with a "flag pole".
$ The pole is almost a cantilever beam. It is a set of 10 beam elements
$ which are pinned at the first and second nodes.
$
$ This sample, also treats the pole in water. It again extracts 2 modes
$ using the same boundary conditions as those used later. The difference
$ between this solution and that obtained in the second samples is that
$ here, the modes are computed without a weight at the top. We compute
$ RAOs of this system, we then apply the weight at the top, and compute
$ them again. The thing here, is that the RAOs for the case with the
$ weight are "sensibly the same as those found in sample 2, even though
$ we are using modes which were created without the top weight.
$
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$*********************************************      basic parameters
$
&dimen -dimen meters k-nts
&device  -g_default device -oecho no -auxin f_pole.dat
$&device  -g_default screen
&set num_mode = 10
&title Flag Pole - Modes Computed in Water - No Weight
$
$*********************************************      read model
$
inmodel
&apply -perc @ 0 f_pole 100
$
$*********************************************      estimates
$
e_estimate
$
$*********************************************      set configuration
$
&instat -locat 0 0 00
$
$*********************************************      define connections
$
medit
   connector c1 ~spr1 *1
   connector c2 ~spr2 *2
end
&equi
&stat force
&event_store 1
$
$*********************************************      get modes
$
struct -init
   modes -num_ev %num_modes% -no_fix
end
strpos
   modes value
   &loop i 1 %num_mode%
      &set cname = M0@0%i%
      modes vector -load %cname%
   &endloop
end
$
$*********************************************      use modes
$
&describe body f_pole -gen 1:%num_mode%
&equi
&stat motion
$
$*********************************************      raos
$
freq_resp
    mor_work 0
    rao -heading 0 -period %liste%
    &subt Motion of Node *%n_nodes% - No Weight
    fr_poin *%n_nodes%
       report
       add_col Energy -inpu %l_drao
       plot 2 3 15 -t_left 'X Motion RAO (m)' \
                   -t_x    'Period (Seconds)'  \
                   -legend 1 'Moses '         \
                   -legend 2 'Energy Estimate'
       end
end
$
$*********************************************      raos
$
&apply -percent @ 100 -force -percent @ 100
&equi
&stat motion
freq_resp
    mor_work %twei -nonl
    rao -heading 0 -period %liste%
    &subt Motion of Node *%n_nodes% - With Weight
    fr_poin *%n_nodes%
       report
       add_col Energy -inpu %l_drao
       plot 2 3 15 -t_left 'X Motion RAO (m)' \
                   -t_x    'Period (Seconds)'  \
                   -legend 1 'Moses '         \
                   -legend 2 'Energy Estimate'
       end
end
$
$*********************************************      all done
$
&fini