一、第二类拉格朗日方程#
一般地,质点系各个质点的矢量坐标可以表示为广义坐标的函数r=ri(q1,q2,⋯,qk;t)。质点速度相应地通过广义速度表示为:
r˙i=l=1∑k∂ql∂riq˙l+∂t∂r˙i上式两边对广义速度q˙j求偏导数,得到消点规则:
∂qj∂ri=∂q˙∂r˙质点速度关于广义坐标的导数为
∂qj∂r˙i=l=1∑k∂ql∂qj∂2riq˙l+∂t∂qj∂2ri而质点矢量坐标关于广义坐标的导数仍为广义坐标的函数,
∂qj∂ri=∂qj∂ri(q1,q2,⋯,qk;t)再将它关于时间求全导数,可得
dtd(∂qj∂ri)=l=1∑k∂ql∂(∂qj∂ri)q˙l+∂t∂(∂qj∂ri)=l=1∑k∂qj∂ql∂2riq˙l+∂qj∂t∂2ri这里,假定矢量坐标 ri 一阶与二阶连续可微,即得另一个恒等关系式,即交换规则:
dtd(∂qj∂ri)=∂qj∂r˙i于是可将广义惯性力表示为:
Qgj=−i=1∑nmidtdr˙i⋅∂qj∂ri=−i=1∑nmidtd(r˙i⋅∂qj∂ri)+i=1∑nmir˙i⋅dtd(∂qj∂ri)即
Qgj=−dtd(i=1∑nmir˙i⋅∂qj∂ri)+i=1∑nmir˙i⋅∂qj∂r˙i(2.10)注意到:
r˙i⋅∂q˙j∂r˙i=∂q˙j∂(21r˙i⋅r˙i)=∂q˙j∂(21vi2)r˙i⋅∂qj∂r˙i=∂qj∂(21r˙i⋅r˙i)=∂qj∂(21vi2)那么广义惯性力可以表示为:
Qgj=−dtd[i=1∑nmi∂q˙j∂(21vi2)]+i=1∑nmi∂qj∂(21vi2)=−dtd[∂q˙j∂(i=1∑n21mivi2)]+∂qj∂(i=1∑n21mivi2)=−dtd(∂q˙j∂T)+∂qj∂T其中,T=∑i=1n21mivi2 是系统的动能。因此,广义惯性力可以表示成动能的导数表达,将广义惯性力代入到广义坐标系中的动力学普遍方程中,得到:
dtd(∂q˙j∂T)−∂qj∂T=Qjj=1,2,⋯,k这就是(第二类)拉格朗日方程
二、广义速度表示的动能#
利用速度表达式,可以将质点系的动能表示为:
T=l=1∑n21mlvl2=l=1∑n21mlr˙l⋅r˙l=l=1∑n21ml(i=1∑k∂qi∂rlq˙i+∂t∂rl)⋅(j=1∑k∂qj∂rlq˙j+∂t∂rl)=21i=1∑kj=1∑kaijq˙iq˙j+i=1∑kbiq˙i+c=T2+T1+T0其中动能关于广义速度的二次项、一次项与零次项部分分别为:
T2=21i=1∑kj=1∑kaijq˙iq˙j,aij=l=1∑nml∂qi∂rl⋅∂qj∂rlT1=i=1∑kbiq˙i,bi=l=1∑nml∂qi∂rl⋅∂t∂rlT0=c=21l=1∑nml∂t∂rl⋅∂t∂rl这里的aij,bi,c都是广义坐标及时间的函数,aij=aji具有对称性。
三、保守系统的拉格朗日方程#
保守系统受到的主动力都是有势力,而有势力相应地广义力可以表示为势能的导数的负数,所以拉格朗日方程成为:
dtd(∂q˙j∂T)−∂qj∂T+∂qj∂V=0(j=1,2,⋯,k)势能仅仅取决于质点系的位形,只是广义坐标的函数,从而:
∂q˙j∂V=0引入拉格朗日函数(或者说动势):
L=T−V拉格朗日函数L=L(q˙1,q˙2,⋯,q˙k;q1,q2,⋯,qk),则拉格朗日方程可以表示为:
dtd(∂q˙j∂L)−∂qj∂L=0(j=1,2,⋯,k)上面括号中的内容叫做广义动,后一项为拉格朗日力。此时,只需计算系统的动能与势能,无需计算广义力。利用动能表达式,可以将拉格朗日函数表示为:
L=T2+T1+T0−V
四、非保守系统的拉格朗日方程#
对于非保守系统,主动力可以分为有势力和非有势力两类:
dtd(∂q˙j∂T)−∂qj∂T=Qj=−∂qj∂V+Q~j系统的拉格朗日方程为:
dtd(∂q˙j∂L)−∂qj∂L=Q~j(j=1,2,...,k)其中Q~是非有势力的广义力。