基于Newton-Raphson算法的径向轴承瞬态热流体动力润滑的研究

建立了径向轴承在载荷和速度突然变化时的三维数学模型，模型中考虑了轴瓦的热变形，在油膜和轴瓦交界面采用热流连续的理想边界条件，数值模拟轴承的瞬态温度场，并对轴承的瞬态性能进行分析。在每一瞬时，用Newton—Raphson算法同时求解Reynolds方程、膜厚方程和轴颈运动方程获得轴承油膜的压力分布和轴颈中心的运动速度，然后数值积分压力分布得到轴承的油膜力，差分运动速度得到轴颈中心位置和运动加速度。用一有效的有限差分法同时求解油膜和轴瓦的温度控制方程。最后将Reynolds方程和能量方程通过节点压力和温度相耦合获得轴承的瞬态三维温度场。结果表明本所介绍的方法收敛快，大大节约计算时间。     

可倾瓦径向滑动轴承绝热瞬态过程分析

为了研究可倾瓦径向轴承在载荷扰动下的绝热瞬态行为.对瞬态油膜压力和温度分布采用了时变的Reynolds方程和能量方程, 建立了轴颈及瓦块动力学方程.应用有限差分及Newton-Raphson法对方程联立求解, 模拟了阶跃载荷冲击工况下可倾瓦径向轴承热瞬态非线性响应过程.给出瞬态过程中油膜最高温度、最小膜厚等参数的变化规律.结论认为瞬态过程中, 油膜温度和膜厚变化量较大, 油膜温度和膜厚响应有一定的超调量, 瞬态过程中有可能因油膜温升过高或膜厚太小而导致失效.      

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

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

一种三维弹性机翼气动优化新方法

基于控制理论方法和有限元分析软件,通过将网格生成、流场计算、有限元分析、粘性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,研究与发展了计及静气动弹性影响的三维机翼气动优化设计方法,其中,雷诺平均Navier-Stokes方程为主控方程,气动载荷和结构静弹性变形量由气动/结构方程的耦合迭代求解得到,目标函数梯度信息由共轭方程数值求解得到。典型大展弦比机翼气动减阻设计结果表明：研究及发展的计及静气动弹性影响的三维机翼气动优化设计方法是有效的,能够有效考虑静气动弹性的影响。     

航空燃气轮机挤压油膜阻尼器油膜温度场的理论分析

油膜温度是挤压油膜阻尼器的一项重要工作参数。本文运用有限元法对挤压油膜阻尼器油膜温度场进行分析,即根据油膜边界条件,联立求解适用于阻尼器的雷诺方程、能量方程、滑油粘温特性方程得到油膜温度分布。文中进一步分析航空发动机阻尼器的工作环境,提出启动阶段,滑油温升由于阻尼减振发热所致;稳态阶段主要由于油膜环向滑油的强迫对流换热,并提出用试凑法来求解稳态滑油工作温度。     

稀薄效应对可倾瓦动压气体轴承性能的影响

以微型三可倾瓦动压气体轴承为研究对象,采用连续模型、一阶滑移模型和WU新滑移模型速度边界条件,建立考虑稀薄效应的滑移修正雷诺方程。结合牛顿迭代法和有限差分法求解修正雷诺方程,计算考虑稀薄效应的三可倾瓦动压气体轴承各瓦块压力分布和承载力。计算结果表明:三可倾瓦动压气体轴承偏心率越大,转速越大,压缩系数越大,轴承承载力越大;考虑稀薄效应后,计算出的可倾瓦动压气体轴承承载力明显下降,且随着克努森数的增大,WU新滑移模型计算得到的承载力明显低于一阶滑移模型。     

行星减速器油膜均载的分析计算

运用流体动力润滑理论分析了油膜均衡行垦齿轮间载荷的机理,解释了共转滑动轴承的动态特性。利用数值方法求解了Ｒｅｖｎｏｌｄｓ方程。计算结果和试验结果有良好的一致性,并表明：共转滑动轴承的油膜柔度随载荷的增大而减小;随轴承间隙的增大而增大;随转速的增大而减小。     

载荷扰动下可倾瓦推力轴承绝热瞬态过程分析

研究了可倾瓦滑动推力轴承在阶跃载荷扰动下的绝热瞬态行为.通过对轴承油膜参数及轴承元件的动力学建模,获得了大载荷扰动工况下可倾瓦推力轴承热瞬态运动参数的非线性响应曲线.由于瞬态分析中大的计算量及油膜力的非线性,计算中采用了Newton-Raphson法对模型进行求解.最终给出了瞬态过程中轴承最高温度、最小膜厚等参数的变化规律.结果表明:如果设计和运行不当,轴承有可能因为瞬态过程中温升过高等原因失效.      

一种关于静压气体轴承节流孔系数的计算方法

基于层流边界层方程的分离变量算法和雷诺方程的解析算法，提出了一种关于单节流孔静压气体止推轴承的节流孔系数的计算方法。该方法通过比较层流边界层方程计算获得的气体轴承的质量流量和雷诺方程计算获得的质量流量计算获得了节流孔系数。将计算获得的节流孔系数和节流孔系数为常数0.8代入单节流孔气体止推轴承的雷诺方程中，计算获得的承载力与分离变量算法求解层流边界层方程获得的承载力进行对比，可以发现，相对于采用节流孔系数为0.8来说， 采用该计算的节流孔系数求解雷诺方程的承载力与分离变量算法求解获得的承载力结果精度最大提高了8%。从而验证了该计算节流孔系数方法的正确性。      

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

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

航空发动机帽罩热气膜防冰的加热特性

针对复合材料帽罩采用的一种热气膜防冰系统,采用计算流体力学方法研究了其内部冲击换热和外部气膜加热效果。发展了二维轴对称计算方法并与全三维计算结果进行了对比,搭建了热气吹风实验台,利用红外测温仪和热电偶测量帽罩温度,验证了数值计算方法。全面研究了不同射流雷诺数、相对冲击距和气膜缝结构参数对热气加热特性的影响,结果表明:增大射流雷诺数有利于提高内部冲击换热效果和外部气膜加热效果;最佳相对冲击距随射流雷诺数的增大而增大,射流雷诺数为10000~40000范围内最佳相对冲击距在5~8内取得;气膜射流能大幅提高缝后外壁面温度,加热效率与气膜缝宽和位置有关,前缘开设的气膜缝还会提升前缘冲击换热效果。      

带气膜出流的旋转叶片冲击冷却的实验研究

用实验方法研究了放大 4倍的涡轮叶片前缘模型在旋转状态下,带气膜出流时叶片内冲击面的换热特性。实验结果表明:雷诺数的增加显著地增强了换热,但叶片前缘和尾缘的增加是不同的,而浮升力对换热的影响是复杂的;在低雷诺数的情况下,旋转对换热影响不明显,而在高雷诺数时,会使换热降低。      

预旋进气位置对转静盘腔换热影响的数值研究

对具有不同预旋进气位置的转静盘腔内的换热进行了数值模拟.在进气流量恒定的条件下, 改变预旋孔在静盘上的径向位置, 研究了旋转雷诺数、预旋角、转静盘间距等参数对进出口气流的无量纲总温降和转盘换热效果的影响.研究表明:旋转雷诺数的增大和转静盘间距的减小都能够增强气流对转盘的换热效果;当预旋孔的径向位置增大时, 无量纲总温降和转盘表面的平均努塞尔数都是先增大后减小, 预旋孔存在一个最优的径向位置.      

液体火箭发动机燃烧室壁液膜冷却的数值模拟

为研究液体火箭发动机的液膜冷却,建立了液膜模型。考虑核心气流与液膜间的对流传热,辐射传热以及壁面与液膜的对流传热分析传热量,由液膜的卷吸和液膜的蒸发计算传质,并由气液界面和液固界面的摩擦力分析流动情况。在400N小发动机内流场数值模拟中采用了该液膜模型,计算得到的壁面温度分布与试验结果符合较好,表明该模型是合理可行的。改变发动机燃烧室半径和圆筒段长度,将数值模拟结果对比分析发现:在一定范围内随着半径和圆筒段长度的增加,液膜长度减小,室壁温度升高,冷却效果变差。研究结果可为发动机的设计提供参考。      

旋转状态气膜换热机理的数值研究

大涡模拟研究了旋转对单孔平板气膜冷却的影响,气膜出流的雷诺数为2600,吹风比为0.5,计算了静止和旋转数为0.02两种状态的流动和换热,从湍流结构演化的角度讨论了旋转状态气膜换热的机理.结果表明:①发卡涡是决定气膜换热的主要湍流结构,旋转状态发卡涡的形状和运动规律主要受哥氏力的影响;②气膜在哥氏力作用下向高半径方向偏...     

叶片弦中区内部气膜孔局部换热特性实验

对叶片弦中区内部有、无冲击射流的气膜出流冷却方式中,冷气侧气膜孔局部换热特性进行了实验研究。无冲击射流时,主要研究了来流雷诺数、气膜出流与横流密流比的变化对气膜孔局部换热特性的影响;研究发现,气膜孔局部的换热均随两者的增加而强化,且孔后的换热要好于孔前的换热,越靠近孔的地方换热越强,并对倾角为30°和90°时气膜孔的“溢流效应”进行了比较。有冲击射流时,改变冲击雷诺数、横流射流密流比等流动参数,通过对数据整理得出了这些参数对其局部换热特性的影响规律。      

前台阶流场的传热实验研究

本文尝试用瞬态实验方法, 在卢德威格风洞实验台上, 对小台阶高度 (H=4.0mm)的前台阶二维流场进行传热实验研究, 测定了当雷诺数为 11800至32200时的流场内热流通量和努谢尔特数的变化规律以及努谢尔特数随雷诺数变化的准则关系式。实验结果是令人满意的。      

基于三维实体建模的刷式密封传热机理数值研究

建立了基于三维实体建模的刷式密封传热特性求解模型,在验证数值模型准确性的基础上,分析了刷式密封流场与温度场分布特性,研究了压比、转速、干涉量和热流密度对刷丝最高温度的影响,揭示了刷式密封的传热机理。结果表明:高温区主要集中在末排刷丝与转子面接触位置,刷丝的最高温度随着压比、转速、干涉量和热流密度增加而增大,其中干涉量对刷丝最高温度的影响最为明显。当干涉量从0.1 mm增至0.7 mm时,刷丝的最高温度上升1.61倍;刷式密封热量的主要来源为刷丝与转子表面摩擦产生的热量,其传热形式包括导热和对流换热,摩擦热量通过导热形式进入刷丝和转子,当刷丝与转子之间的接触力增加时,摩擦热量增大,刷丝的最高温度升高,摩擦热量通过对流换热形式在流体和固体之间进行传递,热量散失主要形式为泄漏气流带走部分热量。      

高旋转雷诺数下预旋进气转-静盘腔流动换热特性

运用RNG k-ε湍流模型对高旋转雷诺数和预旋进口速度下,静盘外缘预旋进气、转盘外缘轴向出流模型的流动和换热过程进行了三维数值模拟,主要研究了冷气流量Cw、旋转雷诺数R ee等参数对转盘对流换热系数和出流口温度分布的影响,并与垂直进气方式进行了对比。研究表明:预旋进气方式与垂直进气相比可降低涡轮叶片冷气入口总温;冷气流量增大以及旋转雷诺数增大均使得转盘平均换热增强;涡轮叶片入口温度随冷气流量增大而降低,随着旋转雷诺数的增大先升高后降低。