$$N_c = \frac{60}{2\pi}\sqrt{\frac{k}{m}}\text{ rpm}$$
不釣り合い応答振幅:
$$A = \frac{e \cdot r^2}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}$$
$r = \omega/\omega_c$(速度比)
質量・軸剛性・減衰比・偏心量を入力して臨界速度Ncを算出。不釣り合い応答曲線と歳差運動アニメーションでロータの共振挙動を直感的に理解しよう。
$$N_c = \frac{60}{2\pi}\sqrt{\frac{k}{m}}\text{ rpm}$$
不釣り合い応答振幅:
$$A = \frac{e \cdot r^2}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}$$
$r = \omega/\omega_c$(速度比)
タービン・発電機:火力や原子力プラントの大型タービンでは、複数の臨界速度が設計範囲内に存在します。起動・停止時にこれらの速度域を如何に安全に通過させるかが運転シミュレーションの重要課題です。
航空機エンジン:ジェットエンジンのコンプレッサーやタービンロータは非常に高速度で回転します。軽量化のために剛性を抑えると臨界速度が下がるため、ほぼ確実に超臨界運転領域に入り、その設計が必須となります。
産業用ポンプ・ファン:化学プラントなどで使われる多段遠心ポンプでは、内部の羽根車(インペラー)がロータを構成します。軸受の劣化による支え剛性の変化が、臨界速度を変え、振動問題を引き起こすことがあります。
自動車のターボチャージャー:エンジンの排気で駆動されるタービンは小型ながら数十万rpmで回転します。この超高回転域で安定して動作させるため、ロータダイナミクス解析により臨界速度を十分に高い位置に設計します。
まず、「臨界速度は絶対に通過してはいけない危険速度だ」という思い込み。確かにその速度で定常運転するのは危険ですが、多くの高速回転機械は「超臨界運転」を前提に設計されています。重要なのは、減衰を適切に設計し、起動・停止時に素早く通過することです。例えば、あるターボ分子ポンプでは、1次の臨界速度が10,000 rpm付近にありますが、運転速度は80,000 rpmです。停止時にこの10,000 rpm域で一瞬振動が大きくなりますが、減衰設計により安全に通過できます。
次に、シミュレーターで使う「等価質量」や「等価剛性」の設定。これは現場で一番悩むところです。例えば、軸の中央に集中した質量(m)を考える時、実際のロータの質量をそのまま代入するわけではありません。軸自体の質量の一部(通常1/3から1/2)を加算する必要があります。剛性(k)も、軸の剛性だけでなく軸受や支持構造の剛性が支配的になるケースが多々あります。ベアリングの種類(玉軸受 vs すべり軸受)やオイルフィルムの剛性が大きく影響するので注意しましょう。
最後に、「偏心量(e)をゼロにすれば振動はなくなる」という理想論。現実では加工・組立誤差、材質の不均一、運転中の熱変形や摩耗で必ず不釣り合いは発生します。ツールでeを変えて振幅がどう変わるかを見るのは、バランス調整の効果を予測するためです。例えば、偏心量を10μmから5μmに削減(バランス修正)すれば、共振点での振幅はほぼ半分にできると定性的に理解できます。完全なゼロを目指すよりも、許容振幅以内に収めるコスト感覚が重要です。
本ツールの既定値(質量m=10kg、軸剛性k=10×10⁵=1×10⁶N/m、偏心e=0.1mm、減衰比ζ=0.05)の場合を計算します。臨界速度Nc=(60/2π)√(k/m)=(60/2π)√(1×10⁶/10)≈3020rpm(ωc=316.2rad/s)と求まり、共振点での振幅ピークe/(2ζ)=0.1/(2×0.05)=1.0mmです。運転速度1000rpmでは速度比r=n/Nc≈0.331(亜臨界)で振幅は約0.012mm。運転速度を4000rpmへ上げると r≈1.325 の超臨界域に入り振幅は約0.229mmとなります。軸剛性kを下げると臨界速度は√k に比例して低下します。