Pochhammer频散方程高阶纵模式波近似解.pdfVIP

Pochhammer频散方程高阶纵模式波的近似解 崔志文张精刘金霞王克协 (吉林大学物理学院声学与微波物理系长春130021) An formulation for order modesinanelastic approximate high longitudinal cylinder CUIZhi-wenZHANGLIUJin—xia Jing WANGKe—xie 1．引言 考虑～半径为a的无限大各向同性的弹性圆柱体。弹性圆柱体的纵、横波速度分别用CL和CT表示。分别 用(-0和k表示角频率和轴向波数。沿应力自由表面的圆柱体轴向传播的纵模式波的频敖方程为”，2l (1) 它是著名的Pochhammer频散方程，是PocllhaIIlm目于1876年首先给出的。其中， c：-k (2) P=√∞2 2，g=、／印2c；一k2 定的实波数k，可以找到频率方程的无限个根，这些根代表在圆柱体中传播波型的频率。由此可绘出频率 波数曲线或频率相速度曲线。关于最低阶纵模式波的研究广泛的出现在研究生教材‘1捌中，很少讨论高阶纵 模式波。最近，Abroad[31研究表明，在速度位于纵横波之间的高阶模式波频散特性可由一个简单的近似公 线间的误差进一步减小。 2．高阶纵模式波的近似解 我们只考虑相速度在纵横波速度之间的实根。用相速度c和无因次波数“表示具有频率国和波数k的 则方程(1)可写为 粼+磊制中。厶∽¨0口)。蟛(1一口：r厶0口)r” (3) 其中埘∞五分别为0阶和1阶第一类变形Bessel函数。假定“，“口，加远大于1，保留到二阶小-量l／u2项 有，。∽)／，。∽)～l+1／(2“卢)+15／(8≯∥)2，代入(3)式有 量厶 0～缸^J砖 心A+咐15盯]，， 4=》耥卜 若-Io，以取保留到二阶小量l／甜2项，我们可得到如下近似的表达式“1 拙卜围[·一南]cos(一41留1 in(一41} 删～压{[1+☆M…三)去c。p詈]} 代入(5)整理得 3(3A 9r)] m =1、X-auln丽』函t百一面一，， ＼ ／“ ＼1二o u-|P，“J 。4～(彳一詈0)去U+陋+器]击tu 方鲫，解近似为铲吉(…署+tan4刁 48A+9y ， 7 蚍2≤蒜3 咖扣p扣毳A 班吉(仃石+三+一，]叫即是F‰妒给出的结果。 3．近似解与数值解的比较 图1是近似解与数值解的比较，其中实线是公式(7)的结果，点线是Ahrnad给出结果，虚线是数值 解。计算中取横波速度为2000m／s，泊松比为0．355。可以看出近似解给出很好的近似。图2表明利用本文 近似解比Ahmad给出结果相速度误差降低。以s，模式为例，可以看出本文给出的结果最大误差控制在2％以 下。模式波序数越高近似程度越好。 参考文献 【1]黎在良，刘殿魁，固体中的波[M】，北京：科学出版社，1995．262-266． J