CN111879400B - System and method for measuring module modal parameters of spacecraft electromechanical product - Google Patents
System and method for measuring module modal parameters of spacecraft electromechanical product Download PDFInfo
- Publication number
- CN111879400B CN111879400B CN202010628832.3A CN202010628832A CN111879400B CN 111879400 B CN111879400 B CN 111879400B CN 202010628832 A CN202010628832 A CN 202010628832A CN 111879400 B CN111879400 B CN 111879400B
- Authority
- CN
- China
- Prior art keywords
- modal
- rotor
- amplitude
- voltage signal
- module
- 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
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000013016 damping Methods 0.000 claims abstract description 41
- 238000005259 measurement Methods 0.000 claims abstract description 24
- 238000004458 analytical method Methods 0.000 claims description 65
- 238000001514 detection method Methods 0.000 claims description 40
- 238000012545 processing Methods 0.000 claims description 37
- 238000012360 testing method Methods 0.000 claims description 21
- 230000001133 acceleration Effects 0.000 claims description 18
- 238000001914 filtration Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 8
- 229920000742 Cotton Polymers 0.000 claims description 5
- 239000000725 suspension Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims 2
- 238000007667 floating Methods 0.000 abstract description 6
- 239000000306 component Substances 0.000 description 20
- 238000010586 diagram Methods 0.000 description 10
- 238000013461 design Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000008358 core component Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- 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
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域technical field
本发明涉及一种航天器机电产品的组件模态参数测量系统及方法,组件为陀螺的高速转子、控制力矩陀螺的高速转子或飞轮的转子,模态参数包括频率、振型和阻尼,属于小质量机械零组件模态测量技术领域。The invention relates to a component modal parameter measurement system and method for aerospace mechanical and electrical products. The component is a high-speed rotor of a gyro, a high-speed rotor of a control torque gyro or a rotor of a flywheel. The technical field of modal measurement of quality mechanical components.
背景技术Background technique
高速转子是航天器惯性执行机构、三浮陀螺、二浮陀螺等产品的核心部件。其模态设计影响整机稳定性和运行精度,而模态特性(频率、振型和阻尼)测量是验证模态设计准确性并进行设计迭代的重要手段。High-speed rotors are the core components of spacecraft inertial actuators, three-floating gyroscopes, and two-floating gyroscopes. Its modal design affects the stability and operation accuracy of the whole machine, and the measurement of modal characteristics (frequency, mode shape and damping) is an important means to verify the accuracy of modal design and carry out design iteration.
转子结构模态是重要设计要素,组件模态分析是整机模态分析的基础与核心之一,组件模态参数的获知精度直接影响整机模态分析的准确度。而小质量组件模态测量与实现更难。The modal of the rotor structure is an important design element, and the modal analysis of the components is one of the foundations and cores of the modal analysis of the whole machine. The accuracy of the modal parameters of the components directly affects the accuracy of the modal analysis of the whole machine. The modal measurement and implementation of small-mass components are more difficult.
CMG、陀螺等由多个零组件装配在一起,多零组件连接建模误差、零部件加工及装配误差均会给整机模态的高精度有限元分析带来困难,在大部分情况下,单纯的有限元分析对模态参数获取仅具有定性作用,这就需要对模态参数的进行精确测试。而测得的模态参数又可给设计分析提供反馈,用以修正有限元分析模型。CMG, gyroscope, etc. are assembled by multiple components. The connection modeling error of multiple components, component processing and assembly errors will bring difficulties to the high-precision finite element analysis of the modal of the whole machine. In most cases, Simple finite element analysis only has a qualitative effect on the acquisition of modal parameters, which requires accurate testing of modal parameters. The measured modal parameters can provide feedback to the design analysis to modify the finite element analysis model.
机械振动测试主要采用接触式加速度计来实现,测试传感器固定在被测件上。目前可获得的振动传感器一般为压电式或电荷式的,重量在10克以上,传感器自身质量会改变被测件结构,对测试精度有较大影响。另外,此类高精度振动传感器需要外配电荷放大器及高速数据采集器,成本极高,实现复杂。The mechanical vibration test is mainly realized by the contact accelerometer, and the test sensor is fixed on the DUT. Currently available vibration sensors are generally piezoelectric or charge-type, with a weight of more than 10 grams. The mass of the sensor itself will change the structure of the tested object, which has a great impact on the test accuracy. In addition, such high-precision vibration sensors require external charge amplifiers and high-speed data collectors, which are extremely costly and complicated to implement.
急需实现一种对被测组件质量分布影响小(特别对于小质量组件),且检测精度高,重量轻,量程高且便于实现的高精度组件模态参数测量方法。There is an urgent need to realize a high-precision component modal parameter measurement method that has little influence on the mass distribution of the component under test (especially for small-mass components), and has high detection accuracy, light weight, high range and easy implementation.
发明内容SUMMARY OF THE INVENTION
本发明解决的技术问题是:克服现有技术的不足,提供了一种航天器机电产品的组件模态参数测量系统及方法。The technical problem solved by the present invention is: to overcome the deficiencies of the prior art, and to provide a component modal parameter measurement system and method of a spacecraft electromechanical product.
本发明的技术解决方案是:The technical solution of the present invention is:
一种航天器机电产品的组件模态参数测量系统,该测量系统包括振动检测模块、信号处理模块、模态振型分析模块和模态阻尼分析模块;A component modal parameter measurement system of a spacecraft electromechanical product, the measurement system comprising a vibration detection module, a signal processing module, a modal mode shape analysis module and a modal damping analysis module;
所述的振动检测模块为基于加速度芯片的PCB式微型振动测量模块,包括加速度芯片、滤波模块和PCB板,加速度芯片、滤波模块固定在PCB板上,加速度芯片选择具有体积重量小、带宽高、量程宽的MEMS加速度芯片,以CMG转子为例,所述振动检测模块通过粘连方式固定在被测转子上,粘贴位置为轴承支承点(两个)和由Ansys预分析确定的转子模态振型的节点、最大值点等,用棉线(或柔软细线)悬挂固定CMG转子,固定点为转子模态振型的节点;The vibration detection module is a PCB-type micro vibration measurement module based on an acceleration chip, including an acceleration chip, a filter module and a PCB board. The acceleration chip and the filter module are fixed on the PCB board. A MEMS acceleration chip with a wide range, taking a CMG rotor as an example, the vibration detection module is fixed on the measured rotor by means of adhesion, and the sticking positions are the bearing support points (two) and the rotor modal vibration shape determined by Ansys pre-analysis The node, maximum point, etc., are fixed with cotton thread (or soft thin thread) to fix the CMG rotor, and the fixed point is the node of the modal vibration shape of the rotor;
当用模态锤敲击CMG转子端部(振型最大值点处)时,位于粘贴点的振动检测模块中的加速度芯片敏感到振动后将振动转化为电压信号后输出给滤波模块,滤波模块对接收到的电压信号进行滤波处理输出滤波后的电压信号给信号处理模块;When the end of the CMG rotor (at the maximum point of the mode shape) is struck with a modal hammer, the acceleration chip in the vibration detection module located at the sticking point is sensitive to the vibration and converts the vibration into a voltage signal and outputs it to the filter module. Filter the received voltage signal and output the filtered voltage signal to the signal processing module;
所述的信号处理模块用于接收振动检测模块输出的滤波后的电压信号,并对接收到的滤波后的电压信号进行FFT分析,得到滤波后的电压信号峰值频率和对应幅值,并将得到的滤波后的电压信号峰值频率和对应幅值输出给模态振型分析模块和模态阻尼分析模块中;The signal processing module is used to receive the filtered voltage signal output by the vibration detection module, and perform FFT analysis on the received filtered voltage signal to obtain the peak frequency and corresponding amplitude of the filtered voltage signal, and will obtain The peak frequency and corresponding amplitude of the filtered voltage signal are output to the modal mode shape analysis module and the modal damping analysis module;
所述的模态振型分析模块用于接收信号处理模块发送的滤波后的电压信号峰值频率和对应幅值,并根据接收到的不同粘贴点的在同一时刻的峰值频率的幅值进行描点和三次样条插值计算,得到转子模态振型;The modal shape analysis module is used to receive the filtered voltage signal peak frequency and corresponding amplitude sent by the signal processing module, and to plot and sum the received peak frequencies of different sticking points at the same time. Cubic spline interpolation calculation to obtain the rotor mode shape;
所述的模态阻尼分析模块用于接收信号处理模块发送的滤波后的电压信号峰值频率和对应幅值,并根据接收到的同一粘贴点在不同时刻的峰值频率的幅值变化进行分析运算,得到转子模态阻尼;The modal damping analysis module is used for receiving the filtered peak frequency and corresponding amplitude of the voltage signal sent by the signal processing module, and performing analysis and calculation according to the amplitude change of the received peak frequency of the same sticking point at different times, get the rotor modal damping;
所述的滤波模块由RC滤波实现,滤波带宽为所测峰值频率频率的5-10倍,加速度芯片感应被测件振动力,并转化为振动电信号;The filtering module is realized by RC filtering, the filtering bandwidth is 5-10 times of the measured peak frequency, and the acceleration chip senses the vibration force of the tested piece and converts it into a vibration electrical signal;
所述的信号处理模块,以四通道示波器为例,以轴承支承点位置粘贴的振动检测模块输出的电压信号U1为基准信号接入示波器通道1,对电压信号U1进行FFT变换,得到转子的模态频率F1的幅值并记为M11,将另一个轴承支承点位置粘贴的振动检测模块输出的电压信号U2为基准信号接入示波器通道2,以转子模态振型的节点位置粘贴的振动检测模块输出的电压信号U3为基准信号接入示波器通道3,以转子模态振型的最大值点位置粘贴的振动检测模块输出的电压信号U4为基准信号接入示波器通道4,并对电压信号U2、电压信号U3和电压信号U4分别进行FFT变换,分别得到转子的模态频率F1的幅值并记为M12、转子的模态频率F1的幅值并记为M13、转子的模态频率F1的幅值并记为M14,然后在时域上对示波器通道2输入的电压信号U2、示波器通道3输入的电压信号U3、示波器通道4输入的电压信号U4分别与示波器通道1输入的电压信号U1在频率F1尺度上的相位进行比较,如果同相则不进行任何处理,如果反相则在对应模态频率的幅值前加负号;The signal processing module, taking the four-channel oscilloscope as an example, takes the voltage signal U1 output by the vibration detection module pasted at the bearing support point as the reference signal and connects to the oscilloscope channel 1, and performs FFT transformation on the voltage signal U1 to obtain the model of the rotor. The amplitude of the modal frequency F1 is recorded as M11, and the voltage signal U2 output by the vibration detection module pasted at the position of another bearing support point is used as the reference signal to connect to the oscilloscope channel 2, and the vibration detection is pasted at the node position of the rotor modal mode shape. The voltage signal U3 output by the module is connected to the oscilloscope channel 3 as the reference signal, and the voltage signal U4 output by the vibration detection module pasted at the maximum point position of the rotor modal mode shape is connected to the oscilloscope channel 4 as the reference signal, and the voltage signal U2 is connected to the oscilloscope channel 4. , the voltage signal U3 and the voltage signal U4 are respectively subjected to FFT transformation to obtain the amplitude of the modal frequency F1 of the rotor and denote it as M12, the amplitude of the modal frequency F1 of the rotor and denote it as M13, and the modal frequency of the rotor F1. The amplitude is recorded as M14, and then in the time domain, the voltage signal U2 input by oscilloscope channel 2, the voltage signal U3 input by oscilloscope channel 3, and the voltage signal U4 input by oscilloscope channel 4 and the voltage signal U1 input by oscilloscope channel 1 are respectively in the time domain. The phases on the scale of frequency F1 are compared, if they are in phase, no processing is performed, and if they are out of phase, a negative sign is added before the amplitude of the corresponding modal frequency;
用如上步骤可测出模态频率F2、F3……的各测点幅值;Using the above steps, the amplitude of each measuring point of the modal frequency F2, F3... can be measured;
若有剩余测点,按如上步骤进行测试,直至测完为止,得到M15、M16……;If there are remaining measuring points, follow the above steps to test until the measurement is completed, and get M15, M16...;
所述的模态振型分析模块首先对信号处理模块所得到的各测试点对应测试模态频率的振动幅值进行归一化处理,将幅值M11参考标准化为1,得到振幅序列[1,M12’,M13’,……],用三次样条插值函数拟合振幅序列[1,M12’,M12’,……],即幅值M11的归一化振幅为1,幅值M12的归一化振幅为M12’,以此类推;The modal shape analysis module firstly normalizes the vibration amplitude of each test point corresponding to the test modal frequency obtained by the signal processing module, and normalizes the amplitude M11 reference to 1 to obtain an amplitude sequence [1, M12', M13', ...], use the cubic spline interpolation function to fit the amplitude sequence [1, M12', M12', ...], that is, the normalized amplitude of the amplitude M11 is 1, and the normalized amplitude of the amplitude M12 is 1. The normalized amplitude is M12', and so on;
转子材料本身应变滞后应力会形成材料内阻,需要进行测试,所述模态阻尼分析模块采用时域方法,测试转子振动衰减速度,得到相应振动模态的阻尼,为消除锤击作用过程的扰动,从振动开始后一段时间开始(>5ms)开始计时,模态振幅随时间变化关系为:The strain hysteresis stress of the rotor material itself will form the material internal resistance, which needs to be tested. The modal damping analysis module adopts the time domain method to test the rotor vibration attenuation speed to obtain the damping of the corresponding vibration mode, in order to eliminate the disturbance of the hammering action process. , the timing starts from a period of time (>5ms) after the vibration starts, and the relationship between the modal amplitude and time is:
其中F1为转子模态频率、ξ为相应模态阻尼,A(t)为在F1模态频率下幅值在t时刻的U(t)。Among them, F1 is the rotor modal frequency, ξ is the corresponding modal damping, and A(t) is the U(t) of the amplitude at time t under the F1 modal frequency.
同理可求出在模态频率为F2、F3……下响应的模态阻尼。In the same way, the modal damping of the response at the modal frequencies F2, F3... can be obtained.
一种航天器机电产品的组件模态参数测量方法,该方法的步骤包括:A method for measuring modal parameters of components of a spacecraft electromechanical product, the steps of the method include:
(1)Ansys预分析确定出转子模态振型的节点、最大值点等。振动检测模块粘贴在轴承支承点和转子模态振型的节点、最大值点处,用棉线(或柔软细线)悬挂固定转子,一端固定点为转子模态振型的节点,一端为悬吊点;(1) Ansys pre-analysis determines the nodes and maximum points of the rotor mode shape. The vibration detection module is pasted on the bearing support point and the node and maximum point of the modal mode shape of the rotor, and the rotor is suspended and fixed with cotton threads (or soft thin wires). point;
(2)用模态锤敲击CMG转子端部(振型最大值点处),位于不同粘贴点的振动检测模块将转子响应转化为的电压信号滤波后同时送给信号处理模块;有两种类型数据,一组为不同粘贴点同一时刻的响应信号;一组为每个粘贴点锤击后随时间变化的响应信号;(2) Hit the end of the CMG rotor (at the maximum point of the mode shape) with a modal hammer, and the vibration detection modules located at different sticking points filter the voltage signal converted by the rotor response and send it to the signal processing module at the same time; there are two kinds of Type data, one group is the response signals of different sticking points at the same time; one group is the response signals that change with time after each sticking point is hammered;
(3)信号处理模块对接收振动到的检测模块输出的滤波后的电压信号,进行FFT分析,得到滤波后的电压信号峰值频率和对应幅值,并将得到的滤波后的电压信号峰值频率和对应幅值输出给模态振型分析模块和模态阻尼分析模块中;(3) The signal processing module performs FFT analysis on the filtered voltage signal output by the received vibration detection module to obtain the peak frequency and corresponding amplitude of the filtered voltage signal, and compares the obtained filtered voltage signal peak frequency and The corresponding amplitude is output to the modal mode shape analysis module and the modal damping analysis module;
(4)模态振型分析模块根据信号处理模块输出的电压信号峰值频率、对应幅值以及接收到的不同粘贴点的在同一时刻的峰值频率的幅值进行描点和三次样条插值计算,得到转子模态振型;(4) The mode shape analysis module performs point tracing and cubic spline interpolation calculation according to the peak frequency of the voltage signal output by the signal processing module, the corresponding amplitude and the received amplitude of the peak frequency of different sticking points at the same time. The modal shape of the rotor;
(5)模态阻尼分析模块根据信号处理模块输出的电压信号峰值频率、对应幅值以及根据接收到的同一粘贴点在不同时刻的峰值频率的幅值变化进行分析运算,通过如下公式得到转子模态阻尼;(5) The modal damping analysis module analyzes and calculates according to the peak frequency of the voltage signal output by the signal processing module, the corresponding amplitude and the amplitude change of the received peak frequency of the same sticking point at different times, and obtains the rotor modulus by the following formula state damping;
其中F1为转子模态频率、ξ为相应模态阻尼,A(t)为在F1模态频率下幅值在t时刻的U(t)。Among them, F1 is the rotor modal frequency, ξ is the corresponding modal damping, and A(t) is the U(t) of the amplitude at time t under the F1 modal frequency.
同理可求出在模态频率为F2、F3……下响应的模态阻尼。In the same way, the modal damping of the response at the modal frequencies F2, F3... can be obtained.
本发明与现有技术相比的优点在于:The advantages of the present invention compared with the prior art are:
(1)本发明的PCB式振动传感器是一种微型、重量轻、量程宽且精度高的mems振动传感器。(1) The PCB type vibration sensor of the present invention is a miniature, light-weight, wide-range and high-precision mems vibration sensor.
(2)本发明的PCB式振动传感器质量小,对被测组件质量分布影响小(特别对于小质量组件),检测精度高。(2) The PCB-type vibration sensor of the present invention has a small mass, has little influence on the mass distribution of the component under test (especially for small-mass components), and has high detection accuracy.
(3)本发明可得到被测结构体的模态频率、振型、阻尼等参数,可实现被测体模态的精确测量。(3) The present invention can obtain the modal frequency, mode shape, damping and other parameters of the measured structure, and can realize the accurate measurement of the modal of the measured body.
(4)本发明使用示波器代替传统加速度计的前置器和分析软件,利用简单数值解算实现获得模态频率、模态振型和模态阻尼,实现简单,大大降低了模态测量成本。(4) The present invention uses an oscilloscope to replace the traditional accelerometer's front end device and analysis software, and uses simple numerical calculation to achieve modal frequency, modal mode shape and modal damping, which is simple to implement and greatly reduces modal measurement costs.
(5)一种微型高精度组件模态参数测量方法,采用MEMS加速度计进行小质量组件如陀螺高速转子的模态测量,降低传感器重量对被测转子质量分布的影响,提高转子模态测量精度。使用示波器代替传统加速度计的前置器和分析软件,利用简单运算实现获得模态频率、模态振型和模态阻尼,大大降低了模态测量成本。可应用于控制力矩陀螺、飞轮、二浮陀螺、三浮陀螺等产品的转子等组件模态参数的精确测量。本发明亦可应用于高速电机转子的模态测量,市场前景广阔。(5) A method for measuring the modal parameters of miniature high-precision components, which uses MEMS accelerometers to measure the modal parameters of small-mass components such as gyro high-speed rotors, reduces the influence of the sensor weight on the mass distribution of the measured rotor, and improves the modal measurement accuracy of the rotor. . The oscilloscope is used to replace the traditional accelerometer's preprocessor and analysis software, and the modal frequency, modal mode shape and modal damping can be obtained by simple operation, which greatly reduces the modal measurement cost. It can be used to accurately measure the modal parameters of components such as rotors of products such as control moment gyroscopes, flywheels, two-floating gyroscopes, and three-floating gyroscopes. The invention can also be applied to the modal measurement of the rotor of the high-speed motor, and has a broad market prospect.
附图说明Description of drawings
图1为本发明振动检测模块粘结示意图;Fig. 1 is the schematic diagram of bonding of vibration detection module of the present invention;
图2为本发明振动检测模块结构框图;Fig. 2 is the structural block diagram of the vibration detection module of the present invention;
图3为本发明信号处理模块示意图;3 is a schematic diagram of a signal processing module of the present invention;
图4为本发明模态振型分析模块示意图;4 is a schematic diagram of a modal mode shape analysis module of the present invention;
图5为本发明模态阻尼分析模块示意图。FIG. 5 is a schematic diagram of a modal damping analysis module of the present invention.
具体实施方式Detailed ways
如图1所示,为本发明结构框图,包括振动检测模块、信号处理模块、模态振型分析模块、和模态阻尼分析模块。振动检测模块为基于加速度芯片的PCB式微型振动测量模块。包括加速度芯片、滤波模块,安装在PCB上。加速度芯片选择具有体积重量小、带宽高、量程宽)的MEMS加速度芯片。以CMG转子为例,所述振动检测模块通过粘连方式固定在被测转子上,粘贴位置为轴承支承点和由Ansys预分析确定的模态振型的节点、最大值点等,用模态锤敲击转子端部(振型最大值点处)。振动检测模块获取粘贴点的转子振动信号,传送到通用示波器中,示波器对振动检测模块所测振动信号进行FFT分析,得到振动信号的峰值频率和对应幅值,送入模态振型分析模块和模态阻尼分析模块中,模态振型分析模块根据信号处理模块所获取的不同测试点在同一时刻的峰值频率的幅值进行描点和三次样条插值计算,得到转子模态振型。模态阻尼分析模块根据信号处理模块所获取的同一测试点在不同时刻的峰值频率的幅值变化进行分析运算,得到转子模态阻尼。As shown in FIG. 1, it is a structural block diagram of the present invention, including a vibration detection module, a signal processing module, a modal mode shape analysis module, and a modal damping analysis module. The vibration detection module is a PCB-type micro vibration measurement module based on an acceleration chip. Including acceleration chip, filter module, mounted on PCB. The acceleration chip selects a MEMS acceleration chip with small volume and weight, high bandwidth and wide range). Taking the CMG rotor as an example, the vibration detection module is fixed on the tested rotor by adhesion, and the sticking position is the bearing support point and the node and maximum point of the modal vibration shape determined by Ansys pre-analysis. Tap the end of the rotor (at the point of maximum mode shape). The vibration detection module obtains the rotor vibration signal at the sticking point and transmits it to the general oscilloscope. The oscilloscope performs FFT analysis on the vibration signal measured by the vibration detection module to obtain the peak frequency and corresponding amplitude of the vibration signal, which are sent to the modal mode shape analysis module and In the modal damping analysis module, the modal vibration analysis module performs point plotting and cubic spline interpolation calculations based on the amplitudes of the peak frequencies of different test points obtained by the signal processing module at the same time, and obtains the rotor modal vibration shape. The modal damping analysis module performs analysis and calculation according to the amplitude changes of the peak frequencies of the same test point obtained by the signal processing module at different times, and obtains the rotor modal damping.
如图2所示,为本发明振动检测模块结构框图,包括加速度芯片、滤波模块,安装在PCB上。加速度芯片选择具有体积重量小、带宽高、量程宽的MEMS加速度芯片。其滤波模块由RC滤波实现,滤波带宽为所测模态频率的5-10倍。加速度芯片感应被测件振动力,并转化为振动电信号,送入到示波器模块。As shown in FIG. 2 , it is a structural block diagram of the vibration detection module of the present invention, including an acceleration chip and a filter module, which are installed on the PCB. The acceleration chip selects a MEMS acceleration chip with small volume and weight, high bandwidth and wide range. Its filtering module is realized by RC filtering, and the filtering bandwidth is 5-10 times of the measured modal frequency. The acceleration chip senses the vibration force of the DUT, converts it into a vibration electrical signal, and sends it to the oscilloscope module.
如图3所示,为本发明信号处理模块示意图,以四通道示波器为例,选择一端部测点的振动信号为基准信号接入示波器通道1,对信号进行FFT变换,从中获得所测模态频率F1的幅值并记为M11。将测点2-4接入示波器通道2-4,逐次对各信号进行FFT变换,从中获得所测模态频率的幅值记为M12-M14。并在时域上对通道2-4信号与通道1信号在频率F1尺度上的相位进行比较。如果反相,在对应模态频率的幅值前加负号。若有剩余测点,按如上步骤进行测试,直至测完为止,得到M15、M16……。用如上步骤可测出模态频率F2、F3……的各测点幅值。As shown in Figure 3, it is a schematic diagram of the signal processing module of the present invention. Taking a four-channel oscilloscope as an example, the vibration signal at one end of the measuring point is selected as the reference signal to connect to the oscilloscope channel 1, and the signal is subjected to FFT transformation to obtain the measured mode. The amplitude of frequency F1 is also recorded as M11. Connect the measuring point 2-4 to the oscilloscope channel 2-4, perform FFT transformation on each signal one by one, and obtain the amplitude of the measured modal frequency, which is recorded as M12-M14. And compare the phase of the channel 2-4 signal and the channel 1 signal on the frequency F1 scale in the time domain. If inverted, add a minus sign before the magnitude of the corresponding modal frequency. If there are remaining measuring points, follow the above steps to test until the measurement is completed, and get M15, M16... . The amplitude of each measuring point of modal frequency F2, F3... can be measured by the above steps.
如图4所示,为本发明模态振型分析模块示意图。首先对信号处理模块所得到的各测试点对应测试模态频率的振动幅值进行归一化处理,以基准信号为参考标准化为1,得到振幅序列[1,M11’,M12’,……],[1,M21’,M22’,……],……。As shown in FIG. 4 , it is a schematic diagram of the modal mode shape analysis module of the present invention. First, normalize the vibration amplitude of each test point corresponding to the test modal frequency obtained by the signal processing module, and normalize it to 1 with the reference signal as a reference, and obtain the amplitude sequence [1, M11', M12',...] , [1, M21', M22', ...], ....
然后,用三次样条插值函数拟合振幅序列[1,M11’,M12’,……],拟合后转子振型与原始数据对比。Then, use the cubic spline interpolation function to fit the amplitude sequence [1, M11', M12', ...], and compare the rotor mode shape with the original data after fitting.
如图5所示,为本发明模态阻尼分析模块示意图。转子材料本身应变滞后应力会形成材料内阻,需要进行测试。在此采用时域方法,测试转子振动衰减速度,得到相应振动模态的阻尼。为消除锤击作用过程的扰动,从振动开始后10mS开始计时。As shown in FIG. 5 , it is a schematic diagram of the modal damping analysis module of the present invention. The strain hysteresis stress of the rotor material itself will form the internal resistance of the material, which needs to be tested. Here, the time domain method is used to test the vibration attenuation speed of the rotor to obtain the damping of the corresponding vibration mode. In order to eliminate the disturbance of the hammering action process, the timing was started 10mS after the vibration started.
模态振幅随时间变化关系为:The modal amplitude changes with time as:
其中F1为转子模态频率、ξ为相应模态阻尼,A(t)为在F1模态频率下幅值在t时刻的U(t)。Among them, F1 is the rotor modal frequency, ξ is the corresponding modal damping, and A(t) is the U(t) of the amplitude at time t under the F1 modal frequency.
转子一阶模态F1振幅在110mS时刻比10mS时刻高11.3dB,由(1)式推导得到下式:The amplitude of the first-order modal F1 of the rotor is 11.3dB higher at the moment of 110mS than that at the moment of 10mS. The following formula is derived from formula (1):
解算得ξ=0.004776。Solved to get ξ=0.004776.
实施例Example
以CMG转子为例,首先通过Ansys预分析确定出转子模态振型的节点和最大值点;将振动检测模块粘贴在如图1所示轴承支承点1、轴承支承点2和转子模态振型的节点、最大值点处。用棉线(或柔软细线)悬挂固定转子,一端固定点为转子模态振型的节点,一端为悬吊点。Taking the CMG rotor as an example, the nodes and maximum points of the rotor modal vibration shape are firstly determined through Ansys pre-analysis; The node of the type, the maximum point. Use cotton thread (or soft thin thread) to suspend and fix the rotor. One end of the fixed point is the node of the rotor's mode shape, and the other end is the suspension point.
用模态锤敲击CMG转子端部(振型最大值点处),位于轴承支承点1的振动检测模块输出U1(t),其他粘贴点的振动检测模块输出U2(t)、U3(t)和U4(t)。将U1(t)、U2(t)、U3(t)和U4(t)同时送给信号处理模块。以轴承支承点1位置粘贴的振动检测模块输出的电压信号U1(t)为基准信号接入示波器通道1,U2(t)、U3(t)和U4(t)接入示波器2、3、4通道。Hit the end of the CMG rotor (at the maximum point of the mode shape) with a modal hammer, the vibration detection module at bearing support point 1 outputs U1(t), and the vibration detection modules at other sticking points output U2(t), U3(t) ) and U4(t). Send U1(t), U2(t), U3(t) and U4(t) to the signal processing module at the same time. Use the voltage signal U1(t) output by the vibration detection module pasted at the bearing support point 1 as the reference signal to connect to the oscilloscope channel 1, and U2(t), U3(t) and U4(t) to connect to the oscilloscope 2, 3, 4 aisle.
信号处理模块对接收到的振动检测模块输出的U1(t)、U2(t)、U3(t)和U4(t),进行FFT分析,得到滤波后的电压信号峰值频率和对应幅值,对U1(t)FFT分析找到幅值最高的频率得到转子的模态频率702Hz;幅值并记为M11=1.2V;对U2(t)FFT分析得到幅值并记为M12=1.2V;对U3(t)FFT分析得到幅值并记为M13=0.1V;对U4(t)FFT分析得到转子的幅值并记为M14=1.4V;然后在时域上对示波器通道2输入的电压信号U2、示波器通道3输入的电压信号U3、示波器通道4输入的电压信号U4分别与示波器通道1输入的电压信号U1在频率702Hz尺度上的相位进行比较,如果同相则不进行任何处理,如果反相则在对应模态频率的幅值前加负号。The signal processing module performs FFT analysis on the received U1(t), U2(t), U3(t) and U4(t) output by the vibration detection module, and obtains the peak frequency and corresponding amplitude of the filtered voltage signal. U1(t) FFT analysis finds the frequency with the highest amplitude to obtain the modal frequency of the rotor 702Hz; the amplitude is recorded as M11=1.2V; for U2(t) FFT analysis, the amplitude is obtained and recorded as M12=1.2V; for U3 (t) The amplitude of the FFT analysis is obtained and recorded as M13=0.1V; the amplitude of the rotor is obtained by the FFT analysis of U4 and recorded as M14=1.4V; then the voltage signal U2 input to the oscilloscope channel 2 in the time domain , The voltage signal U3 input by oscilloscope channel 3 and the voltage signal U4 input by oscilloscope channel 4 are compared with the phase of the voltage signal U1 input by oscilloscope channel 1 on the frequency scale of 702Hz. Add a minus sign before the magnitude of the corresponding modal frequency.
将得到的滤波后的电压信号峰值频率和对应幅值输出给模态振型分析模块和模态阻尼分析模块中;模态振型分析模块根据信号处理模块输出的电压信号U1(t)、U2(t)、U3(t)和U4(t)、峰值频率702Hz以及对应幅值M11=1.2V、M12=1.2V、M13=0.1V、M14=1.4V进行描点和三次样条插值计算,首先对信号处理模块所得到的各测试点对应测试模态频率的振动幅值进行归一化处理,将幅值M11参考标准化为1,得到振幅序列[1,1,0.1/1.2,1.4/1.2],即幅值M11的归一化振幅为1,幅值M12的归一化振幅为1,以此类推,得到转子模态振型。The obtained peak frequency and corresponding amplitude of the filtered voltage signal are output to the modal mode shape analysis module and modal damping analysis module; the modal mode shape analysis module is based on the voltage signals U1(t) and U2 output by the signal processing module. (t), U3(t) and U4(t), peak frequency 702Hz, and corresponding amplitudes M11=1.2V, M12=1.2V, M13=0.1V, M14=1.4V to perform point tracing and cubic spline interpolation calculations, first Normalize the vibration amplitude of each test point corresponding to the test modal frequency obtained by the signal processing module, normalize the amplitude M11 reference to 1, and obtain the amplitude sequence [1, 1, 0.1/1.2, 1.4/1.2] , that is, the normalized amplitude of the amplitude M11 is 1, the normalized amplitude of the amplitude M12 is 1, and so on, to obtain the rotor modal shape.
模态阻尼分析模块根据信号处理模块输出的电压信号峰值频率、对应幅值以及根据接收到的同一粘贴点在不同时刻的峰值频率的幅值变化通过如下公式得到转子模态阻尼;模态振幅随时间变化关系为:The modal damping analysis module obtains the rotor modal damping by the following formula according to the peak frequency and corresponding amplitude of the voltage signal output by the signal processing module, and the amplitude change of the received peak frequency at the same sticking point at different times; The time-varying relationship is:
转子在一阶模态F1=702Hz振幅在110mS时刻比10mS时刻高11.3dB,由(1)式推导得到下式:The rotor's first-order mode F1=702Hz amplitude is 11.3dB higher at 110mS time than 10mS time. The following formula is derived from formula (1):
解算得ξ=0.004776。Solved to get ξ=0.004776.
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。The content not described in detail in the specification of the present invention belongs to the well-known technology of those skilled in the art.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010628832.3A CN111879400B (en) | 2020-07-01 | 2020-07-01 | System and method for measuring module modal parameters of spacecraft electromechanical product |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010628832.3A CN111879400B (en) | 2020-07-01 | 2020-07-01 | System and method for measuring module modal parameters of spacecraft electromechanical product |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111879400A CN111879400A (en) | 2020-11-03 |
CN111879400B true CN111879400B (en) | 2022-07-05 |
Family
ID=73149859
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010628832.3A Active CN111879400B (en) | 2020-07-01 | 2020-07-01 | System and method for measuring module modal parameters of spacecraft electromechanical product |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111879400B (en) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102169046B (en) * | 2010-12-22 | 2012-07-25 | 北京航空航天大学 | Online elastic mode testing system for magnetically suspended electromechanical equipment |
DK201100234A (en) * | 2011-03-30 | 2012-10-01 | Brincker Rune | Method for improved estimation of one or more experimentally obtained mode shapes |
CN103983412B (en) * | 2014-05-30 | 2017-06-06 | 北京航空航天大学 | For vibrating FEM updating avionic device operation mode measuring method |
CN107271126B (en) * | 2017-05-31 | 2020-04-14 | 北京航空航天大学 | On-line test device for rotor modal frequency and mode shape based on magnetic levitation technology |
CN109743856A (en) * | 2018-12-29 | 2019-05-10 | 积高电子(无锡)有限公司 | A kind of packaging technology of sensor |
-
2020
- 2020-07-01 CN CN202010628832.3A patent/CN111879400B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN111879400A (en) | 2020-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102095574B (en) | Joint surface dynamic characteristic parameter testing device of rolling guide rail and testing method thereof | |
CN103256941A (en) | Practical method of high order temperature compensation for MEMS (Micro Electro Mechanical Systems) gyroscope | |
CN106568604A (en) | Calculation method for vibration isolation rate of automotive power assembly suspension system | |
CN108645583A (en) | One kind being fitted to each other face normal direction contact damping high-precision detection device and method | |
CN101566483A (en) | Vibration error compensation method of fiber optic gyro strap-down inertia measurement system | |
CN101158582A (en) | A MEMS gyroscope differential measurement method | |
CN103954304B (en) | A kind of MEMS of being applied to is used to zero inclined shot and long term changing value method of testing of group | |
CN111537170A (en) | Dynamic stiffness testing method for servo actuator | |
CN102175394B (en) | Permanent Calibration Method in Dynamic Unbalance Test of Rigid Rotor with Soft Bearing | |
CN113514666B (en) | A micromechanical accelerometer based on PT symmetric resonator and its detection method | |
CN102538941A (en) | Device and method for conventional balance to measure natural frequency of cantilever supported model in wind tunnel | |
Xu et al. | Dynamic modeling and compensation of robot six-axis wrist force/torque sensor | |
CN111879400B (en) | System and method for measuring module modal parameters of spacecraft electromechanical product | |
CN107291988B (en) | Method for acquiring equivalent excitation force of momentum wheel mounting interface | |
CN203069261U (en) | Modal test system based on single-point laser continuous plane scanning vibration measurement | |
CN105222973A (en) | For the field calibration method of vibration signal order tracking technique | |
CN105157694B (en) | A kind of apparatus and method of accurate assessment optical fibre gyro output delay time | |
CN210271424U (en) | A self-test device for modal analysis | |
CN118583434A (en) | A pulse wind tunnel balance with self-compensation for inertial disturbance and compensation method thereof | |
CN116296237B (en) | Vertical wind load test method for large carrier rocket in low-speed wind tunnel | |
CN109239403B (en) | Single-device virtual accelerometer based on time measurement and implementation method thereof | |
CN116499696B (en) | Method for improving dynamic accuracy of attitude angle of wind tunnel model test model | |
CN107941192A (en) | A kind of angle of rudder reflection tester | |
CN116155143B (en) | Method and device for acquiring motor parameters | |
CN116499665A (en) | System and method for dynamic measurement of nonlinear stiffness of reed return mechanism of low-frequency vibration table |
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 |