CN104537189B - 一种静压转台运动误差建模及计算方法 - Google Patents
一种静压转台运动误差建模及计算方法 Download PDFInfo
- Publication number
- CN104537189B CN104537189B CN201510028962.2A CN201510028962A CN104537189B CN 104537189 B CN104537189 B CN 104537189B CN 201510028962 A CN201510028962 A CN 201510028962A CN 104537189 B CN104537189 B CN 104537189B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msubsup
- turntable
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Sliding-Contact Bearings (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
一种静压转台运动误差建模及计算方法,该方法包括基于雷诺方程的单油垫理论模型计算出预压油垫及支承油垫在考虑运动误差时的压力分布pi(r,)及py(r,),然后后对压力分布进行积分,得到预压油垫及支承油垫的承载力Fy,Fi和Fj;将单个油垫的油膜厚度简化为油垫中心的油膜厚度,考虑转台的加工轮廓误差,并将其拟合成函数,通过各油垫油膜厚度之间的几何关系确定各油垫油膜厚度及速度;建立转台的动力学平衡方程;在Matlab软件中用龙格库塔法求解转台的动力学方程从而得到转台转速、加工轮廓误差等主要参数对运动误差的影响,为静压转台的加工精度提出一个合理的范围。
Description
技术领域
本发明属于静压转台分析领域,涉及一种双圈油垫支承的静压转台运动误差建模与仿真计算方法,更具体是一种考虑了转台转速、加工轮廓误差的静压转台运动误差的建模与仿真计算方法。
背景技术
静压转台(Hydrostatic Rotary Table)是用有压力的流体使有相对运动的两个表面分开并借助流体静压来承载。由于运动副之间完全被油膜隔开,运动副间的摩擦力大大减小,同时其承载能力、运动精度与寿命却大大提高。正因为液体静压支承的诸多优点,所以它广泛的应用于重型机床并成为其关键部件之一。
对于重型静压转台,由于其尺寸较大,就目前现有的加工条件来说,在重型静压转台加工的过程中,转台的加工误差是难以控制到较低水平的,所以其加工误差对转台的运动误差的影响不能忽略。目前,国内外学者对静压支承基础理论及静压转台的一般设计方法进行了广泛而深入的研究,但对转台加工误差的研究并不充分,加工误差及加工误差对运动误差的影响,尚无完备的理论研究来分析。
发明内容
本发明的目的是提供一种静压转台的运动误差的建模方法及计算方法,用以分析加工误差对运动误差的影响。该方法首先计算转台支承油垫及预压油垫的动态承载力,然后建立转台的整体动力学模型,最后通过数值方法计算转台转速、加工轮廓误差等主要参数对运动误差的影响。
本发明是采用以下技术手段实现的:
1、忽略转台自身的变形,即近似认为转台为刚体,同时考虑转台的加工轮廓误差,并将加工轮廓误差拟合成函数,对于转台的导轨面而言其表面的波动波长一般大于转台的直径。所以忽略转台倾角引起的油膜厚度不均匀对油垫承载能力的影响。
2、基于雷诺方程的单油垫理论模型计算出预压油垫及支承油垫在考虑加工轮廓误差时的压力分布及之后对压力分布进行积分得到预压油垫及支承油垫的承载力Fy,Fi和Fj;
3、转台在外载作用时其所有受力及力矩应当平衡,那么根据力平衡及力矩平衡原则可以建立转台的动力学平衡方程;
4、考虑转台的加工轮廓误差,通过各油垫油膜厚度之间的几何关系确定各油膜厚度及速度;
5、转台的动力学方程为二阶非线性微分方程,根据数值分析的相关理论可以用龙格库塔法在matlab软件中编写方程的求解程序最终得到方程的数值解,也就是求出转台转速、加工轮廓误差等主要参数对转台运动误差的影响。
本发明的特点在于基于雷诺方程的单油垫理论模型,建立静压转台的动力学模型并用数值方法求解出求出转台转速、加工轮廓误差等主要参数对运动误差的影响,模型中考虑了转台加工轮廓误差的情况,也考虑了转台转速的影响。发明内容包括四部分。在第一部分中计算出转台各油垫的承载力,在第二部分中,建立转台的平衡方程,在第三部分中计算各油垫的油膜厚度,在第四部分中,采用Matlab软件计算转台转速、加工轮廓误差等主要参数对运动误差的影响。最后通过实例论证本发明提出的模型。
通过下面的描述并结合附图说明,本发明会更加清晰,附图说明用于解释本发明方法及实施例。
附图说明
图1静压转台结构简图
图2液体静压转台两圈导轨轮廓加工误差示意图
图3液体静压转台受力示意图
图4支承油垫(圆形)简图
图5预压油垫(环形)简图
图6转台运动误差程序计算流程图
图7本发明实例中转台转速对运动精度的影响
图8本发明实例中转台内外圈油垫处的波动比值随波长的变化规律
具体实施方式
本发明实施一种考虑加工轮廓误差的静压转台运动误差建模及计算方法,下面结合附图,对本发明的实施进行具体说明。
图1为静压转台结构简图,转台由转台,基座,支承油垫,预压油垫、径向轴承组成,其中支承油垫为双圈圆形油垫,预压油垫为环形预压油垫,在转台的中心部位安装有一个径向轴承。转台各油垫都由定量泵供油,由于加工的原因转台的导轨面会高低不平,用正弦函数来等效描述转台导轨面的起伏。转台的自重为G,其中支承油垫的供油量为Q0,预压油垫的供油量为Q1,径向油垫的供油量为Q2。
图2为液体静压转台两圈导轨轮廓加工误差示意图,内外圈波长呈放射状,两者之比等于导轨半径之比。
步骤(1),相关公式推导
基于雷诺方程的单油垫理论模型,推导出转台支承油垫及预压油垫承载力公式,而后列出转台平衡方程。
1、建立单油垫理论模型
在单油垫内,任意割取半径为dr,高度为dz的微球体,定义该微区中心点坐标为(r,z),在某一瞬时t时,过此微元体中心点的两个分速分别为ur、vz。p为油膜压力,η为粘度系数。根据量纲分析可以将极坐标系下的连续性方程和N-S方程简化为:
由(3)式可认为p与z无关,所以对(2)式两边对z进行二次积分并代入边界条件:z=h,ur=0;z=0,ur=0可得
将(4)式带入(1)式并对z积分,再代入边界条件z=h,ur=0,z=0,ur=0,vz=0可以得到雷诺方程为:
分析对象中,支承油垫为圆形油垫,预紧油垫为环形油垫。
对于圆形油垫有:
r=R1,p=p0;r=R2,p=0;Q(R1)=Q0 (7)
根据(5)、(6)、(7)式可解得油腔压力p0和封油边压力分布p(r):
那么圆形油垫的承载能力为:
即:
对于环形预压压油垫:
对于环形油垫有边界条件为:
r=RC1,p=0;r=RC2,p=p0;r=RC3,p=p0;r=RC4,p=0;Q1=-Q(RC2)+Q(RC3) (11)
根据(5)、(6)、(11)式可解得油腔压力p0和封油边压力分布p(r):
当r∈(RC1,RC2)时压力分布为:
当r∈(RC3,RC4)时压力分布为:
那么根据(12)(13)(14)式可以计算环形油垫的承载力为:
即:
2、转台力平衡方程
图3为转台受力简图,转台由两圈油垫共24个油垫支承,对于第一圈油垫,其支承力分别为Fi(i=1~8),它们与加力点和转台中心的连线的夹角分别为油垫中心到转台中心的距离为RL,对于第二圈支承油垫,其支承力分别为Fj(j=1~16),它们与加力点和转台中心的连线的夹角分别为油垫中心到转台中心的距离为Rs,预压油垫的支承力为Fy。转台的半径为RT,转台的转速为ω。Fω为外载,它与x轴的角度为它到转台中心的距离为b。h为转台所有支承油垫油膜厚度的平均值,它等于转台中心在竖直方向的位移加上油膜的初始厚度h0,θx和θy为在几何轮廓误差作用下转台分别对x轴和y轴的倾角。显然转台的最大倾角为近似认为转台自身不发生变形,所以转台的动力学平衡方程为
其中M是转台质量,J是转台在外载Fω和几何轮廓误差作用下偏转的转动惯量为
由于认为转台为刚体,对于转台的导轨面而言,其表面的波动波长一般大于油垫的直径。所以忽略转台倾角引起的油膜厚度不均匀对油垫承载能力的影响。所以对于每一个油垫而言其油膜厚度及速度的计算为:
其中i=1~8,j=1~16。
步骤(2),计算程序编写
其程序计算流程如图6所示,首先根据(10)、(15)、(16)、(17)式在Matlab软件中编写油膜厚度计算函数,支承油垫承载力计算函数,预压油垫承载力计算函数及根据平衡方程编写转台的各个方向的加速度计算函数,然后在主函数中写入转台各项输入参数包括转台的结构尺寸,各油垫的结构尺寸,油液参数及外载和边界条件,再调用承载力计算函数计算出各油垫的承载力,带入平衡方程,编写平衡方程计算转台各方向的加速度对其进行积分就得到转台的位移与速度,根据龙格库塔法,上述计算需要进行四次修正。之后重复上述过程进行下一时间的计算值到达到时间边界条件计算才终止。
下面将给出一个计算实例,转台的各参数如表一所示。图7分别给出了不同波长下λ=0.78rad,A=0.2h0(a图)λ=1.25rad,A=0.2h0(b图)时运动误差随转台转速的变化规律,纵轴代表转台的运动误差,即波动比值,为油膜厚度的变化量与油膜初始厚度的比值,横轴为转台转速,A为导轨面的波动幅值。当运动误差大于100%时代表油膜对于几何误差有放大作用,当运动误差等于100%时代表转台的运动误差等于转台的几何误差,而当转台的运动误差小于100%时代表转台的油膜对于转台的几何误差有抑制作用,且运动误差越小。可以看出,在转台转速在ω=5~60r/min区间上,随着转速的提高,波动比值越来越小,运动误差越小,当转台转速ω=60r/min时,波动比值是最小的,分别为21.104%,0.078%。并且,当波长λ=1.25rad,A=0.2h0时的波动比值要明显小于波长λ=0.78rad,A=0.2h0时的波动比值。
在内外圈导轨面均存在制造误差的情况下,在Matlab中求解出转台的运动误差,可以得出如图8所示的转台内外圈油垫处的波动比值随波长的变化规律,其中纵轴为波动比值,即油膜厚度的变化量与转台导轨面的波动幅值的比值z,横轴为无量纲化的导轨轮廓波长λ,即导轨轮廓实际波长与油垫直径的比值。虚线表示外圈油垫处的波动比值随波长的变化规律,实线表示内圈油垫处的波动比值随波长的变化规律。在计算中,可以得到对内圈导轨处的油垫来讲,当λ内=3.62时,波动比值最小为82.56%;而对外圈导轨处的油垫来讲,当λ外=3.49时,波动比值最小为22.78%。
通过以上实例分析总结出:本发明方法可以准确的建立转台的动力学模型,快速的计算出转台转速、加工轮廓误差等主要参数对运动误差的影响,为转台的设计提供理论指导,对转台的使用也有一定的借鉴。
表1静压转台结构及供油系统参数
Claims (2)
1.一种静压转台的运动误差的建模及计算方法,其特征在于:该方法包括下述流程,
1)忽略静压转台自身的变形,即近似认为静压转台为刚体,同时考虑静压转台的加工轮廓误差,并将加工轮廓误差拟合成函数,对于静压转台的导轨面而言其表面的波动波长大于静压转台的直径;所以忽略静压转台倾角引起的油膜厚度不均匀对油垫承载能力的影响;
2)基于雷诺方程的单油垫理论模型计算出预压油垫及支承油垫在考虑加工轮廓误差时的压力分布及之后对压力分布进行积分得到预压油垫及支承油垫的承载力Fy,Fi和Fj;
3)静压转台在外载作用时其所有受力及力矩应当平衡,那么根据力平衡及力矩平衡原则可以建立静压转台的动力学平衡方程;
4)考虑静压转台的加工轮廓误差,通过各油垫油膜厚度之间的几何关系确定各油膜厚度及速度;
5)静压转台的动力学方程为二阶非线性微分方程,根据数值分析的相关理论可以用龙格库塔法在matlab软件中编写方程的求解程序最终得到方程的数值解,也就是求出静压转台转速、加工轮廓误差对静压转台运动误差的影响。
2.根据权利要求1所述的一种静压转台的运动误差的建模及计算方法,其特征在于:静压转台由转台,基座,支承油垫,预压油垫、径向轴承组成,其中支承油垫为双圈圆形油垫,预压油垫为环形预压油垫,在转台的中心部位安装有一个径向轴承;转台各油垫都由定量泵供油,由于加工的原因转台的导轨面会高低不平,用正弦函数来等效描述转台导轨面的起伏;转台的自重为G,其中支承油垫的供油量为Q0,预压油垫的供油量为Q1,径向油垫的供油量为Q2;该方法包括下述流程,
步骤(1),相关公式推导
基于雷诺方程的单油垫理论模型,推导出转台支承油垫及预压油垫承载力公式,而后列出转台平衡方程;
1、建立单油垫理论模型
在单油垫内,任意割取半径为dr,高度为dz的微球体,定义该微球体中心点坐标为(r,z),在某一瞬时t时,过此微球体中心点的两个分速分别为ur、vz;p为油膜压力,η为粘度系数;根据量纲分析可以将极坐标系下的连续性方程和N-S方程简化为:
<mrow>
<mfrac>
<mn>1</mn>
<mi>r</mi>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>ru</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>z</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>p</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>&eta;</mi>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>u</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>p</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
由(3)式可认为p与z无关,所以对(2)式两边对z进行二次积分并代入边界条件:z=h,ur=0;z=0,ur=0可得
<mrow>
<msub>
<mi>u</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>p</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
将(4)式带入(1)式并对z积分,再代入边界条件z=h,ur=0,z=0,ur=0,vz=0可以得到雷诺方程为:
<mrow>
<mfrac>
<mn>1</mn>
<mi>r</mi>
</mfrac>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mi>rh</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mn>12</mn>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>p</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>rh</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mn>12</mn>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>p</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
分析对象中,支承油垫为圆形油垫,预紧油垫为环形油垫;
对于圆形油垫有:
r=R1,p=p0;r=R2,p=0;Q(R1)=Q0 (7)
根据(5)、(6)、(7)式可解得油腔压力p0和封油边压力分布p(r):
<mrow>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>6</mn>
<mi>&eta;</mi>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&pi;h</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>Q</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msubsup>
<mi>&pi;R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mn>3</mn>
<mi>&eta;</mi>
</mrow>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<msubsup>
<mi>h</mi>
<mi>T</mi>
<mn>3</mn>
</msubsup>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
那么圆形油垫的承载能力为:
<mrow>
<mi>F</mi>
<mo>=</mo>
<msubsup>
<mi>&pi;R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>2</mn>
<mi>&pi;</mi>
<munderover>
<mo>&Integral;</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</munderover>
<mi>r</mi>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
即:
<mrow>
<mi>F</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>(</mo>
<mn>24</mn>
<msub>
<mi>Q</mi>
<mn>0</mn>
</msub>
<mi>&eta;</mi>
<mo>-</mo>
<mn>12</mn>
<mi>&pi;</mi>
<mi>&eta;</mi>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>8</mn>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
对于环形预压压油垫:
对于环形油垫有边界条件为:
r=RC1,p=0;r=RC2,p=p0;r=RC3,p=p0;r=RC4,p=0;Q1=-Q(RC2)+Q(RC3) (11)
根据(5)、(6)、(11)式可解得油腔压力p0和封油边压力分布p(r):
<mrow>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>6</mn>
<mi>&eta;</mi>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&pi;h</mi>
<mn>3</mn>
</msup>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Q</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
当r∈(RC1,RC2)时压力分布为:
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mn>3</mn>
<mi>&eta;</mi>
</mrow>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
当r∈(RC3,RC4)时压力分布为:
<mrow>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mn>3</mn>
<mi>&eta;</mi>
</mrow>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
那么根据(12)(13)(14)式可以计算环形油垫的承载力为:
<mrow>
<mi>F</mi>
<mo>=</mo>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>2</mn>
<mi>&pi;</mi>
<munderover>
<mo>&Integral;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>c</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>c</mi>
<mn>2</mn>
</mrow>
</msub>
</munderover>
<msub>
<mi>rp</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
<mi>&pi;</mi>
<munderover>
<mo>&Integral;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>c</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>c</mi>
<mn>3</mn>
</mrow>
</msub>
</munderover>
<msub>
<mi>rp</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
即:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mn>3</mn>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>&pi;</mi>
<mi>&eta;</mi>
</mrow>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>C</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
2、转台力平衡方程
转台由两圈油垫共24个油垫支承,对于第一圈油垫,其支承力分别为Fi,i=1~8,它们与加力点和转台中心的连线的夹角分别为i=1~8,油垫中心到转台中心的距离为RL,对于第二圈支承油垫,其支承力分别为Fj,j=1~16,它们与加力点和转台中心的连线的夹角分别为j=1~16,油垫中心到转台中心的距离为Rs,预压油垫的支承力为Fy;转台的半径为RT,转台的转速为ω;Fω为外载,它与x轴的角度为它到转台中心的距离为b;h为转台所有支承油垫油膜厚度的平均值,它等于转台中心在竖直方向的位移加上油膜的初始厚度h0,θx和θy为在几何轮廓误差作用下转台分别对x轴和y轴的倾角;显然转台的最大倾角为近似认为转台自身不发生变形,所以转台的动力学平衡方程为
其中M是转台质量,J是转台在外载Fω和几何轮廓误差作用下偏转的转动惯量为
由于认为转台为刚体,对于转台的导轨面而言,其表面的波动波长一般大于油垫的直径;所以忽略转台倾角引起的油膜厚度不均匀对油垫承载能力的影响;所以对于每一个油垫而言其油膜厚度及速度的计算为:
其中i=1~8,j=1~16;
步骤(2),计算程序编写
首先根据(10)、(15)、(16)、(17)式在Matlab软件中编写油膜厚度计算函数,支承油垫承载力计算函数,预压油垫承载力计算函数及根据平衡方程编写转台的各个方向的加速度计算函数,然后在主函数中写入转台各项输入参数包括转台的结构尺寸,各油垫的结构尺寸,油液参数及外载和边界条件,再调用承载力计算函数计算出各油垫的承载力,带入平衡方程,编写平衡方程计算转台各方向的加速度对其进行积分就得到转台的位移与速度,根据龙格库塔法,上述计算需要进行四次修正;之后重复上述步骤(2)过程进行下一时间的计算值到达到时间边界条件计算才终止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510028962.2A CN104537189B (zh) | 2015-01-21 | 2015-01-21 | 一种静压转台运动误差建模及计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510028962.2A CN104537189B (zh) | 2015-01-21 | 2015-01-21 | 一种静压转台运动误差建模及计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104537189A CN104537189A (zh) | 2015-04-22 |
CN104537189B true CN104537189B (zh) | 2018-02-02 |
Family
ID=52852716
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510028962.2A Active CN104537189B (zh) | 2015-01-21 | 2015-01-21 | 一种静压转台运动误差建模及计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104537189B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105677972A (zh) * | 2016-01-07 | 2016-06-15 | 燕山大学 | 一种液体静压导轨的精确承载力计算方法 |
CN106354987B (zh) * | 2016-11-09 | 2019-09-27 | 北京工业大学 | 一种重载静压转台承载力与油垫温度场分布规律关系计算方法 |
CN107061494A (zh) * | 2016-11-09 | 2017-08-18 | 北京工业大学 | 一种用于非对称工件装夹的静压转台支撑垫铁布局方法 |
CN106767482A (zh) * | 2017-01-08 | 2017-05-31 | 北京工业大学 | 一种小型精密静压实验台台面变形检测装置 |
CN109815548B (zh) * | 2018-12-27 | 2021-01-19 | 西安交通大学 | 一种基于Garlerkin思想的流体膜压力计算方法 |
CN113094832A (zh) * | 2021-04-08 | 2021-07-09 | 西安交通大学苏州研究院 | 径向静压轴承径向回转精度评价方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246786A (zh) * | 2013-05-24 | 2013-08-14 | 哈尔滨理工大学 | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7219747B2 (en) * | 2004-03-04 | 2007-05-22 | Halliburton Energy Services, Inc. | Providing a local response to a local condition in an oil well |
US8100196B2 (en) * | 2005-06-07 | 2012-01-24 | Baker Hughes Incorporated | Method and apparatus for collecting drill bit performance data |
-
2015
- 2015-01-21 CN CN201510028962.2A patent/CN104537189B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246786A (zh) * | 2013-05-24 | 2013-08-14 | 哈尔滨理工大学 | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 |
Non-Patent Citations (3)
Title |
---|
仿真转台误差分析及误差建模;曲智勇 等;《计算机仿真》;20060331;第23卷(第3期);全文 * |
液体静压转台倾斜油膜承载特性分析;熊万里 等;《中国机械工程》;20141231;第25卷(第24期);全文 * |
重型数控机床静压转台承载力分析;赵永胜 等;《高技术通讯》;20130531;第23卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104537189A (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104537189B (zh) | 一种静压转台运动误差建模及计算方法 | |
CN104776998B (zh) | 一种基于动态刚度系数和阻尼系数的转子轴心轨迹求解方法 | |
CN104143026B (zh) | 一种考虑倾斜与热的静压油垫性能计算方法 | |
CN103823940B (zh) | 一种基于整体动力学模型的静压转台动态响应计算方法 | |
CN103226634B (zh) | 基于三维动网格的旋喷泵非定常流场的计算方法 | |
CN105653781B (zh) | 一种复合材料螺旋桨空泡性能的计算方法 | |
CN103632006B (zh) | 一种多级离心泵湿转子临界转速的确定方法 | |
CN107169158B (zh) | 一种基于流固耦合作用的静压滑座工作性能计算方法 | |
Zheng et al. | Generation of noncircular spiral bevel gears by face-milling method | |
CN105608297A (zh) | 一种角接触球轴承动刚度的数值模拟方法 | |
CN102902859B (zh) | 基于相变模型的动静压滑动轴承流固耦合设计方法 | |
CN110096784B (zh) | 一种具有轴向压差的径向滑动轴承的快速计算与设计方法 | |
CN104091021B (zh) | 一种计算倾斜时油垫承载力的方法 | |
CN105005705B (zh) | 一种改进的圆形静压油垫离心力计算方法 | |
CN105181327B (zh) | 一种裂纹轮齿啮合刚度计算的方法 | |
CN107526914B (zh) | 基于结构化动网格的可倾瓦滑动轴承变流域流场计算方法 | |
CN104036079B (zh) | 一种基于灵敏度分析的定量式静压转台优化设计方法 | |
CN106777467A (zh) | 一种考虑流固耦合的静压滑座静平衡计算方法 | |
CN107885699A (zh) | 一种建立立式水泵机组轴系轨迹的轴系振动摆度解析表达式的方法 | |
CN104615801A (zh) | 重型静压轴承临界润滑状态转速值确定方法 | |
CN104091004B (zh) | 一种基于改进的粒子群算法的定量式静压转台优化设计方法 | |
CN109117522A (zh) | 一种基于matlab-ansys软件的静压溜板流固耦合工作性能的计算方法 | |
CN105354364B (zh) | 一种基于云平台的静压支撑系统模型制造方法 | |
CN105095536B (zh) | 一种考虑表面形貌的静压油垫流场特性数值仿真方法 | |
CN102722656B (zh) | 基于混合维模型的精密偶件配合间隙泄漏特性预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |