
 p r o g r a m  C A I S S O N _ V H M  w32 version 00.12_b04 
 ------------------------------------  at 16:58:06 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
    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
    1
 cgeom
 CIRCULAR
 nfhxp  wt_D_side  wt_D_top
  0   0.780E-02   0.780E-02
 iro_lug  r_lug_user
  1    2.25    
 D  L  zlug
   4.50  17.00  11.50
 cfdata: END --- >


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


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


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


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

 datain_caisson: END --- >


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



caisson_vhm: nlines_warnings :    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. chndat: user-supplied EXTRA buoyant weight, DELTA_BW_USER <= 0 kN.
WARNING. ldata__anchor: ILC =    1; THETA_LUG_DEG = 0. [horizontal chain inclination at lug level].
WARNING. ldata__anchor: ILC =   10; THETA_LUG_DEG = 90. [vertical chain inclination at lug level].
WARNING. ldata__anchor: ILC =   10; THETA_LUG_DEG = 90 & IRO_LUG > 0. [ZLUG_CL at infinite depth].

solve_MODE22_sub: < --- BEGIN

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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    6735.449
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500       0.000
   zlug_cL                               [m]:      11.500
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6735.449
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    -598.182    6735.449  -77457.662
   Vstar,Hstar,Mstar                     [-]:      -0.261      11.148      -7.541
   root_Tlug_22_fun=f_VHM_fun               

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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    6704.340
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      10.000
   zlug_cL                               [m]:      11.897
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1164.196
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6602.486
   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]:     566.014    6602.486  -78548.036
   Vstar,Hstar,Mstar                     [-]:       0.247      10.928      -7.647
   root_Tlug_22_fun=f_VHM_fun               

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

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    6292.905
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      20.000
   zlug_cL                               [m]:      12.319
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2152.300
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    5913.397
   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]:    1554.118    5913.397  -72846.736
   Vstar,Hstar,Mstar                     [-]:       0.679       9.787      -7.092
   root_Tlug_22_fun=f_VHM_fun               

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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    5316.321
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      30.000
   zlug_cL                               [m]:      12.799
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2658.161
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4604.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]:    2059.979    4604.069  -58927.659
   Vstar,Hstar,Mstar                     [-]:       0.901       7.620      -5.737
   root_Tlug_22_fun=f_VHM_fun               

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

 > root_Tlug_22_fun                      [-]:       -0.00007
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    4381.277
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      40.000
   zlug_cL                               [m]:      13.388
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2816.231
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3356.253
   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]:    2218.049    3356.253  -44933.431
   Vstar,Hstar,Mstar                     [-]:       0.970       5.555      -4.375
   root_Tlug_22_fun=f_VHM_fun               

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

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    3734.341
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      50.000
   zlug_cL                               [m]:      14.181
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2860.671
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2400.388
   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]:    2262.489    2400.388  -34040.974
   Vstar,Hstar,Mstar                     [-]:       0.989       3.973      -3.314
   root_Tlug_22_fun=f_VHM_fun               

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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    3319.782
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      60.000
   zlug_cL                               [m]:      15.397
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2875.015
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1659.891
   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]:    2276.833    1659.891  -25557.530
   Vstar,Hstar,Mstar                     [-]:       0.995       2.747      -2.488
   root_Tlug_22_fun=f_VHM_fun               

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

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    3064.121
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      70.000
   zlug_cL                               [m]:      17.682
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2879.332
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1047.991
   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]:    2281.150    1047.991  -18530.393
   Vstar,Hstar,Mstar                     [-]:       0.997       1.735      -1.804
   root_Tlug_22_fun=f_VHM_fun               

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

 > root_Tlug_22_fun                      [-]:       -0.00813
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    2926.855
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      80.000
   zlug_cL                               [m]:      24.260
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2882.390
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     508.243
   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]:    2284.208     508.243  -12330.173
   Vstar,Hstar,Mstar                     [-]:       0.999       0.841      -1.200
   root_Tlug_22_fun=f_VHM_fun               

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

 > root_Tlug_22_fun                      [-]:        0.00161
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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]:    2883.004
   L , D                                 [m]:      17.000       4.500
   zlug , theta_lug_star           [m] [deg]:      11.500      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     598.182       0.000
   su_av_L                             [kPa]:       7.898
   Vmax                                 [kN]:    2287.525
   Hfac=L   D su_av_L                   [kN]:     604.202
   Mfac=L^2 D su_av_L                  [kNm]:   10271.428
   aMH , bMH                             [-]:      13.506       1.688
   phiMH , phiMH                 [rad] [deg]:      -0.593     -33.976
   aVH , bVH                             [-]:       4.818       3.065
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2883.004
   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]:    2284.822       0.000   -6486.759
   Vstar,Hstar,Mstar                     [-]:       0.999       0.000      -0.632
   root_Tlug_22_fun=f_VHM_fun               

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

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.54
   eta_Fo, eta_REB                       [-]:  1.000  0.500
   wt_c = D * wt_D_side                  [m]:   0.04
   wt_x = wt_c                           [m]:   0.04
   Di = D - 2 * wt_c                     [m]:   4.43
   su_av_L   = av 0    to L            [kPa]:    7.90    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1233.80100
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1009.01385
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    3.78    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    9.00    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    14.7    
   qv_tip = Nc * su_av_tip             [kPa]     133.    
   Abase = pi * D * D / 4              [m^2]:    15.9    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.492    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.105E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2287.52520
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2242.81484
   Qleaking   = Fo                      [kN]:     1233.80100
 > Vmax_fun                             [kN]:     2287.52520
Vmax_fun : END --- >


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


Bw_plug_fun : < --- BEGIN
 < L,D                                   [m]:  17.00   4.50
 < wt_D_side (= WT / D)                  [-]:  0.008
 < iplug_Vmode                           [-]:  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 =  2
 ilc    = VHM load case #            
 icgeom = caisson geometry #         
 D    = caisson outer diameter       
 L    = caisson embedded length      
 L/D  = caisson embedment ratio      
 Vmax = caisson axial resistance     
 Hmax = caisson FIXED HEAD lateral resistance 
 H0   = caisson FREE  HEAD lateral resistance 
 Bw   = caisson submerged weight (only used for anchor foundations)
 Bw*  = caisson+soil plug submerged weight (only for IPLUG_VMODE=5)
 su_av_L = soil average su between caisson head (seafloor) and tip (L)
 ----------------------------------------------------------------------------------------
  ilc       D      L      L/D       Vmax        Hmax         H0          Bw       su_av_L
   icgeom  [m]    [m]     [-]       [kN]        [kN]        [kN]        [kN]        [kPa]
 ----------------------------------------------------------------------------------------
    1  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
    2  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
    3  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
    4  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
    5  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
    6  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
    7  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
    8  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
    9  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
   10  1   4.50  17.00   3.78      2287.5      6767.3      1825.2       598.2         7.9
 ----------------------------------------------------------------------------------------
  ilc       D      L      L/D       Vmax        Hmax         H0          Bw       su_av_L
   icgeom  [m]    [m]     [-]       [kN]        [kN]        [kN]        [kN]        [kPa]
 ----------------------------------------------------------------------------------------


 ------------------------------------
 solve_MODE22_sub: SUMMARY OF RESULTS
 ilc    = VHM load case #            
 icgeom = caisson geometry #         
 D      = caisson outer diameter           [m] 
 L      = caisson embedded length          [m] 
 zlug   = caisson lug level                [m] 
 theta_lug = load inclination            [deg] 
 Tlug      = UNFACTORED anchor chain load [kN] 
 Bw        = caisson buoyant weight       [kN] 
 iters     = # iterations                  [-] 
 f_VHM     = VHM yield function            [-] 
 -----------------------------------------------------------------------
 ilc        D      L    zlug  theta    Tlug     Bw    ( iters    f_VHM )
  icgeom   [m]    [m]   [m]   [deg]    [kN]    [kN]   ( [-]      [-]   )
 -----------------------------------------------------------------------
    1  1   4.50  17.00  11.50  0.00   6735.4   598.2       5  -0.300E-05   ! ? missing text
    2  1   4.50  17.00  11.50 10.00   6704.3   598.2       5  -0.742E-05   ! ? missing text
    3  1   4.50  17.00  11.50 20.00   6292.9   598.2       6   0.357E-06   ! ? missing text
    4  1   4.50  17.00  11.50 30.00   5316.3   598.2       5  -0.749E-04   ! ? missing text
    5  1   4.50  17.00  11.50 40.00   4381.3   598.2       4   0.312E-03   ! ? missing text
    6  1   4.50  17.00  11.50 50.00   3734.3   598.2       6  -0.103E-03   ! ? missing text
    7  1   4.50  17.00  11.50 60.00   3319.8   598.2      10   0.376E-03   ! ? missing text
    8  1   4.50  17.00  11.50 70.00   3064.1   598.2      10  -0.813E-02   ! ? missing text
    9  1   4.50  17.00  11.50 80.00   2926.9   598.2       8   0.161E-02   ! ? missing text
   10  1   4.50  17.00  11.50 90.00   2883.0   598.2      24  -0.905E-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] =   17.00
  D         [m] =    4.50
  zlug      [m] =   11.50
  iro_lug   [-] =   1: offset + D/2 (lug on outer face)
  m_fac     [-] =   1.300
  Hmax     [kN] =    0.677E+04
  Vmax     [kN] =    0.229E+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.229E+04   0.140E-12   0.229E+04       0.999    0.000        1
    89.0       0.229E+04    39.9       0.229E+04       0.999    0.006        2
    88.0       0.229E+04    79.8       0.229E+04       1.000    0.012        3
    87.0       0.229E+04    120.       0.229E+04       1.000    0.018        4
    86.0       0.229E+04    160.       0.229E+04       1.000    0.024        5
    85.0       0.229E+04    200.       0.230E+04       1.000    0.030        6
    84.0       0.229E+04    240.       0.230E+04       0.999    0.036        7
    83.0       0.229E+04    281.       0.230E+04       0.999    0.041        8
    82.0       0.229E+04    321.       0.231E+04       1.000    0.047        9
    81.0       0.229E+04    362.       0.231E+04       0.999    0.054       10
    80.0       0.229E+04    403.       0.232E+04       0.999    0.060       11
    79.0       0.229E+04    444.       0.233E+04       1.000    0.066       12
    78.0       0.229E+04    486.       0.234E+04       1.000    0.072       13
    77.0       0.229E+04    528.       0.235E+04       0.999    0.078       14
    76.0       0.229E+04    570.       0.236E+04       0.999    0.084       15
    75.0       0.229E+04    612.       0.237E+04       0.999    0.090       16
    74.0       0.229E+04    655.       0.238E+04       0.999    0.097       17
    73.0       0.229E+04    699.       0.239E+04       0.999    0.103       18
    72.0       0.229E+04    743.       0.240E+04       0.999    0.110       19
    71.0       0.229E+04    787.       0.242E+04       0.999    0.116       20
    70.0       0.229E+04    832.       0.243E+04       0.999    0.123       21
    69.0       0.229E+04    877.       0.245E+04       0.999    0.130       22
    68.0       0.229E+04    923.       0.246E+04       0.999    0.136       23
    67.0       0.229E+04    970.       0.248E+04       0.999    0.143       24
    66.0       0.229E+04   0.102E+04   0.250E+04       0.999    0.150       25
    65.0       0.229E+04   0.107E+04   0.252E+04       0.999    0.157       26
    64.0       0.228E+04   0.111E+04   0.254E+04       0.999    0.165       27
    63.0       0.228E+04   0.116E+04   0.256E+04       0.999    0.172       28
    62.0       0.229E+04   0.122E+04   0.259E+04       0.999    0.180       29
    61.0       0.228E+04   0.127E+04   0.261E+04       0.998    0.187       30
    60.0       0.228E+04   0.132E+04   0.264E+04       0.998    0.195       31
    59.0       0.228E+04   0.137E+04   0.267E+04       0.999    0.203       32
    58.0       0.228E+04   0.143E+04   0.269E+04       0.998    0.211       33
    57.0       0.228E+04   0.148E+04   0.272E+04       0.998    0.219       34
    56.0       0.228E+04   0.154E+04   0.275E+04       0.998    0.227       35
    55.0       0.228E+04   0.160E+04   0.279E+04       0.998    0.236       36
    54.0       0.228E+04   0.166E+04   0.282E+04       0.997    0.245       37
    53.0       0.228E+04   0.172E+04   0.286E+04       0.997    0.254       38
    52.0       0.228E+04   0.178E+04   0.289E+04       0.997    0.263       39
    51.0       0.228E+04   0.185E+04   0.293E+04       0.997    0.273       40
    50.0       0.228E+04   0.191E+04   0.298E+04       0.996    0.283       41
    49.0       0.228E+04   0.198E+04   0.302E+04       0.996    0.293       42
    48.0       0.228E+04   0.205E+04   0.306E+04       0.996    0.303       43
    47.0       0.228E+04   0.212E+04   0.311E+04       0.995    0.314       44
    46.0       0.228E+04   0.220E+04   0.316E+04       0.995    0.325       45
    45.0       0.227E+04   0.227E+04   0.322E+04       0.994    0.336       46
    44.0       0.227E+04   0.235E+04   0.327E+04       0.993    0.348       47
    43.0       0.227E+04   0.243E+04   0.333E+04       0.993    0.360       48
    42.0       0.227E+04   0.252E+04   0.339E+04       0.992    0.372       49
    41.0       0.227E+04   0.261E+04   0.345E+04       0.991    0.385       50
    40.0       0.226E+04   0.270E+04   0.352E+04       0.990    0.399       51
    39.0       0.226E+04   0.279E+04   0.359E+04       0.988    0.413       52
    38.0       0.226E+04   0.289E+04   0.367E+04       0.987    0.427       53
    37.0       0.225E+04   0.299E+04   0.374E+04       0.985    0.442       54
    36.0       0.225E+04   0.310E+04   0.383E+04       0.983    0.457       55
    35.0       0.224E+04   0.320E+04   0.391E+04       0.981    0.473       56
    34.0       0.224E+04   0.332E+04   0.400E+04       0.978    0.490       57
    33.0       0.223E+04   0.343E+04   0.410E+04       0.975    0.507       58
    32.0       0.222E+04   0.356E+04   0.419E+04       0.971    0.526       59
    31.0       0.221E+04   0.368E+04   0.430E+04       0.967    0.544       60
    30.0       0.220E+04   0.381E+04   0.440E+04       0.962    0.563       61
    29.0       0.219E+04   0.395E+04   0.451E+04       0.957    0.583       62
    28.0       0.217E+04   0.409E+04   0.463E+04       0.950    0.604       63
    27.0       0.216E+04   0.423E+04   0.475E+04       0.942    0.625       64
    26.0       0.214E+04   0.438E+04   0.487E+04       0.933    0.647       65
    25.0       0.211E+04   0.453E+04   0.500E+04       0.923    0.669       66
    24.0       0.208E+04   0.468E+04   0.512E+04       0.911    0.692       67
    23.0       0.205E+04   0.484E+04   0.525E+04       0.897    0.714       68
    22.0       0.202E+04   0.499E+04   0.538E+04       0.881    0.737       69
    21.0       0.197E+04   0.514E+04   0.551E+04       0.863    0.760       70
    20.0       0.193E+04   0.530E+04   0.564E+04       0.843    0.783       71
    19.0       0.187E+04   0.545E+04   0.576E+04       0.820    0.805       72
    18.0       0.182E+04   0.559E+04   0.588E+04       0.794    0.826       73
    17.0       0.175E+04   0.573E+04   0.599E+04       0.766    0.846       74
    16.0       0.168E+04   0.586E+04   0.610E+04       0.735    0.866       75
    15.0       0.160E+04   0.598E+04   0.619E+04       0.701    0.884       76
    14.0       0.152E+04   0.610E+04   0.628E+04       0.664    0.901       77
    13.0       0.143E+04   0.620E+04   0.636E+04       0.626    0.916       78
    12.0       0.134E+04   0.629E+04   0.644E+04       0.585    0.930       79
    11.0       0.124E+04   0.638E+04   0.650E+04       0.542    0.943       80
    10.0       0.114E+04   0.645E+04   0.655E+04       0.497    0.953       81
    9.00       0.103E+04   0.652E+04   0.660E+04       0.451    0.963       82
    8.00        923.       0.657E+04   0.663E+04       0.404    0.971       83
    7.00        812.       0.661E+04   0.666E+04       0.355    0.977       84
    6.00        699.       0.665E+04   0.669E+04       0.306    0.983       85
    5.00        585.       0.668E+04   0.671E+04       0.256    0.987       86
    4.00        469.       0.671E+04   0.672E+04       0.205    0.991       87
    3.00        352.       0.672E+04   0.673E+04       0.154    0.994       88
    2.00        235.       0.674E+04   0.674E+04       0.103    0.996       89
    1.00        118.       0.675E+04   0.675E+04       0.052    0.997       90
    0.00        0.00       0.676E+04   0.676E+04       0.000    0.999       91
   -1.00       -118.       0.677E+04   0.677E+04      -0.052    1.000       92
   -2.00       -236.       0.677E+04   0.678E+04      -0.103    1.001       93
   -3.00       -355.       0.677E+04   0.678E+04      -0.155    1.001       94
   -4.00       -474.       0.677E+04   0.679E+04      -0.207    1.001       95
   -5.00       -592.       0.676E+04   0.679E+04      -0.259    1.000       96
   -6.00       -709.       0.675E+04   0.679E+04      -0.310    0.997       97
   -7.00       -826.       0.673E+04   0.678E+04      -0.361    0.994       98
   -8.00       -941.       0.670E+04   0.676E+04      -0.411    0.989       99
   -9.00      -0.105E+04   0.665E+04   0.674E+04      -0.461    0.983      100
   -10.0      -0.116E+04   0.660E+04   0.670E+04      -0.509    0.975      101
   -11.0      -0.127E+04   0.654E+04   0.666E+04      -0.555    0.966      102
   -12.0      -0.137E+04   0.646E+04   0.660E+04      -0.600    0.954      103
   -13.0      -0.147E+04   0.637E+04   0.654E+04      -0.643    0.941      104
   -14.0      -0.156E+04   0.626E+04   0.646E+04      -0.683    0.926      105
   -15.0      -0.165E+04   0.615E+04   0.637E+04      -0.720    0.909      106
   -16.0      -0.173E+04   0.602E+04   0.627E+04      -0.755    0.890      107
   -17.0      -0.180E+04   0.589E+04   0.615E+04      -0.787    0.870      108
   -18.0      -0.186E+04   0.574E+04   0.603E+04      -0.815    0.848      109
   -19.0      -0.192E+04   0.559E+04   0.591E+04      -0.841    0.825      110
   -20.0      -0.197E+04   0.543E+04   0.577E+04      -0.863    0.802      111
   -21.0      -0.202E+04   0.526E+04   0.564E+04      -0.883    0.778      112
   -22.0      -0.206E+04   0.510E+04   0.550E+04      -0.900    0.753      113
   -23.0      -0.209E+04   0.493E+04   0.536E+04      -0.915    0.729      114
   -24.0      -0.212E+04   0.477E+04   0.522E+04      -0.928    0.704      115
   -25.0      -0.215E+04   0.460E+04   0.508E+04      -0.938    0.680      116
   -26.0      -0.217E+04   0.444E+04   0.495E+04      -0.948    0.657      117
   -27.0      -0.219E+04   0.429E+04   0.481E+04      -0.955    0.634      118
   -28.0      -0.220E+04   0.414E+04   0.469E+04      -0.962    0.612      119
   -29.0      -0.221E+04   0.399E+04   0.457E+04      -0.967    0.590      120
   -30.0      -0.222E+04   0.385E+04   0.445E+04      -0.972    0.569      121
   -31.0      -0.223E+04   0.372E+04   0.434E+04      -0.976    0.549      122
   -32.0      -0.224E+04   0.359E+04   0.423E+04      -0.979    0.530      123
   -33.0      -0.225E+04   0.346E+04   0.413E+04      -0.982    0.511      124
   -34.0      -0.225E+04   0.334E+04   0.403E+04      -0.985    0.493      125
   -35.0      -0.226E+04   0.322E+04   0.393E+04      -0.987    0.476      126
   -36.0      -0.226E+04   0.311E+04   0.385E+04      -0.988    0.460      127
   -37.0      -0.226E+04   0.300E+04   0.376E+04      -0.990    0.444      128
   -38.0      -0.227E+04   0.290E+04   0.368E+04      -0.991    0.429      129
   -39.0      -0.227E+04   0.280E+04   0.361E+04      -0.992    0.414      130
   -40.0      -0.227E+04   0.271E+04   0.353E+04      -0.993    0.400      131
   -41.0      -0.227E+04   0.262E+04   0.347E+04      -0.994    0.387      132
   -42.0      -0.228E+04   0.253E+04   0.340E+04      -0.995    0.373      133
   -43.0      -0.228E+04   0.244E+04   0.334E+04      -0.995    0.361      134
   -44.0      -0.228E+04   0.236E+04   0.328E+04      -0.996    0.349      135
   -45.0      -0.228E+04   0.228E+04   0.322E+04      -0.996    0.337      136
   -46.0      -0.228E+04   0.220E+04   0.317E+04      -0.996    0.325      137
   -47.0      -0.228E+04   0.213E+04   0.312E+04      -0.997    0.314      138
   -48.0      -0.228E+04   0.205E+04   0.307E+04      -0.997    0.304      139
   -49.0      -0.228E+04   0.198E+04   0.302E+04      -0.997    0.293      140
   -50.0      -0.228E+04   0.192E+04   0.298E+04      -0.998    0.283      141
   -51.0      -0.228E+04   0.185E+04   0.294E+04      -0.998    0.273      142
   -52.0      -0.228E+04   0.178E+04   0.290E+04      -0.998    0.264      143
   -53.0      -0.228E+04   0.172E+04   0.286E+04      -0.998    0.254      144
   -54.0      -0.228E+04   0.166E+04   0.282E+04      -0.998    0.245      145
   -55.0      -0.228E+04   0.160E+04   0.279E+04      -0.999    0.236      146
   -56.0      -0.228E+04   0.154E+04   0.276E+04      -0.999    0.228      147
   -57.0      -0.228E+04   0.148E+04   0.272E+04      -0.999    0.219      148
   -58.0      -0.229E+04   0.143E+04   0.270E+04      -0.999    0.211      149
   -59.0      -0.228E+04   0.137E+04   0.266E+04      -0.998    0.203      150
   -60.0      -0.229E+04   0.132E+04   0.264E+04      -0.999    0.195      151
   -61.0      -0.228E+04   0.127E+04   0.261E+04      -0.999    0.187      152
   -62.0      -0.229E+04   0.122E+04   0.259E+04      -0.999    0.180      153
   -63.0      -0.229E+04   0.116E+04   0.257E+04      -0.999    0.172      154
   -64.0      -0.229E+04   0.111E+04   0.254E+04      -0.999    0.165      155
   -65.0      -0.229E+04   0.107E+04   0.252E+04      -0.999    0.158      156
   -66.0      -0.229E+04   0.102E+04   0.250E+04      -0.999    0.150      157
   -67.0      -0.229E+04    970.       0.248E+04      -0.999    0.143      158
   -68.0      -0.229E+04    924.       0.247E+04      -0.999    0.136      159
   -69.0      -0.229E+04    878.       0.245E+04      -0.999    0.130      160
   -70.0      -0.229E+04    832.       0.243E+04      -0.999    0.123      161
   -71.0      -0.229E+04    787.       0.242E+04      -0.999    0.116      162
   -72.0      -0.229E+04    743.       0.240E+04      -0.999    0.110      163
   -73.0      -0.229E+04    699.       0.239E+04      -0.999    0.103      164
   -74.0      -0.229E+04    656.       0.238E+04      -1.000    0.097      165
   -75.0      -0.229E+04    613.       0.237E+04      -0.999    0.091      166
   -76.0      -0.229E+04    570.       0.236E+04      -0.999    0.084      167
   -77.0      -0.229E+04    528.       0.235E+04      -0.999    0.078      168
   -78.0      -0.229E+04    486.       0.234E+04      -1.000    0.072      169
   -79.0      -0.229E+04    444.       0.233E+04      -0.999    0.066      170
   -80.0      -0.229E+04    404.       0.232E+04      -1.000    0.060      171
   -81.0      -0.229E+04    362.       0.232E+04      -1.000    0.054      172
   -82.0      -0.229E+04    321.       0.231E+04      -1.000    0.047      173
   -83.0      -0.229E+04    281.       0.231E+04      -1.000    0.042      174
   -84.0      -0.229E+04    240.       0.230E+04      -1.000    0.036      175
   -85.0      -0.229E+04    200.       0.229E+04      -0.999    0.030      176
   -86.0      -0.229E+04    160.       0.229E+04      -0.999    0.024      177
   -87.0      -0.229E+04    120.       0.229E+04      -0.999    0.018      178
   -88.0      -0.229E+04    79.8       0.229E+04      -0.999    0.012      179
   -89.0      -0.229E+04    39.9       0.229E+04      -0.999    0.006      180
   -90.0      -0.229E+04   0.140E-12   0.229E+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
