Mooney-Rivilin压缩杆模型的精确解及其动力学性质王婧【摘要】利用动力系统理论中的积分分支法研究Mooney-Rivilin压缩杆模型的行波解,获得了包括尖孤子解、爆破波解、周期波解、亮孤子解和暗孤子解等各种精确行波解.进一步讨论了这些精确行波解随时间演化的动力学行为,并利用Maple 软件绘出了具有代表性的精确行波解随时间演化的坐标图形.与现有文献的结果相比,本文所获得的精确解比较新颖.【期刊名称】《玉溪师范学院学报》【年(卷),期】2017(033)012【总页数】6页(P15-20)【关键词】积分分支法;精确行进波;孤立波解;尖孤子解【作者】王婧【作者单位】重庆师范大学数学科学学院,重庆 401331【正文语种】中文【中图分类】TQ330自然科学中的许多非线性问题,都可以用非线性偏微分方程和发展方程来建模,且非线性偏微分方程和发展方程的解可以用来解释模型的某些物理现象,特别是质点随时间演化的各种动力学行为和动力学性质.有时为了能够更加精准地了解模型所代表的事物内在变化的规律,我们就要去求其相应模型的精确解.目前,在精确求解非线性偏微分方程和发展方程方面已获得的一些有效的方法,如:反散射变换[1,2]、达布变换[3,4]、双线性方法和多元线性方法[5,6]、齐次平衡法[7~9]、Lie 群法[10,11]、平面动力系统分支理论[12]等.在文献[13,14]中,Rui等在平面动力系统分支理论[12]的基础上,提出了一个比较简化的方法叫做积分分支方法,该方法与平面动力系统分支理论中方法相比,绕开了相对复杂的平面像图分析过程,采取因式分解技巧与直接积分相结合的方法来获得非线性偏微分方程和发展方程的各种精确行波解并进一步讨论解的各种动力学行为和动力学性质,以此来解释模型所代表的问题中各种非线性物理现象.本文将利用积分分支法与因式分解相结合的方法来调查下列非线性Mooney-Rivilin压缩杆模型ut-uxxt+3uux-λ(2uxuxx+uuxxx)+βux=0(1)的各种精确解,并进一步讨论它们的动力学行为和动力学性质.方程(1)是戴辉辉从可压缩的超弹性材料中导出的弹性杆模型,又称Mooney-Rivilin压缩杆模型,其中u=u(x,t),而参数λ为材料参数依赖于预应力,γ为一常数[15].特别地,当λ=1时,方程(1)就可约化成著名的潜水波模型Camassa-Holm方程,曾经被许多研究人员研究过.同样,这个可压缩的超弹性杆模型(1)也被很多研究人员研究过它的行波的间断和爆破现象,详细报道见文献[16]中所引用的大量文献.2 模型的精确行波解及动力学性质为了获得方程(1)的各种精确行波解,我们做下列行波变换u(x,t)=φ(ξ),ξ=lx-ct(2)其中l为波数,c为波速度,二者均为非零实常数.将(2)式代入(1)式后整理得:(βl-c)φ′+cl2φ‴+3lφφ′-2λl3φ′φ″-λl3φφ‴=0,(3)其中φ′=dφ/dξ将(3)式积分一次得:g+(βl-c)φ-0.5λl3(φ′)2+(cl2-λl3φ)φ″+1.5lφ2=0,(4)其中g为积分常数.当φ≠c/(λl)时,方程(4)化为φ″=[g+(βl-c)φ-0.5λl3(φ′)2+1.5lφ2]/(λl3φ-cl2).(5)令φ′≡dφ/dξ=y,则方程(5)可化为下列平面动力系统:dφ/dξ=y,dy/dξ=[g+(βl-c)φ-0.5λl3y2+1.5lφ2]/(λl3φ-cl2),(6)系统(6)是一个奇异系统,因为当φ=c/(λl)时,dy/dξ无定义,由此我们称φ=c/(λl)为系统(6)的奇异直线.显然当φ=c/(λl)时,系统(6)与方程(4)不等价,然而,φ=c/(λl)也是方程(4)的一个解(一个常数解).为了获得与方程(4)完全等价的平面动力系统,我们做下列尺度变换dξ=(λl3φ-cl2)dτ,(7)其中τ为参数.则系统(6)在变换(7)下被化成一个规则的二维平面动力系统:=y(λl3φφ-0.5λl3y2+1.5lφ2.(8)显然无论变量φ如何变化,系统(8)始终与方程(4)等价.通过简单计算后不难发现系统(6)和(8)具有相同的首次积分,即0.5lφ3+0.5(βl-c)φ2-0.5λl3y2φ+0.5cl2y2+gφ=h,(9)其中h为新的积分常数.记函数H(φ,y)=0.5lφ3+0.5(βl-c)φ2-0.5λl3y2φ+0.5cl2y2+gφ=h.(10)容易验证函数H(φ,y)满足∂H/∂φ=dy/dτ,∂H/∂y=-dφ/dτ所以系统(8)是一个哈密顿系统,满足哈密顿守恒律.显然O(0,0)原点始终是系统(8)的其中一个平衡点.利用(10)式,不难计算出原点的Hamilton量为h0=H(0,0)=0.根据积分常数和哈密顿量是否为零的情况,下面讨论方程(1)的精确行波解.首先,由方程(9)可解得(11)情形1 当两个积分常数g=0,h=0时,将(11)式代入系统(8)的第一式后化成=±dτ(12)在计算的过程中,为方便起见,记(12)式等号右边的正负号为ε,那么在不同的参数条件下积分(12)后换回变量u,并将积分所得的结果分别代入变换式(7)中再次积分后联立u和ξ两式,我们便得到下列方程(1)的12种参数形式的精确行波解:(13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24)在以上这些参数形式的精确行波解中,解(13)(14)(15)为孤立波解,解(17)(18)(20)为周期波解,解(23)和(24)为有理解.为了直观地显示以上解的动力学行为,以其中部分解为例画出其坐标图,见图1.图1 参数条件g=0,h=0下各种参数型精确行波解的坐标图情形2 当g≠0,h=0,λ>0时,将(11)式代入(7)式和(8)式的第一个式子可得φ(25)其中φ3=c/(λl),φ4=0.当积分常数满足g=(c-βl)2/(8l)条件且c=λβl/(λ-2)时,则有φ1,2=c/(λl),那么(25)式可化成(26)完成(26)式的积分后换回变量,我们获得方程(1)的一个显函数形式的精确行波解:(27)当积分常数满足g=(c-βl)2/8l条件,但c≠λβl/(λ-2)时,则仅有φ1=φ2=(c-βl)/2l,那么(26)可化为φ(28)完成(30)式的积分后换回变量,我们获得方程(1)的一个隐函数形式的精确行波解:(29)为了直观地显示精确行波解(27)和(29)的动力学行为,我们画出了它们的坐标图,如图2所示.情形3 当积分常数g≠0,h≠0时,将(11)式直接代入(8)式中的第一个式子化简可得(30)(a)在特定的参数条件下(30)式可分解为φ(31)图2 由解(27)和(29)定义的U型波当φ1>φ2>φ3≥0>φ,完成(31)式积分后换回相应变量,可获得方程(1)的一种隐函数形式的解:其中为第三类椭圆积分函数且式中其他参数为(b)在特定的参数条件下,(30)又可以分解为φ(33)其中φ1>φ2≥φ>φ3>0.类似地,完成(33)式的积分后换回相应变量,我们得到方程(1)的另一个隐函数形式的精确行波解:(34)其中为第三类椭圆积分函数且式中其他参数为该文由芮伟国教授指导完成.参考文献:[1]V.E.Zakharov and A.B.Shabat,Exact theory of two-dimensional self-focusing and one dimensional self-modulation of waves in nonlinear media[J].JETP,1972(34):62-69.[2]C.S.Garder,J.M.Greene,M.D.Kruskal,R.M.Mirura,Method for solving the Korteweg-de -Vries equation[J].Physical Review Letters,1967(19):1095-1097.[3]C.S.Garder,J.M.Greene,M.D.Kruskal,R.M.Mirura,Korteweg-de Vries equation and generalizations[J].VI.Methods for exactsolution,Communications on Pure and Applied Mathematics,1974,27(1):97-133.[4]V.B.Matveev,M.A.Salle,Darboux transformations andSolitons[M].Berlin:Springer,1991.[5]I.T.Habibullin, Backlund transformation and integrable initial-boundary value problems[J].Mathematical notes of the Academy of Sciences of the USSR,1991,49(4):418-423.[6]R.Hirota,Exact solution of the Korteweg-de Vries equation for multiple collisions of solitons[J].Physical Review Letters,1971(27):1192-1194. [7]A.M.Wazwaz,The tanh method for traveling wave solutions of nonlinear equations[J].Applied Mathematics and Computation,2004,154(3):713-723.[8]M.Wang,Y.Zhou,Z.Li,Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics[J].Physics Letters A,1996,216(1):67-75.[9]E.Fan,H.Zhang,A note on the homogeneous balance method[J].Physics Letters A,1998,246(5):403-406.[10]E.D.Belokolos,A.I.Bobenko,Algebro-geometric approach to nonlinear integrable Equations[M].Berlin:Springer-Verlag,1994.[11]N.C.Freeman,J.J.C.Nimmo,Soliton solutions of the Korteweg de Vries and the Kadomtsev-Petviashvili equations: the Wronskiantechnique[J].Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences.The Royal Society,1977(389):319-329.[12]J.Li,Z.Liu,Smooth and non-smooth travelling waves in nonlinearly dispersive equation[J].Applied Mathematical Modeling,2000(25):41-56. [13]W.Ru,B.He,Y.Long,et al.,The integral bifurcation method and its application for solving a family of third-order dispersive PDEs[J].Nonlinear Analysis: Theory, Methods and Applications,2008,69(4):1256-1267.[14]W.Rui,B.He,S.Xie,et al.Application of the integral bifurcation method for solving modified Camassa-Holm and Degasperis-Procesiequations[J].Nonlinear Analysis:Theory, Methods andApplications,2009,71(7):3459-3470.[15]H.H.Dai,Model equations for nonlinear dispersive waves in a compressible Mooney-Rivlin Rod[J].Acta Mech.1998(127):193-207. [16]F.Ma,Y.Liu,C.Qu,Wave-breaking phenomena for the nonlocal Whitham-type equations[J].Journal of Differential Equations,2016(261):6029-6054.。