跟踪问题最优控制律精细积分

构成有限时间最优跟踪系统的控制律需要求解 Riccati微分方程及外部控制输入向量满足的微分方程 ,前者是非线性矩阵微分方程 ,后者是变系数线性微分方程。在结构力学与最优控制的模拟理论基础上所发展的精细积分方法借鉴了计算结构力学中的算法 ,可以精确有效地求解这些微分方程。这种方法的特点之一在于步长幅度变化较大时 ,Riccati微分方程的数值解仍可以保持很高的精度 ,并且变系数线性微分方程的求解亦可纳入其体系而不必用通常的差分方法。本文介绍了用精细积分方法求解这些方程的过程 ,并给出了数值算例。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

粘性气体中粘性液体射流分裂与雾化机理研究

采用线性稳定性分析的方法对粘性气体中的粘性液体射流的分裂与雾化机理进行了分析,数值计算表明:液体射流分裂与雾化过程中存在一临界气体韦伯数We2c=1,We2＜We2c对应的是射流分裂过程,We2＞We2c对应的是射流雾化过程,射流分裂过程和雾化过程的机理有所不同.当We2＜1时,We2对射流分裂过程具有稳定性的作用;当We2＞1时,We2对射流雾化过程起着不稳定性的作用.液体Reynolds数Re1在整个射流过程中始终起着不稳定性的作用,气体Reynolds数Re2的作用却相反.气液密度比Q,即气动力在整个射流分裂与雾化过程中始终起着不稳定性的作用.     

非结构网格中LU-SGS隐式算法的非平衡性影响

在非结构网格中用LU-SGS隐式算法求解欧拉方程时,两个近似分解因子的项数可能相等,为平衡;也可能不相等,为非平衡.采用对网格重新编号的方法很难达到我们希望的平衡.本文对非平衡性的影响进行了探讨.对二维问题,四边形单元的非结构网格,设计出平衡与非平衡的编号方式.先对标量模型方程分析LU-SGS隐式算法的增长因子,然后通过数值试验来验证这种非平衡性的影响.结果表明,尽管非平衡时也能达到收敛,但平衡却远优于非平衡的情况.     

HTPB复合固体推进剂含损伤和老化本构研究

通过热氧老化实验、声发射损伤检测实验和蠕变柔量测试实验确定HTPB复合固体推进剂非线性本构关系，并与Schapery非线性本构方程进行比较，结果表明，含损伤老化模型较好地解释了HTPB复合固体推进剂的非线性特征，可以应用分析HTPB复合固体推进剂的力学特性。     

坐标系转换矩阵及其微分方程

在坐标单位矢量导数为非零条件下，通过建立飞行器的转换矩阵及其微分方程，将飞行器的动力学、运动学方程表达得更清晰、简洁。所得结论易于理解动静两种坐标系在时域内的联系，其数学表达式方便用MAT—LAB语言实现，有利于飞行器的建模与仿真。     

基于B-P本构模型的涡轮盘应力分析

针对Bodner-Partom统一粘塑性本构方程,采用Euler前差显式积分方法,编制了用户子程序,把该本构方程结合进了大型通用有限元软件ANSYS,进行了一维和二维的算例验证,并针对某涡轮盘进行了三维应力分析。该程序可考虑材料变形的率相关、循环加载等情况,适用于航空发动机高温部件非线性结构分析。      

非壅塞固体火箭冲压发动机二次燃烧试验研究

通过直连式试验研究非壅塞固体火箭冲压发动机二次燃烧影响因素。试验结果表明,当空燃比在一定范围内变化时,若空燃比变大,则燃烧效率升高,当空燃比达到一定程度后再增加,则燃烧效率降低;对于铝镁贫氧推进剂取较小的后置长度时燃烧效率较高;与两股燃气射流向外喷射相比,两股燃气射流向内喷射的燃烧效率明显高;燃气射流与空气流在进气道出口直接撞击不利于燃烧效率的提高。      

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

边条翼布局双垂尾抖振的数值模拟

对边条翼布局的双垂尾抖振进行了较为深入的数值模拟研究。模拟来流马赫数为0.2,迎角为10°~40°。通过非定常Euler方程计算各迎角下的非定常流场及垂尾根部弯矩系数。并将随时间脉动的根部弯矩系数进行计算得到根部弯矩系数均方根值,从而得出根部弯矩响应大小随迎角的变化曲线。结合流场特性对该布局双垂尾抖振的发生机理及抖振响应随迎角的变化规律作了深入分析。结果表明:该边条翼布局双垂尾抖振主要是由于边条破裂涡作用在垂尾上的脉动载荷引起的。最后,将垂尾根部弯矩响应的计算结果与该模型双垂尾抖振的风洞实验结果作了比较,结果符合得较好。