🧑🎓
「臨界速度」って何ですか？タービンとかのカタログで見たことがあります。
🎓
ざっくり言うと、回転体が一番大きく振れ回ってしまう危険な回転速度だよ。例えば、洗濯機が脱水の途中で「ガタガタ！」と揺れるあの現象に近いね。このシミュレーターで「質量」と「剛性」のスライダーを動かすと、計算される臨界速度がどう変わるか、すぐに確認できるよ。
🧑🎓
え、じゃあその速度で運転しちゃダメなんですか？でもジェットエンジンのタービンはもっと速く回ってますよね。
🎓
良いところに気が付いたね。実務では「超臨界運転」といって、危険な速度域を素早く通り抜けて、その上で運転するんだ。シミュレーターで「運転速度」を臨界速度より高く設定して、グラフの山を越えてみて。減衰比を大きくすると、この山が低くなって通りやすくなるのがわかるよ。
🧑🎓
「偏心量」って何ですか？これもスライダーで変えられますが。
🎓
回転体の重心がほんの少し軸の中心からズレている量だ。完全なバランスは現実には不可能なんだ。この「e」の値を大きくすると、グラフの振幅全体が大きくなるでしょう？これが不釣り合い力の大きさで、現場で振動トラブルの原因で一番多いんだよ。
このシミュレーターの基礎は、1919年にJeffcottが提案した古典的で強力なモデルです。軸の中央に質量が集中し、バネとダッシュポットで支えられていると単純化して考えます。これにより導かれる臨界速度 $N_c$ は次の式で計算されます。
$$N_c = \frac{60}{2\pi}\sqrt{\frac{k}{m}}\quad \text{[rpm]}$$
$m$: ロータの等価質量 [kg]
$k$: 軸系の等価剛性 [N/m]
質量が大きく剛性が低い（柔らかい）ロータほど、臨界速度は低くなります。
運転速度が変化した時の振れ回り振幅 $A$ は、以下の「周波数応答関数」で表されます。グラフの曲線はこの式そのものです。
$$A = \frac{e \cdot r^2}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}$$
$e$: 重心の偏心量 [m]
$r$: 速度比 ($r = \omega / \omega_n = n / N_c$)
$\zeta$: 減衰比
$r=1$（運転速度が臨界速度と一致）の時に分母が最小となり、振幅が極大（共振）します。減衰比 $\zeta$ がこのピークの高さを決める鍵です。
よくある誤解と注意点
まず、「臨界速度は絶対に通過してはいけない危険速度だ」という思い込み。確かにその速度で定常運転するのは危険ですが、多くの高速回転機械は「超臨界運転」を前提に設計されています。重要なのは、減衰を適切に設計し、起動・停止時に素早く通過することです。例えば、あるターボ分子ポンプでは、1次の臨界速度が10,000 rpm付近にありますが、運転速度は80,000 rpmです。停止時にこの10,000 rpm域で一瞬振動が大きくなりますが、減衰設計により安全に通過できます。
次に、シミュレーターで使う「等価質量」や「等価剛性」の設定。これは現場で一番悩むところです。例えば、軸の中央に集中した質量（m）を考える時、実際のロータの質量をそのまま代入するわけではありません。軸自体の質量の一部（通常1/3から1/2）を加算する必要があります。剛性（k）も、軸の剛性だけでなく軸受や支持構造の剛性が支配的になるケースが多々あります。ベアリングの種類（玉軸受 vs すべり軸受）やオイルフィルムの剛性が大きく影響するので注意しましょう。
最後に、「偏心量(e)をゼロにすれば振動はなくなる」という理想論。現実では加工・組立誤差、材質の不均一、運転中の熱変形や摩耗で必ず不釣り合いは発生します。ツールでeを変えて振幅がどう変わるかを見るのは、バランス調整の効果を予測するためです。例えば、偏心量を10μmから5μmに削減（バランス修正）すれば、共振点での振幅はほぼ半分にできると定性的に理解できます。完全なゼロを目指すよりも、許容振幅以内に収めるコスト感覚が重要です。
関連する工学分野
この臨界速度計算の考え方は、見方を変えると様々な工学分野の「共振問題」と本質的に同じです。まず真っ先に挙がるのは構造動力学です。ビルや橋梁が風や地震で受ける振動も、固有振動数（臨界速度に相当）と強制振動力（不釣り合い力に相当）の関係で解析します。モード解析という手法は、ロータの複雑な振動モードを解くのにも使われます。
もう一つは制御工学、特に「振動制御」の分野です。シミュレーターで減衰比(ζ)を大きくすると共振ピークが抑えられるのを見ましたね。これは、能動的な制御（アクティブ振動制御）で言う「フィードバックゲインを上げて安定化させる」操作に似ています。実際、磁気軸受を使ったロータでは、この制御系の設計がそのまま減衰特性を決め、臨界速度通過のしやすさを左右します。
さらに材料力学・疲労解析とも深く関連します。共振時の大きな振動振幅は、軸に繰り返し応力を発生させ、材料の疲労寿命を大幅に短縮します。例えば、臨界速度付近で長時間運転されたポンプの軸が、設計寿命より早く破断したという事例は珍しくありません。振動解析から得られた応力振幅とサイクル数は、疲労寿命予測の直接的な入力データになります。
発展的な学習のために
Jeffcottロータのモデルに慣れたら、次のステップは多自由度系への拡張です。実機のロータは複数のディスクが付いた「連続体」に近い。これを複数の質量とバネでモデル化し、行列を使って運動方程式を立てます（$[M]\{\ddot{x}\} + [C]\{\dot{x}\} + [K]\{x\} = \{F\}$）。これを解くことで、第1次、第2次、第3次…と複数の臨界速度とそのときの振動モード（たわみ形）が求められます。これが「キャンベル線図」作成の基礎になります。
数学的には、今扱った周波数応答関数の導出過程を追うことが大切です。これは2階線形常微分方程式の強制振動解です。減衰のない場合の解から始め、次に減衰項を加え、さらに不釣り合い強制力の項を加える…という流れを教科書で確認してみてください。この理解があれば、ツールが出力するグラフの形（なぜr=1でピークが立ち、減衰でピークが丸くなるか）が腑に落ちます。
実務に直結する次のトピックは、「不安定振動」のメカニズムです。臨界速度のような共振は「不安定」とは呼びません。回転体特有の、ある速度域に達すると発散的に振動が大きくなる「オイルホイップ」や「ハイスピード不安定」といった現象があります。これらは、減衰が振動エネルギーを「吸収」から「供給」に変えてしまう自己励起現象で、超臨界運転設計における最大の関門の一つです。まずは、安定性を判別する「レイリーの安定性判別」の概念から学ぶことをお勧めします。