フラッター速度計算 戻る
航空宇宙

フラッター速度計算(V-gダイアグラム)

翼断面の曲げ・ねじり2自由度フラッターを準定常空力理論で解析。各パラメータを変えてフラッター速度とV-gダイアグラムをリアルタイムで確認。

翼断面パラメータ
半コード長 b (m)
m
曲げ固有振動数 ω_h (rad/s)
rad/s
ねじり固有振動数 ω_α (rad/s)
rad/s
質量比 μ = m/(πρb²L)
振動数比 r = ω_h/ω_α
弾性軸位置 a(無次元)
フラッター解析結果
計算結果
V_flutter (m/s)
V* (換算速度)
ω_flutter (rad/s)
f_flutter (Hz)
V-g/V-f 線図

翼断面の曲げ変位h・ねじり角αと弾性軸位置

理論・主要公式
換算速度:$V^* = \dfrac{V}{b\,\omega_\alpha}$

準定常揚力:$L = \pi\rho b^2(b\omega_\alpha)^2 V^{*2}\cdot \alpha$
準定常モーメント:$M = \pi\rho b^3(b\omega_\alpha)^2 V^{*2}\cdot(\tfrac{1}{2}+a)\alpha$

固有値問題:
$\det\!\left[\mathbf{M}^{-1}\mathbf{K}_{aero}(\omega,V^*) - \omega^2\mathbf{I}\right]=0$

フラッター:$\mathrm{Im}(\omega) < 0$

フラッター速度計算(V-gダイアグラム)とは

🙋
フラッターって何ですか?飛行機の翼が壊れるって聞いたことがあります。
🎓
大まかに言うと、構造物が空気の流れの中で自ら激しく振動し始め、ついには壊れてしまう現象だよ。1940年にアメリカのタコマナローズ橋が風で大きく揺れて崩落したのは有名な例だね。このシミュレーターでは、翼の断面モデルでその危険な速度を計算できるんだ。まずは左側の「半コード長 b」のスライダーを動かしてみて。
🙋
V-gダイアグラムのグラフで、線が0を横切るところがあるけど、あれがフラッター速度なんですか?
🎓
その通り!グラフの縦軸gは「ダンピング」、つまり振動を抑える力だ。これが正の間は振動は収まるけど、飛行速度Vが上がってgが0を下回る(負になる)と、振動がどんどん増幅して発散する。この境目の速度がフラッター速度だ。例えば「質量比 μ」を大きく(重く)すると、この線が右に動いてフラッター速度が上がるのが確認できるよ。
🙋
「弾性軸位置 a」ってパラメータもあって、これを変えるとグラフがすごく変わりますね。実務ではどういう意味があるんですか?
🎓
いいところに気づいたね。aは翼の前縁から「弾性軸」という、曲げとねじりの中心線までの距離(無次元化)だ。実務では、主翼の重量配分や燃料タンクの位置でこれが決まる。aを大きくして軸を後ろにずらすと、多くの場合フラッター速度が下がる。設計ではこのパラメータを慎重に選んで、十分な安全マージン(通常15%以上)を確保するんだ。

よくある質問

「g」は構造減衰比を表します。V-gダイアグラムでは、各速度におけるフラッター振動の実効的な減衰比gをプロットしており、gが正から負に転じる速度がフラッター発生速度です。g=0が安定限界を示すため、設計ではgが負にならない速度範囲で運用します。
まず、変更したパラメータがフラッターに影響する物理量(曲げ剛性、ねじり剛性、重心位置、空力中心など)であるか確認してください。特に、質量や剛性の比率が変わらないとフラッター速度は変化しにくいです。また、リアルタイム更新が有効になっているか、ブラウザの再描画を試してください。
いいえ、あくまで教育・初期検討用です。準定常空力理論は非定常効果や3次元流れを簡略化しており、実際のフラッターはより複雑です。設計の目安として活用し、最終的には風洞試験や高精度な非定常CFD解析で検証する必要があります。
一般的に、曲げとねじりの固有振動数が近いほどフラッター速度は低下しやすくなります。設計ではねじり振動数を曲げより20~30%以上高く設定することが多いです。本ツールで振動数比を変えながらV-g線図を確認し、安全なマージンを確保する値を探すことを推奨します。

実世界での応用

航空機の主翼・尾翼設計:最も重要な応用分野です。新規設計した翼のフラッター速度を計算し、想定最大飛行速度を大きく上回ることを確認します。高高度(空気密度ρが低い)でも安全であるか、質量比μの変化を考慮して検証します。

ブレード付き機械(タービン、プロペラ):ジェットエンジンのタービンブレードや大型プロペラは、高速回転による相対風でフラッターが発生するリスクがあります。設計段階でねじり振動数ω_αなどを調整し、危険速度域を運転範囲外に設定します。

長大橋梁・煙突などの風致振動解析:タコマ橋の教訓から、風による自励振動(フラッターの一種)の検討は必須です。橋の桁断面のねじり剛性や質量バランス(弾性軸位置aに相当)が、風の安定性を大きく左右します。

F1マシンのウイングやソーラーカーのボディ:地上を高速走行する車両の空力部品も対象です。軽量化が進むと質量比μが小さくなり、フラッターリスクが高まるため、CFDと連成した空力弾性解析が行われます。

よくある誤解と注意点

このツールを使い始める際、特にCAE初心者が陥りがちな落とし穴がいくつかあります。まず第一に、「パラメータを個別にしか見ていない」という点です。例えば、「質量比μを大きくすればフラッター速度は上がる」と学ぶと、無闇に重くすれば良いと考えがちです。しかし、実際の設計では重量増加は燃費悪化に直結します。さらに、質量を増やすと慣性モーメントも変わり、ねじり振動の特性が変化するため、単純な右シフトにならないこともあります。パラメータは相互に連関していることを常に意識しましょう。

次に、「準定常空力理論」の限界を理解していないことです。このツールの計算はシンプルで強力ですが、非定常な渦の剥離や超音速流れは考慮できません。例えば、遷音速域(マッハ0.8前後)での「バフェット」や「ショックインダクストフラッター」の解析には不向きです。実務では、このような簡易ツールでトレンドを掴んだ後、より高精度な非定常CFD(計算流体力学)と連成解析に進むのが一般的です。

最後に、安全マージンの過小評価です。シミュレーションで「フラッター速度 = 500kt」と出たからといって、最大飛行速度を480ktに設定するのは危険です。材料のばらつき、製造誤差、経年劣化、計算モデルの不確実性を考慮し、通常は15%から20%の安全マージンを見込みます。つまり、この例では設計最大速度を425kt以下に抑える検討が必要になります。

使い方ガイド

  1. 曲げ剛性EI(N・m²)とねじり剛性GJ(N・m²)を入力。アルミ翼の場合EI=8×10⁵、GJ=1.2×10⁵が目安です
  2. 曲げ周波数ωb(rad/s)とねじり周波数ωt(rad/s)を設定。ωb=15~25 rad/s、ωt=30~50 rad/sが一般的です
  3. 翼の質量比μ(無次元)を入力。μ=翼質量/(空気密度×弦長×半幅)で計算します。μ=5~15が標準的です
  4. シミュレータが準定常空力理論でV-gダイアグラムを生成し、フラッター速度Vf(m/s)と換算速度V*を自動計算します

具体的な計算例

航空機の主翼例:EI=1.2×10⁶ N・m²、GJ=8×10⁴ N・m²、ωb=18 rad/s、ωt=35 rad/s、μ=8、巡航高度5000m(ρ=0.736 kg/m³)。曲げ・ねじり連成により計算結果Vf=145 m/s、V*=110 m/sが得られます。この値が設計速度Vd=125 m/sを超えるため安全マージン15%を確保できます

実務での注意点

  1. ωb<ωtの場合、曲げ優先フラッターが発生します。翼設計時にωtを意図的に高める(ねじり剛性向上)がフラッター速度上昇の鍵です
  2. μが小さい(軽量翼)ほどフラッター速度が低下。超音速機はスキン厚さ増加でμを高めています
  3. 高度変化で空気密度が変わるため、V*(換算速度)は高度に無関係ですが、Vfは密度に反比例して変化します
  4. 準定常理論の適用範囲はマッハ数0.3~0.7です。超音速・低速域ではより精密な非定常空力解析が必須です