考虑维修力量影响及载荷动态分配的k/n系统模糊可靠性分析

工程实际中,维修活动开展前往往存在一定时长的准备期,且由于环境时变性、系统长期运行后的劣化累积等因素导致部件状态性能水平存在不确定性,使得系统可靠性建模较为困难。对此,运用模糊数表征系统部件的失效转移率、修复转移率及修理工维修准备率的同时,以Power Law规则刻画部件间的故障相关关系,认为部件承担载荷超过某阈值时才会引发故障相关现象,并考虑了修理工数量与故障件数量之间关系对系统可靠性的影响,研究了载荷动态分配条件下带维修准备期的多修理工n中取k模糊多状态系统模型,建立了状态转移微分方程组,提出用逆向逐层分析的思路建立系统稳态概率系数的递推关系,应用α水平截集及Zadeh扩张原理确定了模糊状态概率的截集区间,得到了系统模糊稳态指标,最后通过算例给出了修理工数量及部件参数模糊程度对系统稳态指标的影响,验证了模型的适用性。     

剪切稀化非牛顿射流撞击液膜破碎直接数值模拟

液态射流撞击是液体火箭推进系统中广泛采用的一种燃料雾化方法,其破碎特征直接影响燃料最终的掺混及燃烧效率。采用直接数值模拟（DNS）工具,研究了低雷诺数（Re_l=41）和中等韦伯数（We_l=163）条件下剪切稀化非牛顿射流撞击液膜破碎的问题,着重分析了对角液膜的三维结构、破碎特征和非牛顿特性等。研究结果表明：在所研究的射流参数下,该非牛顿撞击液膜破碎属于Open Rim类别,破碎过程具有三维特性并伴随液丝与边缘的融合、液丝向液滴的转变等时域流动特征。液体的总表面积随时间不断增长,但单位表面积随液膜破碎的发生而下降,液膜扩张半角随时间逐渐增加并趋于恒定值43°,而后部液膜的长度几乎不随时间发生变化。此外,撞击液膜表现出明显的剪切稀化特性,液体内部最低黏性系数仅为零剪切黏性系数的1/5。     

基于位姿可达空间的太空机械臂容错路径规划

为了使太空机械臂在关节锁定故障后仍能继续完成后续任务,提出一种基于位姿可达空间的太空机械臂容错路径规划方法。基于牛顿-拉夫逊法计算太空机械臂关节人为限位,完成满足任务需求的退化工作空间求解,通过构造姿态可达度指标,在退化工作空间的基础上建立故障机械臂基坐标系下的位姿可达空间。通过在代价函数中增加最小奇异值代价项改进传统A^*算法,基于改进A^*算法在所建立的位姿可达空间内完成太空机械臂容错路径规划。所提方法综合了位姿可达空间与改进A^*算法各自的优势,实现了关节锁定故障太空机械臂同时满足避奇异与位姿可达要求的轨迹搜索。通过建立7自由度太空机械臂运动学模型开展数值仿真研究,仿真结果验证了所提容错路径规划方法的有效性。     

多电极偶接对金属大气腐蚀影响的试验与仿真

搭建薄液膜厚度测量与控制装置,采用微距参比电极后置法组建三电极体系,测量了厚度为100 μm的3.5wt% NaCl液膜下2A12铝合金、TC18钛合金和30CrMnSiNi2A高强钢的极化曲线和耦合体系的电偶电流;建立了基于薄壳电流分布的Comsol腐蚀仿真模型,得到了100 μm液膜下电极电位分布和电流密度分布图,通过对电极表面进行局部电流面积分计算得到了其电偶电流,对不同电极进行了面积的参数化扫描,讨论分析了不同面积比对于电偶电流的影响。结果表明：利用仿真模型得到的电偶电流值与试验值吻合较好,在2A12/30CrMnSiNi2A/TC18多电极体系中2A12充当阳极,TC18/30CrMnSiNi2A充当阴极,且阴阳极极性不随面积比的变化而发生转变,电偶电流主要集中在30CrMnSiNi2A与2A12之间,电偶电流与TC18面积呈线性关系,其对数与2A12和30CrMnSiNi2A面积的对数之间呈线性关系。利用建立的仿真模型对模拟搭接件大气腐蚀进行了仿真研究,仿真发现垫圈、2A12与30CrMnSiNi2A偶接接缝处及2A12与TC18偶接铆钉处电流密度分布更为集中,易发生腐蚀,与暴晒试验结果相符,说明该仿真模型对于多电极偶接的腐蚀分布具有很好的预测作用。     

基于直升机舰面起降动态仿真的风限图计算

建立了基于直升机舰面起降动态仿真的风限图（WOD）计算方法,综合计入舰船尾流、直升机运动和驾驶员操控等多因素的作用提升风限图计算的准确度。基于耦合舰船非定常尾流的飞行动力学模型,发展了适于直升机多轴协同操控的驾驶员模型,建立了舰面起降轨迹的数学描述和生成方法,形成计入舰船尾流、直升机运动和驾驶员操控等多因素综合作用的直升机舰面起降动态仿真方法。在此基础上,总结风限图计算判据,建立风限图计算方法。计算结果表明,某些风况下直升机在舰面起飞和降落过程中受到舰船尾流的干扰远大于在甲板上方悬停时受到的作用。本文方法能够捕捉到不同起降点和起降方式导致舰船尾流时空变化的干扰,与传统计算方法相比,显著提升了风限图计算的准确度。     

内吹式襟翼环量控制翼型升力响应特性

传统尖尾缘翼型通过控制迎角,综合利用襟翼、缝翼来改变升力,升力对迎角变化的时间响应历程可以用Wagner函数来描述,而内吹式襟翼（IBF）主要通过控制分离来拓展最大升力,并在一定范围内通过调节射流强度改变驻点位置和环量来对升力进行有效控制,其升力随吹气动量变化的时间响应尺度是否与传统尖尾缘翼型相同还不是很清楚。本文主要研究内吹式襟翼升力响应过程,并将其与传统尖后缘翼型升力响应特性进行对比。首先通过某襟翼偏角为30°的双圆弧环量控制翼型对数值方法进行验证,再对某最大厚度为18%弦长的亚声速翼型内吹式襟翼定常吹气控制下的流场进行非定常数值模拟,并分析了其中的瞬态特征。结果表明内吹式襟翼环量控制翼型对激励响应的时间依赖特征与Wagner函数有很好的相互关系,并可以用该函数来描述。     

一种基于改进WLBF频谱的异构双基地前视SAR成像算法

异构平台组成的双基地前视合成孔径雷达（SAR）配置灵活,可实现多种构型下的前视成像,具有潜在的应用价值。然而特殊的几何构型使得其回波信号具有新的特性,给频谱模型的推导及成像算法的设计带来不便。针对这一问题,提出了一种基于改进WLBF （Weighted Loffeld's Bistatic Formula）频谱的双基地前视SAR成像算法。首先,引入时域去走动和切比雪夫降阶法,减小了传统频谱模型中相位历程的加权分割误差和二阶展开误差,在此基础上推导了回波信号的改进WLBF频谱表达式。然后,使用切比雪夫多项式将改进WLBF频谱展开并设计了双基地前视SAR成像算法。最后,对频谱模型和成像算法进行了仿真分析。结果表明所提改进WLBF频谱模型与传统频谱模型相比具有较高的精确度,设计的成像算法可满足异构双基地大前视系统的成像要求。     

跨声速层流翼型的混合反设计/优化设计方法

跨声速层流翼型设计须兼顾优良的超临界特性和自然层流特性,因而对设计方法提出了更高的要求。针对现有反设计方法和直接优化设计方法的不足,发展了一种适用于跨声速层流翼型的混合反设计/优化设计方法。该方法引入了基于经验的局部流场特征作为反设计目标,翼型性能指标作为直接优化设计目标,然后加权形成了混合反设计/优化设计总目标,并同时考虑了气动和几何约束。优化算法采用基于自适应并行加点技术的代理优化,流动数值模拟采用耦合基于线性稳定性理论的e^N转捩自动判定的雷诺平均Navier-Stokes（RANS）方程求解器。针对现代中短程民用客机需求,以NPU-LSC-72613翼型为基准,开展了层流翼型减阻的混合反设计/优化设计。分别将局部目标压力分布、总阻力作为反设计和直接优化设计目标,得到了较好的优化结果,验证了方法的有效性。经过2轮优化结果显示混合反设计/优化设计总目标显著下降。所设计翼型吸力面局部压力分布与目标压力分布基本一致,总阻力下降15.5%;吸力面和压力面层流范围均大于55%倍弦长,激波强度显著减弱,说明所设计翼型同时具有优良的超临界和层流特性。将所设计翼型配置到机翼上,通过三维数值模拟进行校验,结果显示所设计跨声速层流机翼升阻比提高了6.64%,在一定升力系数范围内,气动性能均有显著提高,验证了所设计跨声速层流翼型在机翼设计中的适用性。     

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

近地卫星严格回归轨道保持控制

研究了近地卫星基于严格回归参考轨道的轨道保持控制方法：将卫星编队理论引入单星绝对轨道保持控制,提出了"虚拟卫星编队"的概念,分析了卫星轨道相对于参考空间轨迹在轨道摄动情况下的偏离状态及变化趋势,然后根据卫星编队相对运动学,推导出了偏离状态与虚拟卫星编队构形参数之间的对应关系,并设计了以轨道参数超调、偏置及阈值触发为特征的管道保持控制策略。数值仿真结果表明,使用该控制策略能够将卫星轨道保持在以空间参考轨迹为中心的轨道管道内,并且有效减少了因周期性轨道摄动波动造成的管道保持控制量和控制频次。研究成果对于有空间轨迹回归要求的卫星轨道保持控制具有指导意义。