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

例6-3 如图6-10所示，当液压减震器工作时，它的活塞在套筒内作直线往复运动。设活塞的加速度 ( 为活塞的速度, 为比例常数)，初速为 ，求活塞的运动规律。 解: ●建模 活塞作直线运动。取坐标轴 如图所示。 ● Maple程序 restart: #清零。 a:=-k*v: #已知条件。 v:=diff(x(t),t): #活塞的速度。 eq:=diff(x(t),t$2)=a: #活塞的加速度。 dsolve({eq,D(x)(0)=v0,x(0)=x0},x(t));#求解二阶常微分方程。 例6-4 列车沿半径为 的圆弧轨道作匀加速运动。如初速度为零,经过 后,速度达到 。求列车在起点和末点的加速度。 解：●建模 由于列车沿圆弧轨道作匀加速运动，切向加速度 等于恒量。在起点，初速度为零，因此法向加速度等于零，列车只有切向加速度。在末点时，速度不等于零，既有切向加速度，又有法向加速度。 ● Maple程序 restart: #清零。 a[0]:=a[t]: #在起点，列车的加速度。 a[1]:=sqrt(a[t]^2+a[n]^2): #在末点，列车的加速度。 theta:=arctan(a[t]/a[n]): #在末点，列车加速度与法向的夹角。 a[t]:=v[1]/tau: #切向加速度。 a[n]:=v[1]^2/R: #在末点的法向加速度。 v[1]:=54*1000/3600: v[0]:=0: tau:=2*60: R:=800: #已知条件。 a[0]:=evalf(a[0],4); #在起点列车加速度的数值。 a[1]:=evalf(a[1],4); #在末点列车加速度的数值。 theta:=evalf(theta*180/Pi,4); #在末点列车加速度与法向的夹角的数值。 theta0:=trunc(theta); #夹角的度。 theta1:=trunc(frac(theta)*60+0.5); #夹角的分。 例6-5 已知点的运动方程为 ， ， 。求点运动轨迹的曲率半径 。 解：●建模 点运动轨迹的曲率半径可以根据法向加速度的大小公式 求出。其中 ， ， ， 。 ● Maple程序 restart: #清零。 with(linalg): #加载向量运算函数库。 x:=2*sin(4*t): y:=2*cos(4*t): z:=4*t: #已知条件。 V:=vector([diff(x,t),diff(y,t),diff(z,t)]): #速度向量。 A:=vector([diff(x,t$2),diff(y,t$2),diff(z,t$2)]): #加速度向量。 v:=norm(V,2): #速度向量的模。 a:=norm(A,2): #加速度向量的模。 v:=simplify(v,symbolic): #化简根号，不管根号里变量的正负问题。 a:=simplify(a,symbolic): #化简根号，不管根号里变量的正负问题。 at:=diff(v,t): #切向加速度大小。 an:=sqrt(a^2-at^2): #法向加速度大小。 rho:=v^2/an: #运动轨迹的曲率半径 。 rho:=evalf(rho,4); #曲率半径 的数值。 例6-6半径为 的轮子沿直线轨道无滑动地滚动(称为纯滚动)，设轮子转角 为常值)如图6-11a所示，求用直角坐标和弧坐标表示的轮缘上任一点 的运动方程，并求该点速度