
 p r o g r a m  C A I S S O N _ V H M  w32 version 00.12_b04 
 ------------------------------------  at 15:48:59 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.655E-02   0.655E-02
 iro_lug  r_lug_user
  1    3.00    
 D  L  zlug
   6.00  12.00   9.00
 cfdata: END --- >


 ssidat: < --- BEGIN
 alpha_Fo  alpha_Fi  eta_Fo
   0.67   0.67   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]:    4042.918
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000       0.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    4042.918
   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]:    -666.141    4042.918  -36386.259
   Vstar,Hstar,Mstar                     [-]:      -0.292       9.138      -6.854
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  9.138      -6.854      -0.292
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.883
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.110
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.993
 > f_VHM_fun                             [-]: -0.3032915E-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]:  12.00   6.00

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3752.605
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      10.000
   zlug_cL                               [m]:       9.529
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:     651.633
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3695.595
   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]:     -14.508    3695.595  -35215.251
   Vstar,Hstar,Mstar                     [-]:      -0.006       8.353      -6.633
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  8.353      -6.633      -0.006
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.767
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.233
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -1.000
 > f_VHM_fun                             [-]: -0.2471006E-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]:  12.00   6.00

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3507.910
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      20.000
   zlug_cL                               [m]:      10.092
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1199.776
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3296.358
   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]:     533.635    3296.358  -33266.546
   Vstar,Hstar,Mstar                     [-]:       0.234       7.451      -6.266
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  7.451      -6.266       0.234
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.636
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.361
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.997
 > f_VHM_fun                             [-]: -0.5133269E-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]:  12.00   6.00

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3294.529
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      30.000
   zlug_cL                               [m]:      10.732
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1647.265
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2853.146
   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]:     981.124    2853.146  -30620.110
   Vstar,Hstar,Mstar                     [-]:       0.430       6.449      -5.767
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  6.449      -5.767       0.430
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.498
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.473
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.971
 > f_VHM_fun                             [-]: -0.5457665E-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]:  12.00   6.00

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    3096.189
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      40.000
   zlug_cL                               [m]:      11.517
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    1990.192
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    2371.818
   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]:    1324.051    2371.818  -27316.940
   Vstar,Hstar,Mstar                     [-]:       0.580       5.361      -5.145
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  5.361      -5.145       0.580
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.364
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.552
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.915
 > f_VHM_fun                             [-]: -0.3901130E-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]:  12.00   6.00

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2922.712
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      50.000
   zlug_cL                               [m]:      12.575
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2238.927
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1878.683
   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]:    1572.786    1878.683  -23624.928
   Vstar,Hstar,Mstar                     [-]:       0.689       4.246      -4.450
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  4.246      -4.450       0.689
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.245
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.594
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.839
 > f_VHM_fun                             [-]: -0.1397380E-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]:  12.00   6.00

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2788.265
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      60.000
   zlug_cL                               [m]:      14.196
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2414.708
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1394.133
   l_facT, l_fac_BW                      [-]:       1.000       1.000
   Vload=Vlug*l_facT  - Bw*l_fac_Bw
         Hload=Hlug*l_facT
               Mload=-Hload*zlug_cL
   Vstar=Vload/Vmax
         Hstar=Hload/Hfac
               Mstar=Mload/Mfac
   Vload,Hload,Mload         [kN] [kN] [kNm]:    1748.567    1394.133  -19791.319
   Vstar,Hstar,Mstar                     [-]:       0.766       3.151      -3.728
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  3.151      -3.728       0.766
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.150
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.608
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.757
 > f_VHM_fun                             [-]: -0.2205743E-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]:  12.00   6.00

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2701.409
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      70.000
   zlug_cL                               [m]:      17.242
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2538.494
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     923.936
   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]:    1872.353     923.936  -15930.907
   Vstar,Hstar,Mstar                     [-]:       0.820       2.088      -3.001
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  2.088      -3.001       0.820
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.079
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.601
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.679
 > f_VHM_fun                             [-]: -0.2129697E-04
f_VHM_fun : END --- >

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


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2666.534
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      80.000
   zlug_cL                               [m]:      26.014
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2626.023
   Hlug=Tlug*cos(theta_lug_star)        [kN]:     463.039
   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]:    1959.882     463.039  -12045.417
   Vstar,Hstar,Mstar                     [-]:       0.858       1.047      -2.269
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  1.047      -2.269       0.858
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.031
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.580
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.611
 > f_VHM_fun                             [-]: -0.2837511E-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]:  12.00   6.00

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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


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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:    2687.504
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      90.000
   zlug_cL                               [m]:************
   Bw , Bw_plug                         [kN]:     666.141       0.000
   su_av_L                             [kPa]:       6.145
   Vmax                                 [kN]:    2283.406
   Hfac=L   D su_av_L                   [kN]:     442.429
   Mfac=L^2 D su_av_L                  [kNm]:    5309.146
   aMH , bMH                             [-]:      12.141       1.700
   phiMH , phiMH                 [rad] [deg]:      -0.594     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2687.504
   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]:    2021.363       0.000   -8062.512
   Vstar,Hstar,Mstar                     [-]:       0.885       0.000      -1.519
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.000      -1.519       0.885
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.141       1.700     -34.042
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.141       1.700     -34.042
   aVH , bVH                             [-]:       3.533       3.546
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.005
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.548
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.553
 > f_VHM_fun                             [-]:  0.8205130E-04
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:        0.00008
root_Tlug_22_fun: END --- >


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

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

   wt_D_side                             [-]:   0.01
   alpha_Fo                              [-]:   0.67
   alpha_Fi                              [-]:   0.67
   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]:   5.92
   su_av_L   = av 0    to L            [kPa]:    6.14    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:      931.25384
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:      919.05441
   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]:    11.4    
   qv_tip = Nc * su_av_tip             [kPa]     95.6    
   Abase = pi * D * D / 4              [m^2]:    28.3    
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Qtip  = eta_REB * Abase * qv_tip     [kN]:   0.135E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:     2283.40650
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     1850.30825
   Qleaking   = Fo                      [kN]:      931.25384
 > Vmax_fun                             [kN]:     2283.40650
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  2
   wt_D_side                             [-]:  0.007
   wt_D_top                              [-]:  0.007
   wt_c = D * wt_D_side                  [m]:  0.039
   wt_x = wt_c                           [m]:  0.039
   wt_t = D * wt_D_top                   [m]:  0.039
   Di = D - 2 * wt_c                     [m]:  5.921
   Awall = Abase - pi*Dinner*Dinner/4  [m^2]:   0.736    
   Axp   = nfnxp * wt_x * Dinner       [m^2]:    0.00    
   Atop  = pi * D * D / 4              [m^2]:    28.3    
   Vst_w = ( Awall + Axp ) * L         [m^3]:    8.83    
   Vst_t =   Atop          * wt_t      [m^3]:    1.11    
   ETA_VST_WALL                          [-]:    1.00    
   Vst   = Vst_w*eta_Vst_wall + Vst_t  [m^3]:    9.94    
   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]:      666.14120
Bw_caisson_fun : END --- >


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



 -----------------------------------------------
 solve_MODE22_sub: CAISSON RESISTANCES & WEIGHTS
 MODE_VHM, problem solution mode =   22
 IPLUG_VMODE, axial capacity mode  =  1
 ISCOPE_BW, caisson weight BW mode =  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   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
    2  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
    3  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
    4  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
    5  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
    6  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
    7  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
    8  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
    9  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
   10  1   6.00  12.00   2.00      2283.4      4451.2      1343.5       666.1         6.1
 ----------------------------------------------------------------------------------------
  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   6.00  12.00   9.00  0.00   4042.9   666.1       5  -0.303E-05   ! ? missing text
    2  1   6.00  12.00   9.00 10.00   3752.6   666.1       5  -0.247E-05   ! ? missing text
    3  1   6.00  12.00   9.00 20.00   3507.9   666.1       5  -0.513E-05   ! ? missing text
    4  1   6.00  12.00   9.00 30.00   3294.5   666.1       5  -0.546E-05   ! ? missing text
    5  1   6.00  12.00   9.00 40.00   3096.2   666.1       5  -0.390E-05   ! ? missing text
    6  1   6.00  12.00   9.00 50.00   2922.7   666.1       5  -0.140E-05   ! ? missing text
    7  1   6.00  12.00   9.00 60.00   2788.3   666.1       5  -0.221E-06   ! ? missing text
    8  1   6.00  12.00   9.00 70.00   2701.4   666.1       4  -0.213E-04   ! ? missing text
    9  1   6.00  12.00   9.00 80.00   2666.5   666.1       4  -0.284E-05   ! ? missing text
   10  1   6.00  12.00   9.00 90.00   2687.5   666.1       3   0.821E-04   ! ? 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] =   12.00
  D         [m] =    6.00
  zlug      [m] =    9.00
  iro_lug   [-] =   1: offset + D/2 (lug on outer face)
  m_fac     [-] =   1.300
  Hmax     [kN] =    0.445E+04
  Vmax     [kN] =    0.228E+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.217E+04   0.133E-12   0.217E+04       0.950    0.000        1
    89.0       0.217E+04    37.8       0.217E+04       0.949    0.009        2
    88.0       0.217E+04    75.6       0.217E+04       0.948    0.017        3
    87.0       0.216E+04    113.       0.217E+04       0.947    0.025        4
    86.0       0.216E+04    151.       0.217E+04       0.946    0.034        5
    85.0       0.216E+04    189.       0.217E+04       0.945    0.042        6
    84.0       0.215E+04    226.       0.217E+04       0.943    0.051        7
    83.0       0.215E+04    264.       0.217E+04       0.942    0.059        8
    82.0       0.215E+04    302.       0.217E+04       0.941    0.068        9
    81.0       0.214E+04    340.       0.217E+04       0.939    0.076       10
    80.0       0.214E+04    378.       0.217E+04       0.938    0.085       11
    79.0       0.214E+04    416.       0.218E+04       0.936    0.093       12
    78.0       0.213E+04    454.       0.218E+04       0.935    0.102       13
    77.0       0.213E+04    492.       0.219E+04       0.933    0.110       14
    76.0       0.213E+04    530.       0.219E+04       0.931    0.119       15
    75.0       0.212E+04    569.       0.220E+04       0.930    0.128       16
    74.0       0.212E+04    607.       0.220E+04       0.928    0.136       17
    73.0       0.211E+04    646.       0.221E+04       0.926    0.145       18
    72.0       0.211E+04    685.       0.222E+04       0.924    0.154       19
    71.0       0.210E+04    724.       0.223E+04       0.921    0.163       20
    70.0       0.210E+04    764.       0.223E+04       0.919    0.172       21
    69.0       0.209E+04    804.       0.224E+04       0.917    0.181       22
    68.0       0.209E+04    844.       0.225E+04       0.915    0.190       23
    67.0       0.208E+04    884.       0.226E+04       0.912    0.199       24
    66.0       0.208E+04    925.       0.227E+04       0.909    0.208       25
    65.0       0.207E+04    966.       0.228E+04       0.907    0.217       26
    64.0       0.206E+04   0.101E+04   0.230E+04       0.904    0.226       27
    63.0       0.206E+04   0.105E+04   0.231E+04       0.901    0.236       28
    62.0       0.205E+04   0.109E+04   0.232E+04       0.898    0.245       29
    61.0       0.204E+04   0.113E+04   0.234E+04       0.895    0.254       30
    60.0       0.204E+04   0.118E+04   0.235E+04       0.891    0.264       31
    59.0       0.203E+04   0.122E+04   0.237E+04       0.888    0.274       32
    58.0       0.202E+04   0.126E+04   0.238E+04       0.884    0.283       33
    57.0       0.201E+04   0.131E+04   0.240E+04       0.881    0.293       34
    56.0       0.200E+04   0.135E+04   0.241E+04       0.877    0.303       35
    55.0       0.199E+04   0.140E+04   0.243E+04       0.873    0.313       36
    54.0       0.198E+04   0.144E+04   0.245E+04       0.868    0.324       37
    53.0       0.197E+04   0.149E+04   0.247E+04       0.864    0.334       38
    52.0       0.196E+04   0.153E+04   0.249E+04       0.859    0.344       39
    51.0       0.195E+04   0.158E+04   0.251E+04       0.854    0.355       40
    50.0       0.194E+04   0.163E+04   0.253E+04       0.849    0.365       41
    49.0       0.193E+04   0.167E+04   0.255E+04       0.843    0.376       42
    48.0       0.191E+04   0.172E+04   0.257E+04       0.838    0.387       43
    47.0       0.190E+04   0.177E+04   0.260E+04       0.832    0.398       44
    46.0       0.189E+04   0.182E+04   0.262E+04       0.826    0.409       45
    45.0       0.187E+04   0.187E+04   0.265E+04       0.819    0.420       46
    44.0       0.185E+04   0.192E+04   0.267E+04       0.812    0.432       47
    43.0       0.184E+04   0.197E+04   0.270E+04       0.805    0.443       48
    42.0       0.182E+04   0.202E+04   0.272E+04       0.798    0.455       49
    41.0       0.180E+04   0.207E+04   0.275E+04       0.790    0.466       50
    40.0       0.179E+04   0.213E+04   0.278E+04       0.782    0.478       51
    39.0       0.177E+04   0.218E+04   0.281E+04       0.773    0.490       52
    38.0       0.174E+04   0.223E+04   0.283E+04       0.764    0.502       53
    37.0       0.172E+04   0.229E+04   0.286E+04       0.755    0.514       54
    36.0       0.170E+04   0.234E+04   0.289E+04       0.745    0.526       55
    35.0       0.168E+04   0.240E+04   0.292E+04       0.735    0.538       56
    34.0       0.165E+04   0.245E+04   0.296E+04       0.724    0.550       57
    33.0       0.163E+04   0.250E+04   0.299E+04       0.712    0.563       58
    32.0       0.160E+04   0.256E+04   0.302E+04       0.700    0.575       59
    31.0       0.157E+04   0.261E+04   0.305E+04       0.688    0.587       60
    30.0       0.154E+04   0.267E+04   0.308E+04       0.675    0.600       61
    29.0       0.151E+04   0.273E+04   0.312E+04       0.662    0.612       62
    28.0       0.148E+04   0.278E+04   0.315E+04       0.648    0.625       63
    27.0       0.144E+04   0.284E+04   0.318E+04       0.633    0.637       64
    26.0       0.141E+04   0.289E+04   0.322E+04       0.617    0.649       65
    25.0       0.137E+04   0.295E+04   0.325E+04       0.601    0.662       66
    24.0       0.134E+04   0.300E+04   0.328E+04       0.585    0.674       67
    23.0       0.130E+04   0.305E+04   0.332E+04       0.568    0.686       68
    22.0       0.126E+04   0.311E+04   0.335E+04       0.550    0.698       69
    21.0       0.121E+04   0.316E+04   0.338E+04       0.531    0.710       70
    20.0       0.117E+04   0.321E+04   0.342E+04       0.512    0.722       71
    19.0       0.112E+04   0.326E+04   0.345E+04       0.492    0.733       72
    18.0       0.108E+04   0.331E+04   0.348E+04       0.472    0.744       73
    17.0       0.103E+04   0.336E+04   0.352E+04       0.450    0.756       74
    16.0        978.       0.341E+04   0.355E+04       0.428    0.767       75
    15.0        927.       0.346E+04   0.358E+04       0.406    0.777       76
    14.0        874.       0.351E+04   0.361E+04       0.383    0.788       77
    13.0        820.       0.355E+04   0.365E+04       0.359    0.798       78
    12.0        765.       0.360E+04   0.368E+04       0.335    0.808       79
    11.0        708.       0.364E+04   0.371E+04       0.310    0.818       80
    10.0        649.       0.368E+04   0.374E+04       0.284    0.828       81
    9.00        590.       0.372E+04   0.377E+04       0.258    0.837       82
    8.00        529.       0.377E+04   0.380E+04       0.232    0.846       83
    7.00        467.       0.381E+04   0.383E+04       0.205    0.855       84
    6.00        404.       0.384E+04   0.386E+04       0.177    0.864       85
    5.00        340.       0.388E+04   0.390E+04       0.149    0.872       86
    4.00        274.       0.392E+04   0.393E+04       0.120    0.880       87
    3.00        207.       0.395E+04   0.396E+04       0.091    0.888       88
    2.00        139.       0.399E+04   0.399E+04       0.061    0.896       89
    1.00        70.2       0.402E+04   0.402E+04       0.031    0.904       90
    0.00        0.00       0.406E+04   0.406E+04       0.000    0.912       91
   -1.00       -71.4       0.409E+04   0.409E+04      -0.031    0.919       92
   -2.00       -144.       0.412E+04   0.413E+04      -0.063    0.926       93
   -3.00       -218.       0.415E+04   0.416E+04      -0.095    0.933       94
   -4.00       -293.       0.418E+04   0.419E+04      -0.128    0.940       95
   -5.00       -369.       0.421E+04   0.423E+04      -0.161    0.946       96
   -6.00       -445.       0.424E+04   0.426E+04      -0.195    0.952       97
   -7.00       -523.       0.426E+04   0.429E+04      -0.229    0.958       98
   -8.00       -602.       0.428E+04   0.433E+04      -0.264    0.963       99
   -9.00       -682.       0.430E+04   0.436E+04      -0.298    0.967      100
   -10.0       -761.       0.432E+04   0.438E+04      -0.333    0.970      101
   -11.0       -841.       0.433E+04   0.441E+04      -0.368    0.972      102
   -12.0       -921.       0.433E+04   0.443E+04      -0.403    0.974      103
   -13.0      -0.100E+04   0.433E+04   0.445E+04      -0.438    0.974      104
   -14.0      -0.108E+04   0.433E+04   0.446E+04      -0.473    0.972      105
   -15.0      -0.116E+04   0.432E+04   0.447E+04      -0.507    0.970      106
   -16.0      -0.123E+04   0.430E+04   0.447E+04      -0.540    0.966      107
   -17.0      -0.131E+04   0.427E+04   0.447E+04      -0.572    0.960      108
   -18.0      -0.138E+04   0.424E+04   0.446E+04      -0.603    0.952      109
   -19.0      -0.145E+04   0.420E+04   0.444E+04      -0.633    0.943      110
   -20.0      -0.151E+04   0.415E+04   0.442E+04      -0.662    0.932      111
   -21.0      -0.157E+04   0.410E+04   0.439E+04      -0.689    0.920      112
   -22.0      -0.163E+04   0.403E+04   0.435E+04      -0.714    0.906      113
   -23.0      -0.168E+04   0.397E+04   0.431E+04      -0.738    0.891      114
   -24.0      -0.173E+04   0.390E+04   0.426E+04      -0.760    0.875      115
   -25.0      -0.178E+04   0.382E+04   0.421E+04      -0.780    0.858      116
   -26.0      -0.182E+04   0.374E+04   0.416E+04      -0.798    0.840      117
   -27.0      -0.186E+04   0.365E+04   0.410E+04      -0.815    0.821      118
   -28.0      -0.190E+04   0.357E+04   0.404E+04      -0.831    0.802      119
   -29.0      -0.193E+04   0.348E+04   0.398E+04      -0.845    0.782      120
   -30.0      -0.196E+04   0.339E+04   0.392E+04      -0.858    0.762      121
   -31.0      -0.198E+04   0.330E+04   0.385E+04      -0.869    0.742      122
   -32.0      -0.201E+04   0.321E+04   0.379E+04      -0.879    0.722      123
   -33.0      -0.203E+04   0.312E+04   0.373E+04      -0.889    0.702      124
   -34.0      -0.205E+04   0.304E+04   0.366E+04      -0.897    0.682      125
   -35.0      -0.206E+04   0.295E+04   0.360E+04      -0.904    0.663      126
   -36.0      -0.208E+04   0.286E+04   0.354E+04      -0.911    0.643      127
   -37.0      -0.209E+04   0.278E+04   0.348E+04      -0.917    0.624      128
   -38.0      -0.211E+04   0.270E+04   0.342E+04      -0.922    0.606      129
   -39.0      -0.212E+04   0.261E+04   0.336E+04      -0.927    0.587      130
   -40.0      -0.213E+04   0.253E+04   0.331E+04      -0.931    0.569      131
   -41.0      -0.214E+04   0.246E+04   0.325E+04      -0.935    0.552      132
   -42.0      -0.214E+04   0.238E+04   0.320E+04      -0.939    0.535      133
   -43.0      -0.215E+04   0.231E+04   0.315E+04      -0.942    0.518      134
   -44.0      -0.216E+04   0.223E+04   0.310E+04      -0.944    0.502      135
   -45.0      -0.216E+04   0.216E+04   0.306E+04      -0.947    0.486      136
   -46.0      -0.217E+04   0.209E+04   0.301E+04      -0.949    0.470      137
   -47.0      -0.217E+04   0.203E+04   0.297E+04      -0.951    0.455      138
   -48.0      -0.218E+04   0.196E+04   0.293E+04      -0.953    0.440      139
   -49.0      -0.218E+04   0.189E+04   0.289E+04      -0.954    0.426      140
   -50.0      -0.218E+04   0.183E+04   0.285E+04      -0.956    0.411      141
   -51.0      -0.219E+04   0.177E+04   0.281E+04      -0.957    0.398      142
   -52.0      -0.219E+04   0.171E+04   0.278E+04      -0.958    0.384      143
   -53.0      -0.219E+04   0.165E+04   0.274E+04      -0.959    0.371      144
   -54.0      -0.219E+04   0.159E+04   0.271E+04      -0.960    0.358      145
   -55.0      -0.219E+04   0.154E+04   0.268E+04      -0.961    0.345      146
   -56.0      -0.220E+04   0.148E+04   0.265E+04      -0.961    0.333      147
   -57.0      -0.220E+04   0.143E+04   0.262E+04      -0.962    0.320      148
   -58.0      -0.220E+04   0.137E+04   0.259E+04      -0.962    0.308      149
   -59.0      -0.220E+04   0.132E+04   0.256E+04      -0.963    0.297      150
   -60.0      -0.220E+04   0.127E+04   0.254E+04      -0.963    0.285      151
   -61.0      -0.220E+04   0.122E+04   0.251E+04      -0.963    0.274      152
   -62.0      -0.220E+04   0.117E+04   0.249E+04      -0.963    0.263      153
   -63.0      -0.220E+04   0.112E+04   0.247E+04      -0.964    0.252      154
   -64.0      -0.220E+04   0.107E+04   0.245E+04      -0.964    0.241      155
   -65.0      -0.220E+04   0.103E+04   0.243E+04      -0.964    0.231      156
   -66.0      -0.220E+04    980.       0.241E+04      -0.964    0.220      157
   -67.0      -0.220E+04    934.       0.239E+04      -0.964    0.210      158
   -68.0      -0.220E+04    889.       0.237E+04      -0.963    0.200      159
   -69.0      -0.220E+04    844.       0.236E+04      -0.963    0.190      160
   -70.0      -0.220E+04    800.       0.234E+04      -0.963    0.180      161
   -71.0      -0.220E+04    757.       0.233E+04      -0.963    0.170      162
   -72.0      -0.220E+04    714.       0.231E+04      -0.962    0.160      163
   -73.0      -0.220E+04    672.       0.230E+04      -0.962    0.151      164
   -74.0      -0.220E+04    630.       0.228E+04      -0.962    0.141      165
   -75.0      -0.220E+04    588.       0.227E+04      -0.961    0.132      166
   -76.0      -0.219E+04    547.       0.226E+04      -0.961    0.123      167
   -77.0      -0.219E+04    506.       0.225E+04      -0.960    0.114      168
   -78.0      -0.219E+04    466.       0.224E+04      -0.960    0.105      169
   -79.0      -0.219E+04    426.       0.223E+04      -0.959    0.096      170
   -80.0      -0.219E+04    386.       0.222E+04      -0.959    0.087      171
   -81.0      -0.219E+04    346.       0.221E+04      -0.958    0.078      172
   -82.0      -0.219E+04    307.       0.221E+04      -0.957    0.069      173
   -83.0      -0.218E+04    268.       0.220E+04      -0.957    0.060      174
   -84.0      -0.218E+04    229.       0.219E+04      -0.956    0.052      175
   -85.0      -0.218E+04    191.       0.219E+04      -0.955    0.043      176
   -86.0      -0.218E+04    152.       0.218E+04      -0.954    0.034      177
   -87.0      -0.218E+04    114.       0.218E+04      -0.953    0.026      178
   -88.0      -0.217E+04    75.9       0.218E+04      -0.952    0.017      179
   -89.0      -0.217E+04    37.9       0.217E+04      -0.951    0.009      180
   -90.0      -0.217E+04   0.133E-12   0.217E+04      -0.950    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
