極断面2次モーメント:$J = \dfrac{\pi d^4}{32}$
ねじり剛性:$k_t = \dfrac{GJ}{L}$
円板慣性モーメント:$I_d = \dfrac{mR^2}{2}$
固有振動数:$f_n = \dfrac{1}{2\pi}\sqrt{\dfrac{k_t}{I_d}}$
角変位:$\theta(t) = \theta_0\cos(2\pi f_n t)$
軸の材料・寸法と円板の質量・半径を変えて,ねじり振り子の固有振動数と角変位波形をリアルタイムで確認。クランクシャフトやプロペラ軸の振動設計に直結する基礎計算ツール。
極断面2次モーメント:$J = \dfrac{\pi d^4}{32}$
ねじり剛性:$k_t = \dfrac{GJ}{L}$
円板慣性モーメント:$I_d = \dfrac{mR^2}{2}$
固有振動数:$f_n = \dfrac{1}{2\pi}\sqrt{\dfrac{k_t}{I_d}}$
角変位:$\theta(t) = \theta_0\cos(2\pi f_n t)$
自動車のクランクシャフト設計:エンジン内部のクランクシャフトは、ピストンからの周期的な力によりねじり振動を起こします。この固有振動数がエンジンの回転数(特に高回転域)と一致しないように、軸径や形状を設計段階でシミュレーションします。
船舶の推進軸系:エンジンからプロペラに至る長い推進軸は、プロペラの水抵抗による変動トルクでねじり振動を受けます。共振による疲労破壊を防ぐため、軸の寸法や中間にダンパーを設置するなどの対策が講じられます。
産業機械のロボットアーム:高速で繰り返し動作するロボットアームの関節部には、サーボモーターと減速機を繋ぐ軸があります。この軸のねじり剛性は、位置決め精度や応答速度に直接影響するため、精密に設計されます。
風力発電タービンの主軸:風の変動によりブレードに不規則なトルクが加わり、主軸にねじり振動が生じます。巨大な構造物であるため、固有振動数を風況の主要な変動周波数から遠ざけることが耐久性向上の鍵となります。
このシミュレーターを使い始めるとき、いくつか勘違いしやすいポイントがあるよ。まず第一に、「軸の材質は軽い方が振動が速くなる」と思いがちだけど、それは間違い。せん断弾性係数 $G$ は「ねじりのばね定数」の材料側の要素で、値が大きいほど剛性が高く、振動は速くなる。例えば、アルミニウム(G≈26GPa)より鋼(G≈79GPa)の方が、同じ形状なら約3倍ねじり剛性が高く、振動数も約1.7倍になるんだ。軽さ(密度)は、ここでは軸の慣性モーメントを無視しているので影響しない点に注意してね。
第二に、「円板の質量」と「円板の半径」は独立したパラメータではないということ。実際の設計では、質量 $m$ を固定して半径 $R$ を変えると、厚さや材質も変わるのが普通だよね。このツールでは $I_d = \frac{1}{2}m R^2$ で計算するから、半径を2倍にすると慣性モーメントは4倍になり、振動は一気に遅くなる。質量だけを2倍にした場合と、結果が大きく異なるから確認してみて。
最後に、実務での落とし穴。このモデルは「一自由度」の理想的なねじり振り子だから、軸の質量は無視し、円板は剛体と仮定している。でも実際の長いシャフトでは、軸自体にも分布質量があって「連続体」として振動する。すると、基本振動数だけでなく、2次、3次…という高次のねじり振動モードも存在するんだ。ツールで計算した振動数は、あくまで非常に剛性の高い(短い・太い)軸に対する一次近似値だと理解しておこう。
クランクシャフト設計例:鋼製軸(G=80GPa)、軸径d=30mm、長さL=150mm、フライホイール質量m=2.5kg、半径r=80mm。極断面2次モーメントJ=π×30⁴/32≈39,761mm⁴=3.98×10⁻⁵m⁴、ねじり剛性kt=80×10⁹×3.98×10⁻⁵/0.15≈2.12×10⁷N·m/rad、回転慣性Id=2.5×0.08²/2≈0.008kg·m²、固有振動数fn=√(2.12×10⁷/0.008)/(2π)≈65.2Hz、周期T≈0.0153秒。低周波外部振動との共振を避けるため、設計段階で50Hz以上の確保が重要。