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

F[C]:=k*delta[0]: #弹簧对点 的拉力 的大小。 F[G]:=F[C]: #弹簧对点 的拉力 的大小。 eq2:=F[Bx]*delta[xB]+F[C]*delta[yC]-F[G]*delta[yG]+F*delta[yG]=0: # 。 solve({eq2},{F[Bx]}); #解方程。 例15-3 图15-6所示椭圆规机构中，连杆 长为 ，滑块 ， 与杆重均不计，忽略各处摩擦，机构在图示位置平衡时，主动力 和 之间的关系。 解法一:●建模 用虚位移原理求解。椭圆规机构平衡系统受主动力： , 。 ● Maple程序 restart: #清零。 eq1:=F[A]*delta[rA]-F[B]*delta[rB]=0: # 。 eq2:=delta[rB]*cos(phi)=delta[rA]*sin(phi):#虚位移之间的关系。 solve({eq1,eq2},{delta[rA],F[A]}); #解方程组。 F[A]:=F[B]*sin(phi)/cos(phi): #主动力 的大小。 F[A]:=convert(F[A],tan); #将三角函数转换成tan形式。 解法二: ●建模 用解析法求解。椭圆规机构平衡系统受主动力： , 。 ● Maple程序 restart: #清零。 x[B]:=l*cos(phi): #点 的 坐标。 y[A]:=l*sin(phi): #点 的 坐标。 delta[xB]:=diff(x[B],phi)*delta[phi]: #对点 的 坐标的变分。 delta[yA]:=diff(y[A],phi)*delta[phi]: #对点 的 坐标的变分。 eq:=-F[A]*delta[yA]-F[B]*delta[xB]=0: # 。 solve({eq},{F[A]}); #解方程。 F[A]:=F[B]*sin(phi)/cos(phi): #主动力 的大小。 F[A]:=convert(F[A],tan); #将三角函数转换成tan形式。 解法三: ●建模 用虚速度法求解。椭圆规机构平衡系统受主动力： , 。 ● Maple程序 restart: #清零。 eq1:=F[A]*v[A]-F[B]*v[B]=0: #虚功率方程。 eq2:=v[B]*cos(phi)=v[A]*sin(phi): #速度投影定理。 solve({eq1,eq2},{F[A],v[A]}); #解方程组。 F[A]:=F[B]*sin(phi)/cos(phi): #主动力 的大小。 F[A]:=convert(F[A],tan); #将三角函数转换成tan形式。 例15-4 如图15-6所示机构，不计各构件自重与各处摩擦，求机构图示位置平衡时，主动力偶矩 与主动力 之间的关系。 解法一:●建模 用虚位移原理求解。平衡系统受主动力： , 。 ● Maple程序 restart: #清零。 OB:=h/sin(theta): # 几何关系。 delta[re]:=OB*delta[theta]: # 牵连虚位移。 delta[ra]:=delta[re]/sin(theta): # 绝对虚位移。 delta[rC]:=delta[ra]: # 点 的虚位移。 eq:=M*delta[theta]-F*de