打开文本图片集
摘要:不同于传统流体力学,在Lagrange坐标下推导浅水波方程.若将水平位移作为基本变量,则推导出的浅水波数学模型可描述为固体力学的非线性大位移问题.运用不可压缩条件,通过变分原理推导出位移法浅水波方程,给出椭圆函数形式的行波解,并分析孤波解产生的条件.该基础研究建立了在分析结构力学中分析浅水波问题的理论基础,有利于进一步开展水动力学的研究.
关键词:
浅水波; Lagrange坐标; 孤立波; 椭圆函数; 微分代数方程; 保辛; 分析结构力学
中图分类号: O353.2
文献标志码: 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 引 言
中国沿海有广阔的大陆架及浅海区,面积较小的渤海平均水深只有25 m,主要由浅海构成,但现阶段中国总体海洋开发程度和利用率都很低.发展中国的海洋事业就需要进行浅海的开发研究.浅海的波浪直接影响渔业、运输和勘探作业等.因为浅海几十米的深度相对于广阔的海洋延伸尺度而言是很小的,所以其波动主要以浅水波为主.因此,对浅水波进行理论研究与分析对中国近海作业的发展有很大的推动作用.
对浅水波理论和方程[1-3]的研究历史悠久,但主要是在Euler坐标系下进行研究,其中著名的KdV浅水波方程就是在Euler坐标系下推导出的.浅水波的速度分布与水深无关,表明位移也与水深无关,这正满足杆件的刚性横截面假定.在Lagrange坐标系下,利用位移法可以推导出一套与传统的浅水波KdV方程完全不同的浅水波方程.于是,分析力学的变分原理即可运用正则变换、近似解的保辛积分等有效手段使数值求解得到很多方便.
1 KdV方程及其行波解
可以看出:式(27)与KdV方程得出的孤波解式(6)比较相似,只是在相同水深和波高的情况下KdV方程的孤波解的波速大于位移法所给出的波速.还应当指出,纯位移法提供过多的约束,从线性系统本征值的变分原理可知,纯位移约束总是使本征值单向提高,即结构的刚度提高.
3 位移法浅水波方程与传统KdV方程的差别
位移法浅水波方程和KdV方程都存在椭圆余弦波的周期解,而且也都导出孤波解,那么位移法方程与传统KdV方程在推导上存在什么不同,解又有什么差异呢?
首先,线性近似的结果是一致的,都支持水平速度沿高度不变的假定;其次,位移法考虑大变形.根据分析力学Lagrange体系的基本要求,运动学约束的不可压缩条件是严格满足的,然后根据变分原理导出动力方程.传统KdV方程是从速度势φ的Laplace方程div(grad φ)=0,即散度为0,进行φ的小参数展开而推导的.零次项近似支持线性理论,一次项自动为0,二次项导出KdV方程.grad φ=u意味着速度严格保证无旋条件,但无旋条件是动力学的条件,意味着无切力.小参数展开意味着散度为0的条件是近似满足的,也就是说运动学条件被放松了.
于是,很明显2种推导的差别在于:位移法推导严格满足运动学(协调)条件,于是可运用分析结构力学,沿着变分原理的路子走,当然保辛[11-13];进一步对其非线性微分方程的积分可考虑保辛等基本要求,这是其重要的优点;传统KdV方程的推导则放松运动学条件,部分满足动力学条件,是否保辛尚需探讨.
既然位移法浅水波方程的推导和积分运算可以保辛,那么相对传统的KdV方程又存在什么优势?下面就通位移法浅水波方程与KdV方程的孤波行波解所存在的差异讨论2种方程的优缺点.
KdV方程从产生至今已经有120多年的历史,虽然被广泛使用至今,但其孤波解存在2个主要不足:1)KdV方程只有单向运动(向左或向右)的孤子解,而真实的孤波可以存在同时向2个方向传播的多孤波解;2)KdV方程孤波解的波速比实际试验观测到的波速要大.
现在已经有很多方程可以用于描述同时向不同方向传播的孤波解,比如Boussinesq方程.位移法方程同样也可获取不同方向传播的孤波解,克服KdV方程的第一个不足之处,但很多浅水波方程的孤波解依然存在第二个问题.孤波的发现者SCOTT曾用重锤落入水槽的方法进行孤波的试验研究.他还从试验得出,孤波移动速度Vo与水槽中静止水深h和孤波波幅η0之间存在如下关系[6]98-103
这与本文的位移法浅水波的孤波解式(27)完全一致.从而,位移法浅水波方程也可克服KdV方程所存在的第二个不足之处,说明位移法浅水波方程的孤波解较其他方程存在理论优势.
4 结束语
运用Lagrange坐标将一维浅水波表达为固体力学的保守体系杆件纵向振动非线性问题.从而可用分析结构力学[11,13]的方法求解,Hamilton体系保辛积分等许多手段便都可应用.在分析浅水波问题时并不遵循传统Euler坐标的浅水波分析途径,而是另辟蹊径,将分析力学、固体力学与流体力学有机结合,运用保守体系的分析特点,在物质坐标系下完成浅水波方程的推导.位移法浅水波方程和KdV浅水波方程都存在椭圆余弦波的周期解,可利用精细积分进行数值分析计算.2种方程也都存在孤波解,且当波高远小于水深的情况下2种方程所给出的数值结果非常接近.在分析位移法浅水波方程与传统的KdV方程在推导上的差别之后,根据2种方程的孤波解的差异指出位移法浅水波方程的优势所在.
参考文献:
[1]COURANT R, FRIEDRICHS K O. Supersonic flow and shock waves[M]. New York: John Wiley & Sons, 1948.
[2] STOKER J J. Water waves[M]. New York: John Wiley & Sons, 1957.
[3] KINNMARK I. The shallow water wave equations: formulation, analysis and application[M]. Berlin: Springer-Verlag, 1986.
[4] GARDNER C S,GREENE J M,KRUSKAL M D, et al. Method for solving the Korteweg-de Vries Equation[J]. Physic Review Letters, 1967, 19(19): 1095-1097. DOI: 10.1103/PhysRevLett.19.1095.
[5] REMOISSENET M. Waves called solitons: Concepts and Experiments[M]. Berlin: Springer, 1996: 62-64.
[6] 倪皖荪, 魏荣爵. 水槽中的孤波[M]. 上海: 上海科技教育出版社, 1997: 7-8.
[7] 郭柏灵, 庞小峰. 孤立子[M]. 北京: 科学出版社, 1987: 5-25.
[8] 钟万勰, 姚征. 位移法浅水孤立波[J]. 大连理工大学学报, 2006, 46(1): 151-156.
ZHONG W X, YAO Z. Shallow water solitary waves based on displacement method[J]. Journal of Dalian University of Technology, 2006, 46(1): 151-156.
[9] 钟万勰, 姚征. 椭圆函数的精细积分算法[C]//应用力学进展论文集. 北京: 科学出版社, 2004: 106-111.
[10] 姚征, 钟万勰. 椭圆函数的精细积分改进算法[J]. 数值计算与计算机应用, 2008, 29(4): 251-260.
YAO Z, ZHONG W X. The improved precise integration method for elliptic functions[J]. Journal on Numerical Methods and Computer Applications, 2008, 29(4): 251-260.
[11] 钟万勰. 分析结构力学与有限元[J]. 动力学与控制学报, 2004, 2(4): 1-8.
ZHONG W X. Analytical structural mechanics and finite element[J]. Journal of Dynamics and Control, 2004, 2(4): 1-8.
[12] 钟万勰, 高强. WKBJ近似保辛吗?[J]. 计算力学学报, 2005, 22(1): 1-7.
ZHONG W X, GAO Q. Is the WKBJ approximation symplectic conservative?[J]. Chinese Journal of Computational Mechanics, 2005, 22(1): 1-7.
[13] 钟万勰, 高强, 彭海军. 经典力学——辛讲[M]. 大连: 大连理工大学出版社, 2013.
(编辑 武晓英)