$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$     ----------     Flag Pole With Generalized DOF    ---------
$                                                file 2
$
$ 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, treats the pole in water. It again extracts 2 modes
$ using the same boundary conditions as those used later. After the
$ modes are obatined, we use them as generalized degrees of freedom.
$ We find equilibrium of the system, and then compute RAOs of the system.
$ The thing of interest here is that the peak in the RAO curve corresponds
$ directly to the first natural frequency.
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$*********************************************      basic parameters
$
&dimen -dimen meters k-nts
&device  -g_default device -oecho no -auxin f_pole.dat
&set num_mode = 2
&title Flag Pole - Modes Computed in Water - With Weight
$
$*********************************************      read model
$
inmodel
$
$*********************************************      estimates
$
e_estimate
$
$*********************************************      set configuration
$
&instate -locat 0 0 00
&apply -percent @ 100 -force -percent @ 100
$
$*********************************************      define connections
$
medit
   connector c1 ~spr1 *1
   connector c2 ~spr2 *2
end
&equi
&status 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
&status motion
$
$*********************************************      raos
$
freq_resp
    rao -heading 0 -period %liste%
    &subt Motion of Node *%n_nodes%
    fr_poin *%n_nodes%
       report
       end
end
$
$*********************************************      all done
$
&fini