
 p r o g r a m  C A I S S O N _ V H M  w32 version 00.12_b04 
 ------------------------------------  at 22:32:05 on 04-JAN-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: "EX1_MODE10.O01.PDF"
       : program name    : "CAISSON_VHM"

 pidata: rtitle          : "SIMPLIFIED BASE CASE 2 D=6M SU=1.25Z WITH VHM=0, 759, 6071 SHOULD GIVE L=2D = 12M"
       : ptitle          : "VERIFY MODE10"
       : 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        6.00       0.100E-05    50.0    
    40.0    
 m_fac
    1.00    
 gmdata: END --- >


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


 cfdata: < --- BEGIN
 ncgeom
    1
 cgeom
 CIRCULAR
 nfhxp  wt_D_side  wt_D_top
  0   0.100E-01   0.200E-01
 iro_lug  r_lug_user
  1    0.00    
 D  (guess)L
   6.00   9.00
 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.00   0.50    1   0.00
 eta_amh  eta_bmh  eta_phimh
   1.00   1.00   1.00
 ssidat: END --- >


 ldata_support: < --- BEGIN
 nlc   nlfac_lc
  1  1
 lfac_lc(*)
    1.00    
 ldata_support: V H M (l_fac) ***
    0.00        759.       0.607E+04           !
 ldata_support: load_matX; nrows =           1
    0.00        759.       0.607E+04    1.00    

print_VHMloadR_sub:
FACTORED VHM load matrix VHMloadR (# rows/NLC=    1)
 V             = vertical load at seafloor        [kN]
 H             = lateral load at seafloor         [kN]
 M             = moment load at seafloor         [kNm]
 ilc           = load case number                  [-]
 ------------------------------------------!----------------------------------------
    V           H            M          ilc!     V           H          M        ilc
   [kN]        [kN]         [kNm]       [-]!    [kN]        [kN]       [kNm]     [-]
 ------------------------------------------!----------------------------------------
    0.00        759.       0.607E+04     1         0.0       759.0      6071.0     1
 ------------------------------------------!----------------------------------------
    V           H            M          ilc!     V           H          M        ilc
   [kN]        [kN]         [kNm]       [-]!    [kN]        [kN]       [kNm]     [-]
 ------------------------------------------!----------------------------------------
 ldata_support: END --- >


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

 datain_caisson: END --- >


clen_theoretical_sub: < --- BEGIN
 ---------------------------------------------------------------- 
 clen_theoretical_sub: SUMMARY OF RESULTS
 icgeom = caisson geometry #        
 ilc    = VHM load case #           
 loadR  = FACTORED V & H load  [kN]
 D   = caisson outer diameter    [m]
 ?L  = guess caisson embedded length   [m] 
 L_Vmax = L to satisfy Vmax and VloadR [m] 
 L_Hmax = L to satisfy Hmax and HloadR [m] 
 L_H0   = L to satisfy H0   and HloadR [m] 
 ------------------------------------------------------------------
 ilc        D    VloadR       HloadR      ?L     L      L      L   
  icgeom   [m]    [kN]         [kN]       [m]    Vmax   Hmax   H0  
 ------------------------------------------------------------------
    1  1   6.00   0.100E-02    759.       9.00   0.00   4.66   7.39
 ------------------------------------------------------------------
 ilc        D    VloadR       HloadR      ?L     L      L      L   
  icgeom   [m]    [kN]         [kN]       [m]    Vmax   Hmax   H0  
 ------------------------------------------------------------------
clen_theoretical_sub: END --- >



caisson_vhm: nlines_warnings :    0
caisson_vhm: nlines_errors   :    0
solve_MODE10_sub: < --- BEGIN
  VHM load case #=    1; caisson geometry #=    1

  < VHMload          =     0.000   759.000  6071.000
  < D,Lguess         =     6.000     9.000
  > D,Lopt           =     6.000    12.000

f_VHM_LD_fun: < --- BEGIN
 < Vload, Hload, Mload                [kN,m]:       0.000     759.000    6071.000
 < L,D                                   [m]:      12.000       6.000

abVH_VHM6_sub: < --- BEGIN
 < L,D                                   [-]:      12.000       6.000
   L/D                                   [-]:       2.000
   ez/L                                  [-]:       0.667
 ez/L , aVH , bVH                        [-]:       0.250       5.028       2.421
 ez/L , aVH , bVH                        [-]:       0.500       5.500       3.000
 ez/L , aVH , bVH                        [-]:       0.667       2.500       3.833
 ez/L , aVH , bVH                        [-]:       0.750       2.500       3.833

 > aVH , bVH                             [-]:        2.50002        3.83333
abVH_VHM6_sub: END --- >


ab_phiMH_VHM6_sub : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00

Np_FIX_VHM6_fun: < --- BEGIN
 < L/D                                   [-]:       2.000
 < ez/L                                  [-]:       0.667
 ez/L , Np                               [-]:       0.250       8.053
 ez/L , Np                               [-]:       0.500      10.367
 ez/L , Np                               [-]:       0.667       9.900
 ez/L , Np                               [-]:       0.750       9.320
 > Np_FIX                                [-]:        9.90002
Np_FIX_VHM6_fun: END --- >


Np_FREE_VHM6_fun: < --- BEGIN
 < L/D                                   [-]:       2.000
 < ez/L                                  [-]:       0.667
 ez/L , Np                               [-]:       0.250       4.515
 ez/L , Np                               [-]:       0.500       3.653
 ez/L , Np                               [-]:       0.667       2.713
 ez/L , Np                               [-]:       0.750       2.079
 > Np_FREE                               [-]:        2.71300
Np_FREE_VHM6_fun: END --- >


DphiMH_VHM6_fun: < --- BEGIN
 < L/D                                   [-]:       2.000
 < ez/L                                  [-]:       0.667
 ez/L , DphiMH                           [-]:       0.250       3.700
 ez/L , DphiMH                           [-]:       0.500       4.000
 ez/L , DphiMH                           [-]:       0.667       2.000
 ez/L , DphiMH                           [-]:       0.750       2.000
 > DphiMH_VHM6_fun                       [-]:        2.00000
DphiMH_VHM6_fun: END --- >


DabMH_VHM6_sub: < --- BEGIN
 < L/D                                   [-]:       2.000
 < ez/L                                  [-]:       0.667
 ez/L , DaMH , DbMH                      [-]:       0.250       0.000       0.241
 ez/L , DaMH , DbMH                      [-]:       0.500       0.000       0.000
 ez/L , DaMH , DbMH                      [-]:       0.667       0.000       0.000
 ez/L , DaMH , DbMH                      [-]:       0.750       0.000       0.000
 > DaMH , DbMH                           [-]:        0.00000        0.00000
DabMH_VHM6_sub: END --- >

   su_av_L                             [kPa]:  7.500
   L/D , ez/L                            [-]:  2.000  0.667
   Np_FIX, Np_FREE                       [-]:  9.900  2.713
   DphiMH                              [deg]:  2.000
    phiMH = -atan(ez/L)-DphiMH         [deg]:-35.690
   DaMH, DbMH                          [deg]:  0.000  0.000
   aMH =  Np_FIX /cos(phiMH) + DaMH      [-]: 12.189
   bMH = -Np_FREE*sin(phiMH) + DbMH      [-]:  1.583
 > aMH, bMH, phiMH             [-],[-],[deg]:       12.18937        1.58276      -35.69006
ab_phiMH_VHM6_sub : END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: T
 < iplug_vmode                           [-]:  1
 < L,D                                   [m]:  12.00   6.00
   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.06
   wt_x = wt_c                           [m]:   0.06
   Di = D - 2 * wt_c                     [m]:   5.88
   su_av_L   = av 0    to L            [kPa]:    7.50    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1102.72038
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1080.66597
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    2.00    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    8.40    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    16.9    
   qv_tip = Nc * su_av_tip             [kPa]     142.    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:    1.12    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.401E+04
   Qwall =           Awall * qv_tip     [kN]:    159.    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     5110.65270
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2342.10048
   Qleaking   = Fo                      [kN]:     1102.72038
 > Vmax_fun                             [kN]:     5110.65270
Vmax_fun : END --- >

   su_av_L = av 0 to L                 [kPa]:       7.500
   aVH,bVH                               [-]:       2.500       3.833
   aMH,bMH                               [-]:      12.189       1.583
   phiMH,phiMH                   [rad] [deg]:      -0.623     -35.690
   Vload, Hload, Mload                [kN,m]:       0.000     759.000    6071.000
   Vmax                                 [kN]:    5110.653
   Vstar = |Vload/Vmax|                  [-]:       0.000
   Vstar = max(Vstar,1)                  [-]:       0.000
   Hstar = Hload/(  L D su_av_L )        [-]:       1.406
   Mstar = Mload/(L^2 D su_av_L )        [-]:       0.937
   f_VHM_LD_fun = f_VHM_fun(Vstar,Hstar,Mstar,aMH,bMH,phiMH,aVH,bVH)

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00   9.00
   X = Hstar; Y = Mstar; Z = Vstar                  1.406       0.937       0.000
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.189       1.583     -35.690
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.189       1.583     -35.690
   aVH , bVH                             [-]:       2.500       3.833
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.002
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.998
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -1.000
 > f_VHM_fun                             [-]:  0.9886767E-05
f_VHM_fun : END --- >

 > f_VHM                                 [-]:  0.9886767E-05
f_VHM_LD_fun: END --- >


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: T
 < iplug_vmode                           [-]:  1
 < L,D                                   [m]:  12.00   6.00
   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.65
   alpha_Fi                              [-]:   0.65
   eta_Fo, eta_REB                       [-]:  1.000  1.000
   wt_c = D * wt_D_side                  [m]:   0.06
   wt_x = wt_c                           [m]:   0.06
   Di = D - 2 * wt_c                     [m]:   5.88
   su_av_L   = av 0    to L            [kPa]:    7.50    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     1102.72038
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     1080.66597
   Fxp=         2*Di*L*alpha_Fi*su_av_L*nfhxp
                                        [kN]:    0.00    
   L_D = L / D                           [-]:    2.00    
   Nc = min( 6*(1+0.2*L_D) , 9)          [-]:    8.40    
   su_av_tip = av L to L+alpha_d_su*D  [kPa]:    16.9    
   qv_tip = Nc * su_av_tip             [kPa]     142.    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:    1.12    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.401E+04
   Qwall =           Awall * qv_tip     [kN]:    159.    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     5110.65270
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     2342.10048
   Qleaking   = Fo                      [kN]:     1102.72038
 > Vmax_fun                             [kN]:     5110.65270
Vmax_fun : END --- >



 -----------------------------------------------
 solve_MODE10_sub: CAISSON RESISTANCES & WEIGHTS
 MODE_VHM, problem solution mode   = 10
 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* (* NOT USED FOR SUPPORT FOUNDATIONS)
 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   6.00  12.00   2.00      5110.7      5346.1      1465.0      1127.5         7.5
 --------------------------------------------------------------------------------------
  ilc      D      L      L/D       Vmax        Hmax         H0          Bw      su_av_L
   icgeom [m]    [m]     [-]       [kN]        [kN]        [kN]        [kN]       [kPa]
 --------------------------------------------------------------------------------------


 ------------------------------------
 solve_MODE10_sub: SUMMARY OF RESULTS
 ilc    = VHM load case #            
 icgeom = caisson geometry #         
 D    = caisson outer diameter   [m] 
 Lopt = caisson embedded length  [m] 
 L/D  = caisson embedment ratio  [-] 
 iter   = # iterations           [-] 
 f_VHM  = VHM yield function     [-] 
 ------------------------------------------------
 ilc        D      L     L/D   ( iter   f_VHM   )
  icgeom   [m]    [m]    [-]   (  [-]    [-]    )
 ------------------------------------------------
    1  1   6.00  12.00   2.00       6   0.989E-05   !
 ------------------------------------------------
 ilc        D      L     L/D   ( iter   f_VHM   )
  icgeom   [m]    [m]    [-]   (  [-]    [-]    )
 ------------------------------------------------
solve_MODE10_sub: END --- >


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