飞行器空气动力计算.docx
- 文档编号:2050273
- 上传时间:2022-10-26
- 格式:DOCX
- 页数:17
- 大小:431.69KB
飞行器空气动力计算.docx
《飞行器空气动力计算.docx》由会员分享,可在线阅读,更多相关《飞行器空气动力计算.docx(17页珍藏版)》请在冰豆网上搜索。
飞行器空气动力计算
第一章飞行器基本知识
1.1飞行器几何参数
飞行器通常由机翼、机身、尾翼以及动力装置等部件组成。
对于气动正问题及气动分析而言,已知飞行器几何外形,求其气动参数。
要解决这一问题首先要计算出飞行器各部件及组合体的几何参数。
当机翼和机身组合成一体时,机翼中间一部分面积为机身所遮蔽。
它外露在气流中的部分两边合起来,所构成的机翼为外露翼,由下标“”表示
在组合体中把外露翼根部的前后缘向机身内延长并交于机身纵对称面,这样的机翼成为毛机翼。
第2章机翼的气动特性分析
2.1机翼几何参数
2.1.1翼型的几何参数
翼型的前缘点与后缘点的连线称为弦线。
他们之间的距离称为弦长,用符号b表示,是翼型的特征长度。
可以想象翼型是由厚度分布和中弧线分布叠加而成的,对于中等厚度和弯度的翼型,上下翼面方程可以写成
(2—1)
式中的正号用于翼型上表面,负号用于下表面。
,分别为纵、横向无量纲坐标。
相对厚度和相对弯度,。
最大厚度位置和最大弯度位置分别用和或用无量纲量和表示。
翼型前缘的内切圆半径叫做前缘半径,用表示,后缘角τ是翼型上表面和下表面在后缘处的夹角。
2.1.2机翼的几何参数
1.机翼平面形状:
根梢比、展弦比和后掠角
机翼面积S是指机翼在xOz平面上的投影面积,即
(2—2)
式中,b(z)为当地弦长。
几何平均弦长和平均气动弦长分别定义为
(2—3)
(2—4)
显然,是面积和展长都与原机翼相等的当量矩形翼的弦长;而是半翼面心所在的展向位置的弦长,通常取作为纵向力矩的参考长度。
除了上述几何参数外,还有根梢比、梢根比和展弦比。
根梢比和梢根比定义为
,=1/(2—5)
展弦比是机翼展向伸长程度的量度,定义为
(2—6)
梯形后掠翼前缘与z轴的夹角叫做前缘后掠角,用表示,常用的还有1/4弦线、1/2弦线和后缘线的后掠角,分别用,和表示。
如图2—2所示。
2.2翼型的低速气动特性
2.2.1翼型的升力和力矩特性
黏性对失速前翼型升力特性的影响是可以忽略的。
此外,只要翼型相对厚度和相对弯度都很小,并且翼型的迎角也不大,那么翼型表面上压强的合力大小和方向就只受到厚度分布的轻微影响。
对于这样的微弯薄翼,翼型的升力和力矩特性可以用气流绕它的中弧线流动而求得,可以用薄翼理论来计算。
压强和载荷
根据伯努利方程,流动中某点的压强系数与该点的速度有如下关系:
(2—7a)
式中,v=i+j,和为扰动速度,为来流速度。
对于小扰动情况,即,略去二阶小量后式(2—7a)简化为
(2—7b)
弦向点x处下翼面与上翼面的压强与之差为载荷,用符号表示,为
(2—8)
式中为载荷系数,。
对于薄翼,整个翼型是由厚度分布和中弧线叠加而成的,图2—3。
在小迎角情况下,根据线化方程和边界条件,翼型的压强系数可以表示成由厚度和弯度(包括迎角)贡献的叠加,即
式中,为当迎角和弯度时,由厚度产生的压强系数;为中弧线和迎角产生的压强系数。
升力和力矩特性
薄翼理论的结果。
翼型的升力系数和绕翼型前缘的力矩系数为
(2—9)
(2—10)
式中,规定力矩使翼型前缘抬头为正,载荷与环境密度的关系为
(2—11)
由薄翼理论有
(2—12)
由式(2—9)至式(2—12)得
(2—13)
(2—14a)
用升力系数表示的力矩系数可写成
(2—14b)
式(2—12)至式(2—14)中的多项式系数与中弧线方程的关系为
(2—15)
1.翼型的升力特性
将式(2—15)的系数代入式(2—13),改写为
(2—16)
式中,为零升迎角,它代表零升力线与弦线的夹角,图2—4。
它仅与中弧线形状有关。
此式说明翼型的升力系数随几何迎角成线性变化。
将对求导,得薄翼理论的升力线斜率
(2—17)
2.翼型的力矩特性
对于给定的翼型,式(2—14b)等号右边的第二项为常量,故与成线性关系,可将式(2—14b)改为
(2—18)
式中,是零升力矩系数,它与翼型的升力或迎角无关,仅是翼型弯度分布的函数;是力矩系数对升力系数的导数。
如果对翼型的1/4弦点取力矩,并利用式(2—18),可得
(2—19)
显然,对于薄翼理论而言,1/4弦点力矩系数与升力系数(或迎角)无关,它就等于零升力矩系数。
在翼型上有两个重要的特性点,一个是焦点(或称气动中心),另一个是压力中心。
1)翼型上存在这样一个点,该点的力矩系数与升力系数无关,这一点称为翼型的焦点。
焦点的弦向相对量用表示。
既然绕焦点的力矩与升力系数无关,故它是升力增量的作用点。
因此,对于前缘力矩系数又可写成
(2—20)
将式(2—20)与式(2—18)的第一式相比较,可得基于薄翼理论的焦点位置
(2—21)
2)翼型的升力作用线与弦线的交点称为压力中心,压力中心的弦向相对位置用表示。
根据上述定义,将前缘力矩系数除以升力系数,可得
(2—22)
从方程(2—15)知,和都与迎角无关,至取决于中弧线形状,故压力中心将随变化。
对于对称翼型(弯度分布),==0,薄翼理论压力中心与焦点重合,即。
【例2—1】某一翼型的弯度分布,试求该翼型的升力和力矩特性。
解该翼型的弯度分布沿的变化率为
由式(2—15)得
于是根据式(2—16)~式(2—22)有
由最后一个式子可以看出,在零迎角下该翼型的压力中心,当迎角或增大时,它将移向焦点。
2.4.2超声速薄翼型的线性化位流理论
超声速线化速势方程为
(2—31)
式中,。
流动方程式(2—31)的通解为
(2—32)
式中,f和g是自变量为和的任意函数。
可以看到
常量,常量(2—33a)
的两族直线对于x轴的倾角分别为和,因此它们正好代表来流的两族马赫波,如下图所示。
在翼型的上半平面流场中,函数代表翼型上表面所发出的扰动沿马赫线常量向下游传播到流场点(x,y)所产生的扰动速度位;而代表翼型下表面发出的扰动沿马赫线常量向下游传播到流场点(x,y)所产生的扰动速度位。
在超声速流场中,有意义的解是往下游传播的,而且,受到扰动的区域也只局限于前后缘马赫波之间。
所以对上、下半平面的流场的小扰动速度位分别是
,(2—34a)
可见,沿着翼型上表面的马赫波(常量)或沿下表面的马赫波(常量)为常量,而且,流场上沿着马赫波的两扰动速度分量和以及其他流动参数也都是常量。
函数和科根据翼型绕流边界条件确定。
设翼型的上表面方程为,由线化边界条件有
(2—35)
对于上表面令,则有
(2—36)
线化压强系数公式为
(2—37)
联立式(2—35)、(2—36)和(2—37)。
得
(2—38)
类似地,如果翼型下表面方程为,则
(2—39)
根据线化理论,翼型表面上任一点的压强系数与该点翼面对于来流方向的斜率成正比。
由于物面上任一点相对于来流方向倾角都很小,所以该点物面斜率科表示为
这样,式(2—38)和(2—39)科合并成
(2—40)
式(2—38)至式(2—40)就是线化、二维超声速的基本关系。
式(2—40)表明,物面上。
任一点的压强系数与该点相对于来流的倾角成正比。
相对于来流为压缩的物面倾角取正直,为膨胀的物面倾角则取负值。
2.4.3翼型的超声速气动特性
对于薄翼,来流迎角很小,可认为翼型的整个气动力是由厚度、弯度和迎角产生的气动力的代数和。
如图2—19所示,将x轴沿着翼弦方向,y轴与x轴垂直,则上、下翼面相对于来流的倾角可表示成
(2—41)
式中,分别代表厚度()翼型和弯度()翼型表面上某一点的倾角,即
(2—42)
将式(2—41)代入式(2—40),得任意形状翼面上下表面压强分布
(2—43)
由式(2—34)可见,对于任意形状的翼型,它的表面压强系数可认为是厚度分布、弯度分布和来流迎角所产生的压强系数的代数和。
1.升力特性
从图2—19可见,作用在翼型微元上的升力为
式中,和分别为翼型上下表面微元长度。
翼型微元升力系数为
(2—44)
将式(2—43)代入式(2—44),得
(2—45)
由于在翼型前后缘,因此有
(2—46)
将式(2—45)对x从零积分到b,并应用是(2—46),得
(2—47)
从式(2—47)可见,在超声速线化理论中,薄翼型的升力系数与厚度和弯度分布无关,升力系数与迎角成正比。
2.波阻力特性
从图2—19可见,翼型微元上的阻力位
因为对于薄翼,所以上式可写为
(2—48)
利用式(2—46),可以看出式(2—48)中的积分为零,所以翼型的波阻力系数可表示为
(2—49)
对于给定翼型都是知道的,令
(2—50)
式中,分别为翼型的相对厚度和相对弯度;分别是与翼型厚度分布和弯度分布有关的几何常数。
利用式(2—50),式(2—49)可改写为
(2—51)
式(2—51)等号右边第二项与迎角无关,仅与翼型厚度分布和弯度分布有关,对于弯度为零的翼型。
零升波阻力系数为
(2—53)
表2—2给出了几种超音速对称翼型的值,按式(2—50)计算,或由经验给出。
菱形翼型的波阻系数是最小的。
翼剖面
简图
正弦形
/8
四角形
1/4(1-)
六角形
1/(1-a/b)
菱形
1
圆弧或抛物线形
4/3
亚声速翼剖面
2.5—4
3.力矩特性
如果忽略压强弦向分量和其他高阶小量对俯仰力矩的贡献,那么对翼型前缘点的微元俯仰力矩可表示为
将式(2—43)代入上式,得翼型俯仰力矩系数
(2—54)
由式(2—47)和式(2—54)可见,由于厚度产生的压强对翼弦是上下对称的,所以厚度对升力和力矩都无贡献。
力矩系数是迎角和弯度作用的代数和,它们分别是
(2—55)
(2—56)
对于给定的翼型,式(2—56)中只是与翼型的弯度分布函数有关的参数,它的表达式为
(2—57a)
或
(2—57b)
压力中心距前缘的相对距离为
(2—58)
翼型焦点距前缘的相对距离为
(2—59)
在低速时,翼型的焦点;而在超声速时。
由此可见,从低速到超声速焦点位置显著后移了,这是研究飞行器的稳定性与操作性时必须注意的一个问题。
由式(2—58)和(2—59)可见,由线化理论给出的压力中心位置和焦点位置仍然不随马赫数变化而变化,这与低、亚声速是相同的。
【例2—2】有一双凸面的翼型如图2—20所示,该翼型的上下表面方程分别为
和
试用线化理论计算该翼型在时的升阻和力矩系数以及压力中心位置。
解设翼型的迎角为,并由给定的翼型表面方程,根据式(2—41)和式(2—42),上下表面任一点切线与来流方向的倾角分别为
(1)升力系数。
升力系数仅与迎角有关,故
(2)波阻系数。
该翼型的厚度分布函数,弯度分布函数,所以零升波阻系数为
该翼型的总波阻系数
(3)力矩系数。
与该翼型的弯度分布函数相对应的任一点倾角为
将上式代入式(2—56),得到弯度产生的绕前缘力矩系数
该翼型总力矩系数为
压力中心位置
(4)升阻比。
升力Y与阻力X之比。
在位流理论中,翼型的阻力就等于波阻力,故
第三章机身的气动特性分析
3.1机身几何参数
机身的几何参数列出如下:
(1)为机身总长度;为机身头部长度;为机身尾部长度;为机身最大宽度;为机身最大高度;
(2)为旋成体机身最大
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 飞行器 空气 动力 计算