高速发动机油膜惯性对活塞裙润滑影响

利用一种新的迭代算法 ,分析了高速发动机油膜惯性对活塞裙润滑特性的影响 .该迭代方法是利用有限元法和差分法交替求解Navier Stocks方程和雷诺方程 ,据此导出了适用于高速发动机活塞裙的润滑计算的混和润滑雷诺模型 .新的模型借助惯性系数 ,引入了油膜惯性项 ;同时给出求解含有油膜惯性项的迭代步骤和有限元表达式。计算结果表明 :随着惯性系数和活塞裙的长径比的不同 ,油膜惯性会对油膜的摩擦力、压力和承载力产生不同的影响 ,这种影响对承载力尤为明显 .该润滑模型也可用于中、低速发动机的活塞裙润滑计算以及不计入惯性项 (惯性系数置为零 )的某些润滑问题求解 .     

挤压油膜阻尼器的热平衡分析

利用变分有限元法对挤压油膜阻尼器（简称ＳＦＤ）的热平衡问题进行了分析计算。文中根据雷诺边界条件，通过联立求解ＳＦＤ的雷诺方程、能量方程得到了ＳＦＤ的压力分布、滑油流量等参数。为ＳＦＤ的设计与计算提供了更为可靠的理论依据。     

高空台扩压器的数值模拟

从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发,利用时间相关法,采用ＭａｃＣｏｒｍａｃｋ两步显格式,结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,数值模拟了高空台扩压器的内部流场,并计算了扩压器的性能。结果表明,采用求解Ｎ－Ｓ方程的方法计算高空台扩压器性能是一种有效的方法。     

共轴对转螺旋桨的非定常气动干扰

为优化共轴对转螺旋桨的气动性能,研究前后桨之间复杂的非定常气动干扰,使用非定常雷诺平均Navier-Stokes （URANS）方程的数值模拟方法,以及滑移网格技术,对6×6构型的共轴对转螺旋桨进行了三维非定常气动干扰计算,获得了共轴对转螺旋桨的气动特性变化规律,并给出了优化气动外形及运行控制的建议。计算结果表明:来流马赫数为0.38的情况下,由于非定常气动干扰导致的整个共轴对转螺旋桨的拉力系数和功率系数的脉动幅值为1.5%左右,推进效率的脉动幅值在0.1%左右。一个旋转周期内,前后桨均出现了12次周期波动,后桨由于能够吸收一部分前桨产生的切向滑流能量,推进效率比前桨更高一些,而且能够提高整个桨的效率。另外通过压力系数以及轴向速度的径向分布情况看,70%左右的桨半径位置上非定常气动干扰强度相对较大。      

单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N-S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。      

金属橡胶挤压油膜阻尼器阻尼性能分析

为了降低发动机转子高速旋转产生的振动，提高发动机的使用寿命，设计了一种端部安装金属橡胶环的新型挤压油膜阻尼器，并在π油膜假设条件下，利用斯托克斯方程和雷诺方程，给出了新型金属橡胶挤压油膜阻尼器的流量表达式以及非线性油膜力的近似表达式。通过传统短挤压油膜阻尼器和新型金属橡胶挤压油膜阻尼器的流量特性和油膜力特性的对比分析，表明新型金属橡胶挤压油膜阻尼器具有更好的阻尼性能和良好的节流特性。并通过试验证明了新型金属橡胶挤压油膜阻尼器具有更好的阻尼减振效果。     

进口气流角对三维振荡涡轮叶片非定常流动影响的数值模拟研究

进一步发展了非定常雷诺平均N—S方程求解程序,并基于影响系数法,在三种不同振型下就进口气流角对三维涡轮振荡叶片绕流的影响进行了数值模拟研究。数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;进口气流角对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,且在不同模态下,进口气流角对振荡叶片非定常流动的影响规律表现不同。     

静偏心对挤压油膜阻尼器减振特性影响的 数值分析

由于装配误差或重力等原因,挤压油膜阻尼器在实际使用中难以避免产生静偏心。为了研究静偏心对挤压油膜阻尼器减振特性的影响,基于广义雷诺方程推导出静偏心条件下挤压油膜阻尼器雷诺方程并建立了静偏心挤压油膜阻尼器-转子系统的动力学模型,随后运用数值积分获取静偏心下转子系统的非线性响应。研究表明:转子进动半径随着静偏心增大而减小;静偏心显著影响转子轴心轨迹,静偏心时转子轴心轨迹在大多数转速下近似为椭圆而不再是圆,静偏心下突加不平衡响应瞬态振幅远高于无静偏心条件下瞬态振幅。     

基于滑动网格的螺旋桨复杂绕流数值模拟研究

利用SST湍流模型,求解雷诺平均Navier-Stokes方程数值模拟螺旋桨粘性绕流。为了实现桨叶的相对运动,采用了滑动网格技术,即将计算区域分为静止区域和随桨叶一起运动的旋转区域。以课题组设计的1.2 m螺旋桨为例,分别对不同风速下的桨叶粘性流动进行了非定常数值模拟,计算得到的螺旋桨气动特性与实验值吻合良好,证明了方法的正确性。     

湿式摩擦离合器挤压过程及接合特性分析

针对湿式摩擦离合器接合过程中的挤压机理及微凸体接触问题，综合考虑摩擦片表面油槽结构，材料渗透性等因素，根据修正雷诺方程和K-E（Kogut-Etsion）接触模型，对湿式摩擦离合器的挤压过程中的动压承载力、微凸体承载力以及转矩转速变化等建模分析，并采用有限差分法对修正雷诺方程求解，对挤压过程中的油膜压缩速度、油膜厚度变化率和片间承载力等挤压特性和转速、转矩等接合特性展开了仿真分析。并采用SAE#2试验机进行了相关试验获得了转速、转矩、压力、摩擦因数等数据，与仿真分析进行了对比。结果表明：湿式摩擦离合器接合过程可分为3个阶段，通常会在1 s内完成，在0.02 s左右开始从挤压阶段经压紧阶段在0.03 s左右过渡到全微凸体接触阶段，试验与理论分析结果一致。     

欧拉方程计算非定常流动的一种数值方法

提出了利用欧拉方程计算非定常流动的一种数值方法。在利用中心差分的有限体积法和双时间推进法求解非定常欧拉方程过程中，通过在固定物面边界上满足完全的非定常边界条件，计算出非定常流场。算例结果表明，本文方法具有高效、简便的特点，而且能描述激波的非定常运动和变化情况。     

带挤压油膜阻尼器的盘式拉杆转子双稳态振动特性

建立了带挤压油膜阻尼器的盘式拉杆转子动力学模型,利用非线性微分方程多周期解求解方法,对不同系统参数组合下的动力学方程进行求解,研究了带挤压油膜阻尼器的盘式拉杆转子双稳态振动特性的影响规律,比较了带与不带挤压油膜阻尼器系统稳定周期解的动力学特性.通过频谱分析得出:忽略挤压油膜阻尼器时双稳态区主要为基频振动,考虑挤压油膜阻尼器的盘式拉杆转子在双稳态区除了基频振动外,还有倍频的出现.      

变来流速度下旋翼翼型非定常气动特性分析

为分析变来流速度状态下的旋翼翼型气动特性,提出了利用翼型平移来模拟来流速度变化的数值方法.在此方法基础上,采用基于隐式LU-SGS(lower upper symmetric Gauss-Seidal)方法的非定常雷诺平均N-S(Navier-Stokes)(RANS)方程,模拟了SC1095旋翼翼型在定迎角 变来流速度及变迎角 变来流速度状态下的非定常气动特性.通过对比分析发现:翼型在变速度-定迎角状态下会表现出明显的非定常现象,产生了前缘分离涡,气动特性会出现明显的迟滞效应及波动现象,脉动速度越大,非定常效果越明显.并且基准速度越大,翼型气动特性的峰值越大；翼型迎角越大,非定常涡出现的也越早.考虑直升机旋翼翼型实际工作环境,在变速度-动态失速状态下,翼型最大迎角处的气动力会得到一定程度的削弱,在小迎角下的气动力得到一定程度的增强,且脉动速度越大,翼型的非定常特性也越强.      

基于动态边界控制的优化方法在高升阻比气动外形优化中的初步应用

本文是利用动网格技术和有限体积法来研究基于动态边界控制的气动优化的方法.该优化方法是通过在求解非定常流动的过程中计算优化的目标函数随设计变量变化的梯度.而利用参数化的结构动网格技术能有效、快速地生成与优化过程出现的外形对应的网格并进行计算,具有较高的效率.本文以典型的高升阻比外形优化为例,用非定常NS方程为流动控制方程进行优化,取得了初步的结果.     

对转螺旋桨流场气动干扰数值模拟

基于结构网格动态面搭接技术，通过求解三维非定常雷诺平均纳维-斯托克斯(RANS)方程，对单独螺旋桨(SRP)滑流流场进行了数值模拟，结果与风洞试验数据吻合良好，在此基础上深入研究了对转螺旋桨（CRP）的非定常滑流效应及前后桨之间的气动干扰现象。计算结果表明:所采用的非定常RANS方法能够良好刻画对转螺旋桨滑流流场的发展变化特征，适用于分析前后桨气动干扰特性。前后桨桨尖涡的复杂相互作用是导致对转螺旋桨气动干扰的直接原因，后桨桨尖涡对于前桨桨尖涡具有一定耗散作用。在气动干扰影响下，对转螺旋桨气动力会产生周期性波动，波动周期与螺旋桨桨叶数相关, 6×6对转螺旋桨在1个旋转周期内产生了12次波动。气流穿过对转螺旋桨会发生两次加速，因此在相同工况下，对转螺旋桨的拉力系数和功率系数相比于单独螺旋桨分别增大约1倍，效率提升约1.5%。      

长直机翼的颤振及混沌运动分析

采用非定常气动力并考虑几何非线性的影响,建立了长直机翼的气动弹性运动方程。运用伽辽金法对方程进行离散,通过数值模拟研究了机翼的颤振特性及混沌运动。结果表明:考虑几何非线性后,出现极限环振动的初始点与线性预测结果基本一致;不同机翼模型,机翼振动从收敛到混沌的过程不同,可由单个极限环振动经拟周期运动进入混沌,也可以由单个极限环到拟周期运动,再回到单环振动,然后经极限环的周期倍化进入混沌状态。     

振荡叶栅三维非定常流动数值分析

以全三维、非线性、非定常欧拉方程计算作用于振荡叶栅的非定常气动载荷。欧拉方程的数值离散采用了具有时、空二阶精度的MacCormack显式分裂的两步格式,用时间推进方法得到其数值解。应用动网格技术捕获叶片的运动,保证叶片振动过程中网格边界与叶片表面保持一致。以相移周期条件处理具有不同相位的叶片振动问题,同时发展了伪边界形状修正方法以简化周期边界条件处理。     

战斗机大迎角/过失速机动下的进气道气动特性

为掌握战斗机在大迎角和过失速机动飞行时进气道的稳、动态气动特性,采用基于动态嵌套网格的非定常雷诺平均Navier-Stokes （URANS）方程和大迎角风洞试验方法对某战斗机进行了研究,并通过大迎角和过失速机动飞行试验进行了验证。结果表明：大迎角稳态下进气道气动性能随迎角增大逐渐降低,天地相关性吻合良好,而计算仿真和飞行试验均捕捉了眼镜蛇机动下进气道的非定常迟滞效应。通过研究获得了战斗机在大迎角和过失速机动下的进气道气动特性,建立了过失速机动下进气道非定常非线性特性问题的研究方法。     

柔性扑翼气动结构耦合特性数值研究

微型扑翼体积小、重量轻,其柔性变形对气动特性有显著的影响。通过求解雷诺平均N-S方程(ReynoldsAveraged Navier-Stokes,RANS)和结构动力学方程,对微型柔性扑翼飞行器的气动结构耦合特性进行了数值模拟研究。针对微型扑翼的大幅运动,发展了适用于扑翼的气动结构耦合数值计算方法,研究了微型扑翼的气动结构耦合特性。通过求解雷诺平均Navier-Stokes(RANS)方程得到微型扑翼的非定常气动特性;利用哈密顿原理(Hamilton Principle)推导了扑翼的结构动力学方程,采用结构有限元方法对该动力学方程进行离散并求解,得到扑翼的动态结构特性;采用松耦合方法进行迭代。计算结果与风洞实验结果相比吻合良好,验证了所发展方法的有效性。在此基础上研究了惯性力和关键运动参数对柔性扑翼气动及结构特性的影响规律,有助于比较详细、全面地了解微型扑翼的气动机理,为柔性扑翼的设计提供了参考依据。     

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

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