丁光涛
(安徽师范大学物理与电子信息学院,芜湖 241000)
根据给定的系统运动性质,确定系统的受力情况,是动力学一类基本问题,这类问题及其各种变形问题称为动力学逆问题;在质点力学中,已知质点的运动规律求作用其上的力,是已知作用在质点上的力求其运动规律问题的逆问题[1-3].变分法逆问题研究给定的运动微分方程能否从变分原理中导出,即能否构造出对应的Lagrange函数,将给定的方程表示成Lagrange方程形式[4-9].显然,上述两类逆问题可以组合起来,根据给定的运动性质来构造Lagrange函数.本文以质点阻尼运动为例,探讨上述两类逆问题的“灰色”特征,即解的不唯一性,以及两类问题的组合导致的新情况,对一个特定的运动存在多种不同的Lagrange函数,但是,在一般情况下这些函数可能不全部是等效的.
设质点质量为m,运动规律为
这个规律与线性阻尼运动规律一致,本文中的阻尼运动就是指按此规律的质点运动.下面导出作用在质点上的力的函数表达式.将式(1)对时间一次和二次导数,分别得到
利用牛顿运动定律,直接得到结果为
即质点受到有关大小随时间按指数规律衰减的阻力作用,其运动微分方程是
必须指出,导出的力的函数表达式不是唯一的,结合式(2)得到
即质点受到大小与速度成正比、方向与速度相反的阻尼力作用,其运动微分方程是
这是通常所说的阻尼运动微分方程.力还可以表示成坐标x的函数,结合式(1),得到
质点运动微分方程是
利用是(2)和(8),还能够得到
质点运动微分方程是
有些文献中认为质点动力学逆问题比正问题简单,在数学上是微分问题,从运动规律导出作用在质点上的(合)力随时间变化的函数关系即可.但是,从物理方面来看,问题并非如此,实际的作用力是与时间、位置和速度相关的,要求的是力随这些因素变化的函数关系,即得到力的变化规律(力律).但是,如果只是已知某个特定的运动规律,信息是不充分的,得到的解不是唯一的.上述导出的四种不同的作用力F的函数形式和运动微分方程,反映了逆问题的“灰色”特性,四种解存在共同点,都可以导出给定的运动规律,但是,在力律上是有区别的.对此给出进一步说明如下:
1)四种运动微分方程存在共同的特解,即可以导出相同的运动规律式(1).值得注意的是,这四个相同的特解都对应着同一种初始条件:
2)四种不同的作用力F的函数表达式中,式(4)的F1和式(8)的F3存在相同点,即表达式中都包含着初始条件v0;而式(6)的 F2和 式(10)的F4也存在相同点,即表达式中都不包含初始条件v0.这两种情况是有原则区别的,动力学逆问题的一般解应当是与初始条件无关的.
3)式(12)初始条件中,t和x可以任意赋值,可以不为0,但是,不同的初值相当于计时原点和坐标原点的移动,不会带来任何动力学的影响;但是,v不可任意赋值,因为在阻尼运动中,阻尼力的力律是大小与速度成正比,方向与速度相反,而这里的速度实质上是质点相对于媒质的速度,将质点相对于参考系的速度与相对于媒质的速度等同起来的条件是媒质相对于参考系静止,换句话说,v的初值不能通过引入运动参考系的变换来改变的.
4)因此,下面的讨论中,将舍弃力的表达式中参数与运动初始条件赋值紧密相关的方程(5)和(9),而只讨论方程(7)和(11).后两个方程可以看作阻尼运动微分方程,质点受到阻尼力作用,方程(7)中质点只受线性阻尼力作用,方程(11)中质点受到变系数阻尼力和有势力的合力作用,从物理方面来看,这两个方程所代表的系统是不同的.
构造Lagrange函数的方法有很多种,本文仅利用一种从运动方程直接构造Lagrange函数的方法[9].给定一维系统运动微分方程为
设其Lagrange函数为
代入Lagrange方程,展开得到
式中
将方程(14)改写成如下形式
比较方程(13),得到
根据 f(t,x,˙x)的具体形式,引入辅助条件,解方程(16),得到 u(t,x,˙x),v1(t,x)和 v0(t,x)从而构造出(15)式中函数L.应当指出,对一维系统利用规范等效变换,适当选择规范变换函数,总可以使函数v1或v0中的一个变换为零,从而简化求解.
在具体问题中可以对函数u的宗量作6种不同的设定,导出对应的特殊解法.
1.设 u=u(t),由(15)和(16)式得
2.设 u=u(x),由(15)和(16)式得
3.设 u=y(˙x),由(15)和(16)式得
4.设 u=u(t,x),由(15)和(16)式得
5.设 u=u(t,˙x),由(15)和(16)式得
6.设 u=u(x,˙x),由(15)和(16)式得
首先,讨论线性阻尼运动,其运动微分方程为(为简化取m=1)
应用特殊解法1,可得
应用特殊解法3,可得
从φ0=λ,可得到不同的解,如
v"1=,等等.
对应地,一组规范等效的Lagrange函数为
应用特殊解法 5,设 u=u(t,˙x),则有
方程(28)存在一个解族
其中 c0是常数,F=F(ξ)=F(˙xeλt+c0)是对其宗量ξ任意的光滑函数,但要求 d2F/dξ2≠0.与(29)式 u对应的是阻尼运动的一个Lagrange函数族[10,11]
当F取不同函数形式时,就得到一系列不同而等效的Lagrange函数,例如
等等.应当指出(26)式中L′也属于这个函数族.
应用特殊解法 6,设 u=u(x,˙x),可以求解,(25)式中L就是特解之一.应用一般解法,设u=u(t,x,˙x),也可求解,(27)式中 L"就是一个特解.
其次,计算变系数阻尼运动的Lagrange函数,其运动微分方程为(m=1)
这个方程比方程(23)复杂,但是,应用特殊解法1,仍然可以解出
即Lagrange函数为
这个Lagrange函数与线性阻尼运动的Lagrange函数没有等效关系,由它们导出的运动方程的解并不等价,仅仅在特殊的初始条件下有一个共同的特解(1).
本文以运动规律(1)表示的阻尼运动为例,讨论动力学逆问题、变分法逆问题和两者的组合问题的特点.
⑴.这个实例清楚地表明了动力学逆问题的“灰色”特点,已知信息不充分,导致结果不确定.即使对简单的一维运动而言,导出的几种力律也是不同的,运动微分方程不是等价的,对应著物理上不同的系统.
⑵.变分法逆问题的解,即使对同一个运动微分方程,导出的Lagrange函数也可能不是唯一的.线性阻尼运动的Lagrange函数和函数族说明了这个特点,但是,这些Lagrange函数是等效的.
⑶.两类逆问题是相关联的,它们的组合问题产生新的特点.由于动力学逆问题可以导出不相同且不等价的运动微分方程,使得变分法逆问题从不同的微分方程出发,构造得到的Lagrange函数更多,这些函数中有些是等效的,它们对应着同一个方程,有些函数分别对应着不等价的方程,它们之间没有通常的等效关系.
1 梅凤翔.动力学逆问题.北京:国防工业出版社,2009(Mei F X.Inverse problems of dynamics.Beijing:National Defense Industry Press,2009(in Chinese))
2 Rosenberg B M.Analytical dynamics of discrete systems.New York:Pleum Prss,1977
3 梅凤翔.分析力学(下卷).北京:北京理工大学出版社,2013(Mei F X.Analyticalmechanics(II).Beijing:Beijing Institute of Technology Press,2013(in Chinese))
4 Santilli R M.Foundations of theoreticalmechanics I.New York:Springer-Verlag,1978
5 Santilli R M.Foundations of theoreticalmechanics II.New York:Springer-Verlag,1983
6 Lopuszanski J.The inverse variational problem in classical mechanics.Singapore:World Scientific,1999
7 梅凤翔.分析力学专题,北京:北京工业学院出版社,1988(Mei FX.Special problems of analyticalmechanics.Beijing:Beijing Institute Technologe Press,1988(in Chinese))
8 Cielinski J L,Nikiciuk T.A direct approach to the construction of standard and non-standard Lagragians for Dissipative-like dynamical systems with variable coefficients.Journal of Physics A:Mathematical and Theoretical,2000,43:175205
9 丁光涛.从运动方程构造Lagrange函数的直接方法.动力学与控制学报,2010,8:305~310(Ding G T.A direct approach to the construction of Lagrangians from themotion equations.Journal of Dynamics and Control,2010,8:305~310(in Chinese))
10 丁光涛.关于一类Lagrange函数族的存在条件.动力学与控制学报,2011,9(3):219~221(Ding G T.On the existence conditions for a class of the Lagrangians families.Journal of Dynamics and Control,2011,9(3):219~221(in Chinese))
11 丁光涛.一类Painleve方程的Lagrange函数族.物理学报,2012,61:110202(Ding G T.The families of Lagrangians of a Painleve equation.Acta Physica Sinica,2012,61:110202(in Chinese) )