$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$     ----------     Flag Pole With Generalized DOF    ---------
$                                                file 1
$
$ This is the first 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, the first sample, treats the pole in air. It 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 fond equilibrium of the system subject to a horizontal force and
$ compare this with that obtained from a standard structural solution.
$
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$*********************************************      basics
$
&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 Air - No Weight
$
$*********************************************      read model
$
inmodel
$
$*********************************************      set configuration
$
&instate -locat 0 0 100
&apply -percent @ 0 f_pole 100 top_for 100 -force -percent @ 100
$
$*********************************************      define connections
$
medit
   connector c1 ~spr1 *1
   connector c2 ~spr2 *2
end
&equi
&status force
&event_store 1
$
$*********************************************      find deflections
$
struct
   lcase -process 1
   ssolve
end
strpost
    joint disp
    beam  load -detail
end
$
$*********************************************      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
$
$*********************************************      all done
$
&fini
