一种中心型间断有限元MMALE方法

针对多介质可压缩流体动力学问题,提出了一种单元中心型二维MMALE(Multi-Material Arbitrary Lagrangian-Eulerian)方法。在拉氏步,流体力学方程组采用中心型间断有限元方法求解。对于混合网格,采用Tipton压力松弛模型更新物理量,用等参坐标法更新物质中心点坐标。界面重构采用一种健壮的MOF(Moment of Fluid)方法。在重映步提出了基于多边形相交的二阶积分守恒重映方法。该方法分为四个部分:多项式重构、多边形相交、积分和后验校正。多边形相交使用"剪裁投影"算法,显著降低了多边形相交算法的复杂度。后验校正是基于MOOD (Multi-dimensional Optimal Order Detection)限制策略,并做了一些改动以适应多介质的计算。数值算例表明,该方法具有二阶的精度和较好的鲁棒性。     

潜入式喷管喉衬界面间隙优化设计

为了解决固体火箭发动机潜入式喷管喉衬在热载荷与内压联合作用下界面间隙优化问题,数值模拟方法求解了喷管热结构响应。采用流体仿真软件,计算了潜入式喷管流动过程的稳态流场,获取了高温燃气流动参数与固体壁面对流换热系数。并采用三维有限元结构计算平台,编写了非均布壁面压力载荷与非均布对流换热系数子程序,求解了喉衬前后搭接界面间隙在0mm,0.05mm,0.10mm,0.15mm,0.20mm下喷管热结构耦合问题,获取了材料内部温度场与应力场分布。结果表明,喉衬环向压应力与拉应力都随时间增加,先增大后减小。其次,随间隙增大,喉衬拉应力先减小后增大,喉衬压应力先增大后减小。再次,依据界面闭合与喉衬环向受力最小的准则,确定了喉衬前后间隙的相对最优值,前间隙相对最优值0.10mm,后间隙相对最优值0.05mm。本文数值方法可为喷管热防护材料界面搭接与喉衬界面间隙设计提供指导,也可应用于评估喷管热防护与结构强度安全余量。     

无隔道超声速进气道附面层排除特性飞行试验研究

为了分析评估某型歼击机无隔道进气道附面层的排除特性,设计搭建鼓包表面附面层压力梯度测量试验系统,进行了不同飞行高度、马赫数和姿态角等工况下的飞行试验。通过对飞行试验数据的整理、计算和对比分析同型号的缩比模型风洞试验结果,研究了无隔道进气道鼓包表面附面层排除特性。研究结果表明：稳定平飞时,在亚音速范围内,随着飞行高度的增加,鼓包构型对附面层的排除效果增大,而在超音速范围内,变化规律相反;在接近马赫数1.8及以上飞行工况下,鼓包表面附面层的扫除能力有所减弱,附面层气流分离加速,进而会造成较大的进气压力损失和畸变。单纯迎角飞行有利于增强附面层的排除能力;而带侧滑角飞行时,附面层压力系数曲线的拐点沿鼓包中心线平行向“背风面”偏移,偏移量与侧滑角成正比,进气道鼓包表面“迎风面”附面层排除能力增大,而“背风面”受气流分离影响而减弱。     

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

基于多代理模型的航空发动机建模及优化方法

为了进一步提高航空发动机建模及优化方法的性能,本文提出了一种基于多代理模型技术的建模及优化方法。本文首先提出了一种新的代理模型全局误差估计方法,以此建立了新的多代理模型建模方法。然后提出了一种组合模型预测偏差估计方法,以此发展了一种基于多代理模型技术的优化方法。6个不同维度及不同训练集大小的解析测试算例的结果表明,本文所发展的建模方法相较于现有方法精度更高,本文所发展的优化方法相较于经典代理模型优化方法算法收敛性更强。同时变循环发动机稳态性能建模及加速燃油控制规律优化实例表明,本文所发展的方法在处理实际工程问题时,依旧可以表现出良好的算法性能。     

民用飞机第23.867条款要求及适航审定研究

随着越来越多先进复合材料应用在民用飞机上以及民用飞机上使用的电子电气设备重要度和集成度越来越高,飞机的闪电防护设计越来越重要。CCAR-23第23.867条对民用飞机机体结构闪电直接效应防护进行了规定。研究了该条款的制定和修订的背景和历程,明确了条款的实质要求。23部飞机的使用用途和设计特点使得23部飞机的闪电防护设计和符合性验证也有自己的特点,为此,重点研究了23.867条款的基本设计、审定步骤和审查关注要素等,梳理出了飞机机体结构的闪电直接效应防护设计和适航审定的基本流程,给出了可接受的符合性方法,提出了23部飞机闪电直接效应防护设计、审定各流程的主要工作、可指导或可参考的文件、重点关注要素和技术审查要点,为民用飞机的闪电直接效应防护的设计和适航审定提供了参考。     

基于离散反电动势估计的永磁同步电机无速度传感器控制

针对高速永磁同步电机(PMSM)的无速度传感器控制问题,提出了一种基于离散反电动势估计的PMSM无速度传感器控制策略。实施离散反电动势估计的设计有三点:首先,设计了离散dq电流观测器以消除反电动势估计中的电感交叉耦合效应;然后,设计了延迟补偿模型以补偿数模转换引起的电压误差,同时设计了较为精确的离散模型,以克服由数字实现导致的估计反电动势和实际反电动势之间的偏差;最后,开展了高速PMSM驱动试验,试验结果验证了所提出的高速PMSM无速度传感器控制方案的性能。     

基于改进Taguchi方法的航空发动机装配成功率计算方法

考虑到运用传统的Taguchi方法在进行公差分析和装配成功率计算时,在计算精度上的不足,无法满足航空发动机研制过程中公差设计高精度、高效率、短周期、多组成环等要求。针对该问题,文中提出一种基于改进的Taguchi方法的装配成功率计算模型。利用改进的Taguchi方法对统计数据的前四阶矩进行计算,并结合Pearson分布理论建立装配成功率的计算模型。在理论分析的基础上,以某型航空发动机涡轮转子的装配成功率为例进行计算。结果表明:与蒙特卡洛仿真结果对比,改进的Taguchi方法在对装配间隙的峰度进行估计时使相对误差由原来的19.38%降低到0.21%;用本文所提方法得出的装配成功率与蒙特卡洛仿真结果的相对误差仅为0.02%,计算结果验证了该方法的有效性。     

固体火箭发动机零维两相燃烧室压强计算方法研究

为了更准确地预估含金属燃料固体火箭发动机的燃烧室压强,在压强计算中考虑两相流的影响,从一维两相喷管流动的求解出发,通过两相平衡流模型、两相常滞后模型、两相等温流模型、颗粒定温模型等模型的简化,分别推导不同模型下喷管中两相混合物的流量计算公式,再把流量公式应用到发动机零维内弹道理论中,推导并简化得到零维燃烧室平衡压强的计算公式。把压强公式用于HTPB推进剂固体火箭发动机和铝冰固体火箭发动机的燃烧室压强计算,结果表明,当固体推进剂中金属含量较高时(如铝含量为21%的HTPB推进剂发动机),用传统零维燃烧室压强公式预估的压强与实验误差较大,而使用合适的两相流模型和对应的零维燃烧室压强计算方法,在HTPB发动机中,能把压强预估结果与实验的误差降低到6%以内。如果使用多维内流场计算的方法,燃烧室压强预测结果的误差将下降到2.5%以内。结论发现在含金属固体火箭发动机的燃烧室压强计算中,考虑两相流的影响是必要的,而使用两相流修正后的零维燃烧室压强计算公式能够快速、较准确地预估这些发动机的燃烧室压强。     

多维对流通量重构方法的改进

为了研究多维对流通量重构格式中不同形式主方向对激波稳定性、压力场计算方面的影响,分别采用以笛卡尔方向、密度梯度及压力梯度作为主方向进行计算,对结果分析后发现,单纯依赖物理量梯度的主方向会造成压力等值线曲折的非物理解。为了消除这种非物理现象,采用一种基于压力函数的混合方式来确定主方向。单层及展向多层高超声速圆柱绕流问题模拟分析的结果表明:混合方式所确定的主方向不仅可以消除笛卡尔坐标系确定的主方向计算带来的红玉现象,同时还可以消除单纯的密度或压力梯度作为主方向所带来的压力场"zigzag"的非物理现象;通过与结构网格结果的比较表明:网格相关格式驻点处热流的相对误差为29.2%,而混合主方向确定方法所计算驻点热流相对误差只为4%,表明此改进方法计算得到的壁面热流是合理的。