简介这份来自北京航空航天大学精品课程的空气动力学基础教学课件以PDF格式呈现共1个文件压缩包大小19.65MB适合航空航天类专业学生、教师及相关工程技术人员系统学习与参考。课件内容涵盖绪论、流体的基本属性、流体静力学、流体运动学、不可压缩无粘流体平面势流、粘性流体动力学基础、边界层理论及其近似、可压缩高速流动基础等核心章节并对空气动力学的研究方法、分类、发展进程及应用场景做了系统梳理结构完整便于按章节循序渐进地学习。目前已有675人学习使用。学习者可通过本课件建立空气动力学整体知识框架理解飞行器绕流现象与受力分析的基本原理为后续从事航空航天工程实践或深入研究打下扎实基础。1. 一份北航精品课程PDF凭什么值得从头啃一遍拿到《空气动力学基础(北航精品课程)-刘沛清.pdf》很多人第一反应是“又是一堆公式”。但我见过太多做CFD的人能跑通算例却说不清升力从哪来残差曲线一乱只能靠调松弛因子碰运气——这是基础没补好。这门课是北航飞行器类专业的核心课刘沛清老师把流体力学筛了一遍留下的都是直接服务于飞行器受力的内容。它适合刚入行的气动仿真工程师、在考研或读研的学生、做总体设计需要估算气动数据的人。它不教你花哨技巧每一页都在回答“受力从哪来、什么条件下成立、算错了会怎样”。这篇笔记写的是我怎么用它搭起一套可复现的学习框架以及踩过的坑。2. 空气动力学基础的骨架控制方程与无量纲参数这份PDF的目录通常先讲大气环境、流体属性然后进入流动控制方程再按“不可压缩、可压缩、粘性”三大块展开。如果你跳着读很容易被符号淹没。我的建议是先抓骨架控制方程决定物理模型无量纲参数决定相似性流动分界决定简化方向。这三件事不立住后面读翼型、读边界层、读激波都会像是在背字典。2.1 从NS方程到欧拉方程简化条件决定适用范围N-S方程是空气动力学的地板但直接从三维非定常N-S方程推起会劝退很多人。刘沛清老师的讲义做法是先写出守恒形式然后一步步删除项告诉你每一项在什么条件下可以忽略。比如不可压缩、无粘、无热传导情况下N-S方程退化为欧拉方程再叠加无旋条件变成速度势方程这就是全位势方程的来源。这里要记住的不是方程的形式而是“删掉哪一项对应什么假设”。我一般看N-S方程会先看它的四个来源质量守恒、动量守恒三个方向、能量守恒。低速、无热源流动中能量方程可以和动量方程解耦剩下连续性和动量方程就够用了。不可压缩时密度视为常数连续性方程里密度对时间的偏导直接清零动量方程里的粘性应力项在边界层外很小但在物体表面附近又必须保留。所以同一个流动在不同位置可以用不同方程这是初学者最容易困惑的地方——为什么“无粘”还能描述真实机翼外面的流场因为粘性的作用局限在边界层和无粘区域边界的相互作用里。在做实际工程估算时我会先把流动归入某一类低速不可压无粘、低速有粘、高速无粘、高速有粘带激波。这个分类直接决定选什么公式。比如估算翼型升力用势流理论足够但估算阻力必须回到边界层估算超音速激波位置要用斜激波关系式而激波后的边界层问题又得回到N-S。PDF里大量篇幅在讲这些简化条件背后的目的就是让你拿到一个问题时能判断该用哪一层工具。方程/模型速度适用范围忽略内容典型应用拉普拉斯方程低速 Ma0.3 无粘粘性、旋涡、压缩性势流法求升力、压力分布Euler方程跨音速/超音速无粘粘性、热传导激波位置、波阻估算边界层方程近壁区粘性主导法向压力梯度近似为零摩擦阻力、转捩判断完整N-S方程全速域无CFD底层、分离流/湍流直接模拟这个表是我从PDF的公式脉络里抽出来的。它不是为了让你背方程而是告诉你当你跑CFD时求解器内部默认用哪一套方程你就能解释为什么某些工况会发散。比如在分离流大攻角下有些基于欧拉方程加边界层耦合的工具会给出乐观的升力因为把分离后的压力分布算得太“干净”了。2.2 三个绕不开的无量纲参数雷诺数、马赫数、升力系数无量纲参数是空气动力学从“玄学”变科学的底气。三个参数必须形成条件反射雷诺数决定流动状态马赫数决定压缩性升力系数把形状和工况解耦。刘老师的讲义里所有实验数据都以这三个参数为主轴你只要摸清它们的量级和物理意义就能看懂大多数图表坐标轴。雷诺数定义为Re ρVL/μ其中ρ是密度V是特征速度L是特征长度μ是动力粘度。它表示惯性力与粘性力的比值。Re低的时候流动是“粘着的”容易层流Re高的惯性力主导扰动容易放大成湍流。实际工程里小型无人机翼弦0.2m、速度15m/sRe大概20万左右民航客机翼弦3m、速度250m/sRe超过5000万。两个Re差了三个数量级同样一块翼型的最大升力系数和失速特性完全不同。很多风洞数据在Re百万以下测得直接用到飞机上必然翻车。马赫数Ma V/aa是当地声速。Ma0.3可以认为密度不变这是不可压缩近似的标准门槛。从0.3到0.8是亚跨音速翼面上出现局部超声速区但主流仍是亚声速必须考虑压缩性修正Ma1是超声速存在激波。课程里对Ma的讨论不只是分类还会给出一套等熵关系式T0/T 1 (γ-1)/2 · Ma²p0/p (T0/T)^(γ/(γ-1))。这两个公式在计算风洞试验条件、推算飞行器驻点温度时特别常用。比如飞行器Ma2时驻点温度按γ1.4算大约是静温的1.8倍如果不做热防护设计结构强度会直接出问题。升力系数Cl L/(q∞·S)其中q∞是动压S是参考面积。它把升力从速度中解耦出来方便对比不同翼型。薄翼型线性范围内Cl≈2π(α-α0)α0是零升攻角。这个线性关系在真实流体中只在一定攻角内成立失速后Cl急降。学习时要把Cl-α曲线从讲义中抽出来反复对照附着流阶段斜率大约等于2π单位弧度非线性段是分离的开始失速点是最大Cl对应的攻角。这三个阶段对应三种物理行为也对应三种CFD收敛难度。2.3 不可压缩与可压缩流动的分界在哪里分界不是一条严格的线而是一个工程决策点。通常以Ma0.3为界。在这个Ma下密度变化最多5%工程上当作不可压缩处理误差可以接受。超过0.3后密度变化带来的流管面积和速度之间的耦合不能忽略需要使用可压缩质量方程、动量方程和能量方程的联立形式。这个分界在PDF里反复出现也在求解器里对应两套压力修正算法不可压缩用压力泊松方程通过压力修正保证连续可压缩用密度基方法直接解守恒变量。在低速部分课程会先引入流函数和速度势建立无旋流的基本解均匀流、点源、点涡、偶极子。叠加这些基本解可以得到圆柱绕流的势流解但没有升力——这就是著名的达朗贝尔佯谬。加入点涡后通过库塔-茹可夫斯基定理L ρ∞·V∞·Γ得到升力Γ是环量。升力从“环量”来这是理解升力的正确入口而不是单纯靠“上表面快下表面慢”那套直觉。进入可压缩部分后课程会从一维等熵流讲起再讲正激波、斜激波和膨胀波。初学者最容易走神的地方是总压和静压的区别。总压是把气流等熵滞止到零速后得到的压力它是流动“有效做功能力”的量度静压就是当地真实气体分子撞击壁面的压力。皮托管测的是总压静压孔测的是静压。跨过激波时总压下降静压上升熵增加。这些概念在飞机设计里直接决定发动机进气道的总压恢复系数也决定气动加热程度。学习时不要只记公式要在PDF的每个例题里标出哪些量是滞止量哪些是当地量。3. 把PDF读成技能按模块推进的学习路线与笔记方法PDF课件和教材不一样。教材有完整的前后推导课件往往把过程压到只剩关键结论。直接从头翻课件很容易在导数、旋度、斯托克斯定理之间迷路。我的建议是不要按页序读把它当成一个“技能树”来铺设先建立受力直觉再按低速、可压缩、粘性三个模块推进最后把每个结论整理成可检索的知识卡片。3.1 第一步建立坐标系、受力定义和速度三角打开PDF第一遍先找到翼型受力图把坐标系、速度三角形、升力和阻力方向定义抄一遍。注意升力永远垂直于来流方向而不是垂直于翼弦。很多初学者画图时随手一画把升力画成和机翼表面垂直后面所有力矩分析全部错位。正确做法是来流速度V∞水平向右攻角α是翼弦与来流的夹角升力垂直V∞阻力平行V∞气动力矩以某参考点为正。这张图几乎决定后面所有公式的方向和符号。然后去PDF里找“压力中心”和“焦点”的定义。压力中心是气动力合力与弦线的交点它随攻角移动焦点是气动力矩系数不随攻角变化的参考点也叫气动中心。你做纵向静稳定性分析时只需要看焦点相对于重心的位置焦点在重心之后扰动攻角产生的低头力矩会恢复原状即静稳定焦点在前扰动会发散。这个判断不用算压力中心算焦点就够了。我建议你把这三个参考点重心、焦点、压力中心画在一张图上标出典型攻角下各自位置再用一个简单算例验证稳定性。比如一个常规布局飞机焦点约在平均气动弦25%处重心在30%处那么焦点在重心前实际上重心在前才是稳定。画图时别搞反。这一步过关了后面读升力线斜率、下洗流、机翼-尾翼干扰才有坐标感。3.2 第二步按低速、可压缩、粘性三个模块推进学习顺序不要按PDF目录从前到后闷头啃。我推荐倒着安排因为低速无粘模块是基础粘性模块解释阻力和分离可压缩模块最后加宽带压缩性和激波。先低速无粘再把无粘结果和真实结果的差距用边界层补上最后处理高速效应。低速无粘部分重中之重是势流理论和环量。手推至少一遍圆柱绕流加点涡的叠加均匀流加偶极子加点涡通过速度势求表面速度再用伯努利求压力对压力沿圆柱表面积分得到升力为ρ∞V∞Γ阻力为零。这个推导串起速度势、流函数、叠加原理、伯努利、库塔-茹可夫斯基定理五件事。做完这一遍你会突然看懂翼型升力线理论里的“平板替换”是怎么回事。粘性部分核心是边界层概念。刘沛清老师在讲义里会对边界层厚度、位移厚度、动量厚度做严格定义。位移厚度是边界层内因速度亏损而把外流“推开”的等效距离它影响翼型压力分布动量厚度直接进入阻力公式。边界层是否分离主要看逆压梯度大小和雷诺数。一个常见误用是“只要Re大边界层就薄阻力就小”实际上Re大往往是湍流边界层摩擦阻力系数反而相对下降但转捩位置提前可能带来更大的分离趋势。这里要理解层流和湍流在速度分布上的差别。可压缩部分不要急于去背激波公式。先掌握一个一维变截面管流问题缩放喷管里面积变化如何影响速度和压力。然后理解正激波的基本关系比如马赫数越大总压损失越大。用PDF里的表格或图查一道“马赫2气流过正激波后马赫数约为0.577静压比约为4.5”的题。这个数量关系能让你知道超音速进气道设计为什么难损失太大。每个模块学完找一道例题不看答案按PDF给的步骤手算。手算不是要你数值精确而是强迫自己停下来看单位和适用条件。算完后用讲义后面的参考数据对照误差在一个数量级内即可。3.3 第三步用算例和可视化把公式变成直觉公式背得再熟不如看一次流场图。PDF里会有翼面压力分布图升力系数曲线也可能有油流或PIV实验照片。你要做的事情是把公式结果和这些可视化对应起来。比如看到上表面吸力峰附近等压线密集就应该想到那里速度高、压力低、局部静压最低。看到失速状态下吸力峰崩溃就应该想到边界层分离、环量损失。我常用的验证方法是做一个手算小表格给定翼型弦长0.15m速度30m/s空气密度1.225kg/m³粘度1.8e-5Pa·sRe≈30.6万。按照PDF里的经验值查这个大攻角下可能的Cl范围比如18万Re的小型无人机在攻角12°时Cl可能在1.1到1.3之间。然后手算升力LCl·0.5·ρ·V²·S得到单位展向升力大约多少牛。这个数字再去和CFD结果对比差得远了就说明计算域或湍流模型有问题。这个“手算预判”的习惯能让你的数值模拟不再是黑匣子。可视化方面PDF里如果有实验影子图或纹影图先看激波前后流线的折转角再对照斜激波理论预测的折转角。一般偏差都来自粘性边界层对等效楔角的影响。别怕偏差怕的是你完全没有预测值。3.4 笔记怎么记参数、公式、适用条件三栏卡片我当年是吃过亏的拿着PPT从头抄了一遍合上文档什么也没留下。后来改成卡片法每个知识点写一张卡片分三栏参数定义与单位、公式与各符号含义、适用条件与反例。比如“马赫数”卡片定义是V/a公式MaV/a其中asqrt(γRT)适用条件是无所谓但用途是判断压缩性反例是有人把Ma1全叫亚音速其实跨音速已经在翼面局部超音速。公式卡片的重点不是公式本身而是“此公式在什么假设下才成立”。例如伯努利方程p 0.5ρV² ρgh const只沿同一条流线成立且要求定常、无粘、不可压。如果你在边界层里用伯努利就是自找麻烦。把这个“适用条件”写得越严格后面复习越省力。我会把PDF里的所有相关曲线重新画一遍草图升力系数随攻角曲线、极曲线、压力系数分布、激波前后参数表。画图不用精美关键在于自己动手标注每个拐点的物理含义。失速点对应什么极曲线上最小阻力点对应什么攻角这些一旦画过就不会忘。4. 避坑自学空气动力学最容易翻车的五个概念下面这条是通读讲义、同时做了多轮工程估算后总结的。以下每一个坑我都亲身踩过有的在考试里扣分有的在仿真结果里绕了几天才发现。按“现象、原因、解决”三步写出来希望能帮你少走这些弯路。4.1 压力中心与焦点的混淆稳定性算错一整片现象做飞机纵向静稳定性估算时有人用压力中心位置随攻角变化来判定稳定算出焦点在重心之后却把结论写反或者拿两个不同攻角下的压力中心说事得出“压力中心随攻角变化所以不稳定”的错误结论。原因压力中心是气动合力作用点它受攻角影响明显。纵向静稳定性看的是“攻角扰动后气动力矩是否产生恢复效应”这需要用焦点气动中心来定义而不是压力中心。焦点被定义为力矩系数斜率dCm/dα为常数的参考点实际上就是矩心变化不影响斜率的那一点。判定静稳定看焦点与重心的相对位置焦点在前重心在后则不稳定。解决把所有力矩分析都换成焦点。从讲义里找“焦点定义和推导力矩系数”的那一节用一次代数运算验证如果把力矩系数对攻角求导就会发现焦点处的导数与参考点无关。计算时用一个固定参考点比如平均气动弦前缘求出dCm/dα和dCL/dα再用公式x_focus x_ref - (dCm/dα)/(dCL/dα)无因次化求出焦点然后对比重心位置。这一步不能省。4.2 把伯努利方程当作一切升力的解释现象解释升力时只讲“上表面气流快压力低下面慢压力高”然后试图用伯努利反推上表面流速误差很大。还有人用这个去解释平板大攻角分离流算出的升力方向和大小都对不上。原因伯努利是能量方程沿流线的体现它把速度场和压力场联系起来但不解释为什么上表面流速会被加快。升力的根源是流体偏转带来的动量通量变率等价于环绕翼型存在净环量可以用库塔-茹可夫斯基定理表示。只有在附着流、无粘、定常的条件下伯努利才是一个便捷的计算工具把它当物理因果就反了。解决遇到升力问题先画控制体积用动量定理写出升力等于控制面上下压力积分或动量通量积分。然后用自由涡系起动涡附着涡说明环量来源。最后再用伯努利去算压力分布得到和库塔-茹可夫斯基一致的结果。这一套流程走完你才能对“升力从哪来”有完整认识。4.3 忽略雷诺数相似风洞数据和飞行数据对不上现象从文献里查一个翼型的最大升力系数直接用于自己的飞行器设计结果试飞时发现失速速度比设计值高升力不足。原因文献数据对应某个雷诺数可能是Re1e6或3e6而你的实际工况只有Re2e5。低雷诺数下边界层更容易分离最大升力系数明显下降失速攻角提前。雷诺数相似是风洞实验的第一前提少了这个数据只能是装饰。解决看数据时先确认Re范围。如果是小尺寸低速无人机就不要采用大型飞机的Re数据。如果没有对应Re的数据可以降低攻角安全裕度或者用转捩带提前固定转捩点来模拟高Re。CFD也是一样湍流模型和转捩模型对Re敏感边界条件里输入速度、弦长、粘度后先算一下Re是否和参考工况一致。这个自查两分钟能省掉一整晚的重新计算。4.4 把摩擦阻力与压差阻力混在一起调参现象模型阻力偏大有人认为是表面太粗糙把CFD壁面粗糙度调到0结果阻力几乎没变甚至分离区更大。原因气动阻力分解为摩擦阻力和压差阻力。粗糙度主要影响湍流摩擦阻力但压差阻力由逆压梯度引起的分离区大小决定。如果机翼已经在大攻角下有分离压差阻力占总阻力大头抛光表面当然没效果。解决用CFD后处理把阻力拆开看。从几何表面做积分摩擦阻力用壁面剪切应力面积分压差阻力用表压沿法向投影积分。如果压差阻力占比超过70%你的调整方向是翼型弯度、前缘半径、加装涡流发生器或调整攻角而不是改粗糙度。PDF里型阻极性和极曲线的部分会告诉你最小阻力点在哪里尽量在那个攻角附近飞行是最直接的调整。4.5 总压、静压、动压测出来的是什么现象用皮托管测机翼表面压力把读数直接当静压算压力系数结果压力系数偏差极大。原因皮托管对着来流中心孔感受到的是滞止压力即总压侧面孔才测静压。如果探针安装角度不正总压管还可能感受到不完全滞止读数异常。把总压当静压相当于动压被重复计入速度越高压差越大。解决先用一个三孔或五孔探针的校准表确认偏航角和俯仰角对测压的影响。手算时记住总压p0 静压p 动压q。在低速不可压缩流中q 0.5ρV²在可压缩流中需用等熵关系修正q p0·(1 - (p/p0)^((γ-1)/γ))·γ/(γ-1)实际上是总静压比公式的变形。无论做实验还是CFD都要明确每个输出的压力是总压还是静压。CFD里总压常标为p0静压标为p别混着用。5. 一个值得养成的验证习惯动压与马赫数的快速手算这套课程给我留下最深的影响不是背下了哪个公式而是养成了“开工前先手算一个量级”的习惯。现在我拿到一个气动问题无论跑不跑CFD都会先算动压、马赫数和雷诺数再把预期的升力、阻力数量级写在纸上。比如一个飞行器在海平面飞行速度137m/s其中弦长1.2m标准大气密度1.225kg/m³粘度1.79e-5Pa·s。那么动压 q 0.5×1.225×137² ≈ 11500Pa马赫数Ma 137/340 ≈ 0.4已经超过0.3不可压缩近似需要谨慎使用至少要做普朗特-格劳厄特修正。雷诺数Re ≈ 1.225×137×1.2/1.79e-5 ≈ 1.12e7。如果预期Cl0.8单位面积升力就是 0.8×11500 ≈ 9200Pa对应升力比重量超出多少马上能判断设计是否合理。在可压缩修正上有一个足够用的技巧当Ma在0.3到0.7之间用普朗特-格劳厄特修正Cp Cp_incompressible / sqrt(1 - Ma∞²)评估翼型表面压力系数变化。这样即使不做完整CFD也能判断压缩性会不会让吸力峰提前达到临界马赫数。这个细节我曾在一次方案讨论中用来质疑一个高翼载无人机的巡航速度设计避免了一次绿皮初样阶段的返工。超音速工况下我会先查一维等熵流表或直接用公式算总压比、静压比。Ma2.0时静压比p/p0约0.128驻点温度比T0/T约1.8。把这个记在心里看到风洞报告里的压力探针读数就能快速排除一个量级误差。每一次这样的手算都在验证你是否真理解了PDF里的物理模型而不是把CFD结果当黑匣子吞下去。我自己的教训是有次为了赶节点跳过手算直接起了一个跨音速算例结果求解器把速度入口压力给成表压而非绝对压力导致全流场密度差了一倍。这种错靠残差曲线看不出来靠手算一分钟就能发现。现在我只信两条先定参数量级再动仿真软件。希望帮到你。本文还有配套的精品资源点击获取