CN119024002B - 大型超重力旋转机械结构的微振动测量校准方法 - Google Patents
大型超重力旋转机械结构的微振动测量校准方法 Download PDFInfo
- Publication number
- CN119024002B CN119024002B CN202411422948.6A CN202411422948A CN119024002B CN 119024002 B CN119024002 B CN 119024002B CN 202411422948 A CN202411422948 A CN 202411422948A CN 119024002 B CN119024002 B CN 119024002B
- Authority
- CN
- China
- Prior art keywords
- acceleration sensor
- acceleration
- mechanical structure
- frequency domain
- rotating mechanical
- 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
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种大型超重力旋转机械结构的微振动测量校准方法,首先计算电磁干扰引起的频域幅值,然后构建大型超重力旋转机械结构的振动方程模型,求解得到微振动引起的真实加速度,并进行不同高度的微振动真实加速度的校正,最后构建传感器校准方程,求解待校准加速度传感器的真实值。本发明可对大型超重力旋转机械结构吊篮内定位于不同高度的加速度传感器进行测量校准,提高超重力环境下,加速度传感器测量数据的准确性。
Description
技术领域
本发明涉及超重力离心技术领域,具体是一种大型超重力旋转机械结构的微振动测量校准方法。
背景技术
目前,大型超重力旋转机械结构在高速旋转的过程中会产生微振动,主要原因是机械在试验端和配重端无法做到完全的平衡,随着转速的不断增加,不平衡力会越来越大,振动也会越来越大,大型超重力旋转机械结构的机械部件在长时间的运行过程中会出现磨损和老化,导致配合间隙增大或者表面粗糙度增加,引起微振动,这种微振动为机械振动。同时由于大型超重力旋转机械结构电磁力的不均匀分布或电磁干扰,也可能引起转子和定子之间的微振动,这种微振动是电磁干扰引起的振动。
在大型超重力旋转机械结构旋转过程引起的微振动会对两侧吊篮内的一些精密试验操作(如超重力下显微镜操作)或者传感器的测量有着很大的影响。目前,超重力环境下,由于微振动本身产生的加速度对电阻式、电容式、压阻式mems加速度传感器有着很大的影响,且无法直接判别出微振动主要是由于电磁干扰引起的还是机械振动。
发明内容
本发明要解决的技术问题是提供一种大型超重力旋转机械结构的微振动测量校准方法,可对大型超重力旋转机械结构吊篮内定位于不同高度的加速度传感器进行测量校准,提高超重力环境下,加速度传感器测量数据的准确性。
本发明的技术方案为:
大型超重力旋转机械结构的微振动测量校准方法,具体包括有以下步骤:
(1)、首先在大型超重力旋转机械结构的吊篮中设置位于不同水平高度的加速度传感器,然后启动大型超重力旋转机械结构旋转,不同水平高度的加速度传感器实时采集频域幅值信号,由于不同水平高度的加速度传感器采集的频域幅值信号由电磁干扰和机械振动引起,不同水平高度的加速度传感器受到电磁干扰引起的频域幅值变化相等,不同水平高度的加速度传感器受到机械振动引起的频域幅值变化的比值等于高度平方的比值,由此建立方程,求解得到电磁干扰引起的频域幅值Yem(f);
(2)、构建大型超重力旋转机械结构的振动方程模型,具体见下式
(3),根据式(3)求解得到水平方向的真实加速度和垂直方向的真实加速度
式(3)中,m表示大型超重力旋转机械结构吊篮内配重块的质量;表示水平方向的真实速度,由求导得到;表示垂直方向的真实速度,由求导得到;x(t)表示水平方向的真实位移,由求导得到;y(t)表示垂直方向的真实位移,由求导得到;表示大型超重力旋转机械结构的综合阻尼系数;k表示大型超重力旋转机械结构的刚度,取值为107~109;kxy表示水平方向的耦合刚度系数,即y方向对x方向的耦合作用,取值为106~107;kyx表示垂直方向的耦合刚度系数,即x方向对y方向的耦合作用,取值为106~107;mu表示不平衡重量,取值为0~10吨;ω表示大型超重力旋转机械结构的旋转角速度;L表示大型超重力旋转机械结构的转臂长度;e表示大型超重力旋转机械结构的转轴偏心距;g表示重力加速度;Yem(t)表示电磁干扰引起的时域信号,由电磁干扰引起的频域幅值Yem(f)转换得到;
(3)、在大型超重力旋转机械结构的吊篮内放置模型箱,模型箱的高度为H,模型箱内底面处的加速度为α底,α底包括有水平方向的真实加速度和垂直方向的真实加速度构建校正方程求解模型箱内任意高度z处的加速度为α(z),具体见下式(4):
式(4)中,n为微振动随着高度变化的衰减指数,取值为1或2;
(4)、大型超重力旋转机械结构吊篮内待校准加速度传感器的测量值为P测(z,t),构建传感器校准方程求解待校准加速度传感器的真实值P真(z,t),具体见下式(5):
P测(z,t)=P真(z,t)+ρα(z) (5);
式(5)中,ρ为修正系数,反映了待校准加速度传感器对微振动的敏感度,取值为0.5~1。
所述的电磁干扰引起的频域幅值Yem(f)的求解过程为:首先在大型超重力旋转机械结构的吊篮内设置配重块,配重块的顶面上连接L型板,L型板的竖直部分上连接有三个位于不同水平高度的加速度传感器,三个加速度传感器分别为邻近L型板竖直部分底端的第一加速度传感器、邻近L型板竖直部分顶端的第三加速度传感器和位于第一加速度传感器和第三加速度传感器之间的第二加速度传感器,第一加速度传感器距离配重块顶面的垂直距离为L1,第二加速度传感器距离配重块顶面的垂直距离为L2,第三加速度传感器距离配重块顶面的垂直距离为L3,然后将大型超重力旋转机械结构进行旋转,第一加速度传感器、第二加速度传感器和第三加速度传感器实时记录频域信号,设定第一加速度传感器、第二加速度传感器、第三加速度传感器受到电磁干扰的频域幅值均为Yem(f)、受到机械振动的频域幅值分别为Y1(f)、Y2(f)和Y3(f),然后建立方程,具体见下式(1):
式(1)中,Y实1(f)、Y实2(f)和Y实3(f)分别为第一加速度传感器、第二加速度传感器和第三加速度传感器实时记录的频域信号,不同高度的加速度传感器受到机械振动的频域幅值Y1(f)、Y2(f)、Y3(f)满足下式(2);
式(2)中,β为待求解的振动常数;
然后将式(2)带入式(1),求解得到振动常数β和电磁干扰的频域幅值Yem(f)。
所述的第一加速度传感器、第二加速度传感器和第三加速度传感器均为多个,即L型板的竖直部分上连接有三排加速度传感器,Y实1(f)为多个位于同一高度的第一加速度传感器实时记录的频域信号的平均值,Y实2(f)为多个位于同一高度的第二加速度传感器实时记录的频域信号的平均值,Y实3(f)为多个位于同一高度的第三加速度传感器实时记录的频域信号的平均值。
所述的由下式(6)计算得到:
式(6)中,f0表示线性阻尼系数,取值为1~10,c1和c2是非线性阻尼系数,取值范围都为0.1~10。
所述的待校准加速度传感器为IEPE加速度传感器。
本发明的优点:
(1)、对于机械振动,随着大型超重力旋转机械结构吊篮内加速度传感器定位水平高度的增加,振动幅度会放大;而电磁干扰引起的振动则对不同水平高度的加速度传感器,响应幅度基本相同,因为电磁干扰是环境中的干扰信号,对不同位置的加速度传感器影响相同,基于上述原理,建立方程,求解得到电磁干扰的频域幅值和受到机械振动的频域幅值,获知得到主要引起微振动的原因。
(2)、本发明根据大型超重力旋转机械结构转臂模态的不同,构建振动方程模型,求解微振动下的真实加速度,便于对不同高度微振动下的真实加速度进行校正,并最终建立传感器校准方程,实现加速度传感器测量值的校准,从而消除振动对实验测量的干扰,确保加速度传感器读数的准确性。
附图说明
图1是本发明L型板连接于配重块上的结构示意图。
图2是本发明L型板上不同高度加速度传感器的结构示意图。
图3是本发明配重块置于大型超重力旋转机械结构吊篮内的结构示意图。
附图标记:1-吊篮,2-配重块,3-L型板,4-第一加速度传感器,5-第二加速度传感器,6-第三加速度传感器。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
大型超重力旋转机械结构的微振动测量校准方法,具体包括有以下步骤:
(1)、见图1-图3,首先在大型超重力旋转机械结的吊篮1内设置配重块2,采用螺栓在配重块2的顶面上连接L型板3,L型板3的竖直部分上连接有三排位于不同水平高度的加速度传感器,三排加速度传感器分别为邻近L型板3竖直部分底端的一排第一加速度传感器4、邻近L型板3竖直部分顶端的一排第三加速度传感器6和位于一排第一加速度传感器4和一排第三加速度传感器6之间的一排第二加速度传感器5,每排加速度传感器均包括有四个加速度传感器,第一加速度传感器4距离配重块2顶面的垂直距离为L1,第二加速度传感器5距离配重块2顶面的垂直距离为L2,第三加速度传感器6距离配重块2顶面的垂直距离为L3,然后将大型超重力旋转机械结构进行旋转,第一加速度传感器4、第二加速度传感器5和第三加速度传感器6实时记录频域信号,设定第一加速度传感器4、第二加速度传感器5、第三加速度传感器6受到电磁干扰的频域幅值均为Yem(f)、受到机械振动的频域幅值分别为Y1(f)、Y2(f)和Y3(f),然后建立方程,具体见下式(1):
式(1)中,Y实1(f)、Y实2(f)和Y实3(f)分别为第一加速度传感器4、第二加速度传感器5和第三加速度传感器6实时记录的频域信号;Y实1(f)为四个位于同一高度的第一加速度传感器4实时记录的频域信号的平均值,Y实2(f)为四个位于同一高度的第二加速度传感器5实时记录的频域信号的平均值,Y实3(f)为四个位于同一高度的第三加速度传感器6实时记录的频域信号的平均值,满足下式(2);
式(2)中,β为待求解的振动常数;
然后将式(2)带入式(1),求解得到振动常数β和电磁干扰的频域幅值Yem(f)
(2)、构建大型超重力旋转机械结构的振动方程模型,具体见下式
(3),根据式(3)求解得到水平方向的真实加速度和垂直方向的真实加速度
式(3)中,m表示大型超重力旋转机械结构吊篮内配重块的质量;表示水平方向的真实速度,由求导得到;表示垂直方向的真实速度,由求导得到;x(t)表示水平方向的真实位移,由求导得到;y(t)表示垂直方向的真实位移,由求导得到;k表示大型超重力旋转机械结构的刚度,取值为107~109;kxy表示水平方向的耦合刚度系数,即y方向对x方向的耦合作用,取值为106~107;kyx表示垂直方向的耦合刚度系数,即x方向对y方向的耦合作用,取值为106~107;mu表示不平衡重量,取值为0~10吨;ω表示大型超重力旋转机械结构的旋转角速度;L表示大型超重力旋转机械结构的转臂长度;e表示大型超重力旋转机械结构的转轴偏心距;g表示重力加速度;Yem(t)表示电磁干扰引起的时域信号,由电磁干扰引起的频域幅值Yem(f)转换得到;表示大型超重力旋转机械结构的综合阻尼系数,由下式(6)计算得到:
式(6)中,c0表示线性阻尼系数,取值为1~10,c1和c2是非线性阻尼系数,取值范围都为0.1~10;
(3)、在大型超重力旋转机械结构的吊篮内放置模型箱,模型箱的高度为H,模型箱内底面处的加速度为α底,α底包括有水平方向的真实加速度和垂直方向的真实加速度构建校正方程求解模型箱内任意高度z处的加速度为α(z),具体见下式(4):
式(4)中,n为微振动随着高度变化的衰减指数,取值为1或2;
(4)、大型超重力旋转机械结构吊篮内待校准加速度传感器(电阻式、电容式或压阻式)的测量值为P测(z,t),构建传感器校准方程求解待校准加速度传感器的真实值P真(z,t),具体见下式(5):
P测(z,t)=P真(z,t)+ρα(z) (5);
式(5)中,ρ为修正系数,反映了待校准加速度传感器对微振动的敏感度,取值为0.5~1。
实施例
在大型超重力旋转机械结构(超重力离心机)的离心范围150g的情况下,L型板第一加速度传感器距离配重块顶面的垂直距离L1=5cm,第二加速度传感器距离配重块顶面的垂直距离L2=10cm,第三加速度传感器距离配重块顶面的垂直距离L3=15cm,实时测得的Y实1(f)大致为0.07g,Y实2(f)大致为0.22g,Y实3(f)大致为0.47g;
由式(1)和式(2)得到:从而求解出β=0.002,Yem(f)=0.02,即电磁干扰的频域幅值为0.02,不同高度受到机械振动的频域幅值分别为0.05、0.2、0.45,即机械振动的影响大、电磁干扰小。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (5)
1.大型超重力旋转机械结构的微振动测量校准方法,其特征在于:具体包括有以下步骤:
(1)、首先在大型超重力旋转机械结构的吊篮中设置位于不同水平高度的加速度传感器,然后启动大型超重力旋转机械结构旋转,不同水平高度的加速度传感器实时采集频域幅值信号,由于不同水平高度的加速度传感器采集的频域幅值信号由电磁干扰和机械振动引起,不同水平高度的加速度传感器受到电磁干扰引起的频域幅值变化相等,不同水平高度的加速度传感器受到机械振动引起的频域幅值变化的比值等于高度平方的比值,由此建立方程,求解得到电磁干扰引起的频域幅值Yem(f);
(2)、构建大型超重力旋转机械结构的振动方程模型,具体见下式(3),根据式(3)求解得到水平方向的真实加速度和垂直方向的真实加速度
式(3)中,m表示大型超重力旋转机械结构吊篮内配重块的质量;表示水平方向的真实速度,由求导得到;表示垂直方向的真实速度,由求导得到;x(t)表示水平方向的真实位移,由求导得到;y(t)表示垂直方向的真实位移,由求导得到;表示大型超重力旋转机械结构的综合阻尼系数;k表示大型超重力旋转机械结构的刚度,取值为107~109;kxy表示水平方向的耦合刚度系数,即y方向对x方向的耦合作用,取值为106~107;kyx表示垂直方向的耦合刚度系数,即x方向对y方向的耦合作用,取值为106~107;mu表示不平衡重量,取值为0~10吨;ω表示大型超重力旋转机械结构的旋转角速度;L表示大型超重力旋转机械结构的转臂长度;e表示大型超重力旋转机械结构的转轴偏心距;g表示重力加速度;Yem(t)表示电磁干扰引起的时域信号,由电磁干扰引起的频域幅值Yem(f)转换得到;
(3)、在大型超重力旋转机械结构的吊篮内放置模型箱,模型箱的高度为H,模型箱内底面处的加速度为α底,α底包括有水平方向的真实加速度和垂直方向的真实加速度构建校正方程求解模型箱内任意高度z处的加速度为α(z),具体见下式(4):
式(4)中,n为微振动随着高度变化的衰减指数,取值为1或2;
(4)、大型超重力旋转机械结构吊篮内待校准加速度传感器的测量值为P测(z,t),构建传感器校准方程求解待校准加速度传感器的真实值P真(z,t),具体见下式(5):
P测(z,t)=P真(z,t)+ρα(z) (5);
式(5)中,ρ为修正系数,反映了待校准加速度传感器对微振动的敏感度,取值为0.5~1。
2.根据权利要求1所述的大型超重力旋转机械结构的微振动测量校准方法,其特征在于:所述的电磁干扰引起的频域幅值Yem(f)的求解过程为:首先在大型超重力旋转机械结构的吊篮内设置配重块,配重块的顶面上连接L型板,L型板的竖直部分上连接有三个位于不同水平高度的加速度传感器,三个加速度传感器分别为邻近L型板竖直部分底端的第一加速度传感器、邻近L型板竖直部分顶端的第三加速度传感器和位于第一加速度传感器和第三加速度传感器之间的第二加速度传感器,第一加速度传感器距离配重块顶面的垂直距离为L1,第二加速度传感器距离配重块顶面的垂直距离为L2,第三加速度传感器距离配重块顶面的垂直距离为L3,然后将大型超重力旋转机械结构进行旋转,第一加速度传感器、第二加速度传感器和第三加速度传感器实时记录频域信号,设定第一加速度传感器、第二加速度传感器、第三加速度传感器受到电磁干扰的频域幅值均为Yem(f)、受到机械振动的频域幅值分别为Y1(f)、Y2(f)和Y3(f),然后建立方程,具体见下式(1):
式(1)中,Y实1(f)、Y实2(f)和Y实3(f)分别为第一加速度传感器、第二加速度传感器和第三加速度传感器实时记录的频域信号,不同高度的加速度传感器受到机械振动的频域幅值Y1(f)、Y2(f)、Y3(f)满足下式(2);
式(2)中,β为待求解的振动常数;
然后将式(2)代入式(1),求解得到振动常数β和电磁干扰的频域幅值Yem(f)。
3.根据权利要求2所述的大型超重力旋转机械结构的微振动测量校准方法,其特征在于:所述的第一加速度传感器、第二加速度传感器和第三加速度传感器均为多个,即L型板的竖直部分上连接有三排加速度传感器,Y实1(f)为多个位于同一高度的第一加速度传感器实时记录的频域信号的平均值,Y实2(f)为多个位于同一高度的第二加速度传感器实时记录的频域信号的平均值,Y实3(f)为多个位于同一高度的第三加速度传感器实时记录的频域信号的平均值。
4.根据权利要求1所述的大型超重力旋转机械结构的微振动测量校准方法,其特征在于:所述的由下式(6)计算得到:
式(6)中,c0表示线性阻尼系数,取值为1~10,c1和c2是非线性阻尼系数,取值范围都为0.1~10。
5.根据权利要求1所述的大型超重力旋转机械结构的微振动测量校准方法,其特征在于:所述的待校准加速度传感器为IEPE加速度传感器。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411422948.6A CN119024002B (zh) | 2024-10-12 | 2024-10-12 | 大型超重力旋转机械结构的微振动测量校准方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411422948.6A CN119024002B (zh) | 2024-10-12 | 2024-10-12 | 大型超重力旋转机械结构的微振动测量校准方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN119024002A CN119024002A (zh) | 2024-11-26 |
| CN119024002B true CN119024002B (zh) | 2025-02-25 |
Family
ID=93540337
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411422948.6A Active CN119024002B (zh) | 2024-10-12 | 2024-10-12 | 大型超重力旋转机械结构的微振动测量校准方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN119024002B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102589681A (zh) * | 2012-04-05 | 2012-07-18 | 邓昌建 | 一种高可靠旋转设备状态监测的无线振动测量方法与装置 |
| CN116699177A (zh) * | 2023-06-06 | 2023-09-05 | 华中科技大学 | 一种加速度计性能测试装置、方法及系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5806031B2 (ja) * | 2011-07-29 | 2015-11-10 | 株式会社ミツトヨ | 長周期振動感知器、及び長周期振動感知器の出力値補正方法 |
| WO2020160852A1 (en) * | 2019-02-04 | 2020-08-13 | Asml Netherlands B.V. | Electronic system, accelerometer, calibration method, lithographic apparatus and device manufacturing method |
| CN117346877A (zh) * | 2023-10-07 | 2024-01-05 | 清华大学 | 一种用于机械设备转速测量的离散频谱校正方法及装置 |
| CN118624216A (zh) * | 2024-05-29 | 2024-09-10 | 西安热工研究院有限公司 | 一种风力发电机组轴系振动测量装置检测方法 |
-
2024
- 2024-10-12 CN CN202411422948.6A patent/CN119024002B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102589681A (zh) * | 2012-04-05 | 2012-07-18 | 邓昌建 | 一种高可靠旋转设备状态监测的无线振动测量方法与装置 |
| CN116699177A (zh) * | 2023-06-06 | 2023-09-05 | 华中科技大学 | 一种加速度计性能测试装置、方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN119024002A (zh) | 2024-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106839968B (zh) | 转子空间弯曲轴线测试系统及其测试方法 | |
| CN103115726B (zh) | 一种基于应变的旋转零部件动平衡方法 | |
| JPS60214228A (ja) | 振動分析方法 | |
| CN101987718A (zh) | 具有低的热漂移的微机电z轴探测结构 | |
| US11372129B2 (en) | Post-compensation method for motion errors of rotating accelerometer gravity gradiometer | |
| CN113885098B (zh) | 一种重力敏感器低频频率响应误差在线建模及补偿方法 | |
| JP2013079856A (ja) | ジャイロ校正のための2回転軸直交ダブルターンテーブル | |
| CN108931824B (zh) | 一种动基座旋转加速度计重力梯度仪误差增益系数标定方法 | |
| CN106383366A (zh) | 一种旋转簧片式重力梯度仪 | |
| CN115931009B (zh) | 一种基于陀螺仪和激光测距的惯性器件离心测量方法 | |
| CN116678581B (zh) | 一种转子碰摩诱导内共振特性研究系统及其研究方法 | |
| CN119024002A (zh) | 大型超重力旋转机械结构的微振动测量校准方法 | |
| CN102564862B (zh) | 一种摆片刚度的测试方法 | |
| CN111289772A (zh) | 低深宽比的单质量块三轴mems惯性加速度计及其制备方法 | |
| CN115453146A (zh) | 电容式微机械加速度计 | |
| CN111337054B (zh) | 一种光纤陀螺动态特性测量修正方法 | |
| CN117805936A (zh) | 一种ragg线运动误差系数标定方法、补偿方法和系统 | |
| CN115164941B (zh) | 一种惯性导航设备振动测试装置及测试方法 | |
| CN115493747A (zh) | 一种基于应变技术的旋转机械轴系动平衡检测方法 | |
| CN212060315U (zh) | 一种超低频三轴向核电厂地震加速度仪标定台 | |
| CN111337053B (zh) | 一种光纤陀螺动态误差特性测量标定方法及标定系统 | |
| Zhang et al. | A new multi-disk shaft balancing method for rotating machinery using strain gauges | |
| JPS6140331B2 (zh) | ||
| CN216348671U (zh) | 一种改进的用于陀螺频率特性测量的角振动发生装置 | |
| Ameen et al. | Flexible rotor balancing without trial runs using experimentally tuned FE based rotor model |
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 |