固定边矩形膜振动的膜态分析.pdfVIP

固定边矩形膜振动的膜态分析 Page 1 of 8 固定边矩形膜振动的膜态分析 戴玲燕 摘要：膜振动有二维的，三维的，膜边界形状有规则和不规则的，本文是关于对固定边矩形膜振动 的研究。首先写出固定边矩形膜的定解问题，运用数学物理方程中所涉及的方法，例如分离变量法，拉 普拉斯变换法对其进行求解，此处我们用分离变量法求解，得出任一时刻矩形框内任一处膜振动的位移， 由于运用数学物理方程中所涉及的方法，得出精确解的形式不是积分就是级数，我们很难从解的表达式 中直观理解矩形膜振动的特性，所以可以借助图像的直观生动这些特性从而进一步了解矩形膜振动的特 性，这就需要运用 MATLAB 对解析解模拟，可视化。最后对矩形膜振动的振动过程进行讨论。 关键词：膜振动，分离变量，MATLAB，可视化 引言现实生活中膜振动的例子很多，不少实际生产中经常用到膜振动，如电话的话筒和听筒就用到 了膜振动的原理。数学物理方程中，描述膜振动过程的数学模型是膜振动方程。膜振动有二维的，三维 的，膜边界形状有规则和不规则的，定结问题的解与初始条件和边界条件有密切关系。 到目前为止，对于矩形膜的描述已有大量的资料参考，很多资料都导出了膜振动的方程。如哈尔滨 工业大学数学系组编的《数学物理方程》，复旦大学数学系组编的《数学物理方程》等等。对于矩形膜 振动的膜态分析，彭芳麟编写的《数学物理方程的 MATLAB 解法与可视化》用分离变量法的方法做了一 定分析。 本文中我们用分离变量法求解矩形区域内二维膜振动方程，然后用 MATLAB 的方法对其可视化，并 对物理过程进行讨论。 一 用分离变量法求解膜振动的混合问题 （一）长为 c 宽为 b 的矩形膜，其矩形框固定不变，若给定矩形框内任意处特定位移，在不受外界 情况影响下，矩形膜的振动情况，定解问题为   2 utt  a (uxx  uyy ) (0  x , 0  y  b, 0  t)(1)  u |  u |  u |  u |  0(2) x 0 x c y 0 y b    y u |  Ax (x  c) sin ,u |  0(3) t 0 t t 0  b 2 式（1）是膜振动的方程，完整式应是utt  a (uxx  uyy ) +f。其中u 是膜振动过程中，任意 处的横向位移，是x , y ,t 的函数，a 是任意实数，f 也是x , y ,t 的函数。这里 f 取 0。式（2）是矩形的 边界条件，在长为 c 宽为 b 的矩形框的四周，函数值为 0。式（3）是初始条件，在时间t=0 时刻任意  y 处的横向位移 u 的值取为Ax (x  c) sin ，横向位移u 对时间 t 的一阶偏导在时刻 t=0 的函数值取为 b 0。 （二）解定解问题的分离变量法 ①设（1）式有可以分离变量的特解 u(x , y ,t)  T (t)v(x)(4) 将上式带入（1）式，得到 2 2 2  T 2  v  v v(x , y ) 2  a T (t)( 2  2 ) t x y 1 固定边矩形膜振动的膜态分析 Page 2 of 8 2 2 2