无界弦振动的研究

无界弦振动的研究

马玉荣

摘 要 用行波法、积分变换法（傅里叶变换法、拉普拉斯变换法）、分离变量法、格林函数法求解无界弦的自由振动和受迫振动问题。计算和分析表明：对于无界弦的自由振动问题，行波法和傅里叶变换法比较简便，这是常用的求解方法。对于无界弦受迫振动问题，利用叠加原理应用行波法和齐次化原理求解最简便。行波法对于求解无界弦振动问题有其特殊的优点，即，行波法已求出无界弦自由振动问题的达朗贝尔公式，无界弦受迫振动问题的公式，这些公式是通用的，只要把具体问题中初始条件的函数带入计算即可。

关键词 无界弦 行波法 傅里叶变换法 拉普拉斯变换法 分离变量法 格林函数法

一、引言

物理上及工程技术上常需要研究各种各样的振动问题，如弦的振动，杆的振动、膜的振动、体的振动等。弦的振动又有无界弦[1]

的振动、有界弦的振动。其中，研究无界弦的振动问题受到了人们的重视。

通过众多学者的努力，对无界弦振动问题的研究方法越来越多

[2-6]

。比如在运用特征线方法的基

础上利用线积分予以求解[3]

；有学者用分离变量法求解[4]

，将分离变量形式的解代入泛定方程求出泛定方程的特解，再将所有可能的特解线性组合为通解，最后将初始条件代入通解计算各项系数，最后得出定解。分离变量法本来适用于有界问题，作者这里用它求解无界问题，开拓了求解无界弦振动问题的新思路。还有用傅里叶变换法[5]

、行波法[6]

等求解无界弦振动问题。本篇文章将用行波法、傅里叶变换法、拉普拉斯变换法、分离变量法、格林函数法求解无界弦的自由振动和受迫振动问题。通过比较，找出计算比较简便的方法和最佳方法，并且运用Matlab 软件模拟出无界弦自由振动的几个图形，方便大家理解弦的自由振动。

二、无界弦的振动问题

无界弦的振动问题包括无界弦的自由振动和受迫振动。两种问题的方程分别为(I )和(II )它们都由泛定方程[1]

和初始条件[1]

构成。无界弦自由振动的泛定方程为(I )中的(1)式，受迫振动的泛定方程为(II )中的(1)式，两者的初始条件为(2)式和(3)式。其中tt u 是弦的横向加速度；xx u 是u 关于x 的二阶导，质点间的牵连体现在xx u 上；a 是振动在弦上的传播速度，错误！未找到引用源。是t 时刻

作用于x 处单位质量上的横向外力，错误！未找到引用源。是初始位移，)0,(x u t 是初始速度，错误！未找到引用源。和错误！未找到引用源。是任意函数，由具体题目给定。

(I ) ⎪⎩⎪

⎨⎧===-)3()()0,()2()

()0,()

1(02x x u x x u u a u t xx tt ψϕ (II ) ⎪⎩⎪

⎨⎧===-)

3()()0,()2()

()0,()1(),(2x x u x x u t x f u a u t xx tt ψϕ

1、无界弦的自由振动问题

这是一种最简单的情况：一根无限长的均质柔软轻弦在初始条件作用下所引起的自由横向振动在弦中传播的情况。其定解问题为(I )。 （1）行波法

[6]

错误！未找到引用源。式，从而求其通解。②用定解条件确定通解中的任意函数（或常数），从而求得其特解。

⎪⎩⎪

⎨⎧===-)

3()()0,()2()()0,()1(02x x u x x u u a u t

xx tt ψϕ 泛定方程(1)的通解是

)4()

()(),(21at x f at x f t x u -++=

其中1f 和2f 是是两个任意函数，其形式可由初始条件确定。将(4)式代入(2)式和(3)，有：

)

()()()0,()()()()0,('2'121x x af x af x u x x f x f x u t ψϕ=-==+=

即 c d a

x f x f x

+=

-⎰ξξψ021)(1)()( 则⎰++=x c d a x x f 012)(21)(21)(ξξψϕ ⎰--=x c

d a x x f 022

)(21)(21)(ξξψϕ 则⎰++++=+at x c

d a at x at x f 012)(21)(21)(ξξψϕ

⎰-+-+=-at x c

d a at x at x f 022

)(21)(21)(ξξψϕ

所以，由(4)式得：

⎰+-+-++=at

x at

x d a at x at x t x u ξξψϕϕ)(21)]()([21),( （5）

由于大多数偏微分方程的通解难以求得，用定解条件定任意常数或函数也绝非易事，故行波法有很大的局限性，但对于研究波动问题，有其特殊优点。

-10-8-6-4-20246810

-0.5

-0.4

-0.3

-0.2

-0.1

0.1

-10

-8-6-4-20246810

-0.10

0.1

0.2

0.3

0.4

0.5

u

例：现取初始位移⎪⎩

⎪⎨⎧≤≤=)

(0)7473(

7sin )(其余l

x l x

l

x πϕ，初始速度0)(=x ψ，由达朗贝尔公式得

)](7sin )(7[sin 21),(at x l

at x l t x u -++=π

π，用matlab 作出图像如图1：

图1 初位移不为

0初速度为0的达朗贝尔公式的图形

现取初始位移0)(=x ϕ，初速度为⎩⎨⎧≤≤=)(0

)

10(1

)(其余x x ψ，用matlab 作图如下：

图2初位移为0，初速度不为0的达朗贝尔公式的图形

我们把图2分解为图3和图4，图3为开始时)(at x +ψ的波形，图4为开始时)(at x --ψ的波形。

图3 图4

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

-0.05

-0.04-0.03-0.02-0.0100.01

0.020.030.040.05-10

-8

-6

-4

-20246810

-1-0.8-0.6-0.4-0.200.20.40.60.81