在轨服务的相对运动耦合动力学建模与分析

研究在轨服务航天器逼近与捕获目标航天器的相对轨道姿态耦合动力学建模问题。考虑航天器姿态与对接位置的运动耦合,建立目标运行在任意轨道下的相对轨道姿态耦合动力学模型,并对模型中的运动耦合进行深入分析。设计一种非线性的输出反馈姿态控制律,将建立耦合动力学模型与CW方程进行仿真比较,验证轨道与姿态的运动耦合对两航天器对接点之间相对位置的运动影响。     

高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。     

常开洞工业厂房风致内压试验研究

通过一系列不同条件下的单个常开洞单跨厂房刚性模型风洞测压试验,分析了风向角、开洞位置、开洞形状以及开洞大小对内压平均值与脉动值的影响,并与基于稳态理论与非定常理论的内压预测值进行了比较。研究发现,单个常开洞厂房内压的空间分布是均匀的,几乎处处相等,仅在开洞附近测点的风压系数值稍有差别;内压的脉动特性主要来源于来流湍流作用、旋涡脱落与尾流干扰以及赫姆霍兹共振等三个方面的影响;内压的均值与均方根值受开洞条件以及风向角等因素影响较大;稳态理论计算得到的内压系数均值与试验所得非常接近,但是均方根值的预测偏保守。非定常理论推导得到的内压响应方程在风速方向正对开洞时可以给出良好的预测,但由于受旋涡脱落、尾流干扰等影响,斜风向时运用该方法预测内压是不合适的。     

应用自由变形法的翼型几何外形参数化研究

对一种新的基于NURBS方法的自由变形参数化方法进行了研究,该方法继承了NURBS方法的优点并避免了其缺点。利用该方法对两副翼型分别进行了参数化和优化设计、反设计,流动控制方程为N-S方程,目标函数梯度计算方法采用了离散共轭方法。算例表明,自由变形参数化方法能够满足设计需要。     

太阳同步轨道卫星初轨捕获策略研究

从分析轨道控制动力学方程出发,结合工程应用实际约束,针对具有回归、冻结等特性的太阳同步轨道,设计其初轨捕获策略。该策略将初轨捕获任务规划为4个阶段,前3个阶段进行发动机推力效率标校、大控制量的轨道面内修正和轨道面内的精调,完成轨道面内捕获;最后一个阶段完成轨道倾角捕获。最后对某太阳同步轨道卫星的初轨捕获策略进行了详细设计并通过仿真计算验证该策略的合理有效性。     

基于伴随方程和自动微分的雷达散射截面表面灵敏度计算

反隐身技术的发展对军用飞行器的隐身性能提出更高要求,针对经典电磁场伴随方法无法获得表面灵敏度,难以为优化设计提供直观指导,且梯度求解时需反复填充阻抗矩阵,当设计变量、入射角度较多时梯度求解效率下降的问题,提出了基于伴随方程和自动微分的雷达散射截面表面灵敏度计算方法,并结合基函数特点提出了表面灵敏度的稀疏存储和稀疏矩阵-矢量相乘方法,避免了直接计算表面灵敏度时计算量、存储量无法承受的问题。提出的表面灵敏度计算方法可以通过一次阻抗矩阵偏导数求解得到所有网格点的表面灵敏度,避免反复填充阻抗矩阵的问题。当入射角度改变时,求解任意数量设计变量梯度的计算量约为8次矩阵-向量乘,显著提升基于伴随方程的梯度计算效率,为作战飞机气动隐身一体化优化设计提供有力技术支撑。     

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

基于本征方程的柔性机翼非线性颤振分析

柔性机翼在气动载荷作用下常常会产生较大的变形,颤振特性会随之发生变化,针对此问题线性理论常常难以进行合理的预测。以几何精确本征梁模型建立了机翼的运动方程,耦合ONERA-EDlin非线性气动模型,建立了柔性机翼的非线性气动弹性分析模型。利用Newton-Raphson和Backward-Differentiation-Formula(BDF)分别求解机翼的静态变形和动态响应,基于机翼平衡位置附近的线性化方程来判断系统的稳定性,进而确定颤振临界速度。通过算例验证了模型的准确性,并分析了不同刚度、后掠角、机翼安装角等参数对颤振速度的影响。     

基于曲线坐标系浸没边界方法的折流燃烧室模拟

航空发动机折流燃烧室几何结构复杂,其高保真数值模拟需要高效的网格与边界条件处理方法。采用曲线坐标系隐式浸没边界方法结合大涡模拟-概率密度函数输运方程湍流燃烧模型开发自研软件,并实现WP11中折流燃烧室的高保真模拟。流动模拟中准确解析了该燃烧室中的三股主要气流,且三股气流分别约占进口流量的75%,12.5%和12.5%。两相燃烧模拟中针对拉格朗日框架下的液滴运动和欧拉框架下的湍流燃烧采用不同网格标记,模拟得到的出口径向温度分布规律与实验一致,平均相对误差为17.95%,表明基于本方法开发的自研软件能准确模拟折流燃烧室中的两相湍流燃烧现象。     

Computational aerodynamic analysis on perimeter reinforced （PR）-compliant wing

Abstract Implementing the morphing technique on a micro air vehicle （MAV） wing is a very chal- lenging task, due to the MAWs wing size limitation and the complex morphing mechanism. As a result, understanding aerodynamic characteristics and flow configurations, subject to wing structure deformation of a morphing wing MAV has remained obstructed. Thus, this paper presents the investigation of structural deformation, aerodynamics performance and flow formation on a pro- posed twist morphing MAV wing design named perimeter reinforced （PR）-compliant wing. The numerical simulation of two-way fluid structure interaction （FSI） investigation consist of a quasi- static aeroelastic structural analysis coupled with 3D incompressible Reynolds-averaged Navier- Stokes and shear-stress-transport （RANS-SST） solver utilized throughout this study. Verification of numerical method on a rigid rectangular wing achieves a good correlation with available exper- imental results. A comparative aeroelastic study between PR-compliant to PR and rigid wing per- formance is organized to elucidate the morphing wing performances. Structural deformation results show that PR-compliant wing is able to alter the wing＇s geometric twist characteristic, which has directly influenced both the overall aerodynamic performance and flow structure behavior. Despite the superior lift performance result, PR-compliant wing also suffers from massive drag penalty, which has consequently affected the wing efficiency in general. Based on vortices investigation, the results reveal the connection between these aerodynamic performances with vortices formation on PR-comoliant wing.