姚征+鐘萬勰

摘要：不同于傳統(tǒng)流體力學(xué)，在Lagrange坐標(biāo)下推導(dǎo)淺水波方程.若將水平位移作為基本變量，則推導(dǎo)出的淺水波數(shù)學(xué)模型可描述為固體力學(xué)的非線性大位移問題.運(yùn)用不可壓縮條件，通過變分原理推導(dǎo)出位移法淺水波方程，給出橢圓函數(shù)形式的行波解，并分析孤波解產(chǎn)生的條件.該基礎(chǔ)研究建立了在分析結(jié)構(gòu)力學(xué)中分析淺水波問題的理論基礎(chǔ)，有利于進(jìn)一步開展水動(dòng)力學(xué)的研究.

關(guān)鍵詞：

淺水波； Lagrange坐標(biāo)； 孤立波； 橢圓函數(shù)； 微分代數(shù)方程； 保辛； 分析結(jié)構(gòu)力學(xué)

中圖分類號(hào)： O353.2

文獻(xiàn)標(biāo)志碼： A

Abstract：The shallow water wave equation is obtained in Lagrange coordinates， which is different from the traditional hydrodynamics. The horizontal displacement is defined as the basic variable， and the deduced mathematical model of shallow water wave just likes the nonlinear large displacement problem in solid mechanics. Under the incompressible condition， the variational principle for shallow water dynamic equation is given on the basis of the displacement method， and the corresponding elliptic traveling wave solutions are presented the conditions of obtaining solitary wave solutions are also discussed. As a basic research work， the present study establishes the theoretical foundation of analyzing the shallow water waves in analytical structural mechanics and is valuable for the further research in hydrodynamics.

Key words：

shallow water wave； Lagrange coordinate； solitary wave； elliptic function； differentialalgebraic equation； symplectic constraint； analytical structural mechanics

0 引 言

中國(guó)沿海有廣闊的大陸架及淺海區(qū)，面積較小的渤海平均水深只有25 m，主要由淺海構(gòu)成，但現(xiàn)階段中國(guó)總體海洋開發(fā)程度和利用率都很低.發(fā)展中國(guó)的海洋事業(yè)就需要進(jìn)行淺海的開發(fā)研究.淺海的波浪直接影響漁業(yè)、運(yùn)輸和勘探作業(yè)等.因?yàn)闇\海幾十米的深度相對(duì)于廣闊的海洋延伸尺度而言是很小的，所以其波動(dòng)主要以淺水波為主.因此，對(duì)淺水波進(jìn)行理論研究與分析對(duì)中國(guó)近海作業(yè)的發(fā)展有很大的推動(dòng)作用.

對(duì)淺水波理論和方程[1-3]的研究歷史悠久，但主要是在Euler坐標(biāo)系下進(jìn)行研究，其中著名的KdV淺水波方程就是在Euler坐標(biāo)系下推導(dǎo)出的.淺水波的速度分布與水深無關(guān)，表明位移也與水深無關(guān)，這正滿足桿件的剛性橫截面假定.在Lagrange坐標(biāo)系下，利用位移法可以推導(dǎo)出一套與傳統(tǒng)的淺水波KdV方程完全不同的淺水波方程.于是，分析力學(xué)的變分原理即可運(yùn)用正則變換、近似解的保辛積分等有效手段使數(shù)值求解得到很多方便.

1 KdV方程及其行波解

可以看出：式（27）與KdV方程得出的孤波解式（6）比較相似，只是在相同水深和波高的情況下KdV方程的孤波解的波速大于位移法所給出的波速.還應(yīng)當(dāng)指出，純位移法提供過多的約束，從線性系統(tǒng)本征值的變分原理可知，純位移約束總是使本征值單向提高，即結(jié)構(gòu)的剛度提高.

3 位移法淺水波方程與傳統(tǒng)KdV方程的差別

位移法淺水波方程和KdV方程都存在橢圓余弦波的周期解，而且也都導(dǎo)出孤波解，那么位移法方程與傳統(tǒng)KdV方程在推導(dǎo)上存在什么不同，解又有什么差異呢？

首先，線性近似的結(jié)果是一致的，都支持水平速度沿高度不變的假定；其次，位移法考慮大變形.根據(jù)分析力學(xué)Lagrange體系的基本要求，運(yùn)動(dòng)學(xué)約束的不可壓縮條件是嚴(yán)格滿足的，然后根據(jù)變分原理導(dǎo)出動(dòng)力方程.傳統(tǒng)KdV方程是從速度勢(shì)φ的Laplace方程div（grad φ）=0，即散度為0，進(jìn)行φ的小參數(shù)展開而推導(dǎo)的.零次項(xiàng)近似支持線性理論，一次項(xiàng)自動(dòng)為0，二次項(xiàng)導(dǎo)出KdV方程.grad φ=u意味著速度嚴(yán)格保證無旋條件，但無旋條件是動(dòng)力學(xué)的條件，意味著無切力.小參數(shù)展開意味著散度為0的條件是近似滿足的，也就是說運(yùn)動(dòng)學(xué)條件被放松了.

于是，很明顯2種推導(dǎo)的差別在于：位移法推導(dǎo)嚴(yán)格滿足運(yùn)動(dòng)學(xué)（協(xié)調(diào)）條件，于是可運(yùn)用分析結(jié)構(gòu)力學(xué)，沿著變分原理的路子走，當(dāng)然保辛[11-13]；進(jìn)一步對(duì)其非線性微分方程的積分可考慮保辛等基本要求，這是其重要的優(yōu)點(diǎn)；傳統(tǒng)KdV方程的推導(dǎo)則放松運(yùn)動(dòng)學(xué)條件，部分滿足動(dòng)力學(xué)條件，是否保辛尚需探討.

既然位移法淺水波方程的推導(dǎo)和積分運(yùn)算可以保辛，那么相對(duì)傳統(tǒng)的KdV方程又存在什么優(yōu)勢(shì)？下面就通位移法淺水波方程與KdV方程的孤波行波解所存在的差異討論2種方程的優(yōu)缺點(diǎn).