Maple理论力学 作者 李银山 第一部分 第5章.pptVIP

● Maple程序 restart: #清零： eq1:=F[T]-P=0: #重物， 。 eq2:=F[T]*r-F[s]*R=0: #鼓轮， 。 eq3:=F[min]*a+F[s]*c-F[N]*b=0: #杠杆 ， 。 eq4:=F[s]=f[s]*F[N]: #补充方程。 solve({eq1,eq2,eq3,eq4},{F[T],F[s],F[N],F[min]}); #求解方程组。 例 5-4 图5-10所示的匀质木箱重 ，它与地面间的静摩擦因数 。图中 ， 。求：①当D处的拉力 时，木箱是否平衡? ②能保持木箱平衡的最大拉力。 解：●建模 欲保持木箱平衡，必须满足两个条件：一是不发生滑动，即要求静摩擦力 ；二是不绕 点翻倒，这时法向约束力 的作用线应在木箱内，即 。平衡时，木箱受力： 、 、 、 。 为求保持平衡的最大拉力 ，可分别求出木箱将滑动时的临界拉力 和木箱将绕 点翻倒时的临界拉力 。二者中取其较小者，即为所求。将滑动时，木箱受力： 、 、 、 。将翻倒时，木箱受力： 、 、 、 。 ● Maple程序 restart: #清零： eq1:=F[s]-f*cos(theta)=0: #木箱平衡时， 。 eq2:=F[N]-P+f*sin(theta)=0: #木箱平衡时， 。 eq3:=h*f*cos(theta)-P*a/2+F[N]*d=0: #木箱平衡时， 。 solve({eq1,eq2,eq3},{F[s],F[N],d}); #求解方程组。 eq4:=F1[s]-F[slide]*cos(theta)=0: #木箱将滑动时， 。 eq5:=F1[N]-P+F[slide]*sin(theta)=0: #木箱将滑动时， 。 eq6:=F1[s]-fs*F1[N]=0: #木箱将滑动时，补充方程。 solve({eq4,eq5,eq6},{F1[s],F1[N],F[slide]}); #求解方程组。 eq7:=h*F[overturn]*cos(theta)-P*a/2=0: solve({eq7},{F[overturn]}); #求解方程。 #木箱将翻倒时， 。 F[s]:=f*cos(theta): #静摩擦力 的大小。 F[N]:=P-f*sin(theta): #法向约束力 的大小。 d:=1/2*(2*h*f*cos(theta)-P*a)/(-P+f*sin(theta)):# 到 点的距离 。 F[slide]:=fs*P/(fs*sin(theta)+cos(theta)): #滑动时的临界拉力 。 F[overturn]:=1/2*P*a/h/cos(theta): #翻倒时的临界拉力 。 F[smax]:=fs*F[N]: #最大静摩擦力 的大小。 F[max]:=min(F[slide],F[overturn]): # , 二者中取其较小者。 P:=5*10^3: f:=1*10^3: #已知条件。 a:=1: h:=2*a: #已知条件。 theta:=Pi/6: fs:=0.4: #已知条件。 F[s]:=evalf(F[s],4); #静摩擦力