可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

AIRS观测资料研究全球平流层重力波特性

利用2012—2014年1月和7月AIRS（Atmospheric Infrared Sounder）第79通道的观测数据，分析了平流层重力波活动强弱的全球分布以及重力波发生频率的全球分布；分析了重力波活动随纬度和经度的变化特征，给出了重力波活动在全球范围内的热点区域及其活动强度；对比了白天与夜间的重力波活动强度及发生频率.研究表明重力波活动强度呈现出随纬度变化的特征，在低纬度地区（0°—30°），冬季半球重力波活动强度低，夏季半球重力波活动高；在中高纬度地区，冬季半球重力波活动强度高而夏季半球重力波活动强度低.在1月，全球重力波活动有4个突出的热点区域，分别为50°N附近欧洲大陆与大西洋交接地带、北美洲与大西洋交接地，20°S附近南美洲与大西洋交接地区、非洲与印度洋交接地区.在7月，重力波活动突出的地方为巴塔哥尼亚至南极半岛地区，以及50°S和75°E附近的印度洋区域.重力波活动强度在夜间大于白天，但是夜间的强重力波活动区域小于白天.      

考虑燃料消耗均衡性的低轨通信星座在轨星座构型重构方法研究

针对低轨（low earth orbit,LEO）通信星座在轨重构问题，分析了低轨通信星座的特点，提出了4个重构指标，即全球平均覆盖率、燃料消耗均衡性、重构总时间和重构总速度增量，并分析了提高轨道高度的重构方式。采用基于分解的多目标进化算法（multiobjective evolutionary algorithm based on decomposition，MOEA/D）对重构进行建模分析。结果表明，模型优化结果可以得出多组最优解，构成帕累托前沿，并得出燃料消耗均衡性最优的解，针对两种失效模式得出的重构总时间均为3.6×105s，速度增量方差分别为261.4m2/s2和293.4m2/s2，这些解均可恢复星座原有的覆盖性能，并使燃料消耗最为均衡。     

压电周期板中耦合禁带影响规律分析

近年来，利用周期结构中弹性波禁带特性减振的研究思路受到了广泛关注。针对以往周期结构难以实现宽频且可调禁带的问题，设计了一种含电路网络的压电周期结构。该结构中的弯曲波和电路网络中的电波通过压电效应可以产生一个较宽的耦合禁带，且通过调整电路参数就可以达到调节禁带位置的目的。首先，为了高效计算该结构的波动特性，开发了适用于压电周期结构的减缩波有限元算法，该算法可以在保证结果准确性的基础上节约90%以上的计算时间。利用该算法研究了压电材料尺寸和形状对耦合禁带性能的影响。结果表明:相同电学参数下，随着压电片尺寸的增大，耦合禁带向低频移动，且禁带带宽增加；相同面积下含圆形和方形压电片的机电系统内耦合禁带差异较小，即形状对耦合禁带的分布影响不大。其次，针对不同尺寸和形状的机电系统，设计了电学参数使得在同一频率附近产生耦合禁带，并分析了其性能差异。最后，为了证明耦合禁带的减振效果，设计了一种有限压电周期板模型，其强迫响应的结果证明了耦合禁带对结构内弹性波可以进行有效调控。      

基于DOA聚类算法的ARM抗有源诱偏技术研究

有源诱偏技术是雷达对抗反辐射导弹的一种重要手段,可以大大降低反辐射导弹的作战效能。为提高对敌目标的精确打击能力,在分析闪烁诱饵诱偏原理的基础上,针对有源诱偏干扰下被动雷达测角精度与稳定度不高的问题,通过对有源诱偏信号的时域特征进行分析,提出一种基于脉冲前沿检测的DOA聚类分选算法,找出前沿超前的辐射源信号,实现了高性能抗有源诱偏干扰和高精度抗干扰测向。仿真结果表明,该技术可以有效对抗四点源有源诱偏干扰。     

转盘轴承集聚效应分析及其对回转性能的影响

针对航天发射台转盘轴承工作时出现滚动体集聚导致摩擦力矩增大的问题，研究转盘轴承中滚动体之间间隙变化对集聚效应的影响。通过对典型工位处进行有限元分析，对比不同接触类型转盘轴承的滚道变形和滚动体接触作用力的分布规律。在确定滚动体与隔离块的几何关系后，进一步分析了滚道变形和滚动体接触作用力的时变性对滚动体间隙变化的影响，并将间隙变化和集聚效应建立了联系。有限元计算得到的发射台框架变形与实验测量数据一致，验证了有限元计算的可靠性。仿真结果显示，由于下滚道发生轴向变形，滚动体在运动过程中不断交替进行“爬坡”和“下坡”；滚道径向变形导致滚动体之间增大了较多的间隙；各滚动体的驱动摩擦力不断变化，运动状态具有时变性，当转盘轴承中间隙增加时，运动状态不断变化的滚动体导致增加的间隙出现累积，引起滚动体集聚。为消除滚动体集聚效应，探索其对摩擦力矩的影响，进行了对比实验。结果表明:集聚效应导致发射台的摩擦力矩增加，且变化剧烈。      

多等级区域侦察弹性星座设计方法

针对传统侦察星座目标单一、弹性低的问题，提出了多等级区域侦察弹性星座的设计方法。该方法将星座设计过程按区域等级信息分为多个子星座逐步设计，直到整体星座对所有的区域性能满足设计要求。以区域被划分为3个等级为例，首先对星座设计需求、设计指标及设计步骤进行了分析。其次推导了地面最低分辨率和轨道高度的关系并确定了不同子星座的轨道高度。最后考虑轨道倾角、一箭多星发射、光照和升交点漂移同步约束，构建基础星座、子星座1和子星座2的优化模型。最终设计星座为3层混合星座，共8个轨道面和70颗卫星，星座对各等级区域的最大重访时间分别为10937s，12241s和17437s，弹性指数为2213%，2420%和6361%。结果表明该方法设计的星座可实现对区域覆盖和弹性分级的设计要求，证明了方法的有效性。对比Walker星座设计方法，在同等设计要求下，Walker星座所需卫星数为156颗，多等级区域侦察弹性星座所需卫星数远低于Walker星座，结果进一步证明了该星座设计方法的优越性。     

MIMO窄带加宽带非高斯随机振动试验

提出了一种用于多输入多输出窄带加宽带非高斯随机振动试验方法。首先，详细分析了随机相位与非高斯随机信号偏度和峭度的关系，提出了一种可用于快速生成具有指定偏度和峭度的迭代相位调节法。其次，根据控制通道之间的耦合关系，利用各通道之间随机相位的不相关特性，实现了多通道非高斯随机信号生成的解耦。然后，将生成的非高斯随机信号作为参考输出响应，采用时域逆系统方法生成驱动信号，并利用控制算法分别对响应信号的功率谱、偏度和峭度进行控制。最后，通过三轴向振动台试验验证了本文提出的多输入多输出窄带加宽带非高斯随机振动试验方法的可行性。结果表明，控制点上加速度的功率谱被控制在了±3 dB的容差限内，满足预设的容差要求。同时，时域响应信号的峭度和偏度也满足给定的参考要求。     

并行稀疏滤波在轴承声信号下的故障诊断

振动信号是航空发动机故障监测的常用信号。由于航空发动机结构复杂，对振动传感器的布置要求日益严格。声学信号以其非接触式、易布置、低成本的优点，在轴承智能故障诊断中引起了广泛的关注。然而，由于航空发动机内声信号所处的环境噪声较强，传统的轴承故障诊断方法无法实现精确的特征提取。为此，研究有效的特征提取方法实现轴承声信号下的智能故障诊断显得尤为重要。稀疏表示是智能故障诊断中的一个研究热点，在稀疏特征提取方面显示出强大的力量。对强噪声下的声信号进行有效的稀疏特征提取，可为轴承的非接触式故障诊断提供解决路径。提出一种基于并行稀疏滤波的轴承故障诊断方法，能够实现对轴承声信号的稀疏特征提取。并行稀疏滤波通过在传统稀疏滤波的基础上增加另一个归一化方向来实现进一步的稀疏特征提取，然后采用权值归一化方法约束训练得到的权值矩阵。最后，通过仿真和实验数据验证了所提方法的优越性。结果表明，并行稀疏滤波能够实现轴承声信号的有效稀疏特征提取和精准分类，可用于声学信号下的轴承智能故障诊断。     

基于改进冲突搜索的多智能体路径规划算法

多智能体路径规划问题在航空航天领域的多机任务中应用广泛但求解困难。基于改进冲突搜索的算法被设计用来快速求解多智能体路径规划问题。全局路径规划方面，首先设计综合考虑路径代价总和以及最大完工时间的多目标代价函数，其次提出基于唯一最短路径的冲突分类及消解方案，降低多智能体路径规划的计算量。在线冲突消解方面，利用速度障碍法在线检测和消解智能体与动态障碍物间的突发冲突。仿真结果表明，本文算法在全局路径规划方面保留基于冲突搜索算法的最优性并且降低了算法计算量，同时本文算法能够有效实现在线冲突检测与消解。