def test_from_encn452_2013():
"""
Values from HW#7 crib ENCN452 course 2013
"""
length = 6.0 # actually a strip in
width = 3.0
depth = 1.5
phi = 0.0
cohesion = 40.0
unit_dry_weight = 18.0
sl = geofound.create_soil(phi, cohesion, unit_dry_weight)
fd = geofound.create_foundation(length, width, depth)
geofound.capacity_terzaghi_1943(sl, fd, verbose=0)
geofound.capacity_vesics_1975(sl, fd, verbose=0)
assert np.isclose(fd.q_ult, 298.0, rtol=0.001)
def test_terzaghi_again():
"""
values from: Gunaratne, Manjriker. 2006. "Spread Footings: Analysis and Design."
"""
length = 2
width = 2
depth = 2
phi = 15
cohesion = 20.0 # kPa
unit_dry_weight = 17 # kN/m3
sl = geofound.create_soil(phi, cohesion, unit_dry_weight)
fd = geofound.create_foundation(length, width, depth)
geofound.capacity_terzaghi_1943(sl, fd, verbose=0)
print(fd.q_ult)
assert np.isclose(fd.nc_factor, 12.86, rtol=0.001)
assert np.isclose(fd.nq_factor, 4.45, rtol=0.001)
assert np.isclose(fd.ng_factor, 2.168, rtol=0.001)
assert np.isclose(fd.q_ult, 515.0, rtol=0.001)
def test_terzaghi():
"""
values from: Gunaratne, Manjriker. 2006. "Spread Footings: Analysis and Design."
- except qult was not validated
"""
length = 2
width = 2
depth = 1
phi = 35
cohesion = 0
unit_dry_weight = 17
sl = geofound.create_soil(phi, cohesion, unit_dry_weight)
fd = geofound.create_foundation(length, width, depth)
geofound.capacity_terzaghi_1943(sl, fd, verbose=0)
assert np.isclose(fd.nc_factor, 57.8, rtol=0.001)
assert np.isclose(fd.nq_factor, 41.4, rtol=0.001)
assert np.isclose(fd.ng_factor, 47.3, rtol=0.001)
assert np.isclose(fd.q_ult, 1347.0, rtol=0.001)
print("DONE")