CN102042903A - 一种基于有限元模型的旋转设备支承动刚度参数测量方法 - Google Patents
一种基于有限元模型的旋转设备支承动刚度参数测量方法 Download PDFInfo
- Publication number
- CN102042903A CN102042903A CN 201010512019 CN201010512019A CN102042903A CN 102042903 A CN102042903 A CN 102042903A CN 201010512019 CN201010512019 CN 201010512019 CN 201010512019 A CN201010512019 A CN 201010512019A CN 102042903 A CN102042903 A CN 102042903A
- Authority
- CN
- China
- Prior art keywords
- omega
- rotor
- bearing
- formula
- bpx
- 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.)
- Granted
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种基于有限元模型的旋转设备支承动刚度参数测量方法,先在旋转设备主轴上的轴承转子两侧,沿轴线方向选取两点p、q,在p点所在截面分别安装两个电涡流传感器,使两个电涡流传感器与测点p的连线呈90度,在q点的两个电涡流传感器按同样方式配置,四个电涡流传感器的信号输出端分别与工控机相连,再建立转子有限元模型,并计算出不平衡力fu,推导出支承动刚度ZBF,再通过系列计算得到
Description
技术领域
本发明涉及旋转设备支承动刚度参数测试领域,具体涉及一种基于有限元模型的旋转设备支承动刚度参数测量方法。
背景技术
旋转设备是机械工业中广泛应用的设备,旋转设备支承的动态性能的优劣直接影响着机械加工的质量以及生产率,旋转设备的支承动刚度是其动态性能的主要评价指标,它反映了旋转设备抵抗外载荷的能力,对于工件的加工质量以及旋转设备支承寿命都有影响。对于目前的旋转设备设计及制造,普遍的观点是要通过增加刚度来减少震动影响,因此,动刚度参数的辨识十分重要。然而,国内外对于旋转设备动刚度的研究还处于起步阶段,旋转机械动刚度的测量往往依靠经验进行,并没有明确的动刚度准则,存在动刚度的测试精度低,测量成本高的缺点。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种基于有限元模型的旋转设备支承动刚度参数测量方法,提高了动刚度参数的辨识精度,具有测试精度高,操作简便的优点。
为了达到上述目的,本发明采用的技术方案为:
一种基于有限元模型的旋转设备支承动刚度参数测量方法,包括以下步骤:
第一步,在旋转设备主轴5上的轴承转子7两侧,沿轴线方向选取两点p、q,在p点所在截面分别安装第一电涡流传感器1和第二电涡流传感器2,使第一电涡流传感器1、第二电涡流传感器2与测点p的连线呈90度,在q点所在截面分别安装第三电涡流传感器3和第四电涡流传感器4,使第三电涡流传感器3、第四电涡流传感器4与测点q的连线呈90度,四个电涡流传感器的信号输出端分别与安装在机柜6上的工控机8相连,对转子振动情况进行测试,得到p测点两相互垂直位移量分别为rR,px、rR,py,q测点两相互垂直位移量分别为rR,qx、rR,qy,传送至工控机8进行处理,储存;
其中,Z表示刚度,p表示轴承转子上的测点的自由度,b和i分别表示轴承和转子内部的自由度;
展开 可得到:
ZBFrR,b=-fR,b (1)
式(1)即为轴承支承刚度ZBF的表达式,式(1)可表示为:
[KB(ω)+(jω)2MB(ω)+jωCB(ω)]rR,b(ω)=P(ω) (2)
其中:
rR,i为转子内部振动响应,rR,b为转子和轴承间的振动响应,fR,b为转子和轴承间的作用力;
P(ω)=-fR,b (3)
第四步,将旋转设备支承动刚度分解为转子内部、轴承和测点p进行分析,求得转子内部振动响应rR,i,转子和轴承间的振动响应rR,b,转子和轴承间的作用力fR,b:
将旋转设备支承动刚度分解为转子内部、轴承和测点p进行分析,建立刚度方程:
展开式(4)可以得到:
由式(5)可以求出转子内部振动响应rR,i,转子和轴承间的振动响应rR,b,代入式(6)可以求出转子与轴承作用力fR,b;
第五步,通过支承动刚度表达式计算出支承动刚度实部与虚部表达式KMB(ω)与CB(ω);
对式(2)分离实部和虚部有:
由于对于特定旋转设备质量矩阵一定,因此,在转速ω下,[KB(ω)-ω2MB(ω)]可作为整体,以KMB(ω)表示,
其中,
KMB(ω)={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)}T (9)
CB(ω)={cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T (10)
第六步,将第一步所测得的振动量rR,px、rR,py分别依次代入式(5)、式(6)、式(7)、式(8)中,有:
其中:
rR,bpx表示由p点x方向振动响应rR,px计算得到的转子与轴承间x方向振动响应rR,b,rR,bpy表示由p点y方向振动响应rR,py计算得到的转子与轴承间y方向振动响应rR,b;
βB(ω)=[KMB(ω)CB(ω)]T
={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T
(14)
第七步,将第一步所测得的振动量rR,qx、rR,qy,重复第六步操作,有:
其中:
βB(ω)=[KMB(ω)CB(ω)]T
={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T
(19)
第九步,通过计算出的不同工作转速下的旋转设备支承动刚度参数,可以建立标准,对旋转设备的运行进行安全评定。
由于本方法采用了精确的有限元建模方法,同时采用了合理的传感器进行振动测试,使得本方法测试以及计算得到的数据准确性高,提高了动刚度参数的辨识精度,故而具有测试精度高,操作简便的优点。
附图说明
附图为本发明四个电涡流传感器的安装示意图。
具体实施方法
下面结合附图对本发明作进一步详细描述。
一种基于有限元模型的旋转设备支承动刚度参数测量方法,包括以下步骤:
第一步,参照附图,在旋转设备主轴5上的轴承转子7两侧,沿轴线方向选取两点p、q,在p点所在截面分别安装第一电涡流传感器1和第二电涡流传感器2,使第一电涡流传感器1、第二电涡流传感器2与测点p的连线呈90度,即两传感器测得位移信号相位相差90度,在q点所在截面分别安装第三电涡流传感器3和第四电涡流传感器4,使第三电涡流传感器3、第四电涡流传感器4与测点q的连线呈90度,即两传感器测得位移信号相位相差90度,四个电涡流传感器的信号输出端分别与安装在机柜6上的工控机8相连,对转子振动情况进行测试,得到p测点两相互垂直位移量分别为rR,px、rR,py,q测点两相互垂直位移量分别为rR,qx、rR,qy,传送至工控机8进行处理,储存;
其中,p表示轴承转子上的测点的自由度,b和i分别表示轴承和转子内部的自由度;
展开 可得到:
ZBFrR,b=-fR,b (1)
式(1)即为轴承支承刚度ZBF的表达式。式(1)可表示为:
[KB(ω)+(jω)2MB(ω)+jωCB(ω)]rR,b(ω)=P(ω) (2)
其中:
rR,i为转子内部振动响应,rR,b为转子和轴承间的振动响应,fR,b为转子和轴承间的作用力;
P(ω)=-fR,b (3)
第四步,将旋转设备支承动刚度分解为转子内部、轴承和测点p进行分析,求得转子内部振动响应rR,i,转子和轴承间的振动响应rR,b,转子和轴承间的作用力fR,b:
将旋转设备支承动刚度分解为转子内部、轴承和测点p进行分析,建立刚度方程:
展开式(4)可以得到:
由式(5)可以求出转子内部振动响应rR,i,转子和轴承间的振动响应rR,b,代入式(6)可以求出转子与轴承作用力fR,b;
第五步,通过支承动刚度表达式计算出支承动刚度实部与虚部表达式KMB(ω)与CB(ω);
对式(2)分离实部和虚部有:
由于对于特定旋转设备质量矩阵一定,因此,在转速ω下,[KB(ω)-ω2MB(ω)]可作为整体,以KMB(ω)表示。
其中,
KMB(ω)={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)}T (9)
CB(ω)={cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T (10)
第六步,将第一步所测得的振动量rR,px、rR,py分别依次代入式(5)、式(6)、式(7)、式(8)中,有:
其中:
rR,bpx表不由p点x方向振动响应rR,px计算得到的转子与轴承间x方向振动响应rR,b,rR,bpy表示由p点y方向振动响应rR,py计算得到的转子与轴承间y方向振动响应rR,b;
βB(ω)=[KMB(ω)CB(ω)]T
={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T
(14)
第七步,将第一步所测得的振动量rR,qx、rR,qy,重复第六步操作,有:
其中:
βB(ω)=[KMB(ω)CB(ω)]T
={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T
(19)
第九步,通过计算出的不同工作转速下的旋转设备支承动刚度参数,可以建立标准,对旋转设备的运行进行安全评定。
附图中:1为第一电涡流传感器,2为第二电涡流传感器,3为第三电涡流传感器,4为第四电涡流传感器,5为旋转设备主轴,6为机柜,7为轴承转子,8为工控机。
Claims (1)
1.一种基于有限元模型的旋转设备支承动刚度参数测量方法,其特征在于:包括以下步骤:
第一步,在旋转设备主轴(5)上的轴承转子(7)两侧,沿轴线方向选取两点p、q,在p点所在截面分别安装第一电涡流传感器(1)和第二电涡流传感器(2),使第一电涡流传感器(1)、第二电涡流传感器(2)与测点p的连线呈90度,在q点所在截面分别安装第三电涡流传感器(3)和第四电涡流传感器(4),使第三电涡流传感器(3)、第四电涡流传感器(4)与测点q的连线呈90度,四个电涡流传感器的信号输出端分别与安装在机柜(6)上的工控机(8)相连,对转子振动情况进行测试,得到p测点两相互垂直位移量分别为rR,px、rR,py,q测点两相互垂直位移量分别为rR,qx、rR,qy,传送至工控机(8)进行处理,储存;
其中,Z表示刚度,p表示轴承转子上的测点的自由度,b和i分别表示轴承和转子内部的自由度;
展开 可得到:
ZBFrR,b=-fR,b (1)
式(1)即为轴承支承刚度ZBF的表达式,式(1)可表示为:
[KB(ω)+(jω)2MB(ω)+jωCB(ω)]rR,b(ω)=P(ω) (2)
其中:
rR,i为转子内部振动响应,rR,b为转子和轴承间的振动响应,fR,b为转子和轴承间的作用力;
P(ω)=-fR,b (3)
第四步,将旋转设备支承动刚度分解为转子内部、轴承和测点p进行分析,求得转子内部振动响应rR,i,转子和轴承间的振动响应rR,b,转子和轴承间的作用力fR,b:
将旋转设备支承动刚度分解为转子内部、轴承和测点p进行分析,建立刚度方程:
展开式(4)可以得到:
由式(5)可以求出转子内部振动响应rR,i,转子和轴承间的振动响应rR,b,代入式(6)可以求出转子与轴承作用力fR,b;
第五步,通过支承动刚度表达式计算出支承动刚度实部与虚部表达式KMB(ω)与CB(ω);
对式(2)分离实部和虚部有:
由于对于特定旋转设备质量矩阵一定,因此,在转速ω下,[KB(ω)-ω2MB(ω)]可作为整体,以KMB(ω)表示,
其中,
KMB(ω)={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)}T (9)
CB(ω)={cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T (10)
第六步,将第一步所测得的振动量rR,px、rR,py分别依次代入式(5)、式(6)、式(7)、式(8)中,有:
其中:
rR,bpx表示由p点x方向振动响应rR,px计算得到的转子与轴承间x方向振动响应rR,b,rR,bpy表示由p点y方向振动响应rR,py计算得到的转子与轴承间y方向振动响应rR,b;
βB(ω)=[KMB(ω) CB(ω)]T
={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T
(14)
第七步,将第一步所测得的振动量rR,qx、rR,qy,重复第六步操作,有:
其中:
βB(ω)=[KMB(ω)CB(ω)]T
={kmxx(ω)kmxy(ω)kmyx(ω)kmyy(ω)cxx(ω)cxy(ω)cyx(ω)cyy(ω)}T
(19)
第九步,通过计算出的不同工作转速下的旋转设备支承动刚度参数,可以建立标准,对旋转设备的运行进行安全评定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010512019 CN102042903B (zh) | 2010-10-18 | 2010-10-18 | 一种基于有限元模型的旋转设备支承动刚度参数测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010512019 CN102042903B (zh) | 2010-10-18 | 2010-10-18 | 一种基于有限元模型的旋转设备支承动刚度参数测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102042903A true CN102042903A (zh) | 2011-05-04 |
CN102042903B CN102042903B (zh) | 2012-10-17 |
Family
ID=43909251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010512019 Expired - Fee Related CN102042903B (zh) | 2010-10-18 | 2010-10-18 | 一种基于有限元模型的旋转设备支承动刚度参数测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102042903B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706527A (zh) * | 2012-05-30 | 2012-10-03 | 哈尔滨工程大学 | 不平衡激扰力及不平衡激扰力矩联合作用下的设备内源激励载荷定量测试方法 |
CN102737146A (zh) * | 2012-06-12 | 2012-10-17 | 中国科学院工程热物理研究所 | 一种估算转子临界转速的工程方法 |
CN102928222A (zh) * | 2012-09-29 | 2013-02-13 | 广东电网公司电力科学研究院 | 一种滑动轴承动力特性系数试验识别方法 |
CN103344385A (zh) * | 2013-07-04 | 2013-10-09 | 上海工程技术大学 | 一种基于有限元的曲轴动平衡的获取方法 |
CN103530468A (zh) * | 2013-10-23 | 2014-01-22 | 清华大学 | 一种考虑轴承刚度耦合非线性的多支撑轴系有限元方法 |
CN107389267A (zh) * | 2017-07-12 | 2017-11-24 | 哈尔滨工程大学 | 一种转子***动平衡激励识别方法 |
CN111898298A (zh) * | 2020-07-25 | 2020-11-06 | 福建工程学院 | 一种基于有限元模型数值仿真的pcr基座温度控制***的参数优化方法 |
CN111898299A (zh) * | 2020-07-25 | 2020-11-06 | 福建工程学院 | 一种基于有限元模型数值仿真的pcr基座制造参数优化方法 |
CN112240352A (zh) * | 2019-07-16 | 2021-01-19 | 广州汽车集团股份有限公司 | 一种轮毂轴承的有限元分析方法及其*** |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103292958B (zh) * | 2013-05-28 | 2015-12-02 | 西安交通大学 | 一种基于模型的转子无试重失衡参数辨识方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86101634A (zh) * | 1985-03-22 | 1986-09-17 | 通用电气公司 | 具有各向异性刚度的轴承 |
CN1132323A (zh) * | 1995-03-27 | 1996-10-02 | 东北电力学院 | 旋转机械动平衡的动刚度系数方法及其装置 |
-
2010
- 2010-10-18 CN CN 201010512019 patent/CN102042903B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86101634A (zh) * | 1985-03-22 | 1986-09-17 | 通用电气公司 | 具有各向异性刚度的轴承 |
CN1132323A (zh) * | 1995-03-27 | 1996-10-02 | 东北电力学院 | 旋转机械动平衡的动刚度系数方法及其装置 |
Non-Patent Citations (2)
Title |
---|
《中国科学院研究生院学报》 19940531 田炯等 非协调轴对称杂交应力元单元刚度矩阵的准确建立 第39-44页 1 第11卷, 第1期 2 * |
《液压与气动》 20091130 凌玲等 高速旋转静压轴承接头芯轴的静强度和刚度分析 第46-48页 1 , 第11期 2 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706527A (zh) * | 2012-05-30 | 2012-10-03 | 哈尔滨工程大学 | 不平衡激扰力及不平衡激扰力矩联合作用下的设备内源激励载荷定量测试方法 |
CN102706527B (zh) * | 2012-05-30 | 2014-08-06 | 哈尔滨工程大学 | 不平衡激扰力及不平衡激扰力矩联合作用下的设备内源激励载荷定量测试方法 |
CN102737146A (zh) * | 2012-06-12 | 2012-10-17 | 中国科学院工程热物理研究所 | 一种估算转子临界转速的工程方法 |
CN102928222A (zh) * | 2012-09-29 | 2013-02-13 | 广东电网公司电力科学研究院 | 一种滑动轴承动力特性系数试验识别方法 |
CN102928222B (zh) * | 2012-09-29 | 2015-05-20 | 广东电网公司电力科学研究院 | 一种滑动轴承动力特性系数试验识别方法 |
CN103344385A (zh) * | 2013-07-04 | 2013-10-09 | 上海工程技术大学 | 一种基于有限元的曲轴动平衡的获取方法 |
CN103530468A (zh) * | 2013-10-23 | 2014-01-22 | 清华大学 | 一种考虑轴承刚度耦合非线性的多支撑轴系有限元方法 |
CN103530468B (zh) * | 2013-10-23 | 2016-03-30 | 清华大学 | 一种考虑轴承刚度耦合非线性的多支撑轴系有限元方法 |
CN107389267A (zh) * | 2017-07-12 | 2017-11-24 | 哈尔滨工程大学 | 一种转子***动平衡激励识别方法 |
CN107389267B (zh) * | 2017-07-12 | 2019-02-01 | 哈尔滨工程大学 | 一种转子***动平衡激励识别方法 |
CN112240352A (zh) * | 2019-07-16 | 2021-01-19 | 广州汽车集团股份有限公司 | 一种轮毂轴承的有限元分析方法及其*** |
CN112240352B (zh) * | 2019-07-16 | 2022-05-31 | 广州汽车集团股份有限公司 | 一种轮毂轴承的有限元分析方法及其*** |
CN111898298A (zh) * | 2020-07-25 | 2020-11-06 | 福建工程学院 | 一种基于有限元模型数值仿真的pcr基座温度控制***的参数优化方法 |
CN111898299A (zh) * | 2020-07-25 | 2020-11-06 | 福建工程学院 | 一种基于有限元模型数值仿真的pcr基座制造参数优化方法 |
CN111898298B (zh) * | 2020-07-25 | 2022-07-05 | 福建工程学院 | 基于有限元模型的pcr基座温度控制***的参数优化方法 |
CN111898299B (zh) * | 2020-07-25 | 2022-08-09 | 福建工程学院 | 一种基于有限元模型数值仿真的pcr基座制造参数优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102042903B (zh) | 2012-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102042903A (zh) | 一种基于有限元模型的旋转设备支承动刚度参数测量方法 | |
CN102095574B (zh) | 滚动导轨结合面动态特性参数测试装置及其测试方法 | |
CN102928222B (zh) | 一种滑动轴承动力特性系数试验识别方法 | |
CN103364071B (zh) | 单点激光连续扫描测振的薄壁圆筒模态测试***及方法 | |
CN108827536B (zh) | 一种组合回转系装配在线动平衡检测***及检测方法 | |
CN102507205A (zh) | 一种检测航空发动机风扇叶片颤振故障的方法 | |
CN107192488B (zh) | 一种宽频带切削力测量装置 | |
CN110319922A (zh) | 基于叶尖定时原理的叶片扭转振动位移测量方法 | |
CN103292958B (zh) | 一种基于模型的转子无试重失衡参数辨识方法 | |
CN106248328B (zh) | 一种测量轴系结构扭转动柔度的间接方法 | |
CN108827539B (zh) | 一种多转子装配在线动平衡试验装置 | |
CN104101464A (zh) | 一种基于旋转坐标系的多轮盘转子动平衡检测方法 | |
CN102721462A (zh) | 旋转机械启停车过程波德图/奈奎斯特图的快速计算方法 | |
CN113190786B (zh) | 一种大型旋转装备利用多维装配参数的振动预测方法 | |
CN205520745U (zh) | 电主轴机电一体化检测*** | |
CN101105419A (zh) | 基于全柔性振动***的动平衡测量***及其测量方法 | |
CN102928611A (zh) | 基于二倍频获取汽轮机转子临界转速的方法 | |
CN207248474U (zh) | 特高压直流穿墙套管抗震试验装置 | |
CN102706527B (zh) | 不平衡激扰力及不平衡激扰力矩联合作用下的设备内源激励载荷定量测试方法 | |
Xiang et al. | Torsional vibration of a shafting system under electrical disturbances | |
CN102788662A (zh) | 无附加传感器的垂直回转体动平衡检测*** | |
CN107525643A (zh) | 特高压直流穿墙套管抗震试验装置及试验方法 | |
CN113418591A (zh) | 五轴机床工作台轴向振动双位非接触检测装置及预测方法 | |
Lu et al. | Research on vibration characteristics of multistage gears transmission system driven by internal and external excitation | |
CN110046441A (zh) | 旋转机械轴系疲劳强度检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121017 Termination date: 20171018 |
|
CF01 | Termination of patent right due to non-payment of annual fee |