抗力・揚力計算機 戻る
流体工学

抗力・揚力計算機 — 球・円柱・平板の流体抵抗

Reynolds数依存のCd曲線から抗力・揚力・終端速度・カルマン渦周波数をリアルタイム算出。形状・流体・速度を変えて設計に活かそう。

物体形状と流体条件
物体形状
代表長さ D0.10 m
流速 U
物体長さ L(円柱・平板)1.00 m
迎角 α(平板・流線形)
°
流体種別
表面粗さ ε/D
計算結果
計算結果
Reynolds数 Re
抗力係数 Cd
抗力 FD (N)
揚力 FL (N)
終端速度 (m/s)
Strouhal数 / 渦周波数
Cd − Re 曲線(対数スケール)
抗力 FD − 流速 U 曲線
理論・主要公式

$$F_D = C_D \cdot \frac{1}{2}\rho U^2 A$$

$$F_L = C_L \cdot \frac{1}{2}\rho U^2 A$$

Stokes領域(Re < 1):$C_D = 24/Re$

カルマン渦周波数:$f = St \cdot U / D$,$St \approx 0.2$

抗力・揚力計算機とは

🙋
このシミュレーターで計算できる「抗力係数Cd」って何ですか?教科書で見たけど、具体的にどう使うんですか?
🎓
大まかに言うと、物体が受ける抵抗の大きさを表す無次元の数字だよ。例えば、同じ速度で動く球と円柱では、受ける力が大きく異なるよね。その違いをCdが表しているんだ。このツールでは、上の「物体形状」を球から円柱に変えるだけで、Cdが大きく変わって、計算される抗力も変わるのがすぐにわかるよ。
🙋
え、Cdって変わるんですか?「流速U」のスライダーを動かしたら、Cdの値も変わってる!なんで速度が変わると抵抗の係数も変わるんですか?
🎓
良いところに気づいたね!それは「レイノルズ数Re」が変わるからなんだ。Reは流れの状態(なめらかか乱れてるか)を決める重要な数で、これが変わると物体のまわりの流れのパターンが変わる。例えば、球の場合、Reが小さい(ゆっくり)とCdは大きく、Reが大きくなる(速く)とCdは小さくなって、あるところで急にガクンと落ちる「抗力危機」って現象も起こるんだ。シミュレーターで流速をゆっくりから一気に速くしてみて、Cdのグラフがどう変わるか確かめてみよう。
🙋
「終端速度」って計算結果に出てますけど、あれは何に使うんですか?あと、円柱を選ぶと「カルマン渦周波数」も出てきますね。
🎓
終端速度は、例えば雨粒が空から落ちてくる最終的な速度を予測するのに使うよ。重力と抗力が釣り合ってこれ以上速くならない状態だ。シミュレーターで流体を「空気」から「水」に変えてみると、終端速度が大きく変わるのがわかる。カルマン渦は、円柱や橋脚の後ろにできる交互に剥がれる渦で、この渦が発生する周波数が構造物を揺らす原因になる。実務では、煙突や海底パイプラインがこの振動で壊れないように設計する時に、この計算が非常に重要になるんだ。

よくある質問

本ツールでは、レイノルズ数に依存した標準的な実験データ(球は標準抗力曲線、円柱はRe=10^3〜10^6の遷移域を含む)を参照してCdを自動設定しています。形状選択時に適切な曲線が適用されるため、流速やサイズを変えるとリアルタイムにCdが更新されます。
物体の質量と流体の密度・粘度、物体の代表寸法を入力してください。重力と抗力が釣り合う速度を反復計算で求めます。球や円柱ではCdが速度(Re)に依存するため、収束計算が必要ですが、本ツールが自動で行います。
円柱形状を選択し、レイノルズ数が約47〜10^5の範囲にある場合に表示されます。この範囲で規則的な渦放出が発生するため、ストローハル数(約0.2)を用いて周波数を算出します。範囲外では渦放出が不安定または消失するため表示されません。
抗力・揚力はニュートン[N]、終端速度はm/s、渦周波数はHzで表示します。数値は有効数字3桁で丸めています。ただし、Cd曲線は実験近似値のため、実際の流れとは誤差が生じる可能性があります。設計の目安としてご利用ください。

実世界での応用

自動車・航空機の設計:車体や翼の形状を決める際、このツールのような計算で抗力(空気抵抗)と揚力を評価します。特に「迎角α」を変えて揚力係数がどう変わるか調べることは、翼の設計に直結します。

構造物の風荷重・流体力評価:高層ビル、煙突、橋脚、風力発電の支柱などが風や水流から受ける力を算定するために使用されます。円柱形状の場合、カルマン渦による振動(渦励振)の周波数を予測し、共振を避ける設計が必須です。

粒子の沈降・分離プロセス:水処理施設での沈殿池設計や、化学プラントでの粒子分離において、球体粒子の終端速度の計算は粒子の分離効率を決める基礎となります。流体を「水」に設定し、密度や直径を変えて終端速度を比較できます。

スポーツ工学:ゴルフボールのディンプル(表面の小さなくぼみ)は、あえて表面粗さを導入することで「抗力危機」のレイノルズ数領域を早め、抵抗を減らす効果があります。ツールの「表面粗さε/D」パラメータを変えると、この影響をシミュレートできます。

よくある誤解と注意点

この手の計算で最初に陥りがちなのが、「代表面積A」と「代表長さD」の取り方です。ツールでは自動で選んでくれますが、自分で計算する時は要注意。例えば、平板の抗力計算で代表面積を「表面積」にしてしまうと、実際より何倍も大きな抗力が出て大混乱します。正しくは流れに対して垂直な「前面投影面積」です。円柱も同様で、流れに垂直に立てた場合は直径×長さの投影面積、流れに平行に寝かせた場合はまた別の考え方になります。

次に、「係数Cdは形状固有の定数ではない」という点。これは本当に大事。ある教科書で「球のCdは0.47」と覚えて、全ての流速で使ってしまう間違いが多発します。このツールで流速スライダーを動かせば一目瞭然ですが、Cdはレイノルズ数Reに依存してガラリと変わります。実務では、対象がどのRe領域で動作するのかをまず見極め、その領域に対応したCdの値や相関式を使わないと、見積もりが大きく外れます。

最後に、計算結果の現実解釈。ツールは理想化された一様流中の単独物体を想定しています。例えば、車のボンネット上のアンテナの抗力を計算する時、アンテナ単体のCdを使って「終端速度」を出しても、実際は車体による流れの干渉で大きく異なる値になります。また、カルマン渦周波数は「剥離周波数」を理論的に与えますが、実際の構造物では固有振動数と一致した時にのみ大きな振動(共振)が起きるので、両方を比較することが設計の肝になります。

使い方ガイド

  1. 球・円柱・平板から対象形状を選択し、代表寸法D(mm)を入力します。例:球径50mm、円柱直径25mm
  2. 流体種(空気・水・油)を選択するか、密度ρ(kg/m³)と動粘性係数ν(m²/s)を直接指定します。水の場合ρ=1000kg/m³、ν=1.0×10⁻⁶m²/s
  3. 流速U(m/s)と迎え角α(°)を入力すると、Reynolds数Re、抗力係数Cd、抗力FD、揚力FL、カルマン渦周波数をリアルタイム計算します

具体的な計算例

直径50mmの球が20℃の水中で流速3.0m/sで流動される場合:Re=150,000、Cd=0.47、抗力FD≈106N。同じ球が空気中(ρ=1.2kg/m³)で風速20m/sの場合:Re=67,000、Cd=0.41、抗力FD≈2.1N。直径20mmの円柱が垂直配置(α=90°)で同一風速時:Re=26,800、Cd=1.2、抗力FD≈0.86N、Strouhal数St=0.21からカルマン渦周波数f=210Hzが算出されます。

実務での注意点

  1. Re<1(Stokes領域)では式Cd=24/Reを用いて精密計算;Re=100~10⁵ではCd値の非線形性が顕著で、迎え角変化時に揚力が急激に増減するため設計マージンを確保してください
  2. 円柱のカルマン渦周波数は共振を引き起こしやすく、スタックスが低周波(5~50Hz)の場合は機械疲労を検討します
  3. 粘度が大きい油(ν=100×10⁻⁶m²/s)での計算では、層流支配的となりCd値が単調減少するため、乱流係数表を参照せず本ツールの曲線依存設定で検証してください