[发明专利]一种考虑倾斜与热的静压油垫性能计算方法

摘要

本发明公开了一种考虑倾斜于热的静压油垫各项性能计算方法，该方法考虑了倾斜与热的影响并运用有限差分法计算油垫的压力分布及温度分布。其具体实现流程为。该方法首先简化N‑S方程与能量方程并推导二维雷诺方程，再对各方程进行无量纲处理，之后计算油垫的承载力，刚度与阻尼，最后基于有限差分法使用Matlab编写计算程序得到油垫的各项性能参数。本发明的特点在于使用简化的N‑S方程推导雷诺方程，能量方程与雷诺方程通过粘温方程与密温方程相互耦合。耦合的方程使用有限差分法求解，油垫的承载力通过承载面积分得到，而刚度与阻尼通过对承载力求导得到，本发明提供的方法可为转台的使用及设计提供指导。

主权项

一种考虑倾斜与热的静压油垫性能计算方法，其特征在于：该方法的实现步骤如下，(1)该方法首先简化N‑S方程与能量方程并推导二维雷诺方程确定其边界条件；根据油垫的倾角计算倾斜时油垫各处的油膜厚度；(2)对能量方程、雷诺方程、温粘方程方程进行无量纲处理；(3)计算油垫的承载力、刚度与阻尼性能参数；(4)基于有限差分法使用Matlab编写计算程序得到油垫的性能参数；步骤(1)，相关方程的推导1.1雷诺方程的建立油膜比较薄，静压油的流态为层流，油液的密度仅有少量改变，在径向油液的散热以对流为主那么对N‑S方程进行简化有：∂p∂r-η∂2ur∂z2=0---(2)]]>∂p∂z=0---(4)]]>上式中(1)为连续性方程，(2)，(3)，(4)分别为简化后的r向，向及z向动量方程；式中ρ为油液密度，η为油液粘度，ur,vz分别为油液在r向，向及z向的流速；对(3)式两边对z积分并代入边界条件z＝0,z＝hs,有：(2)式对z进行积分并代入边界条件z＝0,ur＝0；z＝hs,ur＝0；有:ur=z2-zhs2η∂p∂r---(6)]]>将(5)、(6)式带入(1)式并对z积分可得到雷诺方程：1.2能量方程油液由摩擦产热，r向和向的散热主要依靠对流，z向散热以传导为主，则可有简化的能量方程为：同时温粘方程与温密方程为：η=η0e-β(T-T0)---(9)]]>ρ＝ρ0(1‑λ(T‑T0))  (10)1.3油膜厚度的计算R1,R2为油垫的内外径，hs是支承油垫的油膜厚度，θ为转台导轨面相对与油垫的倾斜角度，Q0是油垫的供油量；这样当油垫倾斜时油垫封油边的实际油膜厚度为：步骤(2)方程的无量纲处理上式中r,均为极坐标系的坐标参数；雷诺方程及能量方程需要进行无量纲处理，于是令h0＝hs,式中h0为初始油膜厚度，ρ0及η0是常温时油液的密度与粘度，p0为油腔压力，vz为导轨挤压油膜的速度；同时定义平均密度与平均粘度为则无量纲的雷诺方程及油膜计算方程为：无量纲的周向及径向流速为：式中的ur与将按照(5)和(6)式计算；无量纲的能量方程为：无量纲的温密方程及粘温方程为：η‾=e-βT0(T‾-1)---(15)]]>ρ‾=(1-λT0(T‾-1))---(16)]]>雷诺方程(13)与能量方程(14)将使用有限差分法求解，求解的边界条件为：0≤R‾≤R1R2,p‾=1;R‾=1,p‾=0;0≤R‾≤R1R2,T‾=1;R‾=1,∂T‾∂R‾=0;---(17)]]>步骤(3)各项性能的计算对油液径向流速进行积分得到支承油垫的流量为：所以支承油垫的油腔压力为：对油膜压力分布进行积分就得到油垫的承载力为：承载力对油膜厚度求微分就得到油垫的刚度为：Kz=∂wfh0∂H‾---(20)]]>承载力对油垫的挤压速度vz求导就得到油垫的阻尼为：C=∂wfh0∂vz‾---(21)]]>步骤(4)计算流程能量方程与雷诺方程将使用有限差分法求解，收敛的条件为相对误差ζ＜10‑6；而相对误差的计算式为(1)初始化参数包括设置输入参数h0,ρ0,η0,R1,R2,T0及网格的大小nr,nz；(2)计算油垫各处的油膜厚度；(3)计算油膜的压力分布(4)计算油液的径向速度和周向速度(5)根据方程(14)计算油膜各处的温度分布(6)根据(15)，(16)式计算新的油液密度与粘度并更新平均油液密度与粘度(7)计算相对误差ζ；(8)重复第二步到第七步直到相对误差ζ＜10‑6；(9)迭代计算结束，计算油腔压力与承载力参数；上述计算过程将使用matlab编程实现；油垫的参数及油液的参数为：R2＝175mm，h0＝0.1mm，Q0＝1×10‑4m3/s，ρ0＝872kg/m3，η0＝0.091pa·s，T0＝22℃，c＝1845.0J/kgK，β＝0.0372/K，λ＝0.00064737/K，k＝0.145W/(mK)；网格大小设定为根据这些参数按照步骤(4)所述的计算流程计算油膜的压力场与温度场进而计算油垫的承载力、刚度与阻尼。