达朗贝尔原理 (英语:d'Alembert principle )是因其发现者法国 物理学家与数学家让·达朗贝尔 而命名。达朗贝尔原理阐明,对于任意物理系统,所有惯性力 或施加的外力,经过符合约束条件 的虚位移 ,所作的虚功 的总和等于零[ 1] :
∑
i
(
F
i
+
I
i
)
⋅
δ
r
i
=
0
{\displaystyle \sum _{i}\ (\mathbf {F} _{i}+\mathbf {I} _{i})\cdot \delta \mathbf {r} _{i}=0\,\!}
;
本条目中,矢量 与标量 分别用粗体 与斜体 显示。例如,位置矢量通常用
r
{\displaystyle \mathbf {r} \,\!}
表示;而其大小则用
r
{\displaystyle r\,\!}
来表示。
达朗贝尔
其中,
I
i
=
−
m
i
a
i
{\displaystyle \mathbf {I} _{i}=-m_{i}\mathbf {a} _{i}\,\!}
是粒子
P
i
{\displaystyle P_{i}\,\!}
感受到的惯性力,
m
i
{\displaystyle m_{i}\,\!}
和
a
i
{\displaystyle \mathbf {a} _{i}\,\!}
分别是粒子
P
i
{\displaystyle P_{i}\,\!}
的质量 和加速度 ,
F
i
{\displaystyle \mathbf {F} _{i}\,\!}
是施加于粒子
P
i
{\displaystyle P_{i}\,\!}
的外力(不包括约束力 )、
δ
r
i
{\displaystyle \delta \mathbf {r} _{i}\,\!}
是符合系统约束 的虚位移。
静力学 的虚功原理 在动力学 的版本是达朗贝尔原理。假若一个物理系统的每一个约束条件都只约束位置 或时间,而不约束速度 ,则称此物理系统为完整系统 。达朗贝尔原理比哈密顿原理 的适用范围更广阔,可以用于不仅是完整系统。
因为达朗贝尔原理,在一个动力系统里,约束力所作的虚功自动抵消,也就是说,不需要顾虑约束力所作的虚功。
思考由一群粒子构成的一个物理系统。按照牛顿运动定律 [ 2] ,
F
i
(
T
)
=
m
i
a
i
{\displaystyle \mathbf {F} _{i}^{(T)}=m_{i}\mathbf {a} _{i}\,\!}
;
其中,
F
i
(
T
)
{\displaystyle \mathbf {F} _{i}^{(T)}\,\!}
是所有施加于粒子
P
i
{\displaystyle P_{i}\,\!}
的作用力的合力 (包括约束力)。
将方程右边的加速度项目移至左边,
F
i
(
T
)
−
m
i
a
i
=
0
{\displaystyle \mathbf {F} _{i}^{(T)}-m_{i}\mathbf {a} _{i}=\mathbf {0} \,\!}
。
达朗贝尔建议将这加速度项目视为一种因为粒子的运动而产生的作用力,称为惯性力 :
I
i
=
−
m
i
a
i
{\displaystyle \mathbf {I} _{i}=-m_{i}\mathbf {a} _{i}\,\!}
。
这样,施加于每一个粒子的作用力(包括惯性力)的矢量和皆等于零:
F
i
(
T
)
+
I
i
=
0
{\displaystyle \mathbf {F} _{i}^{(T)}+\mathbf {I} _{i}=\mathbf {0} \,\!}
。
采用达朗贝尔这绝顶聪明的建议,这系统内所有的作用力的矢量和变为零,也就是说,这系统达到平衡状态。假若动力系统的动态平衡可以视为静力系统的静态平衡,则所有静力系统内有关于平衡状态的理论都可以适用于动力系统,而这动力系统的运动问题的一大部分也可以当作静力系统的平衡问题来解析。因此,当然也可以将静力学的虚功原理搬迁至动力学里。
对于每一个粒子,经过虚位移
δ
r
i
{\displaystyle \delta \mathbf {r} _{i}\,\!}
,其矢量和所作的虚功等于零:
δ
W
i
=
(
F
i
(
T
)
+
I
i
)
⋅
δ
r
i
=
0
{\displaystyle \delta W_{i}=(\mathbf {F} _{i}^{(T)}+\mathbf {I} _{i})\cdot \delta \mathbf {r} _{i}=0\,\!}
。
作用于每一个粒子的虚功的总和
δ
W
{\displaystyle \delta W\,\!}
等于零:
δ
W
=
∑
i
δ
W
i
=
∑
i
(
F
i
(
T
)
+
I
i
)
⋅
δ
r
i
=
0
{\displaystyle \delta W=\sum _{i}\delta W_{i}=\sum _{i}(\mathbf {F} _{i}^{(T)}+\mathbf {I} _{i})\cdot \delta \mathbf {r} _{i}=0\,\!}
。
将作用于每一个粒子上的合力
F
i
(
T
)
{\displaystyle \mathbf {F} _{i}^{(T)}\,\!}
,细分为外力
F
i
{\displaystyle \mathbf {F} _{i}\,\!}
与约束力
C
i
{\displaystyle \mathbf {C} _{i}}
:
δ
W
=
∑
i
F
i
⋅
δ
r
i
+
∑
i
C
i
⋅
δ
r
i
+
∑
i
I
i
⋅
δ
r
i
=
0
{\displaystyle \delta W=\sum _{i}\ \mathbf {F} _{i}\cdot \delta \mathbf {r} _{i}+\sum _{i}\ \mathbf {C} _{i}\cdot \delta \mathbf {r} _{i}+\sum _{i}\mathbf {I} _{i}\cdot \delta \mathbf {r} _{i}=0\,\!}
。
假设,每一个约束力,因为虚位移,所做的虚功的总和是零[ 3] 。则约束力的项目可以从方程内移去,达朗贝尔原理成立:
δ
W
=
∑
i
(
F
i
+
I
i
)
⋅
δ
r
i
=
0
{\displaystyle \delta W=\sum _{i}(\mathbf {F} _{i}+\mathbf {I} _{i})\cdot \delta \mathbf {r} _{i}=0\,\!}
。(1)
现在,总和内的每一个单独
F
i
−
m
i
a
i
{\displaystyle \mathbf {F} _{i}-m_{i}\mathbf {a} _{i}\,\!}
很可能不等于零。
定义有效力
F
i
e
f
f
{\displaystyle \mathbf {F} _{i}^{eff}\,\!}
为外力加惯性力:
F
i
e
f
f
=
F
i
+
I
i
{\displaystyle \mathbf {F} _{i}^{eff}=\mathbf {F} _{i}+\mathbf {I} _{i}\,\!}
。
达朗贝尔原理又可表达为:对于任意物理系统,所有有效力,经过符合约束条件 的虚位移 ,所作的虚功 的总合等于零,
δ
W
=
∑
i
F
i
e
f
f
⋅
δ
r
i
=
0
{\displaystyle \delta W=\sum _{i}\mathbf {F} _{i}^{eff}\cdot \delta \mathbf {r} _{i}=0\,\!}
。
注意到这推论里的约束力假设。在这里,约束力就是牛顿第三定律 的反作用力 。因此,可以称此假设为反作用力的虚功假设 :所有反作用力所做的符合约束条件的虚功,其总合是零。这是分析力学额外设立的假设,无法从牛顿运动定律 推导出来[ 1] 。
主项目:拉格朗日方程
拉格朗日力学 是对经典力学的一种不同的表述。拉格朗日方程是拉格朗日力学的基要方程,可以用来描述物体的运动,特别适用于理论物理的研究。拉格朗日方程的功能等价于牛顿力学 中的牛顿第二定律 。
从达朗贝尔原理,可以推导出拉格朗日方程[ 3] 。设定粒子
P
i
{\displaystyle P_{i}\,\!}
的位置
r
i
{\displaystyle \mathbf {r} _{i}\,\!}
为广义坐标
q
1
,
q
2
,
⋯
,
q
n
{\displaystyle q_{1},q_{2},\cdots ,q_{n}\,\!}
与时间
t
{\displaystyle t\,\!}
的函数:
r
i
=
r
i
(
q
1
,
q
2
,
⋯
,
q
n
,
t
)
{\displaystyle \mathbf {r} _{i}=\mathbf {r} _{i}(q_{1},q_{2},\cdots ,q_{n},\ t)\,\!}
。
转换为广义坐标的主要的目的,是要除去物体内粒子位置与粒子位置之间的相依性。这问题在后面会有更详细的说明。
虚位移可以表示为
δ
r
i
=
∑
j
∂
r
i
∂
q
j
δ
q
j
{\displaystyle \delta \mathbf {r} _{i}=\sum _{j}\ {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\delta q_{j}\,\!}
。(2)
粒子的速度
v
i
=
v
i
(
q
1
,
q
2
,
⋯
,
q
n
,
q
˙
1
,
q
˙
2
,
⋯
,
q
˙
n
,
t
)
{\displaystyle \mathbf {v} _{i}=\mathbf {v} _{i}(q_{1},q_{2},\cdots ,q_{n},{\dot {q}}_{1},{\dot {q}}_{2},\cdots ,{\dot {q}}_{n},\ t)\,\!}
是
v
i
=
d
r
i
d
t
=
∂
r
i
∂
t
+
∑
j
∂
r
i
∂
q
j
q
˙
j
{\displaystyle \mathbf {v} _{i}={\frac {d\mathbf {r} _{i}}{dt}}={\frac {\partial \mathbf {r} _{i}}{\partial t}}+\sum _{j}{\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}{\dot {q}}_{j}\,\!}
。
取速度对于广义速度的偏微分:
∂
v
i
∂
q
˙
j
=
∂
r
i
∂
q
j
{\displaystyle {\frac {\partial \mathbf {v} _{i}}{\partial {\dot {q}}_{j}}}={\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\,\!}
。(3)
思考方程(1)的加速度项目,将方程(2)代入,
∑
i
m
i
a
i
⋅
δ
r
i
=
∑
i
,
j
m
i
a
i
⋅
∂
r
i
∂
q
j
δ
q
j
{\displaystyle \sum _{i}\ m_{i}\mathbf {a} _{i}\cdot \delta \mathbf {r} _{i}=\sum _{i,j}\ m_{i}\mathbf {a} _{i}\cdot {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\delta q_{j}\,\!}
。
应用乘积法则 ,
∑
i
,
j
m
i
a
i
⋅
∂
r
i
∂
q
j
δ
q
j
=
∑
i
,
j
(
d
d
t
(
m
i
v
i
⋅
∂
r
i
∂
q
j
)
−
m
i
v
i
⋅
d
d
t
(
∂
r
i
∂
q
j
)
)
δ
q
j
{\displaystyle \sum _{i,j}\ m_{i}\mathbf {a} _{i}\cdot {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\delta q_{j}=\sum _{i,j}\left({\frac {d}{dt}}\left(m_{i}\mathbf {v} _{i}\cdot {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\right)-m_{i}\mathbf {v} _{i}\cdot {\frac {d}{dt}}\left({\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\right)\right)\delta q_{j}\,\!}
。
注意到
∂
r
i
∂
q
j
{\displaystyle {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\,\!}
的参数为
q
1
,
q
2
,
⋯
,
q
n
,
t
{\displaystyle q_{1},q_{2},\cdots ,q_{n},\ t\,\!}
,而速度
v
i
{\displaystyle \mathbf {v} _{i}\,\!}
的参数为
q
1
,
q
2
,
⋯
,
q
n
,
q
˙
1
,
q
˙
2
,
⋯
,
q
˙
n
,
t
{\displaystyle q_{1},q_{2},\cdots ,q_{n},{\dot {q}}_{1},{\dot {q}}_{2},\cdots ,{\dot {q}}_{n},\ t\,\!}
,所以,
d
d
t
(
∂
r
i
∂
q
j
)
=
(
∂
∂
t
+
∑
k
q
˙
k
∂
∂
q
k
)
(
∂
r
i
∂
q
j
)
=
∂
2
r
i
∂
q
j
∂
t
+
∑
k
∂
2
r
i
∂
q
j
∂
q
k
q
˙
k
{\displaystyle {\frac {d}{dt}}\left({\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\right)=\left({\frac {\partial }{\partial t}}+\sum _{k}{\dot {q}}_{k}{\frac {\partial }{\partial q_{k}}}\right)\left({\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\right)={\frac {\partial ^{2}\mathbf {r} _{i}}{\partial q_{j}\partial t}}+\sum _{k}{\frac {\partial ^{2}\mathbf {r} _{i}}{\partial q_{j}\partial q_{k}}}{\dot {q}}_{k}\,\!}
、
∂
v
i
∂
q
j
=
∂
∂
q
j
(
∂
r
i
∂
t
+
∑
k
∂
r
i
∂
q
k
q
˙
k
)
=
∂
2
r
i
∂
q
j
∂
t
+
∑
k
∂
2
r
i
∂
q
j
∂
q
k
q
˙
k
{\displaystyle {\frac {\partial \mathbf {v} _{i}}{\partial q_{j}}}={\frac {\partial }{\partial q_{j}}}\left({\frac {\partial \mathbf {r} _{i}}{\partial t}}+\sum _{k}{\frac {\partial \mathbf {r} _{i}}{\partial q_{k}}}{\dot {q}}_{k}\right)={\frac {\partial ^{2}\mathbf {r} _{i}}{\partial q_{j}\partial t}}+\sum _{k}{\frac {\partial ^{2}\mathbf {r} _{i}}{\partial q_{j}\partial q_{k}}}{\dot {q}}_{k}\,\!}
。
因此,以下关系式成立:
d
d
t
(
∂
r
i
∂
q
j
)
=
∂
v
i
∂
q
j
{\displaystyle {\frac {d}{dt}}\left({\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\right)={\frac {\partial \mathbf {v} _{i}}{\partial q_{j}}}\,\!}
。(4)
将方程(3)与(4)代入,加速度项目成为
∑
i
,
j
m
i
a
i
⋅
∂
r
i
∂
q
j
δ
q
j
=
∑
i
,
j
(
d
d
t
(
m
i
v
i
⋅
∂
v
i
∂
q
˙
j
)
−
m
i
v
i
⋅
∂
v
i
∂
q
j
)
δ
q
j
{\displaystyle \sum _{i,j}\ m_{i}\mathbf {a} _{i}\cdot {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\delta q_{j}=\sum _{i,j}\left({\frac {d}{dt}}\left(m_{i}\mathbf {v} _{i}\cdot {\frac {\partial \mathbf {v} _{i}}{\partial {\dot {q}}_{j}}}\right)-m_{i}\mathbf {v} _{i}\cdot {\frac {\partial \mathbf {v} _{i}}{\partial q_{j}}}\right)\delta q_{j}\,\!}
。
思考这个系统的动能
T
{\displaystyle T\,\!}
,
T
=
∑
i
1
2
m
i
v
i
⋅
v
i
{\displaystyle T=\sum _{i}\ {\frac {1}{2}}m_{i}\mathbf {v} _{i}\cdot \mathbf {v} _{i}\,\!}
。
加速度项目与动能的关系为
∑
i
,
j
m
i
a
i
⋅
∂
r
i
∂
q
j
δ
q
j
=
∑
j
(
d
d
t
(
∂
T
∂
q
˙
j
)
−
∂
T
∂
q
j
)
δ
q
j
{\displaystyle \sum _{i,j}\ m_{i}\mathbf {a} _{i}\cdot {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\delta q_{j}=\sum _{j}\ \left({\frac {d}{dt}}\left({\frac {\partial T}{\partial {\dot {q}}_{j}}}\right)-{\frac {\partial T}{\partial q_{j}}}\right)\delta q_{j}\,\!}
。(5)
思考方程(1)的外力项目,将方程(2)代入,
∑
i
F
i
⋅
δ
r
i
=
∑
i
,
j
F
i
⋅
∂
r
i
∂
q
j
δ
q
j
=
∑
j
F
j
δ
q
j
{\displaystyle \sum _{i}\ \mathbf {F} _{i}\cdot \delta \mathbf {r} _{i}=\sum _{i,j}\ \mathbf {F} _{i}\cdot {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\delta q_{j}=\sum _{j}\ {\mathcal {F}}_{j}\delta q_{j}\,\!}
;(6)
这里,
F
{\displaystyle {\boldsymbol {\mathcal {F}}}\,\!}
是广义力 :
F
j
=
∑
i
F
i
⋅
∂
r
i
∂
q
j
{\displaystyle {\mathcal {F}}_{j}=\sum _{i}\ \mathbf {F} _{i}\cdot {\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\,\!}
。
将方程(5)与(6)代入方程(1),会得到
∑
j
(
d
d
t
(
∂
T
∂
q
˙
j
)
−
∂
T
∂
q
j
−
F
j
)
δ
q
j
=
0
{\displaystyle \sum _{j}\ \left({\frac {d}{dt}}\left({\frac {\partial T}{\partial {\dot {q}}_{j}}}\right)-{\frac {\partial T}{\partial q_{j}}}-{\mathcal {F}}_{j}\right)\delta q_{j}=0\,\!}
。(7)
假设所有的广义坐标都相互独立,则所有的广义坐标的虚位移也都相互独立。由于这些虚位移都是任意设定的,只有满足下述方程,才能使方程方程(7)成立:
d
d
t
(
∂
T
∂
q
˙
j
)
−
∂
T
∂
q
j
−
F
j
=
0
{\displaystyle {\frac {d}{dt}}\left({\frac {\partial T}{\partial {\dot {q}}_{j}}}\right)-{\frac {\partial T}{\partial q_{j}}}-{\mathcal {F}}_{j}=0\,\!}
。(8)
假设这系统是单演系统 ,也就是说,这系统的广义力与广义位势
V
{\displaystyle V\,\!}
之间的关系式为
F
j
=
d
d
t
(
∂
V
∂
q
˙
j
)
−
∂
V
∂
q
j
{\displaystyle {\mathcal {F}}_{j}={\frac {d}{dt}}\left({\frac {\partial V}{\partial {\dot {q}}_{j}}}\right)-{\frac {\partial V}{\partial q_{j}}}\,\!}
,
那么,
d
d
t
(
∂
(
T
−
V
)
∂
q
˙
j
)
−
∂
(
T
−
V
)
∂
q
j
=
0
{\displaystyle {\frac {d}{dt}}\left({\frac {\partial (T-V)}{\partial {\dot {q}}_{j}}}\right)-{\frac {\partial (T-V)}{\partial q_{j}}}=0\,\!}
。
广义位势 也是系统的势能 。注意到拉格朗日量
L
{\displaystyle L\,\!}
定义为系统的动能减去势能:
L
=
d
e
f
T
−
V
{\displaystyle L\ {\stackrel {def}{=}}\ T-V\,\!}
,
则可得到拉格朗日方程:
d
d
t
(
∂
L
∂
q
˙
j
)
−
∂
L
∂
q
j
=
0
{\displaystyle {\frac {d}{dt}}\left({\frac {\partial L}{\partial {\dot {q}}_{j}}}\right)-{\frac {\partial L}{\partial q_{j}}}=0\,\!}
。
假设这系统是保守系统 ,也就是说,这系统的广义力与位势
V
{\displaystyle V\,\!}
之间的关系式为
F
j
=
−
∂
V
∂
q
j
{\displaystyle {\mathcal {F}}_{j}=-{\frac {\partial V}{\partial q_{j}}}\,\!}
,
则拉格朗日方程也成立。
根据对于刚体的牛顿第二定律 ,一个运动中的刚体,其运动方程为
∑
i
F
i
=
m
a
{\displaystyle \sum _{i}\ \mathbf {F} _{i}=m\mathbf {a} \,\!}
、
∑
i
M
i
=
I
α
{\displaystyle \sum _{i}\ \mathbf {M} _{i}={\boldsymbol {\mathcal {I}}}{\boldsymbol {\alpha }}\,\!}
;
其中,
F
i
{\displaystyle \mathbf {F} _{i}\,\!}
是施加于刚体的外力,
m
{\displaystyle m\,\!}
是刚体的质量,
a
{\displaystyle \mathbf {a} \,\!}
是刚体质心的加速度,
M
i
{\displaystyle \mathbf {M} _{i}\,\!}
是每一个外力
F
i
{\displaystyle \mathbf {F} _{i}\,\!}
对于刚体质心 的力矩 、
I
{\displaystyle {\boldsymbol {\mathcal {I}}}\,\!}
是对于刚体质心的惯性张量 ,
α
{\displaystyle {\boldsymbol {\alpha }}\,\!}
是刚体的角加速度 。
达朗贝尔建议将加速度项目
−
m
a
{\displaystyle -m\mathbf {a} \,\!}
视为一种因为刚体的运动而产生的作用力,称为惯性力
I
{\displaystyle \mathbf {I} \,\!}
,又将角加速度项目
−
I
α
{\displaystyle -{\mathcal {I}}{\boldsymbol {\alpha }}\,\!}
视为一种因为刚体的运动而产生的力矩,称为惯性力矩
M
{\displaystyle \mathbf {M} \,\!}
:
I
=
−
m
a
{\displaystyle \mathbf {I} =-m\mathbf {a} \,\!}
、
M
=
−
I
α
{\displaystyle \mathbf {M} =-{\boldsymbol {\mathcal {I}}}{\boldsymbol {\alpha }}\,\!}
。
那么,运动方程变为
I
+
∑
i
F
i
=
0
{\displaystyle \mathbf {I} +\sum _{i}\ \mathbf {F} _{i}=\mathbf {0} \,\!}
、
M
+
∑
i
M
i
=
0
{\displaystyle \mathbf {M} +\sum _{i}\ \mathbf {M} _{i}=\mathbf {0} \,\!}
。
在工程力学 里,达朗贝尔惯性力原理 阐明:刚体的惯性力与所有作用于刚体的外力的合力等于零,刚体的惯性力矩与所有作用于刚体的力矩的合力矩等于零。[ 4] 。这原理可以帮助分析正在运动中的某连杆所感受到的作用力。
请注意,惯性力必须作用于质心;而惯性力矩是力偶矩 ,可以作用于物体的任何一位置。靠着达朗贝尔惯性力原理,动力系统可以变为像静力系统一样的解析。这方法的优点是,在等价的静力系统里,可以选择任何一点(不只是质心)来计算力矩。这时常会导至较简易的运算。因为,如果选择出正确的力矩作用点,在计算力矩时,可以忽略许多作用力(这些作用力与选择点同直线)。
Lanczos, Cornelius, The Variational Principles of Mechanics, Dovers Publications, Inc: pp. 90–106, 1970, ISBN 978-0-486-65067-8
Goldstein, Herbert. Classical Mechanics 3rd. United States of America: Addison Wesley. 1980: pp. 18–21. ISBN 0201657023 (英语) .
Beer, Ferdinand; E. Russell Johnston, Jr., William E. Clausen. Vector Mechanics for Engineers 7th. United States of America: Elizabeth A. Jones. 2004: pp. 1029, 1167. ISBN 0-07-230491-X (英语) .