Maple理论力学 作者 李银山 第五部分B 第22章.pptVIP

● Maple程序 v[A]:=Dx: #小车的速度。 v[e]:=v[A]: #重物B的牵连速度 。 v[r]:=l*Dphi: #重物B的相对速度 。 v[a]:=sqrt(v[e]^2+v[r]^2+2*v[e]*v[r]*cos(phi)): #重物B的绝对速度 。 v[B]:=v[a]: #重物B的速度。 T:=1/2*m[1]*v[A]^2+1/2*m[2]*v[B]^2: #系统动能。 V:=1/2*k*x^2+m[2]*g*l*(1-cos(phi)): #系统势能。 L:=T-V: #拉格朗日函数。 Lx:=diff(L,x): # 。 Lx:=subs(x=x(t),Lx): #代换。 LDx:=diff(L,Dx): # 。 LDx:=subs(Dx=Dx(t),Dphi=Dphi(t),phi=phi(t),LDx):#代换。 eq1:=diff(LDx,t)-Lx=0: #由拉格朗日方程 。 eq1:=subs(diff(Dx(t),t)=DDx,diff(Dphi(t),t)=DDphi, diff(phi(t),t)=Dphi,Dphi(t)=Dphi,x(t)=x,phi(t)=phi,eq1): # 。 eq1:=simplify(eq1): #化简。 eq1:=collect(eq1,DDx): #多项式排列。 Lphi:=diff(L,phi): # 。 Lphi:=subs(Dx=Dx(t),Dphi=Dphi(t),phi=phi(t),Lphi): #代换。 LDphi:=diff(L,Dphi): # 。 LDphi:=subs(Dx=Dx(t),Dphi=Dphi(t),phi=phi(t),LDphi):#代换。 eq2:=diff(LDphi,t)-Lphi=0: #由拉格朗日方程 。 eq2:=subs(diff(Dx(t),t)=DDx,diff(Dphi(t),t)=DDphi, diff(phi(t),t)=Dphi,Dx(t)=Dx,Dphi(t)=Dphi,x(t)=x, phi(t)=phi,eq2): #代换。 eq2:=simplify(eq2/(m[2]*l)): #化简。 eq2:=collect(eq2,DDphi): #多项式排列。 eq1:=subs(sin(phi)=phi,cos(phi)=1,Dphi^2=0,eq1): #线性化。 eq2:=subs(sin(phi)=phi,cos(phi)=1,eq2): #线性化。 x:=A*sin(omega*t+theta): #设方程解的形式。 phi:=B*sin(omega*t+theta): #设方程解的形式。 DDx:=diff(x,t$2): # 对时间 求二阶导数。 DDphi:=diff(phi,t$2): # 对时间 求二阶导数。 eq3:=simplify(eq1/sin(omega*t+theta)): #化简。 eq4:=simplify(eq2/sin(omega*t+theta)): #化简。 eq3:=subs(B=A*lambda,eq3): #代换。 eq4:=subs(B=A*lambda,eq4): #代换。 eq3:=expand(eq3/A): #展开。 eq4:=expand(eq4/A): #展开。 m[1]:=2220: m[2]:=2040: g:=9.8: #已知条件。 k:=852.6*10^3: l:=14: v0:=1: #已知条件。 solve({eq3,eq4},{omega,lambda});#解方程组。 omega1:=.8359589802: #第一固有频率。 omega2:=19