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


 datain_caisson: < --- BEGIN

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

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

 pidata: END --- >


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


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


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


 ssidat: < --- BEGIN
 alpha_Fo  alpha_Fi  eta_Fo
   0.62   0.54   1.00
 eta_reb  alpha_d_su  iplug_vmode  vmax_user
   0.50   0.50    1   0.00
 eta_amh  eta_bmh  eta_phimh
   1.00   1.00   1.00
 ssidat: END --- >


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


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


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

 datain_caisson: END --- >


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



caisson_vhm: nlines_warnings :    5
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. 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]:    6935.631
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000       0.000
   zlug_cL                               [m]:      10.000
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6935.631
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    -548.032    6935.631  -69356.307
   Vstar,Hstar,Mstar                     [-]:      -0.222      10.988      -7.325
   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.988      -7.325      -0.222
   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.995
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.001
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.996
 > f_VHM_fun                             [-]: -0.3010790E-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                           [-]:  1
 < L,D                                   [m]:  15.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6912.192
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      10.000
   zlug_cL                               [m]:      10.397
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1200.290
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6807.181
   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]:     652.258    6807.181  -70772.457
   Vstar,Hstar,Mstar                     [-]:       0.265      10.784      -7.475
   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.784      -7.475       0.265
   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.982
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.011
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.994
 > f_VHM_fun                             [-]: -0.7230851E-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                           [-]:  1
 < L,D                                   [m]:  15.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6480.763
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      20.000
   zlug_cL                               [m]:      10.819
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2216.551
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6089.925
   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]:    1668.520    6089.925  -65886.493
   Vstar,Hstar,Mstar                     [-]:       0.677       9.648      -6.959
   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.648      -6.959       0.677
   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.807
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.053
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.859
 > f_VHM_fun                             [-]:  0.2857876E-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                           [-]:  1
 < L,D                                   [m]:  15.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5483.152
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      30.000
   zlug_cL                               [m]:      11.299
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2741.576
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4748.549
   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]:    2193.544    4748.549  -53654.034
   Vstar,Hstar,Mstar                     [-]:       0.890       7.523      -5.667
   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.523      -5.667       0.890
   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.505
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.089
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.594
 > f_VHM_fun                             [-]: -0.5169329E-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                           [-]:  1
 < L,D                                   [m]:  15.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    4538.504
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      40.000
   zlug_cL                               [m]:      11.888
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2917.294
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3476.696
   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]:    2369.262    3476.696  -41330.867
   Vstar,Hstar,Mstar                     [-]:       0.961       5.508      -4.365
   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.508      -4.365       0.961
   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.280
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.106
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.386
 > f_VHM_fun                             [-]: -0.6098825E-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                           [-]:  1
 < L,D                                   [m]:  15.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3880.550
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      50.000
   zlug_cL                               [m]:      12.681
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2972.674
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2494.369
   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]:    2424.642    2494.369  -31632.210
   Vstar,Hstar,Mstar                     [-]:       0.984       3.952      -3.341
   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.952      -3.341       0.984
   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.151
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.114
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.264
 > f_VHM_fun                             [-]:  0.1410586E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00014
root_Tlug_22_fun: END --- >


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

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3455.229
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      60.000
   zlug_cL                               [m]:      13.897
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2992.316
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1727.615
   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]:    2444.284    1727.615  -24008.857
   Vstar,Hstar,Mstar                     [-]:       0.992       2.737      -2.536
   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.737      -2.536       0.992
   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.078
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.118
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.196
 > f_VHM_fun                             [-]: -0.1460455E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00015
root_Tlug_22_fun: END --- >


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

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3192.775
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      70.000
   zlug_cL                               [m]:      16.182
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3000.227
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1091.993
   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]:    2452.195    1091.993  -17670.444
   Vstar,Hstar,Mstar                     [-]:       0.995       1.730      -1.866
   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.730      -1.866       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.035
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.121
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.156
 > f_VHM_fun                             [-]: -0.4807659E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00048
root_Tlug_22_fun: END --- >


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

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3049.754
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      80.000
   zlug_cL                               [m]:      22.760
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3003.421
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     529.584
   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]:    2455.389     529.584  -12053.540
   Vstar,Hstar,Mstar                     [-]:       0.996       0.839      -1.273
   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.839      -1.273       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.123
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.136
 > f_VHM_fun                             [-]: -0.1824628E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00182
root_Tlug_22_fun: END --- >


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

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3004.829
   L , D                                 [m]:      15.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.000      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     548.032       0.000
   su_av_L                             [kPa]:       9.352
   Vmax                                 [kN]:    2464.134
   Hfac=L   D su_av_L                   [kN]:     631.227
   Mfac=L^2 D su_av_L                  [kNm]:    9468.400
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3004.829
   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]:    2456.797       0.000   -6760.864
   Vstar,Hstar,Mstar                     [-]:       0.997       0.000      -0.714
   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.714       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.125
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.126
 > f_VHM_fun                             [-]: -0.3374859E-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                           [-]:  1
 < L,D                                   [m]:  15.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    7845.124
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500       0.000
   zlug_cL                               [m]:      10.500
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    7845.124
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    -581.019    7845.124  -82373.800
   Vstar,Hstar,Mstar                     [-]:      -0.217      11.108      -7.290
   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.108      -7.290      -0.217
   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.987
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.009
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.996
 > f_VHM_fun                             [-]: -0.2972339E-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                           [-]:  1
 < L,D                                   [m]:  16.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    7881.385
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      10.000
   zlug_cL                               [m]:      10.897
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1368.588
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    7761.649
   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]:     787.569    7761.649  -84576.642
   Vstar,Hstar,Mstar                     [-]:       0.294      10.990      -7.485
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  16.00
   X = Hstar; Y = Mstar; Z = Vstar                 10.990      -7.485       0.294
   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.989
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.002
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.991
 > f_VHM_fun                             [-]: -0.8448694E-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                           [-]:  1
 < L,D                                   [m]:  16.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    7337.824
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      20.000
   zlug_cL                               [m]:      11.319
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2509.684
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6895.299
   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]:    1928.665    6895.299  -78047.432
   Vstar,Hstar,Mstar                     [-]:       0.720       9.763      -6.907
   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.763      -6.907       0.720
   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.800
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.027
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.826
 > f_VHM_fun                             [-]:  0.4281164E-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                           [-]:  1
 < L,D                                   [m]:  16.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6081.270
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      30.000
   zlug_cL                               [m]:      11.799
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3040.635
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5266.535
   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]:    2459.616    5266.535  -62140.042
   Vstar,Hstar,Mstar                     [-]:       0.919       7.457      -5.499
   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.457      -5.499       0.919
   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.479
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.055
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.534
 > f_VHM_fun                             [-]: -0.8520862E-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                           [-]:  1
 < L,D                                   [m]:  16.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    4965.184
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      40.000
   zlug_cL                               [m]:      12.388
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3191.559
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3803.552
   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]:    2610.540    3803.552  -47118.298
   Vstar,Hstar,Mstar                     [-]:       0.975       5.386      -4.170
   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.386      -4.170       0.975
   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.258
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.071
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.330
 > f_VHM_fun                             [-]: -0.8362746E-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                           [-]:  1
 < L,D                                   [m]:  16.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    4218.960
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      50.000
   zlug_cL                               [m]:      13.181
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3231.911
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2711.895
   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]:    2650.892    2711.895  -35746.700
   Vstar,Hstar,Mstar                     [-]:       0.990       3.840      -3.163
   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.840      -3.163       0.990
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.137
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.081
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.222
 > f_VHM_fun                             [-]: -0.4317695E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00432
root_Tlug_22_fun: END --- >


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

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3748.643
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      60.000
   zlug_cL                               [m]:      14.397
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3246.420
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1874.321
   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]:    2665.401    1874.321  -26984.818
   Vstar,Hstar,Mstar                     [-]:       0.996       2.654      -2.388
   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.654      -2.388       0.996
   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.070
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.088
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.158
 > f_VHM_fun                             [-]:  0.4438164E-04
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00004
root_Tlug_22_fun: END --- >


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

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3459.766
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      70.000
   zlug_cL                               [m]:      16.682
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3251.117
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1183.310
   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]:    2670.098    1183.310  -19739.766
   Vstar,Hstar,Mstar                     [-]:       0.997       1.676      -1.747
   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.676      -1.747       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.031
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.093
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.126
 > f_VHM_fun                             [-]: -0.2021890E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00202
root_Tlug_22_fun: END --- >


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

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3303.470
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      80.000
   zlug_cL                               [m]:      23.260
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3253.283
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     573.642
   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]:    2672.264     573.642  -13343.123
   Vstar,Hstar,Mstar                     [-]:       0.998       0.812      -1.181
   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.812      -1.181       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.010
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.098
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.108
 > f_VHM_fun                             [-]: -0.7646162E-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                           [-]:  1
 < L,D                                   [m]:  16.00   4.50

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

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


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3253.765
   L , D                                 [m]:      16.000       4.500
   zlug , theta_lug_star           [m] [deg]:      10.500      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     581.019       0.000
   su_av_L                             [kPa]:       9.809
   Vmax                                 [kN]:    2677.051
   Hfac=L   D su_av_L                   [kN]:     706.237
   Mfac=L^2 D su_av_L                  [kNm]:   11299.787
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    3253.765
   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]:    2672.746       0.000   -7320.971
   Vstar,Hstar,Mstar                     [-]:       0.998       0.000      -0.648
   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.648       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.102
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.103
 > f_VHM_fun                             [-]: -0.4291057E-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                           [-]:  1
 < L,D                                   [m]:  16.00   4.50

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

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


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  16.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  1
 > 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  =  1
 ISCOPE_BW, caisson weight BW mode =  1
 ilc    = VHM load case #            
 icgeom = caisson geometry #         
 D    = caisson outer diameter       
 L    = caisson embedded length      
 L/D  = caisson embedment ratio      
 Vmax = caisson axial resistance     
 Hmax = caisson FIXED HEAD lateral resistance 
 H0   = caisson FREE  HEAD lateral resistance 
 Bw   = caisson submerged weight (only used for anchor foundations)
 Bw*  = caisson+soil plug submerged weight (only for IPLUG_VMODE=5)
 su_av_L = soil average su between caisson head (seafloor) and tip (L)
 ----------------------------------------------------------------------------------------
  ilc       D      L      L/D       Vmax        Hmax         H0          Bw       su_av_L
   icgeom  [m]    [m]     [-]       [kN]        [kN]        [kN]        [kN]        [kPa]
 ----------------------------------------------------------------------------------------
    1  1   4.50  15.00   3.33      2464.1      6931.9      1893.3       548.0         9.4
    1  2   4.50  16.00   3.56      2677.1      7831.9      2125.6       581.0         9.8

    2  1   4.50  15.00   3.33      2464.1      6931.9      1893.3       548.0         9.4
    2  2   4.50  16.00   3.56      2677.1      7831.9      2125.6       581.0         9.8

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

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

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

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

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

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

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

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


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

    2  1   4.50  15.00  10.00 10.00   6912.2   548.0       5  -0.723E-05   ! ? missing text
    2  2   4.50  16.00  10.50 10.00   7881.4   581.0       5  -0.845E-05   ! ? missing text

    3  1   4.50  15.00  10.00 20.00   6480.8   548.0       6   0.286E-06   ! ? missing text
    3  2   4.50  16.00  10.50 20.00   7337.8   581.0       6   0.428E-06   ! ? missing text

    4  1   4.50  15.00  10.00 30.00   5483.2   548.0       5  -0.517E-04   ! ? missing text
    4  2   4.50  16.00  10.50 30.00   6081.3   581.0       5  -0.852E-05   ! ? missing text

    5  1   4.50  15.00  10.00 40.00   4538.5   548.0       4  -0.610E-05   ! ? missing text
    5  2   4.50  16.00  10.50 40.00   4965.2   581.0       9  -0.836E-05   ! ? missing text

    6  1   4.50  15.00  10.00 50.00   3880.5   548.0       7   0.141E-03   ! ? missing text
    6  2   4.50  16.00  10.50 50.00   4219.0   581.0       6  -0.432E-02   ! ? missing text

    7  1   4.50  15.00  10.00 60.00   3455.2   548.0       9  -0.146E-03   ! ? missing text
    7  2   4.50  16.00  10.50 60.00   3748.6   581.0       8   0.444E-04   ! ? missing text

    8  1   4.50  15.00  10.00 70.00   3192.8   548.0       9  -0.481E-03   ! ? missing text
    8  2   4.50  16.00  10.50 70.00   3459.8   581.0       8  -0.202E-02   ! ? missing text

    9  1   4.50  15.00  10.00 80.00   3049.8   548.0      11  -0.182E-02   ! ? missing text
    9  2   4.50  16.00  10.50 80.00   3303.5   581.0      24  -0.765E-05   ! ? missing text

   10  1   4.50  15.00  10.00 90.00   3004.8   548.0      24  -0.337E-05   ! ? missing text
   10  2   4.50  16.00  10.50 90.00   3253.8   581.0      24  -0.429E-05   ! ? missing text
 -----------------------------------------------------------------------
 ilc        D      L    zlug  theta    Tlug     Bw    ( iters    f_VHM )
  icgeom   [m]    [m]   [m]   [deg]    [kN]    [kN]   ( [-]      [-]   )
 -----------------------------------------------------------------------

VHlug_resist_sub: < --- BEGIN
  load case #=    1; caisson geometry #=    1 (#rows results= =  181)
  L         [m] =   15.00
  D         [m] =    4.50
  zlug      [m] =   10.00
  iro_lug   [-] =   1: offset + D/2 (lug on outer face)
  m_fac     [-] =   1.000
  Hmax     [kN] =    0.693E+04
  Vmax     [kN] =    0.246E+04
  theta*        = resistance inclination at lug
  Vlug_theta*, Hlug_theta*
                = corresponding VH UNFACTORED resistances
  T_theta*      = sqrt( Vlug_theta*^2 + Hlug_theta*^2)
------------------------------------------------------------------------------
 theta*     Vlug_theta*  Hlug_theta*  T_theta*   Vlug_theta*  Hlug_theta*  row
  [deg]        [kN]         [kN]        [kN]      / Vmax [-]   / Hmax [-]  [-]
------------------------------------------------------------------------------
    90.0       0.246E+04   0.151E-12   0.246E+04       0.998    0.000        1
    89.0       0.246E+04    43.0       0.246E+04       0.999    0.006        2
    88.0       0.246E+04    85.9       0.246E+04       0.999    0.012        3
    87.0       0.246E+04    129.       0.246E+04       0.999    0.019        4
    86.0       0.246E+04    172.       0.247E+04       0.999    0.025        5
    85.0       0.246E+04    215.       0.247E+04       0.999    0.031        6
    84.0       0.246E+04    259.       0.247E+04       0.999    0.037        7
    83.0       0.246E+04    302.       0.248E+04       0.999    0.044        8
    82.0       0.246E+04    346.       0.249E+04       0.999    0.050        9
    81.0       0.246E+04    390.       0.249E+04       0.999    0.056       10
    80.0       0.246E+04    434.       0.250E+04       0.999    0.063       11
    79.0       0.246E+04    478.       0.251E+04       0.999    0.069       12
    78.0       0.246E+04    523.       0.252E+04       0.999    0.075       13
    77.0       0.246E+04    568.       0.253E+04       0.998    0.082       14
    76.0       0.246E+04    613.       0.254E+04       0.998    0.088       15
    75.0       0.246E+04    659.       0.255E+04       0.998    0.095       16
    74.0       0.246E+04    705.       0.256E+04       0.998    0.102       17
    73.0       0.246E+04    752.       0.257E+04       0.998    0.108       18
    72.0       0.246E+04    799.       0.259E+04       0.998    0.115       19
    71.0       0.246E+04    847.       0.260E+04       0.998    0.122       20
    70.0       0.246E+04    895.       0.262E+04       0.998    0.129       21
    69.0       0.246E+04    943.       0.263E+04       0.997    0.136       22
    68.0       0.246E+04    993.       0.265E+04       0.998    0.143       23
    67.0       0.246E+04   0.104E+04   0.267E+04       0.997    0.151       24
    66.0       0.246E+04   0.109E+04   0.269E+04       0.998    0.158       25
    65.0       0.246E+04   0.115E+04   0.271E+04       0.997    0.165       26
    64.0       0.246E+04   0.120E+04   0.273E+04       0.997    0.173       27
    63.0       0.246E+04   0.125E+04   0.276E+04       0.997    0.181       28
    62.0       0.246E+04   0.131E+04   0.278E+04       0.997    0.188       29
    61.0       0.246E+04   0.136E+04   0.281E+04       0.996    0.196       30
    60.0       0.246E+04   0.142E+04   0.284E+04       0.997    0.205       31
    59.0       0.246E+04   0.148E+04   0.286E+04       0.996    0.213       32
    58.0       0.245E+04   0.153E+04   0.289E+04       0.996    0.221       33
    57.0       0.245E+04   0.159E+04   0.293E+04       0.996    0.230       34
    56.0       0.245E+04   0.165E+04   0.296E+04       0.996    0.239       35
    55.0       0.245E+04   0.172E+04   0.299E+04       0.995    0.248       36
    54.0       0.245E+04   0.178E+04   0.303E+04       0.995    0.257       37
    53.0       0.245E+04   0.185E+04   0.307E+04       0.995    0.266       38
    52.0       0.245E+04   0.191E+04   0.311E+04       0.994    0.276       39
    51.0       0.245E+04   0.198E+04   0.315E+04       0.994    0.286       40
    50.0       0.245E+04   0.205E+04   0.320E+04       0.993    0.296       41
    49.0       0.245E+04   0.213E+04   0.324E+04       0.993    0.307       42
    48.0       0.244E+04   0.220E+04   0.329E+04       0.992    0.318       43
    47.0       0.244E+04   0.228E+04   0.334E+04       0.991    0.329       44
    46.0       0.244E+04   0.236E+04   0.339E+04       0.991    0.340       45
    45.0       0.244E+04   0.244E+04   0.345E+04       0.990    0.352       46
    44.0       0.244E+04   0.252E+04   0.351E+04       0.989    0.364       47
    43.0       0.243E+04   0.261E+04   0.357E+04       0.988    0.377       48
    42.0       0.243E+04   0.270E+04   0.363E+04       0.987    0.390       49
    41.0       0.243E+04   0.279E+04   0.370E+04       0.985    0.403       50
    40.0       0.242E+04   0.289E+04   0.377E+04       0.984    0.417       51
    39.0       0.242E+04   0.299E+04   0.384E+04       0.982    0.431       52
    38.0       0.241E+04   0.309E+04   0.392E+04       0.980    0.446       53
    37.0       0.241E+04   0.320E+04   0.400E+04       0.978    0.461       54
    36.0       0.240E+04   0.331E+04   0.409E+04       0.975    0.477       55
    35.0       0.240E+04   0.342E+04   0.418E+04       0.972    0.494       56
    34.0       0.239E+04   0.354E+04   0.427E+04       0.969    0.511       57
    33.0       0.238E+04   0.366E+04   0.437E+04       0.965    0.528       58
    32.0       0.237E+04   0.379E+04   0.447E+04       0.960    0.546       59
    31.0       0.235E+04   0.392E+04   0.457E+04       0.955    0.565       60
    30.0       0.234E+04   0.405E+04   0.468E+04       0.950    0.585       61
    29.0       0.232E+04   0.419E+04   0.479E+04       0.943    0.605       62
    28.0       0.230E+04   0.433E+04   0.491E+04       0.935    0.625       63
    27.0       0.228E+04   0.448E+04   0.503E+04       0.927    0.646       64
    26.0       0.226E+04   0.463E+04   0.515E+04       0.916    0.668       65
    25.0       0.223E+04   0.478E+04   0.528E+04       0.905    0.690       66
    24.0       0.220E+04   0.494E+04   0.540E+04       0.892    0.712       67
    23.0       0.216E+04   0.509E+04   0.553E+04       0.877    0.734       68
    22.0       0.212E+04   0.525E+04   0.566E+04       0.860    0.757       69
    21.0       0.207E+04   0.540E+04   0.578E+04       0.841    0.779       70
    20.0       0.202E+04   0.555E+04   0.591E+04       0.820    0.801       71
    19.0       0.196E+04   0.570E+04   0.603E+04       0.796    0.822       72
    18.0       0.190E+04   0.584E+04   0.614E+04       0.770    0.843       73
    17.0       0.183E+04   0.598E+04   0.625E+04       0.742    0.862       74
    16.0       0.175E+04   0.611E+04   0.635E+04       0.711    0.881       75
    15.0       0.167E+04   0.623E+04   0.645E+04       0.677    0.898       76
    14.0       0.158E+04   0.634E+04   0.653E+04       0.641    0.914       77
    13.0       0.149E+04   0.644E+04   0.661E+04       0.603    0.929       78
    12.0       0.139E+04   0.653E+04   0.668E+04       0.563    0.942       79
    11.0       0.128E+04   0.661E+04   0.673E+04       0.521    0.953       80
    10.0       0.118E+04   0.668E+04   0.678E+04       0.478    0.964       81
    9.00       0.107E+04   0.674E+04   0.682E+04       0.433    0.972       82
    8.00        954.       0.679E+04   0.686E+04       0.387    0.979       83
    7.00        839.       0.683E+04   0.688E+04       0.340    0.985       84
    6.00        721.       0.686E+04   0.690E+04       0.293    0.990       85
    5.00        603.       0.689E+04   0.692E+04       0.245    0.994       86
    4.00        483.       0.691E+04   0.693E+04       0.196    0.997       87
    3.00        363.       0.692E+04   0.693E+04       0.147    0.999       88
    2.00        242.       0.693E+04   0.694E+04       0.098    1.000       89
    1.00        121.       0.694E+04   0.694E+04       0.049    1.002       90
    0.00        0.00       0.695E+04   0.695E+04       0.000    1.002       91
   -1.00       -121.       0.695E+04   0.695E+04      -0.049    1.003       92
   -2.00       -243.       0.695E+04   0.696E+04      -0.099    1.003       93
   -3.00       -364.       0.695E+04   0.696E+04      -0.148    1.003       94
   -4.00       -486.       0.695E+04   0.696E+04      -0.197    1.002       95
   -5.00       -607.       0.694E+04   0.696E+04      -0.246    1.001       96
   -6.00       -727.       0.692E+04   0.696E+04      -0.295    0.998       97
   -7.00       -847.       0.690E+04   0.695E+04      -0.344    0.995       98
   -8.00       -964.       0.686E+04   0.693E+04      -0.391    0.990       99
   -9.00      -0.108E+04   0.682E+04   0.690E+04      -0.438    0.984      100
   -10.0      -0.119E+04   0.677E+04   0.687E+04      -0.484    0.976      101
   -11.0      -0.130E+04   0.670E+04   0.683E+04      -0.529    0.967      102
   -12.0      -0.141E+04   0.663E+04   0.677E+04      -0.572    0.956      103
   -13.0      -0.151E+04   0.654E+04   0.671E+04      -0.613    0.943      104
   -14.0      -0.161E+04   0.644E+04   0.664E+04      -0.651    0.929      105
   -15.0      -0.170E+04   0.633E+04   0.655E+04      -0.688    0.913      106
   -16.0      -0.178E+04   0.620E+04   0.646E+04      -0.722    0.895      107
   -17.0      -0.186E+04   0.607E+04   0.635E+04      -0.754    0.876      108
   -18.0      -0.193E+04   0.593E+04   0.624E+04      -0.782    0.856      109
   -19.0      -0.199E+04   0.579E+04   0.612E+04      -0.808    0.835      110
   -20.0      -0.205E+04   0.563E+04   0.599E+04      -0.832    0.813      111
   -21.0      -0.210E+04   0.548E+04   0.587E+04      -0.853    0.790      112
   -22.0      -0.215E+04   0.532E+04   0.573E+04      -0.872    0.767      113
   -23.0      -0.219E+04   0.515E+04   0.560E+04      -0.888    0.744      114
   -24.0      -0.222E+04   0.499E+04   0.547E+04      -0.902    0.720      115
   -25.0      -0.225E+04   0.483E+04   0.533E+04      -0.915    0.697      116
   -26.0      -0.228E+04   0.468E+04   0.520E+04      -0.926    0.675      117
   -27.0      -0.230E+04   0.452E+04   0.508E+04      -0.935    0.652      118
   -28.0      -0.232E+04   0.437E+04   0.495E+04      -0.943    0.631      119
   -29.0      -0.234E+04   0.423E+04   0.483E+04      -0.950    0.609      120
   -30.0      -0.236E+04   0.408E+04   0.471E+04      -0.957    0.589      121
   -31.0      -0.237E+04   0.394E+04   0.460E+04      -0.962    0.569      122
   -32.0      -0.238E+04   0.381E+04   0.449E+04      -0.966    0.550      123
   -33.0      -0.239E+04   0.368E+04   0.439E+04      -0.970    0.531      124
   -34.0      -0.240E+04   0.356E+04   0.429E+04      -0.974    0.513      125
   -35.0      -0.241E+04   0.344E+04   0.420E+04      -0.977    0.496      126
   -36.0      -0.241E+04   0.332E+04   0.411E+04      -0.979    0.479      127
   -37.0      -0.242E+04   0.321E+04   0.402E+04      -0.981    0.463      128
   -38.0      -0.242E+04   0.310E+04   0.394E+04      -0.983    0.447      129
   -39.0      -0.243E+04   0.300E+04   0.386E+04      -0.985    0.432      130
   -40.0      -0.243E+04   0.290E+04   0.378E+04      -0.987    0.418      131
   -41.0      -0.243E+04   0.280E+04   0.371E+04      -0.988    0.404      132
   -42.0      -0.244E+04   0.271E+04   0.364E+04      -0.989    0.390      133
   -43.0      -0.244E+04   0.262E+04   0.358E+04      -0.990    0.377      134
   -44.0      -0.244E+04   0.253E+04   0.352E+04      -0.991    0.365      135
   -45.0      -0.244E+04   0.244E+04   0.346E+04      -0.992    0.353      136
   -46.0      -0.244E+04   0.236E+04   0.340E+04      -0.992    0.341      137
   -47.0      -0.245E+04   0.228E+04   0.335E+04      -0.993    0.329      138
   -48.0      -0.245E+04   0.220E+04   0.329E+04      -0.994    0.318      139
   -49.0      -0.245E+04   0.213E+04   0.325E+04      -0.994    0.307      140
   -50.0      -0.245E+04   0.206E+04   0.320E+04      -0.995    0.297      141
   -51.0      -0.245E+04   0.199E+04   0.316E+04      -0.995    0.286      142
   -52.0      -0.245E+04   0.192E+04   0.311E+04      -0.996    0.276      143
   -53.0      -0.245E+04   0.185E+04   0.307E+04      -0.996    0.267      144
   -54.0      -0.245E+04   0.178E+04   0.303E+04      -0.996    0.257      145
   -55.0      -0.246E+04   0.172E+04   0.300E+04      -0.996    0.248      146
   -56.0      -0.246E+04   0.166E+04   0.296E+04      -0.997    0.239      147
   -57.0      -0.246E+04   0.160E+04   0.293E+04      -0.997    0.230      148
   -58.0      -0.246E+04   0.154E+04   0.290E+04      -0.997    0.221      149
   -59.0      -0.246E+04   0.148E+04   0.287E+04      -0.997    0.213      150
   -60.0      -0.246E+04   0.142E+04   0.284E+04      -0.997    0.205      151
   -61.0      -0.246E+04   0.136E+04   0.281E+04      -0.997    0.197      152
   -62.0      -0.246E+04   0.131E+04   0.278E+04      -0.998    0.189      153
   -63.0      -0.246E+04   0.125E+04   0.276E+04      -0.998    0.181      154
   -64.0      -0.246E+04   0.120E+04   0.274E+04      -0.998    0.173      155
   -65.0      -0.246E+04   0.115E+04   0.271E+04      -0.998    0.165      156
   -66.0      -0.246E+04   0.109E+04   0.269E+04      -0.997    0.158      157
   -67.0      -0.246E+04   0.104E+04   0.267E+04      -0.998    0.151      158
   -68.0      -0.246E+04    994.       0.265E+04      -0.998    0.143      159
   -69.0      -0.246E+04    944.       0.263E+04      -0.998    0.136      160
   -70.0      -0.246E+04    895.       0.262E+04      -0.998    0.129      161
   -71.0      -0.246E+04    847.       0.260E+04      -0.998    0.122      162
   -72.0      -0.246E+04    799.       0.259E+04      -0.998    0.115      163
   -73.0      -0.246E+04    752.       0.257E+04      -0.998    0.109      164
   -74.0      -0.246E+04    706.       0.256E+04      -0.999    0.102      165
   -75.0      -0.246E+04    659.       0.255E+04      -0.999    0.095      166
   -76.0      -0.246E+04    614.       0.254E+04      -0.999    0.089      167
   -77.0      -0.246E+04    568.       0.253E+04      -0.999    0.082      168
   -78.0      -0.246E+04    523.       0.252E+04      -0.999    0.075      169
   -79.0      -0.246E+04    478.       0.251E+04      -0.999    0.069      170
   -80.0      -0.246E+04    434.       0.250E+04      -0.999    0.063      171
   -81.0      -0.246E+04    390.       0.249E+04      -0.999    0.056      172
   -82.0      -0.246E+04    346.       0.249E+04      -0.999    0.050      173
   -83.0      -0.246E+04    302.       0.248E+04      -0.998    0.044      174
   -84.0      -0.246E+04    259.       0.248E+04      -0.999    0.037      175
   -85.0      -0.246E+04    215.       0.247E+04      -0.999    0.031      176
   -86.0      -0.246E+04    172.       0.247E+04      -0.999    0.025      177
   -87.0      -0.246E+04    129.       0.246E+04      -0.999    0.019      178
   -88.0      -0.246E+04    85.9       0.246E+04      -0.998    0.012      179
   -89.0      -0.246E+04    42.9       0.246E+04      -0.998    0.006      180
   -90.0      -0.246E+04   0.151E-12   0.246E+04      -0.998    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
