一种考虑倾斜与热的静压油垫性能计算方法技术

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

【技术实现步骤摘要】

本专利技术属于静压转台分析领域，涉及一种定量供油的转台支承油垫承载力、刚度及阻尼的计算方法，该方法运用有限差分法计算并考虑了摩擦热及倾斜的作用。
技术介绍
静压转台(Hydrostatic Rotary Table)用有压力的流体使有相对运动的两个表面分开并借助流体静压来承载。由于运动副之间完全被油膜隔开，所以运动副间的摩擦力大大减小，同时其承载能力、运动精度与寿命却大大提高。正因为液体静压支承的诸多优点，所以它广泛的应用于重型机床并成为其关键部件之一。静压油垫是静压转台支承系统的关键部件，通常在其设计时认为其与导轨面永远保持平行，然而这并不现实，事实上由于加工装配误差、转台的局部变形及偏载的影响油垫与导轨面往往难以保持平行。油液的发热也是影响油垫性能的关键因素，它使油液的粘度降低同时也会使转台产生热变形从而影响转台的精度。所以在转台的设计及使用过程中需要充分预估倾斜及热对油垫承载性能的影响。
技术实现思路
本专利技术的目的是提供一种定量供油的转台支承油垫承载力、刚度及阻尼的计算方法，该方法运用有限差分法计算并考虑了摩擦热及倾斜的作用。，该方法首先简化N-S方程与能量方程并推导二维雷诺方程，再对各方程进行无量纲处理，之后计算油垫的承载力，刚度与阻尼，最后基于有限差分法使用Matlab编写计算程序得到油垫的各项性能参数。为实现上述目的，本专利技术采用的技术方案为一种考虑倾斜与热的静压油<br>垫性能计算方法，该方法的实现步骤如下，1、该方法首先简化N-S方程与能量方程并推导二维雷诺方程确定其边界条件。根据油垫的倾角计算倾斜油垫各处的时油膜厚度。2、对能量方程、雷诺方程、温粘方程方程进行无量纲处理。3、计算油垫的承载力、刚度与阻尼性能参数。4、基于有限差分法使用Matlab编写计算程序得到油垫的各项性能参数。与现有技术相比，本专利技术具有如下有益效果。本专利技术的特点在于使用简化的N-S方程推导雷诺方程，能量方程与雷诺方程通过粘温方程与密温方程相互耦合。耦合的方程使用有限差分法求解，油垫的承载力通过承载面积分得到，而刚度与阻尼通过对承载力求导得到，本专利技术提供的方法可为转台的使用及设计提供指导。附图说明图1导轨相对与油垫倾斜的模型图。图2支承油垫结构简图。图3差分法计算流程图。图4油垫各项性能随倾角的变化曲线。图5油垫平均温度随倾角的变化曲线。具体实施方式本专利技术实施一种静压转台偏载时承载力计算分析方法，下面结合附图，对本专利技术的实施进行具体说明。图1为导轨相对与油垫倾斜的模型图，油垫为圆形阶梯油垫由定量泵向其供油供油量为Q0，以油垫的中心点为原点建立极坐标系r,为坐标参数。转台在承受偏载或有局部变形时油膜厚度在油垫上的分布将不再均匀，此时转台的导轨面与油垫面有一个倾角θ(如图1所示)。图2为支承油垫结构简图。图中R1为油垫的内径，R2为油垫的外径，hs为油膜厚度。步骤(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油膜厚度的计算图2为支承油垫结构简图，R1,R2为油垫的内外径，hs是支承油垫的油膜厚度，θ为转台导轨面相对与油垫的倾斜角度，Q0是油垫的供油量。这样当油垫倾斜时油垫封油边的实际油膜厚度为：步骤(2)方程的无量纲处理上式中r,均为极坐标系的坐标参数。雷诺方程及能量方程需要进行无量纲处理，于是令h0＝hs,R‾=rR2,H‾=hsh0,ρ‾=ρρ0,η‾=ηη0,p‾=pp0,z‾=zh0,vz‾=vzh0,T‾=TT0,]]>式中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本文档来自技高网...

【技术保护点】
一种考虑倾斜与热的静压油垫性能计算方法，其特征在于：该方法的实现步骤如下，1、该方法首先简化N‑S方程与能量方程并推导二维雷诺方程确定其边界条件；根据油垫的倾角计算倾斜油垫各处的时油膜厚度；2、对能量方程、雷诺方程、温粘方程方程进行无量纲处理；3、计算油垫的承载力、刚度与阻尼性能参数；4、基于有限差分法使用Matlab编写计算程序得到油垫的性能参数。

【技术特征摘要】
1.一种考虑倾斜与热的静压油垫性能计算方法，其特征在于：该方法的
实现步骤如下，
1、该方法首先简化N-S方程与能量方程并推导二维雷诺方程确定其边界条件；
根据油垫的倾角计算倾斜油垫各处的时油膜厚度；
2、对能量方程、雷诺方程、温粘方程方程进行无量纲处理；
3、计算油垫的承载力、刚度与阻尼性能参数；
4、基于有限差分法使用Matlab编写计算程序得到油垫的性能参数。
2.根据权利要求1所述的一种考虑倾斜与热的静压油垫性能计算方法，
其特征在于：
步骤(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向

\t的流速；对(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,R‾=rR2,H‾=hsh0,ρ‾=ρρ0,η‾=ηη0,p‾=pp0,z‾=zh0,vz‾=vzh0,T‾=TT0,]]>式
中h0为初始油膜厚度，ρ0及η0是常温时油液的粘度与密度，p0为油腔压力，vz为导轨挤压油膜的速度；同时定义平均密...

【专利技术属性】
技术研发人员：刘志峰，湛承鹏，李小燕，马雷，李彦虎，唐振勇，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11