基于六自由度方程的飞机结冰问题仿真

基于六自由度方程对飞机的结冰飞行动力学问题进行了仿真研究。首先针对飞机遭遇结冰时的结冰飞行过程进行仿真,开发了一种结冰过程中飞机结冰严重程度及气动参数的计算方法。通过在各仿真时刻对飞机气动参数进行结冰影响修正以实现对飞机结冰飞行过程的模拟,同时分析了结冰过程对飞行特性的影响。进而对飞机结冰后带冰飞行时不同结冰严重程度下的操纵面阶跃响应进行了仿真,并对仿真结果进行了分析。     

卫星姿态的鲁棒反馈滤波器方法研究

针对传统卫星姿态估计过程中计算量大的问题,提出了基于控制反馈的卫星姿态鲁棒滤波方法,将滤波器的新息作为闭环系统的输入,由此可以分离动态滤波方程中的时变因素(卫星姿态角速度),再将时变因素作为系统的不确定性考虑,最终通过鲁棒滤波方法对线性时不变系统的最终稳态值进行求解。该方法可以通过陀螺与星敏感器的参数信息,提前离线计算出滤波器的相关参数,具有很强的工程可用性。仿真结果表明,在卫星存在旋转角速度时,三轴的估计精度(MSE)比不考虑系统时变因素的控制反馈滤波器分别提升了12%,21. 05%,4%。     

基于事件触发的异构多智能体最优包含控制

针对异构多智能体系统的输出包含控制问题，提出一种基于边的事件触发最优控制协议，保证所有跟随者的输出能进入到由领航者的输出所形成的凸包中。同时，使系统达到最优性能，最小化控制代价。考虑到不是所有的跟随者都可获得领航者的信息，提出一种基于边的分布式事件触发观测器，估计领航者输出形成的凸包内点的轨迹。设计加权代价函数评价包含控制的性能，并将输出包含问题转化为最优状态反馈控制设计问题。利用贝尔曼方程和黎卡提方程，给出异构多智能体系统最优输出包含控制的参数设计。选择不同类型的机器人构成多智能体系统，验证算法的有效性。     

航空发动机燃烧室数字孪生体系关键技术

通过曲线坐标系浸没边界方法（IBM），在设计阶段可以实现对真实航空发动机燃烧室的高保真虚实映射，保留全部几何结构信息。在IBM方法基础上，采用大涡模拟（LES）结合概率密度函数输运方程湍流燃烧模型（TPDF），对双旋流燃烧室、某单头部直流燃烧室以及某折流燃烧室1/10模型进行真实结构仿真，测试数字孪生体系关键技术的有效性。对比预测结果和实验结果，双旋流燃烧室的旋流器出口附近轴向、径向、切向速度平均误差分别为15.7%、23.8%、15.0%；非稳态解析了真实直流燃烧室与折流燃烧室的详细湍流燃烧场，出口温度分布的平均相对误差分别为11.66%和17.95%。因此基于虚实映射得到的燃烧室数字孪生体系具有一定的有效性，该方法具有潜在的工程应用前景。     

基于网格自适应的飞行器防热材料热传导系数辨识

为准确辨识高超声速飞行器防热材料热传导系数，构造了针对热传导偏微分方程系统的约束泛函极值问题，基于表征多项式逼近性能的Weierstrass定理，采用Lagrange多项式对热传导系数进行参数化建模，进而将无穷维的约束泛函极值问题转化为有限维的非线性规划问题，再利用基于动态优化方程的优化方法将此非线性规划问题转化为常微分方程初值问题进行求解。为从较有限的测量数据中准确地“学习”热传导系数，建立采用“过拟合”判别准则的网格自适应迭代算法改善辨识精度。研究表明本文采用的辨识策略与优化方法有效，辨识结果可以较准确地反映材料热传导系数的变化规律。     

基于离散伴随方程的三维雷达散射截面几何敏感度计算

针对有限差分法计算雷达散射截面（RCS）梯度效率低，采用高精度雷达散射截面评估时计算代价高的问题，提出了一种基于麦克斯韦积分方程离散伴随方程的RCS梯度高效计算方法。基于伴随方程的梯度计算可以通过一次雷达散射截面求解、一次伴随方程求解获得RCS关于所有设计变量的梯度。其中麦克斯韦积分方程离散伴随方程的形式与原方程基本一致，可以采用矩量法（MOM）及多层快速多极子算法（MLFMA）求解。伴随方程求解计算量与直接雷达散射截面评估基本一致，存储量在直接雷达散射截面评估的基础上增加不明显。通过双椎体模型、导弹模型对基于矩量法、多层快速多极子算法的伴随梯度进行验证，证明了基于伴随方法的RCS梯度计算可以实现复杂外形中RCS梯度的高效、高精度求解，为基于梯度的高精度气动/隐身一体化优化提供了基础。     

随机恒速运动目标的搜索方程及持续探测概率

根据目标位置的初始分布和速度分布,给出了随机恒速运动目标在任意时刻位置和速度分布密度函数。定义了探测函数、目标存在和探测不成功的联合概率密度函数和目标存活概率函数,建立了基于目标速度分布的搜索方程并给出了其特征迹线解。在目标初始分布和速度分布均为圆正态分布条件下,分析了搜索方程特征迹线的基本特征和持续探测发现概率积分域的形态。给出了一个以直升机吊放声纳探测潜艇为背景的持续探测发现概率的算例,算例表明,对随机恒速运动目标持续探测的发现概率,与对静止目标进行持续探测的发现概率相似,服从时间的指数规律。     

基于VUMAT的三维裂纹动态扩展有限元模拟

文章基于ABAQUS/EXPLICIT平台的VUMAT接口,直接对材料积分点添加自定义本构方程和失效准则,开发了三维裂纹动态扩展子程序.在网格适当加密的情况下就能完成三维裂纹动态扩展模拟.在加载策略上选用适应性较强的增量法,直接在软件中通过增量步控制求解过程和收敛性,计算结果与试验结果一致.本文方法基于商业有限元软件ABAQUS,仅对材料本构方程和积分点失效进行编码,无需调整和重新划分网格,操作简单,计算效率高,适合推广到复杂的本构关系和大型工程模型.     

适用于气动弹性的小波非定常气动力降阶方法

发展非定常气动力模型降阶技术旨在缩减计算耗费并且使得计算流体力学信息能够应用于气动伺服弹性以及设计优化当中。采用小波方法建立基于Volterra级数的非定常气动力降阶模型。在模型识别过程中,激励是光滑连续的且以零阶保持方式进行离散,采样频率选择二的整数幂次。近似的一阶Volterra核基于Haar尺度函数族展开,鉴于Volterra核在系统响应中的衰减特性,可在合适的有限时长截断一阶核。为了获得一阶核的各项展开系数,需要求解由输入/输出数据组成的超定方程,其中涉及到奇异值分解算法。为了验证小波方法的有效性,算例选取了二维的NACA64a010翼型。数值仿真结果表明该方法能够比较准确地预测结构小扰动引起的非定常气动力响应且能描述一定的非线性现象。