
 p r o g r a m  C A I S S O N _ V H M  w32 version 00.12_b04 
 ------------------------------------  at 11:11:57 on 09-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: "EX6_MODE22_BC1.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
    1
  zdep        bden       su_top      su_bot   
   [m]     [kN/m^3]      [kPa]        [kPa]   
    0.00        8.00       0.100E-06    62.5    
    50.0    
 m_fac
    1.00    
 gmdata: END --- >


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


 cfdata: < --- BEGIN
 ncgeom
    1
 cgeom
 CIRCULAR
 nfhxp  wt_D_side  wt_D_top
  0   0.655E-02   0.655E-02
 iro_lug  r_lug_user
  0    0.00    
 D  L  zlug
   5.00   7.50   5.23
 cfdata: END --- >


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


 chndat: < --- BEGIN
 iscope_Bw  l_fac_Bw  eta_Vst_wall  delta_Bw_user  Bw_caisson_user
    1    1.00        1.02        99.9       0.100E+04
 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 # 1; BDEN=8.00 kN/m^3 is inconsistent with SU_TOP=0.00 kPa.
WARNING. ssidat: reduction factor on caisson tip (reverse) end-bearing, ETA_REB > 1;
WARNING. ssidat: typically ETA_REB around 0.8 for anchor foundations.
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].

solve_MODE22_sub: < --- BEGIN

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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    1684.451
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233       0.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1684.451
   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]:    -299.895    1684.451   -8814.733
   Vstar,Hstar,Mstar                     [-]:      -0.131       9.583      -6.686
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  9.583      -6.686      -0.131
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.991
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.009
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -1.000
 > f_VHM_fun                             [-]: -0.1434777E-06
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    1710.690
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      10.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:     297.058
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1684.701
   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]:      -2.837    1684.701   -8816.040
   Vstar,Hstar,Mstar                     [-]:      -0.001       9.584      -6.687
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  9.584      -6.687      -0.001
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.991
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.009
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -1.000
 > f_VHM_fun                             [-]: -0.2033217E-06
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    1792.505
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      20.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:     613.073
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1684.403
   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]:     313.178    1684.403   -8814.483
   Vstar,Hstar,Mstar                     [-]:       0.137       9.582      -6.686
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  9.582      -6.686       0.137
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.991
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.009
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -1.000
 > f_VHM_fun                             [-]: -0.5798943E-06
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    1938.158
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      30.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:     969.079
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1678.494
   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]:     669.184    1678.494   -8783.558
   Vstar,Hstar,Mstar                     [-]:       0.293       9.549      -6.662
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  9.549      -6.662       0.293
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.984
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.009
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.993
 > f_VHM_fun                             [-]: -0.3595378E-05
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2144.035
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      40.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1378.159
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1642.426
   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]:    1078.264    1642.426   -8594.815
   Vstar,Hstar,Mstar                     [-]:       0.472       9.344      -6.519
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  9.344      -6.519       0.472
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.942
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.009
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.950
 > f_VHM_fun                             [-]:  0.1404985E-07
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2359.652
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      50.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1807.598
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1516.755
   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]:    1507.703    1516.755   -7937.179
   Vstar,Hstar,Mstar                     [-]:       0.660       8.629      -6.020
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  8.629      -6.020       0.660
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.803
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.007
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.811
 > f_VHM_fun                             [-]:  0.9469072E-07
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2504.565
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      60.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2169.017
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1252.282
   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]:    1869.122    1252.282   -6553.194
   Vstar,Hstar,Mstar                     [-]:       0.818       7.124      -4.971
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  7.124      -4.971       0.818
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.547
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.005
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.553
 > f_VHM_fun                             [-]: -0.1243965E-05
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2566.021
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      70.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2411.271
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     877.631
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    2111.376     877.631   -4592.642
   Vstar,Hstar,Mstar                     [-]:       0.924       4.993      -3.484
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  4.993      -3.484       0.924
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.269
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.002
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.271
 > f_VHM_fun                             [-]: -0.5113365E-04
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2582.843
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      80.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2543.604
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     448.506
   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]:    2243.709     448.506   -2347.032
   Vstar,Hstar,Mstar                     [-]:       0.982       2.552      -1.780
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  2.552      -1.780       0.982
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.070
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.001
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.071
 > f_VHM_fun                             [-]:  0.4157660E-06
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2585.219
   L , D                                 [m]:       7.500       5.000
   zlug , theta_lug_star           [m] [deg]:       5.233      90.000
   zlug_cL                               [m]:       5.233
   Bw , Bw_plug                         [kN]:     299.895       0.000
   su_av_L                             [kPa]:       4.688
   Vmax                                 [kN]:    2285.325
   Hfac=L   D su_av_L                   [kN]:     175.781
   Mfac=L^2 D su_av_L                  [kNm]:    1318.359
   aMH , bMH                             [-]:      11.739       1.943
   phiMH , phiMH                 [rad] [deg]:      -0.625     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2585.219
   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]:    2285.325       0.000      -0.000
   Vstar,Hstar,Mstar                     [-]:       1.000       0.000      -0.000
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   5.00   7.50
   X = Hstar; Y = Mstar; Z = Vstar                  0.000      -0.000       1.000
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      11.739       1.943     -35.815
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      11.739       1.943     -35.815
   aVH , bVH                             [-]:       2.000       4.000
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.000
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.000
   Z:|Z| (=|Vstar|) > 1; -(1-|Z|)                    [-]:      -0.000
 > f_VHM_fun                             [-]:  0.8828986E-32
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.150
   wt_c = D * wt_D_side                  [m]:   0.03
   wt_x = wt_c                           [m]:   0.03
   Di = D - 2 * wt_c                     [m]:   4.93
   su_av_L   = av 0    to L            [kPa]:    4.69    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      358.95151
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      354.24925
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    1.50    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    7.80    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    10.9    
   qv_tip = Nc * su_av_tip             [kPa]     85.3    
   Abase = pi * D * D / 4              [m^2]:    19.6    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.193E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2285.32460
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:      713.20076
   Qleaking   = Fo                      [kN]:      358.95151
 > Vmax_fun                             [kN]:     2285.32460
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < iscope_Bw                             [-]:  1
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.033
   wt_x = wt_c                           [m]:  0.033
   wt_t = D * wt_D_top                   [m]:  0.033
   Di = D - 2 * wt_c                     [m]:  4.934
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.511    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    19.6    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    3.83    
   Vst_t =   Atop          * wt_t      [m^3]:   0.643    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    4.48    
   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]:      299.89478
Bw_caisson_fun : END --- >


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:   7.50   5.00
 < wt_D_side (= WT / D)                  [-]:  0.007
 < iplug_Vmode                           [-]:  1
 > Bw_plug_fun                          [kN]:        0.00000
Bw_plug_fun : END --- >



 -----------------------------------------------
 solve_MODE22_sub: CAISSON RESISTANCES & WEIGHTS
 MODE_VHM, problem solution mode =   22
 IPLUG_VMODE, axial capacity mode  =  1
 ISCOPE_BW, caisson weight BW mode =  1
 ilc    = VHM load case #            
 icgeom = caisson geometry #         
 D    = caisson outer diameter       
 L    = caisson embedded length      
 L/D  = caisson embedment ratio      
 Vmax = caisson axial resistance     
 Hmax = caisson FIXED HEAD lateral resistance 
 H0   = caisson FREE  HEAD lateral resistance 
 Bw   = caisson submerged weight (only used for anchor foundations)
 Bw*  = caisson+soil plug submerged weight (only for IPLUG_VMODE=5)
 su_av_L = soil average su between caisson head (seafloor) and tip (L)
 ----------------------------------------------------------------------------------------
  ilc       D      L      L/D       Vmax        Hmax         H0          Bw       su_av_L
   icgeom  [m]    [m]     [-]       [kN]        [kN]        [kN]        [kN]        [kPa]
 ----------------------------------------------------------------------------------------
    1  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
    2  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
    3  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
    4  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
    5  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
    6  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
    7  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
    8  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
    9  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
   10  1   5.00   7.50   1.50      2285.3      1673.3       583.6       299.9         4.7
 ----------------------------------------------------------------------------------------
  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   5.00   7.50   5.23  0.00   1684.5   299.9       5  -0.143E-06   ! ? missing text
    2  1   5.00   7.50   5.23 10.00   1710.7   299.9       5  -0.203E-06   ! ? missing text
    3  1   5.00   7.50   5.23 20.00   1792.5   299.9       5  -0.580E-06   ! ? missing text
    4  1   5.00   7.50   5.23 30.00   1938.2   299.9       5  -0.360E-05   ! ? missing text
    5  1   5.00   7.50   5.23 40.00   2144.0   299.9       6   0.140E-07   ! ? missing text
    6  1   5.00   7.50   5.23 50.00   2359.7   299.9       6   0.947E-07   ! ? missing text
    7  1   5.00   7.50   5.23 60.00   2504.6   299.9       5  -0.124E-05   ! ? missing text
    8  1   5.00   7.50   5.23 70.00   2566.0   299.9       4  -0.511E-04   ! ? missing text
    9  1   5.00   7.50   5.23 80.00   2582.8   299.9       7   0.416E-06   ! ? missing text
   10  1   5.00   7.50   5.23 90.00   2585.2   299.9       3   0.883E-32   ! ? 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] =    7.50
  D         [m] =    5.00
  zlug      [m] =    5.23
  iro_lug   [-] =   0: no offset (lug on centreline)
  m_fac     [-] =   1.000
  Hmax     [kN] =    0.167E+04
  Vmax     [kN] =    0.229E+04
  theta*        = resistance inclination at lug
  Vlug_theta*, Hlug_theta*
                = corresponding VH UNFACTORED resistances
  T_theta*      = sqrt( Vlug_theta*^2 + Hlug_theta*^2)
------------------------------------------------------------------------------
 theta*     Vlug_theta*  Hlug_theta*  T_theta*   Vlug_theta*  Hlug_theta*  row
  [deg]        [kN]         [kN]        [kN]      / Vmax [-]   / Hmax [-]  [-]
------------------------------------------------------------------------------
    90.0       0.229E+04   0.140E-12   0.229E+04       1.000    0.000        1
    89.0       0.229E+04    39.9       0.229E+04       1.000    0.024        2
    88.0       0.228E+04    79.8       0.229E+04       0.999    0.048        3
    87.0       0.228E+04    120.       0.229E+04       0.999    0.071        4
    86.0       0.228E+04    159.       0.229E+04       0.998    0.095        5
    85.0       0.228E+04    199.       0.229E+04       0.996    0.119        6
    84.0       0.227E+04    239.       0.229E+04       0.995    0.143        7
    83.0       0.227E+04    279.       0.229E+04       0.993    0.167        8
    82.0       0.226E+04    318.       0.229E+04       0.991    0.190        9
    81.0       0.226E+04    358.       0.229E+04       0.989    0.214       10
    80.0       0.225E+04    397.       0.229E+04       0.986    0.237       11
    79.0       0.225E+04    437.       0.229E+04       0.983    0.261       12
    78.0       0.224E+04    476.       0.229E+04       0.979    0.284       13
    77.0       0.223E+04    515.       0.229E+04       0.976    0.308       14
    76.0       0.222E+04    554.       0.229E+04       0.972    0.331       15
    75.0       0.221E+04    592.       0.229E+04       0.968    0.354       16
    74.0       0.220E+04    631.       0.229E+04       0.963    0.377       17
    73.0       0.219E+04    669.       0.229E+04       0.958    0.400       18
    72.0       0.218E+04    707.       0.229E+04       0.953    0.423       19
    71.0       0.216E+04    745.       0.229E+04       0.947    0.445       20
    70.0       0.215E+04    783.       0.229E+04       0.941    0.468       21
    69.0       0.214E+04    820.       0.229E+04       0.935    0.490       22
    68.0       0.212E+04    857.       0.229E+04       0.928    0.512       23
    67.0       0.210E+04    893.       0.229E+04       0.921    0.534       24
    66.0       0.209E+04    929.       0.228E+04       0.913    0.555       25
    65.0       0.207E+04    965.       0.228E+04       0.905    0.577       26
    64.0       0.205E+04   0.100E+04   0.228E+04       0.897    0.598       27
    63.0       0.203E+04   0.103E+04   0.228E+04       0.888    0.618       28
    62.0       0.201E+04   0.107E+04   0.228E+04       0.879    0.639       29
    61.0       0.199E+04   0.110E+04   0.227E+04       0.870    0.658       30
    60.0       0.196E+04   0.113E+04   0.227E+04       0.860    0.678       31
    59.0       0.194E+04   0.117E+04   0.226E+04       0.849    0.697       32
    58.0       0.192E+04   0.120E+04   0.226E+04       0.839    0.716       33
    57.0       0.189E+04   0.123E+04   0.225E+04       0.827    0.734       34
    56.0       0.186E+04   0.126E+04   0.225E+04       0.816    0.752       35
    55.0       0.184E+04   0.129E+04   0.224E+04       0.804    0.769       36
    54.0       0.181E+04   0.131E+04   0.224E+04       0.791    0.785       37
    53.0       0.178E+04   0.134E+04   0.223E+04       0.778    0.801       38
    52.0       0.175E+04   0.137E+04   0.222E+04       0.765    0.816       39
    51.0       0.172E+04   0.139E+04   0.221E+04       0.751    0.831       40
    50.0       0.169E+04   0.141E+04   0.220E+04       0.737    0.845       41
    49.0       0.165E+04   0.144E+04   0.219E+04       0.723    0.858       42
    48.0       0.162E+04   0.146E+04   0.218E+04       0.708    0.871       43
    47.0       0.158E+04   0.148E+04   0.217E+04       0.693    0.883       44
    46.0       0.155E+04   0.150E+04   0.215E+04       0.678    0.894       45
    45.0       0.151E+04   0.151E+04   0.214E+04       0.662    0.905       46
    44.0       0.148E+04   0.153E+04   0.213E+04       0.647    0.915       47
    43.0       0.144E+04   0.155E+04   0.211E+04       0.631    0.924       48
    42.0       0.140E+04   0.156E+04   0.210E+04       0.615    0.932       49
    41.0       0.137E+04   0.157E+04   0.208E+04       0.598    0.940       50
    40.0       0.133E+04   0.159E+04   0.207E+04       0.582    0.947       51
    39.0       0.129E+04   0.160E+04   0.205E+04       0.566    0.954       52
    38.0       0.125E+04   0.161E+04   0.204E+04       0.549    0.960       53
    37.0       0.122E+04   0.162E+04   0.202E+04       0.533    0.965       54
    36.0       0.118E+04   0.162E+04   0.201E+04       0.516    0.970       55
    35.0       0.114E+04   0.163E+04   0.199E+04       0.500    0.975       56
    34.0       0.110E+04   0.164E+04   0.198E+04       0.483    0.979       57
    33.0       0.107E+04   0.164E+04   0.196E+04       0.467    0.983       58
    32.0       0.103E+04   0.165E+04   0.195E+04       0.451    0.986       59
    31.0        994.       0.165E+04   0.193E+04       0.435    0.989       60
    30.0        958.       0.166E+04   0.192E+04       0.419    0.991       61
    29.0        921.       0.166E+04   0.190E+04       0.403    0.993       62
    28.0        886.       0.167E+04   0.189E+04       0.388    0.995       63
    27.0        850.       0.167E+04   0.187E+04       0.372    0.997       64
    26.0        815.       0.167E+04   0.186E+04       0.357    0.999       65
    25.0        780.       0.167E+04   0.185E+04       0.341    1.000       66
    24.0        746.       0.168E+04   0.183E+04       0.326    1.001       67
    23.0        712.       0.168E+04   0.182E+04       0.311    1.002       68
    22.0        678.       0.168E+04   0.181E+04       0.297    1.003       69
    21.0        645.       0.168E+04   0.180E+04       0.282    1.004       70
    20.0        612.       0.168E+04   0.179E+04       0.268    1.004       71
    19.0        579.       0.168E+04   0.178E+04       0.253    1.005       72
    18.0        546.       0.168E+04   0.177E+04       0.239    1.005       73
    17.0        514.       0.168E+04   0.176E+04       0.225    1.006       74
    16.0        483.       0.168E+04   0.175E+04       0.211    1.006       75
    15.0        451.       0.168E+04   0.174E+04       0.197    1.006       76
    14.0        420.       0.168E+04   0.174E+04       0.184    1.006       77
    13.0        389.       0.168E+04   0.173E+04       0.170    1.006       78
    12.0        358.       0.168E+04   0.172E+04       0.157    1.007       79
    11.0        327.       0.168E+04   0.172E+04       0.143    1.007       80
    10.0        297.       0.168E+04   0.171E+04       0.130    1.007       81
    9.00        267.       0.168E+04   0.171E+04       0.117    1.007       82
    8.00        237.       0.168E+04   0.170E+04       0.104    1.007       83
    7.00        207.       0.168E+04   0.170E+04       0.091    1.007       84
    6.00        177.       0.168E+04   0.169E+04       0.077    1.007       85
    5.00        147.       0.168E+04   0.169E+04       0.064    1.007       86
    4.00        118.       0.168E+04   0.169E+04       0.052    1.007       87
    3.00        88.3       0.168E+04   0.169E+04       0.039    1.007       88
    2.00        58.8       0.168E+04   0.169E+04       0.026    1.007       89
    1.00        29.4       0.168E+04   0.168E+04       0.013    1.007       90
    0.00        0.00       0.168E+04   0.168E+04       0.000    1.007       91
   -1.00       -29.4       0.168E+04   0.168E+04      -0.013    1.007       92
   -2.00       -58.8       0.168E+04   0.169E+04      -0.026    1.007       93
   -3.00       -88.3       0.168E+04   0.169E+04      -0.039    1.007       94
   -4.00       -118.       0.168E+04   0.169E+04      -0.052    1.007       95
   -5.00       -147.       0.168E+04   0.169E+04      -0.064    1.007       96
   -6.00       -177.       0.168E+04   0.169E+04      -0.077    1.007       97
   -7.00       -207.       0.168E+04   0.170E+04      -0.091    1.007       98
   -8.00       -237.       0.168E+04   0.170E+04      -0.104    1.007       99
   -9.00       -267.       0.168E+04   0.171E+04      -0.117    1.007      100
   -10.0       -297.       0.168E+04   0.171E+04      -0.130    1.007      101
   -11.0       -327.       0.168E+04   0.172E+04      -0.143    1.007      102
   -12.0       -358.       0.168E+04   0.172E+04      -0.157    1.007      103
   -13.0       -389.       0.168E+04   0.173E+04      -0.170    1.006      104
   -14.0       -420.       0.168E+04   0.174E+04      -0.184    1.006      105
   -15.0       -451.       0.168E+04   0.174E+04      -0.197    1.006      106
   -16.0       -483.       0.168E+04   0.175E+04      -0.211    1.006      107
   -17.0       -514.       0.168E+04   0.176E+04      -0.225    1.006      108
   -18.0       -546.       0.168E+04   0.177E+04      -0.239    1.005      109
   -19.0       -579.       0.168E+04   0.178E+04      -0.253    1.005      110
   -20.0       -612.       0.168E+04   0.179E+04      -0.268    1.004      111
   -21.0       -645.       0.168E+04   0.180E+04      -0.282    1.004      112
   -22.0       -678.       0.168E+04   0.181E+04      -0.297    1.003      113
   -23.0       -712.       0.168E+04   0.182E+04      -0.311    1.002      114
   -24.0       -746.       0.168E+04   0.183E+04      -0.326    1.001      115
   -25.0       -780.       0.167E+04   0.185E+04      -0.341    1.000      116
   -26.0       -815.       0.167E+04   0.186E+04      -0.357    0.999      117
   -27.0       -850.       0.167E+04   0.187E+04      -0.372    0.997      118
   -28.0       -886.       0.167E+04   0.189E+04      -0.388    0.995      119
   -29.0       -921.       0.166E+04   0.190E+04      -0.403    0.993      120
   -30.0       -958.       0.166E+04   0.192E+04      -0.419    0.991      121
   -31.0       -994.       0.165E+04   0.193E+04      -0.435    0.989      122
   -32.0      -0.103E+04   0.165E+04   0.195E+04      -0.451    0.986      123
   -33.0      -0.107E+04   0.164E+04   0.196E+04      -0.467    0.983      124
   -34.0      -0.110E+04   0.164E+04   0.198E+04      -0.483    0.979      125
   -35.0      -0.114E+04   0.163E+04   0.199E+04      -0.500    0.975      126
   -36.0      -0.118E+04   0.162E+04   0.201E+04      -0.516    0.970      127
   -37.0      -0.122E+04   0.162E+04   0.202E+04      -0.533    0.965      128
   -38.0      -0.125E+04   0.161E+04   0.204E+04      -0.549    0.960      129
   -39.0      -0.129E+04   0.160E+04   0.205E+04      -0.566    0.954      130
   -40.0      -0.133E+04   0.159E+04   0.207E+04      -0.582    0.947      131
   -41.0      -0.137E+04   0.157E+04   0.208E+04      -0.598    0.940      132
   -42.0      -0.140E+04   0.156E+04   0.210E+04      -0.615    0.932      133
   -43.0      -0.144E+04   0.155E+04   0.211E+04      -0.631    0.924      134
   -44.0      -0.148E+04   0.153E+04   0.213E+04      -0.647    0.915      135
   -45.0      -0.151E+04   0.151E+04   0.214E+04      -0.662    0.905      136
   -46.0      -0.155E+04   0.150E+04   0.215E+04      -0.678    0.894      137
   -47.0      -0.158E+04   0.148E+04   0.217E+04      -0.693    0.883      138
   -48.0      -0.162E+04   0.146E+04   0.218E+04      -0.708    0.871      139
   -49.0      -0.165E+04   0.144E+04   0.219E+04      -0.723    0.858      140
   -50.0      -0.169E+04   0.141E+04   0.220E+04      -0.737    0.845      141
   -51.0      -0.172E+04   0.139E+04   0.221E+04      -0.751    0.831      142
   -52.0      -0.175E+04   0.137E+04   0.222E+04      -0.765    0.816      143
   -53.0      -0.178E+04   0.134E+04   0.223E+04      -0.778    0.801      144
   -54.0      -0.181E+04   0.131E+04   0.224E+04      -0.791    0.785      145
   -55.0      -0.184E+04   0.129E+04   0.224E+04      -0.804    0.769      146
   -56.0      -0.186E+04   0.126E+04   0.225E+04      -0.816    0.752      147
   -57.0      -0.189E+04   0.123E+04   0.225E+04      -0.827    0.734      148
   -58.0      -0.192E+04   0.120E+04   0.226E+04      -0.839    0.716      149
   -59.0      -0.194E+04   0.117E+04   0.226E+04      -0.849    0.697      150
   -60.0      -0.196E+04   0.113E+04   0.227E+04      -0.860    0.678      151
   -61.0      -0.199E+04   0.110E+04   0.227E+04      -0.870    0.658      152
   -62.0      -0.201E+04   0.107E+04   0.228E+04      -0.879    0.639      153
   -63.0      -0.203E+04   0.103E+04   0.228E+04      -0.888    0.618      154
   -64.0      -0.205E+04   0.100E+04   0.228E+04      -0.897    0.598      155
   -65.0      -0.207E+04    965.       0.228E+04      -0.905    0.577      156
   -66.0      -0.209E+04    929.       0.228E+04      -0.913    0.555      157
   -67.0      -0.210E+04    893.       0.229E+04      -0.921    0.534      158
   -68.0      -0.212E+04    857.       0.229E+04      -0.928    0.512      159
   -69.0      -0.214E+04    820.       0.229E+04      -0.935    0.490      160
   -70.0      -0.215E+04    783.       0.229E+04      -0.941    0.468      161
   -71.0      -0.216E+04    745.       0.229E+04      -0.947    0.445      162
   -72.0      -0.218E+04    707.       0.229E+04      -0.953    0.423      163
   -73.0      -0.219E+04    669.       0.229E+04      -0.958    0.400      164
   -74.0      -0.220E+04    631.       0.229E+04      -0.963    0.377      165
   -75.0      -0.221E+04    592.       0.229E+04      -0.968    0.354      166
   -76.0      -0.222E+04    554.       0.229E+04      -0.972    0.331      167
   -77.0      -0.223E+04    515.       0.229E+04      -0.976    0.308      168
   -78.0      -0.224E+04    476.       0.229E+04      -0.979    0.284      169
   -79.0      -0.225E+04    437.       0.229E+04      -0.983    0.261      170
   -80.0      -0.225E+04    397.       0.229E+04      -0.986    0.237      171
   -81.0      -0.226E+04    358.       0.229E+04      -0.989    0.214      172
   -82.0      -0.226E+04    318.       0.229E+04      -0.991    0.190      173
   -83.0      -0.227E+04    279.       0.229E+04      -0.993    0.167      174
   -84.0      -0.227E+04    239.       0.229E+04      -0.995    0.143      175
   -85.0      -0.228E+04    199.       0.229E+04      -0.996    0.119      176
   -86.0      -0.228E+04    159.       0.229E+04      -0.998    0.095      177
   -87.0      -0.228E+04    120.       0.229E+04      -0.999    0.071      178
   -88.0      -0.228E+04    79.8       0.229E+04      -0.999    0.048      179
   -89.0      -0.229E+04    39.9       0.229E+04      -1.000    0.024      180
   -90.0      -0.229E+04   0.140E-12   0.229E+04      -1.000    0.000      181
------------------------------------------------------------------------------
 theta*     Vlug_theta*  Hlug_theta*  T_theta*   Vlug_theta*  Hlug_theta*  row
  [deg]        [kN]         [kN]        [kN]      / Vmax [-]   / Hmax [-]  [-]
------------------------------------------------------------------------------
VHlug_resist_sub: END --- >

solve_MODE22_sub: END --- >


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