NACA翼型 是美国国家航空咨询委员会 (NACA)开发的一系列翼型 。每个翼型的代号由“NACA”这四个字母与一串数字组成,将这串数字所描述的几何参数代入特定方程中即可得到翼型的精确形状。
翼型参数:1.零升力线;2.前缘;3.鼻圆;4.最大厚度;5.弯度;6.上表面;7.后缘;8.中弧线;9.下表面
NACA四位数翼型的几何参数为:[ 1]
首位数字代表最大弯度 占弦长 的百分比。
第二位数字代表最大弯度距机翼前缘 的距离占弦长的十分之几。
后两位数字代表机翼最大厚度占弦长的百分比。[ 2]
例如,NACA 2412翼型的最大弯度为弦长的2%,位于距前缘0.4弦长处,而其最大厚度为弦长的12%。四位数翼型默认最大厚度位于距前缘0.3弦长处。
又如,NACA 0015翼型是对称翼型:00表示其没有弯度,15则表示该翼型的最大厚度为弦长的15%。
NACA 0015翼型
对于对称的四位数翼型NACA 00xx,其中“xx”表示最大厚度占弦长的百分比。其形状方程为:
[ 3]
y
t
=
5
t
c
[
0.2969
x
c
+
(
−
0.1260
)
(
x
c
)
+
(
−
0.3516
)
(
x
c
)
2
+
0.2843
(
x
c
)
3
+
(
−
0.1015
)
(
x
c
)
4
]
,
{\displaystyle 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],}
[ 4] [ 5]
其中:
c 表示弦长
x 表示沿弦长的距离(从0到c )
y
t
{\displaystyle y_{t}}
表示x 处厚度的一半
t 表示最大厚度占弦长的比例(100t 即NACA四位数翼型的后两位数字)
需要注意的是,此方程中(x /c )=1处(后缘)的厚度并不完全为零。如果需要零后缘厚度的话,可以改动其中一个系数使得所有系数之和为零。改动最后一个系数对于翼型的整体形状影响最小。
翼型前缘则近似一个半径为
r
=
1.1019
t
2
{\displaystyle r=1.1019t^{2}}
的圆柱。
翼型上下表面的坐标
(
x
U
,
y
U
)
{\displaystyle (x_{U},y_{U})}
与
(
x
L
,
y
L
)
{\displaystyle (x_{L},y_{L})}
分别为:
x
U
=
x
L
=
x
,
y
U
=
+
y
t
,
y
L
=
−
y
t
.
{\displaystyle x_{U}=x_{L}=x,\qquad y_{U}=+y_{t},\qquad y_{L}=-y_{t}.}
NACA 2412翼型(红色为中弧线,紫色为对称0012翼型)
生成非对称的NACA四位数翼型所使用的方程与00xx对称翼型的方程相同,只不过非对称翼型的中弧线是弯曲的。其中弧线方程为:[ 3]
y
c
=
{
m
x
p
2
(
2
p
−
x
c
)
,
0
≤
x
≤
p
c
m
c
−
x
(
1
−
p
)
2
(
1
+
x
c
−
2
p
)
,
p
c
≤
x
≤
c
{\displaystyle 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.}
其中:
m 表示最大弯度(100m 即四位数翼型的首位数字)
p 表示最大弯度的位置(10p 即四位数翼型的第二位数字)
对于有弯度翼型,由于其厚度是垂直于中弧线方向的,所以其上下表面的坐标
(
x
U
,
y
U
)
{\displaystyle (x_{U},y_{U})}
与
(
x
L
,
y
L
)
{\displaystyle (x_{L},y_{L})}
分别为:[ 6]
x
U
=
x
−
y
t
sin
θ
,
y
U
=
y
c
+
y
t
cos
θ
,
x
L
=
x
+
y
t
sin
θ
,
y
L
=
y
c
−
y
t
cos
θ
,
{\displaystyle {\begin{aligned}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{aligned}}}
其中
θ
=
arctan
(
d
y
c
d
x
)
,
{\displaystyle \theta =\arctan {\left({\frac {dy_{c}}{dx}}\right)},}
d
y
c
d
x
=
{
2
m
p
2
(
p
−
x
c
)
,
0
≤
x
≤
p
c
2
m
(
1
−
p
)
2
(
p
−
x
c
)
,
p
c
≤
x
≤
c
{\displaystyle {\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.}
NACA五位数翼型可以描述更为复杂的机翼形状:[ 7]
首位数字乘以0.15表示理想攻角下的设计最佳升力系数。
第二位数字乘以5表示最大弯度距前缘的距离占弦长的百分比。
第三位数字表示翼型是否为负弯度的。
第四、五位数字表示最大厚度占弦长的百分比(与四位数翼型一致)
例如,NACA 23112翼型的设计升力系数为0.3(0.15*2)、最大弯距位于15%弦长处(5*3),其为负弯度翼型,而最大厚度为弦长的12%。
翼型的中弧线方程为[ 8]
y
c
=
{
k
1
6
{
x
3
−
3
m
x
2
+
m
2
(
3
−
m
)
x
}
,
0
<
x
<
p
k
1
m
3
6
(
1
−
x
)
,
p
<
x
<
1
{\displaystyle 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}}}
其中x 与y 已经进行了归一化处理。需选择适当的常数m 使得最大弯度位于
x
=
p
{\displaystyle x=p}
处。例如,对于230xx翼型的中弧度而言,
p
=
0.3
/
2
=
0.15
{\displaystyle p=0.3/2=0.15}
、
m
=
0.2025
{\displaystyle m=0.2025}
。最后,选择适当的常数
k
1
{\displaystyle k_{1}}
可以得到理想升力系数,例如在230xx翼型中
k
1
=
15.957
{\displaystyle k_{1}=15.957}
。