1、第十二章 机 翼 理 论课堂提问:雁群迁徙时为什么呈 ”人字形 ”飞行 ?机翼理论:研究支持飞机升空,水翼船飞腾的机翼理论。在航空,舰船等工程上应用最多,舵、螺旋桨,减摇鳍、水翼、扫雷展开器,研究船舶的操纵性时可以把船体的水下部分看作是一个机翼(短翼) 。此外在风扇,鼓风机,压缩机,水上运动器械如帆板,脚蹼等都与机翼理论有关。本章内容:1. 几何特性2. 流体动力特性3. 有限翼展机翼(三元机翼)本章重点:1. 机翼几何特性。2. 机翼几何特性对流体动力特性的影响。3. 下洗速度形成的概念及计算,自由涡、附着涡形成的概念。4升力线理论的概念。5. 诱导阻力的概念,诱导阻力的计算。6. 展弦比换
2、算的思路及计算。本章难点: 1. 机翼几何特性对流体动力特性的影响。2. 升力线理论的概念。3. 展弦比换算。-机翼的几何特性一、翼型(profile)翼剖面的重要参数:中线(center line) ,翼弦(chord)b,拱度(camber)f,相对拱度 f/b,展长 ,厚度 t,l相对厚度 t/b, (thicheness) ,攻角(angle of attach),翼型面积 S,展弦比 等。根据工程应用的需要,机翼的平面形状多样。展弦比 2lS对于矩形机翼 , 所以 无限翼展机翼:Slb2lb12:短翼:2, 大展弦比机翼: 2 船用舵 , 水翼0.5:57战斗机 ,轰炸24:机 ,风
3、洞试验一般采用标71准机翼 。56机翼的攻角又分为:几何攻角:来流速度 与弦线之间的夹角。 0U基本形状:后缘总是尖的(产生环量)圆前缘:减小形状阻力尖前缘:减小压缩性所引起的激波阻力或自由表面所引起的兴波阻力翼型: 几种常见的翼型翼型(美国国家航空咨询委员会(National Advisori committee for Aeronautics,简称 NACA)设计发表的)目前在舰船的舵、螺旋桨上用得较多的是 NACA 翼型系列。NACA 四组翼型:)NACA 四位数字翼型(12-2)()2)1()22fffff ffff xxxyx 该翼型系列的厚度表达式为(12-3)432507.41.
4、7580.63.0845.( xxxty 翼型系列的 % ,40 %,前缘半径 。翼型系列有九种相对厚度:3tx ,92tr前6%, 8%, 9%, 10% 12%, 15%, 18%, 21%, 24%;有三种相对拱度:0, 1%, 2%。)NACA 五位数字翼型五位数字翼型的厚度分布仍与四位数字翼型相同,都是(-)式,相对厚度有 12%,15%, 18%, 21%, 24%五种; 都是 15%;设计fx升力系数都是 0.3。)NACA 层流翼型,层流翼型在航空工程中早已受到重视。若用于船舶螺旋桨,减阻效益甚小,但对于延迟空泡的产生是有利的。 二、机翼的平面图形机翼的平面图形是多种多样的。对
5、于船用舵,舵高就是翼展。若将整个船体的水下部分看作是一个机翼,则其吃水就是翼展。展弦比:翼展的平方和机翼面积之比称(12-5)Sl2对于矩形机翼, 等于展与弦之比(12-6)bl2风洞试验用标准机翼,,(美国用,原苏联用) 。战斗机 ;轰炸机 ;水翼 ;船用舵 0.51.5。不同展弦比的机翼,其流动特性有很大差别,在理论研究方法上亦有很大不同。通常对 的机翼,称小展弦比机翼;的机翼称大展弦比机翼;,即为二元机翼(翼剖面) 。- 库塔儒可夫斯基定理 一、定理的证明用动量定理来证明该定理:不计质量力。在方向列立动量方程。通过控制面 Cr 的动量为)cossin)(cos(2 rvvV忽略扰动速度
6、Vr 和 Vs 的二阶以上小量,求出积分得 Cr 边界在向动量变化为Vov sr ()作用于控制体边界上方向的力为翼型的反作用力-L 作用于控制体边界 Cr 上流体压力在方向分量的积分为()dprsin20压力可用柏努利方程确定222 1)si()cos(21 oosor VpvVvVp 忽略扰动速度的二阶以上小量得=P 0- vrcos v sin代入压力积分,可得 Cr 上所受方向的力为()osVdpsin20将() , () , ()代入动量方程得v s vs所以 v s (d) 令 Cr 上沿顺时针方向速度环量为 ,则有 v s在无旋流场中,绕周线 Cr 的速度环量 亦即等于绕翼剖面周
7、线的速度环量 ,因此儒可夫斯基定理得证: (12-7)二、机翼绕流环量形成的物理过程静止流场中机翼加速到 Vo 的过程中环量产生的机理:a) 作包围机翼剖面并延伸到充分远的封闭流体周线 CDFE,启动前此封闭流体周线上的速度环量为零。由汤姆逊定理,此流体周线上的环量将始终保持为零。b)机翼突然启动,速度很快达 Vo,流体处处无旋。绕翼型的环量为零。后驻点不在后缘而在处,流体绕过后缘尖点流到翼背上去,T 附近速度很大,压力很低,处速度为零,压力很高,流体由流向时遇到很大逆压梯度,使边界层分离,形成起动涡。起动涡随着流体向下游运动。根据汤姆逊定理,沿流体周线 CDFE 的环量仍应为零,故绕翼剖T面
8、必将产生一速度环量,其大小与起动涡相等方向相反。由于环量的作用,后驻点向后缘点移动。不断有反时针方向的旋涡流向下游,绕机翼的环量 也不断增大,驻点不断向后缘点推移,直到后驻点推移到后缘点为止。当机翼剖面以速度继续飞行,后缘不再有旋涡脱落,环量 也不再变化, 就只与翼剖面的几何形状以及来流的速度大小与方向有关。c) 这时翼剖面上、下两股流体将在翼剖面的后缘处汇合。这时为正常飞行的有利流动图案,流体绕流过机翼时,上面的流线较密,下面的流线较稀,故上面流体的速度大、压力小,下面的流体速度小,压力大,因而产生升力。飞机机翼至少一部分是由流过上表面的空气把它吸起来的,上表面产生的负压对全部升力的贡献比下
9、表面正压力的贡献大。- 机翼的流体动力特性流体动力特性升力 L:绕流物体上、下物面上流动的不对称,引起压力的不对称,在垂直于运动方向产生的压力差。阻力 R:二元机翼的总阻力有摩擦阻力和形状阻力两部分组成。流体动力系数:是无量纲参数,主要有升力系数 ,阻力系数 力矩系数LC,R,MC(12-19)lAVLC201(12-20)D2(12-21)lbM01一、升力系数 为风洞试验求得的升力系数曲线:攻角 的增加,升力系数 按直线比例上升,达到临界攻角 时升力系数达到最大值吸 力 压 力吸 力吸 力 压 力引起的兴波阻力压 力 系 数 分 布 曲 线压 力 系 数 分 布 曲 线 Lmax。失速:攻
10、角增加到某一值升力突然减小并伴随着阻力突然增大。机翼或水翼突然丧失了支承力, 舵失去操纵作用的现象。原因:边界层分离造成。临界攻角:由实验确定,对于翼剖面一般在 1020之间。零升力角 :翼型升力为零时所对应的攻角,零升力线与弦线之间的夹角,一般约为00-2o。在这一攻角附近,机翼的阻力最小。对称翼型 , 。f 的大小在数值上约等于 大小的百分数,即f - 100% (12-22)最大升力系数:主要与翼型的相对拱度 、相对厚度 以及雷诺数有关。 随的增大而ft略有减小,这是由于大将推迟翼剖面边界层分离,从而减小边界层压差阻力的结果。升力系数曲线斜率 :反映升力系数随几何攻角的变化程度。LdC当
11、 2 时,在很大攻角范围内,升力系数为0()()LLad当 =时,升力曲线的斜率的理论值为: (1/弧度) ,()2LdC但试验结果为: (1/弧度) 。()(0.89)2LdC:二、阻力系数翼型阻力:摩擦阻力和压差阻力(亦称形状阻力)两部。摩擦阻力:由翼剖面上流体粘性切应力在翼型运动方向的投影所产生。压差阻力:由翼型表面边界层分离,或无分离情况下流动受粘性排挤产生边界层前后压力差所 t造成。 绝对值的增加而阻力系数增大,零攻角 处取极小值。CD随 的绝对值增大而增加,在 时 D取极小值。三、极曲线 ( C D关系曲线)包括了上面两条曲线的全部内容。从原点到曲线上任一点的矢径就表示了在该对应攻
12、角下的总流体动 力系数的大小和方向。矢径的斜率:该攻角下升力与阻力之比C L/CD,简称升阻比。过原点作极曲线的切线,就得飞机(或机翼)的最大升阻比,显然这是飞机最有利的飞行状态。 四、俯仰力矩系数压力中心 CP:机翼上流体动力合力与弦线的交点,或合力作用点的位置。.俯仰力矩的大小与参考点有关。参考点常有两种取法:前缘为参考点的力矩o 定义为(12-23)lAbVMC2010离前缘 1/4 弦长处的力矩: M 1/利用俯仰力矩曲线,结合相应的升阻力曲线, 可求得压力中心位置(流体合力与翼弦交点) 。图为一 NACA 对称翼型的俯仰力矩曲线,发现该翼型失速前,m 的值恒等于零,并与攻角 和升力系
13、数均无关。对于对称翼型,其压力中心恒在离前缘弦长处。这就是舵杆为什么常安装在弦长处的原因。一个优良的翼型,其压力中心位置随攻角改变移动不能太大,否则机翼的稳定性较差。- 有限翼展机翼一、有限翼展机翼的升力线理论展弦比对机翼的流体动力特性有重要的影响,因此由展弦比将机翼分成两类:大展弦比机翼。:小展弦比机翼或短翼,本节的结果只适用于大展弦比机翼。 实践表明,展弦比 时,机翼的附着涡系可用一根涡丝来代替(它可视为各 形涡的附着涡的迭合) ,这根涡丝通常称为升力线(lift line) 。以升力线为理想模型的计算机翼动力特性的理论称为升力线理论。 三元机翼的流体动力特性主要问题是翼端效应,即机翼上、
14、下翼面出现沿展向的横流,与来流合成产生自由涡,自由涡产生诱导速度(下洗速度) ,使有效来流速度改变了方向,机翼的流体动力合力在无穷远来流方向有了投影,即诱导阻力。 (插入动画附着涡,自由涡说明)二、下滑速度 滑角 导阻力升力线模型:机翼的附着涡系可由一根涡丝替代,这一涡丝称为升力线,用升力线作为机翼的理想模型来研究大展弦比机翼的流体动力特性。绝对攻角(流体动力攻角) : ,零升力线与无穷远来流 之间的夹角。a00U有效攻角 : ,有效来流速度与弦线之间的夹角。k临界攻角 :又称失速角。当机翼攻角达到某一数值时,升力突然减小,阻力急剧增加所cr对应的攻角。其大小与机翼的剖面形状,几何攻角有关。下
15、洗速度: ,翼端绕流引起后缘速度不连续,在翼后缘出现旋涡层,产生诱导2w速度,其方向向下即下洗速度。下洗角 : ,有限翼展机翼的翼端绕流产生自由涡而形成下0LCwtgU:洗速度使有效来流速度与 之间产生的夹角。诱导阻力: ,有限翼展机翼在理想流体中作等速直线运动时,下洗速度使有效iDL:来流速度与无穷远来流方向之间有了夹角。机翼的流体动力合力在无穷远来流方向有了投影,即诱导阻力。诱导阻力系数: 2LDiC三元机翼的总阻力: 又摩擦阻力,形状阻力,诱导阻力三部分组成。引入两点假定:(1)自由涡面是平面,延伸至无穷远而不翻卷成两股大涡,自由涡面旋涡角速度矢量平行来流(2)翼面上横向流动很小,任一剖
16、面处可作平面流动处理,三元效应仅考虑各翼剖面处下洗速度和下洗角的不同。涡丝在升力线上点产生的下洗速度为: 沿展向积分得整个自由涡在 y 处的诱导速度:当 y=, 上式为旁义积分,取主值为: )()()( 2/2/02/ lim ydydyd lyyll因所以而宽度为 dy 的一段机翼的二维升力为按定义升力垂直于来流诱导阻力积分的整个机翼上的升力和诱导阻力1()4ddW2()ll20kVW01()kV0tan()kdLydcosVtan()idDLWyd220()4lllliLVyDd或由此可知,要求出诱导阻力,必须要知道沿翼展的速度环量。下面来求速度环量。三、展弦比的换算:设计机翼或舵时,例如
17、需要将 NACA 翼型中 的某种翼型的 曲线换算成 与所要设计物相同的升力曲线。这就是展弦比的换算方法。设翼型、弦长相同的两个平面形状为矩形的机翼和,由式(12-58)可知,下洗角沿翼展的分布为(12-69)knAy1si)(下洗角的平均值为(12-70)2/)(ldy将(12-54)以及(12-69)式代入式(12-70)(12-71)knkn AlAl 1210si2i因此(12-72)(LC式中(12-73)12Akn对于不同平面形状的机翼,已将 的值列于表 12-1 中。)(机翼和机翼的有效攻角分别为(12-)1(11Lk C74)(12-)1(222 Lk75)要使和具有相同的升力系
18、数 ,必须让它们在相同的有效攻角下工作,即 22()()4lllli ydD几何攻角必须不同,由(12-74)和(12-75)两式相减得(12-76))1(212LC设机翼的 曲线为已知。在该曲线上任取一点,它对应的升力系数为 ,按(12-76)式算出几何攻角之差)1(212L若 1,则自点沿水平方向向右作长度为( 1)的直线段,其端点即为 的机翼的 上的一个相应点。在机翼的 曲线上再取几个点,按相同方法找出机翼的 曲线上另几个相应点。把这些相应点连接起来,就得到机翼的 曲线。作水平线段时,应注意方向。 小的机翼的 曲线总是在 大的 曲线的右边。同样的方法可以对有限翼展机翼和无限翼展机翼的 曲
19、线进行换算。 例 12.1 已知一飞机重 21582N,机翼面积为 20m ,翼展 l=11m,若飞机飞行速度为 280km/h,其飞行方向为水平方向,1.226kg/m 3,求:1) 系数,展弦比;环量;2)若翼型平面形状为矩形, 求诱导阻力系数.50)1(解: 展弦比 .620Al因飞机水平方向飞行, 所以升力与飞机自重力平衡,则升力系数29.01)(6.1582.30221 VLC由升力定理l,可求得环量为sml /57 26.3280诱导阻力系数0.005)1(2LDiC例 12.2 船舵为 NACA0015的对称翼型,对于无限翼展机翼,当翼弦为时,求:1) =0和 15时的升力和诱导
20、阻力各为多少?2)若翼弦不变,翼展为,在相同的几何冲角下,升力和诱导阻力又各为多少?解: 已知 , )/(1097.odCL1) =, 0时 )(0L故 因为无限翼展机翼,无下洗,故 15时 45.1097.)(dCL故 )(272145.2123 NVVblVL)当翼弦=1,翼展=2时, , , ,对称翼型 0.1775 oL radC6.15)(273.05.24.1)( 再根据展弦比换算,求出 时 =2 的升力系数。 OCDABoLoL156.152 )(71320732 23NVblV查表得 0.05,所以=0.085)1(2LDiC所以诱导阻力 (852NVRii本题 ,用有限翼展机翼公式进行计算有误差。例 12-3 一机翼弦长2,展长10,以360的速度在大气中飞行,设机翼中部的环量 20 2/s 两端为零,环量沿翼展呈椭园型分布。试求该机翼的升力系数及诱导阻力系数(1.2 ) 。解:取轴为展向,原点在机翼中点,故环量分布为1)2/(max2ly所以 22max5102)(yyl则升力)(85.1425KNdyVL升力系数 157.02ACL诱导阻力系数 .Di