fenda88 / geofound Goto Github PK

import geofound
import geofound.settlement
import geofound.capacity


length = 21.7
width = 1.5  # metres
depth = 1.5  # metres
phi = 38
cohesion = 0
unit_weight = 18.5  # submerged
youngs_modulus_soil = 30000


fd = geofound.create_foundation(length, width, depth)
sl = geofound.create_soil(phi, cohesion, unit_weight)

sl.unit_sat_weight = 18.5

q_lim = geofound.capacity.capacity_vesics_1975(sl, fd)
p_max = q_lim * length * width

print(' ')
print('Ultimate bearing stress is q_lim = ' + str(round(q_lim,0)) + ' kPa')
print('Ultimate load is Q_lim = ' + str(round(p_max, 0)) + ' kN')

load = 21484
s = geofound.settlement.settlement_schmertmann(sl, fd, load, youngs_modulus_soil)
print(' ')
print('Settlement is si = ' + str(round(s,2)) + ' m')