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

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

各向异性单晶叶片强度寿命研究—第I部分:晶体滑移本构模型及应用

单晶叶片技术是提高航空发动机及地面燃气轮机性能、寿命及可靠性的关键技术之一，但单晶材料机械、力学性能的各向异性特性制约了其发展和应用，对其工程应用及应用的理论基础提出了挑战。课题组开展了各向异性单晶叶片强度分析和寿命预测方面的一些研究工作。包括：建立并验证了弹塑性、蠕变滑移本构模型及蠕变持久寿命预测方法；进行了不同晶体取向DD3单晶在不同温度、不同速率或不同温度、不同应力水平下的拉伸试验及蠕变试验，这些实验数据及由其反映的单晶中、高温各向异性特性对单晶材料的应用具有重要意义。此外还进行了某种单晶叶片的实验研究。作为上述研究的应用，对某发动机单晶涡轮叶片进行了强度分析和寿命预测。本文这一部分介绍本构模型及应用，实验研究将在第II部分介绍。     

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

采用线性稳定性分析的方法对粘性气体中的粘性液体射流的分裂与雾化机理进行了分析,数值计算表明:液体射流分裂与雾化过程中存在一临界气体韦伯数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复合固体推进剂的力学特性。     

聚氨酯泡沫塑料拉伸本构关系及其失效机理的研究

 通过 3种密度硬质聚氨酯泡沫塑料的拉伸实验,研究了它们的应力 -应变特性及其应变率效应。同时,为确定泡沫塑料拉伸失效机理,对材料进行了扫描电镜下的细观拉伸实验,观察了试件表面胞体的变形与失效过程,再结合宏观拉伸试件断口的扫描电镜分析,进一步讨论了泡沫塑料在拉伸加载下的失效机理。基于泡沫塑料的拉伸应力-应变曲线,用数值方法拟合了泡沫塑料的拉伸本构关系。     

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

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

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

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

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

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

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

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