🧑🎓
「1次元FEM固有値解析」って何ですか？棒が震える様子をコンピュータで予測するってことですか？
🎓
その通り！ざっくり言うと、細長い棒や梁を仮想的にバネとおもりの集合体に置き換えて、その「揺れやすいパターン」と「揺れる速さ」を計算する方法だよ。例えば、高層ビルが風でどの高さで大きく揺れるか、エンジンのクランクシャフトがどの回転数で激震するかを予測するのに使われるんだ。このシミュレーターでは、上のスライダーで「要素数」を増やすと、棒がより細かいバネに分割されて計算精度が上がる様子がわかるよ。
🧑🎓
え、バネとおもり？「一貫質量マトリクス」とか聞いたことありますが、それがその「おもり」に対応するんですか？
🎓
鋭いね！実務では「剛性マトリクス」がバネの硬さ、「質量マトリクス」がおもりの重さと慣性を表すんだ。このシミュレーターで「密度」や「断面積」を変えると、質量マトリクスが変わって固有振動数がどう変わるか、すぐに確認できる。軽くすると振動数が高くなる（速く揺れる）のが体感できるはずだよ。
🧑🎓
「境界条件」で結果が大きく変わるみたいですが、例えば片側だけ固定した片持ち梁と、両端を固定した場合では、どれくらい振動数が違うんですか？
🎓
良いところに気づいたね！固定の仕方で振動のしやすさがまるで変わるんだ。例えば同じ棒でも、片持ち（固定-自由）の最低次振動数は、両端固定の約1/4になる。つまり、片持ちの方がはるかにゆっくり揺れるんだ。シミュレーターで「境界条件」を切り替えて、グラフの振動モードの形と数値がどう変わるか、実際に試してみて！
有限要素法では、連続体である棒を有限個の「要素」に分割します。各要素は2つの節点で結ばれた単純な棒要素とみなし、その剛性（変形のしにくさ）を表すのが要素剛性マトリクスです。
$$\mathbf{K}_e = \frac{E A}{L_e}\begin{bmatrix}1 & -1\\-1 & 1\end{bmatrix}$$
ここで、$E$はヤング率（材質の硬さ）、$A$は断面積、$L_e$は要素長さです。この行列は、節点に力を加えた時の変形の関係（フックの法則）を行列形式で表しています。
要素の質量と慣性を考慮するために用いられるのが、一貫質量マトリクスです。要素の質量を節点に適切に分配したモデルです。
$$\mathbf{M}_e = \frac{\rho A L_e}{6}\begin{bmatrix}2 & 1\\1 & 2\end{bmatrix}$$
ここで、$\rho$は密度です。この行列により、要素の運動エネルギーを正確に表現できます。全ての要素の$\mathbf{K}_e$と$\mathbf{M}_e$を組み立てて全体の行列$\mathbf{K}$と$\mathbf{M}$を作り、一般化固有値問題$\mathbf{K}\boldsymbol{\phi}= \omega^2 \mathbf{M}\boldsymbol{\phi}$を解くことで、固有角振動数$\omega$と振動モード$\boldsymbol{\phi}$が求まります。
よくある誤解と注意点
このシミュレーターを使い始めるときに、特に気をつけてほしいポイントがいくつかあります。まず、「要素数を増やせば必ず精度が上がる」というのは誤解です。確かに、理論的には要素を細かくするほど真の解に近づきますが、コンピュータの計算には「丸め誤差」がつきものです。例えば、要素数を100から1000に増やしても、表示される固有振動数がほとんど変わらないどころか、場合によっては数値的な不安定性で結果がおかしくなることもあります。実務では「計算コストと精度のトレードオフ」を考え、必要な精度を得られる最小限の要素数を見極めることが重要です。
次に、パラメータ設定の単位系には一貫性を持たせましょう。このツールでは無次元化されているかもしれませんが、自分で計算する場合、例えばヤング率を[GPa]、密度を[kg/m³]、長さを[m]で入力したなら、断面積は[m²]で統一する必要があります。単位が混在すると、とんでもない数値が出てきて「計算が合わない」という事態になります。例えば、鋼材(E=210 GPa, ρ=7850 kg/m³)で長さ1mの片持ち梁を考えるとき、断面積を100 mm²と入力するなら、0.0001 m²に換算しなければなりません。
最後に、「計算された低次モードだけ見ていればいい」と思わないでください。確かに基本振動数は最も重要ですが、例えばある特定の周波数で強制振動が起きる場合、2次、3次の高次モードが共振して破損に至るケースもあります。このシミュレーターで境界条件を「固定-自由」にし、3次モードまで表示させてみてください。節（振動しない点）の位置がどのように変化するか観察すると、高次モードの複雑さが実感できます。
関連する工学分野
この1次元棒の固有値解析は、一見シンプルですが、実はさまざまな工学分野の基礎として顔を出します。まず挙げるのは音響工学です。楽器の弦や管の中の空気柱の振動は、支配方程式が全く同じ形（波動方程式）をしています。例えば、クラリネットの筒を1次元の管とみなし、一端閉塞（固定端）、一端開放（自由端）の境界条件で固有振動数を求めると、実際の音階（奇数倍音が卓越する）を説明できます。このシミュレーターで「固定-自由」条件に設定した時の振動モードの形が、まさに閉管の定在波そのものなのです。
もう一つの重要な応用先は伝熱工学の分野です。熱が物体中を伝わる「熱伝導」の非定常問題と、構造物の振動問題は、数学的に「同型」です。振動の固有値問題 $\mathbf{K}\boldsymbol{\phi}= \omega^2 \mathbf{M}\boldsymbol{\phi}$ は、熱伝導の時間定数を求める問題に置き換えることができます。つまり、このシミュレーターで計算しているアルゴリズムを理解すれば、エンジン部品がどの速さで暖まるか、電子基板の過熱を防ぐための熱拡散の設計などにも応用が利く考え方を身につけられるのです。
発展的な学習のために
このツールに慣れてきたら、次のステップとして「なぜその数式になるのか？」という背景を学ぶことをお勧めします。まず手に取って欲しいのは、変分法と重み付き残差法の基礎です。有限要素法は「ポテンシャルエネルギー最小化」や「微分方程式を近似的に満たす」というこれらの原理から出発しています。例えば、ここで使われている剛性マトリクスは、「仮想仕事の原理」から導出されているんです。
次に、2次元や3次元への拡張をイメージしましょう。1次元の棒要素が2次元では三角形や四角形の「シェル要素」や「ソリッド要素」に変わります。剛性マトリクスも $2\times2$ から、例えば三角形要素なら $6\times6$ の大きさになります。このシミュレーターで「要素分割」の概念を体感した後で、2次元メッシュの画像を見ると、その意味がぐっと理解しやすくなるはずです。
最後のステップは、減衰や非線形を考慮した現実的な振動解析です。実際の構造物には必ずエネルギーを散逸する減衰（粘性や構造減衰）が存在し、振動振幅が大きくなると材料の特性が変わる「幾何学的非線形」も無視できません。このツールで学んだ「無減衰・線形」の固有値解析は、それら全ての複雑な解析の土台となる、最も純粋で重要な第一歩なのです。