引言:7张图,翼型分析全部完成
上一课我们用单雷诺数Re=210000跑完了NACA 4412翼型分析,3分钟出报告。这一课,我们把分析结果生成的7张图逐张解读,从中读出5个关键设计决策。然后,翼型分析就全部完成了——下一步进入3D机翼分析。
这7张图分为两类:4条核心曲线(直接用于设计决策)和3条辅助曲线(用于深入分析)。
一、7张图总览
| 图序 | 图表类型 | 重要性 | 核心结论 |
|---|---|---|---|
| 第1张 | Cl-Cd 极曲线 | 核心 | 评估整体升阻效率 |
| 第2张 | Cl vs α 升力曲线 | 核心 | 确定失速迎角和巡航迎角 |
| 第3张 | Cm vs α 力矩曲线 | 核心 | 判断稳定性、尾翼设计依据 |
| 第4张 | Cl vs Xtr top 转捩点 | 辅助 | 边界层转捩分析 |
| 第5张 | Hinge Fy vs α 铰链力 | 辅助 | 操纵面铰链载荷 |
| 第6张 | Cd vs α 阻力曲线 | 核心 | 阻力随迎角变化规律 |
| 第7张 | Cl/Cd vs α 升阻比 | 核心 | 确定最佳巡航迎角 |
XFLR5默认只显示5个图表,第6张(Cd vs α)和第7张(Cl/Cd vs α)需要通过右键图表 → Current Graph → Define Graph Setting,修改Variables选项卡中的Y轴变量来调出。改完点OK后,必须按R键刷新才能看到变化。
二、第1张:Cl-Cd 极曲线
极曲线是翼型分析的"名片",横轴是阻力系数Cd,纵轴是升力系数Cl。单雷诺数分析的好处是曲线只有1条绿色线,比12条线叠在一起清晰得多。
关键信息
- 最大Cl约1.4(曲线最高点)
- 失速区域:曲线顶端开始变平,Cl=1.3~1.4
- 最佳效率点:曲线最靠左的点(最小Cd对应的Cl),约Cl=0.3~0.5
当前工况点(Alpha=15°)
| 参数 | 数值 | 说明 |
|---|---|---|
| Cl | 1.39367 | 接近最大升力 |
| Cd | 0.06673 | 阻力较大(已失速) |
| Cl/Cd | 约20.9 | 失速状态,效率低 |
图1:Cl-Cd极曲线(Graph 2),横轴Cd阻力系数0~1.0,纵轴Cl升力系数0~2.0。红色标记点为Alpha=15°工况,Cl=1.39367,Cd=0.06673。曲线左侧最靠前的点即为最小阻力点
三、第2张:Cl vs α 升力曲线(最核心!)
这是最重要的曲线——升力系数随迎角的变化。横轴是迎角α(-5°到15°),纵轴是升力系数Cl(0到1.4)。
曲线三段特征
| 阶段 | 迎角范围 | 曲线表现 | 物理意义 |
|---|---|---|---|
| 线性段 | -5° ~ 8° | 近似直线,斜率约0.1/度 | 气流附着,升力线性增长 |
| 过渡段 | 8° ~ 12° | 斜率减小,曲线弯曲 | 气流开始分离 |
| 失速段 | 12° ~ 15° | 曲线变平,Cl接近最大 | 深度失速 |
关键数据
| 参数 | 数值 | 设计意义 |
|---|---|---|
| 0°时Cl | 约0.4 | 正弯度翼型,0度也有升力 |
| 升力斜率 | 约0.1/度 | 每增1度迎角,Cl增0.1 |
| 失速迎角 | 约12-14° | 超过后升力不再增长 |
| 最大Cl | 约1.39 | 翼型最大升力系数 |
对1公斤飞机的意义:巡航需要Cl约0.7,对应迎角约4度,失速迎角12度,安全余量8度。NACA 4412完全满足要求。最大升力 = 1.39 × 0.25 × 1.225 × 15² ÷ 2 ≈ 4.9kg,远大于1kg重量。
图2:Cl vs α升力曲线(Graph 1),横轴迎角α(-5°~15°),纵轴Cl(0~1.4)。曲线在-5°~8°为线性段,8°~12°斜率减小开始过渡,12°以上变平进入失速。0°时Cl约0.4,证实正弯度翼型特性
四、第3张:Cm vs α 力矩曲线
力矩曲线横轴是迎角,纵轴是俯仰力矩系数Cm(全为负值)。
| 参数 | 数值 | 设计意义 |
|---|---|---|
| Cm0(0°时) | 约-0.09 | 负值,有低头力矩,需要尾翼配平 |
| 线性段斜率 | 接近0 | 翼型自身几乎中性稳定 |
| 失速后Cm | 温和变负 | 自动低头恢复,安全教练机特性 |
设计结论:Cm0约-0.09意味着NACA 4412有弯度翼型天然产生低头力矩,必须有水平尾翼产生抬头力矩来平衡。但Cm0绝对值不大,常规尾翼设计即可配平。失速时Cm温和变化(不急剧变负),这对教练机是安全特性。
图3:Cm vs Cl力矩曲线,横轴Cl(0~11.4),纵轴Cm(-0.10~-0.02)。Cm全为负值,0°时Cm0约-0.09,表示翼型有低头力矩。曲线斜率接近0说明翼型自身中性稳定,需要水平尾翼配平
五、第4张和第5张:辅助曲线
第4张:Cl vs Xtr top(转捩点分析)
横轴是上表面转捩点位置(0~1.0弦长),纵轴是Cl。当前工况点转捩点在弦长3.4%处,非常靠前,说明15度大迎角下气流几乎在前缘就转捩为湍流。这是辅助分析图,不直接用于设计决策。
图4:Cl vs Xtr top转捩点分析,横轴上表面转捩点位置(0~1.0弦长),纵轴Cl(0~1.8)。红色标记点在Xtr top=0.034处,说明15°大迎角下气流在3.4%弦长处就转捩为湍流
第5张:Hinge Fy vs α(铰链力分析)
这是襟翼/操纵面铰链力分析图,用于评估作用在操纵面铰链上的气动载荷。当前分析没有设置襟翼偏角,正常飞行迎角下铰链力接近0,舵机选型不是问题。这条曲线在操纵面设计阶段才有参考价值。
图5:Hinge Fy vs α铰链力分析,横轴迎角α,纵轴铰链力系数。当前只有一个数据点在α=15°处,铰链力接近0。此图用于操纵面设计阶段,目前无襟翼偏角,参考价值有限
六、第6张:Cd vs α 阻力曲线
这条曲线需要通过Graph Settings调出来。横轴是迎角α(-5°到15°),纵轴是阻力系数Cd(0到0.07)。曲线呈U型/抛物线。
| 区域 | 迎角范围 | Cd表现 | 物理意义 |
|---|---|---|---|
| 最小阻力区 | -2° ~ 4° | Cd约0.009-0.015 | 最优工况,摩擦阻力为主 |
| 缓增区 | 4° ~ 10° | Cd缓慢上升 | 压差阻力开始增加 |
| 急增区 | 10° ~ 15° | Cd急升至0.067 | 流动分离,失速前兆 |
关键数据
- 最小Cd:约0.009-0.010,在α≈0°附近
- 巡航Cd(α=4°):约0.015,阻力 = 0.015 × 0.25 × 1.225 × 15² ÷ 2 ≈ 0.52N
- 失速Cd(α=15°):0.06673,阻力暴增7倍
- 阻力发散迎角:约10°,超过后阻力急剧上升
图6:Cd vs α阻力曲线,横轴迎角α(-5°~15°),纵轴Cd(0~0.07)。曲线呈U型,0°~4°附近Cd最小约0.011(阻力斗区),10°以上阻力急升至0.067,对应流动分离
七、第7张:Cl/Cd vs α 升阻比曲线
这是最后一条也是最关键的设计依据曲线。横轴是迎角α(-5°到15°),纵轴是升阻比Cl/Cd(0到75)。曲线呈单峰形状,先升后降。
| 区域 | 迎角范围 | Cl/Cd表现 | 物理意义 |
|---|---|---|---|
| 上升区 | -5° ~ 6° | 从低值快速上升 | 升力增长快于阻力 |
| 峰值区 | 6° ~ 7° | 达到最大约70-75 | 最佳巡航效率 |
| 缓降区 | 7° ~ 12° | 缓慢下降 | 阻力开始增大 |
| 急降区 | 12° ~ 15° | 急降至约20 | 失速,效率暴跌 |
最关键的发现:最佳巡航迎角是6-7度,不是之前根据升力曲线选的4度!4度巡航时Cl/Cd约55,6-7度巡航时Cl/Cd约75,效率提升36%。如果追求最大航程,应该在6-7度巡航。
图7:Cl/Cd vs α升阻比曲线,横轴迎角α(-5°~15°),纵轴Cl/Cd(0~75)。曲线呈单峰形状,6°~7°达到峰值约75(最佳巡航效率),12°以上急降至约20(失速)。这是确定巡航迎角最关键的曲线
八、5个关键设计决策
7张图分析完毕,从中读出5个关键设计决策:
| 决策 | 数据来源 | 结论 |
|---|---|---|
| 1. 巡航迎角 | Cl/Cd vs α + Cl vs α | 4-7度(4度保守,6-7度最高效) |
| 2. 失速迎角 | Cl vs α | 12度,安全余量5-8度 |
| 3. 尾翼设计 | Cm vs α | Cm0≈-0.09,需要水平尾翼配平 |
| 4. 升力余量 | Cl vs α | Clmax≈1.39,最大升力4.9kg >> 1kg |
| 5. 翼型适用性 | 全部曲线 | NACA 4412完全适合1公斤教练机 |
翼型分析总结
NACA 4412在Re=210000下:巡航迎角4-7度,失速迎角12度,最大升力系数1.39(对应4.9kg升力),最小阻力系数0.009,最大升阻比75。Cm0=-0.09需要尾翼配平,失速特性温和。翼型分析全部完成,结论:完全适合1公斤教练机。
九、如何调出隐藏的曲线
XFLR5默认显示5个图表,其中只有3条是核心曲线。另外2条核心曲线(Cd vs α 和 Cl/Cd vs α)需要手动调出:
我试过了,改了以后,点确定,图没有变化。
这是XFLR5的一个交互特性——修改完变量点OK后,图表不会自动刷新,需要手动触发重绘。
正确的操作步骤
- 按数字键1,切到Graph 1单独大图显示
- 在黑色图表区域上点右键
- 右键菜单选 Current Graph → Define Graph Setting(注意不是All Graph Settings)
- 在Variables选项卡中,Y轴选要显示的变量(如Cd或Cl/Cd),X轴选α
- 点OK
- 按R键刷新——这一步最关键!
关键发现:Graph Settings改完点OK后,必须按R键强制刷新,图表才会按新变量重绘。这是XFLR5的已知交互特性,不是软件bug。另外,要用右键菜单的"Current Graph → Define Graph Setting"来改单个图表变量,不要用"All Graph Settings"(那个只改通用显示属性,不改变量)。
十、下一步:从2D翼型到3D机翼
翼型分析完成了,但这只是第一步。翼型是2D截面,告诉我们的是"如果机翼无限长,性能会怎样"。真实机翼有翼展、有翼尖,气流会从翼尖绕过去产生诱导阻力。3D机翼分析就是算这个。
3D机翼分析能得到什么
| 数据 | 意义 |
|---|---|
| 3D升力曲线 | 考虑翼尖效应后的实际升力(Clmax会比2D低) |
| 诱导阻力 | 翼尖涡流产生的额外阻力 |
| 升力沿展向分布 | 看升力在翼展方向怎么分布 |
| 下洗角 | 翼尖涡流导致的气流下偏 |
| 奥斯瓦尔德效率因子 | 翼型效率到机翼效率的折减 |
需要定义的机翼参数
| 参数 | 数值 | 说明 |
|---|---|---|
| 翼展 | 1.6米 | 左右翼尖之间距离 |
| 弦长 | 0.156米 | 翼型前后缘距离 |
| 翼面积 | 0.25㎡ | 翼展×弦长(梯形翼等效) |
| 展弦比 | 10.24 | 翼展²÷翼面积 |
| 翼型 | NACA 4412 | 刚分析完的翼型 |
重要预期:3D机翼的Clmax会比翼型低(约1.1-1.2 vs 翼型的1.39),阻力会比翼型高(多了诱导阻力),但整体性能仍然足够1公斤飞机使用。展弦比10.24属于中等偏大展弦比,诱导阻力相对较小,适合滑翔型教练机。
十一、学习小结
从第五课第一次打开XFLR5,到第十一课翼型分析全部完成,我们经历了:
- 第五课:认识XFLR5界面,生成第一个NACA 4412翼型
- 第六课:理解为什么从翼型开始分析
- 第七课:XFLR5 v6.62实操,生成翼型并启动分析
- 第八课:分析跑完能得到什么——5个设计决策和4条曲线
- 第九课:12组雷诺数的曲线怎么看
- 第十课:只跑1组雷诺数,3分钟出精准报告
- 第十一课(本课):7张图逐张解读,翼型分析全部完成
翼型分析阶段告一段落。下一课开始,我们进入3D机翼定义和分析——把2D翼型数据变成真实机翼的气动特性。这才是航模设计中最接近真实飞行的部分。