高超声速飞行器末段轨迹快速优化

利用高斯伪谱法求解具有路径约束和落角、落速等终端约束条件下的高超声速飞行器末段轨迹优化问题.分析了不同LG节点数目对该问题求解精度和计算时间的影响,以及不同LG节点的拟合多项式与实际积分结果之间的相对误差关系.根据精度需要,选择合适的LG节点进行优化,分析了状态变量与积分结果的相对误差传播情况.仿真结果表明,采用节点分析后的高斯伪谱法在求解该问题上具有较高的计算精度和计算效率.     

高超声速飞行器结构/控制耦合综合研究

针对高超声速飞行器的结构/控制耦合问题,分析了结构/控制耦合效应关系.通过有限元分析方法对高超声速飞行器结构进行静力计算,结合计算结果,从结构弹性变形及弹性振动两方面出发,分析研究结构刚度特性对控制制导特性产生的影响,并对影响因素和耦合效应进行定量或定性的分析及评估,为高超声速飞行器的结构设计以及控制制导系统的设计提供参考.     

飞行器典型热密封结构

对航天飞机、X-38和X-51等飞行器的热密封结构进行了综合评述.分析了机身TPS、机身开口部位及控制面三个安装部位的热密封结构特点,综述了其热密封组件的研究历程和使用现状.     

跨声速风扇叶片反扭影响因素研究

基于结构几何非线性大变形的静态分析和流场分析,使用叶片反扭设计的流固双向耦合的数值模拟方法,得到NASA Rotor67跨声速风扇叶片的冷态加工叶型。研究了材料、气动工况、转速对叶片静态变形和反扭设计参数的影响。结果表明:转速对叶片反扭的影响最显著,气动工况次之,材料的影响最弱;另外,这三种因素和叶片反扭的关系,与其和叶片静变形量的关系有较大相关性。     

高超声速飞行器气动伺服弹性的自适应抑制

针对弹性高超声速飞行器的气动伺服弹性问题，提出一种结合线性自抗扰和自适应陷波器的综合控制方案。针对强耦合和强不确定性问题，采用线性自抗扰控制对总扰动进行快速估计和补偿。为了在最小化对刚体控制性能影响的同时实现对频率未知且时变的弹性模态的抑制，采用能够在线估计弹性频率的自适应陷波器。根据气动伺服弹性抑制问题的特点提炼出对自适应陷波器的性能需求。针对这些需求设计了两种基于递推最大似然法的多频率直接辨识方案——参数单独自适应和同时自适应方案。为了提高辨识算法在各种随机扰动下的鲁棒性，在此基础上提出一种在线有效性监督机制，并通过大量的数字仿真对两种方案进行了性能对比。最后，在弹性高超声速飞行器模型上进行仿真，仿真结果验证了所提控制方案的有效性。     

壁面辐射平衡DSMC方法及其在双锥构型中应用

临近空间高超声速飞行器在高空长时间飞行,受气流加热影响,飞行器表面温度显著升高,依赖地面试验和传统DSMC仿真预测的热流值明显高于飞行观测值,导致飞行器防热系统的保守设计.本文发展了一种基于壁面辐射平衡的DSMC边界模型,通过热流值反算辐射平衡壁面温度,并以此温度作为下一个时间步DSMC计算的边界条件,迭代更新至给出壁面温度的收敛值.基于该温度边界条件,开发了适用于轴对称构型的DSMC求解器,并以钝锥构型对计算模型和求解器进行了验证.重点针对激波风洞试验条件下的双锥构型,开展数值模拟研究,结果表明:该构型恒温冷壁条件得到的壁面压力分布和热流与风洞试验结果吻合,两种温度条件下的压力峰值差异约为15.4％,但是整体气动力特性差异仅约为0.33％;相对于冷壁,辐射平衡计算得到的前缘处热流峰值降低约50％,再附点处的热流峰值降低约三分之二;两种条件相结合,可以给出壁面热流的预测范围.     

超声速横向射流三维流场结构特征

为了研究超声速燃烧室内燃料与空气快速掺混过程的流场特性，基于可压缩Navier-Stokes方程，采用大涡模拟方法和高精度WENO-TCD混合格式对来流马赫数为2.68，喷压比为36的超声速横向射流流场结构进行数值研究。数值结果清晰描述了超声速主流与横向射流相互作用过程的流场结构特征，得到了三维激波形态的演变规律以及它们在强化混合过程中的作用。另外，因桶形激波背风面低压区处的斜压效应，射流气体在桶形激波背风面形成一对螺旋向上的反向涡对，反向涡对的卷吸作用诱导进入壁面边界层的主流向上运动，形成冲击射流。冲击射流以v=557m/s的法向速度向上冲击桶形激波背风面，因而在桶形激波背风面留下类三角锥面凹痕。     

适用于复杂流动的热气动弹性降阶建模方法

针对高超声速巡航类飞行器面临的复杂流动下的热气动弹性问题，发展了一种时变热模态适用的非定常气动力降阶模型，建立了基于流-热-固时空耦合分析策略的热气动弹性分析方法，对高超声速飞行器前体进气压缩面进行了实际飞行加热过程的时变颤振分析。结果表明，在局部高热流载荷下，压缩面固有模态频率和振型随时间发生了大幅度变化，而所建方法能够适应上述变化，在获得高置信度非定常广义气动力的同时，避免了重复性的非定常气动力数值计算，可将热颤振分析中的气动力计算效率提升若干量级；此外，在达到热平衡后，进气压缩面的颤振动压降为初始时刻的0.64%，使得飞行包线大幅收窄。相关方法有效缓和了热气动弹性分析效率与精度的矛盾，提升了热气动弹性问题的工程分析能力。     

非均匀来流下三维激波反问题的微元密切轴对称解法

传统的密切轴对称理论被广泛应用于均匀来流下的三维密切曲面激波反设计，为解决非均匀来流条件下的三维曲面激波反问题，提出了一种微元密切轴对称流场（MOA）求解方法。该方法沿激波面的周向和流向构建一系列微元密切面，在每个微元面内进行三维向二维流动的等效转换，从而突破了传统密切方法中不能有横向波后流动的限制。利用该方法编写设计程序，分别基于带攻角来流条件和外锥型流来流条件重构了标准内锥曲面激波，并与数值仿真结果进行了比较。结果表明，非均匀来流下激波曲面的三维形状均与预设形状完全一致，实现了非均匀来流下曲面激波形状可控。MOA方法在吸气式高超声速推进领域中前体/进气道一体化设计方面有重要应用前景。     

基于网格自适应的飞行器防热材料热传导系数辨识

为准确辨识高超声速飞行器防热材料热传导系数，构造了针对热传导偏微分方程系统的约束泛函极值问题，基于表征多项式逼近性能的Weierstrass定理，采用Lagrange多项式对热传导系数进行参数化建模，进而将无穷维的约束泛函极值问题转化为有限维的非线性规划问题，再利用基于动态优化方程的优化方法将此非线性规划问题转化为常微分方程初值问题进行求解。为从较有限的测量数据中准确地“学习”热传导系数，建立采用“过拟合”判别准则的网格自适应迭代算法改善辨识精度。研究表明本文采用的辨识策略与优化方法有效，辨识结果可以较准确地反映材料热传导系数的变化规律。