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


 datain_caisson: < --- BEGIN

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

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

 pidata: END --- >


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


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


 cfdata: < --- BEGIN
 ncgeom
    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.63   0.54   1.00
 eta_reb  alpha_d_su  iplug_vmode  vmax_user
   1.00   0.50    2   0.00
 eta_amh  eta_bmh  eta_phimh
   1.00   1.00   1.00
 ssidat: END --- >


 chndat: < --- BEGIN
 iscope_Bw  l_fac_Bw  eta_Vst_wall  delta_Bw_user  Bw_caisson_user
    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]:    5316.743
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5316.743
   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    5316.743  -53167.434
   Vstar,Hstar,Mstar                     [-]:      -0.308      10.950      -7.300
   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.950      -7.300      -0.308
   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.988
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.001
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.989
 > f_VHM_fun                             [-]: -0.2966609E-05
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5325.897
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:     924.832
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5244.984
   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]:     376.800    5244.984  -54530.715
   Vstar,Hstar,Mstar                     [-]:       0.212      10.802      -7.487
   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.802      -7.487       0.212
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.985
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.011
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.997
 > f_VHM_fun                             [-]: -0.6897169E-05
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5023.671
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1718.197
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4720.707
   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]:    1170.165    4720.707  -51073.011
   Vstar,Hstar,Mstar                     [-]:       0.658       9.722      -7.012
   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.722      -7.012       0.658
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.819
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.054
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.873
 > f_VHM_fun                             [-]:  0.5425149E-06
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    4247.249
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2123.624
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3678.225
   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]:    1575.592    3678.225  -41560.406
   Vstar,Hstar,Mstar                     [-]:       0.886       7.575      -5.706
   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.575      -5.706       0.886
   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.512
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.090
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.602
 > f_VHM_fun                             [-]: -0.1045477E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00010
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3509.414
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2255.808
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2688.367
   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]:    1707.776    2688.367  -31959.241
   Vstar,Hstar,Mstar                     [-]:       0.961       5.537      -4.388
   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.537      -4.388       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.283
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.107
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.390
 > f_VHM_fun                             [-]: -0.1747379E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00017
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2997.794
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2296.443
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1926.945
   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]:    1748.412    1926.945  -24436.447
   Vstar,Hstar,Mstar                     [-]:       0.983       3.969      -3.355
   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.969      -3.355       0.983
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.238       1.675     -33.949
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.238       1.675     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.152
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.115
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.269
 > f_VHM_fun                             [-]: -0.1981601E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00198
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2669.014
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2311.434
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1334.507
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1763.402    1334.507  -18545.794
   Vstar,Hstar,Mstar                     [-]:       0.992       2.748      -2.546
   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.748      -2.546       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.119
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.197
 > f_VHM_fun                             [-]: -0.1681460E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00017
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2465.686
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2316.987
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     843.314
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1768.955     843.314  -13646.363
   Vstar,Hstar,Mstar                     [-]:       0.995       1.737      -1.874
   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.737      -1.874       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.122
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.160
 > f_VHM_fun                             [-]: -0.2670728E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00267
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2355.456
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2319.672
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     409.021
   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]:    1771.640     409.021   -9309.468
   Vstar,Hstar,Mstar                     [-]:       0.996       0.842      -1.278
   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.842      -1.278       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.124
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.137
 > f_VHM_fun                             [-]: -0.1241160E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00124
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2320.530
   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]:       7.193
   Vmax                                 [kN]:    1778.003
   Hfac=L   D su_av_L                   [kN]:     485.559
   Mfac=L^2 D su_av_L                  [kNm]:    7283.385
   aMH , bMH                             [-]:      13.238       1.675
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.949
   aVH , bVH                             [-]:       4.501       3.182
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2320.530
   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]:    1772.498       0.000   -5221.192
   Vstar,Hstar,Mstar                     [-]:       0.997       0.000      -0.717
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   4.50  15.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.000      -0.717       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.126
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.128
 > f_VHM_fun                             [-]: -0.1200863E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00120
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.19    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      967.12170
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      810.88121
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.33    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    13.3    
   qv_tip = Nc * su_av_tip             [kPa]     119.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.190E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2866.56639
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1778.00291
   Qleaking   = Fo                      [kN]:      967.12170
 > Vmax_fun                             [kN]:     1778.00291
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6017.733
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6017.733
   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    6017.733  -63186.197
   Vstar,Hstar,Mstar                     [-]:      -0.292      11.077      -7.269
   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.077      -7.269      -0.292
   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.982
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.009
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.991
 > f_VHM_fun                             [-]: -0.2950586E-05
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6076.572
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1055.186
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5984.255
   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]:     474.167    5984.255  -65208.847
   Vstar,Hstar,Mstar                     [-]:       0.238      11.015      -7.502
   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.015      -7.502       0.238
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.993
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.002
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.995
 > f_VHM_fun                             [-]: -0.7509203E-05
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    5723.031
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1957.392
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5377.890
   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]:    1376.373    5377.890  -60871.979
   Vstar,Hstar,Mstar                     [-]:       0.692       9.899      -7.003
   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.899      -7.003       0.692
   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.822
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.027
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.849
 > f_VHM_fun                             [-]:  0.6815793E-06
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    4779.633
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2389.817
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4139.284
   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]:    1808.798    4139.284  -48839.569
   Vstar,Hstar,Mstar                     [-]:       0.909       7.619      -5.619
   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.619      -5.619       0.909
   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.500
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.058
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.558
 > f_VHM_fun                             [-]: -0.4878786E-04
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3912.395
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2514.839
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2997.069
   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]:    1933.820    2997.069  -37127.610
   Vstar,Hstar,Mstar                     [-]:       0.972       5.517      -4.271
   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.517      -4.271       0.972
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.271
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.075
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.346
 > f_VHM_fun                             [-]: -0.1678013E-04
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3328.403
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2549.705
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2139.456
   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]:    1968.686    2139.456  -28201.128
   Vstar,Hstar,Mstar                     [-]:       0.990       3.938      -3.244
   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.938      -3.244       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.144
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.085
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.228
 > f_VHM_fun                             [-]:  0.1429273E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00143
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2956.345
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2560.270
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1478.173
   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]:    1979.251    1478.173  -21281.421
   Vstar,Hstar,Mstar                     [-]:       0.995       2.721      -2.448
   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.721      -2.448       0.995
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.073
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.092
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.168
 > f_VHM_fun                             [-]: -0.2172757E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00217
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2728.390
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2563.848
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     933.164
   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]:    1982.829     933.164  -15566.885
   Vstar,Hstar,Mstar                     [-]:       0.997       1.718      -1.791
   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.718      -1.791       0.997
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      13.371       1.682     -33.965
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      13.371       1.682     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.033
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.098
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.139
 > f_VHM_fun                             [-]: -0.8341883E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00834
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2605.142
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2565.564
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     452.378
   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]:    1984.545     452.378  -10522.489
   Vstar,Hstar,Mstar                     [-]:       0.998       0.833      -1.211
   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.833      -1.211       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.103
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.122
 > f_VHM_fun                             [-]: -0.8635608E-02
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00864
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2566.708
   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]:       7.545
   Vmax                                 [kN]:    1989.287
   Hfac=L   D su_av_L                   [kN]:     543.259
   Mfac=L^2 D su_av_L                  [kNm]:    8692.144
   aMH , bMH                             [-]:      13.371       1.682
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.965
   aVH , bVH                             [-]:       4.656       3.125
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2566.708
   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]:    1985.689       0.000   -5775.094
   Vstar,Hstar,Mstar                     [-]:       0.998       0.000      -0.664
   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.664       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.107
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.108
 > f_VHM_fun                             [-]: -0.4920409E-05
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.63
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.55    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1082.04686
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      907.23997
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.56    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.0    
   qv_tip = Nc * su_av_tip             [kPa]     126.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.200E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     3084.28012
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1989.28683
   Qleaking   = Fo                      [kN]:     1082.04686
 > Vmax_fun                             [kN]:     1989.28683
Vmax_fun : END --- >


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


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



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

    2  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
    2  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5

    3  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
    3  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5

    4  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
    4  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5

    5  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
    5  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5

    6  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
    6  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5

    7  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
    7  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5

    8  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
    8  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5

    9  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
    9  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5

   10  1   4.50  15.00   3.33      1778.0      5332.3      1456.4       548.0         7.2
   10  2   4.50  16.00   3.56      1989.3      6024.5      1635.1       581.0         7.5
 ----------------------------------------------------------------------------------------
  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   5316.7   548.0       5  -0.297E-05   ! ? missing text
    1  2   4.50  16.00  10.50  0.00   6017.7   581.0       5  -0.295E-05   ! ? missing text

    2  1   4.50  15.00  10.00 10.00   5325.9   548.0       5  -0.690E-05   ! ? missing text
    2  2   4.50  16.00  10.50 10.00   6076.6   581.0       5  -0.751E-05   ! ? missing text

    3  1   4.50  15.00  10.00 20.00   5023.7   548.0       6   0.543E-06   ! ? missing text
    3  2   4.50  16.00  10.50 20.00   5723.0   581.0       6   0.682E-06   ! ? missing text

    4  1   4.50  15.00  10.00 30.00   4247.2   548.0       5  -0.105E-03   ! ? missing text
    4  2   4.50  16.00  10.50 30.00   4779.6   581.0       5  -0.488E-04   ! ? missing text

    5  1   4.50  15.00  10.00 40.00   3509.4   548.0       4  -0.175E-03   ! ? missing text
    5  2   4.50  16.00  10.50 40.00   3912.4   581.0       6  -0.168E-04   ! ? missing text

    6  1   4.50  15.00  10.00 50.00   2997.8   548.0       6  -0.198E-02   ! ? missing text
    6  2   4.50  16.00  10.50 50.00   3328.4   581.0       6   0.143E-02   ! ? missing text

    7  1   4.50  15.00  10.00 60.00   2669.0   548.0       9  -0.168E-03   ! ? missing text
    7  2   4.50  16.00  10.50 60.00   2956.3   581.0       8  -0.217E-02   ! ? missing text

    8  1   4.50  15.00  10.00 70.00   2465.7   548.0       8  -0.267E-02   ! ? missing text
    8  2   4.50  16.00  10.50 70.00   2728.4   581.0      10  -0.834E-02   ! ? missing text

    9  1   4.50  15.00  10.00 80.00   2355.5   548.0       9  -0.124E-02   ! ? missing text
    9  2   4.50  16.00  10.50 80.00   2605.1   581.0       8  -0.864E-02   ! ? missing text

   10  1   4.50  15.00  10.00 90.00   2320.5   548.0      11  -0.120E-02   ! ? missing text
   10  2   4.50  16.00  10.50 90.00   2566.7   581.0      24  -0.492E-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.300
  Hmax     [kN] =    0.533E+04
  Vmax     [kN] =    0.178E+04
  theta*        = resistance inclination at lug
  Vlug_theta*, Hlug_theta*
                = corresponding VH FACTORED resistances
  T_theta*      = sqrt( Vlug_theta*^2 + Hlug_theta*^2)
------------------------------------------------------------------------------
 theta*     Vlug_theta*  Hlug_theta*  T_theta*   Vlug_theta*  Hlug_theta*  row
  [deg]        [kN]         [kN]        [kN]      / Vmax [-]   / Hmax [-]  [-]
------------------------------------------------------------------------------
    90.0       0.178E+04   0.109E-12   0.178E+04       0.999    0.000        1
    89.0       0.178E+04    31.0       0.178E+04       0.999    0.006        2
    88.0       0.178E+04    62.0       0.178E+04       0.999    0.012        3
    87.0       0.178E+04    93.1       0.178E+04       0.999    0.017        4
    86.0       0.178E+04    124.       0.178E+04       0.999    0.023        5
    85.0       0.178E+04    155.       0.178E+04       0.999    0.029        6
    84.0       0.178E+04    187.       0.179E+04       0.999    0.035        7
    83.0       0.178E+04    218.       0.179E+04       0.999    0.041        8
    82.0       0.178E+04    250.       0.179E+04       0.999    0.047        9
    81.0       0.178E+04    281.       0.180E+04       0.999    0.053       10
    80.0       0.178E+04    313.       0.180E+04       0.999    0.059       11
    79.0       0.178E+04    345.       0.181E+04       0.999    0.065       12
    78.0       0.178E+04    377.       0.182E+04       0.999    0.071       13
    77.0       0.178E+04    410.       0.182E+04       0.999    0.077       14
    76.0       0.178E+04    443.       0.183E+04       0.999    0.083       15
    75.0       0.178E+04    476.       0.184E+04       0.999    0.089       16
    74.0       0.178E+04    509.       0.185E+04       0.999    0.095       17
    73.0       0.178E+04    543.       0.186E+04       0.999    0.102       18
    72.0       0.178E+04    577.       0.187E+04       0.999    0.108       19
    71.0       0.178E+04    611.       0.188E+04       0.999    0.115       20
    70.0       0.178E+04    646.       0.189E+04       0.999    0.121       21
    69.0       0.178E+04    681.       0.190E+04       0.998    0.128       22
    68.0       0.178E+04    717.       0.191E+04       0.998    0.135       23
    67.0       0.178E+04    753.       0.193E+04       0.998    0.141       24
    66.0       0.177E+04    790.       0.194E+04       0.998    0.148       25
    65.0       0.177E+04    828.       0.196E+04       0.998    0.155       26
    64.0       0.177E+04    865.       0.197E+04       0.998    0.162       27
    63.0       0.177E+04    904.       0.199E+04       0.998    0.170       28
    62.0       0.177E+04    944.       0.201E+04       0.998    0.177       29
    61.0       0.177E+04    983.       0.203E+04       0.997    0.184       30
    60.0       0.177E+04   0.102E+04   0.205E+04       0.997    0.192       31
    59.0       0.177E+04   0.107E+04   0.207E+04       0.997    0.200       32
    58.0       0.177E+04   0.111E+04   0.209E+04       0.997    0.208       33
    57.0       0.177E+04   0.115E+04   0.211E+04       0.997    0.216       34
    56.0       0.177E+04   0.120E+04   0.214E+04       0.997    0.224       35
    55.0       0.177E+04   0.124E+04   0.216E+04       0.997    0.233       36
    54.0       0.177E+04   0.129E+04   0.219E+04       0.996    0.241       37
    53.0       0.177E+04   0.133E+04   0.222E+04       0.996    0.250       38
    52.0       0.177E+04   0.138E+04   0.225E+04       0.996    0.259       39
    51.0       0.177E+04   0.143E+04   0.228E+04       0.995    0.269       40
    50.0       0.177E+04   0.148E+04   0.231E+04       0.995    0.278       41
    49.0       0.177E+04   0.154E+04   0.234E+04       0.995    0.288       42
    48.0       0.177E+04   0.159E+04   0.238E+04       0.994    0.298       43
    47.0       0.177E+04   0.165E+04   0.242E+04       0.994    0.309       44
    46.0       0.177E+04   0.170E+04   0.245E+04       0.993    0.320       45
    45.0       0.176E+04   0.176E+04   0.250E+04       0.992    0.331       46
    44.0       0.176E+04   0.183E+04   0.254E+04       0.991    0.342       47
    43.0       0.176E+04   0.189E+04   0.258E+04       0.991    0.354       48
    42.0       0.176E+04   0.195E+04   0.263E+04       0.990    0.367       49
    41.0       0.176E+04   0.202E+04   0.268E+04       0.989    0.379       50
    40.0       0.176E+04   0.209E+04   0.273E+04       0.988    0.392       51
    39.0       0.175E+04   0.217E+04   0.279E+04       0.986    0.406       52
    38.0       0.175E+04   0.224E+04   0.284E+04       0.985    0.420       53
    37.0       0.175E+04   0.232E+04   0.290E+04       0.983    0.435       54
    36.0       0.174E+04   0.240E+04   0.297E+04       0.981    0.450       55
    35.0       0.174E+04   0.248E+04   0.303E+04       0.979    0.466       56
    34.0       0.174E+04   0.257E+04   0.310E+04       0.976    0.482       57
    33.0       0.173E+04   0.266E+04   0.318E+04       0.973    0.500       58
    32.0       0.172E+04   0.276E+04   0.325E+04       0.969    0.517       59
    31.0       0.172E+04   0.286E+04   0.333E+04       0.965    0.536       60
    30.0       0.171E+04   0.296E+04   0.342E+04       0.961    0.555       61
    29.0       0.170E+04   0.306E+04   0.350E+04       0.955    0.575       62
    28.0       0.169E+04   0.317E+04   0.359E+04       0.949    0.595       63
    27.0       0.167E+04   0.329E+04   0.369E+04       0.942    0.616       64
    26.0       0.166E+04   0.340E+04   0.379E+04       0.933    0.638       65
    25.0       0.164E+04   0.352E+04   0.389E+04       0.924    0.660       66
    24.0       0.162E+04   0.364E+04   0.399E+04       0.912    0.683       67
    23.0       0.160E+04   0.377E+04   0.409E+04       0.899    0.706       68
    22.0       0.157E+04   0.389E+04   0.420E+04       0.884    0.730       69
    21.0       0.154E+04   0.402E+04   0.430E+04       0.867    0.753       70
    20.0       0.151E+04   0.414E+04   0.441E+04       0.848    0.777       71
    19.0       0.147E+04   0.427E+04   0.451E+04       0.826    0.800       72
    18.0       0.142E+04   0.439E+04   0.461E+04       0.801    0.822       73
    17.0       0.138E+04   0.450E+04   0.471E+04       0.774    0.844       74
    16.0       0.132E+04   0.461E+04   0.480E+04       0.744    0.865       75
    15.0       0.126E+04   0.472E+04   0.488E+04       0.711    0.884       76
    14.0       0.120E+04   0.481E+04   0.496E+04       0.675    0.903       77
    13.0       0.113E+04   0.490E+04   0.503E+04       0.636    0.919       78
    12.0       0.106E+04   0.498E+04   0.509E+04       0.595    0.934       79
    11.0        982.       0.505E+04   0.515E+04       0.552    0.947       80
    10.0        901.       0.511E+04   0.519E+04       0.507    0.959       81
    9.00        818.       0.516E+04   0.523E+04       0.460    0.969       82
    8.00        732.       0.521E+04   0.526E+04       0.412    0.977       83
    7.00        644.       0.525E+04   0.528E+04       0.362    0.984       84
    6.00        554.       0.527E+04   0.530E+04       0.312    0.989       85
    5.00        463.       0.530E+04   0.532E+04       0.261    0.993       86
    4.00        372.       0.531E+04   0.533E+04       0.209    0.996       87
    3.00        279.       0.533E+04   0.533E+04       0.157    0.999       88
    2.00        186.       0.533E+04   0.534E+04       0.105    1.000       89
    1.00        93.2       0.534E+04   0.534E+04       0.052    1.002       90
    0.00        0.00       0.535E+04   0.535E+04       0.000    1.002       91
   -1.00       -93.4       0.535E+04   0.535E+04      -0.053    1.003       92
   -2.00       -187.       0.535E+04   0.535E+04      -0.105    1.003       93
   -3.00       -280.       0.535E+04   0.536E+04      -0.158    1.003       94
   -4.00       -374.       0.534E+04   0.536E+04      -0.210    1.002       95
   -5.00       -467.       0.533E+04   0.535E+04      -0.262    1.000       96
   -6.00       -559.       0.532E+04   0.535E+04      -0.314    0.997       97
   -7.00       -650.       0.530E+04   0.534E+04      -0.366    0.993       98
   -8.00       -740.       0.527E+04   0.532E+04      -0.416    0.987       99
   -9.00       -828.       0.523E+04   0.529E+04      -0.466    0.980      100
   -10.0       -913.       0.518E+04   0.526E+04      -0.514    0.971      101
   -11.0       -995.       0.512E+04   0.522E+04      -0.560    0.960      102
   -12.0      -0.107E+04   0.505E+04   0.516E+04      -0.604    0.947      103
   -13.0      -0.115E+04   0.497E+04   0.510E+04      -0.646    0.933      104
   -14.0      -0.122E+04   0.489E+04   0.503E+04      -0.685    0.916      105
   -15.0      -0.128E+04   0.479E+04   0.496E+04      -0.721    0.898      106
   -16.0      -0.134E+04   0.468E+04   0.487E+04      -0.755    0.878      107
   -17.0      -0.140E+04   0.457E+04   0.478E+04      -0.785    0.857      108
   -18.0      -0.145E+04   0.445E+04   0.468E+04      -0.813    0.834      109
   -19.0      -0.149E+04   0.432E+04   0.457E+04      -0.837    0.811      110
   -20.0      -0.153E+04   0.420E+04   0.447E+04      -0.859    0.787      111
   -21.0      -0.156E+04   0.407E+04   0.436E+04      -0.878    0.763      112
   -22.0      -0.159E+04   0.394E+04   0.425E+04      -0.895    0.738      113
   -23.0      -0.162E+04   0.381E+04   0.414E+04      -0.909    0.714      114
   -24.0      -0.164E+04   0.368E+04   0.403E+04      -0.921    0.690      115
   -25.0      -0.166E+04   0.355E+04   0.392E+04      -0.932    0.666      116
   -26.0      -0.167E+04   0.343E+04   0.382E+04      -0.941    0.643      117
   -27.0      -0.169E+04   0.331E+04   0.372E+04      -0.949    0.621      118
   -28.0      -0.170E+04   0.320E+04   0.362E+04      -0.956    0.599      119
   -29.0      -0.171E+04   0.308E+04   0.353E+04      -0.961    0.578      120
   -30.0      -0.172E+04   0.298E+04   0.344E+04      -0.966    0.558      121
   -31.0      -0.173E+04   0.287E+04   0.335E+04      -0.970    0.539      122
   -32.0      -0.173E+04   0.277E+04   0.327E+04      -0.974    0.520      123
   -33.0      -0.174E+04   0.268E+04   0.319E+04      -0.977    0.502      124
   -34.0      -0.174E+04   0.258E+04   0.312E+04      -0.980    0.484      125
   -35.0      -0.175E+04   0.249E+04   0.304E+04      -0.982    0.468      126
   -36.0      -0.175E+04   0.241E+04   0.298E+04      -0.984    0.452      127
   -37.0      -0.175E+04   0.233E+04   0.291E+04      -0.986    0.436      128
   -38.0      -0.176E+04   0.225E+04   0.285E+04      -0.987    0.421      129
   -39.0      -0.176E+04   0.217E+04   0.279E+04      -0.989    0.407      130
   -40.0      -0.176E+04   0.210E+04   0.274E+04      -0.990    0.393      131
   -41.0      -0.176E+04   0.203E+04   0.269E+04      -0.991    0.380      132
   -42.0      -0.176E+04   0.196E+04   0.264E+04      -0.992    0.367      133
   -43.0      -0.176E+04   0.189E+04   0.259E+04      -0.992    0.355      134
   -44.0      -0.177E+04   0.183E+04   0.254E+04      -0.993    0.343      135
   -45.0      -0.177E+04   0.177E+04   0.250E+04      -0.994    0.331      136
   -46.0      -0.177E+04   0.171E+04   0.246E+04      -0.994    0.320      137
   -47.0      -0.177E+04   0.165E+04   0.242E+04      -0.995    0.309      138
   -48.0      -0.177E+04   0.159E+04   0.238E+04      -0.995    0.299      139
   -49.0      -0.177E+04   0.154E+04   0.235E+04      -0.996    0.289      140
   -50.0      -0.177E+04   0.149E+04   0.231E+04      -0.996    0.279      141
   -51.0      -0.177E+04   0.143E+04   0.228E+04      -0.996    0.269      142
   -52.0      -0.177E+04   0.138E+04   0.225E+04      -0.997    0.260      143
   -53.0      -0.177E+04   0.134E+04   0.222E+04      -0.997    0.250      144
   -54.0      -0.177E+04   0.129E+04   0.219E+04      -0.997    0.241      145
   -55.0      -0.177E+04   0.124E+04   0.216E+04      -0.997    0.233      146
   -56.0      -0.177E+04   0.120E+04   0.214E+04      -0.997    0.224      147
   -57.0      -0.177E+04   0.115E+04   0.211E+04      -0.997    0.216      148
   -58.0      -0.177E+04   0.111E+04   0.209E+04      -0.998    0.208      149
   -59.0      -0.177E+04   0.107E+04   0.207E+04      -0.998    0.200      150
   -60.0      -0.177E+04   0.102E+04   0.205E+04      -0.997    0.192      151
   -61.0      -0.177E+04    984.       0.203E+04      -0.998    0.184      152
   -62.0      -0.177E+04    943.       0.201E+04      -0.998    0.177      153
   -63.0      -0.177E+04    904.       0.199E+04      -0.998    0.170      154
   -64.0      -0.178E+04    866.       0.198E+04      -0.998    0.162      155
   -65.0      -0.178E+04    828.       0.196E+04      -0.998    0.155      156
   -66.0      -0.178E+04    790.       0.194E+04      -0.999    0.148      157
   -67.0      -0.178E+04    754.       0.193E+04      -0.999    0.141      158
   -68.0      -0.178E+04    717.       0.192E+04      -0.999    0.135      159
   -69.0      -0.178E+04    682.       0.190E+04      -0.999    0.128      160
   -70.0      -0.178E+04    646.       0.189E+04      -0.999    0.121      161
   -71.0      -0.178E+04    611.       0.188E+04      -0.999    0.115      162
   -72.0      -0.178E+04    577.       0.187E+04      -0.999    0.108      163
   -73.0      -0.178E+04    543.       0.186E+04      -0.999    0.102      164
   -74.0      -0.178E+04    509.       0.185E+04      -0.999    0.095      165
   -75.0      -0.178E+04    476.       0.184E+04      -0.999    0.089      166
   -76.0      -0.178E+04    443.       0.183E+04      -0.999    0.083      167
   -77.0      -0.178E+04    410.       0.182E+04      -0.999    0.077      168
   -78.0      -0.178E+04    378.       0.182E+04      -0.999    0.071      169
   -79.0      -0.178E+04    345.       0.181E+04      -0.999    0.065      170
   -80.0      -0.178E+04    313.       0.180E+04      -0.999    0.059      171
   -81.0      -0.178E+04    281.       0.180E+04      -0.999    0.053      172
   -82.0      -0.178E+04    250.       0.179E+04      -0.999    0.047      173
   -83.0      -0.178E+04    218.       0.179E+04      -0.999    0.041      174
   -84.0      -0.178E+04    187.       0.179E+04      -0.999    0.035      175
   -85.0      -0.178E+04    155.       0.178E+04      -0.999    0.029      176
   -86.0      -0.178E+04    124.       0.178E+04      -0.999    0.023      177
   -87.0      -0.178E+04    93.1       0.178E+04      -0.999    0.017      178
   -88.0      -0.178E+04    62.0       0.178E+04      -0.999    0.012      179
   -89.0      -0.178E+04    31.0       0.178E+04      -0.999    0.006      180
   -90.0      -0.178E+04   0.109E-12   0.178E+04      -0.999    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
