斜流下小型无人机螺旋桨的对称性破缺：基于物理的叶素/动量分析

《Symmetry》：Symmetry Breaking of a Small UAV Propeller Under Oblique Inflow: A Physics-Based Blade Element/Momentum Analysis

　　

斜流（oblique inflow）通过引入跨转子盘的方位角相关速度分量，破坏了几何对称螺旋桨的名义旋转周期性。本研究建立了一个透明的降阶叶素/动量（blade element/momentum）框架，用以考察代表性两叶小型无人机螺旋桨由此产生的一周一次气动载荷。该公式体系结合了叶素空气动力学、基于动量的诱导速度、傅里叶分解（Fourier decomposition）以及中等斜流条件下的前进-后退叶片（advancing–retreating blade）诊断。研究特别关注区分直接源自假设谐波结构的量与提供独立气动信息的量。一阶解析展开表明，对于小面内速度，归一化一阶谐波载荷近似与面内速度比成比例，因而与J sin(α)成比例，这与数值模型一致。该分析还报告了绝对一阶谐波载荷幅值、不对称性的展向发展、交叉平面桨毂载荷分量、环量不平衡以及对气动模型假设的敏感性。已发表的轴向和非轴向螺旋桨数据被用于确定降阶趋势可被解释的范围。结果表明，推进比（advance ratio）和入流角控制建模的循环不对称性，而定量幅值仍对截面气动假设和诱导流表示敏感。因此，该框架旨在作为斜流对称性破缺的计算低成本筛选与解释工具，而非相分辨CFD（phase-resolved CFD）或实验测量的替代品。

研究背景方面，小型多旋翼飞行器很少在完美轴向入流下运行，前飞、姿态变化、机动和環境扰动会引入同时垂直于和平行于转子盘的速度分量，使叶片在每转中经历相对速度和有效迎角的周期性变化，导致盘平均推力可能保持平滑而单个叶片承受显著的一周一次载荷，产生俯仰/偏航力矩和偏斜尾流。低雷诺数螺旋桨性能强烈依赖运行点且难以用高雷诺数假设外推，已有研究表明叶素动量/涡格方法可再现斜流中小螺旋桨力的主要趋势但定量精度随非轴向性增加而下降，实验观测到推力及非零俯仰/偏航力矩随盘迎角增大而增大。这些发现促使开展对假设透明的降阶分析：高保真CFD计算昂贵，而降阶叶素/动量模型可快速评估受控运行矩阵并揭示一阶方位谐波的数学来源。本研究的贡献在于一个紧凑的面向对称性的诊断框架，将一阶方位载荷谐波与两叶转子的名义π周期状态相联系，区分归一化谐波含量与其绝对载荷幅值，将一阶不对称源项与其径向载荷贡献相关联，并量化气动模型假设对对称性破缺指标的影响；该框架用于高保真研究之前的快速解释和运行点筛选。

研究人员开展了基于物理的降阶叶素/动量分析，建立了对称性框架，定义了几何对称转子在均匀轴向入流下的名义π周期性以及斜流引入的外部对称性破缺，以第一方位载荷谐波为主要表征，并保留半转误差Eπ作为对称空间表示。研究采用两叶0.254 m直径、名义螺距0.114 m的基线转子，活动叶素从r/R = 0.28开始，弦长规律由三个几何锚点重构，局部几何螺距角遵循等螺距螺旋关系；并设置外侧30%半径施加10°后掠的保对称变体（SLE），以CL,SLE = CL cos2(Λ)表示后掠区法向流修正。运行矩阵在5000 RPM下，取J = 0.15和0.30与α = 0°、10°、20°、30°组合，形成十六个主要模型算例。结论表明，斜流即使在两叶转子保持几何对称时也会产生一周一次叶片载荷分量，一阶模型和数值矩阵均显示在中等等入射角下Λ1,T近似按J sin(α)缩放；归一化和绝对循环载荷必须区分；Eπ对于近正弦载荷历史并非独立气动指标而是与Λ1,T解析关联的对称空间表示；10°外侧后掠算例不能用于得出物理后掠无效的结论；外部文献比较支持预测非轴向趋势的方向但未提供几何匹配验证；该公式体系的实用有效范围为中等斜流入流（J = 0.15–0.30，α ≤ 30°），用于快速趋势筛选和对称性诊断，最强条件J = 0.30、α = 30°仍推荐用于未来相分辨CFD/PIV或测力传感器验证。其重要意义在于提供了一个计算低成本、假设透明的筛选与解释工具，明确了模型有效边界，并为后续高保真验证指明目标运行点。

研究人员为开展研究用到的主要关键技术方法包括：降阶叶素/动量公式体系，结合叶素空气动力学、基于动量的诱导速度迭代求解、傅里叶分解和前进-后退叶片诊断；采用准定常低雷诺数截面极曲线作为闭合，升力斜率a = 2π rad?1、零升迎角α0 = ?2°、CD0 = 0.08、k = 0.06，并施加|CL| ≤ 1.25限制；在轴向悬停参考条件下用两个标量修正因子kT = 1.2933和kQ = 1.8223进行轴向标定，目标系数CT = 0.1065和CQ = 0.01462来自独立执行的5000 RPM悬停瞬态CFD基线（SST k?ω湍流闭合、滑移网格旋转域，约12.0百万单元格参考解）；采用120 × 360径向-方位离散分辨率；并开展外部基准比较和模型形式敏感性分析。样本队列来源为已发表的低雷诺数螺旋桨测量数据（UIUC数据库）以及斜流螺旋桨实验与数值研究文献。

研究结果方面，数值收敛性分析表明所选120 × 360离散在最强算例下接近240 × 720数值参考，CT差异0.16%、CQ差异0.24%、Λ1,T差异0.09%、Eπ差异0.08%、IΓ差异0.28%、差异0.06%，最粗30 × 90网格最大差异约2.0%（IΓ），支持120 × 360作为实用数值分辨率。外部基准结果表明，轴向比较受限于尺度一致性，当前CT在J = 0.15时为0.0833、J = 0.30时为0.0541，落在小型固定螺距螺旋桨预期系数范围内，但几何、截面和雷诺数分布差异排除了逐点误差百分比的可辩护性；斜流入流基准中，基线平均CT在α = 0°至30°之间于J = 0.15时上升5.3%、J = 0.30时上升19.6%，与Serrano等和Theys等的趋势报告一致。方位载荷与一阶谐波增长方面，α = 0°时叶素模型方位均匀，Λ1,T和Eπ均为数值零；斜流引入近似正弦的一周一次分量。基线J = 0.15时Λ1,T从10°的0.046增至20°的0.090和30°的0.129；α = 30°时J从0.15增至0.30使归一化一阶谐波比从0.129增至0.322，但该相对增加部分被平均推力系数同步下降所放大，对应绝对无量纲一阶谐波幅值CT,1 = CTΛ1,T从约0.0113增至0.0208（约84%）。一阶缩放也在数值上得到恢复：J = 0.15时Λ1,T/sin α在α = 10°、20°、30°分别约为0.267、0.263、0.258，该比值的弱变化证实低入射角响应主要由施加的面内速度控制。半转对称性损失与桨毂力矩方面，基线Eπ在J = 0.15、α = 30°时达0.182，在J = 0.30、α = 30°时达0.444，与近正弦波形下Λ1,T的解析关系一致；相应合成交叉平面力矩系数分别为0.00254和0.00468。平均推力响应方面，基线CT在J = 0.15时从α = 0°的0.0833增至30°的0.0878（+5.3%），J = 0.30时从0.0541增至0.0647（+19.6%），该定性行为与Serrano等的实验趋势及基于入射角的动量解释一致；增大α降低轴向分量Va = V∞cos α使螺旋桨移向较低轴向推进条件，同时V‖ = V∞sin α产生前进-后退不对称性，降阶矩阵报告的是组合响应。环量不平衡与近尾流偏斜方面，前进-后退环量不平衡随α和J增大，基线α = 30°时IΓ在J = 0.15为0.0655、J = 0.30为0.1787（增加173%）；近尾流偏斜指标在x/D = 2处随J从0.15增至0.30由0.620增至1.280（α = 30°），这些值描述来自致动盘替代模型的运动学尾流偏斜位移估计，不应读作解析尾流质心或叶尖涡轨迹预测。简化后掠表示的结构敏感性方面，cos2Λ修正主要作为局部力尺度修改而不引入额外方位相关机制，归一化循环指标按构造预期几乎不变，SLE算例仅用于证明该降阶公式体系的结构性质，不用于得出物理前缘后掠在斜流下影响可忽略的结论。局部模型形式敏感性方面，对每个扰动模型重新标定至相同轴向CT和CQ目标后，最强非轴向算例对±10%升力斜率和±20% CD0变化相对不敏感；零升迎角影响最大，α0偏移?1°使Λ1,T降至约0.291，偏移+1°使其增至约0.367，Eπ在同一扰动范围内为0.403至0.502，证实绝对幅值依赖假设的低雷诺数截面极曲线，尽管对α和J的定性依赖在测试扰动内稳健。

讨论部分总结指出，降阶结果区分了常被混淆的两种效应：几何对称性是固体转子的属性，而气动对称性还受外部速度场控制；两叶转子可在180°几何旋转下保持精确不变，却因斜流面内分量不随叶片旋转而产生大的一周一次叶片载荷，Λ1,T和Eπ随α和J的强增长不是数值缺陷而是优选外部方向的预期数学后果。在当前一阶后掠表示内，运行条件效应远大于适度后掠平面形状的建模效应，后掠修正主要修改外侧叶片区域的局部力尺度因而产生平均和循环载荷的相似比例变化，归一化对称指标的近不变性应解释为当前降阶表示的结果而非后掠普遍无效的证据。模型趋势在两个方与已发表非轴向研究一致：固定推进比下平均推力随入射角增大，俯仰/偏航力矩和叶间载荷差异随前进-后退不对称性增强而增大；但结果应解释为趋势级证据，低阶方法在斜流中定量精度下降，故本文不声称高保真尾流预测或实验验证。局部模型形式敏感性提供了实用不确定性尺度：J = 0.30、α = 30°时零升迎角±1°变化使Λ1,T在 nominal 0.322附近从0.291变至0.367（约?10%至+14%），IΓ从0.159变至0.208（约?11%至+16%），运动学尾流偏斜估计变化小得多（约±3%）；此外准定常假设在约化频率k ≈ 0.12时可能高估循环幅值约15%，均匀诱导不再现纵向/横向入流变化相关的相移，这些范围比单一 nominal 值更有用地定义了模型的实际有效边界：该公式体系适用于J = 0.15–0.30和α ≤ 30°的趋势筛选，但不适用于相精确载荷、高入射角流动、解析尾流拓扑或认证级预测。半转量Eπ被保留是因为它以直接与对称性分析相关形式表达名义π周期叶片状态的损失，对于当前模型生成的近正弦历史它解析依赖Λ1,T故不计为独立气动证据，其价值在诊断：偏离Λ1,T–Eπ解析关系将揭示高次谐波或非正弦载荷。同样，IΓ源自建模截面环量，不作为解析尾流可视化的替代。最强建模条件J = 0.30、α = 30°是后续高保真验证活动的自然候选，有针对性的URANS或实验活动可在扩展矩阵前先在该单一运行点测量相分辨推力/力矩历史和近尾流速度。局限性方面，本研究是基于物理的降阶分析而非实验或CFD模拟，叶片几何由重构弦长规律和等螺距螺旋角表示而非完整CAD表面，截面气动采用准定常截断线性升力律和二次阻力极曲线，代表性雷诺数范围低、零升迎角影响大且无精确叶片截面的实测极曲线，未显式包含Prandtl叶尖损失因子，诱导速度在盘上均匀故标准纵向/横向入流变化和动态入流未被解析，运动学尾流偏斜位移在代数外推中使用远尾流动量速度而非解析涡轨迹或尾流质心预测，一阶后掠修正仅表示外侧展向后掠的宏观法向流效应，转捩、雷诺数相关截面极曲线、叶片柔性、电机/机体干扰、旋翼-旋翼干扰、大气湍流和声辐射均超出当前范围，因此结果支持相对趋势、有界对称性诊断和假设生成而非认证级性能预测。研究结论翻译如下：斜流即使两叶转子保持几何对称也会产生一周一次叶片载荷分量，一阶模型和数值矩阵均显示中等人射角下Λ1,T近似按J sin α缩放；归一化和绝对循环载荷必须区分，α = 30°时J从0.15增至0.30使Λ1,T从0.129升至0.322而绝对无量纲一阶谐波幅值从约0.0113升至0.0208（约84%），先前引用的149%归一化增加因此包含平均推力的同步下降；Eπ对于此处产生的近正弦载荷历史不是独立气动指标而是与Λ1,T解析关联的对称空间表示，保留作为π周期性和潜在高次谐波偏离的诊断；10°外侧展向后掠算例不能用于得出物理后掠无效的结论，在当前cos2Λ型力修正下平均和循环载荷被相似缩放，故低于0.5%的归一化差异很大程度上是降阶表示的结构性后果，三维后掠效应需要CFD或实验；外部文献比较支持预测非轴向趋势的方向——推力与交叉平面载荷随入射角上升及尾流组织改变——但未提供Λ1,T、Eπ或尾流偏斜估计的几何匹配验证，局部模型形式敏感性表明最强算例下±1°零升迎角扰动带来Λ1,T约?10%至+14%不确定性，另有k ≈ 0.12时约15%量级的准定常幅值不确定性；当前公式体系的实用有效范围为中等斜流入流（J = 0.15–0.30，α ≤ 30°），用于快速趋势筛选和对称性诊断，最强条件J = 0.30、α = 30°仍推荐作为未来相分辨CFD/PIV或测力传感器验证的要点。