def capacity_terzaghi_1943(sl, fd, round_footing=False, verbose=<NUM_LIT:0>, **kwargs):
if not kwargs.get("<STR_LIT>", False):<EOL><INDENT>models.check_required(sl, ["<STR_LIT>", "<STR_LIT>", "<STR_LIT>"])<EOL>models.check_required(fd, ["<STR_LIT>", "<STR_LIT:width>", "<STR_LIT>"])<EOL><DEDENT>a02 = ((np.exp(np.pi * (<NUM_LIT> - sl.phi / <NUM_LIT>) * np.tan(sl.phi_r))) ** <NUM_LIT:2>)<EOL>a0_check = (np.exp((<NUM_LIT> - sl.phi) / <NUM_LIT> * np.pi * np.tan(sl.phi_r)))<EOL>if (a02 - a0_check) / a02 > <NUM_LIT>:<EOL><INDENT>raise DesignError<EOL><DEDENT>fd.nq_factor = (a02 / (<NUM_LIT:2> * (np.cos((<NUM_LIT> + sl.phi / <NUM_LIT:2>) * np.pi / <NUM_LIT>)) ** <NUM_LIT:2>))<EOL>fd.ng_factor = (<NUM_LIT:2> * (fd.nq_factor + <NUM_LIT:1>) * np.tan(sl.phi_r) / (<NUM_LIT:1> + <NUM_LIT> * np.sin(<NUM_LIT:4> * sl.phi_r)))<EOL>if sl.phi_r == <NUM_LIT:0>:<EOL><INDENT>fd.nc_factor = <NUM_LIT><EOL><DEDENT>else:<EOL><INDENT>fd.nc_factor = (fd.nq_factor - <NUM_LIT:1>) / np.tan(sl.phi_r)<EOL><DEDENT>if round_footing:<EOL><INDENT>s_c = <NUM_LIT><EOL>s_g = <NUM_LIT><EOL><DEDENT>elif fd.length / fd.width < <NUM_LIT:5>:<EOL><INDENT>s_c = <NUM_LIT><EOL>s_g = <NUM_LIT><EOL><DEDENT>else:<EOL><INDENT>s_c = <NUM_LIT:1.0><EOL>s_g = <NUM_LIT:1.0><EOL><DEDENT>s_q = <NUM_LIT:1.0><EOL>q_d = sl.unit_dry_weight * fd.depth<EOL>fd.q_ult = (sl.cohesion * fd.nc_factor * s_c + q_d * fd.nq_factor * s_q + <NUM_LIT:0.5> * fd.width *<EOL>sl.unit_dry_weight * fd.ng_factor * s_g)<EOL>if verbose:<EOL><INDENT>log("<STR_LIT>", fd.nc_factor)<EOL>log("<STR_LIT>", fd.nq_factor)<EOL>log("<STR_LIT>", fd.ng_factor)<EOL>log("<STR_LIT>", s_c)<EOL>log("<STR_LIT>", s_q)<EOL>log("<STR_LIT>", s_g)<EOL>log("<STR_LIT>", fd.q_ult)<EOL><DEDENT>return fd.q_ult<EOL>
Calculates the foundation capacity according Terzaghi (1943)
Ref: http://geo.cv.nctu.edu.tw/foundation/
download/BearingCapacityOfFoundations.pdf

:param sl: Soil object
:param fd: Foundation object
:param round_footing: if True, then foundation is round
:param verbose: verbosity
:return: ultimate bearing stress
Note: the shape factor of 1.3 is used for aspect ratio > 6