概要 中程度の迎角におけるNACA 0012翼型周りの流線 典型的な翼型の揚力曲線と抗力曲線 固定翼航空機 の翼と安定板、およびヘリコプターの ローターブレードは、翼型断面で作られています。翼型は、プロペラ、ファン 、コンプレッサー 、タービンにも見られます。帆も翼型であり、 センターボード 、舵 、キール などのヨットの水中面は、断面が翼型に似ており、翼型と同じ原理で動作します。泳ぐ生き物や飛ぶ生き物、さらには多くの植物や固着生物も翼型/水中翼型を利用しており、一般的な例としては、鳥の翼、魚の体、 タコノマクラ の形状などがあります。翼型の翼は、自動車 やその他の動力車両にダウンフォースを 発生させ、トラクション を向上させることができます。
平板、建物、橋のデッキなどの物体によって風が遮られると、その物体は抗力 と風に垂直な空力を受けます。しかし、これはその物体が翼型であることを意味するものではありません。翼型は非常に効率的な揚力形状であり、同じ面積の平板よりも大きな揚力を発生させることができ、抗力を大幅に少なく抑えることができます。翼型は、航空機、プロペラ、ローターブレード、風力タービンなどの設計や、航空工学の様々な用途で使用されています。
右図は、風洞 試験で得られた揚力と抗力の曲線を示しています。この曲線は、正のキャンバー を持つ翼型を表しており、迎角がゼロのときに揚力が発生します。迎角が増加すると、揚力はほぼ直線的に増加し、これを揚力曲線の傾きと呼びます。約18度でこの翼型は 失速し 、それ以降は揚力が急速に低下します。揚力の低下は、失速角で上面境界層が剥離し、上面上で大きく厚くなることによって説明できます。厚くなった境界層の 変位厚さ によって翼型の実効形状が変化し、特に実効キャンバーが 減少します。これにより、全体の流れ場が変化し、循環と揚力が減少します。また、境界層が厚くなると 圧力抗力 も大きく増加するため、失速点付近および失速点を過ぎると全体の抗力が急激に増加します。
翼型設計は空力学 の重要な側面です。さまざまな翼型は、異なる飛行モードに対応します。非対称翼型は迎角ゼロで揚力を発生させることができますが、対称翼型は曲技飛行機のように頻繁に逆さまになる飛行に適している場合があります。 エルロン 付近や翼端 付近では、対称翼型を使用することで迎角の範囲を広げ、スピン 失速を回避できます。そのため、 境界層剥離を 起こすことなく、広い迎角範囲を使用できます。亜音速翼型は丸い前縁を持ち、迎角の影響を受けにくいのが特徴です。ただし、断面は厳密には円形ではありません。翼が最大厚さに達する前に曲率半径を 大きくすることで、境界層剥離の可能性を最小限に抑えています。これにより翼が長くなり、最大厚さの位置が前縁から後方に移動します。
超音速翼型は 、形状がより角張っており、迎角に非常に敏感な鋭い前縁を持つことができます。超臨界翼型 は、最大厚みが前縁付近にあり、長い翼幅を持つことで、超音速流をゆっくりと亜音速に戻すことができます。一般的に、このような遷音速 翼型と超音速翼型は、抗力発散を 低減するためにキャンバーが小さくなっています。現代の航空機の翼は、翼幅に沿って異なる翼型断面を持つ場合があり、それぞれの翼型断面は翼の各セクションの状況に合わせて最適化されています。
可動式の高揚力装置であるフラップ (場合によってはスラット) は、ほぼすべての航空機の翼型に取り付けられています。後縁フラップは補助翼と同様の働きをしますが、補助翼とは異なり、使用しないときは翼内に部分的に格納することができます。
多くの翼では、上面付近の気流の境界層が最大厚さ点を超えると急速に乱流になり、表面摩擦抗力が増加します。 層流翼で は、最大厚さ点が翼弦に沿って、一般的な翼弦の25%の位置から前縁から60%以上までかなり後方に移動します。[ 7 ] これにより、翼のより広い範囲で滑らかな層流が維持され、抗力が大幅に減少します。ただし、表面汚染により境界層が乱れ、乱流になります。たとえば、昆虫が翼に衝突して付着すると、翼表面のくさび形の層流領域が失われます。これは、離陸速度の速い航空機にとって特に問題です。多くの昆虫は地面付近に生息しているため、飛行中に層流を維持することは困難です。[ 8 ] グライダーは 、低速で低抗力の空力構造が必要なため、層流翼型が広く採用されています。[ 7 ]
1970年代と1980年代のNASAの研究以前は、一般的な製造公差と表面の不完全性のため、層流翼の設計は実用的ではありませんでした。機械加工された金属と複合材料を使用した新しい製造方法が開発されるまでは(例えば、繊維強化プラスチック 製の翼に使用するためにフランツ・ヴォルトマン 教授によって開発された層流翼型など)。1980年代のNASAの研究は、層流翼の設計の実用性と有用性を明らかにし、亜音速の一般航空機から遷音速の大型輸送機、超音速設計に至るまで、現代の実用的な航空機の表面への層流の応用への道を開きました。[ 9 ]
翼型を定義するための様々な方式が考案されており、NACAシステムは その一例である。また、様々な翼型生成システムも用いられている。NACAシステムよりも古くから存在し、広く応用されている汎用翼型の一例として、クラークY 型が挙げられる。今日では、コンピュータプログラムを用いて特定の機能に合わせた翼型を設計することも可能となっている。
翼型用語 翼型命名法 翼型に関連するさまざまな用語は、以下のように定義されます。[ 10 ]
吸引面 (通常は上面)は、一般的に流速が高く、静圧が低い。 圧力面 (通常は下面)は、吸引面よりも比較的高い静圧を持つ。これら2つの面間の圧力勾配が、特定の翼型で発生する揚力に寄与する。 翼型の形状は、さまざまな用語で表現されます 。
前縁 とは、翼型の前方で曲率が最大(半径が最小)になる点のことである。 後縁 とは、翼型上で前縁から最も遠い点のことです。後縁における上面と下面の間の角度を後縁角 といいます。 弦線は 、前縁と後縁を結ぶ直線です。弦の長さ 、または単に弦 、c {\displaystyle c} は翼弦線の長さです。これは翼型断面の基準寸法です。 翼型厚さの異なる定義 ウィングレット 用に設計された翼型(PSU 90-125WL)翼型の形状は、以下の幾何学的パラメータを用いて定義される。
平均キャンバー線 または平均線 とは、上面と下面の中間にある点の軌跡のことです。その形状は、弦に沿った厚みの分布によって決まります。 翼型の厚さは 、翼弦に沿って変化する。厚さの測定方法は2つある。 厚さはキャンバーラインに垂直に測定されます。これは「アメリカ式」と呼ばれることもあります。 厚さは弦線に垂直に測定される。これは「英国慣習」と呼ばれることもある。 翼型の形状を表す重要なパラメータには、キャンバー と厚み があります。例えば、NACA 2415(2-4-15と読みます)のようなNACA 4桁シリーズの翼型は、0.40コード長の位置に0.02コード長のキャンバーがあり、最大厚みが0.15コード長の翼型を表しています。
最後に、流体中を移動する際の翼型の挙動を説明するために用いられる重要な概念は以下のとおりです。
空力中心 とは、ピッチングモーメントが揚力係数と迎角に依存しない、翼弦方向の位置のことである。 圧力中心と は、ピッチングモーメント が瞬間的にゼロになる翼弦方向の位置のことです。キャンバー翼型の場合、圧力中心は迎角や揚力係数の変化に応じて移動するため、固定位置ではありません。無限スパンの一様な翼の周りの二次元流れでは、揚力曲線 の傾きは主に後縁角 によって決まります。傾きは、角度がゼロのときに最大となり、角度が増加するにつれて減少します。[ 15 ] [ 16 ] 有限スパンの翼の場合、翼のアスペクト比 も曲線の傾きに大きく影響します。アスペクト比が減少すると、傾きも減少します。[ 17 ]
薄翼理論 1991年に製造されたこのデニー・キットフォックス 航空機の先端には、翼型断面が展示されている。カモフKa-26 ヘリコプターの下部ローターブレードの翼型プロファイル薄翼理論は、非圧縮性で 非粘性の流れに対する 迎角 と揚力を関連付ける翼の単純な理論です。これは、ドイツの数学者マックス・ムンク によって考案され、1920年代にイギリスの空気力学者ヘルマン・グラウエルト らによってさらに改良されました 。この理論は、翼の周りの流れを薄い翼の周りの2次元の流れとして理想化します。これは、厚さがゼロで翼幅 が無限の翼を扱うものとして考えることができます。
薄翼理論は、二次元非粘性流における翼の以下の重要な特性に対する確固たる理論的基礎を提供したため、当時特に注目されました。
対称翼型では、圧力中心 と空力中心 は一致し、前縁から翼弦長のちょうど4分の1の位置にある。 キャンバー付き 翼型では、空力中心は前縁から翼弦のちょうど4分の1後方に位置するが、迎角が変化すると圧力中心の位置も変化する。揚力係数と迎角の関係を示す 直線の傾きは2 π {\displaystyle 2\pi \!} 単位/ラジアン。(3)の結果として、無限翼幅の薄肉対称翼型の断面揚力係数は次のようになる。
c l = 2 π α {\displaystyle \ c_{l}=2\pi \alpha } どこc l {\displaystyle c_{l}\!} は断面揚力係数です。 α {\displaystyle \alpha \!} は、翼弦線を基準とした 迎角 (ラジアン)です。(上記の表現は、キャンバー付き翼型にも適用可能で、α {\displaystyle \alpha \!} (迎角は、翼弦線ではなく、ゼロ揚力線 を基準として測定される。)
また、(3)の結果として、無限翼幅のキャンバー翼型の断面揚力係数は次のようになる。
c l = c l 0 + 2 π α {\displaystyle \ c_{l}=c_{l_{0}}+2\pi \alpha } どこ c l 0 {\displaystyle \ c_{l_{0}}} は迎角がゼロのときの断面揚力係数です。 薄翼理論では、空気は非粘性流体 であると仮定しているため、翼の失速は考慮されていません。翼の失速 は通常、典型的な翼では迎角が 10° ~ 15° の間で発生します。 [ 21 ] しかし、2000 年代半ばから後半にかけて、前縁失速の開始を予測する理論が、ウォレス J. モリス II によって博士論文で提案されました。 モリスによるその後の改良では、前縁失速現象に関する現在の理論的知識の詳細が示されています。モリスの理論では、前縁失速開始の臨界迎角は、内部流れの解で全体的な剥離領域が予測される条件として予測されます。[ 25 ] モリスの理論は、薄い翼型の周りの亜音速の流れは、翼弦の大部分の周りの外側領域と、ノーズ周辺の内側領域で記述でき、これらは漸近的に一致することを示している。外側領域の流れは古典的な薄い翼型理論によって支配されるため、モリスの方程式は薄い翼型理論の多くの要素を示している。
導出 上から下へ: 層流
乱流
亜音速流
超音速流量
薄翼理論では、(2D)翼の幅は無視できるものとし、翼自体はキャンバー線に沿って迎角 α で向き付けられた1Dブレードに置き換えられます。ブレードに沿った位置をx とし、翼の前縁で0 、後縁で c とします。翼のキャンバーdy / dx は十分に小さいと仮定されるため、 x と胴体に対する位置を区別する必要はありません。 [ 26 ]
翼型を横切る流れは、ブレードの周囲に循環を生成し、これは位置によって強度が変化する渦シート γ ( x ) としてモデル化できます。クッタ条件は γ ( c )=0 を意味しますが、強度はブレード前縁で特異であり、x ≈ 0 の場合、 γ ( x ) ∝ 1 ⁄ √ x となります。主流 V の 密度 がρ の 場合、クッタ-ジュコフスキーの定理 により、総揚力Fは [ 30 ] に比例します。ρ V ∫ 0 c γ ( x ) d x {\displaystyle \rho V\int _{0}^{c}\gamma (x)\,dx} そして、その前縁周りのモーメントMは ρ V ∫ 0 c x γ ( x ) d x 。 {\displaystyle \rho V\int _{0}^{c}x\;\gamma (x)\,dx.}
ビオ・サバールの法則 によれば、渦度γ ( x ) は流れ場を生成する。w ( x ) = 1 2 π ∫ 0 c γ ( x ′ ) x − x ′ d x ′ 、 {\displaystyle w(x)={\frac {1}{2\pi }}\int _{0}^{c}{\frac {\gamma (x')}{xx'}}\,dx'{\text{,}}} x における翼型に垂直な方向 に配置されている。翼型は不透過性の表面 であるため、流れはw ( x ) {\displaystyle w(x)} V からの逆流と釣り合わなければならない。小角近似 により、V は位置 x のブレードに対して角度α - dy ⁄ dx で傾いており、法線成分はそれに応じて( α - dy ⁄ dx ) V となる。したがって、γ ( x )は 畳み込み 方程式 を満たさなければならない。( α − d y d x ) V = − w ( x ) = − 1 2 π ∫ 0 c γ ( x ′ ) x − x ′ d x ′ 、 {\displaystyle \left(\alpha -{\frac {dy}{dx}}\right)V=-w(x)=-{\frac {1}{2\pi }}\int _{0}^{c}{\frac {\gamma (x')}{xx'}}\,dx'{\text{,}}} これは既知の量に関して一意に決定します。
明示的な解は、まず変数変換を 行うことで得られる。x = c ⋅ 1 + コス ( θ ) 2 、 {\displaystyle x=c\cdot {\frac {1+\cos(\theta )}{2}},} そして、dy / dx とγ ( x ) の両方を、修正された先頭項を持つθ の無次元フーリエ級数として展開します。 d y d x = A 0 + A 1 コス ( θ ) + A 2 コス ( 2 θ ) + … γ ( x ) = 2 ( α + A 0 ) ( 罪 θ 1 + コス θ ) + 2 A 1 罪 ( θ ) + 2 A 2 罪 ( 2 θ ) + … 。 {\displaystyle {\begin{aligned}&{\frac {dy}{dx}}=A_{0}+A_{1}\cos(\theta )+A_{2}\cos(2\theta )+\dots \\&\gamma (x)=2(\alpha +A_{0})\left({\frac {\sin \theta }{1+\cos \theta }}\right)+2A_{1}\sin(\theta )+2A_{2}\sin(2\theta )+\dots {\text{.}}\end{aligned}}} 結果として得られる揚力とモーメントは、この級数の最初の数項のみに依存する。
揚力係数は 以下を満たす。C L = 2 π ( α + A 0 + A 1 2 ) = 2 π α + 2 ∫ 0 π d y d x ⋅ ( 1 + コス θ ) d θ {\displaystyle C_{L}=2\pi \left(\alpha +A_{0}+{\frac {A_{1}}{2}}\right)=2\pi \alpha +2\int _{0}^{\pi }{{\frac {dy}{dx}}\cdot (1+\cos \theta )\,d\theta }} モーメント係数C M = − π 2 ( α + A 0 + A 1 − A 2 2 ) = − π 2 α − ∫ 0 π d y d x ⋅ コス ( θ ) ( 1 + コス θ ) d θ 。 {\displaystyle C_{M}=-{\frac {\pi }{2}}\left(\alpha +A_{0}+A_{1}-{\frac {A_{2}}{2}}\right)=-{\frac {\pi }{2}}\alpha -\int _{0}^{\pi }{{\frac {dy}{dx}}\cdot \cos(\theta )(1+\cos \theta )\,d\theta }{\text{.}}} 1/4コードポイント付近の瞬間はC M ( 1 / 4 c ) = − π / 4 ( A 1 − A 2 ) 。 {\displaystyle C_{M}(1/4c)=-\pi /4(A_{1}-A_{2}){\text{.}}} このことから、圧力中心 は「四分の一弦」点0.25 c より後方にあることがわかる。Δ x / c = π / 4 ( ( A 1 − A 2 ) / C L ) 。 {\displaystyle \Delta x/c=\pi /4((A_{1}-A_{2})/C_{L}){\text{.}}} 空力中心 とは、揚力係数の変化によってピッチングモーメントM ′ が変化 しない位置のことである。 ∂ ( C M ′ ) ∂ ( C L ) = 0 。 {\displaystyle {\frac {\partial (C_{M'})}{\partial (C_{L})}}=0{\text{.}}} 薄翼理論によれば、二次元非粘性流においては、空力中心は翼弦の四分の一の位置にある。
参考文献
引用文献 ↑ ハリデー& レスニック 1988 、p. 378 :「翼の効果は、気流に下向きの速度成分を与えることである。偏向された空気塊の反作用力は、翼に等しく反対方向の上向きの成分を与えるために作用しなければならない。」 ↑ Hall, Nancy R. 「流れの向きを変えることによる揚力」。NASAグレン研究センター。2011年7月5日のオリジナルからアーカイブ。2011年6月29日 取得。物体が流れの正味の偏向または向きを変えるように形作られたり、移動されたり、傾いたりすると、局所速度の大きさ、方向、またはその両方が変化する。速度の変化は物体に正味の力を生み出す。 ↑ 「明確な揚力を得るためには、鋭い後縁を持つ翼を使用しなければならないことは、飛行の黎明期から知られていた。」リチャード・フォン・ミーゼス (1945)『飛行理論』 第VIII.2節、179ページ、ドーバー出版、ISBN 0-486-60541-8 ↑ Babinsky 2003 、pp. 497–503 :「流線が湾曲している場合、流線に沿って圧力勾配が存在するはずです。」 1 2 ホートン、エドワード L.、カルザーズ、NB (1982)。 工学系学生のための空気力学 (第3 版)。ロンドン:アーノルド 。ISBN 978-0-7131-3433-9 。↑ Croom, CC; Holmes, BJ (1985-04-01). 層流翼用昆虫汚染防止システムの飛行評価 。 ↑ Holmes, BJ; Obara, CJ; Yip, LP (1984-06-01). "現代の航空機表面における自然層流実験" . NASA Technical Reports : 18592. Bibcode : 1984ntrs.rept18592H . ↑ Hurt, HH Jr. (1965 年 1 月) [1960 年]. 海軍航空士官のための空気力学 . 米国政府印刷局、ワシントン DC: 米国海軍、航空訓練部. pp. 21–22 . NAVWEPS 00-80T-80. ↑ Lyons, DJ、および Bisgood, PL (1945 年 1 月)。アスペクト比の小さい翼型の揚力勾配の解析 。英国航空研究評議会報告書および覚書 No. 2308 ↑ Abbott, IH、および Von Doenhoff, AE (1949)翼断面の理論 、セクション 7.4(b) ↑ Abbott, IH、および Von Doenhoff, AE (1949)翼断面理論 、第 1.3 節 ↑ スコット 2003 :「この式は、中~大アスペクト比の翼を持つ航空機にのみ適用でき、失速角までしか適用できません。一般的な航空機の構成では、失速角は通常10°から15°の間です。」↑ Ramesh, Kiran; Gopalarathnam, Ashok; Granlund, Kenneth; Ol, Michael V.; Edwards, Jack R. (2014年7月)「間欠的な前縁渦放出を伴う非定常翼流に対する新しい放出基準を用いた離散渦法」 Journal of Fluid Mechanics . 751 : 500–538 . Bibcode : 2014JFM...751..500R . doi : 10.1017/jfm.2014.297 . ISSN 0022-1120 . S2CID 121962230 . ↑ Auld & Srinivas 1995 : 「一般的な2次元翼型断面の簡単な解は、厚みの影響を無視し、平均線のみの断面モデルを使用することで得られる。...これはまた、位置の小さな変化が等価であり、 ds ≈ dx であることを意味する。」↑ Acheson, D. J. (1990). Elementary Fluid Dynamics . Oxford Applied Mathematics and Computing Science. Oxford: Clarendon Press (published 2009). pp. 140–141 , 143–145 .
一般情報源 アボット、アイラ・ハーバート ;フォン・ドーンホフ、アルバート・エドワード(1959)。翼型理論、翼型データの概要を含む 。ドーバー出版。ISBN 978-0-486-60586-9 。Auld, Douglass; Srinivas (1995). "2次元薄翼理論" .学生のための空気力学 . シドニー大学。 Babinsky, Holger (2003年11月). 「翼はどのように機能するのか?」(PDF) . Physics Education . 38 (6): 497– 503. Bibcode : 2003PhyEd..38..497B . doi : 10.1088/0031-9120/38/6/001 . S2CID 1657792 . Batchelor, George. K (1967).流体力学入門 . Cambridge UP. pp. 467–471 . Bertin, John J.; Cummings, Russel M. (2009).エンジニアのための空気力学 (第5 版). Pearson Prentice Hall. ISBN 978-0-13-227268-1 。 クランシー、LJ(1975)。空気力学 。ロンドン:ピットマン。ISBN 0-273-01120-0 。 ハリデー、デイビッド、レスニック、ロバート(1988)。物理学の基礎 (第3 版)。ジョン・ワイリー・アンド・サンズ。 Houghton, EL; Carpenter, PW; Collicott, Steven H.; Valentine, Daniel (2012).工学系学生のための空気力学 (第6 版). Elsevier. ISBN 978-0-08-096633-5 。 Morris, Wallace J. II (2009).広範囲のレイノルズ数流れにおける翼型の失速開始の普遍的予測 (博士論文). ハーバード大学. Bibcode : 2009PhDT......146M . Morris, Wallace J.; Rusak, Zvi (2013 年 10 月). 「低~中程度のレイノルズ数流れにおける翼型での失速開始」. Journal of Fluid Mechanics . 733 : 439–472 . Bibcode : 2013JFM...733..439M . doi : 10.1017/jfm.2013.440 . ISSN 0022-1120 . S2CID 122817884 . フィリップス、ウォーレン F. (2004).飛行力学 . ジョン・ワイリー・アンド・サンズ. ISBN 978-0-471-33458-3 。 Scott, Jeff (2003年8月10日). 「質問番号136:揚力係数と薄翼理論」 . ロケット科学者に聞く:空気力学. Aerospaceweb.org. Traub, Lance W. (2016年3月24日). 「翼型ヒステリシスの半経験的予測」 . Aerospace . 3 (2): 9. Bibcode : 2016Aeros...3....9T . doi : 10.3390/aerospace3020009 . ウェルトナー、クラウス;インゲルマン=サンドバーグ、マーティン(1999)。「飛行の物理学 ― 再考」。2011年9月29日時点のオリジナルよりアーカイブ。2021年4月25日 取得。
さらに読む アンダーソン、ジョン、D(2007)。空気力学の基礎 。マグロウヒル。 {{cite book}}: CS1 maint: 複数の名前: 著者リスト (リンク)Ali Kamranpay、Alireza Mehrabadi。「異なる迎角におけるNACA翼型0012の数値解析と空力係数の取得」。Journal of Mechatronics and Automation。2019; 6(3): 8–16p。 Bearman, Matt (2019). 「流れに身を任せる? 1930~1947年の層流研究における英国の貢献」The Aviation Historian (29): 74–87 . ISSN 2051-1930 .
外部リンク UIUC翼型座標データベース 翼型使用法に関する不完全なガイド、イリノイ大学 翼型および水中翼型リファレンスアプリケーション FoilSim NASA の翼型シミュレーター エアフォイル・プレイグラウンド - インタラクティブなウェブアプリ Desktopaero 翼を通過する空気の流れ(ケンブリッジ大学) DesignFOILは、登録不要で利用できる翼型生成・解析ツールです。