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

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

一种基于非Newtonian流变学模型的牵引力计算公式

采用部分弹性流体动力润滑热分析理论及Ree-Eyring流变学模型对滚滑接触状态下航空合成润滑油的牵引特性进行了较为详细的研究,获得了航空合成润滑油牵引力系数随载荷、滚动速度、润滑油粘度及滑滚比变化的规律,并将计算结果与国外相近的实验结果进行了比较。此外,本文还根据研究结果提出了一种适用于航空高速滚动轴承动力学分析的牵引力系数计算公式。     

曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。     

长拉杆-止口连接弯曲刚度损失及对转子系统振动响应影响

航空燃气轮机为了实现高负荷、轻质化的追求,在转子结构设计中,趋向于提高转速和加大长径比。这使得转子系统弯曲模态临界转速降低,转子在工作转速范围内不可避免会产生一定的弯曲变形。转子弯曲变形会影响连接结构界面接触特性的变化,使其连接结构局部弯曲刚度产生损失。因此,对于工作转速靠近弯曲临界的高速转子系统,需要考虑连接结构界面接触状态变化对转子系统振动特性的影响。以高负荷的长拉杆-止口连接转子系统为对象,分析连接界面接触应力分布特性,提出连接结构弯曲刚度损失修正方法,以此为基础建立界面连接转子动力学模型。通过对止口连接三级轴流压气机转子结构动力学特性的仿真和试验研究表明,在靠近弯曲振型临界转速下,转子连接界面接触状态的变化会产生弯曲刚度损失,对转子动力学特性具有显著影响。     

梯形齿车轮月面牵引性能的离散元分析

采用离散元法分析梯形齿车轮在月面软土地形的通过性,提出不规则形状颗粒群系统的离散元建模技术,建立车轮牵引性能参数的细观表达式.基于模拟月壤的直剪实验数据,通过反复调试获得最佳土单元模型参数.轮-土作用的离散元分析结果表明:车轮在月壤地形的牵引性能与滑转率和轮齿结构有关,增加滑转率可获得更高的挂钩牵引力,但同时车轮所需的驱动转矩提高;增加轮齿高度可获得更大的土壤推力,但行驶阻力也相应增大;轮齿数量对车轮牵引性能影响较小,但影响车轮在行驶中的平顺性.     

双流制窄轨电力机车用异步牵引电机设计

将牵引电机应用于双流制窄轨电力机车时,由于窄轨机车轨距小,所以牵引电机安装空间异常狭小,比一般牵引电机在安装空间、电机结构、电磁负荷等方面的要求更加苛刻。针对窄轨双流制电力机车的运用特点和特殊要求进行了分析,提出了牵引电机的设计要求和关键技术难点,并对关键技术难点给出了相应的分析和解决措施。对该牵引电机的制造过程及电机试验测试情况进行了论述。试验结果表明,该牵引电机达到了设计预期,完全能满足设计要求。     

新型航空润滑油油膜拖动力计算研究

在自制的试验装置上测量了新型航空润滑油油膜拖动力;给出了高剪切率弹流润滑下的Ery-ing特征应力、极限剪切应力、剪切弹性模量的新型航空油流变参数拟合计算式;确定了该油本构方程;建立了该油油膜拖动力计算式, 并进行了试验验证.试验结果表明:该油在试验条件下表现为粘弹性;拖动力的理论计算值与试验测量值相当吻合.给出的拖动力计算方法可用于指导润滑分析领域的工程计算.      

Maxwell方程的高阶间断有限元数值解法

采用高精度方法求解时域Maxwell方程,方程的空间离散采用基于计算流体力学(Computational fluid dynamics,CFD)领域的高阶间断有限元格式,非定常时间迭代采用四步龙格-库塔格式。为了提高计算效率,本文采用了Quadrature-free implementation和网格分区并行技术。数值结果表明,采用高阶格式的情况下,采用稀疏网格便可以得到高精度数值解。另外由于本文的方法基于非结构网格,因此非常适合计算复杂外形的情况。     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。