$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$         ----------    JACKET PILE DOCKING COMMANDS ----------
$
$    This set of commands will perform an analysis of a jacket docking
$ with a pile attached to the seafloor.  The pile is simulated as a
$ pair of truselem connectors, in the lateral X and Y directions.  P-Y
$ data is input to provide the stiffness of the pile.  Two time domain
$ analyses are performed.  The first one determines the maximum lateral
$ velocity of the impact point on the jacket, as a result of the
$ applied environment.  The ADD_COL   command in the Disposition Menu
$ is used to find the normal of X and Y  velocities.
$
$    The time domain analysis is then restarted using  the -RESET
$ option, so the connector simulating the pile is recognized as active.
$ Maximum force in this connector is the impact load on the pile.
$ For the sake of brevity,  the time of simulation  is ridicuosly
$ short.  This sample is relatively complete, including a derrick
$ barge with spread mooring,  and sling assemblies for the jacket.
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
&device -g_default file -oecho n
&dimen -dimen feet kips
$
$*********************************************      read model
$
&title 8 Pile Jacket * Pile Docking * TDOM -RESET METHOD *
inmodel
&dcptime Time to read model
$
$*********************************************      Set Param
$
&set time    = 20  1.0
$&set time    = 900 0.5
&set imptime = 8
&set imptste = .1
$
$*********************************************      Set Wgt & Buoy
$
&instate -locate jacket 0 0 -2000
&set jktbuoy = &body(displac jacket)
&set bumar   = &number(real 12377/%jktbuoy)
&instate -locate jacket 0 0 0
&weight -total  jacket 10913 -1.573  -.878   -235.728 135.679      \
             143.317   86.235
&apply -fraction -category str_model #buoy %bumar%
&summary
   &rep_sel -part jacket
   categ
   &rep_sel -part barge
   categ
   &rep_sel -selall
end
$
$*********************************************      Set Initial Cond
$
&instate -locate barge 0 0 -19 0 0 0 -locate jacket 570 0 40 0 0 180
&weight -compute barge 15 41 131 131
$
$*********************************************      Define Connections
$
MEDIT
   &describe body barge
   *bsa   450   64   25
   *bpa   450  -64   25
   *bpf   0    -64   25
   *bsf   0     64   25
   *btip  570   0    30-19+300
   *bpc   425  -15   30
   *bsc   425   15   30
   ~crane tube 20 1
   beam ~crane -nuse @ -use #buoy *btip *bpc
   beam ~crane -nuse @ -use #buoy *btip *bsc
   ~WIRE b_cat 3.0  -depanc 386 -len 4000 -buoydia 0 -wtplen .0125
   CONNECTOR bpa_1  ~WIRE   *bpa -ANC  -45   20
   CONNECTOR bpa_2  ~WIRE   *bpa -ANC  -90   20
   CONNECTOR bpf_1  ~WIRE   *bpf -ANC  -90   20
   CONNECTOR bpf_2  ~WIRE   *bpf -ANC -135   20
   CONNECTOR bsf_1  ~WIRE   *bsf -ANC  135   20
   CONNECTOR bsf_2  ~WIRE   *bsf -ANC   90   20
   CONNECTOR bsa_1  ~WIRE   *bsa -ANC   90   20
   CONNECTOR bsa_2  ~WIRE   *bsa -ANC   45   20
   ~hoist  sling 20 -length 400
   ~sling  sling 4  -length 100
   connector boom ~hoist *btip
   connector sl1381 ~sling  *j1381
   connector sl1401 ~sling  *j1401
   connector sl1389 ~sling  *j1389
   connector sl1409 ~sling  *j1409
   t-h boom sl1381 sl1401 sl1389 sl1409  -initial
$
$*********************************************      Define Pile
$
   &set nodecor = &node(coord *j152 -global)
   &set nodeX   = &token(1 %nodecor%)
   &set nodeY   = &token(2 %nodecor%)
   &set nodeZ   = &token(3 %nodecor%)
   &describe -part ground
   *pile_X -3000+%nodeX%  %nodeY%      %nodeZ%
   *pile_Y       %nodeX%  3000+%nodeY% %nodeZ%
   &set len_X = &node(d_node *j152 *pile_X)
   &set len_Y = &node(d_node *j152 *pile_Y)
   &set py =       0 0,        5 2/12                              \
      10 4/12,    15 7/12,    20 10/12,    25 12.5/12,  30 15/12   \
      35 18.5/12, 40 22/12,   45 26/12,    50 30/12,    55 34/12   \
      60 38/12,   65 42.5/12, 70 48/12,    75 52/12,    80 58/12   \
      85 63/12,   90 69/12
   ~pile_X sl_elem 1.1284*12 -len %len_X% -emodulus 1e10 -x_py %py%
   ~pile_Y sl_elem 1.1284*12 -len %len_Y% -emodulus 1e10 -x_py %py%
   connector pile_X ~pile_X *j152 *pile_x
   connector pile_Y ~pile_Y *j152 *pile_Y
END
$
$*********************************************      Move Anchors
$
&connector pile@    -inactive  bp@  bs@ -A_TENSION 210
&dcptime Time to Connect
$
$*********************************************      Find Equi
$
&compartment -percent legB@ 70 legc@ 100
&weight -compute barge 15 41 131 131
$
$*********************************************      Plots
$
&picture iso
&picture iso -x -100 700 -y -200 200
&picture side
&picture top
$
$*********************************************      PRCOMP
$
hydrodynamics
   g_press barge
end_hydrodynamics
$
$*********************************************      Add Environ
$
&data environment
   c_profile nemba 0 0.52*1.689 386 0.26*1.689
   environment swell1  -sea issc 180-18 4.10 10.0/1.2958  -time %time% \
                       -wind 10 180+10 -current nemba 180+10           \
                       -depth 386 -use_mean y
   environment swell2  -sea issc 180-45 4.10 10.0/1.2958  -time %time% \
                       -wind 10 180+10 -current nemba 180+10           \
                       -depth 386 -use_mean y
   environment swell3  -sea issc 180    4.10 10.0/1.2958  -time %time% \
                       -wind 10 180+10 -current nemba 180+10           \
                       -depth 386 -use_mean y
end_&data
&env swell1
$
$*********************************************      Equi, Set Vari
$
&equi -omega .2
&status tip-hook
&status force -hard
&status f_connect -hard
&set jktloc  = &body(location jacket -global)
&set brgloc  = &body(location barge  -global)
&set jktxor  = &token(1 %jktloc%)
&set jktyor  = &token(2 %jktloc%)
&set jktzor  = &token(3 %jktloc%)
&set jktrzor = &token(6 %jktloc%)
&set brgzor  = &token(3 %brgloc%)
$
$*********************************************      Time Domain
$
&subtitle Before Impact
tdom -newmark
&describe interest -associate *j152
$
$*********************************************      Post Process
$
prcpost
   &if &uniname(bl) .eq. 'feet' &then
      &set l = Ft
      &set v = Ft/Sec
   &else
      &set l = M
      &set v = M/Sec
   &endif
   &set f = &uniname(bf)
$
$*********************************************      Select Max
$
   points *j152
      vlist
      add_column V-XY -norm 10 2
      vlist
      plot 1 4 -rax 5 -no_edit
      plot 1 10 11 35 -no_edit                 \
      -t_x 'Time, Sec'                  \
      -t_left  'Lateral Velocity of Jkt Point, %v%'
      plot 1 12     -no_edit                  \
      -t_x 'Time, Sec'                  \
      -t_left  'Vertical Velocity of Jkt Point, %v%'
      statistic 1   4   5  6 27   -hard
      statistic 1  10  11 12 34 -hard
      extreme   1  10  11 12 34 -hard
      set_variable maxtime  -max 1  35
   end
   cf_magnitude
      report
   end
   trajectory
      report location velocity
   end
end
$
$*********************************************      Impact
$
&subtitle After Impact
&instate -event %maxtime%
&env swell1 -sea issc 180-45 4.10 10.0/1.2958     \
            -time %maxtime%+%imptime% %imptste   \
            -wind 10 180+10 -current nemba 180+10 \
            -depth 386 -use_mean y
&connector pile@ -active  -l_tension 0
&status g_connector
&status f_connector
&type len_X = &node(d_node *j152 *pile_X)
&type len_Y = &node(d_node *j152 *pile_Y)
tdom -newmark -reset  %maxtime%
$
$*********************************************      Post Process
$
prcpost
   cf_magnitude pile@ -event %maxtime% %maxtime%+%imptime%
      add_column pile_xy -rms 2 4
      vlist
      plot 1 2 4 6    -no               \
      -t_x 'Time, Sec'                  \
      -t_left  'Force in Pile, &uniname(bf)'
      extreme   1 2 4 6 -hard
      statistic 1 2 4 6 -hard
      report
   end
end
$
$*********************************************      All Done
$
&finish
