空力弾性・フラッター速度計算機 戻る
空力弾性シミュレーター

空力弾性・フラッター速度計算機

翼断面の曲げ・ねじり2自由度モデルからフラッター速度と発散速度をリアルタイム計算。V-g線図・V-ω線図で安定限界を可視化し、翼断面の振動アニメーションで連成挙動を直感的に確認しよう。

翼断面パラメータ
コード長 c
m
曲げ固有振動数 fh
Hz
ねじり固有振動数 fθ
Hz
空気密度 ρ
kg/m³
密度比 μ = m/(πρb²)
回転半径比 rα
偏心距離 e/b
計算結果
— m/s
フラッター速度 Vf
— m/s
発散速度 Vd
換算振動数 k
周波数比 fh/fθ
翼断面の振動アニメーション(沈下 h + ねじり α の連成)
V-g 線図(減衰比 vs 速度)
V-g 線図(フラッタ速度)
V-ω 線図(振動数 vs 速度)
V-g線図
理論・主要公式
$V_f \approx \omega_\theta \cdot b \cdot r_\alpha \sqrt{\dfrac{\mu}{1+(2\pi/k_c)}}$
$V_{div}= \omega_\theta \cdot b \cdot \sqrt{\dfrac{\mu r_\alpha^2}{2\pi e/b}}$
b = c/2(セミコード), k = ωb/V(換算振動数)

空力弾性・フラッター速度計算機とは

🙋
「フラッター」って何ですか?飛行機の翼が壊れるって聞いたことがあります。
🎓
大まかに言うと、翼の「曲げ」と「ねじれ」が空気力と連動して、ある速度を超えると急激に振動が大きくなる現象だよ。1940年にタコマナローズ橋が崩壊した原因もこれだ。このシミュレーターでは、コード長や固有振動数を変えて、その危険な速度(フラッター速度)がどう変わるか、すぐに計算してグラフで見られるんだ。
🙋
え、橋も関係あるんですか!じゃあ、どうやったらフラッターを防げるんですか?
🎓
実務で多いのは、翼の「ねじり剛性」を高めてねじれ固有振動数を上げることだね。このツールで「ねじり固有振動数 fθ」のスライダーを右に動かしてみて。V-g線図のフラッター速度が右に移動して、より高い速度で発生するようになるのがわかるよ。
🙋
「静的発散速度」ってのも出てきますけど、フラッターとどう違うんですか?
🎓
発散は「静的な」ねじれ不安定で、フラッターは「動的な」振動不安定だ。例えば、発散は速度が上がるにつれて翼がグイグイねじれていって、ある点で復元できなくなる現象。ツールで「偏心距離 e/b」を大きくすると、発散速度が下がって、フラッターより先に危険になる様子が確認できる。両方の危険速度を常に比較するのが設計のポイントだ。

よくある質問

フラッター速度は曲げとねじりの連成振動が発散する速度で、V-g線図の減衰率がゼロを横切る点で判定します。発散速度はねじり剛性のみで決まる静的不安定現象で、空力モーメントが弾性復元力を上回る速度です。本ツールでは両方を自動計算し、グラフ上にマーカー表示します。
V-g線図は横軸に速度、縦軸に減衰率gをプロットし、g=0を超えると不安定(フラッター)です。V-ω線図は速度に対する振動数変化を示し、曲げとねじりのモードが接近・交差する付近でフラッターが発生しやすくなります。両グラフを併せて確認することで安定限界を正確に把握できます。
実際の翼設計値や実験データを参考に、質量、剛性、重心位置(静的不つりあい)、慣性モーメントを入力してください。特に静的不つりあいが大きいほど連成が強まりフラッター速度が低下します。空力係数(揚力傾斜など)は翼型に応じて適切な値を設定することで、より現実的な予測が可能です。
翼断面の曲げ(上下運動)とねじり(回転運動)の連成挙動をアニメーションで可視化します。速度を変化させると振動の位相差や振幅比が変わり、フラッター付近では両モードが同位相で発散する様子を直感的に確認できます。設計上の問題点を視覚的に把握するのに役立ちます。

実世界での応用

航空機の主翼・尾翼設計:フラッター速度は常に最大運用速度(Vmo)よりも十分高いマージンを持って設計されます。特に、高速飛行時に発生する遷音速フラッターは、コンピュータシミュレーション(CFD/CSD連成解析)とともに、このような簡易ツールによるパラメータスタディで初期設計を固めます。

ブレード類の設計(タービン、プロペラ、ヘリコプター):回転するブレードは遠心力による剛性上昇効果がありますが、それでも空力弾性不安定は重大な問題です。特に、ヘリコプターのローターブレードのフラッターは「ストールフラッター」として知られ、設計パラメータの最適化が不可欠です。

長大橋梁・煙突・高層建築物:タコマ橋の教訓後、あらゆる細長い構造物の風致振動評価に空力弾性の考え方が取り入れられました。橋の桁断面の形状(デッキ)は、フラッター発現速度を極力上げるように風洞実験と数値計算で検証されます。

風力発電タービンブレード:大型化が進むブレードは、軽量化と剛性確保のトレードオフの中で設計されます。運転中に発生する乱流や、ブレード自身の変形と空気力の連成(気弾性)を考慮したフラッター解析は、安全性と信頼性を確保するための標準的な設計プロセスです。

よくある誤解と注意点

このツールを使い始める際、特にCAE初心者が陥りがちなポイントがいくつかあるよ。まず大きな誤解は、「計算結果がそのまま安全マージンだ」と思ってしまうこと。このシミュレーターはあくまで2次元断面の最も基本的な理論モデルに基づいている。実際の翼は3次元構造で、複数のモードが複雑に連成する。例えば、ここで出てくるフラッター速度 $V_f$ は、実機ではさらに30%以上の安全マージンを見込んで設計されることが多い。ツールの結果は「傾向を理解し、パラメータの影響を比較する」ためのものと心得よう。

次に、パラメータ設定の落とし穴。密度比 $\mu$ を安易にデフォルト値のままにしないこと。これは「翼の質量と空気の質量の比」で、設計思想を大きく反映する。例えば、軽量なグライダー($\mu$が小さい)と重い戦闘機($\mu$が大きい)では、フラッター特性が大きく異なる。実際のプロジェクトでは、 予想される飛行高度での空気密度 から計算した正確な $\mu$ を使うのが鉄則だ。

最後に、「静的発散とフラッター、どちらか一方だけ見ればいい」という考えは危険。ツールで「偏心距離 e/b」を0.2から0.4に上げてみて。発散速度 $V_{div}$ がガクンと下がり、フラッターより先に危険になることが一目でわかるよね。実務では、全ての飛行条件で両方の危険速度をチェックし、より低い方に注目して設計を進める。これが「空力弾性トリオ(フラッター、発散、操縦反効)」を統合的に考える第一歩だ。

使い方ガイド

  1. コード長(翼弦長)をmm単位で入力します。典型的な小型航空機の主翼は200~400mm、大型旅客機は3000~4000mmです。
  2. 曲げ固有振動数fhと捻じり固有振動数fθをHz単位で設定します。fhは通常1~5Hz、fθはfhより高い周波数(例:fh=3Hz、fθ=8Hz)に設定されます。
  3. 空気密度ρ(kg/m³)を入力します。海面高度では1.225kg/m³、高度10000mでは0.905kg/m³です。
  4. V-g線図(速度-減衰曲線)を確認し、減衰がゼロになる速度(フラッター速度Vf)と静的発散速度(Vd)を読み取ります。

具体的な計算例

3000mm主翼、fh=2.5Hz、fθ=6.8Hz、ρ=1.225kg/m³の場合を想定します。換算振動数k=πfθc/2V(cはコード長)の関係から、V=150m/sにおけるkは約0.284です。古典的Theodorsen理論を適用した2自由度フラッター方程式より、フラッター速度Vfは約142m/s、静的発散速度Vdは約187m/sと計算されます。この値は市街地巡航速度(約100m/s)を上回る安全マージンを確保しています。

実務での注意点

  1. 後退角を持つ翼や胴体接合部の剛性低下は固有振動数を低下させ、フラッター速度を低下させるため、正確な構造解析が必須です。
  2. 飛行高度の増加に伴いρが低下するため、同じ指示速度(IAS)でもVfは高度とともに増加します。高高度飛行時は真対気速度(TAS)の監視が重要です。
  3. 風洞試験や飛行試験で得られた実測減衰値とV-g線図を比較し、Theodorsen近似の妥当性を検証してください。特に高レイノルズ数域では非線形効果が顕著になります。
  4. 橋梁設計では周辺風速、橋長、質量分布により固有振動数が決まるため、風速ハザードマップと照合し安全係数(通常1.5倍以上)を設定します。