
 p r o g r a m  C A I S S O N _ V H M  w32 version 00.12_b04 
 ------------------------------------  at 15:35:38 on 22-AUG-2024
 caisson foundation resistance in undrained soil during VHM load.
 (c) S.Kay Consultant, Leiderdorp, The Netherlands.
 Please read the CAISSON_VHM Users Manual [ref SKA/MAN/003]
         and the CAISSON_VHM FAQ File     [ref SKA/FAQ/003].


 datain_caisson: < --- BEGIN

 pidata: < --- BEGIN
 pidata: input  text file: "JOB.I01"
       : output text file: "JOB.O01"
       : output plot file: "D4.5_PEN15_ZLUG9_SLIDING.O1B.PDF"
       : program name    : "CAISSON_VHM"

 pidata: rtitle          : "RUN 1B: AS RUN 1 BUT NLC=10 (WAS 1) AND T_LUG = 0 THRU 90DEG"
       : ptitle          : "EXAMPLE MODE22"
       : projno          : "SKA/MAN/003"
       : locatn          : "CAISSON_VHM"

 pidata: END --- >


 gmdata: < --- BEGIN
 numsl
    4
  zdep        bden       su_top      su_bot   
   [m]     [kN/m^3]      [kPa]        [kPa]   
    0.00        4.00        1.00        4.00    
   0.700        4.00        4.00        6.50    
    4.60        5.00        6.50        19.0    
    18.0        5.00        19.0        41.0    
    40.0    
 m_fac
    1.00    
 gmdata: END --- >


 psdata: < --- BEGIN
 psdata: vhm_model
 KAY+PALIX_2011
 MODE_VHM
   22
 psdata: END --- >


 cfdata: < --- BEGIN
 ncgeom
    2
 cgeom
 CIRCULAR
 nfhxp  wt_D_side  wt_D_top
  0   0.780E-02   0.111E-01
 iro_lug  r_lug_user
  1    0.00    
 D  L  zlug
   4.50  15.00  10.00
   4.50  16.00  10.50
 cfdata: END --- >


 ssidat: < --- BEGIN
 alpha_Fo  alpha_Fi  eta_Fo
   0.51   0.42   1.00
 eta_reb  alpha_d_su  iplug_vmode  vmax_user
   1.00   0.50    2   0.00
 eta_amh  eta_bmh  eta_phimh
   1.00   1.00   1.00
 ssidat: END --- >


 chndat: < --- BEGIN
 iscope_Bw  l_fac_Bw  eta_Vst_wall  delta_Bw_user  Bw_caisson_user
    1    1.00        1.00        0.00        0.00    
 dtheta_lug_deg  psi_twist_deg
    0.00        0.00    
 chndat: END --- >


 ldata__anchor: < --- BEGIN
 nlc   nlfac_T
 10  1
 ldata__anchor: theta_lug  (l_facT) ***
    0.00        1.00       ! ? missing text
    10.0        1.00       ! ? missing text
    20.0        1.00       ! ? missing text
    30.0        1.00       ! ? missing text
    40.0        1.00       ! ? missing text
    50.0        1.00       ! ? missing text
    60.0        1.00       ! ? missing text
    70.0        1.00       ! ? missing text
    80.0        1.00       ! ? missing text
    90.0        1.00       ! ? missing text
 ldata__anchor: END --- >


 mcdata: < --- BEGIN
 ioutd  ioute
    3   -2
 eps_int  jmax_int  jmin_int
   0.100E-03 15  5
 mcdata: END --- >

 datain_caisson: END --- >


clen_theoretical_sub: < --- BEGIN
clen_theoretical_sub: no theoretical caisson L values for FACTORED VMH load;
clen_theoretical_sub: can't find L without VHload (MODE_VHM = 16, 17, 22 and 23).
clen_theoretical_sub: END --- >



caisson_vhm: nlines_warnings :    7
caisson_vhm: nlines_errors   :    0

WARNING: gmdata. soil layer # 4; BDEN=5.00 kN/m^3 is inconsistent with SU_BOT=41.00 kPa.
WARNING. cfdata: NCGEOM > 1 with offset lug implies variable ZLUG*. [?IRO_LUG=0].
WARNING. ssidat: inner friction factor on peak strength, ALPHA_FI < 0.5.
WARNING. ssidat: typically ALPHA_FI is around 0.65.
WARNING. ldata__anchor: ILC =    1; THETA_LUG_DEG = 0. [horizontal chain inclination at lug level].
WARNING. ldata__anchor: ILC =   10; THETA_LUG_DEG = 90. [vertical chain inclination at lug level].
WARNING. ldata__anchor: ILC =   10; THETA_LUG_DEG = 90 & IRO_LUG > 0. [ZLUG_CL at infinite depth].

solve_MODE22_sub: < --- BEGIN

  VHM load case #=    1; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6915.049
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000       0.000
   zlug_cL                               [m]:      10.000
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6915.049
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    -548.032    6915.049  -69150.486
   Vstar,Hstar,Mstar                     [-]:      -0.299      10.955      -7.303
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                 10.955      -7.303      -0.299
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.989
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.001
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.990
 > f_VHM_fun                             [-]: -0.2817884E-05
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00000
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    2; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6877.865
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      10.000
   zlug_cL                               [m]:      10.397
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1194.329
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6773.374
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:     646.297    6773.374  -70420.983
   Vstar,Hstar,Mstar                     [-]:       0.353      10.730      -7.437
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                 10.730      -7.437       0.353
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.972
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.011
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.984
 > f_VHM_fun                             [-]:  0.4368726E-07
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00000
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    3; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5949.903
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      20.000
   zlug_cL                               [m]:      10.819
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2034.987
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5591.080
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1486.955    5591.080  -60489.522
   Vstar,Hstar,Mstar                     [-]:       0.812       8.857      -6.389
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  8.857      -6.389       0.812
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.680
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.044
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.724
 > f_VHM_fun                             [-]:  0.1386679E-05
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00000
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    4; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    4589.035
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      30.000
   zlug_cL                               [m]:      11.299
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2294.517
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3974.220
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1746.485    3974.220  -44904.869
   Vstar,Hstar,Mstar                     [-]:       0.954       6.296      -4.743
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  6.296      -4.743       0.954
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.354
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.062
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.416
 > f_VHM_fun                             [-]:  0.2328318E-04
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00002
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    5; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3659.672
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      40.000
   zlug_cL                               [m]:      11.888
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2352.392
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2803.472
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1804.360    2803.472  -33327.599
   Vstar,Hstar,Mstar                     [-]:       0.986       4.441      -3.520
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  4.441      -3.520       0.986
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.182
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.069
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.251
 > f_VHM_fun                             [-]: -0.2098971E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00021
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    6; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3091.318
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      50.000
   zlug_cL                               [m]:      12.681
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2368.087
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1987.061
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1820.055    1987.061  -25198.807
   Vstar,Hstar,Mstar                     [-]:       0.994       3.148      -2.661
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  3.148      -2.661       0.994
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.096
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.072
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.168
 > f_VHM_fun                             [-]:  0.2549534E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00025
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    7; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2740.091
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      60.000
   zlug_cL                               [m]:      13.897
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2372.989
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1370.046
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1824.957    1370.046  -19039.680
   Vstar,Hstar,Mstar                     [-]:       0.997       2.170      -2.011
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  2.170      -2.011       0.997
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.049
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.074
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.126
 > f_VHM_fun                             [-]: -0.3586766E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00359
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    8; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2527.869
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      70.000
   zlug_cL                               [m]:      16.182
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2375.420
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     864.582
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1827.388     864.582  -13990.517
   Vstar,Hstar,Mstar                     [-]:       0.998       1.370      -1.478
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  1.370      -1.478       0.998
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.022
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.076
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.098
 > f_VHM_fun                             [-]: -0.2237295E-05
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00000
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    9; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2412.931
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      80.000
   zlug_cL                               [m]:      22.760
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2376.273
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     419.001
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1828.241     419.001   -9536.626
   Vstar,Hstar,Mstar                     [-]:       0.999       0.664      -1.007
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.664      -1.007       0.999
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.007
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.077
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.084
 > f_VHM_fun                             [-]: -0.4248829E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00042
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=   10; caisson geometry #=    1

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2376.302
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    1830.457
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2376.302
   Hlug=Tlug*cos(theta_lug_star)        [kN]:       0.000
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1828.270       0.000   -5346.680
   Vstar,Hstar,Mstar                     [-]:       0.999       0.000      -0.565
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.000      -0.565       0.999
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.001
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.078
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.084
 > f_VHM_fun                             [-]: -0.5234918E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00523
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  15.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.35    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1010.56591
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      819.89100
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    17.3    
   qv_tip = Nc * su_av_tip             [kPa]     155.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.247E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3479.84400
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1830.45691
   Qleaking   = Fo                      [kN]:     1010.56591
 > Vmax_fun                             [kN]:     1830.45691
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.39    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.18    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      548.03197
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    1; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    7826.238
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500       0.000
   zlug_cL                               [m]:      10.500
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    7826.238
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    -581.019    7826.238  -82175.497
   Vstar,Hstar,Mstar                     [-]:      -0.284      11.082      -7.272
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                 11.082      -7.272      -0.284
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.983
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.009
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.992
 > f_VHM_fun                             [-]: -0.2801461E-05
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00000
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    2; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    7834.179
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      10.000
   zlug_cL                               [m]:      10.897
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1360.391
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    7715.160
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:     779.372    7715.160  -84070.061
   Vstar,Hstar,Mstar                     [-]:       0.381      10.924      -7.440
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                 10.924      -7.440       0.381
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.977
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.002
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.979
 > f_VHM_fun                             [-]:  0.5045215E-07
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00000
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    3; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6716.189
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      20.000
   zlug_cL                               [m]:      11.319
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2297.072
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6311.153
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1716.053    6311.153  -71435.520
   Vstar,Hstar,Mstar                     [-]:       0.838       8.936      -6.322
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  8.936      -6.322       0.838
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.670
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.022
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.692
 > f_VHM_fun                             [-]:  0.1568584E-05
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00000
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    4; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5116.338
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      30.000
   zlug_cL                               [m]:      11.799
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2558.169
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4430.878
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1977.150    4430.878  -52280.104
   Vstar,Hstar,Mstar                     [-]:       0.965       6.274      -4.627
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  6.274      -4.627       0.965
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.339
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.039
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.378
 > f_VHM_fun                             [-]: -0.1814900E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00018
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    5; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    4059.378
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      40.000
   zlug_cL                               [m]:      12.388
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2609.318
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3109.664
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    2028.299    3109.664  -38522.439
   Vstar,Hstar,Mstar                     [-]:       0.990       4.403      -3.409
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  4.403      -3.409       0.990
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.173
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.048
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.221
 > f_VHM_fun                             [-]: -0.3830378E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00038
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    6; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3422.593
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      50.000
   zlug_cL                               [m]:      13.181
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2621.859
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2200.001
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    2040.840    2200.001  -28999.189
   Vstar,Hstar,Mstar                     [-]:       0.997       3.115      -2.566
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  3.115      -2.566       0.997
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.090
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.053
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.143
 > f_VHM_fun                             [-]:  0.3098091E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00031
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    7; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3031.752
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      60.000
   zlug_cL                               [m]:      14.397
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2625.574
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1515.876
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    2044.555    1515.876  -21824.242
   Vstar,Hstar,Mstar                     [-]:       0.998       2.146      -1.931
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  2.146      -1.931       0.998
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.046
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.057
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.105
 > f_VHM_fun                             [-]: -0.1439351E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00144
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    8; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2795.698
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      70.000
   zlug_cL                               [m]:      16.682
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2627.097
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     956.185
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    2046.078     956.185  -15950.912
   Vstar,Hstar,Mstar                     [-]:       0.999       1.354      -1.412
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  1.354      -1.412       0.999
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.020
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.061
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.081
 > f_VHM_fun                             [-]: -0.6036946E-05
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00001
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=    9; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2667.827
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      80.000
   zlug_cL                               [m]:      23.260
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2627.297
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     463.263
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    2046.278     463.263  -10775.682
   Vstar,Hstar,Mstar                     [-]:       0.999       0.656      -0.954
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.656      -0.954       0.999
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.006
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.064
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.077
 > f_VHM_fun                             [-]: -0.7180573E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00718
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


  VHM load case #=   10; caisson geometry #=    2

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2627.491
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2047.974
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2627.491
   Hlug=Tlug*cos(theta_lug_star)        [kN]:       0.000
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    2046.472       0.000   -5911.854
   Vstar,Hstar,Mstar                     [-]:       0.999       0.000      -0.523
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.000      -0.523       0.999
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.000
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.067
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.073
 > f_VHM_fun                             [-]: -0.6379166E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00638
root_Tlug_22_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  2
 < L,D                                   [m]:  16.00   4.50

alpha_Fo_twist: < --- BEGIN
 < psi_twist                           [deg]:   0.00
 < alpha_Fo                              [-]:     0.50960
 > alpha_Fo_twist                        [-]:     0.50960
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.51
   alpha_Fi                              [-]:   0.42
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    9.81    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1130.65363
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      917.32041
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    18.2    
   qv_tip = Nc * su_av_tip             [kPa]     164.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.260E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3733.55688
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2047.97405
   Qleaking   = Fo                      [kN]:     1130.65363
 > Vmax_fun                             [kN]:     2047.97405
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.011
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.050
   Di = D - 2 * wt_c                     [m]:  4.430
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    15.9    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    7.88    
   Vst_t =   Atop          * wt_t      [m^3]:   0.794    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.67    
   gamma_sub_steel                  [kN/m^3]:    67.0    
   DELTA_BW_USER                        [kN]:    0.00    
   Bw_caisson = gamma_sub_steel * Vst + delta_Bw_user
 > Bw_caisson_fun                       [kN]:      581.01903
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  2
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >



 -----------------------------------------------
 solve_MODE22_sub: CAISSON RESISTANCES & WEIGHTS
 MODE_VHM, problem solution mode =   22
 IPLUG_VMODE, axial capacity mode  =  2
 ISCOPE_BW, caisson weight BW mode =  1
 ilc    = VHM load case #            
 icgeom = caisson geometry #         
 D    = caisson outer diameter       
 L    = caisson embedded length      
 L/D  = caisson embedment ratio      
 Vmax = caisson axial resistance     
 Hmax = caisson FIXED HEAD lateral resistance 
 H0   = caisson FREE  HEAD lateral resistance 
 Bw   = caisson submerged weight (only used for anchor foundations)
 Bw*  = caisson+soil plug submerged weight (only for IPLUG_VMODE=5)
 su_av_L = soil average su between caisson head (seafloor) and tip (L)
 ----------------------------------------------------------------------------------------
  ilc       D      L      L/D       Vmax        Hmax         H0          Bw       su_av_L
   icgeom  [m]    [m]     [-]       [kN]        [kN]        [kN]        [kN]        [kPa]
 ----------------------------------------------------------------------------------------
    1  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    1  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

    2  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    2  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

    3  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    3  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

    4  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    4  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

    5  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    5  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

    6  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    6  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

    7  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    7  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

    8  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    8  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

    9  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
    9  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8

   10  1   4.50  15.00   3.33      1830.5      6931.9      1893.3       548.0         9.4
   10  2   4.50  16.00   3.56      2048.0      7831.9      2125.6       581.0         9.8
 ----------------------------------------------------------------------------------------
  ilc       D      L      L/D       Vmax        Hmax         H0          Bw       su_av_L
   icgeom  [m]    [m]     [-]       [kN]        [kN]        [kN]        [kN]        [kPa]
 ----------------------------------------------------------------------------------------


 ------------------------------------
 solve_MODE22_sub: SUMMARY OF RESULTS
 ilc    = VHM load case #            
 icgeom = caisson geometry #         
 D      = caisson outer diameter           [m] 
 L      = caisson embedded length          [m] 
 zlug   = caisson lug level                [m] 
 theta_lug = load inclination            [deg] 
 Tlug      = UNFACTORED anchor chain load [kN] 
 Bw        = caisson buoyant weight       [kN] 
 iters     = # iterations                  [-] 
 f_VHM     = VHM yield function            [-] 
 -----------------------------------------------------------------------
 ilc        D      L    zlug  theta    Tlug     Bw    ( iters    f_VHM )
  icgeom   [m]    [m]   [m]   [deg]    [kN]    [kN]   ( [-]      [-]   )
 -----------------------------------------------------------------------
    1  1   4.50  15.00  10.00  0.00   6915.0   548.0       5  -0.282E-05   ! ? missing text
    1  2   4.50  16.00  10.50  0.00   7826.2   581.0       5  -0.280E-05   ! ? missing text

    2  1   4.50  15.00  10.00 10.00   6877.9   548.0       6   0.437E-07   ! ? missing text
    2  2   4.50  16.00  10.50 10.00   7834.2   581.0       6   0.505E-07   ! ? missing text

    3  1   4.50  15.00  10.00 20.00   5949.9   548.0       6   0.139E-05   ! ? missing text
    3  2   4.50  16.00  10.50 20.00   6716.2   581.0       6   0.157E-05   ! ? missing text

    4  1   4.50  15.00  10.00 30.00   4589.0   548.0       5   0.233E-04   ! ? missing text
    4  2   4.50  16.00  10.50 30.00   5116.3   581.0       4  -0.181E-03   ! ? missing text

    5  1   4.50  15.00  10.00 40.00   3659.7   548.0       8  -0.210E-03   ! ? missing text
    5  2   4.50  16.00  10.50 40.00   4059.4   581.0       7  -0.383E-03   ! ? missing text

    6  1   4.50  15.00  10.00 50.00   3091.3   548.0      10   0.255E-03   ! ? missing text
    6  2   4.50  16.00  10.50 50.00   3422.6   581.0      10   0.310E-03   ! ? missing text

    7  1   4.50  15.00  10.00 60.00   2740.1   548.0       6  -0.359E-02   ! ? missing text
    7  2   4.50  16.00  10.50 60.00   3031.8   581.0       6  -0.144E-02   ! ? missing text

    8  1   4.50  15.00  10.00 70.00   2527.9   548.0      24  -0.224E-05   ! ? missing text
    8  2   4.50  16.00  10.50 70.00   2795.7   581.0      24  -0.604E-05   ! ? missing text

    9  1   4.50  15.00  10.00 80.00   2412.9   548.0       9  -0.425E-03   ! ? missing text
    9  2   4.50  16.00  10.50 80.00   2667.8   581.0      10  -0.718E-02   ! ? missing text

   10  1   4.50  15.00  10.00 90.00   2376.3   548.0      11  -0.523E-02   ! ? missing text
   10  2   4.50  16.00  10.50 90.00   2627.5   581.0      10  -0.638E-02   ! ? missing text
 -----------------------------------------------------------------------
 ilc        D      L    zlug  theta    Tlug     Bw    ( iters    f_VHM )
  icgeom   [m]    [m]   [m]   [deg]    [kN]    [kN]   ( [-]      [-]   )
 -----------------------------------------------------------------------

VHlug_resist_sub: < --- BEGIN
  load case #=    1; caisson geometry #=    1 (#rows results= =  181)
  L         [m] =   15.00
  D         [m] =    4.50
  zlug      [m] =   10.00
  iro_lug   [-] =   1: offset + D/2 (lug on outer face)
  m_fac     [-] =   1.000
  Hmax     [kN] =    0.693E+04
  Vmax     [kN] =    0.183E+04
  theta*        = resistance inclination at lug
  Vlug_theta*, Hlug_theta*
                = corresponding VH UNFACTORED resistances
  T_theta*      = sqrt( Vlug_theta*^2 + Hlug_theta*^2)
------------------------------------------------------------------------------
 theta*     Vlug_theta*  Hlug_theta*  T_theta*   Vlug_theta*  Hlug_theta*  row
  [deg]        [kN]         [kN]        [kN]      / Vmax [-]   / Hmax [-]  [-]
------------------------------------------------------------------------------
    90.0       0.183E+04   0.112E-12   0.183E+04       1.000    0.000        1
    89.0       0.183E+04    31.9       0.183E+04       1.000    0.005        2
    88.0       0.183E+04    63.9       0.183E+04       1.000    0.009        3
    87.0       0.183E+04    95.9       0.183E+04       1.000    0.014        4
    86.0       0.183E+04    128.       0.183E+04       1.000    0.018        5
    85.0       0.183E+04    160.       0.184E+04       1.000    0.023        6
    84.0       0.183E+04    192.       0.184E+04       1.000    0.028        7
    83.0       0.183E+04    225.       0.184E+04       1.000    0.032        8
    82.0       0.183E+04    257.       0.185E+04       1.000    0.037        9
    81.0       0.183E+04    290.       0.185E+04       0.999    0.042       10
    80.0       0.183E+04    323.       0.186E+04       1.000    0.047       11
    79.0       0.183E+04    356.       0.186E+04       1.000    0.051       12
    78.0       0.183E+04    389.       0.187E+04       1.000    0.056       13
    77.0       0.183E+04    422.       0.188E+04       1.000    0.061       14
    76.0       0.183E+04    456.       0.189E+04       1.000    0.066       15
    75.0       0.183E+04    490.       0.189E+04       1.000    0.071       16
    74.0       0.183E+04    525.       0.190E+04       1.000    0.076       17
    73.0       0.183E+04    559.       0.191E+04       1.000    0.081       18
    72.0       0.183E+04    594.       0.192E+04       1.000    0.086       19
    71.0       0.183E+04    630.       0.193E+04       1.000    0.091       20
    70.0       0.183E+04    666.       0.195E+04       0.999    0.096       21
    69.0       0.183E+04    702.       0.196E+04       0.999    0.101       22
    68.0       0.183E+04    739.       0.197E+04       0.999    0.107       23
    67.0       0.183E+04    777.       0.199E+04       1.000    0.112       24
    66.0       0.183E+04    814.       0.200E+04       0.999    0.117       25
    65.0       0.183E+04    853.       0.202E+04       0.999    0.123       26
    64.0       0.183E+04    892.       0.204E+04       0.999    0.129       27
    63.0       0.183E+04    932.       0.205E+04       0.999    0.134       28
    62.0       0.183E+04    973.       0.207E+04       0.999    0.140       29
    61.0       0.183E+04   0.101E+04   0.209E+04       0.999    0.146       30
    60.0       0.183E+04   0.106E+04   0.211E+04       0.999    0.152       31
    59.0       0.183E+04   0.110E+04   0.213E+04       0.999    0.158       32
    58.0       0.183E+04   0.114E+04   0.216E+04       0.999    0.165       33
    57.0       0.183E+04   0.119E+04   0.218E+04       0.999    0.171       34
    56.0       0.183E+04   0.123E+04   0.220E+04       0.999    0.178       35
    55.0       0.183E+04   0.128E+04   0.223E+04       0.999    0.185       36
    54.0       0.183E+04   0.133E+04   0.226E+04       0.999    0.192       37
    53.0       0.183E+04   0.138E+04   0.229E+04       0.999    0.199       38
    52.0       0.183E+04   0.143E+04   0.232E+04       0.998    0.206       39
    51.0       0.183E+04   0.148E+04   0.235E+04       0.998    0.213       40
    50.0       0.183E+04   0.153E+04   0.239E+04       0.998    0.221       41
    49.0       0.183E+04   0.159E+04   0.242E+04       0.998    0.229       42
    48.0       0.183E+04   0.164E+04   0.246E+04       0.998    0.237       43
    47.0       0.183E+04   0.170E+04   0.250E+04       0.998    0.246       44
    46.0       0.183E+04   0.176E+04   0.254E+04       0.998    0.254       45
    45.0       0.183E+04   0.183E+04   0.258E+04       0.997    0.263       46
    44.0       0.183E+04   0.189E+04   0.263E+04       0.997    0.273       47
    43.0       0.182E+04   0.196E+04   0.267E+04       0.997    0.282       48
    42.0       0.182E+04   0.203E+04   0.273E+04       0.996    0.292       49
    41.0       0.182E+04   0.210E+04   0.278E+04       0.996    0.303       50
    40.0       0.182E+04   0.217E+04   0.283E+04       0.995    0.313       51
    39.0       0.182E+04   0.225E+04   0.289E+04       0.995    0.324       52
    38.0       0.182E+04   0.233E+04   0.296E+04       0.994    0.336       53
    37.0       0.182E+04   0.241E+04   0.302E+04       0.994    0.348       54
    36.0       0.182E+04   0.250E+04   0.309E+04       0.993    0.361       55
    35.0       0.182E+04   0.259E+04   0.317E+04       0.992    0.374       56
    34.0       0.181E+04   0.269E+04   0.324E+04       0.991    0.388       57
    33.0       0.181E+04   0.279E+04   0.333E+04       0.990    0.403       58
    32.0       0.181E+04   0.290E+04   0.341E+04       0.989    0.418       59
    31.0       0.181E+04   0.301E+04   0.351E+04       0.987    0.434       60
    30.0       0.180E+04   0.312E+04   0.361E+04       0.985    0.450       61
    29.0       0.180E+04   0.325E+04   0.371E+04       0.983    0.468       62
    28.0       0.179E+04   0.337E+04   0.382E+04       0.980    0.487       63
    27.0       0.179E+04   0.351E+04   0.394E+04       0.977    0.506       64
    26.0       0.178E+04   0.365E+04   0.406E+04       0.973    0.527       65
    25.0       0.177E+04   0.380E+04   0.419E+04       0.968    0.548       66
    24.0       0.176E+04   0.396E+04   0.433E+04       0.963    0.571       67
    23.0       0.175E+04   0.412E+04   0.448E+04       0.956    0.595       68
    22.0       0.174E+04   0.430E+04   0.463E+04       0.948    0.620       69
    21.0       0.172E+04   0.448E+04   0.480E+04       0.939    0.646       70
    20.0       0.170E+04   0.466E+04   0.496E+04       0.927    0.673       71
    19.0       0.167E+04   0.486E+04   0.513E+04       0.913    0.700       72
    18.0       0.164E+04   0.505E+04   0.531E+04       0.897    0.729       73
    17.0       0.160E+04   0.525E+04   0.549E+04       0.877    0.757       74
    16.0       0.156E+04   0.545E+04   0.567E+04       0.853    0.786       75
    15.0       0.151E+04   0.564E+04   0.584E+04       0.826    0.814       76
    14.0       0.145E+04   0.583E+04   0.601E+04       0.794    0.841       77
    13.0       0.139E+04   0.601E+04   0.617E+04       0.758    0.867       78
    12.0       0.131E+04   0.618E+04   0.631E+04       0.717    0.891       79
    11.0       0.123E+04   0.633E+04   0.644E+04       0.672    0.913       80
    10.0       0.114E+04   0.646E+04   0.656E+04       0.622    0.932       81
    9.00       0.104E+04   0.658E+04   0.666E+04       0.569    0.949       82
    8.00        938.       0.667E+04   0.674E+04       0.512    0.963       83
    7.00        829.       0.675E+04   0.680E+04       0.453    0.974       84
    6.00        716.       0.681E+04   0.685E+04       0.391    0.983       85
    5.00        600.       0.686E+04   0.689E+04       0.328    0.990       86
    4.00        482.       0.690E+04   0.691E+04       0.263    0.995       87
    3.00        363.       0.692E+04   0.693E+04       0.198    0.998       88
    2.00        242.       0.693E+04   0.694E+04       0.132    1.000       89
    1.00        121.       0.694E+04   0.694E+04       0.066    1.002       90
    0.00        0.00       0.695E+04   0.695E+04       0.000    1.002       91
   -1.00       -121.       0.695E+04   0.695E+04      -0.066    1.003       92
   -2.00       -243.       0.695E+04   0.696E+04      -0.133    1.003       93
   -3.00       -364.       0.695E+04   0.696E+04      -0.199    1.002       94
   -4.00       -485.       0.693E+04   0.695E+04      -0.265    1.000       95
   -5.00       -604.       0.691E+04   0.694E+04      -0.330    0.997       96
   -6.00       -722.       0.687E+04   0.691E+04      -0.395    0.991       97
   -7.00       -837.       0.682E+04   0.687E+04      -0.457    0.983       98
   -8.00       -948.       0.674E+04   0.681E+04      -0.518    0.973       99
   -9.00      -0.105E+04   0.665E+04   0.673E+04      -0.575    0.959      100
   -10.0      -0.115E+04   0.654E+04   0.664E+04      -0.630    0.943      101
   -11.0      -0.124E+04   0.640E+04   0.652E+04      -0.680    0.924      102
   -12.0      -0.133E+04   0.625E+04   0.639E+04      -0.726    0.902      103
   -13.0      -0.140E+04   0.608E+04   0.624E+04      -0.767    0.877      104
   -14.0      -0.147E+04   0.590E+04   0.608E+04      -0.803    0.851      105
   -15.0      -0.153E+04   0.570E+04   0.591E+04      -0.835    0.823      106
   -16.0      -0.158E+04   0.550E+04   0.573E+04      -0.862    0.794      107
   -17.0      -0.162E+04   0.530E+04   0.554E+04      -0.885    0.765      108
   -18.0      -0.166E+04   0.510E+04   0.536E+04      -0.905    0.735      109
   -19.0      -0.169E+04   0.489E+04   0.518E+04      -0.921    0.706      110
   -20.0      -0.171E+04   0.470E+04   0.500E+04      -0.934    0.678      111
   -21.0      -0.173E+04   0.451E+04   0.483E+04      -0.945    0.650      112
   -22.0      -0.175E+04   0.432E+04   0.466E+04      -0.954    0.623      113
   -23.0      -0.176E+04   0.415E+04   0.450E+04      -0.961    0.598      114
   -24.0      -0.177E+04   0.398E+04   0.435E+04      -0.967    0.574      115
   -25.0      -0.178E+04   0.382E+04   0.421E+04      -0.972    0.551      116
   -26.0      -0.179E+04   0.367E+04   0.408E+04      -0.977    0.529      117
   -27.0      -0.179E+04   0.352E+04   0.395E+04      -0.980    0.508      118
   -28.0      -0.180E+04   0.338E+04   0.383E+04      -0.983    0.488      119
   -29.0      -0.180E+04   0.325E+04   0.372E+04      -0.985    0.469      120
   -30.0      -0.181E+04   0.313E+04   0.361E+04      -0.987    0.452      121
   -31.0      -0.181E+04   0.301E+04   0.351E+04      -0.989    0.435      122
   -32.0      -0.181E+04   0.290E+04   0.342E+04      -0.990    0.418      123
   -33.0      -0.181E+04   0.279E+04   0.333E+04      -0.992    0.403      124
   -34.0      -0.182E+04   0.269E+04   0.325E+04      -0.993    0.389      125
   -35.0      -0.182E+04   0.260E+04   0.317E+04      -0.993    0.375      126
   -36.0      -0.182E+04   0.250E+04   0.310E+04      -0.994    0.361      127
   -37.0      -0.182E+04   0.242E+04   0.303E+04      -0.995    0.349      128
   -38.0      -0.182E+04   0.233E+04   0.296E+04      -0.995    0.336      129
   -39.0      -0.182E+04   0.225E+04   0.290E+04      -0.996    0.325      130
   -40.0      -0.182E+04   0.217E+04   0.284E+04      -0.996    0.314      131
   -41.0      -0.182E+04   0.210E+04   0.278E+04      -0.997    0.303      132
   -42.0      -0.183E+04   0.203E+04   0.273E+04      -0.997    0.292      133
   -43.0      -0.183E+04   0.196E+04   0.268E+04      -0.997    0.282      134
   -44.0      -0.183E+04   0.189E+04   0.263E+04      -0.998    0.273      135
   -45.0      -0.183E+04   0.183E+04   0.258E+04      -0.998    0.263      136
   -46.0      -0.183E+04   0.176E+04   0.254E+04      -0.998    0.254      137
   -47.0      -0.183E+04   0.170E+04   0.250E+04      -0.998    0.246      138
   -48.0      -0.183E+04   0.165E+04   0.246E+04      -0.998    0.237      139
   -49.0      -0.183E+04   0.159E+04   0.242E+04      -0.998    0.229      140
   -50.0      -0.183E+04   0.153E+04   0.239E+04      -0.999    0.221      141
   -51.0      -0.183E+04   0.148E+04   0.235E+04      -0.999    0.214      142
   -52.0      -0.183E+04   0.143E+04   0.232E+04      -0.999    0.206      143
   -53.0      -0.183E+04   0.138E+04   0.229E+04      -0.999    0.199      144
   -54.0      -0.183E+04   0.133E+04   0.226E+04      -0.999    0.192      145
   -55.0      -0.183E+04   0.128E+04   0.223E+04      -0.999    0.185      146
   -56.0      -0.183E+04   0.123E+04   0.220E+04      -0.999    0.178      147
   -57.0      -0.183E+04   0.119E+04   0.218E+04      -0.999    0.171      148
   -58.0      -0.183E+04   0.114E+04   0.216E+04      -0.999    0.165      149
   -59.0      -0.183E+04   0.110E+04   0.213E+04      -0.999    0.159      150
   -60.0      -0.183E+04   0.106E+04   0.211E+04      -0.999    0.152      151
   -61.0      -0.183E+04   0.101E+04   0.209E+04      -0.999    0.146      152
   -62.0      -0.183E+04    973.       0.207E+04      -0.999    0.140      153
   -63.0      -0.183E+04    932.       0.205E+04      -0.999    0.134      154
   -64.0      -0.183E+04    892.       0.204E+04      -0.999    0.129      155
   -65.0      -0.183E+04    853.       0.202E+04      -1.000    0.123      156
   -66.0      -0.183E+04    814.       0.200E+04      -0.999    0.117      157
   -67.0      -0.183E+04    777.       0.199E+04      -1.000    0.112      158
   -68.0      -0.183E+04    739.       0.197E+04      -1.000    0.107      159
   -69.0      -0.183E+04    702.       0.196E+04      -1.000    0.101      160
   -70.0      -0.183E+04    666.       0.195E+04      -0.999    0.096      161
   -71.0      -0.183E+04    630.       0.194E+04      -1.000    0.091      162
   -72.0      -0.183E+04    595.       0.192E+04      -1.000    0.086      163
   -73.0      -0.183E+04    559.       0.191E+04      -1.000    0.081      164
   -74.0      -0.183E+04    525.       0.190E+04      -1.000    0.076      165
   -75.0      -0.183E+04    490.       0.189E+04      -1.000    0.071      166
   -76.0      -0.183E+04    456.       0.189E+04      -1.000    0.066      167
   -77.0      -0.183E+04    422.       0.188E+04      -1.000    0.061      168
   -78.0      -0.183E+04    389.       0.187E+04      -1.000    0.056      169
   -79.0      -0.183E+04    356.       0.186E+04      -0.999    0.051      170
   -80.0      -0.183E+04    323.       0.186E+04      -1.000    0.047      171
   -81.0      -0.183E+04    290.       0.185E+04      -1.000    0.042      172
   -82.0      -0.183E+04    257.       0.185E+04      -1.000    0.037      173
   -83.0      -0.183E+04    225.       0.184E+04      -1.000    0.032      174
   -84.0      -0.183E+04    192.       0.184E+04      -1.000    0.028      175
   -85.0      -0.183E+04    160.       0.184E+04      -1.000    0.023      176
   -86.0      -0.183E+04    128.       0.183E+04      -1.000    0.018      177
   -87.0      -0.183E+04    95.9       0.183E+04      -1.000    0.014      178
   -88.0      -0.183E+04    63.9       0.183E+04      -1.000    0.009      179
   -89.0      -0.183E+04    31.9       0.183E+04      -1.000    0.005      180
   -90.0      -0.183E+04   0.112E-12   0.183E+04      -1.000    0.000      181
------------------------------------------------------------------------------
 theta*     Vlug_theta*  Hlug_theta*  T_theta*   Vlug_theta*  Hlug_theta*  row
  [deg]        [kN]         [kN]        [kN]      / Vmax [-]   / Hmax [-]  [-]
------------------------------------------------------------------------------
VHlug_resist_sub: END --- >

solve_MODE22_sub: END --- >


caisson_vhm: input   file: "JOB.I01"
caisson_vhm: output  file: "JOB.O01"
caisson_vhm: normal program end
