降落伞悬挂双锥八面体角反射器的动力学分析

《Aerospace》：Dynamic Analysis of a Parachute-Suspended Bipyramidal Octahedral Corner Reflector

　　

摘要翻译：机载角反射器（ACR）是一种新型雷达无源干扰装置，由于其在降落伞悬挂下可增强干扰覆盖范围，已引起全球研究人员的广泛关注。然而，ACR的定向特性使其效能对飞行姿态动力学高度敏感。研究人员通过分别分析降落伞体和角反射器，提出了一种改进的动力学模型来描述降落伞-载荷系统。关键创新包括：（i）引入11自由度模型进行降落伞安装的双金字塔八面体结构的运动分析，克服了以往方法分析不精确的问题；（ii）明确建立系绳张力和几何约束的模型，以捕获降落伞-载荷耦合机制。稳定下降阶段的数值模拟证明了载荷角速度和欧拉角的收敛性，结果与全尺寸飞行试验数据吻合良好。该模型在计算效率和物理保真度之间取得了良好平衡，特别适用于降落伞下降系统中非轴对称载荷的动力学分析。

论文解读文章

**研究背景与问题**

角反射器是一种经典的雷达无源对抗装置，利用金属几何结构对入射电磁波进行后向反射，从而产生强烈的后向散射信号。传统角反射器主要部署于海面漂浮，近年来，机载角反射器（ACR）通过载体平台投放以实现高空干扰，成为重要的发展方向，如拉斐尔的WIZARD、莱茵金属的OCR、拉克鲁瓦的SEALEM和凯姆林的TORERO等系统。尽管这些系统设计各异，但均依赖降落伞辅助下降以维持驻留时间，并采用八面体反射器作为主要散射结构。现有研究主要关注ACR的电磁散射特性，如张等人研究了ACR的统计雷达散射截面（RCS）特征，吴等人研究了海上机载角反射器的动态多普勒特性，但均未充分解决降落伞下降与载荷姿态之间的动态耦合问题，而这对于实际干扰效果至关重要。以往的运动模型多基于六自由度（6-DOF）公式，可能无法完全捕捉降落伞与载荷之间的耦合运动。此外，传统模型将载荷简化为球体、圆柱体或质点，忽略了复杂非轴对称形状（如多面体角反射器）引起的气动弯矩、力矩-姿态耦合以及非对称气流分离导致的不稳定性。因此，有必要发展一种高保真且计算高效的动力学模型，以准确预测ACR在稳定下降阶段的姿态动力学。

**研究内容与结论**

为此，研究人员发展了一种针对半球形降落伞悬挂双锥八面体载荷的11自由度（11-DOF）受约束动力学模型。该模型将降落伞和载荷视为两个独立刚体，分别建立体固定坐标系，并引入冗余广义坐标向量，通过一个完整约束（系绳长度恒定）将自由度从12降至11。模型的关键创新包括：独立体固定坐标系实现精确解耦的运动学描述；基于运动学约束解析推导系绳张力，避免了经验弹簧-阻尼假设；针对八面体提出了分段解析气动力矩模型，考虑不同偏航角下的非对称投影几何。数值模拟针对稳定下降阶段进行，结果表明载荷的角速度和欧拉角收敛，且与全尺寸飞行试验数据吻合良好。通过与文献中最相似的六自由度模型对比，所提模型在偏航角、俯仰角的均方根误差（RMSE）、平均绝对误差（MAE）和最大绝对误差等指标上均表现更优。研究还发现，所提模型和对比模型均无法复现飞行试验中观察到的滚转角变化，这归因于风场微小变化、局部非对称流动以及模型本身的局限性。该模型在计算效率（在i9 CPU上模拟1000个时间点仅需2.55秒）和物理保真度之间取得了良好平衡，特别适用于非轴对称载荷在降落伞下降系统中的动力学分析。论文发表在《Aerospace》期刊。

**主要关键技术方法**

研究人员采用了以下主要技术方法：首先，基于牛顿第二定律和角动量定理分别建立降落伞和载荷的平动与转动动力学方程，并引入附加质量效应以考虑充气伞衣的方向依赖性；其次，利用运动学约束条件（系绳长度恒定）建立系绳张力模型，通过求解约束方程获得张力，从而耦合两个子系统；第三，针对双锥八面体载荷，根据偏航角范围（0≤φ<π/4和π/4≤β<π/2）分段建立气动力矩模型，利用对称性将迎风面划分为三个典型面（A、B、C），并积分得到总气动力矩。这些方法共同构成了完整的11自由度耦合动力学模型。数值模拟采用自开发的Python 3.10程序，基于四阶龙格-库塔方法进行求解，并与飞行试验数据（来自代表性飞行试验数据集，风速变化较小）进行对比验证。

**研究结果**

**3. 结果**：通过对比所提模型、对比模型（文献[3]中的六自由度模型）以及飞行试验实测数据，发现所提模型获得的偏航角、俯仰角与实测数据更接近，尤其是偏航和俯仰角。然而，两种模型模拟得到的滚转角始终为零，而实测滚转角在零附近有轻微波动，表明两种模型在捕捉滚转角响应方面仍存在局限。定量评估采用RMSE、MAE和最大绝对误差指标，所提模型在所有指标上均优于对比模型。

**4. 讨论**
**4.1 模型精度分析**：模拟与实测俯仰角在前18秒趋势相似，18至31秒模拟结果变化较快，31秒后变化率趋于一致。模拟俯仰角与实测结果合理吻合，但实测数据表现出更快的波动和更大的早期变化。模拟滚转角保持接近零，而实测滚转角在零附近轻微波动。总体而言，模拟结果与现场数据良好一致，偏航和俯仰角变化率相似，滚转角偏差较小。差异主要归因于飞行试验中风速变化。

**4.2 与经典降落伞-载荷模型对比**：所提11自由度模型弥补了简化低保真模型与高保真多体动力学框架之间的差距。与忽略姿态和相对运动的2自由度质点模型相比，该模型明确解析了降落伞和载荷的旋转动力学；与假设刚性连接的6自由度刚体模型相比，基于约束的系绳张力模型避免了经验弹簧-阻尼假设，产生物理一致的力传递。尽管超过12自由度的模型可考虑柔性系绳等更多因素，但通常增加计算复杂性和不确定参数，而该模型在保留关键物理（非轴对称空气动力学、几何约束）的同时维持实时可模拟性。

**4.3 适用性与局限性**：该模型主要适用于降落伞充气后稳定下降阶段，假设伞衣完全充气且忽略大变形，将降落伞和载荷视为刚体，系绳通过几何约束建模。局限性包括：未考虑伞衣充气过程和瞬态展开动力学；忽略伞衣柔性、织物振动和悬挂线弹性；气动系数视为给定参数，未显式解析湍流、尾流干扰和流动分离引起的非定常气动效应；验证基于有限数量的飞行试验记录。因此，该模型最适合于稳定下降阶段的初步轨迹预测、姿态趋势分析和工程评估，而非详细部署模拟或高保真流固耦合分析。

**5. 结论**（翻译）：本研究基于线性动量定理和角动量定理，建立了一个改进的动力学模型，用于描述携带双锥八面体载荷的半球形降落伞。为降落伞和载荷分别引入了独立的体固定坐标系，使其能够独立描述平移运动、旋转运动和姿态演化。此外，分析了两个子系统之间的系绳张力和约束关系，以表征降落伞-载荷系统的动态耦合机制。针对降落伞-载荷系统的稳定下降阶段进行了数值模拟。结果表明，载荷的偏航角速度从初始值逐渐减小并趋于稳定水平，而俯仰角速度呈现振荡收敛。俯仰角也表现出初始快速下降后振荡恢复的趋势。与现场数据对比表明，所提模型能够再现实际系统的主要动态特征。该模型在保持相对简单数学结构的同时，考虑了降落伞-载荷耦合和复杂载荷几何形状的影响。因此，它为机载角反射器和其他降落伞悬挂的非轴对称载荷的动力学建模、轨迹预测和姿态稳定性分析提供了有用的参考。未来工作可能侧重于纳入非定常气动效应、风扰动、柔性系绳动力学以及通过计算流体动力学或风洞试验获得的更详细气动系数，以进一步提高模型在复杂飞行条件下的预测精度。