查看“理论力学”的源代码

<!-- [[文件:Wiki.gif|right]] -->
参考书籍：
#[http://www.amazon.cn/gp/product/B0011C0H0O 朗道，栗弗席兹，《力学（第3版）》（英文影印版），世界图书出版公司]
#;朗道十卷本理论物理教程之首。从分析力学讲起，高屋建瓴，推理严密，内容精炼（170页）。可供仔细研读，若感觉难度太高可结合鞠国兴《朗道力学解读》。
#[http://www.amazon.cn/经典力学/dp/B00475QFVC Goldstein et al.，《经典力学（第3版）》（英文影印版），高等教育出版社]
#;美国普遍采用的理论力学教材，内容丰富，数学推导详尽，有大量习题。内容略多（约630页）。可供选择性阅读，与教材互补。中译本为第2版，内容更加简练，比第3版缺少“混沌”一章。
#[http://www.amazon.cn/经典力学的数学方法-V-I-Arnold/dp/B001174V74 Arnold，《经典力学的数学方法》（英文影印版），世界图书出版公司]
#;经典名著，内容有深度，数学要求略高。可供查阅。
#[http://www.amazon.cn/分析动力学-陈滨/dp/B009SZWSL2 陈滨，《分析动力学（第2版）》，北京大学出版社]
#;力学专业教材，对分析力学进行了细致、深入、大量的讨论。可供参考。

== 质点力学 ==
=== 运动学 ===
----
#运动学变量的矢量表示：坐标 <math> \begin{align}\boldsymbol{r}\end{align} </math>，速度 <math>\boldsymbol{v}=\dot{\boldsymbol{r}}</math>，加速度 <math>\boldsymbol{a}=\dot{\boldsymbol{v}}=\ddot{\boldsymbol{r}}</math>，动量 <math>\boldsymbol{p}=m\boldsymbol{v}</math>，角动量 <math>\boldsymbol{J}=\boldsymbol{r}\times\boldsymbol{p}</math>
#;三种坐标系中运动学变量的具体形式（注意：坐标轴方向是否随时间变化）
##笛卡尔坐标系 <math>\begin{align}\left\{\begin{array}{l}
\boldsymbol{r}=x\boldsymbol{i}+y\boldsymbol{j}+z\boldsymbol{k},\\
\boldsymbol{v}=\dot x\boldsymbol{i}+\dot y\boldsymbol{j}+\dot z\boldsymbol{k},\\
\boldsymbol{a}=\ddot x\boldsymbol{i}+\ddot y\boldsymbol{j}+\ddot z\boldsymbol{k},
\end{array}\right.\end{align}</math>
##平面极坐标系 <math>\begin{align}\left\{\begin{array}{l}
\boldsymbol{r}=r\boldsymbol{i},\\
\boldsymbol{v}=\dot r\boldsymbol{i}+r\dot \theta\boldsymbol{j},\\
\boldsymbol{a}=(\ddot r-r\dot\theta^2)\boldsymbol{i}+\frac{1}{r}\frac{\rm d}{{\rm d}t}(r^2\dot\theta)\boldsymbol{j},
\end{array}\right. 
\hspace{7ex}
\left\{\begin{array}{l}
\frac{{\rm d}\boldsymbol{i}}{{\rm d}t}=\dot\theta \boldsymbol{j},\\
\frac{{\rm d}\boldsymbol{j}}{{\rm d}t}=-\dot\theta \boldsymbol{i},\end{array}
\right.
\end{align}</math>
##自然坐标系 <math>\begin{align}\left\{\begin{array}{l}
\boldsymbol{r}=\boldsymbol{r},\\
\boldsymbol{v}=\dot s\boldsymbol{i},\\
\boldsymbol{a}=\dot v\boldsymbol{i}+\frac{v^2}{\rho}\boldsymbol{j}=a_{\rm t}\boldsymbol{i}+a_{\rm n}\boldsymbol{j},
\end{array}\right.
\hspace{7ex}
\left\{\begin{array}{l}
\frac{{\rm d}\boldsymbol{i}}{{\rm d}t}=\frac{v}{\rho} \boldsymbol{j},\\
\frac{{\rm d}\boldsymbol{j}}{{\rm d}t}=-\frac{v}{\rho} \boldsymbol{i},\end{array}
\right.
\end{align}
\hspace{7ex}
\rho=\left|\frac{{\rm d}s}{{\rm d}\theta}\right|=\left|\frac{y^{\prime\prime}}{\left(1+y^\prime\right)^{3/2}}\right|
</math>
##;三维欧式空间微分几何（参考苏步青《微分几何学》）
##;<math>\boldsymbol{i}=\boldsymbol{e}_{\rm t}</math>：切向单位矢量，t=tangent（切向）；
##;<math>\boldsymbol{j}=\boldsymbol{e}_{\rm n}</math>：法向单位矢量，n=normal（法向）；
##;<math>\boldsymbol{k}=\boldsymbol{e}_{\rm b}</math>：从法向单位矢量，b=binormal（从法向/副法向）。
##;[[文件:教学-理论力学-ch-1-密切平面.png|320px]]
##扩展：球坐标、抛物线坐标、椭圆坐标（见周衍柏第 <math>\S5.10</math> 节）

===质点运动定律===
----
<ol>
<li>牛顿三定律
#质点不受其他物体作用 <math>\Longrightarrow {\it p}=</math> 常数
#<math>\boldsymbol{F}=m\boldsymbol{a}</math>
#<math> \boldsymbol{F}_{12}=-\boldsymbol{F}_{21},\boldsymbol{r}_{12}\times\boldsymbol{F}_{12}=0</math>
</li>
<li>质点运动微分方程：
<ol > <li>自由粒子运动微分方程：  <math>m\ddot{\boldsymbol{r}}=\boldsymbol{F}(\boldsymbol{r},\dot{\boldsymbol{r}},t)</math>。
分量形式为：<br/>
{| border="1" cellspacing="0" cellpadding="5" align="left"
| 笛卡尔系 || 平面极坐标系 || 自然坐标系 
|-
|<math>\begin{align}\left\{\begin{array}{l}
m\ddot{x}=F_x(x,y,z;\dot x,\dot y,\dot z;t)\\
m\ddot{y}=F_y(x,y,z;\dot x,\dot y,\dot z;t)\\
m\ddot{z}=F_z(x,y,z;\dot x,\dot y,\dot z;t)
\end{array}\right.\end{align}</math> 
||
<math>\begin{align}\left\{\begin{array}{l}
m(\ddot r-r{\dot\theta}^2)=F_r(r,\theta;\dot r,\dot\theta;t),\\
m(r\ddot\theta +2\dot r \dot\theta)=F_\theta(r,\theta;\dot r,\dot\theta;t),\\
\end{array}\right.\end{align}</math>
||
<math>\begin{align}\left\{\begin{array}{l}
m\dot v=F_{\rm t},\\
m\frac{v^2}{\rho}=F_{\rm n},\\
0=F_{\rm b}.
\end{array}\right.\end{align}</math>
|}
<br/><br/><br/><br/><br/><br/><br/>
</li>
<li>非自由粒子运动微分方程：<math>\begin{align}\left\{\begin{array}{l}m\ddot{\boldsymbol{r}}=\boldsymbol{F}(\boldsymbol{r},\dot{\boldsymbol{r}},t)+\boldsymbol{R},\\
\varphi(\boldsymbol{r})=0~~~(约束方程)\end{array}\right.\end{align}</math><br/>
分量形式（将<math>\boldsymbol{F}\to \boldsymbol{F}+\boldsymbol{R}</math>）
自然坐标系：<math>\begin{align}\left\{\begin{array}{l}
m\dot v=F_{\rm t},\\
m\frac{v^2}{\rho}=F_{\rm n}+R_{\rm n},\\
0=F_{\rm b}+R_{\rm b},\\
\varphi(\boldsymbol{r})=0.
\end{array}\right.\end{align}</math> 
</li></ol>
</li>
<li>求解质点运动 <math>\to</math> 求解二阶线性/非线性常微分方程
两种解法：（1）数值解法（差分法，适用于线性/非线性、一阶/二阶/高阶、常微分/偏微分等各种方程）；（2）解析法
</li></ol>

===质点运动定理===
----
#保守力：如果 <math>\boldsymbol{F}(\boldsymbol{r})</math> 做功与路径无关，只依赖于出发点和终点坐标，则 <math>\boldsymbol{F}(\boldsymbol{r})</math> 叫保守力。
#;数学表述：曲线积分 <math>\oint_A^B \boldsymbol{F}(\boldsymbol{r})\cdot{\rm d}\boldsymbol{r}</math> 与路径无关，而仅依赖于端点 A 和 B .
#;根据高数知识，必然存在标量函数 <math>V(\boldsymbol{r})</math> 满足 <math>\boldsymbol{F}(\boldsymbol{r})=-\nabla V(r)</math>
#;物理上，标量函数 <math>V(\boldsymbol{r})</math> 为保守力 <math>\boldsymbol{F}(\boldsymbol{r})</math> 对应的势能。
#;保守力：conservative force，原始意义是保持系统能量守恒（conserved）的作用力。(见朗道第 <math>\S6</math> 节)
#;
#质点运动定理和守恒定律
##动量定理：<math>\dot{\boldsymbol{p}}=\boldsymbol{F}</math>或<math>{\rm d}\boldsymbol{p}={\rm d}\boldsymbol{F}{\rm d}t</math>；动量守恒定律：<math>\boldsymbol{F}=0\Longrightarrow\boldsymbol{p}=</math>常数
##角动量定理：<math>\dot{\boldsymbol{J}}=\boldsymbol{M}</math>或<math>{\rm d}\boldsymbol{J}=\boldsymbol{M}{\rm d}t</math>；角动量守恒定律：<math>\boldsymbol{M}=0\Longrightarrow \boldsymbol{J}=</math>常数
##动能定理：<math>\nabla T=\boldsymbol{F}</math>或<math>{\rm d}T=\boldsymbol{F}\cdot{\rm d}\boldsymbol{r}</math>       （动能 <math>T\equiv \frac{1}{2}mv^2</math>）；机械能守恒定律：<math> \boldsymbol{F} </math> 为保守力 <math>\Longrightarrow E=</math>常数       （总能量 <math>E=T+V</math>）

===有心力===

----

#质点在有心力场中运动 <math>\Longrightarrow</math> 角动量守恒（<math>\dot{\boldsymbol{J}}=0</math>）<math>\Longrightarrow</math> 平面运动
#有心力必定为保守力，<math>\boldsymbol{F}=-\nabla U(r)</math>，其中 <math>U(r)</math> 为势能
#质点在有心力场中的运动等价于一维运动：<math>m\ddot r=-\nabla V(r)</math><math>~~~~~\left(V(r)=U(r)+\frac{\boldsymbol{J}^2}{2mr^2}\right)</math>
#轨道有限的条件：有效势能<math>V(r)</math>存在势阱，能量<math>E</math>小于势阱两边势垒的高度，从而质点能够束缚在势阱中运动。         <br />假设势能取幂函数形式，<math>U(r)=cr^\alpha</math>，根据上述条件有<math>V^\prime(r)=0,V^{\prime\prime}(r)>0</math>，从而得到轨道有限的必要条件为<math>\begin{align}\left\{\begin{array}{ll}
c>0,&\alpha>0,\\
c<0,&-2<\alpha<0
\end{array}\right.\end{align}</math>
#质点在有心力场中运动的轨道封闭的条件：<math>\Delta\theta=2\int_{r_{\rm min}}^{r_{\rm max}}\frac{J/r^2{\rm d}r}{\sqrt{2m(E-V(r))}}=\pi \times</math>有理数
#圆形轨道的稳定性：从等效一维运动的角度讨论，稳定性意味着粒子束缚在一个势阱中运动。因此稳定性条件与前边的有限性条件相同。
#万有引力（平方反比引力）作用下运动方程的解：<math>r=\frac{p}{1+e\cos\theta}.\hspace{3ex}\left(p=\frac{J^2}{mk^2}，e=\sqrt{1+\frac{2EJ^2}{mk^4}}\right)</math>
#平方反比斥力：<math>\alpha</math> 粒子的散射，运动方程的解为 <math> r=\frac{p}{e\cos\theta-1}.\hspace{3ex}\left(p=\frac{J^2}{mk^2}，e=\sqrt{1+\frac{2EJ^2}{mk^4}}\right)</math>.
#扩展（SRT和毕业论文候选）
##双恒星系统中行星的轨道？居住在星星上的人的历法会如何发展？

=== 常用数学公式 ===
----
<ol >
<li><math>  \begin{align}
    \int\frac{{\rm d}x}{1+x^2}&=\arctan x+C,& \int\frac{{\rm d}x}{\sqrt{1+x^2}}&={\rm arcsinh}\, x+C,\\
    \int\frac{{\rm d}x}{1-x^2}&={\rm arctanh}\,x+C,& \int\frac{{\rm d}x}{\sqrt{1-x^2}}&=\arcsin x+C.
  \end{align}
</math></li>

<li>二阶常系数齐次线性微分方程 <math> \ddot x+p\dot x+q x=0</math>，解的形式取决于特征方程 <math> r^2+pr+q=0</math> 解的情况。<br/>
<math>\begin{align}x(t)=\left\{\begin{array}{ll}
C_1 e^{r_1 t}+C_2 e^{r_2 t},& (2~个不同实根~ r_1\neq r_2)\\
\left(C_1+C_2x\right)e^{rt},& (2重实根~ r_1=r_2\equiv r)\\
e^{\alpha t}\left(C_1\cos\beta t+C_2\sin\beta t\right),&(2~个共轭复根~  r_{1,2}=\alpha\pm i\beta)
\end{array}\right.\end{align}</math></li>

<li>二阶常系数非齐次线性微分方程 <math> \ddot x+p\dot x+q x=f(x)</math>，通解为 <math>y=y^*+Y</math>，其中 <math>Y</math> 是对应的齐次方程即当 <math>f(x)=0</math> 的时候方程的通解，而 <math> y^* </math> 是非齐次方程的特解。当 <math>f(x)=e^{cx}P_l(x) </math> 时（<math>P_n(x) </math> 为 <math>x</math> 的<math>n</math> 次多项式，<math>c</math> 为复数），特解 <math>y^* </math> 有简单的形式。<br/>
取 <math>f(x)=e^{\lambda x}\left(P_n(x)\cos\omega t+P_l(x)\sin\omega t\right)</math>，则 
<math>\begin{align}y^*=\left\{\begin{array}{l@{~~~~~}l}
   e^{\lambda x}\left(Q_m(x)\cos\omega t+R_m(x)\sin\omega t\right),&\left(\lambda+i\omega\text{不是特征方程}r^2+pr+q=0\text{的根}\right)\\
   xe^{\lambda x}\left(Q_m(x)\cos\omega t+R_m(x)\sin\omega t\right),&\left(\lambda+i\omega\text{是特征方程}r^2+pr+q=0\text{的根}\right)
\end{array}\right.\end{align}</math>
</li>

<li>一阶微分方程 <br/>
<math>\begin{align}\begin{array}{|l|@{~~~}l|@{~~~}l|}\hline
\text{方程类型} & \text{方程形式} & \text{方程通解}\\\hline
\text{可分离变量} & y^\prime=g(x)h(y) & \int\frac{{\rm d}y}{h(y)}=\int g(x){\rm d}x+C\\\hline
\text{齐次微分方程} & y^\prime=\varphi\left(\frac{y}{x}\right) & \int\frac{{\rm d}\tilde y}{\varphi(\tilde y)-\tilde y}=\ln x+C,~~~\left(\tilde y\equiv\frac{y}{x}\right)\\\hline
\text{线性微分方程} & y^\prime+P(x)+Q(y)=0 & y=e^{-\int P(x){\rm d}x}\left[-\int Q(x)e^{\int P(x){\rm d}x}{\rm d}x+C\right]\\\hline
\text{伯努利方程} & y^\prime+P(x)y+Q(y) y^n=0 & \frac{1}{1-n}\tilde y^\prime+P(x)\tilde y^\prime+Q(x)=0,~~~(\tilde y=y^{1-n})\\\hline
\end{array}\end{align}</math>
</li>

</ol>

== 质点组力学 ==
===质点组动力学===
----
#运动学：
#*质点组总动量：<math>\boldsymbol{P}=m\boldsymbol{v}_c+\boldsymbol{0}</math>
#*质点组角动量：<math>\begin{align} \boldsymbol{J}=\boldsymbol{r}\times m\boldsymbol{v}_c+\boldsymbol{J}^\prime\end{align}</math>
#*质点组动能（柯尼希定理）：<math>\begin{align} T=\frac{1}{2}m\boldsymbol{v}_c^2+T^\prime\end{align}</math>，即<math>\begin{align}\frac{1}{2}\sum_{i=1}^n m_i \dot{\boldsymbol{r}}_i^{2}=\frac{1}{2}m\dot{\boldsymbol{r}}_C^2+\frac{1}{2}\sum_{i=1}^n m_i \dot{\boldsymbol{r}}_i^{\prime 2}\end{align}</math>，其中 <math>\begin{align} m=\sum_{i=1}^n m_i \end{align}</math>
#动力学
#*质点组动量定理：<math>\begin{align} {\rm d}\boldsymbol{p}=\left(\sum_{i=1}^n \boldsymbol{F}_i^{(\rm e)}\right){\rm d}t \end{align}</math>
#*质点组角动量定理：<math>\begin{align} {\rm d}\boldsymbol{J}=\left(\sum_{i=1}^n \boldsymbol{M}_i^{(\rm e)}\right){\rm d}t \end{align}</math>
#*质点组动能定理：<math>\begin{align} {\rm d}T=\sum_{i=1}^n \boldsymbol{F}_i^{(\rm e)}\cdot {\rm d}\boldsymbol{r}_i+\sum_{i=1}^n \boldsymbol{F}_i^{(\rm i)}\cdot {\rm d}\boldsymbol{r}_i \end{align}</math>

=== 两体问题 ===
----
两质点组成的系统，如果合外力为 0 即 <math>0=\boldsymbol{F}_1^{(\rm e)}+\boldsymbol{F}_2^{(\rm e)}</math>，对每个质点分别列运动方程：<math>\begin{align}\left\{\begin{array}{l}
m_1\ddot{\boldsymbol{r}}_1=\boldsymbol{F}_1^{(\rm e)}+\boldsymbol{F}_{21}^{(\rm i)},\\
m_2\ddot{\boldsymbol{r}}_2=\boldsymbol{F}_2^{(\rm e)}+\boldsymbol{F}_{12}^{(\rm i)},
\end{array}\right.\end{align}</math>
，其中 <math>\boldsymbol{F}_{21}^{(\rm i)}=-\boldsymbol{F}_{12}^{(\rm i)}</math>。
#两式之和为质心运动微分方程：<math>\begin{align}m\ddot{\boldsymbol{r}}_c=0,\hspace{2ex}\left(m=m_1+m_2\right)\end{align}</math>
#;
#两式分别乘以 <math>m_2</math> 和 <math>m_1</math>，相减得到：<math>m_1m_2\left(\ddot{\boldsymbol{r}}_1-\ddot{\boldsymbol{r}}_2\right)=(m_1+m_2)\left(\boldsymbol{F}_1^{(\rm e)}+\boldsymbol{F}_{21}^{(\rm i)}\right)</math>，即 <math>\mu\ddot{\boldsymbol{r}}=\boldsymbol{F}_1^{(\rm e)}+\boldsymbol{F}_{21}^{(\rm i)}</math>，其中 <math>\boldsymbol{r}\equiv \boldsymbol{r}_1-\boldsymbol{r}_2</math>为两质点的相对位矢， <math>\mu\equiv \frac{m_1m_2}{m_1+m_2}</math> 叫做折合质量。<br />如果两质点不受外力，且两质点相互作用力仅依赖于距离（是保守力），那么 <math>\boldsymbol{F}_{21}^{(\rm i)}=-\nabla U(\boldsymbol{r})</math>，<math>\big(U(\boldsymbol{r})</math> 为两质点相互作用势能<math>\big)</math>，上式可化为有心力场中的单质点运动问题 <math>\begin{align} \mu \ddot{\boldsymbol{r}}=-\nabla U(\boldsymbol{r})\end{align}</math>。

简言之，不受外力的两质点的运动，可以等价为<span style="color: blue">'''质心的匀速直线运动'''</span>与<span style="color: blue">'''具有折合质量的质点的相对运动'''</span>的合成，运动方程分别为
<math>\begin{align}\left\{\begin{array}{l}
m\ddot{\boldsymbol{r}}_c=0,\\
\mu\ddot{\boldsymbol{r}}=-\nabla U(\boldsymbol{r}).
\end{array}\right.\end{align}</math>

=== 位力定理 ===
----
对于满足'''坐标和动量均取有限值'''的运动（包括周期性运动），质点组动能的时间平均值等于位力积时间平均值的负值。即

<math>\overline{T}=-\frac{1}{2}\overline{\sum_{i=1}^n\boldsymbol{F}_i\cdot\boldsymbol{r}_i}</math>

如果势能是坐标的 <math>n</math>次函数，即 <math>V(\boldsymbol{r})=cr^n</math>，则 <math>\boldsymbol{r}\cdot\nabla V(\boldsymbol{r})=nV(\boldsymbol{r})</math>。因此我们得到

<math>\overline{T}=\frac{n}{2}\overline{V}</math>

在统计力学中，位力定理可以导出理想气体状态方程 <math>PV=nRT</math>，也是导出非理想气体的一种系统的方法；在量子力学中位力定理也有重要应用。更详细的讨论见 [https://en.wikipedia.org/wiki/Virial_theorem 位力定理]

== 刚体力学 ==
=== 刚体运动学 ===

----
<ol>

<li>刚体有6个自由度，可以取为固定在刚体上的一点O的位矢 <math>\boldsymbol{r}_O</math> 和转动的角度即三个欧拉角 <math>\theta,~~\varphi,~~~\psi</math>。刚体一般运动为参考点O的平动与绕参考点转动的合成。</li>
<li>取固定在刚体上的三个相互垂直的轴为坐标轴，建立随刚体运动的坐标系<math>o-xyz</math>。（注意：参考系一般选取地面静止系，与坐标系要区分开。）
<br/>刚体的位置可用欧拉角描述，他们分别为进动角 <math>\varphi</math>、自转角 <math>\psi</math>、章动角 <math>\theta</math> ：
<br/> [[文件:欧拉角.png|152px]]           </li>

<li><span style="color: blue">'''欧拉运动学方程'''</span>：<math>\begin{align}\left\{\begin{array}{l}
\omega_x=\dot\varphi\sin\theta\sin\psi+\dot\theta\cos\psi,\\
\omega_y=\dot\varphi\sin\theta\cos\psi-\dot\theta\sin\psi,\\
\omega_z=\dot\varphi\cos\theta+\dot\psi,
\end{array}\right.\end{align}</math> 
<br>上式来源于 <math>\boldsymbol{\omega}=\dot\varphi\boldsymbol{e}_\zeta+\dot\theta\boldsymbol{e}_{ON}+\dot\psi \boldsymbol{e}_z</math>，其中 <math>\zeta</math> 轴和 <math>ON</math> 轴的方位角和极角分别为 <math>\left(\theta,\frac{\pi}{2}-\psi\right)</math> 和 <math>\left(\frac{\pi}{2},-\psi\right)</math> 。刚体转动方向的三个单位矢量在刚体随动系中的表达式为<math>\begin{align}\left\{\begin{array}{l}
\boldsymbol{e}_\zeta=\sin\theta\sin\psi\boldsymbol{i}+\sin\theta\cos\psi\boldsymbol{j}+\cos\theta\boldsymbol{k},\\
\boldsymbol{e}_{ON}=\cos\psi\boldsymbol{i}-\sin\psi\boldsymbol{j},\\
\boldsymbol{e}_z=\boldsymbol{k},
\end{array}\right.\end{align}</math>  <br/>刚体角速度为矢量，按照普通矢量相加 <math>\boldsymbol{\omega}=\boldsymbol{\omega}_1+\boldsymbol{\omega}_2</math></li>

<li> 刚体转动惯量张量为 <math>\begin{align}\hat A=\left(\begin{array}{ccc}
I_{xx} & -I_{xy} & -I_{xz}\\
-I_{xy} & -I_{yy} & -I_{yz}\\
-I_{xz} & -I_{yz} & I_{zz}
\end{array}\right)\end{align}</math>

其中惯量张量矩阵元的一般定义式为
<math>\begin{align}I_{ij}=\int \left(\boldsymbol{r}^2\delta_{ij}-r_ir_j\right)\rho(\boldsymbol{r}){\rm d}^3\boldsymbol{r}\end{align}</math>，即
<math>\begin{align}I_{xx}=\int\left(y^2+z^2\right)\rho(\boldsymbol{r}){\rm d}^3\boldsymbol{r},\\
I_{yy}=\int\left(x^2+z^2\right)\rho(\boldsymbol{r}){\rm d}^3\boldsymbol{r},\\
I_{zz}=\int\left(x^2+y^2\right)\rho(\boldsymbol{r}){\rm d}^3\boldsymbol{r},\end{align}\hspace{5ex}</math>
<math>\begin{align}I_{xy}=\int xy\rho(\boldsymbol{r}){\rm d}^3\boldsymbol{r},\\
I_{yz}=\int yz\rho(\boldsymbol{r}){\rm d}^3\boldsymbol{r},\\
I_{xz}=\int xz\rho(\boldsymbol{r}){\rm d}^3\boldsymbol{r}.\end{align}</math>
</li>

<li>刚体角动量为 <math>\boldsymbol{J}=\hat{\boldsymbol{A}}\cdot\boldsymbol{\omega}</math>，在笛卡尔系中写为 <math>\begin{align}
\left(\begin{array}{l}J_x\\J_y\\J_z\end{array}\right)
=\left(\begin{array}{rrr}I_{xx}&-I_{xy}&-I_{xz}\\-I_{yx}&I_{yy}&-I_{yz}\\-I_{zx}&-I_{zy}&I_{zz}\end{array}\right)
\left(\begin{array}{l}\omega_x\\\omega_y\\\omega_z\end{array}\right)\end{align}</math> <br/>刚体转动动能为 <math>T=\frac{1}{2}\boldsymbol{\omega}\cdot \hat{\boldsymbol{A}}\cdot\boldsymbol{\omega}=\frac{1}{2}I_n\omega^2</math>，其中 <math>\vec\omega=\omega\hat e_n</math>，<math>I_n\equiv \hat e_n\cdot\hat{\boldsymbol{A}}\cdot\hat e_n</math> 为刚体相对于过质心且沿 <math>\hat e_n</math> 方向的轴的转动惯量。<br/>空间转动变换，可以将转动惯量张量 <math>\hat{\boldsymbol{A}}</math> 对角化，此时的坐标轴为惯量主轴。<br/>惯量主轴的求解：惯量椭球 <math>\boldsymbol{r}\cdot\hat{\boldsymbol{A}}\cdot\boldsymbol{r}=1</math> 的三个对称轴。</li>

<li>刚体上任意一点的速度：<math>\boldsymbol{ v}_A=\boldsymbol{ v}_O+\boldsymbol{\omega}\times\boldsymbol{ r}_{OA}</math></li>

<li>刚体的动能有两部分组成：质心动能和转动动能，即<math>T=\frac{1}{2}m\boldsymbol{v}_c^2+\frac{1}{2}I_n\omega^2</math></li>

</ol>

=== 刚体动力学 ===

----
刚体动力学方程包括刚体动量定理、刚体角动量定理（即欧拉动力学方程）和刚体动能定理。

<ol>
<li> 刚体动量定理：<math>m\ddot{\boldsymbol{r}}_c=\boldsymbol{F}</math>，刚体质心运动由刚体受到的合外力决定。</li>

<li> 刚体角动量定理：<math>\dot{\boldsymbol{J}}=\boldsymbol{M}</math>，其中角动量 <math>\boldsymbol{J}=\hat{\boldsymbol{A}}\cdot\boldsymbol{\omega}</math>。为简化计算，欧拉采用了两重简化方式
* 采用固定在刚体上的坐标系，则转动惯量张量 <math>{\rm d}\hat{\boldsymbol{A}}/{\rm d}t=0</math>，角速度 <math>\boldsymbol{\omega}=\omega_x\boldsymbol{i}+\omega_y\boldsymbol{j}+\omega_z\boldsymbol{k}</math>，其中坐标轴的单位矢量随时间变化，满足
<math>\begin{align}\frac{{\rm d}\boldsymbol{j}}{{\rm d}t}=\boldsymbol{\omega}\times\boldsymbol{i},~~~
\frac{{\rm d}\boldsymbol{j}}{{\rm d}t}=\boldsymbol{\omega}\times\boldsymbol{j},~~~
\frac{{\rm d}\boldsymbol{k}}{{\rm d}t}=\boldsymbol{\omega}\times\boldsymbol{k}\end{align},</math>
用<math>\boldsymbol{i,j,k}</math>展开的任意矢量<math>\boldsymbol{G}=G_x\boldsymbol{i}+G_y\boldsymbol{j}+G_z\boldsymbol{k}</math>，对时间的导数为<math>\begin{align}\frac{{\rm d}\boldsymbol{G}}{{\rm d}t}=\frac{{\rm d}^*\boldsymbol{G}}{{\rm d}t}+\boldsymbol{\omega}\times\boldsymbol{G}\end{align}</math>。因此刚体角动量定理变为：<br />
<center><math>\begin{align}\frac{{\rm d}^*\boldsymbol{J}}{{\rm d}t}+\boldsymbol{\omega}\times\boldsymbol{J}=\boldsymbol{M}.\end{align}</math></center>
* 取<math>o-xyz</math>的坐标轴为刚体的惯量主轴，则惯量张量为对角矩阵即<math>\begin{align}\hat{\boldsymbol{A}}=\left(\begin{array}{ccc}
I_1 & 0 & 0\\
0 & I_2 & 0\\
0 & 0 & I_3
\end{array}\right)\end{align}</math>
因此角动量形式为 <math>\begin{align}\boldsymbol{J}=\hat{\boldsymbol{A}}\cdot\boldsymbol{\omega}=I_1\omega_x\boldsymbol{i}+I_2\omega_y\boldsymbol{j}+I_3\omega_z\boldsymbol{k}\end{align}</math>。代入刚体角动量定理，其分量形式即为<span style="color: blue">'''欧拉动力学方程 '''</span>
<center><math>\begin{align}\left\{\begin{array}{l}
I_1\dot\omega_x+(I_3-I_2)\omega_y\omega_z=M_x,\\
I_2\dot\omega_y+(I_1-I_3)\omega_z\omega_x=M_y,\\
I_3\dot\omega_z+(I_2-I_1)\omega_x\omega_y=M_z.
\end{array}\right.\end{align}</math></center></li>

<li> 刚体动能定理：<math>\begin{align}{\rm d}T=\sum_{i=1}^n\boldsymbol{F}^{(\rm e)}_i\cdot{\rm d}\boldsymbol{r}_i\end{align}</math>，其中刚体动能为<math>T=\frac{1}{2}mv_c^2+\frac{1}{2}I_n\omega^2</math></li>
</ol>

=== 刚体运动分类 ===

----
<ol>
<li> 刚体平动：欧拉角为常数，角速度<math>\boldsymbol{\omega}=0</math>。自由度数 <math> s=3</math>，取为刚体质心位矢。 </li>
<li> 刚体定轴转动：刚体上一点固定，且角速度仅沿一个方向。自由度数 <math>s=1</math>，取为自转角。<br>
把角速度方向取为z方向，则<math>\omega_x=\omega_y=0</math>，角动量简化为 <math>\boldsymbol{J}=I_{xz}\omega_z\boldsymbol{i}+I_{yz}\omega_z\boldsymbol{j}+I_{zz}\omega_z\boldsymbol{k}</math>，刚体的角动量定理变为<br />
<center><math>\begin{align}\left\{\begin{array}{r}
I_{xz}\dot\omega_z-I_{yz}\omega_z^2=M_x,\\
I_{yz}\dot\omega_z+I_{xz}\omega_z^2=M_y,\\
I_{zz}\dot\omega_z=M_z,
\end{array}\right.\end{align}</math> </center>  最后一式是动力学方程，前两式可以用来求解对刚体转轴的约束力。
</li>
<li> 刚体平面平行运动：刚体上一点做平面运动，且刚体角速度垂直于这个平面。自由度数 <math> s=3 </math>，取为刚体上一点在平面内的位矢和刚体的转动角度。<br/>
刚体的角动量定理同上，动量定理只涉及x和y方向，即<math>\begin{align}\left\{\begin{array}{l}
m\ddot{\boldsymbol{x}}_c=\boldsymbol{F}_x,\\
m\ddot{\boldsymbol{y}}_c=\boldsymbol{F}_y,
\end{array}\right.\end{align}</math><br/>
平面平行运动的刚体必有转动瞬心<math>O</math>。转动瞬心可以简化刚体的运动描述
# 刚体的总动能可以写为 <math>T=\frac{1}{2}I_O\omega^2</math> ，<math>I_O</math> 是绕瞬心的转动惯量
# 刚体上一点的速度可以写为 <math>\boldsymbol{v}=\boldsymbol{\omega}\times\boldsymbol{r}^\prime</math>，<math>\boldsymbol{r}^\prime</math> 是该点相对于转动瞬心的相对位矢。<br/>
转动瞬心的求法：
# 已知刚体上两点速度，如果速度方向不同，则瞬心为两速度垂线的交点；
# 已知刚体上一点速度和刚体角动量，则瞬心位于过该点且垂直于速度的直线上，且与该点的距离为<math> \frac{v}{\omega}</math>；
# 已知两点的速度，速度方向相同，且过一点速度的垂线穿过第二点，则可以建立相似三角形确定速度；
# 如果两点速度平行且速度的垂涎穿过两点，则刚体的转动瞬心在无穷远处。
<br/>
[[文件:瞬心-1.png|150px]][[文件:瞬心-2.png|135px]]<math>\hspace{3ex}</math>[[文件:瞬心-3.png|165px]]
<br/>
转动瞬心划过的轨迹，在地面静止系中叫做空间极迹；在刚体静止系中叫做本体极迹。
</li>
<li> 刚体定点转动：刚体上有一固定点O，角速度无限制。自由度数 <math>s=3</math>，取为刚体的三个欧拉角。<br/>
刚体定点运动的动力学方程完整的<span style="color: blue">'''欧拉动力学方程'''</span>。<br/>
重刚体：主动力仅有重力的刚体。有解析解的重刚体有三种：
<ol>
<li> 欧拉-潘索情况：外力通过固定点O，<math>\boldsymbol{M}=0</math>，刚体在惯性下运动。欧拉动力学方程简化为
<center><math>\begin{align}\left\{\begin{array}{l}
I_1\dot\omega_x-(I_2-I_3)\omega_y\omega_z=0,\\
I_2\dot\omega_y-(I_3-I_1)\omega_z\omega_x=0,\\
I_3\dot\omega_z-(I_1-I_2)\omega_x\omega_y=0.
\end{array}\right.\end{align}</math></center>
根据上述方程，可以求出<math>\omega_x,\omega_y,\omega_z</math>；再利用欧拉动力学方程，求出欧拉角随时间变化的方程<math>\theta(t),\varphi(t),\psi(t)</math>。这样运动方程就完全解出。<br/>
当<math>I_1=I_2</math>的特殊情况下（比如地球），运动方程可以得到简化。<br/>
潘索引入了几何方法，可以比较简便的说明一般情况下（<math>I_1\neq I_2\neq I_3</math>）刚体的运动图像。</li>
<li> 拉格朗日-泊松情况（对称陀螺）：重力作用点不过固定点O，但有<math>I_1=I_2</math>。欧拉动力学方程化简为
<center><math>\begin{align}\left\{\begin{array}{l}
I_1\dot\omega_x-(I_1-I_3)\omega_y\omega_z=mgl\sin\theta\cos\psi,\\
I_1\dot\omega_y+(I_3-I_1)\omega_z\omega_x=mgl\sin\theta\sin\psi,\\
I_3\dot\omega_z=0.
\end{array}\right.\end{align}</math></center>
根据欧拉运动学方程，易得对称陀螺的拉格朗日函数 <math>\begin{align}L=\frac{1}{2}I_1\left(\dot\theta^2+\sin^2\theta\dot\varphi^2\right)+\frac{1}{2}I_3\left(\cos\theta\dot\varphi+\dot\psi\right)^2-mgl\cos\theta\end{align}</math>。因为拉格朗日函数对<math>\varphi</math>和<math>\psi</math>没有依赖，因此对应的广义动量为一次积分（守恒量）；同时，因为动量是广义速度的二次函数，故能量<math>E=T+V</math>也为守恒量。这样我们就得到对称陀螺的三个一次积分（守恒量）。
<math>\begin{align}\left\{\begin{array}{r}
I_1\sin^2\theta\dot\varphi+I_3\cos\theta\left(\cos\theta\dot\varphi+\dot\psi\right)=c_1,\\
I_3\left(\cos\theta\dot\varphi+\dot\psi\right)=c_2,\\
\frac{1}{2}I_1\left(\dot\theta^2+\sin^2\theta\dot\varphi^2\right)+\frac{1}{2}I_3\left(\cos\theta\dot\varphi+\dot\psi\right)^2+mgl\cos\theta=c_3,
\end{array}\right.\end{align}</math>
</li>
<li> C.B.柯凡律夫斯卡亚（C.B.Kovaleski）情况：<math>I_1=I_2=2I_3</math>，重心在惯量椭球赤道面上。</li>
</ol>
</li>
<li> 拉莫尔进动：匀强磁场中，转动带电物体的角动量随时间的变化满足
<center><math>\begin{align}\frac{{\rm d}\boldsymbol{J}}{{\rm d}t}=-\frac{e\boldsymbol{B}}{2m}\times\boldsymbol{J}\end{align}</math></center>
与刚体角动量随时间的变化公式 
<math>\begin{align} \frac{{\rm d}\boldsymbol{J}}{{\rm d}t}=\frac{{\rm d}^*\boldsymbol{J}}{{\rm d}t}+\boldsymbol{\omega}\times\boldsymbol{J}.\end{align}</math> 
比较，可知带电物体的角动量大小不变，以角速度 <math>\boldsymbol{\omega}_l\equiv \frac{e\boldsymbol{B}}{2m}</math> 转动（绕磁场方向转动，类似于陀螺的进动），称为拉莫尔进动。</li>
</ol>

== 转动参考系 ==
=== 知识要点 ===

----

<ol>
<li>平面转动参考系（坐标原点与静止坐标系相同）： <math>\begin{align}\left\{\begin{array}{l}
\boldsymbol{v}=\boldsymbol{v}^\prime+\boldsymbol{\omega}\times\boldsymbol{r},\\
\boldsymbol{a}=\boldsymbol{a}^\prime+\dot{\boldsymbol{\omega}}\times\boldsymbol{r}-\omega^2\boldsymbol{r}+2\boldsymbol{\omega}\times\boldsymbol{v}^\prime\\
~~=\boldsymbol{a}^\prime+\boldsymbol{a}_{\rm t}+\boldsymbol{a}_{\rm c}
\end{array}\right.\end{align}</math></li>
<li>空间转动参考系（坐标原点与静止坐标系相同）：<math>\begin{align}\left\{\begin{array}{l}
\boldsymbol{v}=\boldsymbol{v}^\prime+\boldsymbol{\omega}\times\boldsymbol{r},\\
\boldsymbol{a}=\boldsymbol{a}^\prime+\dot{\boldsymbol{\omega}}\times\boldsymbol{r}+\boldsymbol{\omega}(\boldsymbol{\omega}\cdot\boldsymbol{r})-\omega^2\boldsymbol{r}+2\boldsymbol{\omega}\times\boldsymbol{v}^\prime\\
~~=\boldsymbol{a}^\prime+\boldsymbol{a}_{\rm t}+\boldsymbol{a}_{\rm c}
\end{array}\right.\end{align}</math><br/>
 <math>\frac{{\rm d}\boldsymbol{G}}{{\rm d}t}=\frac{{\rm d}^*\boldsymbol{G}}{{\rm d}t}+\boldsymbol{\omega}\times\boldsymbol{G},</math></li>
<li>非惯性系动力学：<math>\boldsymbol{a}_{\rm t}</math> 和 <math>\boldsymbol{a}_{\rm c}</math> 的影响等价于 <math>-m\boldsymbol{a}_{\rm t}</math>（惯性力）和 <math>-m\boldsymbol{a}_{\rm c}</math>（科里奥利力）的作用，即
<center><math>\begin{align}\left\{\begin{array}{ll}
m\boldsymbol{a}^\prime=\boldsymbol{F}-m\dot{\boldsymbol{\omega}}\times\boldsymbol{r}+m\omega^2\boldsymbol{r}-2m\boldsymbol{\omega}\times\boldsymbol{v}^\prime,&（平面转动参考系）\\
m\boldsymbol{a}^\prime=\boldsymbol{F}-m\dot{\boldsymbol{\omega}}\times\boldsymbol{r}+m\omega^2\boldsymbol{R}-2m\boldsymbol{\omega}\times\boldsymbol{v}^\prime,&（空间转动参考系）
\end{array}\right.\end{align}</math></center>
</li>
<li>应用
<ol>
<li>地球自转对地表物体运动造成的影响</li>
<li>傅科摆</li>
<li>地漏涡旋方向与地球自转有无关联？</li>
<li>考虑表面完全被水覆盖的行星，自转对行星表面的水波有无影响？对行星尺度的水波有无影响？</li>
<li>木星是巨型气体行星，表面有大红斑。试分析木星自转在大红斑形成和演化过程中的影响？</li></ol></li>
</ol>

== 分析力学 ==
=== 变分法 ===
变分：函数整体的微小变化，<math>f\to f+\delta f</math><br/>
变分问题例子：
* 最速降线问题： <math>\begin{align}I=\int_{x_A}^{x_B}\frac{\sqrt{1+y^{\prime 2}}}{\sqrt{2gy}}{\rm d}x\end{align}</math>
* 最小曲面问题：<math>\begin{align}I=\int_{x_A}^{x_B}y\sqrt{1+y^{\prime 2}}{\rm d}x \end{align}</math>
* 等周问题，等等。<br/>
变分问题：给定关于函数 <math>y(x)</math> 的定积分 <math>\begin{align}I=\int_{x_{A}}^{x_{B}} f(y,y^\prime){\rm d}x\end{align}</math>，求解函数形式变化时积分的极小值。
类比于普通函数极小值的求法，<math>I</math> 取极小值的条件为变分值为 0，即 <math>y\to y+\delta y</math>，有 <math>\delta I=0</math>。由此可以得到<span style="color: blue">'''变分极值满足的条件——欧拉-拉格朗日方程'''</span>：
<center><math>\begin{align} \frac{\rm d}{{\rm d}x}\frac{\partial f}{\partial y^\prime}+\frac{\partial f}{\partial y}=0.\end{align}</math></center>
这里有一些[[变分计算例子]]可供练习。

=== 虚功原理 ===
<ol>
<li> 力学系统平衡条件为 <math>\begin{align} \delta W=\sum_{i=1}^n \boldsymbol{F_i}\cdot\delta\boldsymbol{r}_i\end{align}</math>，亦可表示为 <math> \begin{align} Q_\alpha =0,\hspace{1ex}\left(\alpha=1,\ldots,s\right)\end{align}</math>，其中 <math>s=n-m</math> 是系统的自由度数，<math>m</math> 是约束的个数。</li>

<li>拉格朗日未定乘子法：在约束 <math>\begin{align} f_\beta(x_1,\ldots,x_n)=0,(\beta=1,2,\ldots,m)\end{align}</math> 限制下求解具有 <math>n</math> 个自变量的函数 <math>g(x_1,\ldots,x_n)</math> 的极值，等价于求解具有 <math>n+m</math> 个自变量的函数
<center><math>\begin{align}\tilde g(x_1,\ldots,x_n;\lambda_1,\ldots,\lambda_m)=g(x_1,\ldots,x_n)+\sum_{\beta=1}^m \lambda_\beta f_\beta(x_1,\ldots,x_n)\end{align}</math></center>
的极值。极值的必要条件，即对所有变量的一阶偏导为0，为<br/>
<math>\begin{align}\left\{\begin{array}{l}
\frac{\partial \tilde g}{\partial x_\alpha}=\frac{\partial g}{\partial x_\alpha}+\sum_{\beta=1}^m \lambda_\beta \frac{\partial f_\beta}{\partial x_\alpha}=0,\hspace{3ex}(\alpha=1,\ldots,n,)\\
\frac{\partial\tilde g}{\partial\lambda_\beta}=f_\beta=0,\hspace{3ex}(\beta=1,\ldots,m)
\end{array}\right.\end{align}</math>
</li>
</ol>

=== 拉格朗日方程 ===
<ol>

<li>哈密顿原理：两点之间粒子的运动轨迹是作用量取极值的轨迹，即<math>\delta S=0</math><br/>
上式中作用量 <math>\begin{align}S\equiv\int_{t_1}^{t_2}L(x,\dot x,t){\rm d}t=0\end{align}</math>，拉格朗日函数<math>L=T-V</math>.<br/>
根据哈密顿原理，
<center><math>\begin{align}
\delta S=&\int_{t_A}^{t_B} \delta L(x,\dot x,t){\rm d}t\\
=&\int_{t_A}^{t_B} \left[\frac{\partial L}{\partial x}\delta x+\frac{\partial L}{\partial\dot x}\delta\dot x\right]{\rm d}t\\
=&\int_{t_A}^{t_B} \left[\frac{\partial L}{\partial x}\delta x+\frac{\partial L}{\partial\dot x}\frac{\rm d}{{\rm d}t}\delta x\right]{\rm d}t\\
=&\int_{t_A}^{t_B} \left[\frac{\partial L}{\partial x}\delta x+\frac{\rm d}{{\rm d}t}\left(\frac{\partial L}{\partial\dot x}\delta x\right)-\frac{\rm d}{{\rm d}t}\left(\frac{\partial L}{\partial\dot x}\right)\delta x\right]{\rm d}t\\
=&\int_{t_A}^{t_B} \left[\frac{\partial L}{\partial x}-\frac{\rm d}{{\rm d}t}\left(\frac{\partial L}{\partial\dot x}\right)\right]\delta x{\rm d}t+\left.\left(\frac{\partial L}{\partial\dot x}\delta x\right)\right|_{t_A}^{t_B}\\
=&\int_{t_A}^{t_B} \left[\frac{\partial L}{\partial x}-\frac{\rm d}{{\rm d}t}\left(\frac{\partial L}{\partial\dot x}\right)\right]\delta x{\rm d}t
\end{align}</math></center>
因为<math>\delta x</math>是任意满足<math>\delta x(t_A)=\delta x(t_B)=0</math>的光滑函数，因此他与另一函数的积分值为零意味着这个函数处处为0，因此我们得到欧拉-拉格朗日方程
<center><math>\begin{align}\frac{\rm d}{{\rm d}t}\left(\frac{\partial L}{\partial \dot x}\right)-\frac{\partial L}{\partial x}=0.\end{align}</math></center></li>

<li>拉格朗日方程：如果系统运动由 <math>s</math> 个独立的广义坐标 <math> q_1,\ldots,q_s</math> 描述，则拉格朗日函数为<math> L(q_1,\ldots,q_s;\dot q_1,\ldots,\dot q_s;t)</math>，应用哈密顿原理，根据广义坐标的线性无关性，得到 <math>s</math> 个拉格朗日方程。此为保守系的拉格朗日方程：
<center><math>\begin{align}\frac{\rm d}{{\rm d}t}\left(\frac{\partial L}{\partial q_\alpha}\right)-\frac{\partial L}{\partial q_\alpha}=0,\hspace{4ex}\left(\alpha=1,\ldots,s\right)\end{align}</math></center>
上式中广义坐标 <math> q_\alpha</math>、广义速度 <math>\dot q_\alpha</math>。另外广义动量 <math>p_\alpha=\begin{align}\frac{\partial T}{\partial \dot q_\alpha}\end{align}</math>、广义力 <math>\begin{align} Q_\alpha=-\frac{\partial V}{\partial q_\alpha}=\sum_{i=1}^n\boldsymbol{F}_i\cdot\frac{\partial\boldsymbol{r}_i}{\partial q_\alpha}\end{align}</math>。将<math>L=T-V</math>带入保守系的拉格朗日方程，得到基本形式的拉格朗日方程：
<center><math>\begin{align}\frac{\rm d}{{\rm d}t}\left(\frac{\partial T}{\partial \dot q_\alpha}\right)-\frac{\partial T}{\partial q_\alpha}=Q_\alpha,~~~~~(\alpha=1,2,\ldots,s)\end{align}</math></center>
基本形式的拉格朗日方程可以应用到非保守系。
</li>

<li>拉格朗日方程的运动积分
# 如果 <math>\frac{\partial L}{\partial q_\beta}=0</math>，则<math>q_\beta</math>是系统的循环坐标，对应的广义动量是守恒量即<math>p_\beta\equiv \frac{\partial L}{\partial \dot q_\beta}=C</math>
# <math>T=T_0+T_2</math>，两项分别是广义速度的零次幂和二次幂。如果拉格朗日函数不显含时间，即 <math>\frac{\partial L}{\partial t}=0</math>，那么我们有能量积分 <math> T_2-T_0+V=h </math>，其中 <math>h</math> 不代表总能量；进一步，如果 <math>T_0=0</math>，则有能量积分 <math>T+V=E</math>，能量守恒。
上式分别代表广义的动量守恒和能量守恒律。</li>
</ol>

=== 小振动 ===
力学系统在平衡位置附近做微小振动，运动方程必为简谐振动或若干简谐振动的叠加。
# 一维情况，平衡位置意味着势函数在 <math>x_0</math> 处有 <math>V^\prime(x_0)=0,~~~V^{\prime\prime}(x_0)>0 </math>。将势函数对小参数 <math>\xi\equiv x-x_0</math> 作泰勒展开，得到势函数的近似形式 <math>V(x)\approx V_0+\frac{1}{2}k\xi^2,~~~(k\equiv V^{\prime\prime}(x_0)>0)</math>，因此运动方程变为 <math> m\ddot \xi+k\xi=0</math>，即简谐振动。
# 多维运动的质点，平衡位置附近，有 <math>\begin{align}\left\{\begin{array}{l}
V=\frac{1}{2}\sum_{\alpha,\beta=1}^s c_{\alpha\beta }q_\alpha q_\beta,\\
T=\frac{1}{2}\sum_{\alpha,\beta=1}^s a_{\alpha\beta} \dot q_\alpha\dot q_\beta,
\end{array}\right.\end{align}</math>，解出<math>\begin{align} q_\beta=\sum_{l=1}^s \Delta_{1\beta}(-\nu_l^2)\cos(\nu_l t+\varepsilon_l)\end{align}</math>
#;简正坐标：坐标变换，将 <math>V</math> 和 <math>T</math >对角化，解出的独立振动的坐标 <math>\xi_l=C_l\cos(\nu_lt+\varepsilon_l)</math> 叫做简正坐标，对应的角频率 <math>\nu_l</math>叫做简正频率。
#物理学中大量讨论的是系统在平衡位置附近的微振动，因此均可以应用上述小振动近似，简谐振动<span style="color:blue">'''普适的'''</span>存在于物理学研究的各种系统中。

=== 哈密顿方程 ===
<ol>
<li> 勒让德变换：<math>f(x)</math> 对 <math>x</math> 的勒让德变换为 <math>\bar f(u)=ux(u)-f\big(x(u)\big),~~~\left(u=\frac{{\rm d} f(x)}{{\rm d}x}\Longrightarrow x=x(u)\right)</math><br/>
哈密顿函数：<math>\begin{align} H(p,q,t)=\sum_{i=1}^s p_\alpha\dot q_\alpha-L(q_1,\ldots,q_s;\dot q_1,\ldots,\dot q_s;t)\end{align}</math><br/>
哈密顿正则方程：<math>\begin{align}\left\{\begin{array}{l}
\dot q_\alpha=\frac{\partial H}{\partial p_\alpha},\\
\dot p_\alpha=-\frac{\partial H}{\partial q_\alpha},
\end{array}\right.~~~~(\alpha=1,\ldots,s)\end{align}</math>
式中<math>q_\alpha,p_\alpha</math>叫做正则变量。</li>

<li> 泊松括号：<math>\begin{align} [A,B]\equiv\sum_{\alpha=1}^s\left(\frac{\partial A}{\partial q_\alpha}\frac{\partial B}{\partial p_\alpha}-\frac{\partial A}{\partial p_\alpha}\frac{\partial B}{\partial q_\alpha}\right)\end{align}</math>。<br/>
对于依赖于正则变量的函数 <math>\varphi(q_1,\ldots,q_s;p_1,\ldots,p_s;t)</math>，利用正则方程可以将其随时间的变化率表示为 <math>\begin{align}\frac{{\rm d}\varphi}{{\rm d}t}=\frac{\partial \varphi}{\partial t}+[\varphi,H]\end{align}</math><br/>
关于运动积分：
* 如果哈密顿函数不依赖于某个坐标 <math>q_\beta</math>，即 <math>\frac{\partial H}{\partial q_\beta}=0</math>，那么对应的正则动量守恒 <math>p_\beta=c_\beta,</math>；
* 如果哈密顿函数不显含时间 <math>\frac{\partial H}{\partial t}=0</math>，取 <math>\varphi=H</math>，根据 <math>\frac{{\rm d}\varphi}{{\rm d}t}</math> 的表达式，易见哈密顿函数守恒，即 <math>H=h</math>是守恒量；
* 对于一般的函数<math>\varphi</math>，如果它不显含时间即<math>\frac{\partial \varphi}{\partial t}=0</math>，且他与哈密顿函数的泊松括号为0即 <math>\left[\varphi, H\right]=0</math>，则此函数也是守恒量 <math>\varphi=C</math>
* 泊松定理：如果 <math>\varphi</math> 和 <math>\psi</math> 是正则方程的积分，则 <math>[\varphi,\psi]</math> 也是正则方程的积分。</li>

<li> 例：电磁场中的带电粒子，哈密顿函数为 <math>L=\frac{1}{2}mv^2-q\varphi+q\boldsymbol{A}\cdot\boldsymbol{v}</math>，正则动量依定义有 <math>\boldsymbol{p}\equiv\frac{\partial L}{\partial \boldsymbol{v}}=m\boldsymbol{v}+q\boldsymbol{A}</math>，勒让德变换得到哈密顿函数 <math>H=\boldsymbol{p}\cdot\boldsymbol{v}-L=\frac{1}{2}mv^2+q\varphi=\frac{1}{2m}\left(\boldsymbol{p}-q\boldsymbol{A}\right)^2+q\varphi</math>。<br/>
注意区分两种动量：普通的<span style="color:blue">'''机械动量为'''</span> <math>m\dot{\boldsymbol{r}}=\boldsymbol{p}-q\boldsymbol{A}</math>，仅与质点的运动方程 <math>\boldsymbol{r}(t)</math> 有关；<span style="color:blue">'''正则动量'''</span>为 <math>\boldsymbol{p}=m\dot{\boldsymbol{r}}+q\boldsymbol{A}</math>，与质点的运动方程 <math>\boldsymbol{r}(t)</math> 和电磁场矢量势 <math>\boldsymbol{A}(\boldsymbol{r},t)</math> 均有关。
</ol>

=== 哈密顿方程的一般求解方法：正则变换 & 哈密顿-雅可比理论 ===

<ol>

<li> 正则变换：保证哈密顿方程形式不变的坐标变换；<br/>
正则变换条件为 
<center><math>\begin{align} \sum_{\alpha=1}^s\left(p_\alpha{\rm d}q_\alpha-P_\alpha{\rm d}Q_\alpha\right)+(H^*-H){\rm d}t={\rm d}U_1\end{align}</math> </center>
上式也可以变形为
<center><math>\begin{align} \sum_{\alpha=1}^s\left(-q_\alpha{\rm d}p_\alpha-P_\alpha{\rm d}Q_\alpha\right)+(H^*-H){\rm d}t={\rm d}U_2,\hspace{3ex}\left(U_2=U_1-\sum_{\alpha=1}^s \left(p_\alpha q_\alpha\right)\right)\end{align}</math> </center>
<center><math>\begin{align} \sum_{\alpha=1}^s\left(p_\alpha{\rm d}p_\alpha+Q_\alpha{\rm d}P_\alpha\right)+(H^*-H){\rm d}t={\rm d}U_3,\hspace{3ex}\left(U_3=U_1+\sum_{\alpha=1}^s \left(P_\alpha Q_\alpha\right)\right)\end{align}</math> </center>
上式中<math>U_2,U_3</math>均为<math>U_1</math>的勒让德变换。<br/>
如果正则形式和 <math>U</math> 均不显含 <math>t</math>，则正则变换条件变为
<center> <math>\begin{align} \sum_{\alpha=1}^s\left(p_\alpha{\rm d}p_\alpha-P_\alpha{\rm d}Q_\alpha\right)={\rm d}U\end{align}</math>，且 <math>H^*=H</math>。</center>
正则变换的关键在于确保变换后的力学系统含有尽可能多的循环积分。<br/>
最理想的情况，<math>\begin{align} H^*=H^*(P_1,P_2,\ldots,P_\alpha,t) \end{align}</math>，则方程完全解出：<math>P_\alpha=</math>常数，<math>Q_\alpha=\int\frac{\partial H^*}{\partial P_\alpha}{\rm d}t,~~~(\alpha=1,2,\ldots,s)</math></li>

<li>'''哈密顿-雅可比方程'''：选取母函数 <math>S=U</math> 使得变换后的哈密顿函数 <math>H^*=0</math>，则母函数的满足的方程为 <math>\begin{align}\frac{\partial S}{\partial t}+H\left(t;q_1,q_2,\ldots,q_s;\frac{\partial S}{\partial q_1},\frac{\partial S}{\partial q_2},\ldots,\frac{\partial S}{\partial q_s}\right)=0\end{align}</math><br/>
如果哈密顿函数不显含 <math>t</math> 且约束为稳定的，则 <math>\begin{align}S=-Et+W(q_1,q_2,\ldots,q_s,\alpha_1,\alpha_2,\ldots,\alpha_s,E)+C\end{align}</math>，其中 <math>W</math> 叫做哈密顿-雅可比特性函数。哈密顿-雅可比方程变为
<center><math>\begin{align}H\left(q_1,q_2,\ldots,q_s,\frac{\partial W}{\partial q_1},\frac{\partial W}{\partial q_2},\ldots,\frac{\partial W}{\partial q_s}\right)=E,\end{align}</math></center></li>

<li>分离变量法求解哈密顿-雅可比方程：如果哈密顿函数<math>H</math>对坐标<math>q_\alpha</math>及偏导数<math>\frac{\partial W}{\partial q_\alpha}</math>的依赖可以放到一起，即
<center><math>\begin{align}H\left(q_1,\ldots,q_{\alpha-1},q_{\alpha+1},\ldots,q_s;\frac{\partial W}{\partial q_1},\ldots,\frac{\partial W}{\partial q_{\alpha-1}},\frac{\partial W}{\partial q_{\alpha+1}},\ldots,\frac{\partial W}{\partial q_s};\varphi_\alpha\left(q_\alpha,\frac{\partial W}{\partial q_\alpha}\right)\right)=E\end{align}</math></center>
则我们可以得到
<center><math>\begin{align}\left\{\begin{array}{l}
\varphi_\alpha\left(q_\alpha,\frac{\partial W}{\partial q_\alpha}\right)=c_\alpha,\\
 H\left(q_1,\ldots,q_{\alpha-1},q_{\alpha+1},\ldots,q_s;\frac{\partial W}{\partial q_1},\ldots,\frac{\partial W}{\partial q_{\alpha-1}},\frac{\partial W}{\partial q_{\alpha+1}},\ldots,\frac{\partial W}{\partial q_s};c_\alpha \right)=E
\end{array}\right.\end{align}</math></center>
<math>c_\alpha</math>是特征值。继续对哈密顿函数的其他变量应用分离变量法，可以依次解出方程。
</li>
<li> 哈密顿-雅可比方程是经典力学与量子力学的桥梁。薛定谔在建立量子力学的基本方程——哈密顿方程——的时候，就是从哈密顿-雅可比方程出发进行讨论的。详见薛定谔的论文集《Collected papers on wave mechanics》.
</ol>

=== 刘维尔定理 ===
质点的相空间分布函数<math>\rho(q_1,\ldots,q_s,p_1,\ldots,p_s,t)</math>随时间的变化率满足
<center><math>\begin{align}\frac{{\rm d}\rho}{{\rm d}t}=\frac{\partial\rho}{\partial t}+\left[\rho,H\right]=0.\end{align}</math></center>

=== 弯曲空间的拉格朗日方程与广义相对论 ===

<ol>

<li> 三维欧式空间中曲面和曲线上质点的运动：<br/>
曲面上，质点位矢参数方程为 <math>\boldsymbol{r}(u,v)</math>，速度为 <math>\boldsymbol{v}=\dot{\boldsymbol{r}}=\boldsymbol{r}_u\dot u+\boldsymbol{r}_v\dot v </math>，动能为 
<center><math>\begin{align} T=\frac{1}{2}m\boldsymbol{v}\cdot\boldsymbol{v}=\frac{1}{2}m\left[\boldsymbol{r}_u\cdot\boldsymbol{r}_u\dot u^2+2\boldsymbol{r}_u\cdot \boldsymbol{r}_v\dot u\dot v+\boldsymbol{r}_v\cdot \boldsymbol{r}_v\dot v^2\right]=\frac{1}{2}m\sum_{\alpha,\beta=1}^2 g_{ij}\dot q^i \dot q^j,\end{align} </math></center>
其中 <math>\begin{align}g_{ij}=\left(\begin{array}{cc}
\boldsymbol{r}_u\cdot \boldsymbol{r}_u & \boldsymbol{r}_u\cdot \boldsymbol{r}_v\\
\boldsymbol{r}_u\cdot \boldsymbol{r}_v & \boldsymbol{r}_v\cdot  \boldsymbol{r}_v
\end{array}\right)\end{align}</math>，是弯曲空间的度规张量；<math>q^1=u,~~~q^2=v</math>，是弯曲空间的坐标。<br/>
弯曲空间中两矢量的点积定义为 <math>\begin{align}\boldsymbol{A}\cdot\boldsymbol{B}=\sum_{i,j=1}^2 g_{ij}A^i B^j \end{align}</math>，因此动能可以写为
<center><math>\begin{align} T=\frac{1}{2}m\dot{\boldsymbol{q}}\cdot\dot{\boldsymbol{q}}=\frac{1}{2}m\sum_{i,j=1}^2 g_{ij}\dot q^i \dot q^j \end{align}</math></center>
如果质点除约束力外不受其他外力，则势能 <math>V=0</math>。将拉格朗日函数 <math>L=T</math> 代入拉格朗日方程<math>\begin{align}\frac{\rm d}{{\rm d}t}\left(\frac{\partial L}{\partial\dot q^\alpha}\right)-\frac{\partial L}{\partial q^\alpha}=0\end{align}</math>中，得到运动方程
<center><math>\begin{align} 
\frac{\rm d}{{\rm d}t}\left(\sum_{j=1}^2 g_{\alpha j}\dot q^j\right)-\frac{1}{2}\sum_{i,j=1}^2\frac{\partial }{\partial q^\alpha}\left(g_{ij}\right)\dot q^i \dot q^j=0,\\
2\sum_{i,j=1}^2\frac{\partial g_{\alpha j}}{\partial q^i}\dot q^i\dot q^j+2 \sum_{j=1}^2 g_{\alpha j}\ddot q^j-\sum_{i,j=1}^2\dot q^i\dot q^j\frac{\partial g_{ij}}{\partial  q^\alpha}=0\\
 \sum_{i,j=1}^2\left(\frac{\partial g_{\alpha j}}{\partial q^i}\dot q^i+\frac{\partial g_{\alpha i}}{\partial q^j}\dot q^i\right)\dot q^j+2\sum_{j=1}^2g_{\alpha j}\ddot q^j-\sum_{i,j=1}^2\dot q^i\dot q^j\frac{\partial g_{ij}}{\partial q^\alpha}=0\\
  \sum_{i,j=1}^2\left(\frac{\partial g_{\alpha j}}{\partial q^i}+\frac{\partial g_{\alpha i}}{\partial q^j}-\frac{\partial g_{ij}}{\partial  q^\alpha}\right)\dot q^i\dot q^j+2\sum_{j=1}^2g_{\alpha j}\ddot q^j=0
\end{align}</math></center>
定义 <math>g_{ij,k}\equiv\frac{\partial g_{ij}}{\partial q^k}</math>，且两边同乘以 <math>\frac{1}{2}g^{\alpha\beta}</math>，上式变为
<center><math>\begin{align}
  \sum_{i,j,\alpha=1}^2\frac{1}{2}g^{\alpha\beta}\left(g_{\alpha j,i}+ g_{\alpha i,j}- g_{ij,\alpha}\right)\dot q^i\dot q^j+\sum_{j,\alpha=1}^2g^{\alpha\beta}g_{\alpha j}\ddot q^j=0.
\end{align}</math></center>
根据联络的定义 <math>\begin{align}
  \Gamma^\alpha_{\mu\nu}\equiv\frac{1}{2}g^{\alpha\beta}(g_{\mu\beta,\nu}+g_{\nu\beta,\mu}-g_{\mu\nu,\beta}),
\end{align}</math>
以及度规张量的性质
<math>\begin{align}
  \sum_{\nu=1}^2g^{\mu\nu}g_{\nu\rho}=\delta^\mu_\rho.
\end{align}</math>
上式变为
<center><math>\begin{align}
  \sum_{i,j=1}^2\Gamma^\beta_{ij}\dot q^i\dot q^j+\ddot q^\beta=0,
\end{align}</math></center>
此即“测地线方程”，是不受主动力的质点约束在平面上自由运动的质点的运动方程。
</li>

<li>四维弯曲时空中质点的运动方程：<br>
广义相对论中，质点在四维弯曲时空中运动，按照二维曲面上质点运动的讨论，质点动能仍然可以写为
<center><math>\begin{align} T=\frac{1}{2}m\dot{\boldsymbol{q}}\cdot\dot{\boldsymbol{q}}=\frac{1}{2}m\sum_{\mu,\nu=1}^4 g_{\mu\nu}\dot q^\mu \dot q^\nu \end{align}</math></center>
上式中度规张量 <math> g_{\mu\nu}</math> 是四维时空的度规张量，根据广义相对论的基本方程，它的取值是由空间物质分布决定的。广义相对论将引力的作用效果等价于时空弯曲，质点在引力场中的运动等价于质点在弯曲时空中的<span style="color:blue">'''自由运动'''</span>，因此物体的运动方程仍为<span style="color:blue">'''测地线方程</span>，即
<center><math>\begin{align}
  \sum_{\mu,\nu=1}^4\Gamma^\beta_{\mu\nu}\dot q^\mu\dot q^\nu+\ddot q^\beta=0,
\end{align}</math></center>

</ol>

== 解题常见注意事项 ==
<ol>

<li> 一般解题步骤为：
<ol> <li> 想象运动图像 （理解题意、分析受力、应用守恒律）</li> <li> 根据运动三定理（动量定理、角动量定理、动能定理），列动力学方程 </li> <li> 求解方程</li></ol> </li>

<li> 无质量轻杆（或轻绳），如果仅受2个力作用，则两力必然沿杆（或绳）方向；搭在无质量滑轮上的绳子，两端作用力相等。</li>

<li> 质点或刚体在平衡位置附近作的微小振动，必为简谐振动或几个简谐振动的叠加。</li>

<li> 刚体上一点 B 的速度可以用刚体上另一点 A 的速度求出，即 <math> \boldsymbol{v}_B=\boldsymbol{v}_A+\boldsymbol{\omega}\times \boldsymbol{r}_{BA}</math>，其中 <math>\boldsymbol{r}_{BA}\equiv \boldsymbol{r}_B-\boldsymbol{r}_A</math>；刚体动能为 <math>T=\frac{1}{2}mv_c^2+\frac{1}{2}I_c\omega^2=\frac{1}{2}I_o\omega^2</math>，其中 <math>I_c</math> 和 <math>I_o</math> 分别为刚体绕质心和转动瞬心的转动惯量。</li>

<li> 区分<span style="color:blue">'''参考系'''</span>与<span style="color:blue">'''坐标系'''</span>
* 参照物（参考系）选定，质点的速度/角速度矢量大小和方向完全确定；坐标系仅仅决定了矢量要分解到哪几个方向上。
* 静止参考系中角速度与平动参考系中角速度相同，与转动参考系中角速度不同。例如，如果飞机有角速度，则飞机为参考物时角速度与地面为参照物时角速度不同。
* 固定在参照物上的坐标系（静止坐标系）与转动坐标系中，角速度相同。如刚体的角速度，一旦选取地面为参照物，刚体角速度就固定。固定在地面上的坐标系<math>o^\prime{\rm -}\xi\eta\zeta</math>和固定在刚体上的坐标系<math>o{\rm -}xyz</math>中研究物体的运动，区别仅仅是角速度做不同的坐标轴坐标轴<math>\hat{\boldsymbol{e}}_\xi,\hat{\boldsymbol{e}}_\eta,\hat{\boldsymbol{e}}_\zeta</math> 和 <math>\hat{\boldsymbol{i}},\hat{\boldsymbol{j}},\hat{\boldsymbol{k}}</math> 上分解，对应的分量不一样，即 <math>\boldsymbol{\omega}=\omega_\xi \hat{\boldsymbol{e}}_\xi+\omega_\eta\hat{\boldsymbol{e}}_\eta+\omega_\zeta\hat{\boldsymbol{e}}_\zeta=\omega_x\hat{\boldsymbol{i}}+\omega_y\hat{\boldsymbol{j}}+\omega_z\hat{\boldsymbol{k}}</math>.</li>

<li></li>
</ol>