CN109063404B - 一种空气静压止推轴承刚度的动态预测方法 - Google Patents
一种空气静压止推轴承刚度的动态预测方法 Download PDFInfo
- Publication number
- CN109063404B CN109063404B CN201811231695.9A CN201811231695A CN109063404B CN 109063404 B CN109063404 B CN 109063404B CN 201811231695 A CN201811231695 A CN 201811231695A CN 109063404 B CN109063404 B CN 109063404B
- Authority
- CN
- China
- Prior art keywords
- bearing
- dynamic
- stiffness
- static
- air film
- 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
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
本发明公开了一种空气静压止推轴承刚度的动态预测方法,该方法为一种基于流体计算软件Fluent的动网格建模技术求解轴承静刚度的方法,该方法包括下列步骤:根据空气静压轴承的结构及原理,利用通用CFD流体计算软件ANSYS Fluent的动网格建模技术建立空气静压轴承非稳态计算的数值模型;利用傅里叶变换对CFD计算得到的轴承非稳态时动态承载力的变化数据进行处理,得到轴承气膜的动态刚度曲线;借助数学分析和数值试验,可以建立轴承气膜的动、静刚度在低频段几乎相等的理论判断,进而求得轴承气膜的静刚度。
Description
技术领域
本发明涉及一种动、静刚度的动态预测方法,适用于空气静压止推轴承在不同工况下动、静刚度的求解和预测,实现静压气体轴承动、静刚度的精确预测。
背景技术
超精密制造装备代表着精密制造业的最高发展水平,随着国家精密制造实力的增强,空气静压轴承逐渐取代传统机械接触或液体润滑式轴承,成为制造装备实现高精度、高可靠性的有效手段。空气静压轴承由于其近零摩擦、无发热、有效隔离振动传递等优点,常被空间飞行器三轴气浮仿真装置、惯导测试设备和超精密机床主轴部件选作精密支承,实现设备所需要的回转精度。轴承气膜刚度是静压气体轴承最主要的技术指标之一,因而其动静态刚度的测试成为轴承研制的一个关键技术问题。
目前常用的静压气体轴承气膜静刚度的测试方法,即通过测量工作点处载荷变化量与相应的气膜高度的变化量,再求其比值得到相应的气膜高度下的轴承静刚度。实际测试中由于测量误差的存在,气膜高度变化量Δh不可能取得很小,因此这种测试方法的测试精度并不高,且必须借助于实验装置,有一定的局限性。对于静压气体轴承动态刚度的测试,现有的测试装置比较复杂,且价格昂贵,适用性不强。
机械振动理论表明,线性系统低频时的动刚度数值等于其静刚度。与气动位置伺服系统的非线性类似,气体润滑膜和负载质量组成的等效弹簧阻尼质量系统,也是一个本质非线性的系统,显然,线性系统低频时的动刚度数值等于其静刚度这一结论并不能立即适用于气膜质量这个本质非线性系统。对于气膜质量系统而言,其静态刚度和动态刚度之间是否仍存在某种联系?搞清楚本质非线性气膜质量系统的这些基本物理特性,对于静压气体轴承的研制和测试,具有重要的现实意义。
发明内容
针对空气静压轴承的动静态刚度的求解问题,本发明提供一种空气静压轴承刚度的动态预测方法,本发明以空气静压止推轴承为例,利用有限元软件建立空气静压轴承非稳态计算的有限元模型;借助数学分析和数值实验,求得轴承动静态刚度数值,并从理论上揭示了轴承气膜静刚度和动刚度之间的联系,实现轴承动、静刚度的求解。
本发明采用的技术方案为一种空气静压止推轴承刚度的动态预测方法,该方法包括以下步骤:
(1)建立空气静压轴承非稳态计算的有限元模型;
采用小参数摄动法建立空气静压轴承非稳态气膜润滑力学方程,对其动态特性进行数学描述,如下式
式中,x和y为力学方程对应的横坐标和纵坐标分量,h0为轴承静平衡位置气膜厚度,p0为轴承静平衡位置气膜内压力,Δh为轴承受到扰动时支承绕静态平衡点的简谐小扰动,Δp为Δh引起的气膜压力变化量,μ为气体动力粘度,q为节流孔处流量,δ为Kronecker数,R为气体常数,T为温度,p为轴承气膜内压力,j为虚数单位,ω为扰动频率。
基于上述力学方程,利用通用CFD流体计算软件ANSYS Fluent的动网格建模方法建立空气静压轴承非稳态计算的数值模型,建立空气静压轴承气膜的三维网格模型如图1所示。
(2)轴承动态承载力的数据处理方法,如下式
式中,F(ω)为动态承载力的傅里叶变换结果,S(ω)为扰动频率的傅里叶变换结果,Ki(w)为复刚度,复刚度的实部即为动态刚度。对建立的计算空气静压轴承动态特性的动网格CFD模型进行仿真分析,得到动态承载力的变化数据如图2所示。对动态承载力数据利用傅里叶变换进行处理后,得到该扰动频率下的动刚度数值。
(3)轴承动刚度变化曲线;
逐一改变扰动频率来进行仿真,计算各扰动频率处的动刚度值,得到轴承系统在整个测试频率范围内的动刚度曲线,如图3所示。
(4)轴承动刚度与静刚度之间联系
气膜静态承载力、静刚度的数学表达式为
式中,F0为气膜静态承载力,pa为大气压力,h为气膜厚度,k为气膜静刚度。
上述给出了轴承动静刚度的数学定义,但从中并不能直接观察出二者之间有什么样的联系,进一步的分析需要借助于数值实验。采用质量块对主轴系统进行加载,利用电感测微仪测量不同相应载荷下空气静压轴承的气膜间隙的变化。根据空气静压轴承加载载荷与电感测微仪检测的相应的位移变化量,采用刚度的计算公式,得出实验测量的空气静压轴承的静承载力及静刚度曲线,图4a为静承载力随气膜高度的变化曲线,图4b为对应的静刚度变化曲线。比较试验所得静刚度数值与仿真求得的动刚度变化曲线,发现静刚度数值与低频段动刚度数值几乎相等,经过大量的数值计算表明,任选其他轴承,都可以得到同样的结论。静压气体轴承气膜的这个特征为求解轴承的动静态刚度提供了一种新的、准确的方法。
与现有技术相比,本发明具有以下优点:
本发明以空气静压止推轴承为例,基于流体计算动力学软件Fluent的动网格建模技术建立轴承动静态刚度的有限元模型;利用傅里叶变换对仿真得到的动态承载力进行数学处理,可以得到静压轴承动刚度随扰动频率的变化曲线;借助数学分析和数值试验,从理论上揭示了轴承气膜静刚度和动刚度之间的联系,指出静压气体轴承静刚度的求解可以通过低频下气膜动刚度来代替。实测结果表明,静压气体轴承刚度的动态求解方法是正确可行的,气膜与承载质量构成的本质非线性二阶振动系统,其低频段的动态刚度数值仍然等于其静刚度。
目前常用加载实验来测量静压气体轴承的气膜静刚度,但由于测量误差的存在,气膜高度的变化量不可能取得很小,因此这种测试方法的测试精度并不高。本发明基于仿真分析和数值实验,可以提高静刚度的求解精度和效率。
本发明可以采用较为简单的方法得到轴承的动静态刚度,这对于气体轴承稳定性及精密轴系动态性能的研究具有重要的意义,是轴承的加载实验所无法实现的。
附图说明
图1为空气静压轴承气膜的三维网格模型。
图2为动态承载力的变化数据。
图3为动刚度曲线
图4a为静承载力随气膜高度的变化曲线,
图4b为对应的静刚度变化曲线。
具体实施方式
本发明附所述方法均由流体计算软件Fluent以及MATLAB软件编程程序实现。
本发明所述方法具体包括以下步骤:
步骤1,利用通用CFD流体计算软件ANSYSFluent的动网格建模技术建立空气静压轴承非稳态计算的仿真模型,定义边界为正弦扰动并由UDF定义,图1为其三维网格模型。
步骤2,利用傅里叶变换对仿真求得的动态承载力的变化数据进行处理,从图2可以看出动态承载力的变化数据近似正弦变化,因此可采用正弦函数进行曲线拟合,进而得到动态承载力的函数,从而可计算得到复刚度,复刚度的实部即为动态刚度。
步骤3,逐一改变扰动频率来进行仿真,计算各扰动频率处的动刚度值,可以绘制出轴承系统在整个测试频率范围内的动刚度曲线,图3表明,在小于50Hz的低频段,气膜的动态刚度是一条水平直线,此时动态刚度的数值与频率无关。超过50Hz后,气膜的动刚度曲线发生转折,随频率的增加动刚度数值快速增大。
步骤4,利用加载实验,求得静态承载力随气膜高度的变化曲线,再由静承载力和气膜高度之间的增量比来获得轴承的静刚度。由图4a可知空气静压轴承的承载力随气膜厚度增加而减小。由图4b可以看出气膜厚度为15μm的静刚度约为60N/μm,显然与图3中低频段的动刚度数值非常接近。这表明轴承的静刚度与低频段动刚度数值几乎相等,这个特性不仅在上述气膜高度的轴承处成立,大量的数值计算表明,任选其他气膜厚度的轴承,都可以得到同样的结论。借助于数值试验,可以建立轴承气膜的动静态刚度在低频段几乎相等的理论判断。因此可以利用流体计算软件Fluent的动网格建模技术求得轴承的动静态刚度。
Claims (2)
1.一种空气静压止推轴承刚度的动态预测方法,其特征在于:该方法包括以下步骤:
(1)建立空气静压轴承非稳态计算的有限元模型;
采用小参数摄动法建立空气静压轴承非稳态气膜润滑力学方程,对其动态特性进行数学描述,如下式
式中,x和y为力学方程对应的横坐标和纵坐标分量,h0为轴承静平衡位置气膜厚度,p0为轴承静平衡位置气膜内压力,Δh为轴承受到扰动时支承绕静态平衡点的简谐小扰动,Δp为Δh引起的气膜压力变化量,μ为气体动力粘度,q为节流孔处流量,δ为Kronecker数,R为气体常数,T为温度,p为轴承气膜内压力,j为虚数单位,ω为扰动频率;
基于力学方程,利用通用CFD流体计算软件ANSYS Fluent的动网格建模方法建立空气静压轴承非稳态计算的数值模型,建立空气静压轴承气膜的三维网格模型;
(2)轴承动态承载力的数据处理方法,如下式
式中,F(ω)为动态承载力的傅里叶变换结果,S(ω)为扰动频率的傅里叶变换结果,Ki(w)为复刚度,复刚度的实部即为动态刚度;对建立的计算空气静压轴承动态特性的动网格CFD模型进行仿真分析,得到动态承载力的变化数据;对动态承载力数据利用傅里叶变换进行处理后,得到该扰动频率下的动刚度数值;
(3)轴承动刚度变化曲线;
逐一改变扰动频率来进行仿真,计算各扰动频率处的动刚度值,得到轴承系统在整个测试频率范围内的动刚度曲线;
(4)轴承动刚度与静刚度之间联系
气膜静态承载力、静刚度的数学表达式为
式中,F0为气膜静态承载力,pa为大气压力,h为气膜厚度,k为气膜静刚度;
采用质量块对主轴系统进行加载,利用电感测微仪测量不同相应载荷下空气静压轴承的气膜间隙的变化;根据空气静压轴承加载载荷与电感测微仪检测的相应的位移变化量,采用刚度的计算公式,得出实验测量的空气静压轴承的静承载力及静刚度曲线。
2.一种空气静压止推轴承刚度的动态预测方法,其特征在于:比较试验所得静刚度数值与仿真求得的动刚度变化曲线,静刚度数值与低频段动刚度数值几乎相等。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811231695.9A CN109063404B (zh) | 2018-10-22 | 2018-10-22 | 一种空气静压止推轴承刚度的动态预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811231695.9A CN109063404B (zh) | 2018-10-22 | 2018-10-22 | 一种空气静压止推轴承刚度的动态预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109063404A CN109063404A (zh) | 2018-12-21 |
CN109063404B true CN109063404B (zh) | 2023-04-07 |
Family
ID=64765177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811231695.9A Active CN109063404B (zh) | 2018-10-22 | 2018-10-22 | 一种空气静压止推轴承刚度的动态预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109063404B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110276136B (zh) * | 2019-06-25 | 2023-06-06 | 江苏集萃精凯高端装备技术有限公司 | 一种止推板、止推板设计方法及空气静压轴承 |
CN111144034B (zh) * | 2020-01-17 | 2022-06-24 | 西南科技大学 | 基于Fluent被动网格计算闭式气浮轴承动态的方法 |
CN111666644A (zh) * | 2020-06-18 | 2020-09-15 | 哈尔滨工业大学 | 一种参数化cfd气体静压轴承性能分析方法 |
CN111751220B (zh) * | 2020-06-22 | 2023-06-02 | 哈尔滨工业大学 | 一种考虑流固耦合的气体静压轴承性能计算方法 |
CN112231852B (zh) * | 2020-09-09 | 2023-12-22 | 东南大学 | 静压圆柱导轨五自由度静刚度建模与计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880765A (zh) * | 2012-10-18 | 2013-01-16 | 哈尔滨工业大学 | 用于气体静压主轴动力学预测的仿真方法 |
CN105608297A (zh) * | 2016-02-17 | 2016-05-25 | 中国科学院沈阳自动化研究所 | 一种角接触球轴承动刚度的数值模拟方法 |
CN106644475A (zh) * | 2016-10-12 | 2017-05-10 | 河南科技大学 | 一种微型气体轴承静动态特性及稳定性测试方法 |
CN108287947A (zh) * | 2017-12-30 | 2018-07-17 | 北京工业大学 | 一种空气静压主轴径向回转误差预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI589791B (zh) * | 2015-10-28 | 2017-07-01 | 財團法人工業技術研究院 | 液靜壓軸承監測系統及其方法 |
-
2018
- 2018-10-22 CN CN201811231695.9A patent/CN109063404B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880765A (zh) * | 2012-10-18 | 2013-01-16 | 哈尔滨工业大学 | 用于气体静压主轴动力学预测的仿真方法 |
CN105608297A (zh) * | 2016-02-17 | 2016-05-25 | 中国科学院沈阳自动化研究所 | 一种角接触球轴承动刚度的数值模拟方法 |
CN106644475A (zh) * | 2016-10-12 | 2017-05-10 | 河南科技大学 | 一种微型气体轴承静动态特性及稳定性测试方法 |
CN108287947A (zh) * | 2017-12-30 | 2018-07-17 | 北京工业大学 | 一种空气静压主轴径向回转误差预测方法 |
Non-Patent Citations (4)
Title |
---|
A prediction and evaluation system of the impact factors on the performance of the aerostatic slider;CHEN Dongju等;《Robotics and Computer-Integrated Manufacturing》;20180430;第213-221页 * |
Investigation the gas film in micro scale induced error on the performance of the aerostatic spindle in ultra-precision machining;CHEN Dongju等;《Mechanical Systems and Signal Processing》;20180531;第488-501页 * |
静压止推气浮轴承动特性的微扰动法研究;卢志伟等;《兵工学报》;20160115(第01期);第149-154页 * |
静压气体轴承静刚度的动态测试新方法;郭良斌 等;《机械工程学报》;20070430;第25-30页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109063404A (zh) | 2018-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109063404B (zh) | 一种空气静压止推轴承刚度的动态预测方法 | |
Zhang et al. | Elastic ring deformation and pedestal contact status analysis of elastic ring squeeze film damper | |
Li et al. | Comparisons of rotordynamic characteristics predictions for annular gas seals using the transient computational fluid dynamic method based on different single-frequency and multifrequency rotor whirling models | |
Machado et al. | Discussion about nonlinear boundaries for hydrodynamic forces in journal bearing | |
Meybodi et al. | Numerical analysis of a rigid rotor supported by aerodynamic four-lobe journal bearing system with mass unbalance | |
Hahn et al. | Harmonic balance analysis of general squeeze film damped multidegree-of-freedom rotor bearing systems | |
Merelli et al. | Dynamic coefficients of finite length journal bearing. Evaluation using a regular perturbation method | |
Patel et al. | Nonlinear dynamic response of cylindrical roller bearing–rotor system with 9 degree of freedom model having a combined localized defect at inner–outer races of bearing | |
Zuo et al. | Comparative performance analysis of conical hydrostatic bearings compensated by variable slot and fixed slot | |
Chasalevris | Analytical evaluation of the static and dynamic characteristics of three-lobe journal bearings with finite length | |
Yoon et al. | Introduction to rotor dynamics | |
Delgado et al. | Dynamic force coefficients of hydrostatic gas films for recessed flat plates: experimental identification and numerical predictions | |
CN106354974A (zh) | 一种滚动轴承等效刚度与等效阻尼的计算方法 | |
San Andres et al. | Turbocharger nonlinear response with engine-induced excitations: predictions and test data | |
Andrés et al. | On the experimental dynamic force performance of a squeeze film damper supplied through a check valve and sealed with O-Rings | |
Zhuang et al. | Effect of surface waviness on the performances of an aerostatic thrust bearing with orifice-type restrictor | |
Yu et al. | Thermal-fluid-solid coupling deformation of hydrostatic thrust bearing friction pairs | |
Du et al. | Nonlinear dynamic analysis of a rigid rotor supported by a spiral-grooved opposed-hemisphere gas bearing | |
Bensouilah et al. | Elasto‐aerodynamic lubrication analysis of a self‐acting air foil journal bearing | |
CN110276136B (zh) | 一种止推板、止推板设计方法及空气静压轴承 | |
Du et al. | Dynamic modeling of hydrostatic guideway considering compressibility and inertia effect | |
Andrés et al. | On the force coefficients of a flooded, open ends short length squeeze film damper: from theory to practice (and back) | |
Han et al. | Dynamical characteristic analysis of elastic ring squeeze film damper in rotor system | |
Bert et al. | Transient analysis of gas-lubricated journal bearing systems by differential quadrature | |
Li et al. | Investigation on the stability of periodic motions of a flexible rotor-bearing system with two unbalanced disks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |