- 分析动力学(第二版)
- 陈滨
- 3375字
- 2024-11-02 18:17:30
§2.1 广义坐标
2.1.1 完整约束组的区分
考虑N个质点组成的力学系统.设系统受有K个独立的一阶约束
缺秩为零.假定其中有L个约束是一阶线性的.不失一般性,可假定约束组(1.1)的前L个约束是一阶线性的.将这L个约束写成Pfaff形式:
假定我们在约束组(1.2)中找到了一个由d个约束组成的完全可积组.例如,假定
是一个完全可积组,即它满足1.2.3小节中的Frobenius定理的条件.于是这组约束等价(注:此处的“等价”是就约束组(1.3)和约束组(1.4)整组而言的,可参见1.2.5小节.)于如下的约束组
其中f1(u1,u2,…,u3N,t),f2(u1,u2,…,u3N,t),…,fd(u1,u2,…,u3N,t)是完全可积组(1.3)的积分函数.由此可见,(1.3)约束组等价于如下的d个独立的几何约束
(1.5)其中c1,c2,…,cd是由系统初始条件决定的常数.在规定系统一定要在确定的几何约束曲面上时,c1,c2,…,cd就是确定的常数.这时,系统的初始条件是受限制的,应该满足这些具有确定常数的约束方程组(1.5).因此,这和由初始条件来确定常数是一致的.我们以下研究的正是这种常见的情况.对于系统的完全可积约束组(1.3)仅为一阶约束组而无初始条件限制的情况,只是c1,c2,…,cd应看做未定的常参数,其他和常见的情况并无不同.
几何约束组(1.5)独立的条件是:在所关心的区域上Jacobi矩阵
的缺铁为零;亦即在所关心区域的每个子域上,都能找到点,使矩阵(1.6)的某d阶行列式不等于零.
以上叙述的是一般情况.比较简单的情形有:
(1)系统除完全可积的约束组(1.3)外,别无其他的约束了.这就是系统为完整系统的特殊情况.
(2)系统除完全可积的约束组(1.3)外,还有某些一阶线性约束,但没有非线性约束.这就是系统的约束仅为一阶线性约束的特殊情况.
2.1.2 广义坐标
我们从一般性的意义上来建立广义坐标.
不失一般性,假定在我们所关心的区域上,函数组(1.5)就是对u1,u2,…,ud这一组d个坐标满足
此时,我们引入Descartes位形空间C:{c=[u1,u2,…,u3N]T}到另一新的空间X:{x=[x1,x2,…,x3N]T}之间的变换,其变换关系式为
xs=fs(u1,u2,…,u3N,t),s=1,2,…,3N, (1.8)
其中前d个函数就是(1.5)式的几何约束表达式,而后3N-d个函数可以任选,只要求保证函数组(1.8)在我们关心的区域上是一个无关组,亦即使
根据隐函数定理,我们可以从关系式(1.8)反解得到
us=us(x1,…,xd,xd+1,…,x3N,t),s=1,2,…,3N. (1.12)
这样,关系式(1.8)和(1.12)就建立了C空间和X空间之间在我们所关心的区域上互相变换的一一映射关系.按照这样的变换,将已知的几何约束组(1.5)变到X空间里,并成为系统在X空间里必须满足的约束关系式,即
x1=c1,x2=c2,…,xd=cd. (1.13)
在满足(1.13)式这组几何约束的条件下,变换关系式(1.12)成为
us=us(c1,…,cd,xd+1,…,x3N,t),s=1,2,…,3N, (1.14)
其中c1,c2,…,cd为常数.我们记
xd+1=q1,xd+2=q2,…,x3N=qn, (1.15)
其中n=3N-d.这样,由公式(1.14)就得到在满足约束条件(1.5),亦即满足约束(1.13)的情况下,系统的3N个Descartes位形变元u1,u2,…,u3N用n个变元q1,q2,…,qn表达的关系式
us=us(q1,…,qn,t),s=1,2,…,3N. (1.16)
注意到,对于满足几何约束条件(1.5)这一点来说,变换关系式(1.16)是自动地达到,而不必再对q1,q2,…,qn有其他的限制了.这是变数q1,q2,…,qn的重要特点.同时,变换式(1.16)的Jacobi矩阵
(注:①注意,这个论断只在“我们所关心的区域”上成立.).
如上所得到的n个变数q1,q2,…,qn,我们称之为Lagrange广义坐标.对于广义坐标,除如上从变换的意义来理解而外,还可以这样来解释它的几何意义:在Descartes位形空间C中看,所有满足几何约束(1.5)的位形点c实际上组成了系统在(1.5)约束下的可达子空间.在一组确定的常数c1,c2,…,cd情况下,这个子空间是3N-d=n维的.公式(1.16)表明,广义坐标q1,q2,…,qn正是这n维可达子空间某个局部区域上的n个独立的描述参数,或者说是张在这n维可达子空间某个局部区域上的n维曲线坐标.无论从哪一个观点,系统满足几何约束的任意位形(限制在某个局部区域内)都可以用这n个广义坐标表示出来.反之,任意一组广义坐标数值也对应着系统满足几何约束的一个确定位形.由此,类似于Descartes位形空间C一样,我们记
q=[q1,q2,…,qn]T,
并引入由这n个广义坐标张成的n维空间(注:有关这个空间的度量我们将在需要的地方再赋予.它可能是欧氏的,但更一般的情况则是Riemann的.),来表现系统的位形(指满足几何约束并在某局部区域内的位形;以下未声明者,均同此义).我们称这个空间为广义坐标位形空间Cq,简称为广义坐标空间.
在这里,我们要着重作一个说明.对于一个力学系统,选定一组确定的广义坐标参数往往会导致位形描述不再具有全局性的限制.这个限制反映在要求[Jq]缺秩为零这项条件上.实际上,由于可达位形子空间具有一般微分流形的特征,一组确定的广义坐标参数一般都不能保证在Cq整个空间里处处有[Jq]缺秩为零.这时可能有奇点(或奇线,奇面),而奇点的邻域就是这组广义坐标失效的区域.我们可以举一个例子来说明.
一质点m约束在以O点为心,r为半径的球面上运动.系统的Descartes位形为
图 2.1
c=[u1,u2,u3]T=[x,y,z]T.
这是一个完整系统,受到的约束为
f=[x2+y2+z2-r2]/2=0.
系统的位形可达子空间是由f=0所决定的球面,这是一个二维的子空间.当我们想转移到广义坐标空间来做研究时,我们可以选用球面上的广义坐标q1,q2来描述系统的可达位形.例如,通常选q1为经度角φ,q2为纬度角θ,如图2.1所示.这时,变换关系(1.16)为
u1=x=rcosθcosφ,u2=y=rcosθsinφ,u3=z=rsinθ,
因而
由此不难看到,若要[Jq]的缺秩为零,其充要条件是θ≠±π/2;也就是说,对于上述选定的经纬度坐标来说,球面的上、下两个极点是奇点,它们的邻域是这组广义坐标失效的区域.
图 2.2
球面上经纬度坐标出现的上述缺陷并不是偶然的,也不是我们选择不当所引起的.实际上,由于球面和二维欧氏空间的不同胚性,球面上任何一组广义坐标都必然不是全局的,而只能是局部的.当然,奇点的数目和位置分布会随着广义坐标的选择而变化.例如,上述经纬度坐标是双奇点的.但是,在球面上可以找到单奇点的广义坐标系统.例如,图2.2所示的坐标网就是球面上单奇点的广义坐标.虽然由这组广义坐标决定的变换关系式要比经纬度坐标复杂,但由于奇点数目减少,这种坐标系统在地球球面导航的应用中是有价值的.
上述的广义坐标局部性质并没有在实用上引起重大的困难.在作了上述说明之后,今后我们在谈到整个广义坐标空间的时候,就不再对这些不大的奇点邻域作一一的说明了.
力学系统任一满足几何约束的运动都表现为q=[q1,q2,…,qn]T随时间变化的一维轨道.反过来,Cq空间里任一以时间为变元的一维轨道也对应系统满足几何约束的一种“可能运动”.这种一维轨道叫做系统运动的q轨迹.从这里可以看到,仅从满足几何约束来说,系统运动的q轨迹在Cq空间里是完全自由的了.
2.1.3 广义速度与广义加速度
在考虑了系统的完整约束组(1.5)之后,引入广义坐标q1,q2,…,qn,并且有表达式
us=us(q1,q2,…,qn,t),s=1,2,…,3N. (1.21)
系统的任一运动可表达为广义坐标q1,q2,…,qn随时间而变化,即有
qi=qi(t),i=1,2,…,n,
2.1.4 其他约束
按照2.1.1小节中的假定,对于我们所研究的系统除已经考虑的几何约束组(1.5)以外,还有两个独立的约束组:
一个是一阶线性约束组
在没有考虑这两组约束时,系统的运动在Cq空间里是自由的.现在要考虑这两组约束,我们需要把这两组约束条件转移到Cq空间里去.首先考虑(1.28)式的线性约束组.由(1.22)式得到
将其代入到(1.29)式中,立即得到
这样,我们就得到了附加在广义坐标轨道上两组约束条件(1.31)和(1.34),其中(1.31)式仍然是一阶线性约束,(1.34)式仍然是一般性的一阶约束,而且其独立性不变,即
的缺秩为零.
2.1.5 微变线性空间的变换
当应用(1.8)式的变换
xs=fs(u1,u2,…,u3N,t),s=1,2,…,3N
或者(1.12)式所表达的逆变换式
us=us(x1,x2,…,x3N,t),s=1,2,…,3N,
将C空间变到X空间时,C空间里原来的约束组(1.1)也变换到X空间里去,成为相应的X空间里的约束组,即
在假定约束组(1.1)中有一个几何约束组(1.5)情况下,根据2.1.2小节中的分析,已知前d个约束在X空间里的表达式为
x1=c1,x2=c2,…,xd=cd,
因而有
(1.45)这K-d个限制方程可分为两组,一组是由线性约束(1.31)式产生的,为
其中us和qi的关系式就是C空间到Cq空间的变换式us=us(q1,q2,…,qn,t),s=1,2,…,3N.