基于代理模型全局优化的自适应参数化方法

对于翼型气动隐身设计问题，设计变量的配置对设计结果影响很大，而简单地增加设计变量不能保证得到理想的结果。提出一种适用于代理模型全局优化的自适应参数化方法：利用全局敏感性分析方法——基本效应法，得到设计空间关于目标函数的敏感区域信息，并以此为根据增加设计变量；利用节点插入算法将低维样本在高维空间内进行重构，避免了重新取样的工作量。相对于传统固定设计空间维度方法，自适应参数化方法在设计空间的敏感区域扩展维度，能够更加精准地描述外形并反映目标的变化趋势。通过飞翼布局翼型的气动隐身优化算例，证实自适应参数化方法可以大幅提高优化设计质量和效率。     

最小数据丢失量的地月中继卫星任务调度研究

月球背面的探测器必须依靠地月中继卫星进行数据传输与通信。地月中继任务包括实时性任务和延迟容忍类任务，如数传任务。当探测器等待传输的数据量超出用户存储容量时，延迟容忍类任务会由于探测器本地存储资源不足和地月中继卫星天线资源受限而无法完成，导致任务数据丢失，所以需要设计一种合理的地月中继任务调度策略，提高地月中继卫星的资源利用率，减少数据的丢失。对地月中继卫星任务调度进行了研究，在分析地月中继卫星数传任务的特点及用户的存储限制的基础上，以最小化数据丢失量为优化目标，建立了地月中继卫星任务调度模型，并设计了一种基于离散烟花算法（DFWA）的地月中继卫星任务调度算法。仿真数据分析表明，基于离散烟花算法的地月中继卫星任务调度算法在求解结果上优于遗传算法，是一种合理、有效的调度方法。     

不确定条件下航空不安全事件灵敏度分析的Monte-Carlo方法

为解决不确定条件下航空不安全事件灵敏度分析的难题，基于Bow-tie模型提出了多模式下航空安全性指标及其灵敏度测度，以轮胎爆破事件为例，采用Monte-Carlo方法计算得到安全性指标、基本事件的全局灵敏度及其分布参数的局部灵敏度。根据轮胎爆破事件仿真结果，两类灵敏度指标均随着飞行时间的增加而发生变化，尤其是在500~600 h时变化最为显著，但灵敏度重要性排序保持不变；基本事件类型是影响灵敏度的一个主要因素，电子类基本事件灵敏度测度远远小于机械类基本事件；同类型基本事件中，平均故障间隔时间不是主导因素，灵敏度大小还与失效传递的逻辑关系密切相关。算例结果表明:航空安全水平随着飞行时间逐步下降，应重点关注飞行时间为500~600 h时航空组件失效导致事故发生的可能性；航空组件的灵敏度重要性不会随着飞行时间变化，提高灵敏度较高的基本事件的可靠性水平是防范航空事故的关键。      

数值模拟方法对NASA CRM模型阻力预测的影响

在复杂工程外形的数值模拟中，网格类型、规模和分布、湍流模型、数值格式等都会不同程度影响模拟结果。如何评估这些模拟方法对模拟结果的影响，并识别对模拟结果影响较显著的因素，对于CFD的发展方向有积极的借鉴指导意义。为了综合研究不同因素对商用运输机外形阻力预测的影响，以AIAA第六届阻力预测会议外形NASA CRM为研究对象，同时考虑了网格、湍流模型、无黏通量格式和体心梯度求解方法等因素对阻力预测的影响。分析中采用枚举法和正交试验设计两种方法，并使用了基于聚类分析的定性敏感性分析方法和基于Mckay主影响分析的定量方法，识别出对阻力预测影响较大的因素，这为数值模拟方法的发展指明了方向。     

液氮贮箱自增压实验与数值仿真

为研究低温推进剂在常温下的自增压过程，设计了以液氮为模拟介质可视化低温玻璃贮箱自增压实验系统，研究了自增压过程压力和温度的变化规律及体积充填率对压力和温度变化的影响。实验结果表明：气枕区和液体区存在显著的轴向温度分层，液体区温度的上升速率低于压力引起饱和温度的上升速率。压力上升分为有典型意义的三段：初始段、过渡段和稳定段，稳定段的压力上升速率随体积充填率增加而增加。液体区的对流运动在自增压过程受到抑制，气液界面逐渐进入准静止状态。并以实验测得温度作为边界条件，采用流体体积（VOF）模型对整个自增压过程进行了175 s的数值仿真。仿真得到的压力曲线变化规律与实验结果基本一致，稳定段的压力上升速率是实验值的1.58倍。本文得到的自增压物理参数变化规律，为低温推进剂的贮存和贮箱的热防护设计提供参考。     

光敏推进剂激光烧蚀过程建模与仿真

相较于传统大卫星,微小卫星具有结构紧凑、质量轻便和成本低廉的特点。然而,受功率和质量负载的限制,微小卫星一般不装备推进系统,其航线也局限于近地轨道。为扩展微小卫星的功能,满足日益复杂的任务需求,需给其配备合适的微推进系统。固体推进系统具有结构简单、寿命长、可靠性高的优点,但无法重复启动。为得到可重复启动的固体微推进系统,设计了一种非自持燃烧的光敏推进剂,采用激光控制其燃烧。在背压为大气压的环境下,利用高速摄像机拍摄燃烧过程并记录燃速。之后,对光敏推进剂的激光烧蚀过程进行建模。分析结果表明:激光可控制光敏推进剂的燃烧,燃速与激光强度成线性关系;该光敏推进剂的最小激光点火强度为0.28 W/mm~2;燃速计算值与实测值的误差在10%以内,证明该数学模型具备工程应用价值。     

时间序列法在临近空间大气风场预报中的应用

受多种因素影响，临近空间大气环境要素复杂多变，预报难度很大.本文采用时间序列法中的自回归滑动平均（ARMA）模型对临近空间大气风场开展统计预报方法研究，基于廊坊（39.4°N，116.7°W）中频雷达在88km高度的大气纬向风数据开展预报试验.本次预报试验的样本数据为2015年9月24日至10月24日风场数据，利用过去7天数据对未来第8天风场数据进行预报.试验结果显示，ARMA模型对临近空间大气风场预报有一定的适用性.当风场变化规律性较强，即样本数据风场呈现出比较显著的24h周期性变化时，ARMA模型预报效果较好；当风场发生突变时，预报效果变差.与实测数据的对比结果表明，ARMA模型预报结果的误差在9~27m·s-1，预报效果优于同阶自回归（AR）模型，略优于高阶AR模型.      

热防护系统损伤容限分析

建立了含初始矩形损伤的热防护系统(Thermal protection system,TPS)气动热分析的CFD数值模型,分析结果表明损伤区域侧壁出现了很高的热流密度峰值,并且迎风面侧壁峰值高于背风面,而损伤区域底部热流密度却很低。利用分析获得的热流密度建立了含损伤和无损伤TPS的有限元传热分析模型。分析结果表明:损伤的存在导致防热瓦最高温度急剧上升,超过其材料能承受的极限温度(1 500℃),防热瓦首先失效,而损伤对机体最高温度影响较小。最后进行了TPS损伤容限分析,在防热瓦极限温度约束下,外部热流密度最大值从100kW/m2增加到140kW/m2,矩形损伤宽度最大容许值从22.7mm减小到12.6mm,而弧形损伤宽度最大容许值从34.6mm减小到25.1mm,即随着外部热流密度最大值增加,损伤宽度的最大容许值降低,并且相同外部热流密度水平下弧形损伤宽度的最大容许值大于矩形损伤。     

直升机惯性传感器结构的模态优化

提出了一种基于迭代抽样和径向基插值的自适应代理模型方法。这种自适应方法以减少仿真计算数量和提高代理模型自适应能力为目的,使用多岛遗传算法选择新增样本点并使新增样本点位于设计空间的稀疏区域,使得所有的样本点均匀分布于设计空间。标准误差用来判断代理模型的精度大小以决定是否对代理模型进行更新。这种自适应代理模型结合多岛遗传算法对直升机的惯性传感器结构模态进行优化。用拉丁超立方抽样方法选择10个样本点构建初始的代理模型,自适应代理模型的计算结果表明2%的误差条件下需要额外增加7个样本点。优化结果表明不同的权重系数对最优模态特性的影响很大,惯性传感器结构的一至六阶模态值更加远离直升机的激励频率。     

华南地区台风全局化结构及风场特征实测研究

基于探空气球、雷达风廓线仪等设备长期观测数据,对影响华南地区热带气旋的全局化结构及风场特征进行了实测研究.以典型台风为例,考察了系统水平及垂直主结构特征,揭示了台风双眼壁结构、眼壁置换现象以及登陆后结构的轴非对称变化特征.提出了台风气压场水平轴非对称模型和垂直剖线模型,分析了台风气压场两关键参数(最大风速半径及Holland-B)数值分布特征及其与其他参数间的内在关系.采用复合分析技术,对不同来流状况下台风边界层风场特征进行研究,考察了风速剖线低空急流和风向随高度偏转的特征及其随径向距离和来流地貌的变化规律,建立了以低空急流为特征的台风垂直剖线经验模型,为精确化评估高层建筑风荷载提供数据和技术支撑.相关研究成果有助于深化对华南地区热带气旋全局化结构及风场特征的认识,并为合理评估该地区台风风灾害及构筑物台风风效应提供参考.