$ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $ $ $ ---------- 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%)-1) &instate -locate jacket 0 0 0 &weight -total jacket 10913 -1.573 -.878 -235.728 135.679 \ 143.317 86.235 -margin jacket 0 %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 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 swell2 $ $********************************************* Equi, Set Vari $ &equi -tolerance 3e-4 5e-4 -iter_max 150 &status tip-hook &status force -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 8 9 35 -no_edit \ -t_x 'Time, Sec' \ -t_left 'Lateral Velocity of Jkt Point, %v%' plot 1 10 -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 35 -hard extreme 1 10 11 12 35 -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 swell2 -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 &stat g_c &stat f_c &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