$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$       Sample Problem Comaparing a Longitudinal Strength Analysis
$                      to a Structual Analysis
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$*********************************************      define units
$
&dimen -save -dimen meters knts
$
$*********************************************      define body
$
&describe body barge -section 113e11 0 75.333 
pgen barge
   plane 0   5 10   15  20  25  30  35  40  45 \
        50  55 60   65  70  75  80  85  90  95 \
       100 105 110 115 120 125 130 135 140 145 \
       150 155 160 165 170 175 180 185 190 195 \ 
       200 205 210 215 220 225 230 235 240 245 \
       250 255 260 265 270 275 280 285 290 295 \ 
       300 305 310 315 320 325 330 335 340 345 \
       350 355 360 365 370 375 380 385 390 395 \ 
       400   \                       
                -rect 0 15 100
end pgen
$
$*********************************************      generate beams
$
&describe body barge 
&default -nuse @
$
$ define element class
$
$  Conversion to units for -section option 
$ 1 MPa = nt/mm^2
$ EI = nt/mm^2 * mm^4 = nt-mm^2
$ EI = nt-mm^2* m^2/1000^2 mm^2* k/1000 = knts-m^2
$
$ SM = I/c = 
$
~bardum pri (21000*2) 15000 -section 4.398e8 5.65e16 7.854e16 9.6304e16 
$
&set list      =          $ initialize variable to collect node names
&set num_beams = 80       $ number of beams
&set dist = 5             $ length of beams
&loop lll 1 &number(real %(num_beams)+1) 1
    &set kkk = &number(integer -1+%lll)      $ beam number
    &set xl = &number(integer %dist*%kkk)    $ longitudinal location of end of beam
$
    &set name = &string(o_number n00 %(kkk)) $ setup name of beam and node
    *%(name)   %xl   0   0                   $ define end node
    &set list = %list *%(name)               $ collect node names
    &if %lll .gt. 1 &then
         beam %(name)  ~bardum  %(last) *%(name)   $ define beams
    &endif
    &set last =  *%(name)
$
    &if %lll .gt. 11 .and. %lll .lt. 72 &then
       #elat %(name) 16070 0 0 0 0   -cat elat  $ define element weight
    &endif
$
&endloop
&dimen -remember