$
$*********************************************      basics
$
&dimen -save -dimen meters m-tons
&device -oecho n -g_defaultmary device
&title fixtube
$
$*********************************************      read model
$
inmodel
&instate -locate fixtube 0 0 0
medit
  ~fix fix
  connector ca ~fix *a
  connector cb ~fix *b
  connector cc ~fix *c
end_medit
&parameter -slam_both yes
&set dir = 0
&set height = 19.5
&set period = 12.0
&env -ramp 0 -sea regular %dir %height %period -time %period%*10 0.5
tdom
prcpost
  ldg tub *mid -force @  -even 0 900 .1
  report -even 0 9000 .5
  set_variable maxtotx -maximum 93 93
  set_variable maxtotz -maximum 95 95
  vlist
  plot 1 16 30 37 44 79 93 \
    -t_main 'Forces in global X-direction - H = %height m, T = %period sec, Dir = %dir deg' \
    -t_sub 'Max total = %maxtotx m-tons (&number(real %maxtotx*9.81) kN)' \
    -t_left 'Force [m-tons]' \
    -t_x 'Time [sec]' \
    -no_edit
  plot 1 18 32 39 46 81 95 \
    -t_main 'Forces in global Z-direction - H = %height m, T = %period sec, Dir = %dir deg' \
    -t_sub 'Max total = %maxtotz m-tons (&number(real %maxtotz*9.81) kN)' \
    -t_left '\Force [m-tons]' \
    -t_x 'Time [sec]' \
    -no_edit
  end $plot
end $prcpost
$
$*********************************************      all done
$
&fini
