🧑🎓
「液体スロッシング」って何ですか？タンクの中の水が揺れるだけなら、そんなに問題になるんですか？
🎓
ざっくり言うと、タンク内の液体が大きく揺れ動く現象だね。これが厄介なのは、地震や車の急ブレーキなどの振動と「共振」して、揺れがどんどん大きくなってしまうからなんだ。例えば、石油タンクでこれが起こると、タンクの屋根が破損したり、液体が溢れ出したりする大事故につながるよ。上のシミュレーターで「タンク幅」を変えてみると、揺れの周期（固有振動数）がどう変わるか、すぐに体感できるぞ。
🧑🎓
え、そうなんですか！共振を避けるために、揺れやすい周期を事前に知りたいわけですね。でも、計算は難しそう…。
🎓
そこで役立つのが「ノルトン則」に基づくこの計算式だ。実は、タンクの幅と液深ささえ決まれば、最初の5つの揺れ方（モード）の固有振動数がバーグラフでパッと出るんだ。シミュレーターで「液深」のスライダーを動かしてみて。液深が深くなるほど、振動数が上限に近づいていくのがわかるだろう？これが実務でよく使う「深い液体」の近似だ。
🧑🎓
「等価振り子長」って表示されてますけど、これは何に使うんですか？振り子と液体の揺れって関係あるんですか？
🎓
良いところに気が付いたね！一番基本的な揺れ方（第1モード）は、ある長さの振り子の動きで近似できるんだ。この「等価振り子長」がわかると、液体の揺れによる力を振り子の重りが壁を押す力として簡単に計算できる。例えば、燃料タンクの耐震設計で「スロッシング荷重」を見積もる時、この考え方がとても役立つよ。「液体密度」を変えて、等価振り子長がどう変わるかも確認してみよう。
矩形タンク内の液体の微小振幅振動を仮定した時、n次モードの固有角振動数は以下の分散関係式で与えられます。
$$ \omega_n^2 = g k_n \tanh(k_n h) $$
ここで、$f_n = \omega_n / 2\pi$ が固有振動数[Hz]です。$k_n = n\pi / a$ は波数[$m^{-1}$]、$a$はタンク幅[m]、$h$は液深[m]、$g$は重力加速度[$m/s^2$]、$n$はモード次数（1,2,3...）です。$\tanh$関数は液深の影響を表し、$h$が大きい（深い）と1に近づきます。
最も重要な基本モード（n=1）の動きを理解するために、等価振り子モデルが用いられます。この振り子の長さは以下の式で定義されます。
$$ l_{eq}= \frac{h \tanh(k_1 h)}{k_1 h}$$
$l_{eq}$は等価振り子長[m]です。この長さの振り子の周期が、液体スロッシングの第1モードの周期と一致します。液深$h$が浅い場合、$l_{eq}\approx h$（液深そのもの）に、深い場合、$l_{eq} \approx 1/k_1 = a/\pi$（タンク幅に依存）に収束するのが特徴です。
よくある誤解と注意点
このツールを使い始める際、特にCAE初心者が陥りがちなポイントがいくつかあるよ。まず大きな誤解は、「計算結果がそのまま実設計値」と思ってしまうこと。このシミュレーターは「微小振幅」を仮定した理想的な矩形タンクの理論解だ。実際のタンクは円筒形が多く、液面の揺れが大きくなると非線形現象が無視できなくなる。例えば、計算で出た固有振動数が1.0Hzだったとしても、実機では0.8Hz〜1.2Hzくらいのバラつきを見込んでおく必要があるんだ。
次に、パラメータ入力の落とし穴。液深「h」は静止時の深さだよね？でも、タンクが加速している最中は液面が傾くから、実効的な深さが変わる。車の燃料タンクで考えてみよう。急ブレーキで液面が前方に寄れば、後方の壁面にかかる圧力は計算値より小さくなる。逆に、旋回時は片側の壁面に想定以上の荷重がかかる可能性がある。シミュレーション結果は「一つの基準状態」と捉え、最も過酷な液面姿勢を想定した複数ケースの検討が実務では不可欠だ。
最後に、「第1モードだけ見ていればいい」という思い込み。確かに基本モードが最もエネルギーは大きいけど、高次モードも無視できない場面がある。例えば、タンク内部の細かい構造物（計測機器の取り付け支柱など）は、局所的に第3モードや第5モードの波の節（動かない点）や腹（大きく動く点）に位置するかもしれない。ツールで5モードまで可視化している意味を理解して、自らの設計対象がどのモードの影響を受けるか、考えるクセをつけよう。
関連する工学分野
液体スロッシングの計算式は、一見特殊だけど、実は波動現象を扱う様々な工学分野と根っこが同じなんだ。まず真っ先に挙がるのは「船舶工学」。船体という動く容器中の水の揺れ（バルキング）は、スロッシングと数学的に非常に近い問題だ。さらに、タンク内の表面波の分散関係式 $$ \omega_n^2 = g k_n \tanh(k_n h) $$ は、「海岸工学」で使う沖波から岸辺に近づく波の変形を記述する式と本質的に同一。液深「h」が浅くなるほど波速が遅くなる性質（tanhの挙動）は共通しているよ。
もう一つの大きな繋がりは「構造力学」、特に流体構造連成(FSI)解析だ。スロッシングで発生する動水圧はタンク壁を変形させ、その変形がさらに液体の揺れ方を変える…という相互作用が起こる。このツールで求まる固有振動数は「壁が剛体」という仮定での値だが、FSIを考慮すると振動数が少し低くなることが多い。また、等価振り子長の概念は、「耐震工学」で建物の屋上タンクなどを「振り子型ダンパー」としてモデル化する考え方にも通じる。液体の揺れを単純化した質量モデルとして扱うことで、全体の地震応答解析が格段に楽になるんだ。
発展的な学習のために
このツールの背後にある理論に興味が湧いたら、次のステップに進んでみよう。まず数学的背景を深めるなら、支配方程式である「ラプラス方程式」と境界条件（自由表面条件、底面条件、側面条件）の立て方を学ぶことから始めるのがオススメ。これにより、なぜ解が三角関数と双曲線関数の組み合わせ $$ \phi(x,z) \propto \cos(k_n x) \frac{\cosh(k_n z)}{\cosh(k_n h)} $$ になるのか、イメージが掴めるはずだ。
実務的な次の一歩は、「円筒形タンク」のスロッシングを調べることだ。矩形と違い、解にはベッセル関数が現れる。多くの設計基準（例えばAPI 650）は円筒タンクを対象としているので、ここを押さえると一気に実務に近づく。学習ステップとしては、1) このツールで矩形の感覚を掴む、2) 円筒では半径が幅の役割をし、モードは軸対称と非対称に分かれることを理解する、3) 防波板（バッフル）の効果を定性的に考える（波動を遮断・反射する）、という流れがいいね。
最終的には、非線形スロッシングや減衰の世界に挑戦してみよう。液面が大きく揺れると波打ち際でエネルギーが散逸するし、粘性の影響も無視できなくなる。これらの現象は、このツールの基礎にある線形理論を超える部分で、現在も研究が活発な領域だ。まずは、このシンプルなツールで「幅」「深さ」「密度」が振動数に与える影響を骨の髄まで体感しておくことが、それら複雑な現象を理解する最強の土台になるんだ。