基于matlab的飞行器并行MDO方法的研究与实现

在Parmatlab的基础上开发出适合网络通讯的maflab工具箱MD0504 toolbox，此工具箱提供了在matlab环境下进行并行计算的函数。在此工具箱的基础上开发出进行多学科设计优化的软件GSEOPT，GSEOPT在并行环境下实现了基于全局灵敏度方程（GSE）的单级优化算法，用并行的方法求得了全局灵敏度，在系统层进行优化。最后用此软件对某飞行器进行了并行优化设计，给出了计算结果，并与参考文献结果进行了比较，得到了更好的结果。     

粒子群算法在APU转子平衡工艺中的应用

在辅助动力装置（APU）装配工作中,转子组件平衡是极其重要的步骤。由于转子组件零部件的动力学参数未知,为达到维修工艺所需的剩余动态不平衡量,需要对各零部件的装配角度进行多次匹配尝试,但现装配工作依赖工作者经验而缺乏理论性指导,所测数据无法得到有效利用,造成极大的工时浪费。本文运用欧拉公式对不平衡量的表达进行简化,并使用粒子群算法求解,编写计算机应用程序应用于APU转子平衡工艺,可提高装配质量与效率。     

随机场TPDF在非结构网格平台的开发及检验

为实现对燃烧室内部燃烧过程的高精度数值模拟,增强航空发动机燃烧室研究手段,采用开源的非结构网格平台Saturne,开发可以耦合详细化学反应机理,高精度解析湍流燃烧过程的概率密度函数输运方程（TPDF）湍流燃烧模型。采用随机场方法求解TPDF方程,在原有程序基础上研发了TPDF程序模块及配合的加速算法等模拟单元,发展了针对燃烧室燃烧过程模拟的软件功能。采用射流火焰Flame D和旋流火焰TECFLAM对新软件进行了测试,结果表明：模拟结果与试验数据趋势一致,精度尚可。对某型燃烧室性能模拟验证中,计算的出口温度分布与试验结果一致,平均温度误差小于3.8%,新的软件可进一步应用于燃烧室性能的研究。     

柔性空间机械臂在轨操作仿真与分析

针对挠性航天器利用柔性空间机械臂在轨操作目标进行分析.首先利用Kane方程和假设模态法对挠性航天器上安装有柔性空间机械臂的系统进行动力学建模.其次,采用修正的罗德里格斯参数描述机械臂末端相对服务航天器的姿态,利用五次多项式对机械臂末端的相对位置与姿态进行规划,并将目标航天器的相对运动进行补偿,基于雅克比矩阵的广义逆求解机械臂关节运动规律.然后,将反馈控制与扩张状态观测器结合,分别设计了航天器姿态稳定控制器和机械臂轨迹跟踪控制器.最后,对柔性空间机械臂捕获目标航天器以及安装模块的过程进行闭环数值仿真,结果表明,所设计的控制器能够使机械臂跟踪期望轨迹,同时使得航天器姿态趋于稳定,机械臂可以较高精度完成在轨操作.     

结冰飞机着陆阶段飞行安全包线确定及操纵应对策略

结冰会导致飞行安全包线收缩、严重威胁飞行安全,研究结冰后安全包线变化对于操纵应对策略设计及提高飞行安全具有重要意义。以美国国家航空航天局（NASA）的项目飞机GTM（Generic Transport Model）为研究对象,通过对飞机的气动参数进行多项式拟合,建立了结冰飞机纵向通道的动力学模型。为了得到能随结冰程度变化的安全包线,将可达性分析理论引入到对结冰飞机着陆过程的安全性分析。提出将正向可达集与反向可达集的交集作为飞行安全包线,其中可达集的确定是基于水平集方法求解哈密尔顿-雅克比方程的最优解。最后针对不同程度的结冰条件进行了操纵时域验证,并提出了相应的操纵控制策略。研究结果表明,轻度结冰对安全包线影响较小,整个着陆过程的飞行状态始终能在最优控制指导下保持在安全包线以内;但对于重度结冰,飞行安全包线收缩严重,常规操纵已经很难使飞机达到着陆要求,需要进行飞行状态改出处理。研究结果为指导飞行操纵及实时包线保护打下基础。     

径向磁气组合轴承系统的动态性能

为探究磁轴承失效时能否减轻转子跌落带来的损害,将人字槽径向动压气体轴承引入磁轴承转子系统,研究磁气组合轴承对系统动态性能的影响和动压气体轴承的支承特性。采用有限差分法和小扰动法求解气膜厚度方程和雷诺方程,研究动压气体轴承的静动态特性。对磁轴承电磁力进行了分析,采用磁轴承不完全微分PID（比例-积分-微分）控制策略,对系统进行了理论模态分析和试验模态分析,完成了系统高速旋转试验,测试了动压气体轴承在不同偏心率和转速下的承载力。结果表明,引入动压气体轴承可提高系统的动态性能,在磁轴承失效造成高速转子跌落瞬间,偏心率趋近于1,两个径向动压气体轴承能够产生较大承载力,减轻转子跌落造成的损害。      

基于超椭圆方程和序列响应面法的回转壳开孔形状优化

含有多种功能开孔的轻质回转壳结构是航空航天结构系统内的常用支撑构件,其开孔形状直接影响结构的静动态性能。以回转壳结构为对象,基于超椭圆方程和坐标映射变换推导了回转壳开孔边界的参数化表达,开展了开孔形状动力学优化研究。为提高结构优化计算的精度、效率和收敛性,提出了准等弧长方法和基于均匀设计的序列响应面近似建模方法(SRSM),以分别实现空间超椭圆曲线的精确逼近、减少结构有限元分析成本和加快迭代收敛。以非支配排序遗传算法Ⅱ(NSGA-Ⅱ)作为响应面模型求解算法,结合有限元分析构建了回转壳开孔形状优化设计流程,开展了最大化结构一二阶频率带隙的典型回转壳结构开孔形状优化设计。结果表明,基于超椭圆方程和序列响应面法的开孔优化方法获得了有效改进结构动态特性的回转壳开孔形状,对开展计算耗时工程结构形状优化设计具有一定应用价值。     

基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.      

基于Deform软件二次开发和BP神经网络的TA15多向锻造微观组织预报

通过热压缩实验构建了TA15合金的热变形真应力-应变曲线,以此为基础分别建立了合金双相区及单相区温度区间的热变形本构方程;基于热压缩试样动态再结晶的统计数据建立了TA15合金的动态再结晶模型。借助Deform提供的二次开发功能实现相关数学模型的程序化,制定正交实验方案,实现了TA15合金多向锻造变形的微观组织仿真。通过正交实验分析得出各项因子的影响对象及强弱差异,提出了双相区及单相区温度区间内的多向锻造最佳因子组合。建立了TA15合金多向锻造变形微观组织的BP （Back Propagation）神经网络预测模型,将预测结果与有限元仿真结果进行比较,结果表明两种方法的预测结果基本一致,但神经网络具备有限元仿真难以实现的良好细节预测能力,能更为细致地实现对微观组织分布状态的划分。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.