NACA0012翼型的截面与升力阻力曲线图
飞机升力与阻力详解图文
飞行基础知识①升力与阻力详解图文升力是怎样产生的任何航空器都必须产生大于自身重力的升力才能升空飞行,这是航空器飞行的基本原理.前面我们提到,航空器可分为轻于空气的航空器和重于空气的航空器两大类,轻于空气的航空器如气球、飞艇等,其主要部分是一个大大的气囊,中间充以比空气密度小的气体如热空气、氢气等,这样就如同我们小时候的玩具氢气球一样,可以依靠空气的静浮力升上空中.远在一千多年以前,我们的祖先便发明了孔明灯这种借助热气升空的精巧器具,可以算得上是轻于空气的航空器的鼻祖了.然而,对于重于空气的航空器如飞机,又是靠什么力量飞上天空的呢相信大家小时候都玩过风筝或是竹蜻蜓,这两种小小的玩意构造十分简单,但却蕴含着深刻的飞行原理.飞机的机翼包括固定翼和旋翼两种,风筝的升空原理与滑翔机有一些类似,都是靠迎面气流吹动而产生向上的升力,但与固定翼的飞机有一定的差别；而旋翼机与竹蜻蜓却有着异曲同工之妙,都是靠旋翼旋转产生向上的升力.机翼是怎样产生升力的呢让我们先来做一个小小的试验：手持一张白纸的一端,由于重力的作用,白纸的另一端会自然垂下,现在我们将白纸拿到嘴前,沿着水平方向吹气,看看会发生什么样的情况.哈,白纸不但没有被吹开,垂下的一端反而飘了起来,这是什么原因呢流体力学的基本原理告诉我们,流动慢的大气压强较大,而流动快的大气压强较小,白纸上面的空气被吹动,流动较快,压强比白纸下面不动的空气小,因此将白纸托了起来.这一基本原理在足球运动中也得到了体现.大家可能都听说过足球比赛中的“香蕉球”,在发角球时,脚法好的队员可以使足球绕过球门框和守门员,直接飞入球门,由于足球的飞行路线是弯曲的,形似一只香蕉,因此叫做“香蕉球”.这股使足球偏转的神秘力量也来自于空气的压力差,由于足球在踢出后向前飞行的同时还绕自身的轴线旋转,因此在足球的两个侧面相对于空气的运动速度不同,所受到的空气的压力也不同,是空气的压力差蒙蔽了守门员.对于固定翼的飞机,当它在空气中以一定的速度飞行时,根据相对运动的原理,机翼相对于空气的运动可以看作是机翼不动,而空气气流以一定的速度流过机翼.空气的流动在日常生活中是看不见的,但低速气流的流动却与水流有较大的相似性.日常的生活经验告诉我们,当水流以一个相对稳定的流量流过河床时,在河面较宽的地方流速慢,在河面较窄的地方流速快.流过机翼的气流与河床中的流水类似,由于机翼一般是不对称的,上表面比较凸,而下表面比较平,流过机翼上表面的气流就类似于较窄地方的流水,流速较快,而流过机翼下表面的气流正好相反,类似于较宽地方的流水,流速较上表面的气流慢.根据流体力学的基本原理,流动慢的大气压强较大,而流动快的大气压强较小,这样机翼下表面的压强就比上表面的压强高,换一句话说,就是大气施加与机翼下表面的压力方向向上比施加于机翼上表面的压力方向向下大,二者的压力差便形成了飞机的升力.当飞机的机翼为对称形状,气流沿着机翼对称轴流动时,由于机翼两个表面的形状一样,因而气流速度一样,所产生的压力也一样,此时机翼不产生升力.但是当对称机翼以一定的倾斜角称为攻角或迎角在空气中运动时,就会出现与非对称机翼类似的流动现象,使得上下表面的压力不一致,从而也会产生升力.飞机的阻力凡是懂得物理知识的人都知道,飞机在飞行的过程中,机体上所受的力是平衡的.飞机的重力与飞机产生的升力平衡,而飞机的发动机的作用则是克服飞机所受的阻力,推动飞机前进,使得飞机相对于空气运动,从而产生升力.大家肯定要想,飞机发动机的功率那么大,难道飞机上所受的阻力有那么大吗的确,飞机在高速飞行的同时,会因为不同原因受到非常大的阻力.从产生阻力的不同原因来说,飞机所受的阻力可以分为摩擦阻力、压差阻力、诱导阻力、干扰阻力、激波阻力等.摩擦阻力当两个物体相互滑动的时候,在两个物体上就会产生与运动方向相反的力,阻止两个物体的运动,这就是物体之间的摩擦阻力.当飞机在空气中飞行时,飞机也会受到空气的摩擦阻力,飞机的摩擦阻力是因为空气的粘性造成的.当气流流过物体时,由于粘性,空气微团与物体表面发生摩擦,阻滞了气流的流动,这就是物体对空气的摩擦阻力,反之,空气对物体也给予了摩擦阻力.摩擦阻力是在边界层中产生的.所谓边界层就是紧贴物体表面,流速由外部流体的自由流速逐渐降低到零的那一层薄薄的空气层.边界层中气流的流动情况是不同的.一般机翼大约在最大厚度之前,边界层的气流各层不相混杂而成层地流动,这部分叫做“层流边界层”.在这之后,气流的活动转变为杂乱无章,并且出现了漩涡和横向流动,这部分叫做“紊流边界层”.从“层流边界层”转变为“紊流边界层”的那一点叫做“转捩点”. 边界层中的摩擦阻力大小与流动情况有很大关系,从大量的实践证明,对于层流流动,物体表面受到的摩擦阻力小,而紊流流动对物面的摩擦阻力大的多.在普通的机翼表面,既有层流边界层,又有紊流边界层,所以为了减小摩擦阻力,人们就千方百计地使物体表面的流动保持层流状态,例如通过在机翼表面上钻孔,吸除紊流边界层,这样就可以达到减阻的目的.另外,提高加工精度,使层流边界层尽量的长,延缓转捩点的出现,甚至抑制它的出现,也可以起到很好的效果.这些都是飞机设计中的层流机翼的概念.物体表面受到的摩擦阻力还跟物体的表面积有关系,面积越大,阻力也越大.因此在人们试图减小飞行阻力的时候,减小飞机的尾翼或者机翼的面积也是一个有效的方法.当然前提条件是保证产生足够的升力和控制力.例如使用推力矢量技术的飞机,由于有了发动机推力直接用于飞行控制,这样飞机的尾翼就可以减小或者去除,这样就可以大大的减小摩擦阻力.诱导阻力机翼同一般物体相似,也有摩擦阻力和压差阻力.对于机翼而言,这二者合称“翼型阻力”.机翼上除翼型阻力外,还有“诱导阻力”又叫“感应阻力”.这是机翼所独有的一种阻力.因为这种阻力是伴随着机翼上举力的产生而产生的.也许可以说它是为了产生举力而付出的一种代价.如果有一架飞机以某一正迎角a作水平飞行,它的机翼上面的压强将降低,而下面的压强将增高,加上空气摩擦力,于是产生了举力Y.这是气流作用到机翼上的力,根据作用和反作用定律,必然有一个反作用力即负举刀力-Y,由机翼作用到气流上,它的方向向下,所以使气流向下转折一个角度a,这一角度叫“下洗角”.随着下洗角的出现,同时出现了气流向下的速度.这一速度叫做“下洗速w”.下洗的存在还可由风洞实验观察出来.由实验可知：当飞机飞行时,下翼面压强大、上翼面压强小.由于翼展的长度是有限的,所以上下翼面的压强差使得气流从下翼面绕过两端翼尖,向上翼面流动.当气流绕流过翼尖时,在翼尖那儿不断形成旋涡.旋涡就是旋转的空气团.随着飞机向前方飞行,旋涡就从翼尖向后方流动,并产生了向下的下洗速w.下洗速在两个翼尖处最大,向中心逐渐减小,在中心处减到最小.这是因为旋涡可以诱导四周的空气随之旋转,而这又是由于空气粘性所起的作用.空气在旋转时,越靠内圈,旋转得越快,越靠外圈,旋转得越慢.因此,离翼尖越远,气流垂直向下的下洗速就越小.图示的就是某一个翼剖面上的下洗速度.它与原来相对速度v组成了合速度u .u与v的夹角就是下洗角a1.下洗角使得原来的冲角a减小了.根据举力Y原来的函义,它应与相对速度v垂直,可是气流流过机翼以后,由于下洗速w的作用,使v的方向改变,向下转折一个下洗角a1,而成为u.因此,举力Y也应当偏转一角度a1,而与u垂直成为y 1.此处下洗角很小,因而y与y1一般可看成相等.回这时飞机仍沿原来v的方向前进.y1既不同原来的速度v垂直,必然在其上有一投影为Q；.它的方向与飞机飞行方向相反,所起的作用是阻拦飞机的前进.实际上是一种阻力.这种阻力是由举力的诱导而产生的,因此叫做“诱导阻力”.它是由于气流下洗使原来的举力偏转而引起的附加阻力,并不包含在翼型阻力之内.图中机翼前面的一排小箭头表示原来的流速,后面的一排小箭头则表示流过机翼后偏转一个角度的流速.诱导阻力同机翼的平面形状,翼剖面形状,展弦比,特别是同举力有关.压差阻力“压差阻力”的产生是由于运动着的物体前后所形成的压强差所形成的.压强差所产生的阻力、就是“压差阻力”.压差阻力同物体的迎风面积、形状和在气流中的位置都有很大的关系.用刀把一个物体从当中剖开,正对着迎风吹来的气流的那块面积就叫做“迎风面积”.如果这块面积是从物体最粗的地方剖开的,这就是最大迎风面积.从经验和实验都不难证明：形状相同的物体的最大迎风面积越大,压差阻力也就越大.物体形状对压差阻力也有很大的作用.把一块圆形的平板,垂直地放在气流中.它的前后会形成很大的压差阻力.平板后面会产生大量的涡流,而造成气流分离现象.如果在圆形平板的前面加上一个圆锥体,它的迎风面积并没有改变,但形状却变了.平板前面的高压区,这时被圆锥体填满了.气流可以平滑地流过,压强不会急剧升高,显然这时平板后面仍有气流分离,低压区仍然存在,但是前后的压强差却大为减少,因而压差阻力降低到原来平板压差阻力的大约五分之一.如果在平板后面再加上一个细长的圆锥体,把充满旋涡的低压区也填满,使得物体后面只出现很少的旋涡,那么实验证明压差阻力将会进一步降低到原来平板的大约二十到二十五分之象这样前端圆纯、后面尖细,象水滴或雨点似的物体,叫做“流线形物体”,简称“流线体”.在迎风面积相同的条件下,它的压差阻力最小.这时阻力的大部分是摩擦阻力.除了物体的迎风面积和形状外,物体在气流中的位置也影响到压差阻力的大小.物体上的摩擦阻力和压差阻力合起来叫做“迎面阻力”.一个物体,究竟哪一种阻力占主要部分,这要取决于物体的形状和位置.如果是流线体,那么它的迎面阻力中主要部分是摩擦阻力.如果形状远离流线体的式样,那么压差阻力占主要部分,摩擦阻力则居次要位置,而且总的迎面阻力也较大.干扰阻力飞机上除了摩擦阻力,压差阻力和诱导阻力以外,还有一种“干扰阻力”值得我们注意,实践表明,飞机的各个部件,如机翼、机身、尾翼等,单独放在气流中所产生的阻力的总和并不等于、而是往往小于把它们组成一个整体时所产生的阻力.所谓“干扰阻力”就是飞机各部分之间由于气流相互干扰而产生的一种额外阻力.如图所示,气流流过机翼和机身的连接处,由于机翼和机身二者形状的关系,在这里形成了一个气流的通道.在A处气流通道的截面积比较大,到C点翼面最圆拱的地方,气流通道收缩到最小,随后到B处又逐渐扩大.根据流体的连续性定理和伯努利定理,C处的速度大而压强小,B处的速度小而压强大,所以在CB一段通道中,气流有从高压区B回流到低压区C 的趋势.这就形成了一股逆流.但飞机前进不断有气流沿通道向后流,遇到了后面的这股逆流就形成了气流的阻塞现象,使得气流开始分离,而产生了很多旋涡.这些旋涡表明气流的动能有了消耗,因而产生了一种额外的阻力,这一阻力是气流互相干扰而产生的,所以叫做“干扰阻力”.不但在机翼和机身之间可能产生干扰阻力,而且在机身和尾翼连接处,机翼和发动机短舱连接处,也都可能产生.从干扰阻力产生的原因来看,它显然和飞机不同部件之间的相对位置有关.如果在设计飞机时,仔细考虑它们的相对位置,使得它们压强的增加不大也不急剧,干扰阻力就可减小.另外,还可以采取在不同部件的连接处加装流线型的“整流片”的办法,使连接处圆滑过渡,尽可能减少漩涡的产生,也可减少“干扰阻力”.激波阻力飞机在空气中飞行时,前端对空气产生扰动,这个扰动以扰动波的形式以音速传播,当飞机的速度小于音速时,扰动波的传播速度大于飞机前进速度,因此它的传播方式为四面八方；而当物体以音速或超音速运动时,扰动波的传播速度等于或小于飞机前进速度,这样,后续时间的扰动就会同已有的扰动波叠加在一起,形成较强的波,空气遭到强烈的压缩、而形成了激波.空气在通过激波时,受到薄薄一层稠密空气的阻滞,使得气流速度急骤降低,由阻滞产生的热量来不及散布,于是加热了空气.加热所需的能量由消耗的动能而来.在这里,能量发生了转化--由动能变为热能.动能的消耗表示产生了一种特别的阻力.这一阻力由于随激波的形成而来,所以就叫做"波阻".从能量的观点来看,波阻就是这样产生的.从机翼上压强分布的观点来看,波阻产生的情况大致如下；根据对机翼所作的实验,在超音速飞行时,机翼上的压强分布如图所示.在亚音速飞行情况下,机翼上只有摩擦阻力、压差阻力和诱导阻力.它的压力分布如图中虚线所示.对图中两种不同的飞行情况压强分布加以比较,可以看出：在亚音速飞行情况下,最大稀薄度靠前,压强分布沿着与飞行相反的方向上的合力,不是很大,即阻力不是很大,其中包括翼型阻力和诱导阻力.可是在超音速飞行情况下,压强分布变化非常大,最大稀薄度向后远远地移动到尾部,而且向后倾斜得很厉害,同时它的绝对值也有增加.因此,如果不考虑机翼头部压强的升高,那么压强分布沿与飞行相反方向的合力,急剧增大,使得整个机翼的总阻力相应有很大的增加.这附加部分的阻力就是波阻.由于它来自机翼前后的压力差,所以波阻实际上是一种压差阻力.当然,如果飞机或机翼的任何一点上的气流速度不接过音速,是不会产生激波和波阻的.阻力对于飞机的飞行性能有很大的影响,特别是在高速飞行时,激波和波阻的产生,对飞机的飞行性能的影响更大.这是因为波阻的数值很大,能够消耗发动机一大部分动力.例如当飞行速度在音速附近时,根据计算,波阻可能消耗发动机大约全部动力的四分之三.这时阻力系数Cx急骤地增长好几倍.这就是由于飞机上出现了激波和波阻的缘故.由上面所说的看来,波阻的大小显然同激波的形状有关,而激波的形状在飞行M数不变的情况下；又主要决定于物体或飞机的形状,特别是头部的形状.按相对于飞行速度或气流速度成垂直或成偏斜的状态,有正激波和斜激波两种不同的形状.成垂直的是正激波,成偏斜的是斜激波.在飞行M数超过1时例如M等于2,如果物体的头部尖削,象矛头或刀刃似的,形成的是斜激波；如果物体的头部是方楞的或圆钝的,在物体的前面形成的则是正激波.正激波沿着上下两端逐渐倾斜,而在远处成为斜激波,最后逐渐减弱成为弱扰动的边界波.斜激波的情况也是一样的,到末端也逐渐减弱而转化为边界波.在正激波之后的一小块空间,气流穿过正激波,消耗的动能很大,总是由超音速降低到亚音速,在这里形成一个亚音速区.M数的大小也对激波的形状有影响.当M数等于1或稍大于1例如M＝时,在尖头如炮弹物体前面形成的是正激波.如果M数超过1相当多例如M＝,形成的则是斜激波.正激波的波阻要比斜激波大,因为在正激波下,空气被压缩得很厉害,激波后的空气压强和密度上升的最高,激波的强度最大,当超音速气流通过时,空气微团受到的阻滞最强烈,速度大大降低,动能消耗很大,这表明产生的波阻很大；相反的,斜激波对气流的阻滞较小,气流速度降低不多,动能的消耗也较小,因而波阻也较小.斜激波倾斜的越厉害,波阻就越小.加强图仔细看啊从机翼截面观察气体流场状态下图所示为：飞机在转弯时的受力情况.假定飞机的飞行方向是由外飞进屏幕里,即飞机是在做左转弯.此时飞行员向左侧压杆,使左侧副翼上翻、右侧副翼下翻,在左翼上产生向下的力Fa、右翼上产生向上的力Fb,此二力以机身重心为中心形成一滚动力矩,使飞机向左翻滚.而从整架飞机来考虑,机翼左翻也使总升力 F 向左翻.在竖直和水平方向上将其分解,其竖直分力F1 与飞机重力G 维持平衡,保持飞机的飞行高度；水平分力F2 提供做圆周运动所需的向心力,使飞机转弯. 同理可得,飞机在俯冲时,飞行员向前推杆使平尾上的升降舵下翻,产生向上的力抬起机尾,机头向下形成俯冲姿态；爬升时向后拉杆,升降舵上翻,产生压力压下机尾,使机头向上形成爬升姿态；蹬右踏板使方向舵右翻,产生水平向左的推力推动机尾向左,使机头向右,同理,蹬左踏板使飞机向左.综上所述,如果把操纵杆向左推再向后拉,会使飞机左侧翻时做一个爬升动作,即一个左急转.其实,再复杂的机动动作也是由这么几块操纵面完成的,也就是操纵杆前后左右推拉以及不同高度、速度的排列组合,看起来开飞机好象不那么复杂吧,不过这只是在游戏里,要换成真的,光身体素质这一项就没几个人过得了关了.小常识在电子传动技术被广泛运用于航空领域之前,飞机的操控一直是依靠机械传动的,即所有操纵面的转动都要靠飞行员的体力来完成,在完成一个高过载机动时,机翼承受的加速度往往是七八个重力加速度,甚至更高,飞行员要付出的体力的巨大是可想而知的有力回馈摇杆的玩家都能体会到.而电子传动技术则彻底把飞行员从“力气活”里解放了出来,飞行变得更轻松了,也更注重技巧了,各种高难度的机动动作也诞生了,其难度也更多地反映出飞机的机动性能,而不是飞行员的身体素质.。
飞行器翼型设计
1、翼型的定义与研究发展在飞机的各种飞行状态下，机翼是飞机承受升力的主要部件，而立尾和平尾是飞机保持安定性和操纵性的气动部件。
一般飞机都有对称面，如果平行于对称面在机翼展向任意位置切一刀，切下来的机翼剖面称作为翼剖面或翼型。
翼型是机翼和尾翼成形重要组成部分，其直接影响到飞机的气动性能和飞行品质。
通常飞机设计要求，机翼和尾翼的尽可能升力大、阻力小、并有小的零升俯仰力矩。
因此，对于不同的飞行速度，机翼的翼型形状是不同的。
对于低亚声速飞机，为了提高升力系数，翼型形状为圆头尖尾形；对于高亚声速飞机，为了提高阻力发散Ma数，采用超临界翼型，其特点是前缘丰满、上翼面平坦、后缘向下凹；对于超声速飞机，为了减小激波阻力，采用尖头、尖尾形翼型。
3、NACA翼型编号NACA四位数翼族：其中第一位数代表f，是弦长的百分数；第二位数代表p，是弦长的十分数；最后两位数代表厚度，是弦长的百分数。
例如NACA 0012是一个无弯度、厚12%的对称翼型。
有现成实验数据的NACA四位数翼族的翼型有6%、8%、9%、10%、12%、15%、18%、21%、24五位数翼族的厚度分布与四位数翼型相同。
不同的是中弧线。
具体的数码意义如下：第一位数表示弯度，但不是一个直接的几何参数，而是通过设计升力系数来表达的，这个数乘以3/2就等于设计升力系数的十倍。
第二、第三两位数是2p，以弦长的百分数来表示。
最后两位数仍是百分厚度。
例如NACA 23012这种翼型，它的设计升力系数是（2）×3/20=0.30；p=30/2,即中弧线最高点的弦向位置在15%弦长处，厚度仍为12%。
一般情况下的五位数编号意义如下有现成实验数据的五位数翼族都是230-系列的，设计升力系数都是0.30，中弧线最高点的弦向位置p都在15%弦长处，厚度有12%、15%、18%、21%、24%五种。
其它改型的五位数翼型在此就不介绍了。
1、低速翼型绕流图画低速圆头翼型在小迎角时，其绕流图画如下图示。
飞机升力与阻力详解(图文)
飞行基础知识①升力与阻力详解(图文)升力是怎样产生的任何航空器都必须产生大于自身重力的升力才能升空飞行，这是航空器飞行的基本原理。
前面我们提到，航空器可分为轻于空气的航空器和重于空气的航空器两大类，轻于空气的航空器如气球、飞艇等，其主要部分是一个大大的气囊，中间充以比空气密度小的气体（如热空气、氢气等），这样就如同我们小时候的玩具氢气球一样，可以依靠空气的静浮力升上空中。
远在一千多年以前，我们的祖先便发明了孔明灯这种借助热气升空的精巧器具，可以算得上是轻于空气的航空器的鼻祖了。
然而，对于重于空气的航空器如飞机，又是靠什么力量飞上天空的呢？相信大家小时候都玩过风筝或是竹蜻蜓，这两种小小的玩意构造十分简单，但却蕴含着深刻的飞行原理。
飞机的机翼包括固定翼和旋翼两种，风筝的升空原理与滑翔机有一些类似，都是靠迎面气流吹动而产生向上的升力，但与固定翼的飞机有一定的差别；而旋翼机与竹蜻蜓却有着异曲同工之妙，都是靠旋翼旋转产生向上的升力。
机翼是怎样产生升力的呢？让我们先来做一个小小的试验：手持一张白纸的一端，由于重力的作用，白纸的另一端会自然垂下，现在我们将白纸拿到嘴前，沿着水平方向吹气，看看会发生什么样的情况。
哈，白纸不但没有被吹开，垂下的一端反而飘了起来，这是什么原因呢？流体力学的基本原理告诉我们，流动慢的大气压强较大，而流动快的大气压强较小，白纸上面的空气被吹动，流动较快，压强比白纸下面不动的空气小，因此将白纸托了起来。
这一基本原理在足球运动中也得到了体现。
大家可能都听说过足球比赛中的“香蕉球”，在发角球时，脚法好的队员可以使足球绕过球门框和守门员，直接飞入球门，由于足球的飞行路线是弯曲的，形似一只香蕉，因此叫做“香蕉球”。
这股使足球偏转的神秘力量也来自于空气的压力差，由于足球在踢出后向前飞行的同时还绕自身的轴线旋转，因此在足球的两个侧面相对于空气的运动速度不同，所受到的空气的压力也不同，是空气的压力差蒙蔽了守门员。
(完整版)NACA0012翼型俯仰振荡实例进行讲解动网格
以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程；首先需要声明的是，这个例子也是来源于网络，原作者不详，在此向他表示感谢。
1、问题描述本例是想对作简谐振荡运动的NACA0012翼型的气动特性（升力系数，阻力系数和力矩系数）进行数值计算，来流速度为V，攻角的变化规律为：Alpha（t）=A/2*sin(omega*t)其中，A=10度，omega=10*pi 弧度/秒。
2、该例需要使用动网格来实现，首先需要编写刚体运动UDF实现翼型的俯仰运动，由于在FLUENT的UDF中只能指定速度，角速度；所以，需要将攻角对时间求导，得到转动角速度的规律：D（alpha）/dt=A*omega/2*cos(omega*t)编写的UDF在附件中。
3、由于本例只是为了讲述动网格的实现，至于其他方面的设置及分析就不再讨论；这里详细讲述下动网格的建立以及动网格的预览的结果。
步骤如下：1）将mesh文件读入到FLUENT中，Grid：check，scale…，Smooth/Swap…；Display Grid；2）定义求解器为，Define：pressure-based，2D，unsteady，Implicit，Green-Gauss Node Based（因本例使用的是三角形单元）.3）编译UDF，Define->User-Defined->Functions->Complied…此时打开了Complied UDFs的窗口，Add…在选择UDF的对话框中找到NACA0012DM文件夹中的airfoil.c文件，选中，ok；此时返回到Complied UDFs的窗口点击Build，FLUENT开始进行编译，可以在FLUENT窗口看到编译的一些过程提示；等编译完成，点击Load；就将已经编译好的UDF加载到FLUENT中了。
4）定义动网格参数，Define->Dynamic Mesh（选勾，激活动网格模型）->Parameters…此时打开了Dynamic Mesh Parameters 窗口，在Models中只选取Dynamic Mesh，本例的网格类型为三角形单元，要实现的运动为小幅度的转动，因此选用的动网格更新方法为Smoothing+Remeshing；开始依次对这两种更新方法进行参数设定：Smoothing中的参数设定：Spring Constant Factor（弹簧倔强系数），该值设定为一个较小的值，在0.01到0.1之间，本例选取0.08；Boundary Node Relaxation（边界节点松弛），设定为0.5；Convergence Tolerance（收敛判据），保持默认的0.001；Number of Iterations（迭代次数），保持默认的20；Remeshing中的参数设定：为了得到较好的网格更新，本例在使用局部网格重新划分方法时，使用尺寸函数，也就是Remeshing+Must Improve Skewness+Size Function的策略。
NACA0012翼型的截面与升力阻力曲线图教学文案
NACA0012翼型的截面与升力阻力曲线图NACA0012翼型是垂直轴风力机的优选翼型，这里根据美国网站提供的技术资料绘制了翼型截面图与升力曲线图，还有根据技术书籍描绘的升力阻力曲线图，供大家参考。
NACA0012翼型的截面图由于NACA0012是对称翼型，在下图左侧数据表中仅列出了单边的数据，表中c是弦长（弦长为1.00）；x是弦长坐标（单位是x/c）；y是对应x位置的翼面与弦的距离（单位是y/c）。
图1 NACA0012翼型数据与截面图NACA0012翼型的升力曲线图图2是根据美国的技术资料数据绘制的NACA0012翼型的升力曲线图，在这张图中有多根升力曲线，显示了当雷诺数不同时翼型的最大升力系数与失速攻角都有较大的变化。
图2 NACA0012翼型升力系数曲线图NACA0012翼型的大攻角升力、阻力曲线图图3是根据“21世纪能源与动力工程类创新型应用人才培养规划教材“风力机空气动力学””一书图7.29描绘的NACA0012翼型的大攻角升力、阻力曲线图，供大家参考。
图3 NACA0012翼型大攻角升力、阻力曲线图该图显示了翼型攻角从-5度到180度的升力与阻力系数的变化，攻角在0度至10度升力系数随攻角增大而增大，阻力系数很小；超过12度时升力系数下降，阻力系数上升；攻角到40度后升力与阻力系数先是相同，然后阻力系数继续上升，升力系数下降。
NACA0012翼型的升/阻比变化曲线在图3中显示了翼型的升力与阻力随攻角的变化曲线，在实际应用中更关心的是在正常工作时的阻力情况，图4是根据“风能技术”（美Tony Burton，武鑫译）一书图3.100描绘，图中显示了NACA0012翼型在不同攻角时的升力/阻力系数比变化曲线。
图3 NACA0012翼型升/阻比变化曲线可见在失速前有最大的升力/阻力系数比值，升力约为阻力的50多倍。
当然这是一个光滑的翼型在较高雷诺数时的状态，多数情况会比该比值低一些。
附:叶片雷诺数计算示例从图2看到翼型的升力曲线受雷诺数影响较大，下面给出了叶片雷诺数的简单计算方法。
NACA 23012型翼型升阻力实验.doc
NACA 23012型翼型升阻力实验.doc摘要：本实验使用风洞测试，以NACA 23012型翼型为对象，测量不同攻角下的升力系数和阻力系数，并计算出该翼型的升阻比。
通过实验数据处理及理论计算，可以探讨影响翼型升阻比的因素以及优化翼型设计的方法。
关键词：实验原理：在飞行中，飞行器所受到的空气阻力分为两部分，一部分是法向于速度方向的，称为升力，另一部分则是平行于速度方向的，称为阻力。
升阻比是翼型的重要性能参数之一，它反映了翼型的升力水平与阻力水平的关系。
升阻比越高，则翼型的升力效率越高，相对阻力越小，这对于提高飞机的性能和经济性都是十分重要的。
在风洞实验中，我们可以通过测量翼型在不同攻角下的升力系数和阻力系数，来计算出该翼型的升阻比。
其中，升力系数和阻力系数的定义如下：$$C_L=\frac{L}{\frac{1}{2}\rho v^2S}$$其中，$C_L$为升力系数，$C_D$为阻力系数，$L$为升力，$D$为阻力，$\rho$为空气密度，$v$为风速，$S$为翼面积。
理论上，翼型的升阻比可以用L/D来表示。
但是由于实验中升力和阻力很难同时测量，因此常常采用C_L/C_D来近似表示翼型的升阻比。
实验过程：1.准备工作（1）校准称重仪器，确保称重的准确性。
（2）打开风洞，调整风速为预定值，并记录下风速的值。
（3）根据实验要求选择相应的翼型，并测量其翼展、翼弦等参数，计算出翼面积。
2.实验操作（1）将所选翼型安装在风洞的测试段上，并调整翼型的安装角度，使其与气流的流向平行。
（2）通过调整风洞流量调节阀，使风洞风速达到预定值。
（3）调整支架高度，使翼型处于所需要的攻角。
（4）在稳定的状态下，记录下此时风洞中的压力数据，包括总压力、静压力等。
（5）逐步调节支架高度，改变攻角的大小，并记录下不同攻角下的压力数据。
（6）将实验数据处理，并计算出翼型在不同攻角下的升力系数、阻力系数、升阻比等参数。
实验结果：通过实验，我们得到了NACA 23012型翼型在不同攻角下的升力系数、阻力系数的数据，计算出了该翼型在不同攻角下的升阻比。
基于CFD的机翼颤振分析
freeplay position and friction in the freeplay on aeroelastic response is analyzed. They have a
Key words: aeroelasticity ,stall flutter ,UDF，Fluent ,dynamic stall,freeplay nonlinearity
ii
南京航空航天大学硕士学位论文
图清单
图 1.1 气动弹性力三角形 .................................................................................................. 1 图 2.1CFD 流程图 .............................................................................................................. 9 图 2.2 基于弹簧光滑节点开始状况 ................................................................................ 17 图 2.3 基于弹簧光滑节点结束状况 ................................................................................ 17 图 2.4 二维网格数据结构示意图 .................................................................................... 19 图 2.5 三维网格数据结构示意图 ................................................................................... 20 图 3.1 第一套网格 ............................................................................................................ 23 图 3.2 第二套网格 ............................................................................................................ 23 图 3.3 第一套网格升力系数曲线 .................................................................................... 24 图 3.4 第二套网格升力系数曲线 .................................................................................... 24 图 3.5 阻力系数曲线比较 ................................................................................................ 24 图 3.6 失速机翼周围的流场速度分布 ............................................................................ 24 图 3.7 α 0 = 5° 时升力系数迟滞曲线和力矩系数迟滞曲线 ............................................. 25 图 3.8 α 0 = 10° 时升力系数迟滞曲线和力矩系数迟滞曲线 ........................................... 26 图 3.9 α 0 = 12° 时升力系数迟滞曲线和力矩系数迟滞曲线 ........................................... 26 图 3.10 α 0 = 15° 时升力系数迟滞曲线和力矩系数迟滞曲线 ......................................... 26 图 3.11 深度失速时（ α 0 = 12° ）机翼周围流场的速度分布 ........................................ 28 图 3.12 α1 = 2° 时升力系数迟滞曲线和力矩系数迟滞曲线 ........................................... 29 图 3.13 α1 = 5° 时升力系数迟滞曲线和力矩系数迟滞曲线 ........................................... 29 图 3.14 α1 = 10° 时升力系数迟滞曲线和力矩系数迟滞曲线.......................................... 29 图 3.15 α1 = 15° 时升力系数迟滞曲线和力矩系数迟滞曲线.......................................... 30 图 3.16 k = 0.05 ，不同雷诺数下的非定常特性比较..................................................... 30 图 3.17 k = 0.1 ，不同雷诺数下的非定常特性比较 ....................................................... 31 图 3.18 k = 0.15 ，不同雷诺数下的非定常特性比较..................................................... 31 图 3.19 k = 0.2 ，不同雷诺数下的非定常特性比较....................................................... 31 图 3.20 k = 0.4 ，不同雷诺数下的非定常特性比较....................................................... 32 图 4.2 具有 2 个自由度的翼型示意图 ............................................................................ 36 图 4.3 复合材料夹层板结构机翼模型 ............................................................................ 38 图 4.4V＝40m/s，二维翼型的颤振响应 ........................................................................ 39 图 4.5V＝46.75m/s，二维翼型的颤振响应 ................................................................... 39
NACA0012翼型的截面与升力阻力曲线图
NACA0012翼型的截面与升力阻力曲线图NACA0012翼型是垂直轴风力机的优选翼型，这里根据美国网站提供的技术资料绘制了翼型截面图与升力曲线图，还有根据技术书籍描绘的升力阻力曲线图，供大家参考。
NACA0012翼型的截面图由于NACA0012是对称翼型，在下图左侧数据表中仅列出了单边的数据，表中c 是弦长（弦长为1.00）；x是弦长坐标（单位是x/c）；y是对应x位置的翼面与弦的距离（单位是y/c）。
图1 NACA0012翼型数据与截面图NACA0012翼型的升力曲线图图2是根据美国的技术资料数据绘制的NACA0012翼型的升力曲线图，在这张图中有多根升力曲线，显示了当雷诺数不同时翼型的最大升力系数与失速攻角都有较大的变化。
图2 NACA0012翼型升力系数曲线图NACA0012翼型的大攻角升力、阻力曲线图图3是根据“21世纪能源与动力工程类创新型应用人才培养规划教材“风力机空气动力学”” 一书图7.29描绘的NACA0012翼型的大攻角升力、阻力曲线图，供大家参考。
图3 NACA0012翼型全攻角升力、阻力曲线图该图显示了翼型攻角从-5度到180度的升力与阻力系数的变化，攻角在0度至10度升力系数随攻角增大而增大，阻力系数很小；超过12度时升力系数下降，阻力系数上升；攻角到40度后升力与阻力系数先是相同，然后阻力系数继续上升，升力系数下降。
NACA0012翼型的升/阻比变化曲线在图4中显示了翼型的升力与阻力随攻角的变化曲线，在实际应用中更关心的是在正常工作时的阻力情况，图4是根据“风能技术”（美Tony Burton，武鑫译）一书图3.100描绘，图中显示了NACA0012翼型在不同攻角时的升力/阻力系数比变化曲线。
图4 NACA0012翼型升/阻比变化曲线可见在失速前有最大的升力/阻力系数比值，升力约为阻力的50多倍。
当然这是一个光滑的翼型在较高雷诺数时的状态，多数情况会比该比值低一些。
NACA0012翼型Gurney襟翼增升特性及其机理实验研究
——兰些堕：：！二堂堡！翌！兰垒墨型鱼！翌型堡堡望盐堑苎！！垒些壅鐾堕壅ａ速度矢量陶ｆ无襟翼，ｔ－－Ｏ）一＝》ｏ４ｘ／ｃｏ０１０２０．３ｏ４ｘ／ｃｃ速度矢量图（６‰，ｔ＝Ｏ）ｅ速度矢量图（６‰，仁Ｏ＋４Ｇｔ）ｂ速度矢量图ｆ无谗现，ｉ－，０＋２Ａｄ速度矢量图（６‰，ｗＯ＋２Ｇｔ）Ｏ￡速度矢量图（６‰，ｔ＝Ｏ＋６Ａｔ）ｇ流线圈（６‰，ｔ－－ｏ）ｈ．流线图（６‰，ｔ＝Ｏ＋２Ａｔ）ｉ染色液流动显示国ｆ６％ｃ）图１８ＮＡＣＡ００１２翼型加装Ｇｕｒｎｅｙ襟冀后的瞬时图（Ⅱ：２．５。
）０４ｘ／ｃｏ＾彩ｃ“●卧●上李亚拒、王晋军：ＮＡＣＡ００１２翼型Ｇｕｍｅｙ襟翼增升特性及其机理实验研究１６场、流线和染色液流动显示图，从图ａ、ｂ可知，不加装Ｇｕｍｅｙ襟翼时，虽然尾流也有上Ｆ波动，但并不明显，也没有明显的涡结构存在，此时流场中未见明显分离：而在图ｃ．ｆ中则明显示出，加装６％ｃＧｕｍｅｙ襟翼后，尾流包含结构非常明显的交替脱落的卡门涡街，周期性很强，位于Ｇｕｍｅｙ襟翼的正后方。
考察速度大小可以发现，翼型尾缘处上下翼面附近速度由于此交替脱落的涡街的存在而变得明显增大，且流经上翼面的流体速度方向明显下偏。
图ｇ、ｈ分别为对应图ｃ、ｄ的瞬时流线图，可以看出集中涡的位置上下变化，表明涡脱落是上下摆动的，同时，两图中涡的旋向分别为顺时针和逆时针。
从图ｉ给出的染色液瞬时图也可以看出，在Ｇｕｒｎｅｙ襟翼的最上方和最下方交替有涡脱落且呈现周期性。
圈１９给出的是ＮＡＣＡ００１２翼型在攻角５。
时的瞬时流线图，分别为不加襟翼和加装６％ｃＧｕｍｅｙ襟翼。
不加襟翼时，流体在流经上翼面尾缘处时明显发生分离，不再能够沿着流体表面流动，由此造成流经箍个翼型的流体在经过翼型尾缘后明显向上偏转，尾流区内火部分处于分离状态；虽然尾流也有摆动，但并不存在明显的旋涡结构。
加装６％ｃＧｕｍｅｙ襟翼后，上翼面尾缘处的分离明显减弱，甚至消失，尾流中一童存在两个反向旋转的集中涡，并呈现交替脱落状态。
NACA0012翼型的截面与升力阻力曲线图
NACA0012翼型的截面与升力阻力曲线图NACA0012翼型是垂直轴风力机的优选翼型，这里根据美国网站提供的技术资料绘制了翼型截面图与升力曲线图，还有根据技术书籍描绘的升力阻力曲线图，供大家参考。
NACA0012翼型的截面图由于NACA0012是对称翼型，在下图左侧数据表中仅列出了单边的数据，表中c 是弦长（弦长为1.00）；x是弦长坐标（单位是x/c）；y是对应x位置的翼面与弦的距离（单位是y/c）。
图1 NACA0012翼型数据与截面图NACA0012翼型的升力曲线图图2是根据美国的技术资料数据绘制的NACA0012翼型的升力曲线图，在这张图中有多根升力曲线，显示了当雷诺数不同时翼型的最大升力系数与失速攻角都有较大的变化。
图2 NACA0012翼型升力系数曲线图NACA0012翼型的大攻角升力、阻力曲线图图3是根据“21世纪能源与动力工程类创新型应用人才培养规划教材“风力机空气动力学”” 一书图7.29描绘的NACA0012翼型的大攻角升力、阻力曲线图，供大家参考。
图3 NACA0012翼型全攻角升力、阻力曲线图该图显示了翼型攻角从-5度到180度的升力与阻力系数的变化，攻角在0度至10度升力系数随攻角增大而增大，阻力系数很小；超过12度时升力系数下降，阻力系数上升；攻角到40度后升力与阻力系数先是相同，然后阻力系数继续上升，升力系数下降。
NACA0012翼型的升/阻比变化曲线在图4中显示了翼型的升力与阻力随攻角的变化曲线，在实际应用中更关心的是在正常工作时的阻力情况，图4是根据“风能技术”（美Tony Burton，武鑫译）一书图3.100描绘，图中显示了NACA0012翼型在不同攻角时的升力/阻力系数比变化曲线。
图4 NACA0012翼型升/阻比变化曲线可见在失速前有最大的升力/阻力系数比值，升力约为阻力的50多倍。
当然这是一个光滑的翼型在较高雷诺数时的状态，多数情况会比该比值低一些。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
NACA0012翼型的截面与升力阻力曲线图NACA0012翼型是垂直轴风力机的优选翼型，这里根据美国网站提供的技术资料绘制了翼型截面图与
升力曲线图，还有根据技术书籍描绘的升力阻力曲线图，供大家参考。
NACA0012翼型的截面图
由于NACA0012是对称翼型，在下图左侧数据表中仅列出了单边的数据，表中c是弦长（弦长为1.00）；x是弦长坐标（单位是x/c）；y是对应x位置的翼面与弦的距离（单位是y/c）。
图1 NACA0012翼型数据与截面图
NACA0012翼型的升力曲线图
图2是根据美国的技术资料数据绘制的NACA0012翼型的升力曲线图，在这张图中有多根升力曲线，显示了当雷诺数不同时翼型的最大升力系数与失速攻角都有较大的变化。
图2 NACA0012翼型升力系数曲线图
NACA0012翼型的大攻角升力、阻力曲线图
图3是根据“21世纪能源与动力工程类创新型应用人才培养规划教材“风力机空气动力学””一书图7.29描绘的NACA0012翼型的大攻角升力、阻力曲线图，供大家参考。
图3 NACA0012翼型大攻角升力、阻力曲线图
该图显示了翼型攻角从-5度到180度的升力与阻力系数的变化，攻角在0度至10度升力系数随攻角
增大而增大，阻力系数很小；超过12度时升力系数下降，阻力系数上升；攻角到40度后升力与阻力系
数先是相同，然后阻力系数继续上升，升力系数下降。
NACA0012翼型的升/阻比变化曲线
在图3中显示了翼型的升力与阻力随攻角的变化曲线，在实际应用中更关心的是在正常工作时的阻力情况，图4是根据“风能技术”（美Tony Burton，武鑫译）一书图3.100描绘，图中显示了NACA0012翼型在不同攻角时的升力/阻力系数比变化曲线。
图3 NACA0012翼型升/阻比变化曲线
可见在失速前有最大的升力/阻力系数比值，升力约为阻力的50多倍。
当然这是一个光滑的翼型在较
高雷诺数时的状态，多数情况会比该比值低一些。
附:叶片雷诺数计算示例
从图2看到翼型的升力曲线受雷诺数影响较大，下面给出了叶片雷诺数的简单计算方法。
叶片的雷诺数可直接用公式Re=（ρ/μ）（vl），
将有关数据代入计算即可。
ρ与μ随气温气压变化较大，但在固定的环境里气压变化较小、气温变化较大，我们选在标准大气压下，气温为0度、10度、20度、30度四种情况来计算。
通过查阅相关手册，计算ρ/μ的数值：算得：0度时ρ/μ=75187 10度时ρ/μ=70077
20度时ρ/μ=66186 30度时ρ/μ=48193
如果叶片宽度为1m，叶片与空气的相对速度为30m，代入Re=ρ/μ（vl）计算得：
0度时Re=2255610 10度时Re=2102310
20度时Re=1985580 30度时Re=1445790
如果叶片宽度为0.2m，叶片与空气的相对速度为10m，代入Re=ρ/μ（vl）计算得：
0度时Re=150374 10度时Re=140154
20度时Re=132372 30度时Re=96386
前面的翼型的升力曲线图未提供雷诺数在160000以下的数据，可能原数据主要围绕飞行器使用，没有较小翼型低速运行时的数据。
上述0.2m宽翼型、10m/s速度的雷诺数均小于160000，在该图上就没有了。
由于许多翼型的雷诺数在低于60000后，升力明显下降，建议叶片宽度应在0.1m以上，较宽的翼型才能有较高的运行效率，翼型的雷诺数能在1000000以上是比较理想的状态。
以上计算结果是表面光滑的叶片，如果叶片表面粗糙或有粘覆物雷诺数与计算结果会有差别。