删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于CFD的斜盘/滑靴副油膜特性分析

本站小编 Free考研考试/2021-12-25

?油膜特性与柱塞泵性能、可靠性、寿命密切相关,是影响柱塞泵机械效率、容积效率的关键参数。经过良好设计的滑靴结构,能够使斜盘/滑靴副(以下称滑靴副)在工作状态下保持合适的油膜厚度,既不会导致过多的磨损和机械效率降低,也不会引起过大的泄漏、动态不稳定和容积效率降低。除此之外,油膜特性还是影响柱塞泵滑靴副磨损的最重要参数。
国内外关于滑靴副的研究,多集中在分析结构参数对油膜的影响。浙江大学针对某型具体柱塞泵[1, 2],在滑靴副的力/力矩平衡分析基础上,提出了基于雷诺方程的滑靴副油膜厚度的分析计算流程。Hooke和Li[3]介绍了滑靴副油膜厚度的分析流程,通过实验对比分析了几种不同结构尺寸的滑靴对油膜厚度的影响,并提出了滑靴副油膜厚度的计算流程,考虑了滑靴倾覆力矩的影响。Koc和Hooke[4]分析了滑靴副的平衡情况及影响滑靴平衡的因素,在其系列研究中,对滑靴副最小油膜厚度和压力损失的情况进行了分析。
关于滑靴副油膜特性的研究,早期主要采用实验方法。Koc和Hooke[4]通过实验手段对比分析了节流孔尺寸、斜盘倾角及滑靴结构尺寸对滑靴副润滑情况的影响,然后基于实验结论对滑靴副的设计优化给出了建议。Manring等[5]用实验的方法对不同滑靴球窝接触形式进行了相关研究,以分析不同形式对滑靴副油膜的影响。Canbulut[6]也是通过实验的方法研究分析了滑靴表面粗糙度、柱塞阻尼孔直径、柱塞腔压力和转速对滑靴副的泄漏和摩擦力矩的影响。Fisher[7]把解析分析结果与实验结果进行了对比分析,验证了解析方法的适用性。
国内也有****针对滑靴副开展相关实验研究。庄欠伟等[8, 9]开发了滑靴副润滑特性实验台,该实验台可实现高压、高转速和大流量柱塞泵的滑靴副油膜厚度测量。华中科技大学开发了纯水柱塞泵的滑靴副实验台,分析了滑靴副摩擦副特性[10]。北京航空航天大学针对液压柱塞泵滑靴副磨损及影响特性进行了分析与研究[11, 12]
基于实验的滑靴副油膜特性研究具有直观性。然而实验台并不能完全真实地模拟滑靴在柱塞泵中的工作状态,而且一些关键内部变量的测量也非常困难,如油膜厚度、分离力和流场形式等,所以仅依据实验结果并不能全面地分析柱塞泵中滑靴副的油膜特性。浙江大学设想了滑靴副油膜特性实验平台[2],但仍然是通过对比不同工况下的滑靴底面磨损特性,间接验证了解析分析方法的合理性。
滑靴副油膜特性的理论分析方法很好地弥补了实验方法的缺陷,其主要包括解析方法和数值方法。
解析方法的优点是简单直观。文献[7, 13, 14, 15]利用解析方法表达了滑靴副的结构与工况参数对油膜特性的影响。因为解析分析具有显示表达的优势,所以比较容易发现关键因素并进行设计优化。在柱塞泵油膜特性研究领域,普渡大学Maha实验室取得了大量的研究成果[16],基于其开发的CASPAR软件,可以模拟各摩擦副油膜的详细特性。
在采用数值方法对滑靴副油膜特性进行相关的研究中,尤其是油膜问题的研究中,20世纪以来都是使用雷诺方程,并且在大多数情况下,都能得到比较合理的结果[17, 18, 19, 20]。目前,越来越多的研究者开始利用Navier-Stokes方程解决油膜特性仿真问题[21, 22, 23, 24],主要有以下3个方面的原因:①在一定条件下,雷诺方程是无效的,尤其是流体的惯性、紊流和空化现象十分显著的情况[25]。②越来越多油膜特性的研究不再单纯针对摩擦副的压紧区域,而是研究整个润滑系统[26],尤其是在复杂几何形状的油膜润滑特性分析中,雷诺方程由于自身的局限性,很多问题无法适用,而基于Navier-Stokes方程的数值方法则游刃有余[27]。对于滑靴副而言,雷诺方程更多关注的是封油带内的压力、流场的分布。若采用Navier-Stokes方程,则可以直接将滑靴油池、平衡槽和封油带等结构,甚至柱塞的阻尼孔和球铰油池都当作一个整体来仿真。③越来越多基于Navier-Stokes方程的商业仿真软件(ANSYS、FLUINT和COMSOL)与高速计算机的出现,为基于Navier-Stokes方程的计算流体力学(Computational Fluid Dynamic,CFD)仿真在润滑领域中的应用铺平了道路。
本文研究基于三维Navier-Stokes方程和任意拉格朗日-欧拉(Arbitrary Lagrangian-Eulerian,ALE)方法的滑靴副油膜特性的CFD仿真分析方法。首先介绍了油膜厚度的分析计算流程。然后分别阐述了欧拉描述与ALE描述下的控制方程与边界条件,基于不同工况下的仿真结果,依据滑靴受力平衡原则,通过调节油膜厚度得出对应工况下的油膜厚度。着重分析了滑靴副油膜厚度的影响因素,包括柱塞泵转速、温度和阻尼孔入口压力的影响。最后研究了滑靴副油膜特性分析方法,提出一种基于CFD的滑靴副油膜特性分析方法,并结合滑靴副流量守恒原则,给出了滑靴副油膜厚度计算的解析公式。基于此解析公式,可以直观地分析油膜厚度的影响因素,为滑靴的优化设计提供支持。
1 滑靴副油膜特性分析流程本文研究的轴向柱塞泵含有9个柱塞,单个柱塞/滑靴装配件的结构如图 1所示。柱塞与滑靴借助球铰链接。滑靴中心有阻尼孔,柱塞腔内的高压介质通过阻尼孔流入滑靴底面。此外,滑靴底面还包括了凹槽与油池的结构形式。这些凹槽与油池具有增加润滑和保持滑靴稳定的作用。
图 1 柱塞/滑靴装配件 Fig. 1 Piston/slipper assembly
图选项


滑靴副油膜特性分析流程见图 2。在CFD仿真之前,需要明确滑靴副的结构参数和工况参数。然后基于CFD方法进行滑靴副的油膜厚度分析。通过不断调节油膜厚度,基于三维Navier-Stokes方程和仿真工具,计算在特定结构和工况下的滑靴副受力状态,保证滑靴/斜盘处于受力平衡状态。
图 2 滑靴副油膜厚度分析流程 Fig. 2 Oil film thickness analysis procedure of slipper pair
图选项


图 2所述流程适用于特定结构、特定工况下的油膜厚度分析。在CFD仿真结果基础上,结合滑靴副的流量守恒原则,本文提出了基于液动力参数描述方程的油膜厚度解析公式。该解析公式综合考虑了滑靴副的结构参数和工况参数,可以定量地描述油膜厚度与工况参数的函数关系。
2 基于CFD的油膜厚度计算2.1 计算域描述及控制方程本文研究的柱塞泵滑靴副结构如图 3所示。该滑靴含有非封油带、封油带、通油槽和平衡槽等结构。图中:wd为通油槽宽度;dz为柱塞直径;l为阻尼孔长度;d为阻尼孔直径;h为油膜厚度;h1为平衡槽高度;h2为阻尼孔高度;r1、r2、r3r4为滑靴半径。
图 3 滑靴副结构示意图 Fig. 3 Computational domain of slipper structure
图选项


非封油带的引入增加了与斜盘接触的面积,改善了滑靴中心凸起的现象,提高了滑靴副的抗磨能力。封油带的作用是密封液压油和形成油膜。通油槽的作用是增加滑靴的润滑。平衡槽的作用是调整滑靴姿态。此外,滑靴底面中心内置阻尼孔,阻尼孔将柱塞腔内的高压油液引入滑靴油池,油池借助通油槽与平衡槽相连。
研究采用三维Navier-Stokes方程描述滑靴副的流体力学特性,计算域Ω包括阻尼孔、球铰油池、滑靴油池、通油槽、平衡槽、封油带以及非封油带。与雷诺方程仅描述滑靴副压紧面的局部特性相比,Navier-Stokes方程可以更全面地分析惯性力、复杂几何形状对油膜特性的影响。
针对滑靴副的油膜特性进行CFD仿真研究,有如下基本假设:①滑靴副内的流动状态是层流;②滑靴副油膜厚度是一致的;③滑靴副的工况是稳定的;④CFD仿真考虑滑靴的旋转速度。
欧拉描述下的Navier-Stokes方程如式(1)和式(2)[28]所示。
式中:u为流速;η为动力黏度;ρ为介质密度;p为压强;f为离心体积力。
本文考虑了由于柱塞泵高速旋转产生的离心体积力f,f的计算公式为f=ρω2r,其中,ω为滑靴旋转角速度,r为计算单元距离滑靴选择中心的距离。并且引入了介质的体积模量,考虑到了液压介质在高压下的密度变化。此外,式(2)中动力黏度η受温度和压力的综合影响,但在低压(0~50 MPa)范围内,压力对黏度的影响相对较小。本文采用15号航空液压油,只考虑温度对黏度的影响,当温度为40℃时,黏度η取值0.012 75 Pa·s;当温度为80℃时,黏度η取值0.003 434 Pa·s[29]
式(1)和式(2)的计算域为Ω(t)×[0,T](即Ut∈[0,T]Ω(t)×{t}的简写),该计算域随时间变化,用以表征滑靴在斜盘上的旋转运动。
计算域Ω(t)×[0,T]的边界为Г{ГABin,Гout},对应的边界条件如下。
计算域Ω(t)的初始速度:
u(x,y,z,0)=w(x,y,z,0)
滑靴边界ГA的速度:
uA(x,y,z,t)=w(x,y,z,t)
斜盘边界ГB的速度:
uB(x,y,z,t)=0
阻尼孔入口Гin处压力:
pinin(x,y,z,t)·n
滑靴副出口Гout处压力:
poutout(x,y,z,t)·n
其中,w(x,y,z,t)满足:
w(x,y,z,t)=dr/dt,r·w=0,r×w=rω
式中:σinσout分别为滑靴副入口处和出口处的应力张量。
2.2 ALE描述下的计算域与驱动方程ALE方法融合了拉格朗日描述和欧拉描述[30],是处理自由边界和动边界问题的有效方法[27, 31]。由于滑靴副CFD仿真考虑了柱塞泵的旋转等诸多因素,导致滑靴副计算域的边界条件较为复杂。通过ALE方法可将动边界问题转化为静止边界问题,使问题简化。
g(x,y,z,t)是欧拉描述下的方程,ALE描述方法将方程从欧拉描述的坐标系(x,y,z)转换到参考坐标系(xr,yr,zr)中,ALE转换见式(3)、式(4)[27]
式中:v′为迁移速度,即滑靴副的运动速度。
借助转换式(3)和式(4),可将欧拉描述下的Navier-Stokes方程式(1)、方程式(2)转化为ALE描述下的Navier-Stokes方程式(5)、方程式(6)。CFD仿真的计算域Ω也从运动域Ω(t)×[0,T]转换到固定域Ω(0)×[0,T]中。
ALE描述下滑靴副的边界条件如下。
计算域Ω(0)的初始速度:
u(x,y,z,0)=0
滑靴边界ГA的速度:
uA(xr,yr,zr,t)=0
斜盘边界ГB的速度:
uB(x,y,z,t)=wr(xr,yr,zr,t)
阻尼孔入口Гin处压力:
pinin(xr,yr,zr,t)·n
滑靴副出口Гin处压力:
poutout(xr,yr,zr,t)·n
其中,wr(xr,yr,zr,t)满足:
wr(xr,yr,zr,t)=dr/dt
r·wr=0,r×wr=rω
2.3 解算方法本文采用有限元方法求解滑靴副的CFD仿真问题。对于有限元方法,式(5)和式(6)的弱形式是必不可少的。Codina和Zienkiewicz[32]将Navier-Stokes方程式(1)和方程式(2)的弱形式涵盖在一个变分方程中,如式(7)所示。

式中:v和q为试函数。 u(t+,-n)=limε→0+,-u(tn+ε)
因为Navier-Stokes方程是非线性的对流扩散方程,若直接采用Galerkin有限元方法对Navier-Stokes方程离散化得到的方程式(7)是不稳定的,所以需要相关的数值稳定方法,如SUPG、GBS和GLS等方法。本文采用GLS方法建立ALE描述下的Navier-Stokes方程稳定离散形式如式(8)和式(9)所示。
式中:τgls为数值参数;R为残差。Codina和Zienkiewicz[32]给出了τglsR的详细取值方法。
2.4 滑靴副网格划分与CFD仿真图 4展示了滑靴副几何模型的网格划分。滑靴副采用了四面体的网格划分方法。图 4中滑靴副的几何模型中最大的尺寸有十几毫米,最薄的油膜厚度仅有几微米至十几微米。为模拟油膜中介质的流动特性,对垂直于油膜方向的网格施行密集化处理,仅有几微米厚的油膜内部就划分了多层网格。
图 4 滑靴副CFD仿真的几何模型网格划分 Fig. 4 Geometrical mesh of CFD model of slipper pair
图选项


针对滑靴副的CFD仿真,本文采用有限元分析软件COMSOL。COMSOL的CFD模块可以精确地模拟包括基于ALE方法的动边界流动在内的多种流动。设滑靴副阻尼孔的入口压力为21 MPa,转速为4 000 r/min,介质黏度为0.012 75 Pa·s,油膜厚度取值11 μm时,滑靴副压力的三维分布见图 5(a),滑靴底面的压力分布见图 5(b)。此时有限元方法的自由度为259 115。关于解算时间,在给定工况下,COMSOL软件的仿真时间为1 602 s。仿真机的配置为:CPU为Intel (R) Core (TM) i5-2410M,内存为4.00 GB,操作系统为win7家庭普通版。
图 5 滑靴副压力分布 Fig. 5 Pressure distribution of slipper pair
图选项


2.5 油膜厚度的计算结果根据滑靴副的CFD仿真结果可以得到压力的分布。图 6展示了2种工况下滑靴副C-C剖面(见图 3)的油膜压力分布。图 6的最大特点是转速、阻尼孔入口压力和油膜厚度对滑靴副的压力分布形式都不会产生显著的影响,改变的只是最大压力。需要指出的是,从图 6中可看出C-C剖面中存在压力差,这是旋转体积力在滑靴副中分布不均产生的结果。当转速Rs为3 000 r/min,阻尼孔入口压力为10 MPa时,压力差约为0.31 MPa。仿真结果也与Kumar等[23]的仿真结果类似,证明了本文提出的描述方程及数值方法的有效性。
图 6 工况对滑靴副油膜压力分布的影响 Fig. 6 Effect of running condition on oil film pressure distribution of slipper pair
图选项


针对一种工况,基于CFD仿真结果可以分析滑靴副的受力状态,通过不断调节油膜厚度,并保证滑靴/斜盘处于受力平衡,可得到这种工况下的油膜厚度。滑靴主要受到压紧力和分离力的作用,其中分离力由滑靴副内的高压介质产生,可从CFD仿真中得到。压紧力Fb主要来源于柱塞底部受到的液压力,如式(10)所示。本文并未考虑斜盘倾角γ对油膜特性的影响,取γ=0°。
若改变工况,则需要重新调整油膜厚度,使滑靴/斜盘达到受力平衡。图 7展示了24种工况(包括2个温度等级、3个压力等级和4转速等级)下压紧力Fb与CFD仿真得到的分离力。对比图 7中的24组分离力与3种压力状态下的压紧力Fb,可知两者具有一致性。
图 7 24种工况下的滑靴分离力 Fig. 7 Separation force on slipper under 24 different running conditions
图选项


图 8展示了CFD仿真得到的24种工况下的滑靴副油膜厚度。由图 8可以看出:①随转速变快,油膜变厚。这是因为转速增大引起滑靴副的动压增大,需滑靴副静压承担的分离力减少,所以油膜增厚。②阻尼孔入口压力pin增大导致油膜厚度变薄。增大入口压力pin使得式(10)中滑靴压紧力Fb增大,将滑靴副油膜压薄。③温度升高导致油膜厚度变薄。这是因为温度升高致使介质的黏度下降,降低了油膜的承载能力。
图 8 油膜厚度与工况关系 Fig. 8 Relationship between oil film thickness and running condition
图选项


3 油膜厚度的影响因素分析3.1 液动力参数在文献[1, 2, 3, 4, 5, 6, 7]中,分析滑靴副润滑特性时都提及了液动力参数,但是都没有对这个参数的物理含义进行分析。液动力参数与滑靴油池压力、介质黏度、柱塞泵转速、滑靴结构尺寸以及滑靴副的油膜厚度均相关,可以描述滑靴的综合工况。液动力参数G可以用来描述滑靴副的油膜特性,如式(11)所示。
式中:U为滑靴斜盘相对运动速度;ps为滑靴油池压力。
针对图 8展示的24种工况下的油膜厚度结果,结合式(11)可以计算液动力参数G,计算结果见图 9。
图 9 24种工况下液动力参数G Fig. 9 Hydrodynamic parameter G under 24 different running conditions
图选项


从图 9可以看出,在固定滑靴结构参数情况下,温度(介质黏度)对液动力参数G的影响最为显著,而阻尼孔入口压力和转速的影响可以忽略。这也与文献[4]的结论一致。
3.2 基于液动力参数的油膜厚度分析在温度不变的情况下,结合液动力参数表达式(11)和滑靴副流量守恒原则,可以得到滑靴副油膜厚度的解析公式。
滑靴油池/柱塞腔间阻尼孔的流量:
滑靴与斜盘可变缝隙的流量:
流量Q1Q2相等,得到
将式(14)代入到式(11)中,可得到滑靴副油膜厚度h的解析表达式。
式(15)通过流量守恒原则考虑了滑靴副阻尼孔的作用,又通过液动力参数G考虑了诸如温度、柱塞泵转速和滑靴半径等参数,几乎囊括了滑靴副CFD仿真的所有关键影响因素。
图 10对比了40℃下CFD仿真与解析式(15)计算的滑靴副油膜厚度h。其中,解析式(15)中的液动力参数G参考了图 9中的结果,取值为0.3。观察图 10可知,解析式(15)与CFD仿真得到的油膜厚度相符得很好。
图 10 40℃时CFD与式(15)计算的油膜厚度h Fig. 10 Oil film thickness h obtained from CFD and Eq. (15) (temperature=40℃)
图选项


4 结 论1) 基于CFD和数值仿真方法开展柱塞泵滑靴副油膜特性分析,可以综合考虑滑靴副结构形式、尺寸参数和环境工况对油膜的综合影响,具有全面性。其中基于液动力参数的油膜厚度分析方法,通过解析手段定量直观描述油膜厚度与工况参数之间的函数关系。
2) 在滑靴副结构尺寸固定的情况下,滑靴副油膜厚度与环境工况之间的关系可以通过液动力参数进行定量描述。针对特定的滑靴副,液动力参数只与介质温度(黏度)相关。
3) 鉴于通过实验手段直接测量滑靴副油膜特性难度较大,本文关于油膜厚度的分析结论可以通过对比不同工况下的磨损情况(和油膜厚度密切相关)进行验证。
参考文献
[1] 张斌. 轴向柱塞泵的虚拟样机及油膜压力特性研究[D].杭州:浙江大学,2009:75-97. ZHANG B.Study on virtual prototype and pressure characteristics of oil film for axial piston pump[D].Hangzhou:Zhejiang University,2009:75-97(in Chinese).
Cited By in Cnki (43)
[2] 李迎兵. 轴向柱塞泵滑靴油膜特性研究[D].杭州:浙江大学,2011:23-35. LI Y B.Study on characteristics of oil film for slipper pair of axial piston pumps[D].Hangzhou:Zhejiang University,2011:23-35(in Chinese).
Cited By in Cnki (10)
[3] HOOKE C J,LI K Y.The lubrication of slippers in axial piston pumps and motors-the effect of tilting couples[J].Journal of Mechanical Engineering Science,1989,203(53):343-350.
Click to display the text
[4] KOC E,HOOKE C J.Considerations in the design of partially hydrostatic slipper bearings[J].Tribology International,1997,30(11):815-823.
Click to display the text
[5] MANRING N D,WRAY C L, DONG Z.Experimental studies on the performance of slipper bearings within axial piston pumps[J].Journal of Tribology,2004,126(3):511-518.
Click to display the text
[6] CANBULUT F. The experimental analyses of the effects of the geometric and working parameters on the circular hydrostatic thrust bearings[J].JSME International Journal,2005,48(4):715-722.
Click to display the text
[7] FISHER M J. A theoretical determination of some characteristics of a tilted hydrostatic slipper bearing[R].[S.l.]:British Hydromechanics Research Association,1962.
[8] 庄欠伟,周华,艾青林.轴向柱塞泵滑靴副润滑特性实验台的研制[J].机床与液压,2005(3):113-115. ZHUANG Q W,ZHOU H,AI Q L.Development of experiment rig for lubricate film of the slipper pad in axial piston pumps[J].Machine Tools & Hydraulics,2005(3):113-115(in Chinese).
Cited By in Cnki (11)
[9] 艾青林,周华,张增猛,等.轴向柱塞泵配流副与滑靴副润滑特性试验系统的研制[J].液压与气动,2004(11):22-25. AI Q L,ZHOU H,ZHANG Z M, et al.Development of lubricating characteristics testing system of friction pair in axial piston pump[J].Chinese Hydraulics & Pneumatics,2004(11):22-25(in Chinese).
Cited By in Cnki (30)
[10] NIE L, HUANG G H,LI Y R.Tribological study on hydrostatic slipper bearing with annular orifice damper for water hydraulic axial piston motor[J].Tribology International,2006,39(11):1342-1354.
Click to display the text
[11] 马纪明,李齐林,任春宇,等.轴向柱塞泵/滑靴副润滑磨损的影响因素分析[J].北京航空航天大学学报,2015,41(3):405-410. MA J M,LI Q L,REN C Y,et al.Influence factors analysis on wear of hydraulic axial piston pump/slipper pair[J].Journal of Beijing University of Aeronautics and Astronautics,2015,41(3):405-410(in Chinese).
Cited By in Cnki (2)
[12] MA J M,CHEN J,LI J.Wear analysis of swash plate/slipper pair of axis piston hydraulic pump[J].Tribology International,2015,90:467-472.
Click to display the text
[13] IBOSHI N,YAMAGUCHI A.Characteristics of a slipper bearing for swash plate type axial piston pumps and motors,theoretical analysis[J].Bulletin of the JSME,1982,25(428):1921-1930.
[14] IBOSHI N. Characteristics of a slipper bearing for swash plate type axial piston pumps and motors,design method for a slipper with a minimum power loss in fluid lubrication[J].Bulletin of the JSME,1986,129(254):2529-2538.
[15] BERGADA J M,KUMAR S,DAVIES D L,et al.A complete analysis of axial piston pump leakage and output flow ripples[J].Applied Mathematical Modeling,2012,36(4):1731-1751.
Click to display the text
[16] SCHENK A,IVANTYSYNOVA M.A transient fluid structure interaction model for lubrication between the slipper and swash plate in axial piston machines[C]//Proceedings of the 9th International Fluid Power Conference(9IFK),2014,1:398-409.
[17] KAZAMA T,YAMAGUCHI A.Application of a mixed lubrication model for hydrostatic equipment[J].Journal of Tribology,1993,115(4):686-691.
Click to display the text
[18] FANG Y,SHIRAKASHI M.Mixed lubrication characteristics between the piston and cylinder in hydraulic piston pump-motor[J].Journal of Tribology,1995,117(1):80-85.
Click to display the text
[19] BERGADA J M,WATTON J,KUMMAR S.Pressure,flow,force,and torque between the barrel and port plate in an axial piston pump[J].Journal of Dynamic Systems Measurement and Control,2008,130(1):141-148.
Click to display the text
[20] 王彬,周华,杨华勇.轴向柱塞泵配流副油膜试验原理及控制特性[J].机械工程学报,2009,45(11):113-118. WANG B,ZHOU H,YANG H Y.Principle of oil film test on port pair of axial piston pump and control characteristics[J].Journal of Mechanical Engineering,2009,45(11):113-118(in Chinese).
Cited By in Cnki (7)
[21] BRAIDIC-MITIDIERI P,GOSMAN A D,IOANNIDES E,et al.CFD analysis of a low friction pocketed pad bearing[J].Journal of Tribology,2005,127(4):803-812.
Click to display the text
[22] HOUZEAUX G,CODINA R.A finite element method for the solution of rotary pumps[J].Computers & Fluids,2007,36(4):667-679.
Click to display the text
[23] KUMAR S,BERGADA J M,WATTON J.Axial piston pump grooved slipper analysis by CFD simulation of three-dimensional NVS equation in cylindrical coordinates[J].Computers & Fluids,2009,38(3):648-663.
Click to display the text
[24] KUMAR S. CFD analysis of an axial piston pump[D].Barcelona:Universidad Politecnica de Catalunya,2010:45-86.
[25] ALMQVIST T,ALMQVIST A,LARSSON R.A comparison between computational fluid dynamic and Reynolds approaches for simulating transient EHL line contact[J].Tribology International,2004,37(1):61-69.
Click to display the text
[26] GUO Z,HIRANO T,KIRK R G.Application of CFD analysis for rotating machinery,Part 1:Hydrodynamic,hydrostatic bearings and squeeze film damper[J].ASME,2003,4(1):651-659.
Click to display the text
[27] DUARTE F,GORMAZ R,NATESAN S.Arbitrary Lagrangian-Eulerian method for Navier-Stokes equations with moving boundaries[J].Computer Methods in Applied Mechanics and Engineering,2004,193(45-47):4819-4836.
Click to display the text
[28] 张兆顺. 流体力学[M].北京:清华大学出版社,2006:53-55. ZHANG Z S.Fluid mechanics[M].Beijing:Tsinghua University Press,2006:53-55(in Chinese).
[29] 李玉琳. 液压元件与系统设计[M].北京:北京航空航天大学出版社,1991:13-14. LI Y L.Hydraulic components and system design[M].Beijing:Beihang University Press,1991:13-14(in Chinese).
[30] HIRT C W,AMSDEN A A,COOK J L.An arbitrary Lagrangian-Eulerian computing method for all flow speeds[J].Journal of Computational Physics,1974,14(2):227-253.
Click to display the text
[31] SOULAIMANI A,SAADB Y.An arbitrary Lagrangian-Eulerian finite element method for solving three-dimensional free surface flows[J].Computer Methods in Applied Mechanics and Engineering,1998,162(1-4):79-106.
Click to display the text
[32] CODINA R,ZIENKIEWICZ O C.CBS versus GLS stabilization of the incompressible Navier-Stokes equations and the role of the time step as stabilization parameter[J].Communications in Numerical Methods in Engineering,2002,18(2):99-112.
Click to display the text


相关话题/计算 结构 动力 实验 油膜

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于小扰动理论的桨叶叶素气动载荷计算方法
    ?直升机依靠旋翼桨叶产生推进力、升力和操纵力,靠尾桨产生偏航力矩[1],可以垂直起降而不受地形限制,无需大面积机场,可以向任意方向进行机动飞行。由于旋翼的气动性能直接影响到直升机的整体飞行性能,因此气动性能的准确模拟就显得尤为重要。研究直升机旋翼的非定常气动特性,就是研究旋翼诱导速度场、非定常气动力 ...
    本站小编 Free考研考试 2021-12-25
  • 磁悬浮控制力矩陀螺热-结构耦合分析与研究
    ?磁悬浮控制力矩陀螺(MSCMG)是空间飞行器实现高精度姿态机动通常采用的惯性执行机构,具有转速高、输出力矩大以及寿命长等特点。磁悬浮控制力矩陀螺转子依靠磁轴承支承[1,2],并且在高速电机的驱动下持续高速旋转,陀螺在工作过程中电机与磁轴承会产生铁损与铜损,这些损耗最终转化为热量传递到陀螺其他结构部 ...
    本站小编 Free考研考试 2021-12-25
  • 确定飞机结构WFD平均行为的寿命升降法
    自1988年Aloha航空事故[1]以来,广布疲劳损伤(WidespreadFatigueDamage,WFD)问题引起了学术界和航空界的广泛关注。WFD是指多个相似或相同细节中同时出现尺寸和密度足够的裂纹,使结构的剩余强度不再满足美国联邦航空条例§25.571(b)的要求[2]。WFD包括2类情况 ...
    本站小编 Free考研考试 2021-12-25
  • 前置隔板对圆柱绕流影响的实验分析
    在一定雷诺数下,直匀流在通过圆柱等钝体后会产生周期性的卡门涡街,这种非对称脱落的旋涡可能会导致圆柱失稳甚至诱发振动[1],1940年美国华盛顿州的一处悬索桥在风速不大的情况下坍塌就是典型的案例。为了防止类似事故的发生,研究者采用各种控制方法,对流动分离进行人为干预。目前,流动控制主要分为被动控制和主 ...
    本站小编 Free考研考试 2021-12-25
  • 利用符号计算方法研究生物系统全时滞稳定性
    时滞微分系统的稳定性研究在理论和应用上都有其重要的意义,特别地,从控制理论的角度看,生物系统全时滞稳定即表明该系统对于时滞具有很好的鲁棒性和可靠性.长期以来,人们一直致力于寻找时滞微分系统全时滞稳定的代数判据,已取得了不少进展.秦元勋在文献[1]中第1次将单滞后多维系统的全时滞稳定判据由超越形式的检 ...
    本站小编 Free考研考试 2021-12-25
  • 结构偏差对二维连续地月载荷转移系统动力学影响
    绳系卫星作为一种新型的卫星概念,经过数十年的研究历程,在航天领域扮演着十分重要的角色.绳系卫星具有诸多特点[1,2]:能够实现动量(能量)的传递[3,4,5],能够借助于电动绳的原理为卫星系统提供可再生的能源,能够实现对太空垃圾的捕获等[6,7].动量交换作为绳系的一个重要功能,其基本思想就是通过系 ...
    本站小编 Free考研考试 2021-12-25
  • 航天员受银河宇宙线辐射的剂量计算
    空间中的质子、电子、重离子、中子和X射线等辐射是威胁航天员健康和安全的环境要素.近20年,国际和我国的载人航天事业均有了很大发展,更多航天员参与空间活动,航天员在空间停留的时间也更长,接受越来越多的空间辐射照射.不同于以X射线、γ射线和中子为主的地面辐射,空间辐射含有高传能线密度(LET)的银河宇宙 ...
    本站小编 Free考研考试 2021-12-25
  • 周期结构电磁特性在高频真空器件中的应用
    在行波类真空电子器件中,电子注和电磁波之间产生持续互作用的一个条件就是电磁波的相速小于光速,通常有两种方法来降低相速,一种是采用介质材料加载的波导,另一种就是采用周期性的金属或者介质结构.微波真空器件通常采用金属周期结构,其传输损耗比较小,可以维持较高的射频电场,还可以接收高频散焦等杂散电子注电流. ...
    本站小编 Free考研考试 2021-12-25
  • 基于可编程控制网关的一体化网络体系结构
    随着Internet、移动通信和卫星通信技术的高速发展,涌现出了大量不同类型的通信网络,其中典型的网络有互联网、地面无线通信网络以及各种专用网络.其中互联网和地面无线通信网络以IP(InternetProtocol)技术为基础逐步走向融合.而非IP技术的专用网络,如战术数据链网络、基于国际空间数据系 ...
    本站小编 Free考研考试 2021-12-25
  • 基于自由尾迹/面元法的旋翼下洗干扰计算和直升机配平
    受旋翼涡系的影响,直升机工作在十分复杂的流场当中,旋翼的下洗流对机身振动、桨叶气动载荷和直升机的飞行性能都有十分重要的影响.长期以来,科研人员从数值计算和试验两方面对涡动力学方法进行了大量的研究[1,2,3],对旋翼/机身的复杂气动干扰的理解也愈发加深.当直升机在吊挂飞行、舰载和沙漠等复杂条件起降时 ...
    本站小编 Free考研考试 2021-12-25