
 p r o g r a m  C A I S S O N _ V H M  w32 version 00.12_b04 
 ------------------------------------  at 15:22:53 on 19-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.30    
 gmdata: END --- >


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


 cfdata: < --- BEGIN
 ncgeom
    3
 cgeom
 CIRCULAR
 nfhxp  wt_D_side  wt_D_top
  0   0.780E-02   0.780E-02
 iro_lug  r_lug_user
  1    2.25    
 D  L  zlug
   4.50  17.00  11.50
   4.50  16.00  10.50
   4.50  15.00  10.00
 cfdata: END --- >


 ssidat: < --- BEGIN
 alpha_Fo  alpha_Fi  eta_Fo
   0.65   0.54   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
    2    1.00        1.00        0.00        665.    
 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 :    6
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. chndat: user-supplied EXTRA buoyant weight, DELTA_BW_USER <= 0 kN.
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]:    6733.996
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500       0.000
   zlug_cL                               [m]:      11.500
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6733.996
   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]:    -598.182    6733.996  -77440.957
   Vstar,Hstar,Mstar                     [-]:      -0.267      11.145      -7.539
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                 11.145      -7.539      -0.267
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.993
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.000
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.993
 > f_VHM_fun                             [-]: -0.2986644E-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]:  17.00   4.50

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    6703.145
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      10.000
   zlug_cL                               [m]:      11.897
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1163.989
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6601.309
   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]:     565.807    6601.309  -78534.029
   Vstar,Hstar,Mstar                     [-]:       0.252      10.926      -7.646
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                 10.926      -7.646       0.252
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.975
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.019
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.994
 > f_VHM_fun                             [-]: -0.7860802E-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]:  17.00   4.50

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    6267.315
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      20.000
   zlug_cL                               [m]:      12.319
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2143.548
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5889.350
   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]:    1545.366    5889.350  -72550.505
   Vstar,Hstar,Mstar                     [-]:       0.689       9.747      -7.063
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                  9.747      -7.063       0.689
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.793
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.059
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.852
 > f_VHM_fun                             [-]:  0.4081913E-06
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]:  17.00   4.50

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    5260.447
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      30.000
   zlug_cL                               [m]:      12.799
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2630.223
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4555.681
   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]:    2032.041    4555.681  -58308.331
   Vstar,Hstar,Mstar                     [-]:       0.906       7.540      -5.677
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                  7.540      -5.677       0.906
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.487
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.086
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.573
 > f_VHM_fun                             [-]: -0.5487518E-04
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00005
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    4320.706
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      40.000
   zlug_cL                               [m]:      13.388
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2777.296
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3309.853
   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]:    2179.114    3309.853  -44312.224
   Vstar,Hstar,Mstar                     [-]:       0.972       5.478      -4.314
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                  5.478      -4.314       0.972
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.265
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.094
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.359
 > f_VHM_fun                             [-]: -0.4203294E-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]:  17.00   4.50

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    3678.666
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      50.000
   zlug_cL                               [m]:      14.181
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2818.021
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2364.601
   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]:    2219.839    2364.601  -33533.456
   Vstar,Hstar,Mstar                     [-]:       0.990       3.914      -3.265
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                  3.914      -3.265       0.990
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.141
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.095
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.237
 > f_VHM_fun                             [-]: -0.7813161E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00078
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    3268.992
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      60.000
   zlug_cL                               [m]:      15.397
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2831.030
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1634.496
   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]:    2232.848    1634.496  -25166.523
   Vstar,Hstar,Mstar                     [-]:       0.996       2.705      -2.450
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                  2.705      -2.450       0.996
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.072
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.095
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.168
 > f_VHM_fun                             [-]: -0.1300046E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00130
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    3017.367
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      70.000
   zlug_cL                               [m]:      17.682
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2835.397
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1032.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]:    2237.215    1032.000  -18247.646
   Vstar,Hstar,Mstar                     [-]:       0.998       1.708      -1.777
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                  1.708      -1.777       0.998
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.032
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.094
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.132
 > f_VHM_fun                             [-]: -0.5918779E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00592
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    2881.646
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      80.000
   zlug_cL                               [m]:      24.260
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2837.867
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     500.393
   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]:    2239.685     500.393  -12139.715
   Vstar,Hstar,Mstar                     [-]:       0.999       0.828      -1.182
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.828      -1.182       0.999
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.010
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.094
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.104
 > f_VHM_fun                             [-]: -0.5650790E-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]:  17.00   4.50

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    2838.538
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2242.815
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2838.538
   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]:    2240.356       0.000   -6386.710
   Vstar,Hstar,Mstar                     [-]:       0.999       0.000      -0.622
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  17.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.000      -0.622       0.999
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.506       1.688     -33.976
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.506       1.688     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.001
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.093
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.094
 > f_VHM_fun                             [-]: -0.6065223E-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]:  17.00   4.50

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   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.211E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3341.24941
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2242.81484
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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]:    8.37    
   Vst_t =   Atop          * wt_t      [m^3]:   0.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.93    
   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]:      598.18209
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.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]:    6021.143
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    6021.143
   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]:    -565.195    6021.143  -63221.998
   Vstar,Hstar,Mstar                     [-]:      -0.280      11.083      -7.273
   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.083      -7.273      -0.280
   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.2959725E-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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    6075.664
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    1055.028
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5983.361
   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]:     489.833    5983.361  -65199.102
   Vstar,Hstar,Mstar                     [-]:       0.243      11.014      -7.501
   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.014      -7.501       0.243
   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.993
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.002
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.995
 > f_VHM_fun                             [-]: -0.7413217E-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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    5725.750
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    1958.322
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5380.445
   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]:    1393.127    5380.445  -60900.895
   Vstar,Hstar,Mstar                     [-]:       0.691       9.904      -7.006
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  9.904      -7.006       0.691
   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.823
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.027
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.850
 > f_VHM_fun                             [-]:  0.6039861E-06
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    4792.512
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    2396.256
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4150.437
   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]:    1831.061    4150.437  -48971.169
   Vstar,Hstar,Mstar                     [-]:       0.908       7.640      -5.634
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  7.640      -5.634       0.908
   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.503
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.058
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.561
 > f_VHM_fun                             [-]: -0.5106058E-04
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00005
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    3927.301
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    2524.420
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3008.487
   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]:    1959.225    3008.487  -37269.058
   Vstar,Hstar,Mstar                     [-]:       0.972       5.538      -4.288
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                  5.538      -4.288       0.972
   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.273
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.076
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.349
 > f_VHM_fun                             [-]: -0.2827192E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00028
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    3342.353
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    2560.391
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2148.423
   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]:    1995.196    2148.423  -28319.322
   Vstar,Hstar,Mstar                     [-]:       0.989       3.955      -3.258
   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.955      -3.258       0.989
   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.145
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.086
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.230
 > f_VHM_fun                             [-]:  0.9882723E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00099
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    2968.942
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    2571.179
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1484.471
   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]:    2005.984    1484.471  -21372.095
   Vstar,Hstar,Mstar                     [-]:       0.995       2.733      -2.459
   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.733      -2.459       0.995
   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.074
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.093
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.171
 > f_VHM_fun                             [-]: -0.3799355E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00380
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    2741.156
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    2575.844
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     937.531
   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]:    2010.649     937.531  -15639.721
   Vstar,Hstar,Mstar                     [-]:       0.997       1.726      -1.799
   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.726      -1.799       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.033
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.099
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.133
 > f_VHM_fun                             [-]: -0.1499221E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00150
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    2617.434
   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]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    2577.670
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     454.513
   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]:    2012.475     454.513  -10572.140
   Vstar,Hstar,Mstar                     [-]:       0.998       0.837      -1.216
   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.837      -1.216       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.011
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.104
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.114
 > f_VHM_fun                             [-]:  0.2293059E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00023
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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]:    2578.066
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     565.195       0.000
   su_av_L                             [kPa]:       7.545
   Vmax                                 [kN]:    2016.594
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   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]:    2578.066
   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]:    2012.871       0.000   -5800.649
   Vstar,Hstar,Mstar                     [-]:       0.998       0.000      -0.667
   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.667       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.001
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.108
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.109
 > f_VHM_fun                             [-]: -0.5127306E-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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1109.35403
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   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]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   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.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3111.58730
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2016.59400
   Qleaking   = Fo                      [kN]:     1109.35403
 > Vmax_fun                             [kN]:     2016.59400
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    8.44    
   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]:      565.19503
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 #=    1; caisson geometry #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5320.390
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    5320.390
   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]:    -532.208    5320.390  -53203.898
   Vstar,Hstar,Mstar                     [-]:      -0.295      10.957      -7.305
   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.957      -7.305      -0.295
   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.990
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.001
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.991
 > f_VHM_fun                             [-]: -0.2975645E-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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5325.125
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:     924.698
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5244.224
   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]:     392.490    5244.224  -54522.811
   Vstar,Hstar,Mstar                     [-]:       0.218      10.800      -7.486
   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.800      -7.486       0.218
   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.985
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.011
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.997
 > f_VHM_fun                             [-]: -0.6858782E-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]:  15.00   4.50

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5023.859
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    1718.261
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4720.883
   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]:    1186.053    4720.883  -51074.919
   Vstar,Hstar,Mstar                     [-]:       0.658       9.723      -7.013
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  9.723      -7.013       0.658
   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.819
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.054
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.873
 > f_VHM_fun                             [-]:  0.4809562E-06
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    4254.986
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    2127.493
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3684.926
   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]:    1595.285    3684.926  -41636.123
   Vstar,Hstar,Mstar                     [-]:       0.885       7.589      -5.717
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  7.589      -5.717       0.885
   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.514
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.091
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.604
 > f_VHM_fun                             [-]: -0.1009085E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00010
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3519.695
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    2262.416
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2696.242
   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]:    1730.208    2696.242  -32052.861
   Vstar,Hstar,Mstar                     [-]:       0.960       5.553      -4.401
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  5.553      -4.401       0.960
   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.285
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.108
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.393
 > f_VHM_fun                             [-]: -0.2307044E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00023
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3008.513
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    2304.654
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1933.835
   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]:    1772.447    1933.835  -24523.820
   Vstar,Hstar,Mstar                     [-]:       0.983       3.983      -3.367
   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.983      -3.367       0.983
   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.153
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.115
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.269
 > f_VHM_fun                             [-]:  0.6420457E-04
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00006
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2677.318
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    2318.625
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1338.659
   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]:    1786.417    1338.659  -18603.496
   Vstar,Hstar,Mstar                     [-]:       0.991       2.757      -2.554
   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.757      -2.554       0.991
   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.079
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.120
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.204
 > f_VHM_fun                             [-]: -0.5773618E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00577
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2474.456
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    2325.228
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     846.314
   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]:    1793.020     846.314  -13694.902
   Vstar,Hstar,Mstar                     [-]:       0.995       1.743      -1.880
   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.743      -1.880       0.995
   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.036
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.123
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.161
 > f_VHM_fun                             [-]: -0.3213832E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00321
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2364.014
   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]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    2328.099
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     410.507
   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]:    1795.891     410.507   -9343.290
   Vstar,Hstar,Mstar                     [-]:       0.996       0.845      -1.283
   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.845      -1.283       0.996
   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.011
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.125
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.137
 > f_VHM_fun                             [-]: -0.8929707E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00089
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 #=    3

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2329.471
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     532.208       0.000
   su_av_L                             [kPa]:       7.193
   Vmax                                 [kN]:    1802.410
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   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]:    2329.471
   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]:    1797.263       0.000   -5241.310
   Vstar,Hstar,Mstar                     [-]:       0.997       0.000      -0.720
   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.720       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.001
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.127
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.124
 > f_VHM_fun                             [-]:  0.4291379E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00429
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.65000
 > alpha_Fo_twist                        [-]:     0.65000
alpha_Fo_twist: END --- >

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   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]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      991.52856
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   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]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   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.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2890.97325
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1802.40977
   Qleaking   = Fo                      [kN]:      991.52856
 > Vmax_fun                             [kN]:     1802.40977
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  15.00   4.50
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.008
   wt_D_top                              [-]:  0.008
   wt_c = D * wt_D_side                  [m]:  0.035
   wt_x = wt_c                           [m]:  0.035
   wt_t = D * wt_D_top                   [m]:  0.035
   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.558    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    7.94    
   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]:      532.20798
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 --- >



 -----------------------------------------------
 solve_MODE22_sub: CAISSON RESISTANCES & WEIGHTS
 MODE_VHM, problem solution mode =   22
 IPLUG_VMODE, axial capacity mode  =  2
 ISCOPE_BW, caisson weight BW mode =  2
 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  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    1  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    1  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

    2  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    2  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    2  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

    3  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    3  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    3  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

    4  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    4  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    4  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

    5  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    5  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    5  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

    6  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    6  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    6  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

    7  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    7  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    7  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

    8  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    8  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    8  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

    9  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
    9  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
    9  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2

   10  1   4.50  17.00   3.78      2242.8      6767.3      1825.2       598.2         7.9
   10  2   4.50  16.00   3.56      2016.6      6024.5      1635.1       565.2         7.5
   10  3   4.50  15.00   3.33      1802.4      5332.3      1456.4       532.2         7.2
 ----------------------------------------------------------------------------------------
  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  17.00  11.50  0.00   6734.0   598.2       5  -0.299E-05   ! ? missing text
    1  2   4.50  16.00  10.50  0.00   6021.1   565.2       5  -0.296E-05   ! ? missing text
    1  3   4.50  15.00  10.00  0.00   5320.4   532.2       5  -0.298E-05   ! ? missing text

    2  1   4.50  17.00  11.50 10.00   6703.1   598.2       5  -0.786E-05   ! ? missing text
    2  2   4.50  16.00  10.50 10.00   6075.7   565.2       5  -0.741E-05   ! ? missing text
    2  3   4.50  15.00  10.00 10.00   5325.1   532.2       5  -0.686E-05   ! ? missing text

    3  1   4.50  17.00  11.50 20.00   6267.3   598.2       6   0.408E-06   ! ? missing text
    3  2   4.50  16.00  10.50 20.00   5725.7   565.2       6   0.604E-06   ! ? missing text
    3  3   4.50  15.00  10.00 20.00   5023.9   532.2       6   0.481E-06   ! ? missing text

    4  1   4.50  17.00  11.50 30.00   5260.4   598.2       5  -0.549E-04   ! ? missing text
    4  2   4.50  16.00  10.50 30.00   4792.5   565.2       5  -0.511E-04   ! ? missing text
    4  3   4.50  15.00  10.00 30.00   4255.0   532.2       5  -0.101E-03   ! ? missing text

    5  1   4.50  17.00  11.50 40.00   4320.7   598.2       6  -0.420E-05   ! ? missing text
    5  2   4.50  16.00  10.50 40.00   3927.3   565.2       5  -0.283E-03   ! ? missing text
    5  3   4.50  15.00  10.00 40.00   3519.7   532.2       4  -0.231E-03   ! ? missing text

    6  1   4.50  17.00  11.50 50.00   3678.7   598.2       6  -0.781E-03   ! ? missing text
    6  2   4.50  16.00  10.50 50.00   3342.4   565.2       6   0.988E-03   ! ? missing text
    6  3   4.50  15.00  10.00 50.00   3008.5   532.2       7   0.642E-04   ! ? missing text

    7  1   4.50  17.00  11.50 60.00   3269.0   598.2       9  -0.130E-02   ! ? missing text
    7  2   4.50  16.00  10.50 60.00   2968.9   565.2       8  -0.380E-02   ! ? missing text
    7  3   4.50  15.00  10.00 60.00   2677.3   532.2       8  -0.577E-02   ! ? missing text

    8  1   4.50  17.00  11.50 70.00   3017.4   598.2      10  -0.592E-02   ! ? missing text
    8  2   4.50  16.00  10.50 70.00   2741.2   565.2      10  -0.150E-02   ! ? missing text
    8  3   4.50  15.00  10.00 70.00   2474.5   532.2       8  -0.321E-02   ! ? missing text

    9  1   4.50  17.00  11.50 80.00   2881.6   598.2      24  -0.565E-05   ! ? missing text
    9  2   4.50  16.00  10.50 80.00   2617.4   565.2       8   0.229E-03   ! ? missing text
    9  3   4.50  15.00  10.00 80.00   2364.0   532.2       9  -0.893E-03   ! ? missing text

   10  1   4.50  17.00  11.50 90.00   2838.5   598.2      24  -0.607E-05   ! ? missing text
   10  2   4.50  16.00  10.50 90.00   2578.1   565.2      24  -0.513E-05   ! ? missing text
   10  3   4.50  15.00  10.00 90.00   2329.5   532.2      10   0.429E-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] =   17.00
  D         [m] =    4.50
  zlug      [m] =   11.50
  iro_lug   [-] =   1: offset + D/2 (lug on outer face)
  m_fac     [-] =   1.300
  Hmax     [kN] =    0.677E+04
  Vmax     [kN] =    0.224E+04
  theta*        = resistance inclination at lug
  Vlug_theta*, Hlug_theta*
                = corresponding VH FACTORED 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.224E+04   0.137E-12   0.224E+04       1.000    0.000        1
    89.0       0.224E+04    39.2       0.224E+04       1.000    0.006        2
    88.0       0.224E+04    78.3       0.224E+04       1.000    0.012        3
    87.0       0.224E+04    117.       0.224E+04       0.999    0.017        4
    86.0       0.224E+04    157.       0.225E+04       1.000    0.023        5
    85.0       0.224E+04    196.       0.225E+04       0.999    0.029        6
    84.0       0.224E+04    236.       0.225E+04       1.000    0.035        7
    83.0       0.224E+04    275.       0.226E+04       1.000    0.041        8
    82.0       0.224E+04    315.       0.226E+04       0.999    0.047        9
    81.0       0.224E+04    355.       0.227E+04       1.000    0.052       10
    80.0       0.224E+04    395.       0.228E+04       1.000    0.058       11
    79.0       0.224E+04    436.       0.228E+04       0.999    0.064       12
    78.0       0.224E+04    476.       0.229E+04       0.999    0.070       13
    77.0       0.224E+04    517.       0.230E+04       0.999    0.076       14
    76.0       0.224E+04    559.       0.231E+04       1.000    0.083       15
    75.0       0.224E+04    601.       0.232E+04       1.000    0.089       16
    74.0       0.224E+04    643.       0.233E+04       0.999    0.095       17
    73.0       0.224E+04    685.       0.234E+04       1.000    0.101       18
    72.0       0.224E+04    728.       0.236E+04       0.999    0.108       19
    71.0       0.224E+04    772.       0.237E+04       0.999    0.114       20
    70.0       0.224E+04    816.       0.239E+04       0.999    0.121       21
    69.0       0.224E+04    861.       0.240E+04       1.000    0.127       22
    68.0       0.224E+04    905.       0.242E+04       0.999    0.134       23
    67.0       0.224E+04    951.       0.243E+04       0.999    0.141       24
    66.0       0.224E+04    998.       0.245E+04       0.999    0.147       25
    65.0       0.224E+04   0.104E+04   0.247E+04       0.999    0.154       26
    64.0       0.224E+04   0.109E+04   0.249E+04       0.999    0.161       27
    63.0       0.224E+04   0.114E+04   0.251E+04       0.999    0.169       28
    62.0       0.224E+04   0.119E+04   0.254E+04       0.999    0.176       29
    61.0       0.224E+04   0.124E+04   0.256E+04       0.999    0.183       30
    60.0       0.224E+04   0.129E+04   0.259E+04       0.998    0.191       31
    59.0       0.224E+04   0.134E+04   0.261E+04       0.998    0.199       32
    58.0       0.224E+04   0.140E+04   0.264E+04       0.998    0.207       33
    57.0       0.224E+04   0.145E+04   0.267E+04       0.998    0.215       34
    56.0       0.224E+04   0.151E+04   0.270E+04       0.998    0.223       35
    55.0       0.224E+04   0.157E+04   0.273E+04       0.998    0.232       36
    54.0       0.224E+04   0.163E+04   0.277E+04       0.998    0.240       37
    53.0       0.224E+04   0.169E+04   0.280E+04       0.997    0.249       38
    52.0       0.224E+04   0.175E+04   0.284E+04       0.997    0.258       39
    51.0       0.224E+04   0.181E+04   0.288E+04       0.997    0.268       40
    50.0       0.224E+04   0.188E+04   0.292E+04       0.997    0.277       41
    49.0       0.223E+04   0.194E+04   0.296E+04       0.996    0.287       42
    48.0       0.223E+04   0.201E+04   0.301E+04       0.996    0.297       43
    47.0       0.223E+04   0.208E+04   0.305E+04       0.995    0.308       44
    46.0       0.223E+04   0.216E+04   0.310E+04       0.995    0.318       45
    45.0       0.223E+04   0.223E+04   0.315E+04       0.995    0.330       46
    44.0       0.223E+04   0.231E+04   0.321E+04       0.994    0.341       47
    43.0       0.223E+04   0.239E+04   0.327E+04       0.993    0.353       48
    42.0       0.223E+04   0.247E+04   0.333E+04       0.992    0.365       49
    41.0       0.222E+04   0.256E+04   0.339E+04       0.992    0.378       50
    40.0       0.222E+04   0.265E+04   0.346E+04       0.991    0.391       51
    39.0       0.222E+04   0.274E+04   0.353E+04       0.989    0.405       52
    38.0       0.222E+04   0.284E+04   0.360E+04       0.988    0.419       53
    37.0       0.221E+04   0.294E+04   0.368E+04       0.986    0.434       54
    36.0       0.221E+04   0.304E+04   0.376E+04       0.984    0.449       55
    35.0       0.220E+04   0.315E+04   0.384E+04       0.982    0.465       56
    34.0       0.220E+04   0.326E+04   0.393E+04       0.980    0.481       57
    33.0       0.219E+04   0.337E+04   0.402E+04       0.977    0.499       58
    32.0       0.218E+04   0.350E+04   0.412E+04       0.974    0.516       59
    31.0       0.218E+04   0.362E+04   0.422E+04       0.970    0.535       60
    30.0       0.217E+04   0.375E+04   0.433E+04       0.965    0.554       61
    29.0       0.215E+04   0.388E+04   0.444E+04       0.960    0.574       62
    28.0       0.214E+04   0.402E+04   0.456E+04       0.954    0.595       63
    27.0       0.212E+04   0.417E+04   0.468E+04       0.947    0.616       64
    26.0       0.210E+04   0.431E+04   0.480E+04       0.938    0.638       65
    25.0       0.208E+04   0.447E+04   0.493E+04       0.928    0.660       66
    24.0       0.206E+04   0.462E+04   0.506E+04       0.917    0.683       67
    23.0       0.203E+04   0.478E+04   0.519E+04       0.904    0.706       68
    22.0       0.199E+04   0.493E+04   0.532E+04       0.889    0.729       69
    21.0       0.195E+04   0.509E+04   0.545E+04       0.871    0.752       70
    20.0       0.191E+04   0.525E+04   0.558E+04       0.851    0.775       71
    19.0       0.186E+04   0.540E+04   0.571E+04       0.829    0.798       72
    18.0       0.180E+04   0.555E+04   0.583E+04       0.804    0.820       73
    17.0       0.174E+04   0.569E+04   0.595E+04       0.776    0.841       74
    16.0       0.167E+04   0.583E+04   0.606E+04       0.745    0.861       75
    15.0       0.160E+04   0.595E+04   0.616E+04       0.711    0.880       76
    14.0       0.151E+04   0.607E+04   0.626E+04       0.675    0.897       77
    13.0       0.143E+04   0.618E+04   0.634E+04       0.636    0.913       78
    12.0       0.133E+04   0.628E+04   0.642E+04       0.595    0.928       79
    11.0       0.124E+04   0.637E+04   0.648E+04       0.552    0.941       80
    10.0       0.114E+04   0.644E+04   0.654E+04       0.506    0.952       81
    9.00       0.103E+04   0.651E+04   0.659E+04       0.460    0.962       82
    8.00        923.       0.656E+04   0.663E+04       0.411    0.970       83
    7.00        812.       0.661E+04   0.666E+04       0.362    0.977       84
    6.00        699.       0.665E+04   0.669E+04       0.312    0.983       85
    5.00        584.       0.668E+04   0.671E+04       0.261    0.987       86
    4.00        469.       0.670E+04   0.672E+04       0.209    0.991       87
    3.00        352.       0.672E+04   0.673E+04       0.157    0.994       88
    2.00        235.       0.674E+04   0.674E+04       0.105    0.996       89
    1.00        118.       0.675E+04   0.675E+04       0.053    0.997       90
    0.00        0.00       0.676E+04   0.676E+04       0.000    0.999       91
   -1.00       -118.       0.677E+04   0.677E+04      -0.053    1.000       92
   -2.00       -236.       0.677E+04   0.678E+04      -0.105    1.001       93
   -3.00       -355.       0.677E+04   0.678E+04      -0.158    1.001       94
   -4.00       -474.       0.677E+04   0.679E+04      -0.211    1.001       95
   -5.00       -592.       0.676E+04   0.679E+04      -0.264    0.999       96
   -6.00       -709.       0.675E+04   0.679E+04      -0.316    0.997       97
   -7.00       -826.       0.672E+04   0.677E+04      -0.368    0.994       98
   -8.00       -940.       0.669E+04   0.676E+04      -0.419    0.989       99
   -9.00      -0.105E+04   0.665E+04   0.673E+04      -0.469    0.982      100
   -10.0      -0.116E+04   0.659E+04   0.669E+04      -0.518    0.974      101
   -11.0      -0.127E+04   0.652E+04   0.664E+04      -0.565    0.964      102
   -12.0      -0.137E+04   0.644E+04   0.658E+04      -0.610    0.952      103
   -13.0      -0.146E+04   0.634E+04   0.651E+04      -0.653    0.938      104
   -14.0      -0.155E+04   0.624E+04   0.643E+04      -0.693    0.922      105
   -15.0      -0.164E+04   0.612E+04   0.633E+04      -0.731    0.904      106
   -16.0      -0.172E+04   0.598E+04   0.623E+04      -0.765    0.884      107
   -17.0      -0.179E+04   0.584E+04   0.611E+04      -0.796    0.863      108
   -18.0      -0.185E+04   0.569E+04   0.598E+04      -0.824    0.841      109
   -19.0      -0.190E+04   0.553E+04   0.585E+04      -0.849    0.817      110
   -20.0      -0.195E+04   0.537E+04   0.571E+04      -0.871    0.793      111
   -21.0      -0.200E+04   0.520E+04   0.557E+04      -0.890    0.769      112
   -22.0      -0.203E+04   0.503E+04   0.543E+04      -0.907    0.744      113
   -23.0      -0.207E+04   0.487E+04   0.529E+04      -0.921    0.719      114
   -24.0      -0.209E+04   0.470E+04   0.514E+04      -0.933    0.694      115
   -25.0      -0.212E+04   0.454E+04   0.501E+04      -0.943    0.670      116
   -26.0      -0.213E+04   0.438E+04   0.487E+04      -0.952    0.647      117
   -27.0      -0.215E+04   0.422E+04   0.474E+04      -0.959    0.624      118
   -28.0      -0.216E+04   0.407E+04   0.461E+04      -0.965    0.602      119
   -29.0      -0.218E+04   0.393E+04   0.449E+04      -0.970    0.580      120
   -30.0      -0.219E+04   0.379E+04   0.437E+04      -0.974    0.559      121
   -31.0      -0.219E+04   0.365E+04   0.426E+04      -0.978    0.539      122
   -32.0      -0.220E+04   0.352E+04   0.415E+04      -0.981    0.520      123
   -33.0      -0.221E+04   0.340E+04   0.405E+04      -0.984    0.502      124
   -34.0      -0.221E+04   0.328E+04   0.395E+04      -0.986    0.484      125
   -35.0      -0.222E+04   0.316E+04   0.386E+04      -0.988    0.467      126
   -36.0      -0.222E+04   0.305E+04   0.377E+04      -0.989    0.451      127
   -37.0      -0.222E+04   0.295E+04   0.369E+04      -0.991    0.436      128
   -38.0      -0.222E+04   0.285E+04   0.361E+04      -0.992    0.421      129
   -39.0      -0.223E+04   0.275E+04   0.354E+04      -0.993    0.406      130
   -40.0      -0.223E+04   0.266E+04   0.347E+04      -0.994    0.393      131
   -41.0      -0.223E+04   0.257E+04   0.340E+04      -0.995    0.379      132
   -42.0      -0.223E+04   0.248E+04   0.334E+04      -0.995    0.366      133
   -43.0      -0.223E+04   0.239E+04   0.327E+04      -0.996    0.354      134
   -44.0      -0.223E+04   0.231E+04   0.322E+04      -0.996    0.342      135
   -45.0      -0.223E+04   0.223E+04   0.316E+04      -0.996    0.330      136
   -46.0      -0.224E+04   0.216E+04   0.311E+04      -0.997    0.319      137
   -47.0      -0.224E+04   0.209E+04   0.306E+04      -0.997    0.308      138
   -48.0      -0.224E+04   0.201E+04   0.301E+04      -0.997    0.298      139
   -49.0      -0.224E+04   0.195E+04   0.296E+04      -0.998    0.287      140
   -50.0      -0.224E+04   0.188E+04   0.292E+04      -0.998    0.278      141
   -51.0      -0.224E+04   0.181E+04   0.288E+04      -0.998    0.268      142
   -52.0      -0.224E+04   0.175E+04   0.284E+04      -0.998    0.258      143
   -53.0      -0.224E+04   0.169E+04   0.280E+04      -0.998    0.249      144
   -54.0      -0.224E+04   0.163E+04   0.277E+04      -0.999    0.240      145
   -55.0      -0.224E+04   0.157E+04   0.273E+04      -0.999    0.232      146
   -56.0      -0.224E+04   0.151E+04   0.270E+04      -0.998    0.223      147
   -57.0      -0.224E+04   0.145E+04   0.267E+04      -0.999    0.215      148
   -58.0      -0.224E+04   0.140E+04   0.264E+04      -0.999    0.207      149
   -59.0      -0.224E+04   0.135E+04   0.261E+04      -0.999    0.199      150
   -60.0      -0.224E+04   0.129E+04   0.259E+04      -0.999    0.191      151
   -61.0      -0.224E+04   0.124E+04   0.256E+04      -0.999    0.184      152
   -62.0      -0.224E+04   0.119E+04   0.254E+04      -0.999    0.176      153
   -63.0      -0.224E+04   0.114E+04   0.252E+04      -0.999    0.169      154
   -64.0      -0.224E+04   0.109E+04   0.249E+04      -0.999    0.162      155
   -65.0      -0.224E+04   0.105E+04   0.247E+04      -0.999    0.154      156
   -66.0      -0.224E+04    998.       0.245E+04      -1.000    0.148      157
   -67.0      -0.224E+04    951.       0.244E+04      -0.999    0.141      158
   -68.0      -0.224E+04    906.       0.242E+04      -1.000    0.134      159
   -69.0      -0.224E+04    860.       0.240E+04      -0.999    0.127      160
   -70.0      -0.224E+04    816.       0.239E+04      -0.999    0.121      161
   -71.0      -0.224E+04    772.       0.237E+04      -1.000    0.114      162
   -72.0      -0.224E+04    728.       0.236E+04      -0.999    0.108      163
   -73.0      -0.224E+04    685.       0.234E+04      -1.000    0.101      164
   -74.0      -0.224E+04    643.       0.233E+04      -1.000    0.095      165
   -75.0      -0.224E+04    601.       0.232E+04      -0.999    0.089      166
   -76.0      -0.224E+04    559.       0.231E+04      -0.999    0.083      167
   -77.0      -0.224E+04    518.       0.230E+04      -1.000    0.077      168
   -78.0      -0.224E+04    476.       0.229E+04      -0.999    0.070      169
   -79.0      -0.224E+04    436.       0.228E+04      -1.000    0.064      170
   -80.0      -0.224E+04    395.       0.228E+04      -1.000    0.058      171
   -81.0      -0.224E+04    355.       0.227E+04      -1.000    0.052      172
   -82.0      -0.224E+04    315.       0.226E+04      -1.000    0.047      173
   -83.0      -0.224E+04    275.       0.226E+04      -1.000    0.041      174
   -84.0      -0.224E+04    236.       0.225E+04      -0.999    0.035      175
   -85.0      -0.224E+04    196.       0.225E+04      -1.000    0.029      176
   -86.0      -0.224E+04    157.       0.225E+04      -1.000    0.023      177
   -87.0      -0.224E+04    117.       0.224E+04      -1.000    0.017      178
   -88.0      -0.224E+04    78.3       0.224E+04      -0.999    0.012      179
   -89.0      -0.224E+04    39.1       0.224E+04      -1.000    0.006      180
   -90.0      -0.224E+04   0.137E-12   0.224E+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
