温度对磁屏蔽霍尔推力器磁场构型的影响研究

霍尔推力器磁路设计主要通过常温静态磁场仿真得到，并实测推力器非工作状态常温磁场进行复核。大功率霍尔推力器将面临更为严峻的热问题，推力器工作时磁路系统受高温影响，因此在常温下仿真得到的磁场位形会因温度升高而产生偏移，不能反映推力器真实工作时的磁场情况。为研究霍尔推力器工作时热量对磁路系统的影响，通过热磁耦合仿真对10kW磁屏蔽霍尔推力器的热态磁场分布进行研究，并对热态、常温仿真结果进行了对比，发现在阳极附近的径向磁感应强度Br的差异比放电室出口更大。常温设计的磁屏蔽构型在热态时偏离磁屏蔽，磁场和壁面最大不符合度达到13%，通过陶瓷出口型面修正后重新获得磁屏蔽效果，使最大不符合度降低到4.8%以下。合理热设计有助于降低热载荷，热仿真得到磁路系统最高温度低于500℃，低于0.78倍的居里温度Tc磁性急剧转变点，不会出现磁性能急剧下降，但热量对磁屏蔽霍尔推力器磁场构型的影响是应该考虑的。     

柔性薄壁大部件数字化装配调姿算法研究

采用定位器托架混合调姿定位平台,提出了一种飞机薄壁柔性大部件数字化调姿方法。首先,采用奇异值分解法(SVD)计算机翼位姿,通过运动学逆解得到定位器各轴的运动状态。其次,根据飞机大部件易变形的特点,提出了一种消除位姿计算误差的建立飞机部件坐标系的方法,实现快速高精度调姿。最后,将该方法应用于飞机机翼的调姿,经测量调姿精度满足装配要求。     

导弹配平飞行数值模拟及控制优化

传统基于建模仿真的控制系统检验方法受建模质量的影响,在非定常效应较强时很难得出准确结果。为此,采用基于CFD的虚拟飞行数值模拟方法检验和优化控制系统。首先介绍了数值模拟方法,进行了纵向开环控制飞行模拟,结果与气动力建模仿真结果、风洞试验结果吻合较好,共同验证了计算方法和气动力模型。针对迎角拉起设计了PID控制器,阶跃响应仿真结果表明控制器具有良好的动态性能,然后进行了数值模拟验证,结果表明迎角超调量超过50%、调节时间1.7s。通过试探法对控制参数进行了优化改进,将超调量减小到14%,调节时间减小到1.18s,动态性能显著提高。分析表明,受控时尾舵角速度达到200°/s,非定常效应强,线性气动力模型不能精确描述这些过程,造成了仿真结果的差异;虚拟飞行数值模拟技术在控制律检验和优化上可以起到关键作用。     

适用于非定常流模拟的分布式并行GMRES方法

为提高计算流体力学方法的收敛性和对高性能并行计算机的适应性,发展了适用于非定常流模拟的GMRES并行全隐式方法,并开展了相应的收敛和并行特性研究。采用变子空间数GMRES方法,减小重启过程计算时间;通过分区并行和Hybrid LU-SGS预处理算子实现方法的分布式并行化;采用鲁棒的Negative-SA湍流模型获得更大CFL数,采取计算和存储雅可比矩阵、网格重排序方法提高计算效率。利用这套方法完成了平面流、NACA0012翼型扰流、翼身组合体扰流、F-16战斗机非定常气动弹性和旋翼前飞流场的数值模拟。结果表明其计算效率较LU-SGS方法提高20%～200%;适用于当代高性能计算机分布式并行结构,并行效率非常高,在240个计算核心上出现了加速比的超线性。     

滑阀副零位内泄漏量分布模型与参数灵敏度分析

针对尺寸参数摄动作用下滑阀副零位泄漏量的不确定性，建立了考虑参数摄动的滑阀副零位泄漏量数学模型，获得了尺寸参数摄动下的滑阀副零位泄漏量分布规律。进一步地，建立了工艺参数与零位泄漏量分布特征参数之间映射关系的Kriging代理模型，并基于代理模型利用Sobol法开展了零位泄漏量分布特性对工艺参数的全局灵敏度分析，最终为减小内泄漏量并提高一致性提供了理论依据。研究结果表明，尺寸参数摄动下的零位泄漏量服从正态分布规律，工艺参数之间的交互作用对滑阀副零位泄漏量分布特性的影响可以忽略不计。同时，相比轴向尺寸工艺参数，滑阀副零位泄漏量分布特性对径向尺寸工艺参数和节流边圆角半径工艺参数的变化更为敏感。滑阀副径向间隙最小值对滑阀副零位泄漏量平均值影响最大，而阀芯的节流边圆角半径最大值对滑阀副零位泄漏量的一致性影响最大。滑阀副内泄漏量实验结果表明零位泄漏量呈现出一定的波动性，与理论模型总体预测趋势相一致。     

非结构CFD软件MPI+OpenMP混合并行及超大规模非定常并行计算的应用

常规工程应用中，非定常数值模拟（如多体分离）的计算量十分巨大，如果为了达到更高的计算精度，加密网格或者采用高精度方法将会使得计算量进一步增大，导致非定常数值模拟在CFD工程应用中成为十分耗时和昂贵的工作，因此，提高非定常数值模拟的可扩展性和计算效率十分必要。为充分发挥既有分布内存又有共享内存的多核处理器的性能和效率优势，对作者团队开发的非结构网格二阶精度有限体积CFD软件（HyperFLOW）进行了混合并行改造，在计算节点间采用MPI消息传递机制，在节点内采用OpenMP共享内存的MPI+OpenMP混合并行策略。首先分别实现了两种粒度（粗粒度和细粒度）的混合并行，并基于国产in-house集群采用CRM标模（约4 000万网格单元）定常湍流算例对两种混合并行模式进行了测试和比较。结果表明，粗粒度在进程数和分区数较少的小规模并行时具有效率优势，16线程时效率较高；而细粒度混合并行在大规模并行计算时具有优势，8线程时效率较高。其次，验证了混合并行在非定常计算情况下的可扩展性，采用机翼外挂物投放标模算例，分别生成3.6亿和28.8亿非结构重叠网格，采用对等的（P2P）网格读入模式和优化的重叠网格隐式装配策略，网格读入和重叠网格装配耗时仅需数十秒；采用3.6亿网格，完成了非定常状态效率测试及非定常分离过程的湍流流场计算，在in-house集群上12 288核并行效率达到90%（以768核为基准），在天河2号上12 288核并行效率达到70%（以384核为基准），数值模拟结果与试验结果符合良好。最后，在in-house集群上采用28.8亿非结构重叠网格进行了4.9万核的并行效率测试，结果显示，4.9万核并行效率达到55.3%（以4 096核为基准）。     

采用伪速度冲击响应谱评估蜂窝夹层板破坏边界

采用伪速度冲击响应谱(PVSRS)对航天器蜂窝夹层结构进行了冲击破坏边界评估。目前航天领域常用的冲击响应谱(SRS)是绝对加速度谱(AASRS)，相较于绝对加速度谱，伪速度谱能更清楚地显示出结构在低频、中频和高频段不同的冲击响应特征。分析了弹性力和惯性力在低频、中频和高频段对结构安全性的影响，据此在对数四坐标伪速度谱上给出了相应的结构破坏边界。以一典型蜂窝夹层结构为研究对象，通过有限元分析计算得到的结构破坏边界与伪速度谱给出的结构破坏边界符合得较好。在航天器结构冲击安全性研究领域，本文采用的伪速度谱结构冲击破坏边界是对目前常用的加速度谱方法的一个有益补充。     

机器人旋转超声钻削铝合金叠层构件毛刺特性

针对铝合金叠层构件机器人制孔容易产生钻削毛刺，严重影响飞机装配精度和效率等问题，提出基于弱刚度环境下机器人旋转超声制孔毛刺高度的计算方法。首先，通过实验验证了高频的振动冲击对毛刺高度的影响规律，建立了机器人旋转超声钻削铝合金叠层构件的钻削力经验公式。然后，结合经典薄板弯曲理论和能量法，分析了超声振动及钻削位置刚性对机器人钻削铝合金叠层板毛刺形成的影响机制。实验结果显示:所提计算方法的相对误差在13%以内，具有较高的精度。      

火箭发动机端面密封静环热压工艺及故障分析

静环是端面密封的主要组成部分,石墨环热压环节是整个静环生产的核心,热压成品率直接关系端面密封的质量可靠性。结合生产工艺过程,深入分析了热压工艺机理,建立了静环热压的数学模型,明确了镶嵌应力、过盈量校核等一系列参数。在此基础上,形成了一套完备的热压工艺方法。结合某型号试车端面密封故障问题,对静环进行了仿真分析和试验验证。结果表明,在热反浸过程端面密封受到高温作用,由于石墨材料线性膨胀系数发生变化,导致石墨环和静环座间的过盈量减小,过盈连接部位局部存在微缝隙,石墨环受高压介质挤压不均匀脱出,石墨环脱出不均匀导致石墨端面变形,产生泄漏。