def test_nzs_vm4_load_case_5():
"""
values from: NZ Building code Clause B1 VM4 example in Appendix C
-Retain wall example load case 5
***looks like there is an error in the calculation of d_c
***small discrepancy between Nc, i_c and i_q
"""
length = 10000 # Actually should be a strip
width = 2.65
depth = 0.4
phi = 25.0
cohesion = 12.5 # kPa
unit_dry_weight = 8.2 # kN/m3
sl = geofound.create_soil(phi, cohesion, unit_dry_weight)
fd = geofound.create_foundation(length, width, depth)
h_b = 70.45 * length
vertical_load = 144.48 * length
h_eff_b = 1.44
loc_v_b = 0.813
geofound.capacity_nzs_vm4_2011(sl,
fd,
h_b=h_b,
vertical_load=vertical_load,
h_eff_b=h_eff_b,
loc_v_b=loc_v_b,
verbose=0)
assert np.isclose(fd.nc_factor, 20.72, rtol=0.001)
assert np.isclose(fd.nq_factor, 10.66, rtol=0.01)
assert np.isclose(fd.ng_factor, 9.01, rtol=0.01)
assert np.isclose(fd.q_ult, 145.02, rtol=0.001) # 152.70?
def test_nzs_vm4_load_case_3():
"""
values from: NZ Building code Clause B1 VM4 example in Appendix C
-Retain wall example load case 3
"""
length = 10000 # Actually should be a strip
width = 2.65
depth = 0.4
phi = 0
cohesion = 75.0 # kPa
unit_dry_weight = 18 # kN/m3
sl = geofound.create_soil(phi, cohesion, unit_dry_weight)
fd = geofound.create_foundation(length, width, depth)
h_b = 93.27 * length
vertical_load = 154.87 * length
h_eff_b = 1.78
loc_v_b = 0.854
geofound.capacity_nzs_vm4_2011(sl,
fd,
h_b=h_b,
vertical_load=vertical_load,
h_eff_b=h_eff_b,
loc_v_b=loc_v_b,
verbose=0)
assert np.isclose(fd.nc_factor, 5.14, rtol=0.001)
assert np.isclose(fd.nq_factor, 1.0, rtol=0.01)
assert np.isclose(fd.ng_factor, 0.0, rtol=0.01)
assert np.isclose(fd.q_ult, 301.68, rtol=0.001)