CN106503391B - 一种考虑流固耦合的矩形静压油垫快捷计算方法 - Google Patents

一种考虑流固耦合的矩形静压油垫快捷计算方法 Download PDF

Info

Publication number
CN106503391B
CN106503391B CN201610987511.6A CN201610987511A CN106503391B CN 106503391 B CN106503391 B CN 106503391B CN 201610987511 A CN201610987511 A CN 201610987511A CN 106503391 B CN106503391 B CN 106503391B
Authority
CN
China
Prior art keywords
guide pass
deformation
oil film
node
film thickness
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
Application number
CN201610987511.6A
Other languages
English (en)
Other versions
CN106503391A (zh
Inventor
赵永胜
张海明
刘志峰
蔡力钢
程强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201610987511.6A priority Critical patent/CN106503391B/zh
Publication of CN106503391A publication Critical patent/CN106503391A/zh
Application granted granted Critical
Publication of CN106503391B publication Critical patent/CN106503391B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Bearings And Hydrostatic Bearings (AREA)

Abstract

本发明公开了一种考虑流固耦合的矩形静压油垫快捷计算方法，更具体是一种考虑了静压油垫导轨面变形和油垫压力分布相互影响的矩形承载力快捷计算方法。该方法首先进行网格划分，计算并储存导轨面标准变形，再依据雷诺方程和变形方程相互迭代求解油膜压力分布，最后积分求解油膜承载能力。其特点是在计算静压油膜承载力过程中，充分考虑了流固耦合作用对承载力的影响。并且提出了一种导轨面变形的快捷计算方法，在计算导轨面变形时不必每次计算导轨面变形都计算一次二重积分，极大的提高设计效率。

Description

一种考虑流固耦合的矩形静压油垫快捷计算方法
技术领域
本发明属于静压支承进给系统设计分析领域，涉及一种考虑了流固耦合的矩形油垫承载力计算方法，更具体是一种考虑了静压油垫导轨面变形和油垫压力分布相互影响的矩形承载力快捷计算方法。
背景技术
静压支承进给系统是重型高档数控机床的关键部件，其承载原理是用有压力的流体使有相对运动的两个表面分开并借助流体静压来承载。由于运动副之间完全被油膜隔开，所以运动副间的摩擦力大大减小，同时其承载能力、运动精度与寿命却大大提高。由于其高负载，低摩擦，高精度等众多油垫，静压支承进给系统被广泛运用于重型高档数控机床中存在相对运动的部件之间，如：滑座与导轨间、立柱与横梁间、横梁与溜板间、溜板与滑枕间。静压支承进给系统中油膜厚度通常在0.1毫米左右，由于其高承载，不可避免的会引起与油膜接触的导轨面发生变形，导轨面发生变形后，油膜厚度同时发生变化，油膜厚度变化后必然影响油膜压力分布，进而油膜对导轨面的作用力发生变化，导轨面变形再次发生变化，如此循环。通过静压支承基本原理可知，油膜厚度的变化与油垫承载力的变化之间是三次方关系，油膜厚度的一点变化会极大的影响静压支撑系统的承载能力。因此在静压支承系统设计分析过程中必须考虑流固耦合即油膜压力分布于导轨面变形的相互影响。目前计算导轨面变形过程时，通常利用弹性力学基本原理，每次计算导轨面变形时需要计算一次二重积分，由于重型高档数控机床中油垫尺寸巨大，计算节点多，计算速度十分缓慢，极大的影响设计效率。本发明旨在提出一种考虑流固耦合的矩形油垫快捷计算方法，不必每次计算导轨面变形都计算一次二重积分，极大的提高设计效率。
发明内容
本发明旨在提供一种矩形油垫承载力的快捷计算方法。该方法的主要特点是考虑了流固耦合作用，但在迭代过程中计算导轨面变形时不必每次都计算二重积分，极大地减少了计算量。
本发明是采用以下技术手段实现的：
1、根据导轨面基本尺寸对导轨面进行网格划分，确定计算节点数目。
2、计算导轨面标准变形，计算每个节点只在该节点受单位力作用，其余节点所受作用力为零时，导轨面的变形情况。并储存每个节点的变形结果。
3、根据初始条件(油垫基本尺寸及油膜厚度)及雷诺方程，计算不考虑流固耦合时，油垫各计算节点的压力。求解雷诺方程时，计算节点与初始划分的节点一一对应。
4、计算各节点作用力乘以与之对应的导轨面标准变形中导轨面变形矩阵，并将所有矩阵叠加，得到该作用力下导轨面最终变形。
5、将最后一次计算的导轨面最终变形矩阵与前一次导轨面变形矩阵相减，两者之间差值在设定的误差允许范围内则终止计算，若不在最终标准内，则将最后一次导轨面变形矩阵加到初始油膜厚度上，利用新油膜厚度及雷诺方程，重新计算油垫各计算节点压力。之后重复步骤4。
6、将各节点作用力进行积分，得到考虑了流固耦合作用的矩形静压油垫承载力。
本发明的特点在于计算节点变形时，不必每次都计算一次二重积分，将积分计算转化为数和矩阵相乘计算，极大的提高了计算速度减少了计算量。
附图说明
图1静压支承系统结构简图
具体实施方式
步骤(1)划分网格
根据计算精度需要，对静压导轨面进行网格划分，得到i*j个计算节点。其中i为导轨面X方向上划分节点数，j为导轨面Y方向上划分的节点数。
步骤(2)计算标准变形
根据弹性力学基本原理，利用四边兼职矩形薄板的纳维(Navier)解，求解导轨面在单位集中力W作用下的变形表达式。集中力W的作用点坐标为(m,n)，导轨面X方向上边长为a，Y方向上边长为b，对应坐标分别为x、y。则在坐标为(m,n)处受单位力W作用时，导轨面(x,y)处的变形ω计算公式为：
其中j，k为迭代计算次数，次数越高迭代精度越高，在保证精度前提下，尽量减少计算次数，j，k的值选取为50。其中D为薄板弯曲刚度，计算公式为：
其中E为导轨面弹性模量，μ为导轨面材料泊松比，H为导轨面厚度。通过循环计算，计算出各个节点的变形情况，再分别计算出单位力W在各个节点处作用时，导轨面的变形情况，并储存备用。最终得到i*j个变形矩阵，每个矩阵B0(i，j)包含i*j个数据，与节点一一对应。
步骤(3)计算初始压力分布
静压油膜符合薄膜润滑理论，依据雷诺方程基本假设：流体在界面上无滑动，即贴于表面的流体流速与表面速度相同；在沿润滑膜厚度方向不计压力的变化；忽略油膜曲率的影响,并用平移速度代替转动速度；润滑剂是牛顿流体；流动为层流,油膜中不存在涡流和湍流；与粘性力比较,可忽略惯性力的影响；沿润滑膜厚度方向粘度数值不变。雷诺方程为：
对方程进行无量纲化，并通过有限差分法求解雷诺方程，对方程进行无量纲化：
式中：p为压强；p0为油兜内压强；Ux为导轨X方向移动速度；h为油膜厚度；η为油液粘度。为无量纲压力；为无量纲长度；为无量纲宽度；为无量纲厚度为无量纲导轨移动速度；为无量纲油膜厚度。其中油膜厚度h为矩阵h(i，j)，对应表示各节点处的油膜厚度。在已知初始油膜厚度h的情况下，利用有限差分法求解压强分布P(x,y)。在节点区域对压强分布进行积分，求得各节点处压力分布pf1(x,y)。
步骤(4)计算初始变形
将压力分布pf(x,y)每个节点处压力值对应乘以只在该点处受单位力，其余节点受力为0时的变形矩阵，并将所有节点叠加，得到导轨面变形总矩阵A1(i,j)。
步骤(5)循环终止判断
变形总矩阵Az(i,j)与前一次计算的变形矩阵Az-1(i,j)各项对应相减，挑选其中差值最大点，与终止条件δ＝0.0000000001比较，若最大值大于终止条件值，则进行步骤(6)，若小于终止条件值，则跳出循环。
步骤(6)计算新压力分布和新变形
求将原始油膜厚度与导轨面产生变形量相加得到新油膜厚度h1(i,j)，并利用新油膜厚度重复步骤(3)求解新压力分布pf2，利用新压力分布pf2重复步骤(4)求解导轨面新变形A2(i,j)。
h1(i,j)＝h0(i,j)+A(i,j)
步骤(7)重新进行终止判断，重复步骤(5)。最终油膜厚度分布，导轨面变形分布和对最终压力分布pf1(x,y)在整个油垫范围内进行积分可得油膜承载能力值
F＝∫∫pf1(x,y)dxdy

Claims (2)

1.一种考虑流固耦合的矩形静压油垫快捷计算方法，其特征在于：
S1、根据导轨面基本尺寸对导轨面进行网格划分，确定计算节点数目；
S2、计算导轨面标准变形，计算每个节点只在该节点受单位力作用，其余节点所受作用力为零时，导轨面的变形情况；并储存每个节点的变形结果；
S3、根据油垫基本尺寸、油膜厚度初始条件及雷诺方程，计算不考虑流固耦合时，油垫各计算节点的压力；求解雷诺方程时，计算节点与初始划分的节点一一对应；
S4、计算各节点作用力乘以与之对应的导轨面标准变形中导轨面变形矩阵，并将所有矩阵叠加，得到该作用力下导轨面最终变形；
S5、将计算的导轨面最终变形矩阵与前一次导轨面变形矩阵相减，两者之间差值在设定的误差允许范围内则终止计算，若不在最终标准内，则将最后一次导轨面变形矩阵加到初始油膜厚度上，利用新油膜厚度及雷诺方程，重新计算油垫各计算节点压力；之后重复S4；
S6、将各节点作用力进行积分，得到考虑了流固耦合作用的矩形静压油垫承载力。
2.根据权利要求1所述的一种考虑流固耦合的矩形静压油垫快捷计算方法，其特征在于：
步骤(1)划分网格；
根据计算精度需要，对静压导轨面进行网格划分，得到i*j个计算节点；其中i为导轨面X方向上划分节点数，j为导轨面Y方向上划分的节点数；
步骤(2)计算标准变形；
根据弹性力学基本原理，利用四边兼职矩形薄板的纳维解，求解导轨面在单位集中力W作用下的变形表达式；集中力W的作用点坐标为(m,n)，导轨面X方向上边长为a，Y方向上边长为b，对应坐标分别为x、y；则在坐标为(m,n)处受单位力W作用时，导轨面(x,y)处的变形ω计算公式为：
其中j，k为迭代计算次数，次数越高迭代精度越高，在保证精度前提下，尽量减少计算次数，j，k的值选取为50；D为薄板弯曲刚度，计算公式为：
其中E为导轨面弹性模量，μ为导轨面材料泊松比，H为导轨面厚度；通过循环计算，计算出各个节点的变形情况，再分别计算出单位力W在各个节点处作用时，导轨面的变形情况，并储存备用；最终得到i*j个变形矩阵，每个矩阵B0(i，j)包含i*j个数据，与节点一一对应；
步骤(3)计算初始压力分布；
静压油膜符合薄膜润滑理论，依据雷诺方程基本假设：流体在界面上无滑动，即贴于表面的流体流速与表面速度相同；在沿润滑膜厚度方向不计压力的变化；忽略油膜曲率的影响,并用平移速度代替转动速度；润滑剂是牛顿流体；流动为层流,油膜中不存在涡流和湍流；与粘性力比较，忽略惯性力的影响；沿润滑膜厚度方向粘度数值不变；雷诺方程为：
对方程进行无量纲化，并通过有限差分法求解雷诺方程，对方程进行无量纲化：
式中：p为压强；p0为油兜内压强；Ux为导轨X方向移动速度；h为油膜厚度；η为油液粘度；为无量纲压力；为无量纲长度；为无量纲宽度；为无量纲导轨移动速度；为无量纲油膜厚度；其中油膜厚度h为矩阵h(i，j)，对应表示各节点处的油膜厚度；在已知初始油膜厚度h的情况下，利用有限差分法求解压强分布P(x,y)；在节点区域对压强分布进行积分，求得各节点处压力分布pf1(x,y)；
步骤(4)计算初始变形；
将压力分布pf1(x,y)每个节点处压力值对应乘以只在该点处受单位力，其余节点受力为0时的变形矩阵，并将所有节点叠加，得到导轨面变形总矩阵A1(i,j)；
步骤(5)循环终止判断；
变形总矩阵Az(i,j)与计算的变形矩阵Az-1(i,j)各项对应相减，挑选其中差值最大点，与终止条件δ＝0.0000000001比较，若最大值大于终止条件值，则进行步骤(6)，若小于终止条件值，则跳出循环；
步骤(6)计算新压力分布和新变形；
求将原始油膜厚度与导轨面产生变形量相加得到新油膜厚度h1(i,j)，并利用新油膜厚度重复步骤(3)求解新压力分布pf2，利用新压力分布pf2重复步骤(4)求解导轨面新变形A2(i,j)；
h1(i,j)＝h0(i,j)+A(i,j)
步骤(7)重新进行终止判断，重复步骤(5)；最终油膜厚度分布，导轨面变形分布和对最终压力分布pf1(x,y)在整个油垫范围内进行积分可得油膜承载能力值
F＝∫∫pf1(x,y)dxdy。
CN201610987511.6A 2016-11-09 2016-11-09 一种考虑流固耦合的矩形静压油垫快捷计算方法 Active CN106503391B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610987511.6A CN106503391B (zh) 2016-11-09 2016-11-09 一种考虑流固耦合的矩形静压油垫快捷计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610987511.6A CN106503391B (zh) 2016-11-09 2016-11-09 一种考虑流固耦合的矩形静压油垫快捷计算方法

Publications (2)

Publication Number Publication Date
CN106503391A CN106503391A (zh) 2017-03-15
CN106503391B true CN106503391B (zh) 2019-10-29

Family

ID=58323757

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610987511.6A Active CN106503391B (zh) 2016-11-09 2016-11-09 一种考虑流固耦合的矩形静压油垫快捷计算方法

Country Status (1)

Country Link
CN (1) CN106503391B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107169158B (zh) * 2017-04-09 2021-01-01 北京工业大学 一种基于流固耦合作用的静压滑座工作性能计算方法
CN110378018B (zh) * 2019-07-18 2023-12-26 上海理工大学 一种液体动静压球轴承的稳态承载能力的计算方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462752A (zh) * 2014-10-31 2015-03-25 北京工业大学 一种考虑导轨面变形的静压导轨承载能力分析方法
CN105005705A (zh) * 2015-07-31 2015-10-28 北京工业大学 一种改进的圆形静压油垫离心力计算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462752A (zh) * 2014-10-31 2015-03-25 北京工业大学 一种考虑导轨面变形的静压导轨承载能力分析方法
CN105005705A (zh) * 2015-07-31 2015-10-28 北京工业大学 一种改进的圆形静压油垫离心力计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
偏载液体静压转台旋转工况下承载力及倾覆力矩动网格计算方法;王少力 等;《机械工程学报》;20141231;第50卷(第23期);第66-74页 *
定量式静压转台动态特性建模与影响因素分析;刘志峰 等;《机械工程学报》;20151030;第51卷(第19期);第75-83页 *

Also Published As

Publication number Publication date
CN106503391A (zh) 2017-03-15

Similar Documents

Publication Publication Date Title
Kango et al. Thermal analysis of microtextured journal bearing using non-Newtonian rheology of lubricant and JFO boundary conditions
Hashemi et al. Multibody dynamics of pivot slipper pad thrust bearing in axial piston machines incorporating thermal elastohydrodynamics and mixed lubrication model
Lai et al. Numerical simulation of moving contact lines with surfactant by immersed boundary method
CN107169158B (zh) 一种基于流固耦合作用的静压滑座工作性能计算方法
CN106503391B (zh) 一种考虑流固耦合的矩形静压油垫快捷计算方法
Farooq et al. Steady Poiseuille flow and heat transfer of couple stress fluids between two parallel inclined plates with variable viscosity
Sajid et al. Stretching flows with general slip boundary condition
CN106777467B (zh) 一种考虑流固耦合的静压滑座静平衡计算方法
CN104091021B (zh) 一种计算倾斜时油垫承载力的方法
CN110096784B (zh) 一种具有轴向压差的径向滑动轴承的快速计算与设计方法
CN109117522A (zh) 一种基于matlab-ansys软件的静压溜板流固耦合工作性能的计算方法
Behr On the application of slip boundary condition on curved boundaries
Sharma et al. Effect of non-Newtonian behaviour of lubricant and bearing flexibility on the performance of slot-entry journal bearing
CN107526914B (zh) 基于结构化动网格的可倾瓦滑动轴承变流域流场计算方法
Zhang et al. Numerical and experimental investigation on the performances of a liquid metal bearing with spiral groove structures
CN104462752A (zh) 一种考虑导轨面变形的静压导轨承载能力分析方法
Liu et al. Thermal simulation modeling of a hydrostatic machine feed platform
Yangong et al. Application of direct coupling method to design problems of hydrostatic guideways accounting for fluid-structure interactions
CN105095536B (zh) 一种考虑表面形貌的静压油垫流场特性数值仿真方法
Xing et al. Elastohydrodynamic lubrication analysis with eccentric errors for planetary roller screw mechanism
CN109531409B (zh) 一种复合微织构导轨及其制备方法
CN109002569A (zh) 一种建立静压油膜运动层边界条件的方法
Weinstein et al. Scale dependence of contact line computations
CN106020115A (zh) 一种静压转台载荷误差辨识及其对机床加工误差的影响分析方法
Chen et al. Numerical simulation of superoscillations of a Triton-bearing drop in microgravity

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant