$
$****************************************       macros
$
&macro lc_doit name case num = 1
   &env %name
   &apply -force -fract @ 0 %case  %num
   lcase -static %name
&endmacro
$
$****************************************       dimensions
$
&dimen -dim feet kips
$
$****************************************       describe body
$
&describe body model
&default -emo 29000 -nuse @
$
*j0102 50 50 0
*j0202 41+8/12 43+4/12 50
&describe body ground
*bop -145 -ref *j0102 *j0202
&describe body model
$
$
~bem tube 45 1.25
~pil tube 45 1.25 -soil sol1 -refi 50
$
beam ~bem *j0102 *j0202
beam ~pil *j0102 *bop
$
&data soil sol1
   mpy 45/30
   mtz 3.1415*45*.0139
   mqw 1440*.139
   depth 0
   py 0,0 .002,.5 .002,11
   tz 0,0 .25,.012 .5,.03 .75,.048 1,.12 1,.24 .95,1 .9,2 .9,60
   qw 0,-60 0,0 .25,.012 .5,.03 .75,.048 1,.12 1,.24 .95,1 .9,2 .9,60
   depth 5.8
   py 0,0 .033,.5 .033,1.5 .034,2 .034,50
   depth 6.2
   py 0,0 .255,.5 .384,1 .461,1.5 .513,2 .578,3 .617,4 .662,6 .687,8 +
      .709,11 .725,15 .737,20 .759,50
   depth 10
   py 0,0 .797,.5 1.293,1 1.631,1.5 1.876,2 2.208,3 2.423,4 2.683,6 +
      2.835,8 2.974,11 3.080,15 3.158,20 3.309,50
   depth 16.5
   py 0,0 .878,.5 1.52,1 2.009,1.5 2.394,2 2.962,3 3.361,4 3.884,6 4.212,8 +
      4.524,11 4.776,15 4.966,20 5.349,50
   depth 301
   py 0,0 .878,.5 1.52,1 2.009,1.5 2.394,2 2.962,3 3.361,4 3.884,6 4.212,8 +
      4.524,11 4.776,15 4.966,20 5.349,50
   tz 0,0 .25,.012 .5,.03 .75,.048 1,.12 1,.24 .95,1 .9,2 .9,60
   qw 0,-60 0,0 .25,.012 .5,.03 .75,.048 1,.12 1,.24 .95,1 .9,2 .9,60
end
$
$ #load *j0202 *j0102 -1400 -local -len 0  -b 0
#load *j0202 *j0102 426.927 290.999 -1236.561 +
                    11075.672/12 -19449.881/12 -726.637/12 -len 0 -b 0