CN109884985A - The measurement method of numerically-controlled machine tool complete machine machining state dynamic characteristic - Google Patents

The measurement method of numerically-controlled machine tool complete machine machining state dynamic characteristic Download PDF

Info

Publication number
CN109884985A
CN109884985A CN201910202404.1A CN201910202404A CN109884985A CN 109884985 A CN109884985 A CN 109884985A CN 201910202404 A CN201910202404 A CN 201910202404A CN 109884985 A CN109884985 A CN 109884985A
Authority
CN
China
Prior art keywords
machine tool
vibration
response
mode
numerical control
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.)
Pending
Application number
CN201910202404.1A
Other languages
Chinese (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.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and 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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201910202404.1A priority Critical patent/CN109884985A/en
Publication of CN109884985A publication Critical patent/CN109884985A/en
Pending legal-status Critical Current

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The present invention relates to a kind of measurement methods of numerically-controlled machine tool complete machine machining state dynamic characteristic, the specific steps are as follows: step 1: being respectively mounted vibration acceleration flowmeter sensor on each main component of numerically-controlled machine tool;Step 2: under machining state, the vibration generated using lathe foundation vibration and operation acquires the three-phase vibratory response of each component of numerically-controlled machine tool by multi-channel data acquisition board as excitation;Step 3: using Bayes's operational modal method and Fast Fourier Transform (FFT), obtains numerically-controlled machine tool each channel frequency response function curve under off working state environment, establishes Modal Parameter Identification function;Step 4: the dynamic characteristic parameters such as each rank intrinsic frequency of numerically-controlled machine tool, damping ratio and Mode Shape are obtained by analyzing each channel frequency response function.The present invention solves the problems, such as machining state complete machine tool dynamic characteristic measuring.The characteristics of this method is to be not necessarily to artificial excitation, and can get the dynamic characterization measurement result of complete machine tool machining state.

Description

Method for measuring dynamic characteristics of complete machine processing state of numerical control machine tool
Technical Field
The invention relates to a method for testing a numerical control machine tool and evaluating the dynamic characteristics of a machine tool structure, in particular to a method for measuring the dynamic characteristics of the complete machine processing state of the numerical control machine tool.
Background
The dynamic characteristics of the machine tool refer to the characteristics of a machine tool system in a vibration state, namely the characteristics that the amplitude and the phase of the machine tool change along with the vibration frequency under a certain vibration exciting force, and main indexes of the dynamic characteristics of the machine tool comprise natural frequency, damping ratio, modal shape, dynamic stiffness and the like. The dynamic characteristics of the machine tool determine the cutting characteristics of the machine tool and directly determine performance indexes such as processing stability, cutting capability and precision of the numerical control machine tool. Therefore, the measurement and evaluation of the dynamic characteristics of the machine tool have important significance for the design and the use of the numerical control machine tool. For measuring the dynamic characteristics of the machine tool, a hammering method is generally adopted, as shown in fig. 1(a) and (B), a pulse excitation signal is given to the machine tool by hammering a certain part a of the machine tool, then a plurality of vibration sensors B are arranged at different positions to measure the response signals of the pulse excitation signal, and the dynamic characteristics of the machine tool are obtained through a signal processing theory. However, the traditional experimental mode method based on the hammering method has the disadvantages that the single hammering of a complex system such as a machine tool is difficult to excite the output response, and the dynamic characteristic test and analysis of the complete machine machining state of the machine tool cannot be realized.
In the machining process of the numerical control machine tool, different cutting technological parameters can generate different vibration frequencies of the numerical control machine tool, and all parts of the numerical control machine tool have different vibration inherent frequencies. In the actual machining process, vibration excited by cutting technological parameters often occurs to cause resonance of key functional parts of the machine tool, so that the machining precision of the machine tool is greatly influenced. Therefore, the natural frequency of each order of the numerical control machine tool can be mastered by measuring the dynamic characteristics of the numerical control machine tool, and the vibration phenomenon of the numerical control machine tool in the machining process is avoided. On the other hand, in the machining state, the dynamic characteristics of the whole machine tool in the machining state are obtained, the dynamic design defects of the machine tool are positioned and evaluated, and the design improvement of the numerical control machine tool plays an important role. The method is also a precondition for evaluating structural design optimization of each main part of the machine tool and monitoring the performance of the machine tool.
Disclosure of Invention
The invention provides a method for measuring the dynamic characteristics of the complete machine processing state of a numerical control machine tool, which solves the problem of measuring the dynamic characteristics of the complete machine of the machine tool in the processing state by using the vibration of a machine tool foundation and the vibration generated by operation as excitation according to a modal parameter identification technology combining a Bayesian theory and operation modal analysis. The method is characterized in that manual excitation is not needed, and the dynamic characteristic test result of the complete machine processing state of the machine tool can be obtained, wherein the dynamic characteristic test result comprises main parameters such as natural frequency, damping ratio, vibration mode, signal-to-noise ratio and the like.
The technical scheme of the invention is as follows:
a method for measuring the dynamic characteristics of the complete machine processing state of a numerical control machine tool comprises the following specific steps:
the method comprises the following steps: a vibration accelerometer sensor is arranged on each main part of the numerical control machine tool;
step two: in a machining state, the vibration of a machine tool foundation and the vibration generated by operation are used as excitation, and the three-phase vibration response of each part of the numerical control machine tool is acquired through a multi-channel data acquisition card;
step three: obtaining frequency response function curves of all channels of the numerical control machine tool in a non-working state environment by using a Bayesian operation modal method and fast Fourier transform, and establishing a modal parameter identification function;
step four: and obtaining dynamic characteristic parameters of the numerical control machine tool, such as inherent frequency, damping ratio, modal shape and the like of each order by analyzing frequency response functions of each channel.
The principle and the method for acquiring the three-phase vibration response of each part of the numerical control machine tool through the multi-channel data acquisition card in the machining state are as follows:
the linear vibration system of the numerical control machine tool with n degrees of freedom meets the equation of motion:
wherein M is equal to Rn×n、C∈Rn×n、K∈Rn×nRespectively a mass matrix, a damping ratio matrix and a rigidity matrix of the system; the structural response of the multi-degree-of-freedom structure satisfying linear classical damping is assumed to be:
wherein,is a global mode shape vector of the j order;the mode response is the j order mode response, the mode order is m, and the mode response satisfies the decoupling mode power equation:
wherein, ω isj,ζjAnd Wj(t) are scalar quantities representing the natural frequency, damping ratio and modal force of the j-th order mode of the machine tool, respectively. The main parameters of the mode identification are natural frequency, damping ratio and vibration mode;
the numerical control machine tool generates weak vibration under the condition of random working condition excitationN is the sampling number in each acquisition channel for the vibration acceleration response data obtained by testing; the acceleration response obtained by the hypothesis test is composed of the structural environment vibration test signal and the prediction error of the model, i.e.
Wherein, { εj∈RnThe "prediction error" is the error between the model response and the test response due to the test noise and the model error. Carrying out Fast Fourier transform (Fast Fourier transform, FFT) on the acquired machine tool acceleration signal D, namely
In the formula i2=-1;FkAnd GkRespectively representThe real and imaginary parts of (c); Δ t is the sampling step, corresponding to k being 2,3, …, Nq,fkGiven (k-1)/(N Δ t), the augmentation matrix is assumedAn amplification matrix { Z obtained by FFT of acceleration signals acquired from a machining statekStarting from the parameter, identifying the dynamic characteristic modal parameter theta of the machine tool, wherein the parameter is f, zeta, Se(ii) a Φ, where f ═ fj:j=1,…,m},ζ={ζjJ-1, …, m-and Φ - Φ1,…,Φm]∈Rn×mRespectively representing the natural frequency and the damping ratio in the m-order natural mode and the corresponding integral vibration mode of the machine tool; s is belonged to Cm×mAnd SeRespectively representing the spectral density matrix of the state excitation and the prediction error spectral density;
the method can be obtained by the Bayesian theory principle:
p(θ|D)∝p(D|θ) (6)
p (θ | D) represents the posterior probability of θA density function; p (D | theta) is a likelihood function. From the stochastic process analysis we can see that when high sampling rate and sufficient data length are met, the FFT results are approximately independent and satisfy a gaussian distribution, so ZkSubject to a zero mean Gaussian combined distribution to obtain a likelihood function of
In the formula, det (-) represents determinant; ckIs ZkFor the purpose of analysis and calculation, the theoretical covariance matrix of (a) is a negative log-likelihood function, expressed as L (θ), and p (D | θ) ═ exp [ -L (θ) for analysis]Wherein
Thus, the optimized value of the modal parameter θ can be obtained by maximizing the posterior probability density function, i.e., minimizing the negative log likelihood function L (θ).
Compared with the prior art, the invention has the following beneficial effects:
the dynamic characteristics of the machine tool are generally measured by adopting a hammering method, a pulse excitation signal is given to the machine tool by hammering a certain part of the machine tool, then a plurality of vibration sensors are arranged at different positions to measure response signals of the pulse excitation signal, and the dynamic characteristics of the machine tool are obtained by a signal processing theory. However, the traditional experimental mode method based on the hammering method has the disadvantages that the single hammering of a complex system such as a machine tool is difficult to excite the output response, and the dynamic characteristic test and analysis of the complete machine machining state of the machine tool cannot be realized.
The invention can obtain the dynamic test data of the machine tool in the real machining state, and ensures the accuracy of the dynamic parameters of the machine tool.
Drawings
FIG. 1 is a schematic view of a conventional method for measuring dynamic characteristics of a numerically controlled machine tool;
wherein, (a) a schematic diagram and (b) an example diagram;
FIG. 2 is a schematic view of the method for measuring dynamic characteristics of a numerically controlled machine tool according to the present invention;
wherein, (a) a schematic diagram and (b) an example diagram;
FIG. 3 is a right view of the measurement point arrangement of a certain type of machine tool;
FIG. 4 is a left side view of a measuring point arrangement of a certain type of machine tool;
FIG. 5 shows a first-order vibration mode of a machine tool of a certain type;
FIG. 6 is a second-order vibration mode of a machine tool of a certain type;
FIG. 7 shows a third vibration mode of a machine tool of a certain type.
Detailed Description
The invention is further described with reference to the following figures and examples.
As shown in fig. 2(a) and (b), a method for measuring the dynamic characteristics of the complete machine processing state of a numerical control machine tool specifically comprises the following steps:
the method comprises the following steps: a vibration accelerometer sensor B is arranged on each main part of the numerical control machine tool;
step two: in a machining state, the three-phase vibration response of each part of the numerical control machine tool is acquired through a multi-channel data acquisition card by using the machine tool foundation vibration C and the vibration D generated by operation as excitation;
step three: obtaining frequency response function curves of all channels of the numerical control machine tool in a non-working state environment by using a Bayesian operation modal method and fast Fourier transform, and establishing a modal parameter identification function;
step four: and obtaining dynamic characteristic parameters of the numerical control machine tool, such as inherent frequency, damping ratio, modal shape and the like of each order by analyzing frequency response functions of each channel.
The measurement principle of the invention is as follows:
the linear vibration system of the numerical control machine tool with n degrees of freedom meets the equation of motion:
wherein M is equal to Rn×n、C∈Rn×n、K∈Rn×nRespectively a mass matrix, a damping ratio matrix and a rigidity matrix of the system. The structural response of the multi-degree-of-freedom structure satisfying linear classical damping is assumed to be:
wherein,is a global mode shape vector of the j order;the mode response is the j order mode response, the mode order is m, and the mode response satisfies the decoupling mode power equation:
wherein, ω isj,ζjAnd Wj(t) are scalar quantities representing the natural frequency, damping ratio and modal force of the j-th order mode of the machine tool, respectively. The main parameters of the mode identification are the natural frequency, the damping ratio and the mode shape.
Numerical control machine toolCan generate weak vibration under the excitation condition of random working conditionsAnd N is the sampling number in each acquisition channel for testing the obtained vibration acceleration response data. The acceleration response obtained by the hypothesis test is composed of the structural environment vibration test signal and the prediction error of the model, i.e.
Wherein, { εj∈RnThe "prediction error" is the error between the model response and the test response due to the test noise and the model error. Performing Fast Fourier Transform (FFT) on the acquired machine tool acceleration signal D, namely
In the formula i2=-1;FkAnd GkRespectively representThe real and imaginary parts of (c); Δ t is the sampling step. Corresponding to k 2,3, …, Nq,fk(k-1)/(N Δ t). Hypothesis augmenting matrixThe content of the study is an augmentation matrix { Z ] obtained by FFT of acceleration signals acquired from a processing statekStarting from the parameter, identifying the dynamic characteristic modal parameter theta of the machine tool, wherein the parameter is f, zeta, Se(ii) a Φ, where f ═ fj:j=1,…,m},ζ={ζjJ-1, …, m-and Φ - Φ1,…,Φm]∈Rn×mRespectively representing the natural frequency and the damping ratio in the m-order natural mode and the corresponding integral vibration mode of the machine tool; s is belonged to Cm×mAnd SeRespectively representing the spectral density matrix of the state excitation and the prediction error spectral density.
The method can be obtained by the Bayesian theory principle:
p(θ|D)∝p(D|θ) (6)
p (θ | D) represents the posterior probability density function of θ; p (D | theta) is a likelihood function. From the stochastic process analysis we can see that when high sampling rate and sufficient data length are met, the FFT results are approximately independent and satisfy a gaussian distribution, so ZkSubject to a zero mean Gaussian combined distribution to obtain a likelihood function of
In the formula, det (-) represents determinant; ckIs ZkThe theoretical covariance matrix of (1). Further for the convenience of analysis and calculation, a negative log-likelihood function is extracted for analysis, and expressed as L (theta), then p (D | theta) ═ exp [ -L (theta) ]]Wherein
Thus, the optimized value of the modal parameter θ can be obtained by maximizing the posterior probability density function, i.e., minimizing the negative log likelihood function L (θ).
The specific embodiment is as follows:
a. a model number machine tool complete machine is selected as a structure to be tested, a high-precision three-axis acceleration sensor with the measuring range of +/-2 g and the sensitivity of 2000 mu mv is adopted in an experiment, and a modal test system is formed by a multi-channel data acquisition card and a computer to record data. The distribution of the relevant measuring points E is shown in FIG. 3 and FIG. 4;
b. and (3) obtaining the first ten-order natural frequency of the complete machine tool and the corresponding complete machine tool vibration mode by using a modal parameter identification technology combining Bayes theory and operation modal analysis, as shown in table 1. 5-7 show three-dimensional mode diagrams corresponding to the first three-order frequency of the machine tool respectively; finally, a complete machine tool complete machine dynamic characteristic test report is formed.
TABLE 1

Claims (2)

1. A method for measuring the dynamic characteristics of the complete machine processing state of a numerical control machine tool is characterized by comprising the following specific steps:
the method comprises the following steps: a vibration accelerometer sensor is arranged on each main part of the numerical control machine tool;
step two: in a machining state, the vibration of a machine tool foundation and the vibration generated by operation are used as excitation, and the three-phase vibration response of each part of the numerical control machine tool is acquired through a multi-channel data acquisition card;
step three: obtaining frequency response function curves of all channels of the numerical control machine tool in a non-working state environment by using a Bayesian operation modal method and fast Fourier transform, and establishing a modal parameter identification function;
step four: and obtaining dynamic characteristic parameters of the numerical control machine tool, such as inherent frequency, damping ratio, modal shape and the like of each order by analyzing frequency response functions of each channel.
2. The method for measuring the dynamic characteristics of the complete machine processing state of the numerical control machine according to claim 1, wherein: the principle and the method for acquiring the three-phase vibration response of each part of the numerical control machine tool through the multi-channel data acquisition card in the machining state are as follows:
the linear vibration system of the numerical control machine tool with n degrees of freedom meets the equation of motion:
wherein M is equal to Rn×n、C∈Rn×n、K∈Rn×nRespectively a mass matrix, a damping ratio matrix and a rigidity matrix of the system; the structural response of the multi-degree-of-freedom structure satisfying linear classical damping is assumed to be:
wherein,is a global mode shape vector of the j order;the mode response is the j order mode response, the mode order is m, and the mode response satisfies the decoupling mode power equation:
wherein, ω isj,ζjAnd Wj(t) are tables, respectivelyA scalar showing the natural frequency, damping ratio and modal force of the j-th order mode of the machine. The main parameters of the mode identification are natural frequency, damping ratio and vibration mode;
the numerical control machine tool generates weak vibration under the condition of random working condition excitationN is the sampling number in each acquisition channel for the vibration acceleration response data obtained by testing; the acceleration response obtained by the hypothesis test is composed of the structural environment vibration test signal and the prediction error of the model, i.e.
Wherein, { εj∈RnThe "prediction error" is the error between the model response and the test response due to the test noise and the model error. Carrying out Fast Fourier transform (Fast Fourier transform, FFT) on the acquired machine tool acceleration signal D, namely
In the formula i2=-1;FkAnd GkRespectively representThe real and imaginary parts of (c); Δ t is the sampling step, corresponding to k being 2,3, …, Nq,fkGiven (k-1)/(N Δ t), the augmentation matrix is assumedAn amplification matrix { Z obtained by FFT of acceleration signals acquired from a machining statekStarting from the parameter, identifying the dynamic characteristic modal parameter theta of the machine tool, wherein the parameter is f, zeta, Se(ii) a Φ, where f ═ fj:j=1,…,m},ζ={ζjJ-1, …, m-and Φ - Φ1,…,Φm]∈Rn×mRespectively representing the natural frequency and the damping ratio in the m-order natural mode and the corresponding integral vibration mode of the machine tool; s is belonged to Cm×mAnd SeRespectively representing the spectral density matrix of the state excitation and the prediction error spectral density;
the method can be obtained by the Bayesian theory principle:
p(θ|D)∝p(D|θ) (6)
p (θ | D) represents the posterior probability density function of θ; p (D | theta) is a likelihood function. From the stochastic process analysis we can see that when high sampling rate and sufficient data length are met, the FFT results are approximately independent and satisfy a gaussian distribution, so ZkSubject to a zero mean Gaussian combined distribution to obtain a likelihood function of
In the formula, det (-) represents determinant; ckIs ZkFor the purpose of analysis and calculation, the theoretical covariance matrix of (a) is a negative log-likelihood function, expressed as L (θ), and p (D | θ) ═ exp [ -L (θ) for analysis]Wherein
Thus, the optimized value of the modal parameter θ can be obtained by maximizing the posterior probability density function, i.e., minimizing the negative log likelihood function L (θ).
CN201910202404.1A 2019-03-11 2019-03-11 The measurement method of numerically-controlled machine tool complete machine machining state dynamic characteristic Pending CN109884985A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910202404.1A CN109884985A (en) 2019-03-11 2019-03-11 The measurement method of numerically-controlled machine tool complete machine machining state dynamic characteristic

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910202404.1A CN109884985A (en) 2019-03-11 2019-03-11 The measurement method of numerically-controlled machine tool complete machine machining state dynamic characteristic

Publications (1)

Publication Number Publication Date
CN109884985A true CN109884985A (en) 2019-06-14

Family

ID=66932801

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910202404.1A Pending CN109884985A (en) 2019-03-11 2019-03-11 The measurement method of numerically-controlled machine tool complete machine machining state dynamic characteristic

Country Status (1)

Country Link
CN (1) CN109884985A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110991126A (en) * 2019-12-05 2020-04-10 齐鲁工业大学 Cutting machining robot dynamic stiffness modeling method based on modal analysis
CN111142375A (en) * 2019-11-25 2020-05-12 中国航空工业集团公司洛阳电光设备研究所 Mechanism natural frequency testing method for improving control stability margin
CN113370203A (en) * 2020-03-10 2021-09-10 固高科技(深圳)有限公司 Robot control method, robot control device, computer device, and storage medium
CN113901740A (en) * 2021-10-18 2022-01-07 华北电力大学 Flexible shell modal shape recognition factor and additional mass calculation method and application thereof
CN116050914A (en) * 2023-01-19 2023-05-02 上海理工大学 Quantitative evaluation method for installation consistency of machine tool equipment
TWI849858B (en) * 2023-04-21 2024-07-21 財團法人精密機械研究發展中心 Machine dynamic characteristics based intelligent controller processing parameters adjusting method and system

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005007351A1 (en) * 2003-07-15 2005-01-27 Wacker Construction Equipment Ag Working tool with damped handle
CN101718613A (en) * 2009-11-12 2010-06-02 东莞华中科技大学制造工程研究院 Experimental modal analysis method of numerical control equipment
CN101942547A (en) * 2010-07-13 2011-01-12 北京航空航天大学 Ultrasonic elliptical vibration extrusion device and vibration extrusion processing method for carrying out surface finishing of part by using same
CN103196643A (en) * 2013-03-04 2013-07-10 同济大学 Main shaft-knife handle joint surface nonlinear dynamic characteristic parameter identification method
CN103506299A (en) * 2013-10-18 2014-01-15 西南科技大学 Vibrational excitation method and device capable of automatically tracing natural frequency
CN103823406A (en) * 2014-03-11 2014-05-28 华中科技大学 Numerical control machine tool sensitive-link identification method based on modal mass distribution matrix
CN103970066A (en) * 2014-04-30 2014-08-06 华中科技大学 Numerical-control machine tool frequency-response function obtaining method based on different structure states of machine tool
CN106891306A (en) * 2017-04-25 2017-06-27 西安理工大学 Magnetic auxiliary excitation precision actuation workbench based on variation rigidity flexible structure
CN106897717A (en) * 2017-02-09 2017-06-27 同济大学 Bayesian model modification method under multiple test based on environmental excitation data
EP3261791A1 (en) * 2015-02-27 2018-01-03 Rattunde & Co GmbH Method for reducing the regenerative chatter of chip-removal machines
EP2740563B1 (en) * 2012-12-05 2018-04-04 TRUMPF Werkzeugmaschinen GmbH & Co. KG Processing device, processing machine and method for moving a machining head
CN108031870A (en) * 2017-12-04 2018-05-15 上海理工大学 A kind of main shaft of numerical control machine tool loading performance test device and test evaluation method
CN108469784A (en) * 2018-03-07 2018-08-31 上海理工大学 The measuring device and method of modal parameter suitable for numerically-controlled machine tool machining state
CN108595789A (en) * 2018-04-04 2018-09-28 上海理工大学 A kind of constrained damping structure vibration radiation acoustical power hierarchy optimization design method
CN108733899A (en) * 2018-05-02 2018-11-02 上海理工大学 The precision machine tool Dynamic performance Optimization method that frequency domain response calculates
CN109254321A (en) * 2018-07-27 2019-01-22 同济大学 Quick Bayes's Modal Parameters Identification under a kind of seismic stimulation

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005007351A1 (en) * 2003-07-15 2005-01-27 Wacker Construction Equipment Ag Working tool with damped handle
CN101718613A (en) * 2009-11-12 2010-06-02 东莞华中科技大学制造工程研究院 Experimental modal analysis method of numerical control equipment
CN101942547A (en) * 2010-07-13 2011-01-12 北京航空航天大学 Ultrasonic elliptical vibration extrusion device and vibration extrusion processing method for carrying out surface finishing of part by using same
EP2740563B1 (en) * 2012-12-05 2018-04-04 TRUMPF Werkzeugmaschinen GmbH & Co. KG Processing device, processing machine and method for moving a machining head
CN103196643A (en) * 2013-03-04 2013-07-10 同济大学 Main shaft-knife handle joint surface nonlinear dynamic characteristic parameter identification method
CN103506299A (en) * 2013-10-18 2014-01-15 西南科技大学 Vibrational excitation method and device capable of automatically tracing natural frequency
CN103823406A (en) * 2014-03-11 2014-05-28 华中科技大学 Numerical control machine tool sensitive-link identification method based on modal mass distribution matrix
CN103970066A (en) * 2014-04-30 2014-08-06 华中科技大学 Numerical-control machine tool frequency-response function obtaining method based on different structure states of machine tool
CN103970066B (en) * 2014-04-30 2017-01-11 华中科技大学 Numerical-control machine tool frequency-response function obtaining method based on different structure states of machine tool
EP3261791A1 (en) * 2015-02-27 2018-01-03 Rattunde & Co GmbH Method for reducing the regenerative chatter of chip-removal machines
CN106897717A (en) * 2017-02-09 2017-06-27 同济大学 Bayesian model modification method under multiple test based on environmental excitation data
CN106891306A (en) * 2017-04-25 2017-06-27 西安理工大学 Magnetic auxiliary excitation precision actuation workbench based on variation rigidity flexible structure
CN108031870A (en) * 2017-12-04 2018-05-15 上海理工大学 A kind of main shaft of numerical control machine tool loading performance test device and test evaluation method
CN108469784A (en) * 2018-03-07 2018-08-31 上海理工大学 The measuring device and method of modal parameter suitable for numerically-controlled machine tool machining state
CN108595789A (en) * 2018-04-04 2018-09-28 上海理工大学 A kind of constrained damping structure vibration radiation acoustical power hierarchy optimization design method
CN108733899A (en) * 2018-05-02 2018-11-02 上海理工大学 The precision machine tool Dynamic performance Optimization method that frequency domain response calculates
CN109254321A (en) * 2018-07-27 2019-01-22 同济大学 Quick Bayes's Modal Parameters Identification under a kind of seismic stimulation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YUCN K V: ""Bayesian Fast Fourier Transform Approach for Modal Updating Using Ambient Data"", 《A DVANCCS IN STRUCTURAL ENGINEERING》 *
林剑峰: ""数控机床动态特性测试与分析研究", 《机械制造》 *
胡育佳: ""基于贝叶斯运行模态分析法的砂轮架动态特性分析"", 《中国机械工程》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111142375A (en) * 2019-11-25 2020-05-12 中国航空工业集团公司洛阳电光设备研究所 Mechanism natural frequency testing method for improving control stability margin
CN110991126A (en) * 2019-12-05 2020-04-10 齐鲁工业大学 Cutting machining robot dynamic stiffness modeling method based on modal analysis
CN110991126B (en) * 2019-12-05 2023-05-26 齐鲁工业大学 Cutting machining robot dynamic stiffness modeling method based on modal analysis
CN113370203A (en) * 2020-03-10 2021-09-10 固高科技(深圳)有限公司 Robot control method, robot control device, computer device, and storage medium
CN113901740A (en) * 2021-10-18 2022-01-07 华北电力大学 Flexible shell modal shape recognition factor and additional mass calculation method and application thereof
CN116050914A (en) * 2023-01-19 2023-05-02 上海理工大学 Quantitative evaluation method for installation consistency of machine tool equipment
CN116050914B (en) * 2023-01-19 2024-07-02 上海理工大学 Quantitative evaluation method for installation consistency of machine tool equipment
TWI849858B (en) * 2023-04-21 2024-07-21 財團法人精密機械研究發展中心 Machine dynamic characteristics based intelligent controller processing parameters adjusting method and system

Similar Documents

Publication Publication Date Title
CN109884985A (en) The measurement method of numerically-controlled machine tool complete machine machining state dynamic characteristic
CN106960068B (en) Rapid modal damping ratio calculation method based on pulse excitation response frequency spectrum
CN100561162C (en) A kind of virtual oscillating table detection signal processing method and equipment thereof
El-Sheimy et al. Analysis and modeling of inertial sensors using Allan variance
CN106525226B (en) Evaluation method and system based on-site vibration load recognition
CN109902408B (en) Load identification method based on numerical operation and improved regularization algorithm
CN103823406A (en) Numerical control machine tool sensitive-link identification method based on modal mass distribution matrix
CN108629864A (en) A kind of electro spindle radial accuracy characterizing method and its system based on vibration
Guan et al. The measurement of unsteady surface pressure using a remote microphone probe
CN100498229C (en) Method for processing periodic error in inertial components
CN112733298B (en) Machining performance evaluation method of series-parallel robot at different poses based on spiral hole milling
Prato et al. A reliable sampling method to reduce large sets of measurements: a case study on the calibration of digital 3-axis MEMS accelerometers
CN112199874A (en) Modal test optimal excitation point identification method
US4031744A (en) Method and apparatus for analyzing a damped structural specimen
Herranen et al. Acceleration data acquisition and processing system for structural health monitoring
Cheng et al. AR model-based crosstalk cancellation method for operational transfer path analysis
Chang et al. Analysis of the dynamic characteristics of pressure sensors using ARX system identification
CN109029334B (en) Rock mass structural plane roughness coefficient size effect global search measurement method
Pedersen Identification techniques in composite laminates
Tavares et al. Automated damage localization for lightweight plates
Sims Dynamics diagnostics: Methods, equipment and analysis tools
CN116952354B (en) Data optimization acquisition method for driving measurement sensor
CN115495958B (en) Combined recognition method for bearing force and shafting parameters of ship propeller
CN113029570B (en) Harmonic bearing fault sample generation model and diagnosis method
Koh et al. Estimation of Milling Forces From Compliant Sensors Using a Harmonic Force Model and Kalman Filter

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190614