参数设置
模态序号 1=(0,1)、2=(1,1)、3=(2,1)、4=(0,2)、5=(3,1)、6=(1,2)。缩略图上的黄色方框标记当前模态。
膜面位移 u(r,θ)
从上往下看的膜面图。蓝色＝位移向下／红色＝位移向上／白色＝节线（|u|<ε）。每个模态都有自己的节圆与直径节线。
6 个低阶模态比较
并排显示 (0,1)〜(1,2) 6 个低阶模态。每格标注 (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}$ 是 $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}$。低阶零点：(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×、1.594×、2.135×、2.296×、2.653×、2.917×… 都是无理数比。所以鼓声更像「砰」「咚」，音高感觉模糊，是典型的打击乐音色。
🙋
α_mn = 2.4048 这种奇怪的数字是从哪里来的？
🎓
它是贝塞尔函数 $J_m(x)$ 的零点，是没有解析表达的超越数。边界条件 u(R)=0 要求 $J_m(kR)=0$，所以 kR 必须等于这些零点 $\alpha_{mn}$。$J_0$ 的第一零点是 2.4048…，第二是 5.5201…，依此类推。和弦上 $\sin(kL)=0$ 产生 $kL=n\pi$ 的逻辑完全相同，只是圆形几何把 sin 换成了贝塞尔。
🎓
完全吻合。$f_{mn} = (\alpha_{mn}/2\pi R)\sqrt{T/\rho_s}$，T 翻倍则 $\sqrt{T}$ 提升约 1.414 倍。这就是鼓皮越拧紧音越高的物理依据，也和弦的 $f=(n/2L)\sqrt{T/\mu}$ 同型。定音鼓踏板就是连续改变 T，从而平滑地扫过音高范围。
🎓
两种都是「位移为零」的节线，但来自不同因子。直线节来自方位因子 $\cos(m\theta)$ 的零点，m 越大直线越多（共 2m 条腹相间排列）。圆节来自径向因子 $J_m(\alpha_{mn} r/R)$ 的零点，n 越大内部圆节越多。在膜上撒砂并以共振频率激励，颗粒会聚集到这些节线上，与克拉尼图形完全相同的原理。
常见问题
弦满足一维 Helmholtz 方程 $u_{tt}=c^2 u_{xx}$，边界 $\sin(kL)=0$ 给出 $k_n=n\pi/L$ 这一整数等差谱，从而产生「整数泛音＝清晰音高」。圆膜则是二维极坐标方程，径向变成贝塞尔微分方程，$J_m(x)$ 的零点是非整数比的超越数，所以圆膜频谱本质上是非整数谱。这就是为何鼓被定义为「无明确音高的打击乐」。定音鼓与塔布拉通过质量加载、空气耦合或鼓皮局部阻尼，从这些非整数泛音中挑出近似整数的子集，所以能演奏旋律。
定音鼓采用三大设计：在半径约 1/4 处敲击、鼓腔空气耦合、鼓皮中央质量微调。三者共同抑制轴对称的 (0,1) 模态，保留 m≥1 的 (1,1)、(2,1)、(3,1)、(4,1) 模态，对应 α = 3.832、5.136、6.380、7.588。以 (1,1) 为基准归一化，得到 1.000、1.340、1.665、1.981，最后一项几乎是 2 倍频，1.340 接近完全四度。物理上正是通过这些精心调节实现了「有音高的打击乐」。
本工具的 $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)$ 频率相同（简并），任何线性组合也是本征函数。本工具仅显示余弦项，但真实膜中两者的混合会因敲击位置和不均匀性而显得节线在旋转。当鼓皮被压成椭圆时简并被打破，分裂为两个相近频率，听感上会产生缓慢拍音。m=0（轴对称）模态不简并，所以 (0,1)、(0,2) 为单一频率。
实际应用
打击乐设计（定音鼓、塔布拉、康加鼓）：定音鼓通过脚踏板连续改变张力 T，按 $\sqrt{T}$ 比例扫过近一个八度。塔布拉（印度鼓）在中央粘黑色质量加载片（萨亚伊），抑制特定模态形成近似整数泛音，从而能演奏旋律。康加鼓与非洲鼓借助木腔空气共鸣，强化低频模态。本工具通过调节张力、密度、半径，即可重现这些乐器的调音原理。
扬声器振膜：动圈式扬声器锥盆与球顶高音振膜可视为圆形（或锥形）膜。可听频率范围内的共振模态会引起音色染色，设计上需将 (0,1)、(1,1)、(2,1) 共振频率推到工作频段之外，方法包括调整锥盆几何、密度分布、边缘悬挂刚度。本工具足以用于球顶高音的一次估算。
MEMS 压力传感器与麦克风：电容式 MEMS 压力传感器通过电容感测半径 100µm〜1mm 的硅薄膜在压力下的位移。工作带宽上限由 (0,1) 模态决定，设计上调节半径、张力、厚度，将其推至 30kHz 以上。在公式中加入板弯曲刚度修正，本工具即可用于 MEMS 设计的初步估算。
鼓声录音与扩声：麦克风放在鼓皮中央上方主要拾取 (0,1) 模态（最大位移）；靠近边缘则拾取较高的 m≥1 模态。本工具的节线图能从「拾取的模态不同」角度解释中心与边缘拾音的音色差异。专业录音师正是基于此为底鼓选择中央拾音，为军鼓选择边缘拾音。
常见误解与注意事项
最常见的误解是「鼓的泛音和弦一样为整数比 f、2f、3f、4f」。事实上圆膜频谱以 α_mn 比例展开，比例序列 1、1.594、2.135、2.296、2.653、2.917… 都是无理数。这正是弦乐与打击乐音色差异的物理本质，也是鼓声「音高模糊」的根源。本工具的 6 模态比较视图能直接读出相对频率 f/f_01。定音鼓和塔布拉之所以能演奏旋律，是因为通过物理设计从这些非整数泛音中挑出了近似整数的子集。
其次是认为「α_mn 有简单的闭式公式」。贝塞尔函数零点没有闭式解，必须用牛顿迭代或截断级数等数值方法。本工具硬编码 6 个低阶零点，并用 60 项幂级数计算 $J_m$ 形成膜面图。工程上一般查表（Abramowitz–Stegun）或使用 Boost.Math、SciPy 等库。
最后，「实际鼓边界完全固定」也是误区。本工具采用理想化 u(R)=0 边界条件，但真实鼓皮压在有限刚度的鼓圈上，会使固有频率略微下降。生产环境的 CAE/FEM 模型会在边界添加有限弹簧常数以更接近实际。教育上，本工具给出的是「理想化的上限值」，可作为实测值的参照。