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


 datain_caisson: < --- BEGIN

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

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

 pidata: END --- >


 gmdata: < --- BEGIN
 numsl
    1
  zdep        bden       su_top      su_bot   
   [m]     [kN/m^3]      [kPa]        [kPa]   
    0.00        10.0        2.00        125.    
    50.0    
 m_fac
    1.00    
 gmdata: END --- >


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


 cfdata: < --- BEGIN
 ncgeom
    1
 cgeom
 CIRCULAR
 nfhxp  wt_D_side  wt_D_top
  0   0.655E-02   0.655E-02
 iro_lug  r_lug_user
  0    0.00    
 D  L  zlug
   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
   1.00   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
    0    1.00        1.00        0.00        0.00    
 dtheta_lug_deg  psi_twist_deg
    0.00        0.00    
 chndat: END --- >


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


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

 datain_caisson: END --- >


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



caisson_vhm: nlines_warnings :    3
caisson_vhm: nlines_errors   :    0

WARNING: gmdata. soil layer # 1; BDEN=10.00 kN/m^3 is inconsistent with SU_TOP=2.00 kPa.
WARNING. ldata__anchor: ILC =    1; THETA_LUG_DEG = 0. [horizontal chain inclination at lug level].
WARNING. ldata__anchor: ILC =   10; THETA_LUG_DEG = 90. [vertical chain inclination at lug level].

solve_MODE22_sub: < --- BEGIN

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

root_Tlug_22_fun: < --- BEGIN
 < Tlug                                 [kN]:   11465.383
   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]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:       0.000
   Hlug=Tlug*cos(theta_lug_star)        [kN]:   11465.383
   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]:       0.000   11465.383 -103188.443
   Vstar,Hstar,Mstar                     [-]:       0.000       9.501      -7.126
   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.501      -7.126       0.000
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.951
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.049
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -1.000
 > f_VHM_fun                             [-]: -0.1603061E-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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   11634.725
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      10.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    2020.349
   Hlug=Tlug*cos(theta_lug_star)        [kN]:   11457.968
   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]:    2020.349   11457.968 -103121.708
   Vstar,Hstar,Mstar                     [-]:       0.185       9.495      -7.121
   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.495      -7.121       0.185
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.950
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.049
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.999
 > f_VHM_fun                             [-]: -0.1913590E-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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   12085.967
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      20.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    4133.644
   Hlug=Tlug*cos(theta_lug_star)        [kN]:   11357.094
   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]:    4133.644   11357.094 -102213.849
   Vstar,Hstar,Mstar                     [-]:       0.379       9.412      -7.059
   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.412      -7.059       0.379
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.933
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.048
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.981
 > f_VHM_fun                             [-]: -0.7133960E-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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   12611.996
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      30.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    6305.998
   Hlug=Tlug*cos(theta_lug_star)        [kN]:   10922.309
   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]:    6305.998   10922.309  -98300.781
   Vstar,Hstar,Mstar                     [-]:       0.578       9.051      -6.788
   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.051      -6.788       0.578
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.863
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.045
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.908
 > f_VHM_fun                             [-]: -0.1820294E-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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   12844.185
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      40.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    8256.083
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    9839.216
   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]:    8256.083    9839.216  -88552.948
   Vstar,Hstar,Mstar                     [-]:       0.757       8.154      -6.115
   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.154      -6.115       0.757
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.700
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.036
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.736
 > f_VHM_fun                             [-]: -0.2253866E-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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   12581.954
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      50.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:    9638.336
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    8087.524
   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]:    9638.336    8087.524  -72787.715
   Vstar,Hstar,Mstar                     [-]:       0.884       6.702      -5.027
   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.702      -5.027       0.884
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.473
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.024
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.498
 > f_VHM_fun                             [-]: -0.1982954E-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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   12022.581
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      60.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:   10411.861
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    6011.291
   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]:   10411.861    6011.291  -54101.616
   Vstar,Hstar,Mstar                     [-]:       0.955       4.982      -3.736
   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.982      -3.736       0.955
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.261
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.013
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.275
 > f_VHM_fun                             [-]: -0.5631925E-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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   11455.585
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      70.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:   10764.728
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    3918.041
   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]:   10764.728    3918.041  -35262.366
   Vstar,Hstar,Mstar                     [-]:       0.987       3.247      -2.435
   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.247      -2.435       0.987
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.111
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.006
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.117
 > f_VHM_fun                             [-]:  0.1747576E-03
f_VHM_fun : END --- >

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


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   11051.916
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      80.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:   10884.013
   Hlug=Tlug*cos(theta_lug_star)        [kN]:    1919.145
   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]:   10884.013    1919.145  -17272.306
   Vstar,Hstar,Mstar                     [-]:       0.998       1.590      -1.193
   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.590      -1.193       0.998
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.027
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.001
   Z: - (1-(|Z|**bVH))**(2.0d0/aVH)                  [-]:      -0.029
 > f_VHM_fun                             [-]: -0.8927242E-03
f_VHM_fun : END --- >

 > root_Tlug_22_fun                      [-]:       -0.00089
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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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]:   10902.515
   L , D                                 [m]:      12.000       6.000
   zlug , theta_lug_star           [m] [deg]:       9.000      90.000
   zlug_cL                               [m]:       9.000
   Bw , Bw_plug                         [kN]:       0.000       0.000
   su_av_L                             [kPa]:      16.760
   Vmax                                 [kN]:   10902.515
   Hfac=L   D su_av_L                   [kN]:    1206.720
   Mfac=L^2 D su_av_L                  [kNm]:   14480.640
   aMH , bMH                             [-]:      12.173       1.626
   phiMH , phiMH                 [rad] [deg]:      -0.613     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   Vlug=Tlug*sin(theta_lug_star)        [kN]:   10902.515
   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]:   10902.515       0.000      -0.000
   Vstar,Hstar,Mstar                     [-]:       1.000       0.000      -0.000
   root_Tlug_22_fun=f_VHM_fun               

f_VHM_fun : < --- BEGIN
 < L,D                                   [m]:   6.00  12.00
   X = Hstar; Y = Mstar; Z = Vstar                  0.000      -0.000       1.000
   eta_aMH , eta_bMH , eta_phiMH         [-]:       1.000       1.000       1.000
   aMH,bMH,phiMH               [-],[-],[deg]:      12.173       1.626     -35.133
   aMH_new,bMH_new,phiMH_new   [-],[-],[deg]:      12.173       1.626     -35.133
   aVH , bVH                             [-]:       2.858       3.734
   XY:((+cos(phiMH)*X+sin(phiMH)*Y)/aMH)^2           [-]:       0.000
   XY:((-sin(phiMH)*X+cos(phiMH)*Y)/bMH)^2           [-]:       0.000
   Z:|Z| (=|Vstar|) > 1; -(1-|Z|)                    [-]:      -0.000
 > f_VHM_fun                             [-]:  0.3390294E-32
f_VHM_fun : END --- >

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


Vmax_fun : < --- BEGIN
 < support_foundation                  [T/F]: F
 < iplug_vmode                           [-]:  1
 < L,D                                   [m]:  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  1.000
   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]:    16.8    
   Fo = eta_Fo*pi*D *L*alpha_Fo*su_av_L [kN]:     2539.98520
   Fi =        pi*Di*L*alpha_Fi*su_av_L [kN]:     2506.71139
   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]:    35.2    
   qv_tip = Nc * su_av_tip             [kPa]     296.    
   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.836E+04
   Qwall =           Awall * qv_tip     [kN]:    0.00    
   Qxp   =           Axp   * qv_tip     [kN]:    0.00    
   Qplugged   = Fo+Qtip                 [kN]:    10902.51529
   Qunplugged = Fo+Fi+Fxp+Qwall+Qxp     [kN]:     5046.69659
   Qleaking   = Fo                      [kN]:     2539.98520
 > Vmax_fun                             [kN]:    10902.51529
Vmax_fun : END --- >


Bw_caisson_fun : < --- BEGIN
 < L,D                                   [m]:  12.00   6.00
 < iscope_Bw                             [-]:  0
 > Bw_caisson_fun                       [kN]:        0.00000
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 =  0
 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     10902.5     12013.8      3409.2         0.0        16.8
    2  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
    3  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
    4  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
    5  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
    6  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
    7  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
    8  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
    9  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
   10  1   6.00  12.00   2.00     10902.5     12013.8      3409.2         0.0        16.8
 ----------------------------------------------------------------------------------------
  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  11465.4     0.0       5  -0.160E-05   ! ? missing text
    2  1   6.00  12.00   9.00 10.00  11634.7     0.0       5  -0.191E-05   ! ? missing text
    3  1   6.00  12.00   9.00 20.00  12086.0     0.0       5  -0.713E-05   ! ? missing text
    4  1   6.00  12.00   9.00 30.00  12612.0     0.0       5  -0.182E-04   ! ? missing text
    5  1   6.00  12.00   9.00 40.00  12844.2     0.0       5  -0.225E-04   ! ? missing text
    6  1   6.00  12.00   9.00 50.00  12582.0     0.0       4  -0.198E-04   ! ? missing text
    7  1   6.00  12.00   9.00 60.00  12022.6     0.0       5  -0.563E-05   ! ? missing text
    8  1   6.00  12.00   9.00 70.00  11455.6     0.0       7   0.175E-03   ! ? missing text
    9  1   6.00  12.00   9.00 80.00  11051.9     0.0      11  -0.893E-03   ! ? missing text
   10  1   6.00  12.00   9.00 90.00  10902.5     0.0       3   0.339E-32   ! ? missing text
 -----------------------------------------------------------------------
 ilc        D      L    zlug  theta    Tlug     Bw    ( iters    f_VHM )
  icgeom   [m]    [m]   [m]   [deg]    [kN]    [kN]   ( [-]      [-]   )
 -----------------------------------------------------------------------

VHlug_resist_sub: < --- BEGIN
  load case #=    1; caisson geometry #=    1 (#rows results= =  181)
  L         [m] =   12.00
  D         [m] =    6.00
  zlug      [m] =    9.00
  iro_lug   [-] =   0: no offset (lug on centreline)
  m_fac     [-] =   1.000
  Hmax     [kN] =    0.120E+05
  Vmax     [kN] =    0.109E+05
  theta*        = resistance inclination at lug
  Vlug_theta*, Hlug_theta*
                = corresponding VH UNFACTORED resistances
  T_theta*      = sqrt( Vlug_theta*^2 + Hlug_theta*^2)
------------------------------------------------------------------------------
 theta*     Vlug_theta*  Hlug_theta*  T_theta*   Vlug_theta*  Hlug_theta*  row
  [deg]        [kN]         [kN]        [kN]      / Vmax [-]   / Hmax [-]  [-]
------------------------------------------------------------------------------
    90.0       0.109E+05   0.668E-12   0.109E+05       1.000    0.000        1
    89.0       0.109E+05    190.       0.109E+05       1.001    0.016        2
    88.0       0.109E+05    381.       0.109E+05       1.002    0.032        3
    87.0       0.109E+05    572.       0.109E+05       1.002    0.048        4
    86.0       0.109E+05    764.       0.109E+05       1.002    0.064        5
    85.0       0.109E+05    954.       0.109E+05       1.000    0.079        6
    84.0       0.109E+05   0.115E+04   0.110E+05       0.999    0.095        7
    83.0       0.109E+05   0.134E+04   0.110E+05       1.000    0.111        8
    82.0       0.109E+05   0.153E+04   0.110E+05       0.999    0.127        9
    81.0       0.109E+05   0.172E+04   0.110E+05       0.998    0.143       10
    80.0       0.109E+05   0.192E+04   0.111E+05       0.998    0.160       11
    79.0       0.109E+05   0.211E+04   0.111E+05       0.998    0.176       12
    78.0       0.109E+05   0.231E+04   0.111E+05       0.997    0.192       13
    77.0       0.109E+05   0.251E+04   0.111E+05       0.996    0.209       14
    76.0       0.109E+05   0.271E+04   0.112E+05       0.996    0.225       15
    75.0       0.108E+05   0.291E+04   0.112E+05       0.995    0.242       16
    74.0       0.108E+05   0.311E+04   0.113E+05       0.994    0.259       17
    73.0       0.108E+05   0.331E+04   0.113E+05       0.992    0.275       18
    72.0       0.108E+05   0.351E+04   0.114E+05       0.991    0.292       19
    71.0       0.108E+05   0.371E+04   0.114E+05       0.989    0.309       20
    70.0       0.108E+05   0.392E+04   0.115E+05       0.987    0.326       21
    69.0       0.107E+05   0.412E+04   0.115E+05       0.985    0.343       22
    68.0       0.107E+05   0.433E+04   0.116E+05       0.983    0.360       23
    67.0       0.107E+05   0.454E+04   0.116E+05       0.981    0.378       24
    66.0       0.107E+05   0.475E+04   0.117E+05       0.978    0.395       25
    65.0       0.106E+05   0.496E+04   0.117E+05       0.975    0.413       26
    64.0       0.106E+05   0.517E+04   0.118E+05       0.971    0.430       27
    63.0       0.106E+05   0.538E+04   0.118E+05       0.968    0.448       28
    62.0       0.105E+05   0.559E+04   0.119E+05       0.964    0.465       29
    61.0       0.105E+05   0.580E+04   0.120E+05       0.960    0.483       30
    60.0       0.104E+05   0.601E+04   0.120E+05       0.955    0.500       31
    59.0       0.104E+05   0.622E+04   0.121E+05       0.950    0.518       32
    58.0       0.103E+05   0.643E+04   0.121E+05       0.945    0.536       33
    57.0       0.102E+05   0.665E+04   0.122E+05       0.939    0.553       34
    56.0       0.102E+05   0.686E+04   0.123E+05       0.932    0.571       35
    55.0       0.101E+05   0.707E+04   0.123E+05       0.926    0.588       36
    54.0       0.100E+05   0.727E+04   0.124E+05       0.918    0.605       37
    53.0       0.993E+04   0.748E+04   0.124E+05       0.911    0.623       38
    52.0       0.984E+04   0.769E+04   0.125E+05       0.902    0.640       39
    51.0       0.974E+04   0.789E+04   0.125E+05       0.893    0.657       40
    50.0       0.964E+04   0.809E+04   0.126E+05       0.884    0.673       41
    49.0       0.953E+04   0.828E+04   0.126E+05       0.874    0.690       42
    48.0       0.941E+04   0.848E+04   0.127E+05       0.864    0.706       43
    47.0       0.929E+04   0.867E+04   0.127E+05       0.852    0.721       44
    46.0       0.916E+04   0.885E+04   0.127E+05       0.841    0.737       45
    45.0       0.903E+04   0.903E+04   0.128E+05       0.828    0.752       46
    44.0       0.889E+04   0.920E+04   0.128E+05       0.815    0.766       47
    43.0       0.874E+04   0.937E+04   0.128E+05       0.802    0.780       48
    42.0       0.858E+04   0.953E+04   0.128E+05       0.787    0.794       49
    41.0       0.842E+04   0.969E+04   0.128E+05       0.773    0.807       50
    40.0       0.826E+04   0.984E+04   0.128E+05       0.757    0.819       51
    39.0       0.808E+04   0.998E+04   0.128E+05       0.741    0.831       52
    38.0       0.790E+04   0.101E+05   0.128E+05       0.725    0.842       53
    37.0       0.772E+04   0.102E+05   0.128E+05       0.708    0.853       54
    36.0       0.753E+04   0.104E+05   0.128E+05       0.691    0.863       55
    35.0       0.733E+04   0.105E+05   0.128E+05       0.673    0.872       56
    34.0       0.714E+04   0.106E+05   0.128E+05       0.654    0.881       57
    33.0       0.693E+04   0.107E+05   0.127E+05       0.636    0.889       58
    32.0       0.673E+04   0.108E+05   0.127E+05       0.617    0.896       59
    31.0       0.652E+04   0.108E+05   0.127E+05       0.598    0.903       60
    30.0       0.631E+04   0.109E+05   0.126E+05       0.578    0.909       61
    29.0       0.609E+04   0.110E+05   0.126E+05       0.559    0.915       62
    28.0       0.588E+04   0.111E+05   0.125E+05       0.539    0.920       63
    27.0       0.566E+04   0.111E+05   0.125E+05       0.519    0.925       64
    26.0       0.544E+04   0.112E+05   0.124E+05       0.499    0.929       65
    25.0       0.522E+04   0.112E+05   0.124E+05       0.479    0.932       66
    24.0       0.501E+04   0.112E+05   0.123E+05       0.459    0.936       67
    23.0       0.479E+04   0.113E+05   0.123E+05       0.439    0.939       68
    22.0       0.457E+04   0.113E+05   0.122E+05       0.419    0.941       69
    21.0       0.435E+04   0.113E+05   0.121E+05       0.399    0.943       70
    20.0       0.413E+04   0.114E+05   0.121E+05       0.379    0.945       71
    19.0       0.392E+04   0.114E+05   0.120E+05       0.359    0.947       72
    18.0       0.370E+04   0.114E+05   0.120E+05       0.340    0.948       73
    17.0       0.349E+04   0.114E+05   0.119E+05       0.320    0.950       74
    16.0       0.327E+04   0.114E+05   0.119E+05       0.300    0.951       75
    15.0       0.306E+04   0.114E+05   0.118E+05       0.281    0.951       76
    14.0       0.285E+04   0.114E+05   0.118E+05       0.262    0.952       77
    13.0       0.264E+04   0.114E+05   0.117E+05       0.242    0.953       78
    12.0       0.243E+04   0.115E+05   0.117E+05       0.223    0.953       79
    11.0       0.223E+04   0.115E+05   0.117E+05       0.204    0.953       80
    10.0       0.202E+04   0.115E+05   0.116E+05       0.185    0.954       81
    9.00       0.182E+04   0.115E+05   0.116E+05       0.166    0.954       82
    8.00       0.161E+04   0.115E+05   0.116E+05       0.148    0.954       83
    7.00       0.141E+04   0.115E+05   0.115E+05       0.129    0.954       84
    6.00       0.120E+04   0.115E+05   0.115E+05       0.111    0.954       85
    5.00       0.100E+04   0.115E+05   0.115E+05       0.092    0.954       86
    4.00        802.       0.115E+05   0.115E+05       0.074    0.954       87
    3.00        601.       0.115E+05   0.115E+05       0.055    0.954       88
    2.00        400.       0.115E+05   0.115E+05       0.037    0.954       89
    1.00        200.       0.115E+05   0.115E+05       0.018    0.954       90
    0.00        0.00       0.115E+05   0.115E+05       0.000    0.954       91
   -1.00       -200.       0.115E+05   0.115E+05      -0.018    0.954       92
   -2.00       -400.       0.115E+05   0.115E+05      -0.037    0.954       93
   -3.00       -601.       0.115E+05   0.115E+05      -0.055    0.954       94
   -4.00       -802.       0.115E+05   0.115E+05      -0.074    0.954       95
   -5.00      -0.100E+04   0.115E+05   0.115E+05      -0.092    0.954       96
   -6.00      -0.120E+04   0.115E+05   0.115E+05      -0.111    0.954       97
   -7.00      -0.141E+04   0.115E+05   0.115E+05      -0.129    0.954       98
   -8.00      -0.161E+04   0.115E+05   0.116E+05      -0.148    0.954       99
   -9.00      -0.182E+04   0.115E+05   0.116E+05      -0.166    0.954      100
   -10.0      -0.202E+04   0.115E+05   0.116E+05      -0.185    0.954      101
   -11.0      -0.223E+04   0.115E+05   0.117E+05      -0.204    0.953      102
   -12.0      -0.243E+04   0.115E+05   0.117E+05      -0.223    0.953      103
   -13.0      -0.264E+04   0.114E+05   0.117E+05      -0.242    0.953      104
   -14.0      -0.285E+04   0.114E+05   0.118E+05      -0.262    0.952      105
   -15.0      -0.306E+04   0.114E+05   0.118E+05      -0.281    0.951      106
   -16.0      -0.327E+04   0.114E+05   0.119E+05      -0.300    0.951      107
   -17.0      -0.349E+04   0.114E+05   0.119E+05      -0.320    0.950      108
   -18.0      -0.370E+04   0.114E+05   0.120E+05      -0.340    0.948      109
   -19.0      -0.392E+04   0.114E+05   0.120E+05      -0.359    0.947      110
   -20.0      -0.413E+04   0.114E+05   0.121E+05      -0.379    0.945      111
   -21.0      -0.435E+04   0.113E+05   0.121E+05      -0.399    0.943      112
   -22.0      -0.457E+04   0.113E+05   0.122E+05      -0.419    0.941      113
   -23.0      -0.479E+04   0.113E+05   0.123E+05      -0.439    0.939      114
   -24.0      -0.501E+04   0.112E+05   0.123E+05      -0.459    0.936      115
   -25.0      -0.522E+04   0.112E+05   0.124E+05      -0.479    0.932      116
   -26.0      -0.544E+04   0.112E+05   0.124E+05      -0.499    0.929      117
   -27.0      -0.566E+04   0.111E+05   0.125E+05      -0.519    0.925      118
   -28.0      -0.588E+04   0.111E+05   0.125E+05      -0.539    0.920      119
   -29.0      -0.609E+04   0.110E+05   0.126E+05      -0.559    0.915      120
   -30.0      -0.631E+04   0.109E+05   0.126E+05      -0.578    0.909      121
   -31.0      -0.652E+04   0.108E+05   0.127E+05      -0.598    0.903      122
   -32.0      -0.673E+04   0.108E+05   0.127E+05      -0.617    0.896      123
   -33.0      -0.693E+04   0.107E+05   0.127E+05      -0.636    0.889      124
   -34.0      -0.714E+04   0.106E+05   0.128E+05      -0.654    0.881      125
   -35.0      -0.733E+04   0.105E+05   0.128E+05      -0.673    0.872      126
   -36.0      -0.753E+04   0.104E+05   0.128E+05      -0.691    0.863      127
   -37.0      -0.772E+04   0.102E+05   0.128E+05      -0.708    0.853      128
   -38.0      -0.790E+04   0.101E+05   0.128E+05      -0.725    0.842      129
   -39.0      -0.808E+04   0.998E+04   0.128E+05      -0.741    0.831      130
   -40.0      -0.826E+04   0.984E+04   0.128E+05      -0.757    0.819      131
   -41.0      -0.842E+04   0.969E+04   0.128E+05      -0.773    0.807      132
   -42.0      -0.858E+04   0.953E+04   0.128E+05      -0.787    0.794      133
   -43.0      -0.874E+04   0.937E+04   0.128E+05      -0.802    0.780      134
   -44.0      -0.889E+04   0.920E+04   0.128E+05      -0.815    0.766      135
   -45.0      -0.903E+04   0.903E+04   0.128E+05      -0.828    0.752      136
   -46.0      -0.916E+04   0.885E+04   0.127E+05      -0.841    0.737      137
   -47.0      -0.929E+04   0.867E+04   0.127E+05      -0.852    0.721      138
   -48.0      -0.941E+04   0.848E+04   0.127E+05      -0.864    0.706      139
   -49.0      -0.953E+04   0.828E+04   0.126E+05      -0.874    0.690      140
   -50.0      -0.964E+04   0.809E+04   0.126E+05      -0.884    0.673      141
   -51.0      -0.974E+04   0.789E+04   0.125E+05      -0.893    0.657      142
   -52.0      -0.984E+04   0.769E+04   0.125E+05      -0.902    0.640      143
   -53.0      -0.993E+04   0.748E+04   0.124E+05      -0.911    0.623      144
   -54.0      -0.100E+05   0.727E+04   0.124E+05      -0.918    0.605      145
   -55.0      -0.101E+05   0.707E+04   0.123E+05      -0.926    0.588      146
   -56.0      -0.102E+05   0.686E+04   0.123E+05      -0.932    0.571      147
   -57.0      -0.102E+05   0.665E+04   0.122E+05      -0.939    0.553      148
   -58.0      -0.103E+05   0.643E+04   0.121E+05      -0.945    0.536      149
   -59.0      -0.104E+05   0.622E+04   0.121E+05      -0.950    0.518      150
   -60.0      -0.104E+05   0.601E+04   0.120E+05      -0.955    0.500      151
   -61.0      -0.105E+05   0.580E+04   0.120E+05      -0.960    0.483      152
   -62.0      -0.105E+05   0.559E+04   0.119E+05      -0.964    0.465      153
   -63.0      -0.106E+05   0.538E+04   0.118E+05      -0.968    0.448      154
   -64.0      -0.106E+05   0.517E+04   0.118E+05      -0.971    0.430      155
   -65.0      -0.106E+05   0.496E+04   0.117E+05      -0.975    0.413      156
   -66.0      -0.107E+05   0.475E+04   0.117E+05      -0.978    0.395      157
   -67.0      -0.107E+05   0.454E+04   0.116E+05      -0.981    0.378      158
   -68.0      -0.107E+05   0.433E+04   0.116E+05      -0.983    0.360      159
   -69.0      -0.107E+05   0.412E+04   0.115E+05      -0.985    0.343      160
   -70.0      -0.108E+05   0.392E+04   0.115E+05      -0.987    0.326      161
   -71.0      -0.108E+05   0.371E+04   0.114E+05      -0.989    0.309      162
   -72.0      -0.108E+05   0.351E+04   0.114E+05      -0.991    0.292      163
   -73.0      -0.108E+05   0.331E+04   0.113E+05      -0.992    0.275      164
   -74.0      -0.108E+05   0.311E+04   0.113E+05      -0.994    0.259      165
   -75.0      -0.108E+05   0.291E+04   0.112E+05      -0.995    0.242      166
   -76.0      -0.109E+05   0.271E+04   0.112E+05      -0.996    0.225      167
   -77.0      -0.109E+05   0.251E+04   0.111E+05      -0.996    0.209      168
   -78.0      -0.109E+05   0.231E+04   0.111E+05      -0.997    0.192      169
   -79.0      -0.109E+05   0.211E+04   0.111E+05      -0.998    0.176      170
   -80.0      -0.109E+05   0.192E+04   0.111E+05      -0.998    0.160      171
   -81.0      -0.109E+05   0.172E+04   0.110E+05      -0.998    0.143      172
   -82.0      -0.109E+05   0.153E+04   0.110E+05      -0.999    0.127      173
   -83.0      -0.109E+05   0.134E+04   0.110E+05      -1.000    0.111      174
   -84.0      -0.109E+05   0.115E+04   0.110E+05      -0.999    0.095      175
   -85.0      -0.109E+05    954.       0.109E+05      -1.000    0.079      176
   -86.0      -0.109E+05    764.       0.109E+05      -1.002    0.064      177
   -87.0      -0.109E+05    572.       0.109E+05      -1.002    0.048      178
   -88.0      -0.109E+05    381.       0.109E+05      -1.002    0.032      179
   -89.0      -0.109E+05    190.       0.109E+05      -1.001    0.016      180
   -90.0      -0.109E+05   0.668E-12   0.109E+05      -1.000    0.000      181
------------------------------------------------------------------------------
 theta*     Vlug_theta*  Hlug_theta*  T_theta*   Vlug_theta*  Hlug_theta*  row
  [deg]        [kN]         [kN]        [kN]      / Vmax [-]   / Hmax [-]  [-]
------------------------------------------------------------------------------
VHlug_resist_sub: END --- >

solve_MODE22_sub: END --- >


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