上一篇文章中,我们完成了水平尾翼的定义和Type 2稳定性分析的设置:
设置完成后,点击Analyze按钮准备运行分析——但曲线没出来。这是因为还缺少一个关键前置步骤。
3D粘性分析需要2D翼型极线数据作为输入。机翼用的NACA 4415翼型在之前的步骤中已经完成了2D分析,但尾翼用的NACA 0012还没有跑过2D分析——XFLR5找不到对应的粘性数据,3D分析就无法生成曲线。
解决方法:切到2D模块(Ctrl+1),对NACA 0012翼型运行多线程批量分析(Multi-threaded Batch Analysis):
| 参数 | 值 | 说明 |
|---|---|---|
| 翼型 | NACA 0012 | 尾翼对称翼型 |
| 雷诺数范围 | 100,000 ~ 250,000 | 覆盖尾翼实际工作雷诺数 |
| 雷诺数步长 | 50,000 | 4个雷诺数点 |
| 迎角范围 | -5° ~ 15° | 覆盖可能的工作范围 |
| 迎角步长 | 0.5° | 精细分辨率 |
| 分析类型 | Type 1 (Fixed Speed) | 2D翼型分析 |
2D分析完成后,切回3D模块(Ctrl+3),重新运行Analysis → Analyze,这次就能正常生成曲线了。
在Plane Editor中确认T2-VLM2分析对象已选中,检查右侧Plane Analysis面板的参数设置:
| 参数 | 值 |
|---|---|
| 分析类型 | Type 2 (Fixed Lift) |
| 方法 | VLM2 (Ring vortex) + Viscous |
| 迎角范围 | -2° ~ 12° |
| 步长 | 0.5° |
| Sequence | 勾选(批量运行) |
| Store OpPoint | 勾选(保存 Operating Points) |
以上参数都在Analysis Definition(分析定义)对话框中设置。下面三个截图分别展示Polar Type、Analysis和Inertia三个标签页的实际配置:
点击Analyze按钮,等待计算完成。计算进度会在底部状态栏显示,760个面板的VLM2分析通常几十秒内完成。
分析完成后,需要从3D View切换到Polar View才能看到曲线:
如果切换后看不到曲线,可能是分析没有成功运行。检查方法:在对象浏览器中展开T2-VLM2,看下面有没有Operating Point数据点。如果没有,需要重新点击Analyze按钮运行分析。
Cm-Alpha曲线(俯仰力矩系数-迎角曲线)是判断俯仰稳定性的最直接证据。
判断标准:Cm-Alpha曲线斜率为负 = 俯仰稳定。
物理意义解读:
从曲线数据可以读取:
| 迎角 α | Cm(俯仰力矩系数) | 物理意义 |
|---|---|---|
| -2° | ≈ -0.015 | 低迎角,低头力矩较小 |
| 0° | ≈ -0.04 | 零迎角,已有低头力矩 |
| 3° | ≈ -0.09 | 巡航迎角附近,低头力矩增大 |
曲线斜率(Cm_α)约为 -0.020/°,负值确认了飞机具备俯仰静稳定性。
光知道"稳定"还不够,还需要知道"有多稳定"——这就要用到静稳定裕度(Static Margin, SM)。
SM的计算公式:
SM = -Cm_α / CL_α
其中Cm_α是俯仰力矩曲线斜率,CL_α是升力曲线斜率。需要切换Y轴从Cm改为CL,查看CL-Alpha曲线:
从CL-Alpha曲线读取:
| 参数 | 值 | 来源 |
|---|---|---|
| CL_α(升力曲线斜率) | ≈ 0.097/° | CL-Alpha曲线线性段斜率 |
| 零升力迎角 α₀ | ≈ -2.5° | CL=0对应的迎角 |
| 线性段范围 | -2° ~ 12° | 失速前的线性区域 |
将两个斜率代入公式:
SM = -(-0.020) / 0.097 ≈ 0.21 = 21%
| SM范围 | 稳定性评价 | 操纵性评价 |
|---|---|---|
| 5% ~ 15% | 适中 | 良好 |
| 15% ~ 25% | 偏稳定 | 稍迟钝 |
| > 25% | 过于稳定 | 操纵困难 |
| < 5% | 不足 | 过于灵敏 |
| 负值 | 不稳定 | 危险 |
| 21%(本文) | 偏稳定 | 可接受 |
SM=21%在可接受范围内,飞机俯仰稳定。对于少年组的限时载运飞机来说,偏稳定是好事——飞起来更"听话",不容易失稳。如果后续试飞觉得操纵太迟钝,可以微调重心位置(把重心往后挪一点),将SM降到15%左右。
在3D视图中检查整机布局时,发现水平尾翼的位置似乎偏前了。根据3D模型测量,尾翼前缘到重心应为0.8m左右,但初始设置的尾翼x位置(x=0.600m)导致尾臂偏短。
将尾翼从x=0.600m后移到x=0.900m,修正后的关键参数:
修正后尾翼参数:
| 参数 | 修正前 | 修正后 |
|---|---|---|
| 尾翼x位置 | 0.600 m | 0.900 m |
| 尾翼前缘到CG | 0.51 m | 0.81 m |
| 尾臂(CG到尾翼AC) | 0.555 m | 0.855 m |
| 尾容系数 V_h | 0.34 | 0.52 |
尾容系数从0.34提升到0.52,在0.3~0.6的合理范围内,尾翼面积和力臂的搭配更加合理,能够提供充足的恢复力矩。
| 检查项 | 结果 | 评价 |
|---|---|---|
| Cm-Alpha斜率 | 负值(-0.020/°) | 俯仰稳定 ✓ |
| 静稳定裕度 SM | ≈ 21% | 偏稳定,可接受 |
| 尾容系数 V_h | 0.52 | 合理范围 (0.3~0.6) ✓ |
| 尾翼位置 | x=0.900m(已修正) | 尾臂0.855m ✓ |
| CL-Alpha线性段 | -2° ~ 12° | 工作范围充足 ✓ |
稳定性分析确认飞机"能自动恢复",但还有一个问题没解决——配平(Trim)。从Cm-Alpha曲线可以看到,在巡航迎角3°时Cm≈-0.09,这不是零,意味着飞机在巡航时有低头趋势,需要升降舵偏转来抵消。
下一步的工作:
系列文章导航:
系列一:设计步骤与20块载运图 |
系列二:双箱载运系统 |
系列三:比例优化 |
系列四:仿真分析全景图 |
系列五:从升力公式到雷诺数 |
系列六:NACA 4415翼型分析结果 |
系列七:3D机翼定义教程 |
系列八:定义分析 |
系列九:机翼定义修正 |
系列十:机翼气动分析结果 |
系列十一:水平尾翼定义与俯仰稳定性设置 |
系列十二:水平尾翼稳定性分析——Cm曲线验证俯仰稳定(本文)