🧑🎓
「ねじり振動」って、普通の横揺れの振動と何が違うんですか？
🎓
ざっくり言うと、回転方向の「ねじれ」による振動だね。例えば、エンジンのクランクシャフトが回転トルクの変動で、ギュルギュルとねじれる現象さ。シミュレーターで「慣性モーメント」のスライダーを動かすと、ディスクの大きさが変わって、このねじれやすさがどう変わるかすぐにわかるよ。
🧑🎓
え、じゃあこの振動が大きくなるとまずいんですか？
🎓
そう、共振するとシャフトが折れることもあるんだ。実務では、エンジンの常用回転域に固有振動数が入らないように設計するのが重要。このツールの「キャンベル線図」を見てみて。回転数を上げた時に、エンジンの点火次数（例えば4気筒なら2次）の線と固有振動数の線が交わるところが危険な共振点だ。
🧑🎓
「モード形状」ってグラフで、節点って言われるゼロの点がありますけど、あれは何を意味してるんですか？
🎓
いいところに気が付いたね。節点は「ねじれない点」。その点を境に左右でねじれる向きが逆になるんだ。例えば、1次モードの節点の位置は、シャフトの中で最も大きなねじり応力がかかる場所の目安になる。ツールでシャフトの「ねじり剛性」を変えると、この節点の位置がどう動くか試してみよう。
このシミュレーターの核心は、多自由度系のねじり振動を表す運動方程式を行列形式で立て、固有値問題として解くことです。振動しない状態からの変位（ねじれ角）を考えます。
$$ M \ddot{\boldsymbol{\theta}}+ K \boldsymbol{\theta}= \boldsymbol{0}$$
ここで、$\boldsymbol{\theta}$ は各慣性体のねじれ角を並べたベクトル、$M$は対角成分に慣性モーメント$I_i$が並ぶ慣性行列、$K$はシャフトの剛性$k_i$から決まる剛性行列です。ドットは時間微分を表します。
振動解 $\boldsymbol{\theta}= \boldsymbol{\Theta}e^{j\omega t}$ を仮定して代入すると、以下の固有値問題が得られます。これを解くことで、系が持つ固有振動数 $\omega_n$ と、その時の振動の形（モード形状）$\boldsymbol{\Theta}_n$ が求まります。
$$ [K - \omega^2 M] \{\Theta\}= \{0\}$$
$\omega$: 固有角振動数 [rad/s], $\boldsymbol{\Theta}$: モード形状ベクトル（各点の相対的なねじれ角）, $n$個の慣性体があれば、$n$個の固有値（振動数）と固有ベクトル（モード）が得られます（1つは剛体モード）。
よくある誤解と注意点
まず、「慣性モーメントが大きいほど振動しにくい」は半分正解で半分間違いです。確かに単純な1自由度系ではそうですが、多慣性系では全体のバランスが命。例えば、4慣性系で両端のディスクだけを極端に重くすると、中間の軽いディスクが激しく振動する「局部振動モード」が現れ、かえって問題になることがあります。ツールでI1とI4だけを大きくしてモード形状を見てみると、I2とI3が大きく揺れる様子が確認できますよ。
次に、剛性は「強いほど安全」とは限らないという落とし穴。シャフトを太くしてねじり剛性を上げると、固有振動数は確かに高くなります。しかし、エンジンの高回転域で使う場合、かえって危険な高次のモード（例えば4次モード）が常用回転域に降りてきてしまうことがあるんです。キャンベル線図で、高剛性にした時に高次モードの線が左にシフトするのを確認してみてください。
最後に、減衰（ダンピング）を無視した結果の解釈。このシミュレーターは保守系（減衰なし）の固有値解析なので、共振点の「危険度」まではわかりません。現実にはダンパーがついていたり、材料内部でエネルギーが散逸します。つまり、計算上は共振点があっても、減衰が十分なら振幅は抑えられ、実害がないケースも多い。このツールの役割は「潜在的な危険箇所を洗い出す」ことまでで、最終判断には減衰を考慮した過渡応答解析が必要だと覚えておきましょう。
関連する工学分野
このツールの計算ロジックは、構造物全体の振動解析（モーダル解析）と完全に同じです。ビルや橋梁の地震時の揺れ、飛行機の翼のフラッター現象も、質量行列と剛性行列を立てて固有値問題を解く点では一緒。つまり、このねじり振動の理解は、より複雑な構造振動へと繋がる重要な第一歩なんです。
また、制御工学の「状態空間モデル」とも深く関係しています。今回の運動方程式 $$ M \ddot{\boldsymbol{\theta}}+ K \boldsymbol{\theta}= \boldsymbol{0}$$ は、状態変数を $\boldsymbol{x} = [\boldsymbol{\theta}, \dot{\boldsymbol{\theta}}]^T$ と置き直せば、$\dot{\boldsymbol{x}} = A \boldsymbol{x}$ という標準形で書けます。この行列Aの固有値が振動数に対応。つまり、振動解析は制御系の安定性解析の親戚のようなもの。回転機械の振動を積極的に抑制する「アクティブ振動制御」は、この二つの分野が融合した技術です。
さらに応用として、音響工学（特に室内音響や騒音解析）が挙げられます。閉じた空間の空気の振動（音波）も、波動方程式を離散化すれば同様の行列固有値問題に帰着し、その空間の「固有モード（定在波）」と「固有周波数」を求めることができます。機械のねじり振動が騒音の原因となることも多く、両分野は切っても切れない関係にあるんです。
発展的な学習のために
まず次の一歩は、「連続体としてのシャフト」を考えてみること。このツールは「集中質量モデル」で、慣性はディスクに、弾性はシャフトの一点に集中しています。でも実際のシャフトは質量と剛性が連続的に分布していますよね。これを学ぶと、無限個の固有振動数があることや、ねじり波動方程式 $\frac{\partial^2 \theta}{\partial t^2} = c^2 \frac{\partial^2 \theta}{\partial x^2}$ にたどり着きます。集中質量モデルはその近似だと理解できると、視野が一気に広がります。
数学的には、行列の固有値問題の数値解法を少し深掘りするのがおすすめ。ツールが一瞬で答えを出す裏側では、QR法やランチョス法といったアルゴリズムが働いています。特に、対称行列の固有値問題は物理工学の至る所で現れる最重要課題。その性質（固有値は実数、固有ベクトルは直交する等）を学ぶと、計算結果への信頼感が全く違ってきます。
実務に近づくなら、「過渡応答解析」と「周波数応答解析」が次のトピック。今回の固有値解析は「自由振動」がテーマでしたが、実際の機械はエンジン点火や路面衝撃といった「外力」で振動します。過渡応答解析はその時間歴を、周波数応答関数（FRF）は外力の周波数に対する応答の大きさを明らかにします。これらを学べば、「この共振ピークはどのモードに起因するのか？」と、現象とモードを結びつけて考える、真の振動トラブルシューティングができるようになるでしょう。