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)