查看“︁NACA翼型”︁的源代码

[[File:Airfoil geometry.svg|thumb|right|350px|翼型参数：1.零升力线；2.前缘；3.鼻圆；4.最大厚度；5.弯度；6.上表面；7.后缘；8.中弧线；9.下表面]]
'''NACA翼型'''是[[美国国家航空咨询委员会]]（NACA）开发的一系列[[翼型]]。每个翼型的代号由“NACA”这四个字母与一串数字组成，将这串数字所描述的几何参数代入特定方程中即可得到翼型的精确形状。

== 四位数翼型族 ==
NACA四位数翼型的几何参数为：<ref>E.N. Jacobs, K.E. Ward, & R.M. Pinkerton. [http://hdl.handle.net/2060/19930091108 NACA Report No. 460, "The characteristics of 78 related airfoil sections from tests in the variable-density wind tunnel"] {{Wayback|url=http://hdl.handle.net/2060/19930091108 |date=20191018162335 }}. NACA, 1933.</ref>
* 首位数字代表最大[[弯度]]占[[翼弦|弦长]]的百分比。
* 第二位数字代表最大弯度距机翼[[前缘]]的距离占弦长的十分之几。
* 后两位数字代表机翼最大厚度占弦长的百分比。<ref>"Fundamentals of aerodynamics", John D. Anderson,Jr, third ed, chap 4</ref>

例如，NACA 2412翼型的最大弯度为弦长的2%，位于距前缘0.4弦长处，而其最大厚度为弦长的12%。四位数翼型默认最大厚度位于距前缘0.3弦长处。

又如，NACA 0015翼型是对称翼型：00表示其没有弯度，15则表示该翼型的最大厚度为弦长的15%。

=== 对称四位数翼型方程 ===
[[File:NACA0015 a.png|thumb|right|400px|NACA 0015翼型]]
对于对称的四位数翼型NACA 00xx，其中“xx”表示最大厚度占弦长的百分比。其形状方程为：
<ref name="Moran_7">{{Cite book 
 | title = An introduction to theoretical and computational aerodynamics
 | first = Jack 
 | last = Moran
 | publisher = Dover 
 | year = 2003
 | isbn = 0-486-42879-6
 | page = [http://books.google.com/books?id=4eVP3yWZ1LgC&pg=PA7 7]
 |url=http://books.google.com/?id=4eVP3yWZ1LgC}}</ref>

:<math> y_t = 5tc\, \left[ 0.2969 \sqrt{\frac{x}{c}} + (-0.1260) \left(\frac{x}{c}\right) + (-0.3516) \left(\frac{x}{c}\right)^2 + 0.2843 \left(\frac{x}{c}\right)^3 + (-0.1015) \left( \frac{x}{c} \right)^4 \right],</math><ref>{{Cite web |url=http://www.aerospaceweb.org/question/airfoils/q0041.shtml |title=Aerospaceweb.org {{!}} Ask Us - NACA Airfoil Series<!-- Bot generated title --> |accessdate=2015-02-17 |archive-date=2020-08-16 |archive-url=https://web.archive.org/web/20200816064455/http://www.aerospaceweb.org/question/airfoils/q0041.shtml |dead-url=no }}</ref><ref>{{Cite web |url=http://www.fges.demon.co.uk/cfd/naca.html#07 |title=存档副本 |accessdate=2015-02-17 |archive-date=2010-02-03 |archive-url=https://web.archive.org/web/20100203103956/http://www.fges.demon.co.uk/cfd/naca.html#07 |dead-url=no }}</ref>

其中：
* ''c''表示[[翼弦|弦长]]
* ''x''表示沿弦长的距离（从0到''c''）
* <math>y_t</math>表示''x''处厚度的一半
* ''t''表示最大厚度占弦长的比例（100''t''即NACA四位数翼型的后两位数字）

需要注意的是，此方程中(''x''/''c'')=1处（后缘）的厚度并不完全为零。如果需要零后缘厚度的话，可以改动其中一个系数使得所有系数之和为零。改动最后一个系数对于翼型的整体形状影响最小。

翼型前缘则近似一个半径为<math>r = 1.1019 t^2</math>的圆柱。

翼型上下表面的坐标<math>(x_U,y_U)</math>与<math>(x_L,y_L)</math>分别为：
:<math>x_U = x_L = x, \qquad y_U = +y_t, \qquad y_L = -y_t.</math>

=== 有弯度四位数翼型方程 ===
[[File:NACA 2412.png|thumb|right|400px|NACA 2412翼型（红色为中弧线，紫色为对称0012翼型）]]
生成非对称的NACA四位数翼型所使用的方程与00xx对称翼型的方程相同，只不过非对称翼型的中弧线是弯曲的。其中弧线方程为：<ref name="Moran_7"/>

:<math>y_c = \left\{\begin{array}{ll}
\displaystyle{m\, \frac{x  }{ p     ^2} \left( 2\, p\, - \frac{x}{c}   \right)}, & 0 \leq x \leq pc \\
\\
\displaystyle{m\, \frac{c-x}{(1 - p)^2} \left( 1 + \frac{x}{c} - 2\, p \right)}, & pc \leq x \leq c
\end{array} \right. </math>

其中：
* ''m''表示最大弯度（100''m''即四位数翼型的首位数字）
* ''p''表示最大弯度的位置（10''p''即四位数翼型的第二位数字）

对于有弯度翼型，由于其厚度是垂直于中弧线方向的，所以其上下表面的坐标<math>(x_U,y_U)</math>与<math>(x_L,y_L)</math>分别为：<ref>{{cite web | first=Pier | last=Marzocca | url=http://people.clarkson.edu/~pmarzocc/AE429/The%20NACA%20airfoil%20series.pdf | format=PDF | title=The NACA airfoil series | publisher=Clarkson University | accessdate=07-03-2009 | archive-date=2018-11-23 | archive-url=https://web.archive.org/web/20181123153242/http://people.clarkson.edu/~pmarzocc/AE429/The%20NACA%20airfoil%20series.pdf | dead-url=no }}</ref>

:<math>
  \begin{align}
    x_U &= x   - y_t\, \sin \theta, \qquad &
    y_U &= y_c + y_t\, \cos \theta, \\
    x_L &= x   + y_t\, \sin \theta, &
    y_L &= y_c - y_t\, \cos \theta,
  \end{align}
</math>
其中
:<math>
 \theta = \arctan{ \left( \frac{dy_c}{dx} \right)},
</math>

:<math>\frac{dy_c}{dx} = \left\{\begin{array}{ll}
\displaystyle{\frac{2m}{p^2} \left(p - \frac{x}{c} \right)}, & 0 \leq x \leq pc \\
\\
\displaystyle{\frac{2m}{(1 - p)^2} \left(p - \frac{x}{c}\right)}, & pc \leq x \leq c
\end{array} \right. </math>

== 五位数翼型族==
NACA五位数翼型可以描述更为复杂的机翼形状：<ref>E. N. Jacobs & R. M. Pinkerton 1936 Test in the variable-density wind tunnel of related airfoils having the maximum camber unusually far forward, [http://hdl.handle.net/2060/19930091610 NACA Report No. 537].</ref>
* 首位数字乘以0.15表示理想攻角下的设计最佳升力系数。
* 第二位数字乘以5表示最大弯度距前缘的距离占弦长的百分比。
* 第三位数字表示翼型是否为负弯度的。
* 第四、五位数字表示最大厚度占弦长的百分比（与四位数翼型一致）

例如，NACA 23112翼型的设计升力系数为0.3（0.15*2）、最大弯距位于15%弦长处（5*3），其为负弯度翼型，而最大厚度为弦长的12%。

翼型的中弧线方程为<ref>{{cite book|last=Abbott|first=Ira|title=Theory of Wing Sections: Including a Summary of Airfoil Data|url=https://archive.org/details/theorywingsectio00abbo|year=1959|publisher=Dover Publications|location=New York|isbn=978-0486605869|page=[https://archive.org/details/theorywingsectio00abbo/page/115 115]}}</ref>

:<math>
y_c = \begin{cases}\frac{k_1}{6}\left\{x^3-3mx^2+m^2(3-m)x\right\}, & 0<x<p \\
\frac{k_1 m^3}{6}(1-x), & p<x<1\end{cases}
</math>

其中''x''与''y''已经进行了归一化处理。需选择适当的常数''m''使得最大弯度位于<math>x=p</math>处。例如，对于230xx翼型的中弧度而言，<math>p=0.3/2=0.15</math>、<math>m=0.2025</math>。最后，选择适当的常数<math>k_1</math>可以得到理想升力系数，例如在230xx翼型中<math>k_1=15.957</math>。

== 参考文献 ==
{{reflist}}

{{飞机结构}}

[[Category:空气动力学]]
[[Category:机翼]]