圆周运动、一般曲线运动、阿基米德螺旋线
乔宕一
2021年08月01日 18:21
收录于文集
共1篇

恒长旋转向量的导数

一个恒长旋转向量求导后得到的向量的方向与原向量相比,逆时针旋转了  ,而求导后得到的向量的长度与旋转角速度有关。

证明

例如 

1、对  求导,

%5Cfrac%7Bd%5Cvec%7Ba%7D%7D%7Bd%5Ctheta%7D%3D(-sin%5C%20%5Ctheta%2C%20%5Cquad%20cos%20%5C%20%5Ctheta)%3D%5Bcos(%5Ctheta%2B%5Cfrac%7B%5Cpi%7D%7B2%7D)%2C%20%5Cquad%20sin(%5Ctheta%2B%5Cfrac%7B%5Cpi%7D%7B2%7D)%5D

结论:

 与  相比,在方向上逆时针旋转了 90°,长度不变


2、对 t 求导

 是关于 t 的函数, 或 前者角速度恒定,后者角速度是关于时间的函数。

%5Cfrac%7Bd%5Cvec%7Ba%7D%7D%7Bdt%7D%3D%5Cfrac%7Bd%5Cvec%7Ba%7D%7D%7Bd%5Ctheta%7D%20%5Ccdot%20%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D%3D%5B%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D%20%5Ccdot%20cos(%5Ctheta%2B%5Cfrac%7B%5Cpi%7D%7B2%7D)%2C%20%5Cquad%20%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D%20%5Ccdot%20sin(%5Ctheta%2B%5Cfrac%7B%5Cpi%7D%7B2%7D)%5D

结论:

%5Cfrac%7Bd%5Cvec%7Ba%7D%7D%7Bdt%7D 与 %5Cvec%7Ba%7D 相比,在方向上逆时针旋转了 90° ,在长度上变为原来的 w%3D%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D 倍。如果 %5Cvec%7Ba%7D 作变速圆周运动,%5Cfrac%7Bd%5Cvec%7Ba%7D%7D%7Bdt%7D 的长度将会随着时间变化。


圆周运动中角速度和线速度的关系

一个质点以原点为圆心作圆周运动,设它的位移由  变成 ,如下图所示

 是位移的变化量, , 是弧度增量。

 很小时,  根据  ,可得

使用标量表示:

 ,方程两边同时除以  ,得  ,则得到  ,

这只是大小关系,考虑方向,使用叉乘, 

阿基米德螺旋线

阿基米德螺旋线方程

阿基米德螺旋线的极坐标方程为: 这表示极径与  是线性关系,成正比。

它的参数方程为:

%5Cbegin%7Bcases%7D%0Ax%3Dr%20%5Ccdot%20cos%20%5C%20%5Ctheta%3Da%5Ctheta%20%5Ccdot%20cos%20%5C%20%5Ctheta%20%5C%5C%0Ay%3Dr%20%5Ccdot%20sin%20%5C%20%5Ctheta%3Da%5Ctheta%20%5Ccdot%20sin%20%5C%20%5Ctheta%0A%5Cend%7Bcases%7D

试图消去参数

x%5E2%2By%5E2%3D(a%5Ctheta)%5E2%20%5Ccdot%20(sin%5E2%20%5Ctheta%20%2B%20cos%5E2%20%5Ctheta)%3D(a%5Ctheta)%5E2%20%5Cquad%20%E5%85%B6%E4%B8%AD%EF%BC%8C%5Ctheta%3Darctan(%5Cfrac%7By%7D%7Bx%7D)

这个方程无法化为显函数的形式,所以,最好用极坐标方程或参数方程来表示。我们以  为参数,用 matlab 画出它的图像。

代码块
clike
自动换行
复制代码
syms theta
a=1;
r=a*theta;
x=r*cos(theta);
y=r*sin(theta);
fplot(x,y,[0,50],'LineWidth',1.5);
grid on
axis square
复制成功

螺旋线的时间导数

 为螺旋线的位移, ,可见, 对时间求导只是对  求导后乘上一个系数  而已,所以我们研究对  的导数,而不是对 t 的导数。


对螺旋线的参数方程求导

代码块
clike
自动换行
复制代码
syms a theta
x=a*theta*cos(theta);
y=a*theta*sin(theta);
dx=diff(x,theta)
dy=diff(y,theta)
复制成功

结果如下

代码块
clike
自动换行
复制代码
 
dx =
 
a*cos(theta) - a*theta*sin(theta)
 
 
dy =
 
a*sin(theta) + a*theta*cos(theta)
 
>> 
复制成功

这可以看成是两个向量的合成,分别为 

(a%5Ccdot%20cos%5C%20%5Ctheta%2C%5Cquad%20a%20%5Ccdot%20sin%5C%20%5Ctheta) 和 

(-a%5Ctheta%20%5Ccdot%20sin%5C%20%5Ctheta%2C%5Cquad%20a%5Ctheta%20%5Ccdot%20cos%20%5C%20%5Ctheta)%3D%20%5Ba%5Ctheta%20%5Ccdot%20cos(%5Ctheta%2B%20%5Cfrac%7B%5Cpi%7D%7B2%7D)%2C%20%5Cquad%20a%5Ctheta%20%5Ccdot%20sin(%5Ctheta%2B%5Cfrac%7B%5Cpi%7D%7B2%7D)%5D

分别令

%5Cbegin%7Bcases%7D%0A%5Cvec%7Be_t%7D%3D%5Ba%5Ctheta%20%5Ccdot%20cos(%5Ctheta%2B%20%5Cfrac%7B%5Cpi%7D%7B2%7D)%2C%20%5Cquad%20a%5Ctheta%20%5Ccdot%20sin(%5Ctheta%2B%5Cfrac%7B%5Cpi%7D%7B2%7D)%5D%20%5C%5C%0A%5Cvec%7Be_n%7D%3D(a%5Ccdot%20cos%5C%20%5Ctheta%2C%5Cquad%20a%20%5Ccdot%20sin%5C%20%5Ctheta)%20%5C%5C%0A%5Cend%7Bcases%7D

%5Cbegin%7Bcases%7D%0A%5Cvec%7Bv_t%7D%3D%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D%20%5Ccdot%20%5Cvec%7Be_t%7D%3Dw%5Ccdot%20%5Cvec%7Be_t%7D%20%5C%5C%0A%5Cvec%7Bv_n%7D%3D%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D%20%5Ccdot%20%5Cvec%7Be_n%7D%20%3Dw%5Ccdot%20%5Cvec%7Be_n%7D%5C%5C%0A%0A%5Cend%7Bcases%7D

其中, 的方向沿着半径向外,是法向速度,  沿着切线,与  垂直,总是指向逆时针方向,是切向速度。

可以看出  恒定时,切向速度随着  的增大而增大,法向速度恒定,此时质点一方面绕着原点作圆周运动,线速度越来越大;另一方面,又以恒定的速度远离原点。


切向速度的规律

这个切向速度与原来的螺旋线方程比起来有什么规律呢?

1. 螺旋线参数方程

%5Cbegin%7Bcases%7D%0Ax%3Dr%20%5Ccdot%20cos%20%5C%20%5Ctheta%3Da%5Ctheta%20%5Ccdot%20cos%20%5C%20%5Ctheta%20%5C%5C%0Ay%3Dr%20%5Ccdot%20sin%20%5C%20%5Ctheta%3Da%5Ctheta%20%5Ccdot%20sin%20%5C%20%5Ctheta%0A%5Cend%7Bcases%7D

2. 切向速度

%5Cvec%7Bv_t%7D%3D%5Ba%5Ctheta%20%5Ccdot%20cos(%5Ctheta%2B%20%5Cfrac%7B%5Cpi%7D%7B2%7D)%2C%20%5Cquad%20a%5Ctheta%20%5Ccdot%20sin(%5Ctheta%2B%5Cfrac%7B%5Cpi%7D%7B2%7D)%5D

可以发现,切向速度与螺旋线方程比起来就是逆时针旋转了  而已。

使用 matlab 画出切向速度的图像

代码块
clike
自动换行
复制代码
clc;
clear;
close all;

syms theta
a=1;
x=a*theta*cos(theta+pi/2);
y=a*theta*sin(theta+pi/2);
fplot(x,y,[0,pi*2*5],'LineWidth',1.5,'Color','b');%画出切向速度,蓝色
grid on
axis square
hold on
x=a*theta*cos(theta);
y=a*theta*sin(theta);
fplot(x,y,[0,pi*2*5],'LineWidth',1.5,'Color','r');%画出阿基米德螺旋线,红色
复制成功


切向速度可以使用   的公式。(合速度不行)

1、使用方程组得到  ( 的标量,即它的大小)

%5Cbegin%7Bcases%7D%0A%5Cvec%7Be_t%7D%3D%5Ba%5Ctheta%20%5Ccdot%20cos(%5Ctheta%2B%20%5Cfrac%7B%5Cpi%7D%7B2%7D)%2C%20%5Cquad%20a%5Ctheta%20%5Ccdot%20sin(%5Ctheta%2B%5Cfrac%7B%5Cpi%7D%7B2%7D)%5D%20%5C%5C%0A%5Cvec%7Be_n%7D%3D(a%5Ccdot%20cos%5C%20%5Ctheta%2C%5Cquad%20a%20%5Ccdot%20sin%5C%20%5Ctheta)%20%5C%5C%0A%5Cend%7Bcases%7D

%5Cbegin%7Bcases%7D%0A%5Cvec%7Bv_t%7D%3D%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D%20%5Ccdot%20%5Cvec%7Be_t%7D%3Dw%5Ccdot%20%5Cvec%7Be_t%7D%20%5C%5C%0A%5Cvec%7Bv_n%7D%3D%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D%20%5Ccdot%20%5Cvec%7Be_n%7D%20%3Dw%5Ccdot%20%5Cvec%7Be_n%7D%5C%5C%0A%0A%5Cend%7Bcases%7D

则标量方程为:


2、使用    得到  

v_t%3Dwr%3Dw%5Ccdot%20a%5Ctheta%3Dwa%5Ccdot%20wt%3Dw%5E2at

这两种方式得到的  一样,可见,对于任意的曲线运动,它的切向速度总是适用  ,但是它的合速度不能盲目地使用这个公式。


题外话——蚊香线

两条互差 180° 得阿基米德螺旋线就是蚊香得形状

代码块
JavaScript
自动换行
复制代码
clc;
clear;
close all;

syms theta
a=1;
x=a*theta*cos(theta+pi);
y=a*theta*sin(theta+pi);
fplot(x,y,[0,pi*2*5],'LineWidth',1.5,'Color','b');
grid on
axis square
hold on
x=a*theta*cos(theta);
y=a*theta*sin(theta);
fplot(x,y,[0,pi*2*5],'LineWidth',1.5,'Color','r');
复制成功

在一个平面上切割两条阿基米德螺旋线,然后横着切两刀,剩下的部分丢掉,把要的部分扒开就是两盘蚊香了。