飞行器空气动力计算

1、第一章 飞行器基本知识1.1飞行器几何参数飞行器通常由机翼、机身、尾翼以及动力装置等部件组成。对于气动正问题及气动分析而言，已知飞行器几何外形，求其气动参数。要解决这一问题首先要计算出飞行器各部件及组合体的几何参数。当机翼和机身组合成一体时，机翼中间一部分面积为机身所遮蔽。它外露在气流中的部分两边合起来，所构成的机翼为外露翼，由下标“”表示在组合体中把外露翼根部的前后缘向机身内延长并交于机身纵对称面，这样的机翼成为毛机翼。第2章 机翼的气动特性分析2.1机翼几何参数2.1.1 翼型的几何参数 翼型的前缘点与后缘点的连线称为弦线。他们之间的距离称为弦长，用符号b表示，是翼型的特征长度。可以想象翼

2、型是由厚度分布和中弧线分布叠加而成的，对于中等厚度和弯度的翼型，上下翼面方程可以写成 （21）式中的正号用于翼型上表面，负号用于下表面。，分别为纵、横向无量纲坐标。相对厚度和相对弯度，。最大厚度位置和最大弯度位置分别用和或用无量纲量和表示。翼型前缘的内切圆半径叫做前缘半径，用表示，后缘角是翼型上表面和下表面在后缘处的夹角。2.1.2 机翼的几何参数1.机翼平面形状：根梢比、展弦比和后掠角机翼面积S是指机翼在xOz平面上的投影面积，即 （22）式中，b（z）为当地弦长。几何平均弦长和平均气动弦长分别定义为 （23） （24）显然，是面积和展长都与原机翼相等的当量矩形翼的弦长；而是半翼面心所在的展

3、向位置的弦长，通常取作为纵向力矩的参考长度。除了上述几何参数外，还有根梢比、梢根比和展弦比。根梢比和梢根比定义为 ，=1/ （25）展弦比是机翼展向伸长程度的量度，定义为 （26） 梯形后掠翼前缘与z轴的夹角叫做前缘后掠角，用表示，常用的还有1/4弦线、1/2弦线和后缘线的后掠角，分别用，和表示。如图22所示。2.2 翼型的低速气动特性2.2.1 翼型的升力和力矩特性 黏性对失速前翼型升力特性的影响是可以忽略的。此外，只要翼型相对厚度和相对弯度都很小，并且翼型的迎角也不大，那么翼型表面上压强的合力大小和方向就只受到厚度分布的轻微影响。对于这样的微弯薄翼，翼型的升力和力矩特性可以用气流绕它的中弧

4、线流动而求得，可以用薄翼理论来计算。2.2.1.1 压强和载荷根据伯努利方程，流动中某点的压强系数与该点的速度有如下关系： （27a）式中，v=i+j，和为扰动速度，为来流速度。对于小扰动情况，即，略去二阶小量后式（27a）简化为 （27b）弦向点x处下翼面与上翼面的压强与之差为载荷，用符号表示，为 （28）式中为载荷系数，。对于薄翼，整个翼型是由厚度分布和中弧线叠加而成的，图23。在小迎角情况下，根据线化方程和边界条件，翼型的压强系数可以表示成由厚度和弯度（包括迎角）贡献的叠加，即 式中，为当迎角和弯度时，由厚度产生的压强系数；为中弧线和迎角产生的压强系数。2.2.1.2 升力和力矩特性 薄

5、翼理论的结果。翼型的升力系数和绕翼型前缘的力矩系数为 （29） （210)式中，规定力矩使翼型前缘抬头为正，载荷与环境密度的关系为 (211）由薄翼理论有 （212）由式（29）至式（212）得 （213） （214a）用升力系数表示的力矩系数可写成 （214b）式（212）至式（214）中的多项式系数与中弧线方程的关系为 （215） 1. 翼型的升力特性将式（215）的系数代入式（213），改写为 （216）式中，为零升迎角，它代表零升力线与弦线的夹角，图24。它仅与中弧线形状有关。此式说明翼型的升力系数随几何迎角成线性变化。将对求导，得薄翼理论的升力线斜率 （217）2. 翼型的力矩特性对

6、于给定的翼型，式（214b）等号右边的第二项为常量，故与成线性关系，可将式（214b）改为 （218）式中，是零升力矩系数，它与翼型的升力或迎角无关，仅是翼型弯度分布的函数；是力矩系数对升力系数的导数。如果对翼型的1/4弦点取力矩，并利用式（218），可得 （219）显然，对于薄翼理论而言，1/4弦点力矩系数与升力系数（或迎角）无关，它就等于零升力矩系数。在翼型上有两个重要的特性点，一个是焦点（或称气动中心），另一个是压力中心。1） 翼型上存在这样一个点，该点的力矩系数与升力系数无关，这一点称为翼型的焦点。焦点的弦向相对量用表示。既然绕焦点的力矩与升力系数无关，故它是升力增量的作用点。因此，对

7、于前缘力矩系数又可写成 （220）将式（220）与式（218）的第一式相比较，可得基于薄翼理论的焦点位置 （221）2）翼型的升力作用线与弦线的交点称为压力中心，压力中心的弦向相对位置用表示。根据上述定义，将前缘力矩系数除以升力系数，可得 （222）从方程（215）知，和都与迎角无关，至取决于中弧线形状，故压力中心将随变化。对于对称翼型（弯度分布），=0，薄翼理论 压力中心与焦点重合，即。【例21】 某一翼型的弯度分布，试求该翼型的升力和力矩特性。解 该翼型的弯度分布沿的变化率为 由式（215）得 于是根据式（216）式（222）有 由最后一个式子可以看出，在零迎角下该翼型的压力中心，当迎角或

8、增大时，它将移向焦点。2.4.2 超声速薄翼型的线性化位流理论超声速线化速势方程为 （231）式中，。流动方程式（231）的通解为 （232）式中，f和g是自变量为和的任意函数。可以看到 常量，常量 （233a）的两族直线对于x轴的倾角分别为和，因此它们正好代表来流的两族马赫波，如下图所示。在翼型的上半平面流场中，函数代表翼型上表面所发出的扰动沿马赫线常量向下游传播到流场点（x，y）所产生的扰动速度位；而代表翼型下表面发出的扰动沿马赫线常量 向下游传播到流场点（x，y）所产生的扰动速度位。在超声速流场中，有意义的解是往下游传播的，而且，受到扰动的区域也只局限于前后缘马赫波之间。所以对上、下半平

9、面的流场的小扰动速度位分别是 ， （234a）可见，沿着翼型上表面的马赫波（常量）或沿下表面的马赫波（常量）为常量，而且，流场上沿着马赫波的两扰动速度分量和以及其他流动参数也都是常量。函数和科根据翼型绕流边界条件确定。设翼型的上表面方程为，由线化边界条件有 （235）对于上表面令，则有 （236）线化压强系数公式为 （237）联立式（235）、（236）和（237）。得 （238）类似地，如果翼型下表面方程为，则 （239）根据线化理论，翼型表面上任一点的压强系数与该点翼面对于来流方向的斜率成正比。由于物面上任一点相对于来流方向 倾角都很小，所以该点物面斜率科表示为 这样，式（238）和（23

10、9）科合并成 （240）式（238）至式（240）就是线化、二维超声速的基本关系。式（240）表明，物面上。任一点的压强系数与该点相对于来流的倾角成正比。相对于来流为压缩的物面倾角取正直，为膨胀的物面倾角则取负值。2.4.3 翼型的超声速气动特性对于薄翼，来流迎角很小，可认为翼型的整个气动力是由厚度、弯度和迎角产生的气动力的代数和。如图219所示，将x轴沿着翼弦方向，y轴与x轴垂直，则上、下翼面相对于来流的倾角可表示成 （241）式中，分别代表厚度（）翼型和弯度（）翼型表面上某一点的倾角，即 （242）将式（241）代入式（240），得任意形状翼面上下表面压强分布 （243）由式（234）可见

11、，对于任意形状的翼型，它的表面压强系数可认为是厚度分布、弯度分布和来流迎角所产生的压强系数的代数和。 1.升力特性从图219可见，作用在翼型微元上的升力为式中，和分别为翼型上下表面微元长度。翼型微元升力系数为 （244）将式（243）代入式（244），得 （245）由于在翼型前后缘，因此有 （246）将式（245）对x从零积分到b，并应用是（246），得 （247）从式（247）可见，在超声速线化理论中，薄翼型的升力系数与厚度和弯度分布无关，升力系数与迎角成正比。2.波阻力特性从图219可见，翼型微元上的阻力位 因为对于薄翼，所以上式可写为 （248）利用式（246），可以看出式（248）中的

12、积分为零，所以翼型的波阻力系数可表示为 （249）对于给定翼型都是知道的，令 （250）式中，分别为翼型的相对厚度和相对弯度；分别是与翼型厚度分布和弯度分布有关的几何常数。利用式（250），式（249）可改写为 （251） 式（251）等号右边第二项与迎角无关，仅与翼型厚度分布和弯度分布有关，对于弯度为零的翼型。零升波阻力系数为 （253）表22给出了几种超音速对称翼型的值，按式（250）计算，或由经验给出。菱形翼型的波阻系数是最小的。翼剖面简 图 正弦形 /8 四角形1/4(1-) 六角形 1/(1-a/b) 菱形 1 圆弧或抛物线形 4/3 亚声速翼剖面 2.54 3. 力矩特性如果忽略压

13、强弦向分量和其他高阶小量对俯仰力矩的贡献，那么对翼型前缘点的微元俯仰力矩可表示为 将式（243）代入上式，得翼型俯仰力矩系数 （254）由式（247）和式（254）可见，由于厚度产生的压强对翼弦是上下对称的，所以厚度对升力和力矩都无贡献。力矩系数是迎角和弯度作用的代数和，它们分别是 （255） （256）对于给定的翼型 ，式（256）中只是与翼型的弯度分布函数有关的参数，它的表达式为 （257a）或 （257b）压力中心距前缘的相对距离为 （258）翼型焦点距前缘的相对距离为 （259）在低速时，翼型的焦点；而在超声速时。由此可见，从低速到超声速焦点位置显著后移了，这是研究飞行器的稳定性与操作

14、性时必须注意的一个问题。由式（258）和（259）可见，由线化理论给出的压力中心位置和焦点位置仍然不随马赫数变化而变化，这与低、亚声速是相同的。【例22】 有一双凸面的翼型如图220所示，该翼型的上下表面方程分别为 和 试用线化理论计算该翼型在时的升阻和力矩系数以及压力中心位置。解 设翼型的迎角为，并由给定的翼型表面方程，根据式（241）和式（242），上下表面任一点切线与来流方向的倾角分别为 （1） 升力系数。升力系数仅与迎角有关，故 （2） 波阻系数。该翼型的厚度分布函数，弯度分布函数，所以零升波阻系数为 该翼型的总波阻系数 （3） 力矩系数。与该翼型的弯度分布函数相对应的任一点倾角为 将

15、上式代入式（256），得到弯度产生的绕前缘力矩系数该翼型总力矩系数为压力中心位置（4） 升阻比。升力Y与阻力X之比。在位流理论中，翼型的阻力就等于波阻力，故 第三章 机身的气动特性分析3.1 机身几何参数机身的几何参数列出如下：（1）为机身总长度；为机身头部长度；为机身尾部长度；为机身最大宽度；为机身最大高度；（2）为旋成体机身最大直径；为旋成体半径；（3）是旋成体长细比；为旋成体头部长细比；为旋成体尾部长细比；是旋成体尾部收缩比；（4）为机身相对高度；为机身相对宽度；为机身中轴线距机身纵轴线的距离。旋成体头部母线有各种形式，常用的头部母线方程介绍如下：（1） 锥形头部 （31）式中，为头部半

16、顶角。（2） 蛋形头部 （32）（3） 抛物线头部 （33）（4） （34）式中，（5） 哈克形头部 （35）（6）指数形头部 （36）以上介绍的母线方程是弹箭常用的头部形状，也是导弹和旋成体机身头部常用的形状。除上述几何参数外，机身的体积和机身表面积也是气动力计算中常用到的，对于旋成体机身来说，其体积可表示为 （37）其表面积为 （38）例如对于旋成抛物体积分可得 （39）式中，为机身最大横截面积。旋成体机身表面积与最大横截面积之比可近似表示为 （310）式中，为机身头部有进气道时引入的头部收缩比，为头部进气道口所在处的机身直径。 3.3 细长旋成体小迎角气动特性3.3.1 压强分布1.旋成

17、体公式对于像机翼那样扁平物体，在一阶近似的情况下，物面压强系数表示为 （312）式中，为x方向扰动速度，它与来流方向平行，是扰动速度位。对于旋成体，若用式（312），则不精确，还必须考虑扰动速度的平方项，以柱坐标表示为 （313）对于细长旋成体，式中x方向扰动速度比r方向及方向扰动速度，要小，故常略去，式（313）可简化为 （314）该式是在速度坐标系内的，在实用中常用机体坐标系内的关系式，以（x，y，z），（x，r，）表示体轴坐标，以（），（）表示风轴坐标系，两者之间的换算关系为两个坐标系下扰动速度的关系式为式中，。请注意式中带撇量为风轴系，不带撇量为体轴系。将以上转换关系代入式（314）中

18、，并注意式（314）中各物理量，此时都应该理解为是带撇量，最后可得 将上式配方，加上，则可写成（315）式（315）表明，只要求得旋成体绕流流场中的三个扰动分速度和，就可以求得任一下的。式（315）还表明，旋成体与扰动速度是非线性关系，一般情况下，求流场中某点，不能应用叠加原理，而只有对十分细长体求物面的特殊情况，才存在有叠加性，后面将会介绍。理论上处理旋成体绕流计算采用小扰动假设，将来流以小迎角绕旋成体流动看成是两种流动的叠加：一是来流以速度零迎角绕旋成体的流动，称之为轴向流；另一是来流以垂直于旋成体轴线绕旋成体流动，称之为横向流。2. 细长旋成体轴向流的物面压强分布对于轴向流旋成体（即）物面上所有母线上的分布是一致的，这时是轴对称流动，。根据细长旋成体理论可以求得细长旋成体轴向亚声速流动的扰动速度和，因而可求得轴向绕流时的压强分布。（316）式中，为旋成体横截面面积S沿轴线x的二阶导数；R是旋成体半径；为旋成体半径对x的导数；L是旋成体长度；。式（