パラメータ設定
モード番号 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$ が出るのと同じ仕組みで、円の場合に sin の代わりに Bessel が出てくる、と覚えればいい。
🙋
スライダーで張力を 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}$ と同じスケーリングで、太鼓のチューニング（皮の張り具合）でピッチが変わる物理的根拠になっている。ティンパニ奏者がペダルを踏むと R は変わらず T だけ動かすから √T で連続的にピッチが滑らかに変わるんだ。
🙋
直径方向の「白い線」と同心円の「白い線」、両方ありますね。何が違うんですか？
🎓
どちらも「変位がゼロになる節線」だけど起源が違う。直径方向の線は周方向因子 $\cos(m\theta)$ の零点で、$m$ が大きいほど 2m 本の直線になる。同心円の線は半径方向因子 $J_m(\alpha_{mn} r/R)$ の零点で、$n$ が大きいほど内側に節円が増える。砂粒を膜の上に乗せて低周波で加振すると、ちょうどこの白い線の上に粒子が集まる Chladni 図形と同じ現象が観察できるよ。
よくある質問
弦の波動方程式は 1 次元の Helmholtz 方程式 $u_{tt}=c^2 u_{xx}$ で、境界条件 $\sin(kL)=0$ から $k_n = n\pi/L$ と整数比の固有値列が出ます。これが「整数倍音」=明瞭なピッチをもたらします。一方で円形膜は 2 次元極座標で $u_{tt} = c^2 \nabla^2 u$ となり、半径方向に Bessel 微分方程式が現れます。Bessel 関数 $J_m(x)$ の零点は整数比ではなく超越数の数列なので、固有振動数の比も非整数となります。これが「太鼓は弦と違ってピッチが曖昧な打楽器」である根本的な物理的理由です。タブラやティンパニのように特定モードを抑制する工夫で「歌う打楽器」を作ることもできます。
理想的な円形膜は非整数倍音ですが、ティンパニは「ヘッドの中心から半径の 1/4 付近を叩く」「シェル下の空気共鳴」「中心付近の質量分布の微調整」によって基本軸対称モード (0,1) を抑制し、(1,1)・(2,1)・(3,1)・(4,1) という対称性 m≥1 のモードのみを残します。これらの α は 3.832, 5.136, 6.380, 7.588 で、(1,1) を基準とすると 1.000, 1.340, 1.665, 1.981 となり、最後の 1.981 はほぼ 2 倍音、5.136/3.832=1.340 は完全 4 度に近い音程比になります。物理的にこれらの調整で「明確に音程の聞こえる打楽器」を実現しているのです。
本ツールの $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)$ の 2 つで縮退（同じ振動数）しており、両者の任意の線形結合も固有関数になります。本シミュレーターでは $\cos(m\theta)$ のみを表示していますが、実物では膜の不均一性や打撃位置によって両者が混合し、節線が回転して見えることがあります。これは縮退モードの「方向選択」と呼ばれる現象で、楕円形にゆがんだ太鼓では縮退が解けて 2 つの異なる振動数に分裂します。$m=0$ のモード（軸対称モード）は縮退しないため (0,1)・(0,2) は単独の振動数となります。
実世界での応用
打楽器設計（ティンパニ・タブラ・コンガ）：ティンパニはペダル機構で張力 T を連続的に変化させ、$\sqrt{T}$ に比例して基本振動数を 1 オクターブ近く動かします。タブラ（インドの太鼓）では中央に黒い円形パッチ（シャイ）を貼って質量を局所的に増やし、特定モードを抑制して整数倍音を作り出します。コンガやジェンベは木製シェルの空気共鳴と組み合わせて低音モードを強調しています。本ツールで張力・面密度を変えると、これらの楽器のチューニング原理が直感的に理解できます。
スピーカーの振動板設計：ダイナミックスピーカーのコーン紙やドーム型ツイーターは、エッジで固定された円形（または円錐形）膜の振動として扱えます。共振モードが可聴域内に入ると音色が変色するため、設計では (0,1)・(1,1)・(2,1) 等の固有振動数を狙った範囲外に追い出します。コーンの形状・密度分布・縁の支持剛性を調整して、基本モード以外の振動が音圧周波数応答に山谷を作らないようにします。本ツールは初期検討の解析モデルとして十分な精度をもちます。
MEMS 圧力センサー・マイク：シリコン MEMS の圧力センサー（コンデンサ型）は、半径 100µm〜1mm の薄い円形ダイヤフラムが圧力で変形する量を電気容量で読み取ります。動作周波数帯（音響用は 20Hz〜20kHz）の上限を決めるのが基本振動モード (0,1) で、ダイヤフラム半径とテンション・厚さを設計パラメータとして共振周波数を 30kHz 以上に押し上げます。本ツールの式に膜厚を加えれば MEMS 設計の一次見積もりに使えます。
太鼓の音響工学（録音・PA）：マイクロホンを太鼓のヘッド上方に置くとき、ヘッド中央（最大変位）に近づけると低音 (0,1) モードが、ヘッド端に寄せると高次モードが拾われます。本ツールの節円・腹のパターンを見ると、マイク位置による音色の違いが「拾えるモードの違い」として説明できます。プロのエンジニアはこれを利用してキックドラムでは中央付近、スネアではエッジ寄り、というマイキングを使い分けています。
よくある誤解と注意点
最も多い誤解は、「弦と同じように、太鼓の倍音も f, 2f, 3f, 4f と整数比で並ぶ」と考えることです。実際には円形膜の固有振動数は α_mn の比で決まる非整数比（1, 1.594, 2.135, 2.296, 2.653, 2.917, ...）になります。これが弦楽器と打楽器の音色を分ける本質的な物理であり、太鼓に「明確なピッチ」を感じにくい理由です。本ツールで 6 モードを並べて表示すると、各モードの相対振動数 f/f_01 が直接読み取れます。「打楽器なのに音程が聞こえる」ティンパニやタブラは、これらの非整数倍音から整数比に近い組を意図的に選び出している特殊な楽器です。
次に多いのが、「α_mn は何かの単純な公式で計算できる」と思い込むことです。実際には Bessel 関数の零点は閉形式の解析解を持たず、数値計算（Newton 法や級数の打ち切り）で求めるしかありません。本ツールでは低次の 6 個を定数として埋め込み、波動方程式の計算には Bessel 関数の冪級数展開を 60 項まで使っています。エンジニアリングではテーブル（Abramowitz-Stegun 等）や Boost.Math 等のライブラリを参照するのが標準的です。
最後に、「縁の固定がゆるいと α が変わる」と思わないことも重要です。本ツールは「縁で変位ゼロ（クランプ）」という理想境界条件の解析解で、実物の太鼓は皮を木枠に張り付ける部分の剛性が有限なため、僅かに振動数が下がります。CAE/FEM 解析では境界条件を「弾性支持」としてバネ定数を加えるとより現実的な値が得られます。教育的には本ツールの値が「理想化された上限値」として、実測値との比較に役立ちます。