非线性数值模拟柱面气膜密封动态特性

针对气膜密封研究中密封气膜动态特性分析难点问题,采用非线性数值模拟法对柱面气膜密封系统的动态特性进行研究[D].建立了柱面气膜密封系统的动力学模型,给出了基于有限元数值方法的瞬态气体雷诺方程（含时间项）与系统动力学方程的耦合求解流程和方法,数值仿真给出了柱面密封环的中心轨迹图、确定了一定转速下的临界无量纲质量、不同密封环无量纲质量时等量位移扰动下系统的响应图,分析了密封环无量纲质量对系统稳定性的影响.结果表明:光滑柱面气膜密封系统存在临界无量纲质量,临界时系统受扰后处于半频涡动状态,随着密封环无量纲质量的增加,受扰后系统收敛时间变长,稳定性变差.      

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

超燃冲压发动机内外流场三维燃烧数值模拟

以某型煤油燃料超燃冲压发动机整机模型为研究对象,对发动机外流场的流动与发动机内的流动燃烧过程开展了耦合三维数值模拟研究.燃烧室壁面静压计算结果与实验数据吻合良好,平均相对误差为4.9%,验证了数值方法的有效性.结果表明:内流场的流动燃烧状况既受到外流场捕获空气流量的限制,又会反过来影响着外流场的激波结构;泄压孔对保持发动机中各部件流量匹配起到关键作用;由于进气道内激波与边界层的相互作用以及泄压孔诱发的斜激波的影响,流场在第1处燃油喷射孔附近呈现出燃烧流场不对称特征.      

固体火箭发动机转角涡脱落分类

应用三维大涡模拟(LES)方法对C1x模拟试验器进行了数值模拟,以研究转角涡脱落引起的声涡耦合现象.结果表明:与二维和轴对称LES相比,三维LES方法可以准确捕获转角涡脱落运动的三维特性及锁频现象.根据转角涡脱落是否在腔内形成声反馈循环过程可将其分为两类问题,第1类形成声反馈循环涡;第2类形成自耗散涡.燃烧室的台阶扩张比、台阶前后长径比是影响两类转角涡脱落现象的重要结构参数.     

某型直升机旋翼载荷模拟系统力耦合性分析

直升机主桨载荷模拟系统是一个多输入多输出液压伺服加载系统,该系统助力器和加载缸之间存在复杂的非线性耦合关系。通过建立系统的数学模型分析了该系统存在的力耦合关系,发现其造成的影响在系统允许的误差范围之内,因此可忽略这部分耦合因素。从而为简化系统控制策略提供依据。     

高机动导弹气动／运动／控制耦合的风洞虚拟飞行试验技术

解决先进飞行器大迎角高机动飞行时的气动／运动非线性耦合问题,需要发展基于非线性理论的风洞试验技术,即风洞虚拟飞行试验技术。该试验能够实现较为逼真的模拟飞行器机动运动过程,气动和运动参数的实时同步测量,以及飞行控制律的集成验证与优化,从而达到探索气动／运动耦合特性和机理的目的。本文介绍了风洞虚拟飞行试验的模拟方法、关键技术及其解决措施,并针对典型导弹模型开展了虚拟飞行验证试验。试验结果表明：目前已经初步具备适用于导弹模型跨声速气动／运动／飞行控制一体化研究的风洞虚拟飞行试验能力。     

轴向柱塞泵滑靴副的热结构耦合特性

为了改善轴向柱塞泵滑靴副的耐磨损性能,建立了滑靴与斜盘摩擦副的瞬态热结构耦合模型,分析压力冲击条件下滑靴的表面温度、应力以及变形的变化规律.研究结果表明:某型轴向柱塞泵中滑靴温度随柱塞腔压力呈周期性变化,滑靴温度范围为45.5~49.8℃,且滑靴的最高温度出现在泵的吸排油过渡区.当滑靴处于泵的排油区时,滑靴的最大轴向应力为250MPa,集中在滑靴油腔与密封带之间的边缘区域.滑靴的轴向应力分层显著,引起滑靴的变形分化,其变形量为12.5~15μm,出现在滑靴的边缘.由于滑靴的输入热流密度增强磨粒的剪切力,加剧滑靴表面的微切削和挤压变形,导致滑靴表面出现条状剥落和凹坑磨损,呈现出黏着和磨粒磨损特征.      

非结构变形网格和离散几何守恒律

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉（ALE）方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格（MSA）及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数（RBFs）和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律（DGCL）概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。     

PID控制器与CFD的耦合模拟技术研究及应用

飞行控制系统（FCS）与计算流体力学（CFD）的耦合求解是一个崭新的研究领域。传统的飞行控制系统的工程仿真方法依靠气动力模型或气动力数据库得到不同飞行姿态的气动力;而当前方法通过耦合求解Navier-Stokes方程和刚体动力学方程（RBD）以获取飞行器运动过程实时流场和非定常气动力。由于充分反映了气动力的非定常、非线性效应,因而从根本上保证了飞行控制系统仿真的精度。以方形截面导弹俯仰姿态控制为例,首先给出了系统的传递函数,并基于系统在单位阶跃舵偏操纵下的开环响应特性,提出了传递函数的修正方法,进而设计了该外形俯仰姿态控制的PID控制器。数值模拟了不同控制参数时,P控制器、PD控制器和PID控制器的控制效果。针对不同的控制指令,根据建立的控制律,数值模拟了飞行器在PID控制器作用下的实时响应过程,最终成功实现了对飞行器的俯仰姿态控制。研究发现,当飞行器作慢速机动时,工程仿真与CFD数值计算的结果吻合很好,两种方法可以互相验证;但快速机动时,两种方法给出的结果差异明显,基于CFD的耦合模拟方法由于模拟了飞行器运动和舵面偏转导致的非定常流动过程,其结果比基于静态气动力的工程方法的可靠性更高。在大攻角和快速机动等非定常效应较强时,采用CFD方法评估和验证飞行控制系统是很有必要的。     

多级叶盘结构动力特性分析的一种方法

因各级转子叶片数的不同,基于传统循环对称法的通用有限元程序无法直接用于多级叶盘结构动力特性分析.发展了一种方法,分别对每一级扇区模型施加相应的循环对称边界条件,级间则通过施加自由度耦合来保证位移场的连续.若对每一级叶盘结构选取合理的扇区数目进行建模,最大程度地考虑级间耦合效应,则精度可进一步提高.由于循环对称边界条件和级间耦合通过对位移自由度施加约束方程来实现,因此该方法易于在通用有限元程序系统中实现.对两级叶盘结构分别采用该方法和全环模型进行动力特性分析,两种方法的频率差控制在3%内.探讨了多级叶盘结构中级间耦合对模态的影响,在低节径模态两级叶盘结构耦合效应较为明显,随着节径数的增加,耦合效应降低.