西南交大路基课设重力式挡土墙

⑵剪应力检算，检算合格3.4利用失算法对计算结果进行校核作用在挡土墙墙背的主动土压力，一般按库仑主动土压力公式计算。当破裂面交于路基面时，破裂棱体的面积S随着挡土墙及破裂面位置而变化，但都可归纳为一个表达式：式中，当边界条件确定后，A0、B0为常数，并可从破裂棱体的几何关系求得。附表“各种边界条件下的库仑主动土压力公式”给出了不同边界条件下的库仑主动土压力计算公式。在具体计算时，由于无法预知破裂面的位置，一般是先假设破裂面位置，然后按此情况计算出破裂角θ，再根据几何关系来校核假设是否正确。若假设不合理，则需选用另外的破裂面位置重新计算，直至校核合理。最后可根据附表中公式计算土压力的大小，方向和作用点位置。然后将此结果与穷举方法所得结果进行对比校核。

附录代码如下：ifyi<m1N=B*H*23+Eay;Kc=0.5*(N+Eax*tan(j10))/(Eax-N*sin(j10));endifyi<m2&&yi>m1N=B*H*23+Eay+((H+a)*tan(j6max)-b-k-H*tan(j1))*h0*e;Kc=0.5*(N+Eax*tan(j10))/(Eax-N*sin(j10));endifyi<m3&&yi>m2N=B*H*23+Eay+l0*h0*e;Kc=0.5*(N+Eax*tan(j10))/(Eax-N*sin(j10));endifyi<m4&&yi>m3N=B*H*23+Eay+l0*h0*e+((H+a)*tan(j6max)-b-k-l0-D-H*tan(j1))*h0*e;Kc=0.5*(N+Eax*tan(j10))/(Eax-N*sin(j10));endifyi>m4N=B*H*23+2*l0*h0*e;Kc=0.5*(N+Eax*tan(j10))/(Eax-N*sin(j10));endK0=(B*H*23*(H*tan(j1)+B)/2+Eay*Zymax)/(Eax*Zxmax);Zn=((B*H*23*(H*tan(j1)+B)/2+Eay*Zymax)-(Eax*Zxmax))/(B*H*23+Eay);e0=B/2-Zn;delta1=(B*H*23+Eay+l0*h0*e)/B*(1+6*e0/B);delta2=(B*H*23+Eay+l0*h0*e)/B*(1-6*e0/B);ifdelta1>delta2deltamax=delta1;elsedeltamax=delta2;endfprintf('滑动稳定性系数：%f;\n倾覆稳定性系数：%f\n\n',Kc,K0);ifKc>1.3&&Kc<1.5fprintf('滑动稳定性检算通过\n');elsefprintf('滑动稳定性检算失败\n');endifK0>1.6&&K0<1.8fprintf('倾覆稳定性检算通过\n');elsefprintf('倾覆稳定性检算失败\n');endife0<=B/6fprintf('基底偏心距检算通过\n');ifdeltamax<400fprintf('基底应力检算通过\n');elsefprintf('基底应力检算失败\n');endelsefprintf('基底偏心距检算失败\n');fprintf('基底应力偏心距检算失败\n');endifdeltamax<=600fprintf('墙体法向应力检算通过\n');elsefprintf('墙体法向应力检算失败\n');endifEax/B<=600fprintf('墙体剪应力检算通过\n');fprintf('====================================================\n');elsefprintf('墙体剪应力检算失败\n');fprintf('====================================================\n');endfprintf('检算完成\n\n'