$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$       Sample Problem Comaparing a Longitudinal Strength Analysis
$                      to a Structual Analysis
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$*********************************************      basic settings
$
&device -g_default file -oecho no
$
$*********************************************      dimensions
$
&dimen -dimen meter knts
$
$*********************************************      read model
$
inmodel
$
$*********************************************      set condition
$
&instate -condition barge 12 0 0
$
$
$  fix ends
$
medit
   ~spr1 spr x 1e4 y 1e4 z 1e4
   ~spr2 spr x 1e4 y 1e4 z 1e4 rx 1e8
   rest ~spr1 *n00
   rest ~spr2 *n80
end
$
$*********************************************      plot the model
$
&picture iso
$
$*********************************************      compute weight
$
&status config
&status force
&weight -compute barge 12.5 32 .29*400 .29*400
&equi
&status force
&status b_w
&status cat
$
$*********************************************      summary
$
&summary
   class_sum
   compart_sum long
end
$
$*********************************************      Traditional long. Str.
$
hstatics barge
   moment
      plot 1 2 -rax 3 -t_sub " Traditional Results"
      report
      set shr  -column 1 2  $ capture location and shears into a variable
      set bmom -column 1 3  $ capture location and moment into a variable
      set defl -column 1 5  $ capture location and deflection into a variable
      store 1 2 3 5         $ store values in a csv format file
   end
end
$
$*********************************************      structural analysis
$
&env none
struct
   lcase -static none
   ssolve
end
$
$*********************************************      change sign of traditional results to match
$      coordinate system of structural analysis
$
&set n_shr =                        $ initialize variable for transformed values
&loop ppp 1 &token(n %(shr)) 2
   &set chg = &token(1+%ppp %shr)
   &set chg = &number(real -1*%chg) $ change sign
   &set xval = &token(%ppp %shr)
   &set n_shr = %n_shr %xval %chg   $ save location & transformed value
&endloop
$
&set n_bmom =                         $ initialize variable for transformed values
&loop ppp 1 &token(n %(bmom)) 2
   &set chg = &token(1+%ppp %bmom)
   &set chg = &number(real -1*%chg)   $ change sign
   &set xval = &token(%ppp %bmom)
   &set n_bmom = %n_bmom %xval %chg   $ save location & transformed value
&endloop
$
strpost
   joint disp -local yes -file yes    $ find, report, store joint displacements
   rest -detail
   beam load -file yes
$
$ REMEMBER: the load is distributed but the buoyancy is lumped at the nodes
$           also the variable 'list' was defined in the dat file
$
   bmom_shr %(list)
      add_column  trad_sh -input 1 %n_shr
      add_column  trad_mb -input 1 %n_bmom
      add_column  trad_ds -input 1 %defl
      vlist
      store 1 4 8 6 9                      $ store all variables in a separate file
      plot 1 4 8 -t_sub "Shear Comparison" -t_x "Long. Location" \
                 -t_left "Shear"                                 \
                 -legend 1 Structural -legend 2 Traditional
      plot 1 6 9 -t_sub "Moment Comparison" -t_x "Long. Location" \
                 -t_left "Moment"                                 \
                 -legend 1 Structural -legend 2 Traditional
   end
end
&finish
