喵宅苑 MewoGarden × 技术宅社区II | Z站 Z Station 棒棒哒纯文字二次元技术社区

正文

作者:浑狱弥

[i=s] 本帖最后由浑狱弥于 2012-11-17 12:24 编辑今天上课讲到欧拉—班索的特殊情况，而举出的实例关于纬度变迁引出了同学们的不少疑问。什么是纬度变迁？为什么[tex]\vec\omega[/tex]和[tex]\vec k[/tex]夹角不变，但纬度会发生变化？对于这些疑问，老师并没有做明确解答，下来之后鄙人查了一下有关纬度变迁的资料并结合欧拉—班索情况求解。首先是有关于纬度变迁（latitude variation）的具体定义。在以恒星为标准测定地面纬度，地球上同一点的纬度会随时间变化。而引起这种变化的原因大致可分为两类：钱德勒颤动（Chandler wobble )和蒙气差波动(Irregularities in atmospheric refraction)。目前关于此现象的解释已经非常充分。可参照一篇大众科普文（纬度的变迁）。由此可见，此现象是很复杂且伴随不确定性的周期现象。回到本问题的初始，或者是参照关于Chandler wobble--wiki，还是不难发现这是与欧拉—班索情况相关的。起初欧拉就从公式中推证出其周期为305天，不过论文鄙人不可考，于是自己来重新推证。本现象中假设地球为赤道规则呈圆、两极略短的重刚体，且只受重力，过中心转动。即地球做无外力矩的惯性转动，且[tex]I_x=I_y \ne I_z[/tex]。取惯量主轴为动系坐标系，由欧拉动力学方程得： [tex]\begin{cases} I_x\dot{\omega_x}-\left(I_x-I_y\right)\omega_y\omega_z=0 \\ I_y\dot{\omega_y}-\left(I_z-I_x\right)\omega_z\omega_x=0 \\ I_z\dot{\omega_z}=0 \\ \end{cases}[/tex] 由此可知，[tex]\omega_z=\Omega\left(const\right)[/tex]，即为常数。令：[tex]n=\fac{I_z-I_x}{I_x}\Omega[/tex] 得 [tex]\begin{cases} \dot{\omega_x}=-n\omega_y \\ \dot{\omega_y}=n\omega_x \\ \end{cases}[/tex] 解得此方程组为： [tex]\begin{cases} \omega_x=\omega_0\cos\left(nt+\delta\right) \\ \omega_y=\omega_0\sin\left(nt+\delta\right) \\ \end{cases}[/tex] 因此可知：[tex]\vec\omega^2=\left(\omega_x^2+\omega_y^2+\omega_z^2\right)^{\frac{1}{2}}=\left(\omega_0^2+\Omega^2\right)^{\frac{1}{2}}=const[/tex] 即角速度的大小不变，方向随时间变化。取角动量[tex]\vec J[/tex]为静系[tex]\zeta[/tex]轴方向，且[tex]\vec\dot{\phi}\uparrow\uparrow\vec J[/tex] 由欧拉运动学方程： [tex]\vec\dot{\phi}=\dot{\phi}\sin \theta\sin \psi\vec i+\dot{\phi}\sin \theta\cos \psi\vec j+\dot{\phi}\cos \theta\vec k[/tex] [tex]\vec J=J\sin \theta\sin \psi\vec i+J\sin \theta\cos \psi\vec j+J\cos \theta\vec k=J_x\vec i+J_y\vec j+J_z\vec k[/tex] 则，[tex]J\cos \theta=I_z\omega_z[/tex] 由此可依次代入求解，得到：章动：[tex]\theta=\cos^{-1}\left(\frac{I_z\Omega}{J}\right)=\theta_0\left(const\right),\dot{\theta}=0[/tex] 自转：[tex]\psi=\frac{\pi}{2}-\left(nt+\delta \right),\dot{\psi}=-n[/tex] 进动：[tex]\phi=\left(\Omega+nt\right)\sec \theta_0+\phi_0,\dot{\phi}=\frac{\Omega-\dot{\psi}}{\cos \theta}=\sec \theta_0\left(\Omega+nt\right)[/tex] 而其中 [tex]\begin{cases} \omega_x=\omega_0\cos\left(nt+\delta\right) \\ \omega_y=\omega_0\sin\left(nt+\delta\right) \\ \omega_z=\Omega\left(const\right) \end{cases}[/tex] 其中，[tex]\omega_z=\vec\omega \cdot \vec k=\Omega,n=\frac{I_z-I_x}{I_x}\Omega=\frac{1}{300}\Omega,\Omega=\omega\cos{\left(\vec \omega \cdot \vec k\right)\approx\omega[/tex] 于是[tex]T=\frac{2\pi}{n}=305[/tex]天而，将本设定改为[tex]I_x=I_y=I_z[/tex]，则章动、自转均为0，不会形成纬度偏移。欧拉—班索情况的3D图像（顺带吐槽学院所谓的全国唯一）

就是对于

作者:Ekino

浑狱弥发表于 2012-11-17 13:44 就是对于[tex]I_x[/tex]来说，是[tex]\left(I_y-I_z\right)\omega_y\omega_z[/tex]

…………原来是把列维西维塔算符都省了且用了爱因斯坦约定么？太高级了点没反应过来啊…………

查看回复

章动角是啥

作者:浑狱弥

Ekino 发表于 2012-11-17 13:38 章动角是啥？L和w的夹角么？话说御宅网的回复还支持Latex功能貌似很先进呢~~ ...

章动角就是静系轴（进动的主轴）与惯量椭球自转轴的夹角~ 不知道您过去用的啥教材，反正鄙人用的学校自编教材（1楼的链接就是学校理论物理的网站）上面全是点点，都是照着老师的PPT学的~ 见到LaTex大好，站长好像是做数模出生，这也难怪~

查看回复

和我的书上写的不太一样呢

作者:浑狱弥

Ekino 发表于 2012-11-17 13:36 咦，和我的书上写的不太一样呢，话说[tex]I_j_k[/tex]是啥？

就是对于[tex]I_x[/tex]来说，是[tex]\left(I_y-I_z\right)\omega_y\omega_z[/tex]

查看回复

本帖最后由

作者:Ekino

[i=s] 本帖最后由 Ekino 于 2012-11-17 13:40 编辑

浑狱弥发表于 2012-11-17 12:23 有章动角（大小取决于初始条件），但无章动

章动角是啥？L和w的夹角么？话说御宅网的回复还支持Latex功能貌似很先进呢~~

查看回复

本帖最后由

作者:Ekino

[i=s] 本帖最后由 Ekino 于 2012-11-17 13:43 编辑

浑狱弥发表于 2012-11-17 12:17 此时的[tex]I_x=I_y[/tex]所以第一项的[tex]I_x-I_z=I_y-I_z[/tex] 欧拉运动学方程在第二次化简（取惯量主 ...

咦，和我的书上写的不太一样呢，话说[tex]I_j_k[/tex]是啥？（弱弱的提醒一下原文写的是[tex]I_x-I_y[/tex]是个笔误~）

查看回复

没有外力矩的情况下会有章动么

作者:浑狱弥

Ekino 发表于 2012-11-17 10:03 没有外力矩的情况下会有章动么？

有章动角（大小取决于初始条件），但无章动

查看回复

本帖最后由

作者:浑狱弥

[i=s] 本帖最后由浑狱弥于 2012-11-17 12:21 编辑

Ekino 发表于 2012-11-17 09:50 他的第一个欧拉运动方程貌似不对啦，欧拉方程是 I_i*d(w_i)/dt+e_ijk*w_j*L_k=N_i，所以第一项转动惯量项应 ...

此时的[tex]I_x=I_y[/tex]所以第一项的[tex]I_x-I_z=I_y-I_z[/tex] 欧拉运动学方程在第二次化简（取惯量主轴为坐标轴）之后就是[tex]I_i{\operatorname{d}\!\omega\over\operatorname{d}\!x}-I_{jk} {\omega}_{kj}=M[/tex]

查看回复

没有外力矩的情况下会有章动么

作者:Ekino

没有外力矩的情况下会有章动么？

查看回复