参数设置
模式编号 1=(0,1)、2=(1,1)、3=(2,1)、4=(0,2)、5=(3,1)、6=(1,2)。缩略图列表中黄色框为当前模式。
膜的变位图案 u(r,θ)
膜的俯视图。蓝色=变位向下/红色=变位向上/白色=节线（|u|<ε）。同心节圆和放射状直径节线因模式而异。
6 个低次模式对比
排列 6 个低次模式 (0,1)～(1,2)。每个单元显示 (m,n)、α 值和相对振动频率 f/f_01。当前模式用黄色框高亮。点击可切换到该模式。
理论与主要公式
对于半径 $R$、张力 $T$ [N/m]、面密度 $\rho_s$ [kg/m²] 的圆形膜，波动方程通过变量分离求解，固有函数为
$$u_{mn}(r,\theta,t) = J_m\!\left(\alpha_{mn}\,\frac{r}{R}\right)\cos(m\theta)\,\cos(2\pi f_{mn} t)$$
其中 $\alpha_{mn}$ 是 Bessel 函数 $J_m(x)$ 的第 $n$ 个零点。固有振动频率为
$$f_{mn} = \frac{\alpha_{mn}}{2\pi R}\sqrt{\frac{T}{\rho_s}}$$
波数为 $k_{mn}=\alpha_{mn}/R$，周期为 $T_p = 1/f_{mn}$。低次模式的 $\alpha_{mn}$ 值为 (0,1)=2.4048、(1,1)=3.8317、(2,1)=5.1356、(0,2)=5.5201、(3,1)=6.3802、(1,2)=7.0156，互不成整数比，因此鼓具有非整数倍音的音色。
圆形膜振动模拟器简介
🙋
鼓的音色和吉他不一样，音高不清晰，这是为什么呢？
🎓
问得很好。原因是"倍音不是整数比"。吉他弦的频率是 f、2f、3f、4f...（整数倍），所以人耳会把基音 f 识别为音高。而圆形膜以 f_01 为基准，后续倍音是 1.594×、2.135×、2.296×、2.653×、2.917×...（非整数比）。因此听不出明确的"哆、唻、咪"音高，只是"砰、嗵"的打击乐声。
🙋
为什么 α_mn = 2.4048 这样的半整数会出现？
🎓
这是 Bessel 函数 $J_m(x)$ 的零点，是一个超越数，无法用初等公式求解。边界条件是"膜边缘处变位为零"，写成 $J_m(kR)=0$，其解 $kR = \alpha_{mn}$ 就被量子化。$J_0(x)=0$ 的第一个解是 x=2.4048...，第二个是 5.5201...。这和弦的 $\sin(kL)=0$ 得到 $kL=n\pi$ 是一个道理，只是圆形时用 Bessel 代替了 sin。
🙋
把张力改成 2 倍时，频率增加了 √2 倍左右，这和公式一致吗？
🎓
完全一致。从 $f_{mn} = (\alpha_{mn}/2\pi R)\sqrt{T/\rho_s}$ 可知，张力 T 增加 2 倍时，$\sqrt{T}$ 增加 √2 倍，频率就增加 1.414 倍。这和弦的 $f = (n/2L)\sqrt{T/\mu}$ 是同样的标度规律。定音鼓演奏者踩踏板改变 T，就靠这个 √T 关系连续地改变音高，所以音可以平滑地滑动。
🙋
白线有两种吧？直径方向的和圆形的，它们分别是什么？
🎓
都是"变位为零的节线"，但来源不同。直径方向的线来自周向因子 $\cos(m\theta)$ 的零点，m 越大就有越多条直线（2m 条）。圆形的节圆来自径向因子 $J_m(\alpha_{mn} r/R)$ 的零点，n 越大就有越多的同心圆（n−1 条在内侧）。如果往膜上撒砂粒，用低频加振，砂粒就会聚集在这些白线上，形成克拉德尼图形，是可视化振动模式的经典方法。
常见问题解答
弦的波动方程是一维的 Helmholtz 方程 $u_{tt}=c^2 u_{xx}$，边界条件 $\sin(kL)=0$ 给出 $k_n = n\pi/L$ 这样的整数比固有值列，产生"整数倍音"=清晰音高。而圆形膜是二维极坐标问题 $u_{tt} = c^2 \nabla^2 u$，径向方向出现 Bessel 微分方程。Bessel 函数 $J_m(x)$ 的零点不成整数比，而是 2.405、3.832、5.136、5.520... 这样的超越数序列，导致固有振动频率的比也是非整数。这是"鼓作为打击乐器、音高模糊"的根本物理原因。塔布拉鼓或定音鼓通过抑制某些模式的方法，能使某些频率比接近整数，因而能发出"有音高的打击乐"。
理想圆形膜确实是非整数倍音，但定音鼓通过几项调整实现了清晰音高：（1）击打位置通常在鼓头中心外 1/4 半径处，不是中心，能抑制基本轴对称模式 (0,1)；（2）鼓壳下方的空气共鸣；（3）鼓头内侧特殊的张力分布。这样就让 (1,1)、(2,1)、(3,1)、(4,1) 等模式 (m≥1) 主导，它们的 α 值为 3.832、5.136、6.380、7.588，相对比为 1.000、1.340、1.665、1.981。最后一个接近 2 倍音，形成"近似整数倍"，耳朵就能感知出清晰的音高。
本工具的 $u_{mn}(r,\theta) = J_m(\alpha_{mn} r/R)\cos(m\theta)$ 正是圆形膜波动方程的解析解，即固有函数。用 Ansys/Abaqus/Nastran 做圆形薄板固有值分析，得出的模态形状正是这个形式（板的情况加上弯曲刚性 D，变成 Kirchhoff 板）。触摸屏、扬声器振动膜、MEMS 传感器膜、电子乐器传感器等圆形薄膜结构的共振设计，都需要这种分析。本模拟器是一个教学性的例子，说明"手工计算也能重现 FEM 结果"。
$m≥1$ 的模式存在 $\cos(m\theta)$ 和 $\sin(m\theta)$ 两种形式的简并（频率相同），它们的任意线性组合也是特征函数。本工具只显示 $\cos(m\theta)$ 形式，但实际上膜的不均一性或打击位置会混合这两种形式，节线会显得旋转。椭圆形膜则会解除简并，分裂成两个不同频率。$m=0$ 的模式（轴对称）不简并，(0,1) 和 (0,2) 都是单独的振动频率。
实际应用
打击乐器设计（定音鼓、塔布拉鼓、康加鼓）：定音鼓通过脚踏机制连续改变张力 T，使 $\sqrt{T}$ 与频率成正比，可以在跨越约一个八度的音域内调整。塔布拉（印度鼓）在膜中心贴一块黑色圆形贴片（叫做 Shaya），局部增加质量，抑制特定模式，从而产生整数倍的谐波。康加鼓和杰姆贝利用木制壳体的空气共鸣，强化低音模式。本工具改变张力和面密度，能直观理解这些乐器的调音原理。
扬声器振动膜设计：动圈扬声器的圆锥纸或球顶高音喇叭是边缘固定的圆形膜。如果共振模式 (0,1)、(1,1)、(2,1) 等落在可闻频域内，音频响应曲线就会有山谷，导致音色变色。设计时需要改变锥纸的形状、密度分布、边缘的支持刚性，把这些模式推到可闻带外。本工具提供了初期估算的模型精度。
MEMS 压力传感器和麦克风：硅基 MEMS 压力传感器（电容式）用半径 100µm～1mm 的薄硅膜来感应压力变化。工作频率范围（音频传感器是 20Hz～20kHz）的上限由基模 (0,1) 决定。设计时通过调整膜半径、张力（由夹持方式确定）、厚度，让共振频率在 30kHz 以上。如果考虑膜厚，本公式可用于初步估算。
鼓的录音与扩声技术：麦克风放在鼓头上方不同位置时拾音不同。靠近中心会拾到低频的 (0,1) 模式，靠近边缘会拾到高次模式。从本工具的节圆和腹的分布可以看出，"麦克风位置导致的音色差异"本质上就是"拾音模式范围不同"。专业录音师利用这一点，bass drum 中心附近放麦克风获得低音，snare 的边缘放麦克风获得明亮音色。
常见误解与注意点
最常见的误解是以为"鼓的倍音也像弦一样是 f、2f、3f、4f 的整数比"。实际上圆形膜的固有频率由 α_mn 的比值决定，是非整数的 1、1.594、2.135、2.296、2.653、2.917...，这正是弦乐和打击乐音色根本不同的地方。本工具将 6 种模式并排显示，每个模式的相对频率 f/f_01 可直接读出。打击乐里唯一的"有音高的"就是定音鼓和塔布拉，它们通过物理手段刻意选出接近整数比的频率组合。
另一个常见误解是"α_mn 可以用某个简单公式计算"。实际上 Bessel 函数零点没有闭式解，只能用数值方法（牛顿法或级数截断）求得。本工具把低次的 6 个值（α_01～α_12）硬编码为常数，计算 Bessel 函数时用冪级数展开到 60 项。工程上通常查表（如 Abramowitz-Stegun）或调用数学库（Boost.Math）。
最后，不要以为"边缘固定松动时 α 会改变"。本工具假设边界条件"边缘完全夹紧"（位移为零），给出理想化的解析解。现实的鼓皮用钉子或绳子固定在木框上，夹持刚性有限，所以实际频率会比计算值低 10～20%。CAE/FEM 分析时可以将边界设为"弹性支持"，加入弹簧常数，得到更接近实测的值。教学意义上，本工具的值是"理想情况的上限值"，对比实测数据很有启发。
具体计算案例
分析一个直径 360mm、张力 3500 N/m、面密度 1.5 kg/m² 的音乐会低音鼓：基础模式 (m=0,n=1) 的 α_01≈2.405，固有频率 f₁=(2.405/2π)×√(3500/1.5)/0.18≈82.4Hz，与实测值 80～85Hz 吻合。第一倍频 (m=0,n=2) 的 α_02≈5.52，f₂≈190Hz，为鼓的余音形成。增加张力至 5000 N/m 时，f₁≈101Hz，音高更加突出。