ねじり振り子シミュレーター 戻る
振動解析ツール

ねじり振り子シミュレーター

軸の材料・寸法と円板の質量・半径を変えて,ねじり振り子の固有振動数と角変位波形をリアルタイムで確認。クランクシャフトやプロペラ軸の振動設計に直結する基礎計算ツール。

パラメータ設定
せん断弾性率 G
GPa
軸長さ L
m
軸直径 d
mm
円板質量 m
kg
円板半径 R
m
初期振幅 θ₀
°
計算結果
計算結果
J (m⁴)
kt (N·m/rad)
Id (kg·m²)
fn (Hz)
周期 T (s)
ディスク回転アニメーション
角変位波形 θ(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)$

ねじり振り子とは

🙋
「ねじり振り子」って、普通の振り子と何が違うんですか?
🎓
大まかに言うと、振り子が「揺れる」のに対して、ねじり振り子は「ひねる」振動だね。例えば、天井からワイヤーで円盤を吊るして、円盤をひねって手を離すと、元の位置に戻ろうとするねじれの力で振動するんだ。このシミュレーターでは、軸の長さや太さ、円盤の大きさを変えて、その振動の速さ(固有振動数)がどう変わるか試せるよ。
🙋
へー!じゃあ、軸の直径「d」を大きくすると、振動は速くなるんですか?遅くなるんですか?
🎓
いいところに気づいたね。軸が太くなると、ねじれにくくなる、つまり「ねじり剛性」が高くなるんだ。ばねが硬くなるイメージだね。だから、振動は速くなるよ。逆に軸を長く「L」を大きくすると、ねじれやすくなって振動は遅くなる。実際に上のスライダーで「d」を2倍にしてみると、計算される固有振動数が一気に跳ね上がるのがわかるはずだよ。
🙋
なるほど!で、この計算って、実務ではどこで使われるんですか?
🎓
自動車のエンジンや船舶のプロペラ軸の設計で特に重要だよ。エンジンのクランクシャフトは複雑なねじり振動を起こすんだけど、その固有振動数がエンジンの回転数と一致すると「共振」して、軸が折れる大事故につながる。このツールで学ぶ基礎が、そのような危険な共振を避ける設計の第一歩なんだ。円板の半径「R」を変えて慣性モーメントの影響を確認してみて。

よくある質問

固有振動数は軸の直径の2乗に比例して増加します。ねじり剛性k_tが直径の4乗に比例するため、直径を2倍にすると剛性は16倍になり、固有振動数は約4倍になります。設計時はこの感度の高さに注意してください。
円板の質量を大きくすると慣性モーメントが増加し、同じ外力に対する角変位の振幅は小さくなります。ただし、固有振動数は低下するため、共振条件が変わる点に注意が必要です。実際の設計では振幅と固有振動数のトレードオフを確認してください。
本ツールで軸のねじり剛性と円板の慣性モーメントから求めた固有振動数が、エンジンの回転数域での共振周波数と一致しないかを確認できます。共振を避けるため、軸径や円板質量を調整する際の基礎データとして活用してください。
一般的な鋼材ではG≈80GPa、アルミニウム合金ではG≈26GPa、真鍮ではG≈37GPaです。正確な値は材料メーカーの技術データシートやJIS規格表を参照してください。本ツールでは代表的な材料の標準値をプリセットとして用意しています。

実世界での応用

自動車のクランクシャフト設計:エンジン内部のクランクシャフトは、ピストンからの周期的な力によりねじり振動を起こします。この固有振動数がエンジンの回転数(特に高回転域)と一致しないように、軸径や形状を設計段階でシミュレーションします。

船舶の推進軸系:エンジンからプロペラに至る長い推進軸は、プロペラの水抵抗による変動トルクでねじり振動を受けます。共振による疲労破壊を防ぐため、軸の寸法や中間にダンパーを設置するなどの対策が講じられます。

産業機械のロボットアーム:高速で繰り返し動作するロボットアームの関節部には、サーボモーターと減速機を繋ぐ軸があります。この軸のねじり剛性は、位置決め精度や応答速度に直接影響するため、精密に設計されます。

風力発電タービンの主軸:風の変動によりブレードに不規則なトルクが加わり、主軸にねじり振動が生じます。巨大な構造物であるため、固有振動数を風況の主要な変動周波数から遠ざけることが耐久性向上の鍵となります。

よくある誤解と注意点

このシミュレーターを使い始めるとき、いくつか勘違いしやすいポイントがあるよ。まず第一に、「軸の材質は軽い方が振動が速くなる」と思いがちだけど、それは間違い。せん断弾性係数 $G$ は「ねじりのばね定数」の材料側の要素で、値が大きいほど剛性が高く、振動は速くなる。例えば、アルミニウム(G≈26GPa)より鋼(G≈79GPa)の方が、同じ形状なら約3倍ねじり剛性が高く、振動数も約1.7倍になるんだ。軽さ(密度)は、ここでは軸の慣性モーメントを無視しているので影響しない点に注意してね。

第二に、「円板の質量」と「円板の半径」は独立したパラメータではないということ。実際の設計では、質量 $m$ を固定して半径 $R$ を変えると、厚さや材質も変わるのが普通だよね。このツールでは $I_d = \frac{1}{2}m R^2$ で計算するから、半径を2倍にすると慣性モーメントは4倍になり、振動は一気に遅くなる。質量だけを2倍にした場合と、結果が大きく異なるから確認してみて。

最後に、実務での落とし穴。このモデルは「一自由度」の理想的なねじり振り子だから、軸の質量は無視し、円板は剛体と仮定している。でも実際の長いシャフトでは、軸自体にも分布質量があって「連続体」として振動する。すると、基本振動数だけでなく、2次、3次…という高次のねじり振動モードも存在するんだ。ツールで計算した振動数は、あくまで非常に剛性の高い(短い・太い)軸に対する一次近似値だと理解しておこう。

使い方ガイド

  1. 軸材料を選択(鋼:G=80GPa、アルミ:G=28GPa)し、軸径d(mm)と長さL(mm)を入力して極断面2次モーメントJ=πd⁴/32を自動計算
  2. 円板の質量m(kg)と半径r(mm)を入力し、回転慣性モーメントId=m·r²/2を算出
  3. ねじり剛性kt=G·J/Lから固有振動数fn=√(kt/Id)/(2π)をHz単位で求め、周期T=1/fnを確認
  4. 初期角変位θ₀(度)を設定するとθ(t)=θ₀·cos(2πfn·t)の時間波形がリアルタイム表示される

具体的な計算例

クランクシャフト設計例:鋼製軸(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以上の確保が重要。

実務での注意点

  1. プロペラ軸の場合、軸の回転質量も無視できず、有効慣性モーメントは円板単体より増加するため、Idに軸の寄与分0.1~0.2kg·m²を加算して安全評価する
  2. 減衰がない理想系での計算のため、実機では材料内部摩擦や軸受隙間による減衰で振幅は減少。Q値≈5~10程度見込んで、共振幅を制限
  3. 温度上昇でG値は低下(鋼は100℃で約1%低下)し、fnが減少するため、高温環境ではG値補正が必須
  4. 軸の応力集中係数Kt=1.5~2.0が存在する場合、実際のフィレット部ねじり応力τ=M·r/J×Ktで許容値を超過しないか検証