CN112214844A - Method and system for calculating transmission path parameters of equipment operation conditions - Google Patents

Method and system for calculating transmission path parameters of equipment operation conditions Download PDF

Info

Publication number
CN112214844A
CN112214844A CN202010851261.XA CN202010851261A CN112214844A CN 112214844 A CN112214844 A CN 112214844A CN 202010851261 A CN202010851261 A CN 202010851261A CN 112214844 A CN112214844 A CN 112214844A
Authority
CN
China
Prior art keywords
reference point
transmission path
condition
parameters
vibration response
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
Application number
CN202010851261.XA
Other languages
Chinese (zh)
Other versions
CN112214844B (en
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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202010851261.XA priority Critical patent/CN112214844B/en
Publication of CN112214844A publication Critical patent/CN112214844A/en
Application granted granted Critical
Publication of CN112214844B publication Critical patent/CN112214844B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10—Noise analysis or noise optimisation

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a method and a system for calculating parameters of a transmission path of equipment operation conditions, wherein the position of an excitation source and the transmission path are determined according to the mechanical structure characteristics of an object to be tested, then the number of reference points and the optimal installation position are determined by using a reference point optimization method, under the condition of ensuring the test working condition parameters with different linear irrelevant working conditions, acquiring the vibration response signals of the final reference point installation position and the target point installation position under the stable excitation of the test working condition parameters, and performing crosstalk elimination processing on the obtained vibration response signal, and performing least square support vector regression training on the operation condition transmission path analysis model by using the vibration response signal subjected to crosstalk elimination processing to obtain a transmission function matrix, so that the matrix ill-condition problem of the traditional operation condition transmission path analysis method is avoided, the working condition data can be completely retained, and the good generalization performance can be still kept under the test environment with complex environment.

Description

Method and system for calculating transmission path parameters of equipment operation conditions
Technical Field
The invention belongs to the field of vibration and noise reduction of mechanical equipment, and particularly relates to a method and a system for calculating parameters of a transmission path of operating conditions of equipment.
Background
The high-end equipment of the gas turbine and the like has the characteristics of numerous internal parts and complex coupling connection, and the vibration generated by the active part in the operation process is transmitted to the whole equipment through the connecting part, so that the problems of mechanical structure deformation, main part abrasion aggravation and the like are caused, the performance and the service life of the equipment are further influenced, and therefore, the vibration active control measure for the high-end equipment of the gas turbine and the like is imperative. The vibration active control strategy can be divided into vibration source control and transmission path control, wherein the vibration source control requires design optimization of a high-end equipment active component and is difficult to realize on the premise of not influencing the performance; and the transmission path control focuses on the connecting piece, and the vibration control is realized by the vibration damping device under the condition of not changing the structure of the vibration damping device, so that the operability is higher. The precondition for controlling the vibration transmission Path is to acquire the transmission Path characteristics, and the current main technical approaches include Transmission Path Analysis (TPA) and Operational Transmission Path Analysis (OTPA).
The basis of the transmission path analysis is frequency response function test, and the mechanical structure needs to be disassembled, so that time and labor are wasted; the operation condition transmission path analysis adopts the condition data to replace the frequency response function, so that the integrity of a mechanical system can be ensured in the test, and the time and the efficiency are saved. When high-end equipment such as a gas turbine and the like is used for analyzing the transmission path of the operation working condition, the number of required reference points is large due to large overall dimension and complex mechanical structure, and the serious ill-condition problem of the input matrix is easily caused. According to the traditional method, the matrix ill-condition degree is reduced through a regularization means, so that part of working condition information is omitted in the solving process, and the result cannot completely reflect the real contribution of each transmission path, so that the control precision of the transmission path is influenced, and the vibration active control precision of high-end equipment such as a gas turbine cannot be effectively improved.
Disclosure of Invention
The invention aims to provide a method and a system for calculating parameters of a transmission path of equipment operating conditions, so as to overcome the defects of the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for calculating parameters of a transmission path of equipment operation conditions comprises the following steps:
step 1), determining an excitation source position and a transmission path according to the mechanical structure characteristics of an object to be tested, determining a final reference point installation position and a final reference point installation number around the excitation source position and the transmission path by a reference point optimization method, then determining a target point installation position, and establishing an operation condition transmission path analysis model according to the reference point installation position, the reference point installation number and the target point installation position;
step 2), acquiring a vibration response signal of a final reference point installation position and a target point installation position under stable excitation of the test working condition parameters under the condition of ensuring the test working condition parameters with different working condition linearity irrelevant;
step 3), carrying out crosstalk elimination treatment on the obtained vibration response signals, and carrying out least square support vector regression training on the operation condition transmission path analysis model by using the vibration response signals subjected to crosstalk elimination treatment to obtain a transmission function matrix;
and 4) taking the vibration response signals of the final reference point mounting positions on different transmission paths under the actual working condition as test data, and multiplying the test data by the transfer function matrix obtained in the step 3) to obtain vibration energy occupation ratios of the different transmission paths, wherein the transmission path corresponding to the maximum value of the vibration energy occupation ratio is the vibration interference path.
Furthermore, a plurality of reference point candidate positions are preset for a single transmission path, the reference point signal matrix condition number is used as a judgment standard, after all reference point candidate position combinations are traversed, the reference point candidate position combination with the most complete retained information is selected as the final reference point installation position and the reference point installation number on the transmission path.
Further, the number of the test working condition parameters in the step 2) is larger than the installation number of the reference points.
Further, in step 3), the vibration response signal of the final reference point installation position after the crosstalk elimination processing is used as input data, the vibration response signal of the target point installation position after the crosstalk elimination processing is used as output data, and least square support vector regression training is performed on the operation condition transfer path analysis model to obtain an optimal linear regression function, wherein the regression function is the solved transfer function matrix.
Further, the transfer function matrix is:
Figure BDA0002644805400000031
t is a transfer function matrix, r is the number of parameters of the test working condition, alphaiIs Lagrange multiplier, xiInput data, b is the offset.
Further, specifically, the operation condition transfer path analysis model may be expressed as follows according to the relationship between the reference point signal, the target point signal and the transfer function matrix: y is XT (1),
namely, it is
Figure BDA0002644805400000032
In the formula, X is a reference point signal matrix, Y is a target point signal matrix, T is a transfer function matrix, r is the number of test working condition parameters, and n is the number of reference points; the r groups of working condition data can form a training set
Figure BDA0002644805400000033
Wherein the input data xi∈RnOutput data yie.R, the linear regression function is:
Figure BDA0002644805400000034
wherein, w is a weight vector,
Figure BDA0002644805400000035
is a high-dimensional space mapping function, and b is an offset;
calculating by adopting a secondary penalty function:
Figure BDA0002644805400000041
the constraint conditions are as follows:
Figure BDA0002644805400000042
where γ is the normalization parameter, eiIs a relaxation variable;
using lagrange multiplier alphaiAnd carrying out dual solution:
Figure BDA0002644805400000043
under the KKT condition, each variable respectively calculates partial derivatives of Lagrangian functions, and the optimal solution can be obtained as follows:
Figure BDA0002644805400000044
for ease of solution, the above equation is expressed in the form of a linear system of equations:
Figure BDA0002644805400000045
wherein the content of the first and second substances,
Figure BDA0002644805400000046
e＝[e1,e2,…,er]T，y＝[y1,y2,…,yr]T，Iis an identity matrix, alpha ═ alpha1,α2,…,αr]T，
Figure BDA0002644805400000047
Using Mercer conditions one can obtain:
Figure BDA0002644805400000048
by substituting formula (9) for formula (8) while eliminating w and e, formula (8) can be converted to:
Figure BDA0002644805400000049
solving the system of linear equations yields:
Figure BDA0002644805400000051
Figure BDA0002644805400000052
the transfer rate function matrix can be obtained by substituting equations (7) and (9) for equation (3):
Figure BDA0002644805400000053
further, under the vibration response signals obtained by r groups of test working condition parameters, r-1 groups are selected as a training set to carry out r times of parameter training, and then the average value of the r times of training results is taken as a training result.
A system for calculating parameters of a transmission path of equipment operation conditions comprises a data acquisition module, a data preprocessing module and a data processing module;
the data acquisition module is used for acquiring a vibration response signal of a final reference point installation position, a vibration response signal of a target point installation position under stable excitation of test working condition parameters and a vibration response signal of a final reference point installation position on different transmission paths under actual working conditions; the data acquisition module transmits vibration response signals of the final reference point mounting positions on different transmission paths under actual working conditions to the data processing module;
the data acquisition module transmits the acquired vibration response signal under the stable excitation of the test working condition parameters to the data preprocessing module for data preprocessing, and the preprocessed vibration response signal is transmitted to the data processing module;
the data processing module is used for storing an operation condition transmission path analysis model established according to the reference point installation positions, the reference point installation number and the target point installation positions, optimizing and training the operation condition transmission path analysis model according to the preprocessed vibration response signals to obtain a transmission function matrix, and multiplying the vibration response signals of the final reference point installation positions on different transmission paths under the actual condition by the transmission function matrix to obtain vibration energy ratio results of different transmission paths and outputting the vibration energy ratio results.
Further, the position of an excitation source and the transmission path are determined according to the mechanical structure characteristics of the object to be tested, the final installation position of reference points and the installation number of the reference points are determined around the position of the excitation source and the transmission path by a reference point optimization method, and then the installation position of the target points is determined, and an operation condition transmission path analysis model is established in the data processing module.
Compared with the prior art, the invention has the following beneficial technical effects:
the invention relates to a method for calculating the parameters of a transmission path of equipment operating conditions, which determines the position of an excitation source and the transmission path according to the mechanical structure characteristics of an object to be tested, then determines the number of reference points and the optimal installation position by using a reference point optimization method, acquires vibration response signals of the final reference point installation position and the target point installation position under the stable excitation of the parameters of the test conditions under the condition of ensuring the test conditions which are linearly uncorrelated with different conditions, carries out crosstalk elimination processing on the acquired vibration response signals, carries out least square support vector regression training on an analysis model of the transmission path of the operating conditions by using the vibration response signals after the crosstalk elimination processing to obtain a transmission function matrix, avoids the matrix ill-condition problem of the traditional analysis method of the transmission path of the operating conditions, can completely reserve the data of the operating conditions, takes the reference point response signals under, and solving contribution amounts of all paths, sequencing, and guiding high-end equipment to take vibration reduction measures by using a path contribution analysis result, so that good generalization performance can be still maintained in a test environment with a complex environment.
Furthermore, when a transfer function matrix is obtained, fitting is directly carried out by using least square support vector regression, matrix ill-condition problems caused by excessive reference points are avoided, working condition information is completely retained, random interference brought by the external environment is eliminated to the maximum extent through multi-working condition training, efficient and accurate identification of vibration transfer paths of high-end equipment such as a gas turbine is achieved, and vibration reduction measures are guided to be carried out.
Furthermore, under the vibration response signals obtained by r groups of test working condition parameters, r-1 groups are selected as a training set to carry out r times of parameter training, and then the average value of the r times of training results is taken as a training result, so that the working condition information is fully utilized, and random interference becomes the only factor influencing the model parameters.
Drawings
Fig. 1 is a schematic diagram of simulation setup in an embodiment of the present invention.
Fig. 2(a) is a result diagram of a reference point a response signal after crosstalk is removed in the embodiment of the present invention, fig. 2(B) is a result diagram of a reference point B response signal after crosstalk is removed in the embodiment of the present invention, and fig. 2(C) is a result diagram of a target point C response signal in the embodiment of the present invention.
Fig. 3 is a comparison graph of the results of the OTPA method based on the least squares support vector regression, the OTPA method based on the TSVD, and the theoretical values in the embodiment of the present invention.
Detailed Description
The invention is described in further detail below with reference to the accompanying drawings:
a method for calculating parameters of a transmission path of equipment operation conditions comprises the following steps:
step 1), determining an excitation source position and a transmission path according to the mechanical structure characteristics of an object to be tested, determining a final reference point installation position and a final reference point installation number around the excitation source position and the transmission path by a reference point optimization method, then determining a target point installation position, and establishing an operation condition transmission path analysis model according to the reference point installation position, the reference point installation number and the target point installation position;
specifically, analyzing the mechanical structure characteristics of a test object, determining the position of an excitation source according to the mechanical structure characteristics of the test object, dividing a transmission path, and selecting a reference point installation position suitable for installing a sensor around the excitation source and the transmission path; and selecting the position of a target point suitable for additionally mounting the sensor at a position to be analyzed far away from the excitation source, and according to a reference point optimization method, namely traversing all the reference point mounting position combinations by using a reference point signal matrix condition number as a judgment standard, and then taking the reference point mounting position combination with the most retained information as the selected final reference point mounting position and the reference point mounting number.
Specifically, a plurality of reference point candidate positions are preset for a single transmission path, the reference point signal matrix condition number is used as a judgment standard, and after traversing all reference point candidate position combinations, the reference point candidate position combination with the most complete retained information is selected as the final reference point installation position and the reference point installation number on the transmission path.
Step 2), acquiring a vibration response signal of a final reference point installation position and a target point installation position under stable excitation of the test working condition parameters under the condition of ensuring the test working condition parameters with different working condition linearity irrelevant;
specifically, the number of the test working condition parameters is greater than the installation number of the reference points; the vibration response signal comprises a reference point signal and a target point signal;
step 3), carrying out crosstalk elimination treatment on the obtained vibration response signals, and carrying out least square support vector regression training on the operation condition transmission path analysis model by using the vibration response signals subjected to crosstalk elimination treatment to obtain a transmission function matrix;
when the operation condition transmission path analysis model is trained, under the vibration response signals obtained by r groups of test condition parameters, r-1 groups are selected as a training set to carry out r times of parameter training, and then the average value of the r times of training results is taken as a training result, so that the condition information is fully utilized, and random interference becomes the only factor influencing the model parameters.
Specifically, crosstalk elimination processing is carried out on a final reference point installation position obtained under test working condition parameters which ensure that different working conditions are linear and irrelevant and a vibration response signal obtained at a target point installation position, then the vibration response signal of the final reference point installation position after the crosstalk elimination processing is used as input data, the vibration response signal of the target point installation position after the crosstalk elimination processing is used as output data, least square support vector regression training is carried out on an operation working condition transmission path analysis model, and an optimal linear regression function is obtained, wherein the regression function is a solved transmission function matrix; the specific process is as follows:
the OTPA model (operation condition transfer path analysis model) can be expressed as follows according to the relationship between the reference point signal, the target point signal and the transfer function matrix:
Y＝XT (1)
Figure BDA0002644805400000091
in the formula, X is a reference point signal matrix, Y is a target point signal matrix, T is a transfer function matrix, r is the number of test working condition parameters, and n is the number of reference points; the r groups of working condition data can form a training set
Figure BDA0002644805400000092
Wherein the input data xi∈RnOutput data yie.R, the linear regression function is:
Figure BDA0002644805400000093
wherein, w is a weight vector,
Figure BDA0002644805400000094
is a high-dimensional space mapping function, and b is an offset;
calculating a regression problem by adopting a secondary penalty function, and converting the regression problem into a quadratic programming problem:
Figure BDA0002644805400000095
the constraint conditions are as follows:
Figure BDA0002644805400000096
where γ is the normalization parameter, eiIs a relaxation variable, J is a penalty function, and the optimization goal is to minimize the penalty function;
the difficulty of directly solving the quadratic programming problem is high, and a Lagrange multiplier alpha is introducediAnd converting the quadratic programming problem into a dual form:
Figure BDA0002644805400000097
under the KKT condition, each variable respectively calculates partial derivatives of Lagrangian functions, and the optimal solution can be obtained as follows:
Figure BDA0002644805400000101
for ease of solution, the above equation is expressed in the form of a linear system of equations:
Figure BDA0002644805400000102
wherein the content of the first and second substances,
Figure BDA0002644805400000103
e＝[e1,e2,…,er]T，y＝[y1,y2,…,yr]Ti is an identity matrix, α ═ α1,α2,…,αr]T，
Figure BDA0002644805400000104
Using Mercer conditions one can obtain:
Figure BDA0002644805400000105
wherein Ω is a kernel function satisfying the Mercer condition;
by substituting formula (9) for formula (8) while eliminating w and e, formula (8) can be converted to:
Figure BDA0002644805400000106
solving the system of linear equations yields:
Figure BDA0002644805400000107
Figure BDA0002644805400000108
the transfer rate function matrix can be obtained by substituting equations (7) and (9) for equation (3):
Figure BDA0002644805400000109
wherein r is the number of parameters of the test working condition, alphaiIs Lagrange multiplier, xiFor the i-th input data, b is the offset, ψ is the kernel function, and a gaussian kernel is generally taken.
Step 4), specifically, selecting a working condition with serious abnormal vibration as an actual working condition, and collecting a vibration response signal of the installation position of the reference point under the actual working condition after starting up; taking vibration response signals of final reference point installation positions on different transmission paths under actual working conditions as test data, multiplying the test data by the transfer function matrix obtained in the step 3) to obtain vibration energy occupation ratios of the different transmission paths, namely contribution amounts, sequencing all the transmission paths according to the contribution amounts, wherein the transmission paths with high contribution amounts are vibration interference paths with obvious vibration, vibration reduction measures need to be taken, the vibration reduction measures can be taken from high to low according to actual conditions, in the process of the vibration reduction measures, the vibration reduction measures can be taken to influence the contribution amounts of other transmission paths after the vibration reduction measures with high contribution amounts, and the steps are repeated until the transmission path with the highest contribution amount does not need to be subjected to vibration reduction, and the vibration reduction is within a meeting requirement range.
According to the method, when the reference point is selected, an empirical method is not simply relied on, but a condition number is utilized to screen out a reference point position combination which retains most transmission path information; on the other hand, the regularization method is not adopted when the transfer function matrix is solved, so that the ill-conditioned problem involved in inverse problem solving is avoided, the working condition information is completely reserved, and the accuracy and the reliability of transfer path analysis are effectively improved.
Example 1:
1) a metal flat plate is used as a simulation object, and the material is Q235, the length is 1m, and the thickness is 5 mm. 2 excitation sources are preset A, B, 4 candidate reference points are selected around each excitation source, and a target point C is preset at a position far away from the excitation sources, and the position is shown in fig. 1. Under a given working condition, traversing all the alternative reference point combinations to obtain the reference point signal matrix condition number average value, as shown in the following table:
condition of matrix A1 A2 A3 A4
B1 128.5 127.2 127.0 129.1
B2 127.1 128.5 126.9 129.7
B3 127.2 126.5 128.7 129.4
B4 129.3 129.9 127.3 130.4
Obtaining optimal reference point combinations of A4 and B4;
2) selecting 5 groups of test condition parameters to ensure linear independence among different conditions, wherein reference point response signals after crosstalk elimination are shown in figures 2(a) and 2(b), and target point response signals are shown in figure 2(c) as shown in the following table:
vibration source A frequency (Hz) Vibration source B frequency (Hz)
Working condition 1 10 17
Working condition 2 12 15
Working condition 3 14 10
Working condition 4 16 20
Working condition 5 18 13
3) Performing crosstalk elimination processing on reference point and target point response signals under a test condition, respectively using the processed reference point and target point response signals as input data and output data, training a least square support vector regression model, and constructing an optimal linear regression function, wherein the regression function is a transfer function matrix;
4) selecting an actual working condition of a vibration source A of 20Hz and a vibration source B of 25Hz, and collecting reference point and target point signals;
5) and taking the response signals of each reference point under the actual working condition as test data, and multiplying the test data by the transfer rate function matrix obtained in the step 3) respectively to obtain the vibration energy ratio of each path and obtain the path contribution value. Comparing the contribution quantity value obtained by the invention with a theoretical value, and comparing the contribution quantity value obtained by the OTPA method based on TSVD with the theoretical value, as shown in FIG. 3, calculating the path contribution relative error, as shown in the following table:
Figure BDA0002644805400000121
the method can effectively reduce the error value of the path contribution amount, completely reserve the working condition data and improve the accuracy.

Claims (9)

1. A method for calculating parameters of a transmission path of equipment operation conditions is characterized by comprising the following steps:
step 1), determining an excitation source position and a transmission path according to the mechanical structure characteristics of an object to be tested, determining a final reference point installation position and a final reference point installation number around the excitation source position and the transmission path by a reference point optimization method, then determining a target point installation position, and establishing an operation condition transmission path analysis model according to the reference point installation position, the reference point installation number and the target point installation position;
step 2), acquiring a vibration response signal of a final reference point installation position and a target point installation position under stable excitation of the test working condition parameters under the condition of ensuring the test working condition parameters with different working condition linearity irrelevant;
step 3), carrying out crosstalk elimination treatment on the obtained vibration response signals, and carrying out least square support vector regression training on the operation condition transmission path analysis model by using the vibration response signals subjected to crosstalk elimination treatment to obtain a transmission function matrix;
and 4) taking the vibration response signals of the final reference point mounting positions on different transmission paths under the actual working condition as test data, and multiplying the test data by the transfer function matrix obtained in the step 3) to obtain vibration energy occupation ratios of the different transmission paths, wherein the transmission path corresponding to the maximum value of the vibration energy occupation ratio is the vibration interference path.
2. The method for calculating the parameters of the transmission path of the equipment operation condition according to claim 1, wherein a plurality of reference point candidate positions are preset for a single transmission path, a reference point signal matrix condition number is used as a judgment standard, and after traversing all the reference point candidate position combinations, a reference point candidate position combination with the most complete retained information is selected as a final reference point installation position and a reference point installation number on the transmission path.
3. The method for calculating the transmission path parameter of the operating condition of the equipment according to claim 1, wherein the number of the test operating condition parameters in the step 2) is larger than the installation number of the reference points.
4. The method for calculating the parameters of the transmission path of the equipment operating condition according to claim 1, wherein in the step 3), the vibration response signal of the final reference point installation position after crosstalk elimination is used as input data, the vibration response signal of the target point installation position after crosstalk elimination is used as output data, and least square support vector regression training is performed on the operating condition transmission path analysis model to obtain an optimal linear regression function, wherein the regression function is the solved transmission function matrix.
5. The method for calculating the parameters of the transmission path of the equipment operation condition according to claim 4, wherein the matrix of the transmission function is as follows:
Figure FDA0002644805390000021
t is a transfer function matrix, r is the number of parameters of the test working condition, alphaiIs Lagrange multiplier, xiFor the i-th entry, b is the offset and ψ is the kernel function.
6. The method for calculating the parameters of the transmission path of the equipment operating condition according to claim 5, wherein specifically, the relation of the analysis model of the transmission path of the operating condition according to the reference point signal, the target point signal and the transmission function matrix is expressed as: y is XT (1),
namely, it is
Figure FDA0002644805390000022
In the formula, X is a reference point signal matrix, Y is a target point signal matrix, T is a transfer function matrix, r is the number of test working condition parameters, and n is the number of reference points; r groups of working condition data can form a training set { xi,yi}1 rWherein data x is inputi∈RnOutput data yie.R, the linear regression function is:
Figure FDA0002644805390000023
wherein, w is a weight vector,
Figure FDA0002644805390000024
is a high-dimensional space mapping function, and b is an offset;
calculating by adopting a secondary penalty function:
Figure FDA0002644805390000031
the constraint conditions are as follows:
Figure FDA0002644805390000032
where γ is the normalization parameter, eiIs a relaxation variable, J is a penalty function, and the optimization goal is to minimize the penalty function;
using lagrange multiplier alphaiAnd carrying out dual solution:
Figure FDA0002644805390000033
under the KKT condition, each variable respectively calculates partial derivatives of Lagrangian functions, and the optimal solution can be obtained as follows:
Figure FDA0002644805390000034
for ease of solution, the above equation is expressed in the form of a linear system of equations:
Figure FDA0002644805390000035
wherein the content of the first and second substances,
Figure FDA0002644805390000036
e＝[e1,e2,…,er]T，y＝[y1,y2,…,yr]Ti is an identity matrix, α ═ α1,α2,…,αr]T，
Figure FDA0002644805390000037
Using Mercer conditions one can obtain:
Figure FDA0002644805390000038
wherein Ω is a kernel function satisfying the Mercer condition;
by substituting formula (9) for formula (8) while eliminating w and e, formula (8) can be converted to:
Figure FDA0002644805390000041
solving the system of linear equations yields:
Figure FDA0002644805390000042
the transfer rate function matrix can be obtained by substituting equations (7) and (9) for equation (3):
Figure FDA0002644805390000043
7. the method for calculating the parameters of the transmission path of the equipment running condition according to claim 1, wherein r-1 groups are selected as training sets to perform r times of parameter training under the vibration response signals obtained by r groups of test condition parameters, and then the average value of the r times of training results is taken as the training result.
8. A system for calculating parameters of a transmission path of equipment operation conditions is characterized by comprising a data acquisition module, a data preprocessing module and a data processing module;
the data acquisition module is used for acquiring a vibration response signal of a final reference point installation position, a vibration response signal of a target point installation position under stable excitation of test working condition parameters and a vibration response signal of a final reference point installation position on different transmission paths under actual working conditions; the data acquisition module transmits vibration response signals of the final reference point mounting positions on different transmission paths under actual working conditions to the data processing module;
the data acquisition module transmits the acquired vibration response signal under the stable excitation of the test working condition parameters to the data preprocessing module for data preprocessing, and the preprocessed vibration response signal is transmitted to the data processing module;
the data processing module is used for storing an operation condition transmission path analysis model established according to the reference point installation positions, the reference point installation number and the target point installation positions, optimizing and training the operation condition transmission path analysis model according to the preprocessed vibration response signals to obtain a transmission function matrix, and multiplying the vibration response signals of the final reference point installation positions on different transmission paths under the actual condition by the transmission function matrix to obtain vibration energy ratio results of different transmission paths and outputting the vibration energy ratio results.
9. The system of claim 8, wherein the excitation source location and the transmission path are determined according to the mechanical structure characteristics of the object to be tested, the final reference point installation location and the final reference point installation number are determined around the excitation source location and the transmission path by a reference point optimization method, and then the target point installation location is determined to establish the operation condition transmission path analysis model in the data processing module.
CN202010851261.XA 2020-08-21 2020-08-21 Method and system for calculating transmission path parameters of equipment operation conditions Active CN112214844B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010851261.XA CN112214844B (en) 2020-08-21 2020-08-21 Method and system for calculating transmission path parameters of equipment operation conditions

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010851261.XA CN112214844B (en) 2020-08-21 2020-08-21 Method and system for calculating transmission path parameters of equipment operation conditions

Publications (2)

Publication Number Publication Date
CN112214844A true CN112214844A (en) 2021-01-12
CN112214844B CN112214844B (en) 2022-12-09

Family

ID=74058701

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010851261.XA Active CN112214844B (en) 2020-08-21 2020-08-21 Method and system for calculating transmission path parameters of equipment operation conditions

Country Status (1)

Country Link
CN (1) CN112214844B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113405711A (en) * 2021-05-28 2021-09-17 上海电机学院 Motor operation condition force testing method and device

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106706118A (en) * 2016-12-09 2017-05-24 西安交通大学 Operating condition transmission path analysis method based on blind source separation technology
CN107290048A (en) * 2017-06-30 2017-10-24 西安交通大学 The operating condition Transfer Path Analysis Method of Automobile of technology is cut down based on crosstalk

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106706118A (en) * 2016-12-09 2017-05-24 西安交通大学 Operating condition transmission path analysis method based on blind source separation technology
CN107290048A (en) * 2017-06-30 2017-10-24 西安交通大学 The operating condition Transfer Path Analysis Method of Automobile of technology is cut down based on crosstalk

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
成玮等: "Tikhonov正则化在运行工况传递路径分析的应用", 《振动.测试与诊断》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113405711A (en) * 2021-05-28 2021-09-17 上海电机学院 Motor operation condition force testing method and device

Also Published As

Publication number Publication date
CN112214844B (en) 2022-12-09

Similar Documents

Publication Publication Date Title
CN107733539B (en) Testing method and system for satellite-borne multi-path microwave receiver/frequency converter
CN111638494A (en) Multi-channel amplitude and phase correction method for digital beam synthesis system
CN107807346A (en) Adaptive WKNN outdoor positionings method based on OTT Yu MR data
CN101951634B (en) Method and system for automatic power control calibration of TD-LTE terminal
CN112214844B (en) Method and system for calculating transmission path parameters of equipment operation conditions
US6839613B2 (en) Remote tuning for gas turbines
CN111478737A (en) Radio frequency test method, device, electronic device and storage medium
CN112903091A (en) Elevator guide shoe fault diagnosis method based on VMD-SVD and SSA-SVM
CN108632849A (en) A kind of method of adjustment of antenna-feed parameter, device and equipment
CN206523632U (en) A kind of noise source test system
CN109444571B (en) Moonlet communication load electromagnetic compatibility prediction method
CN117407782A (en) Signal flow graph-based transmission path analysis and vibration contribution evaluation method
CN114279554A (en) Multi-place synchronous self-adaptive performance testing method and system of low-temperature flutter sensor
CN113189513A (en) Ripple-based redundant power supply current sharing state identification method
CN109946010B (en) Device and method for testing stress of net sheet
CN115459868B (en) Millimeter wave communication performance evaluation method and system in complex environment
CN114896932A (en) Analog simulation early-stage evaluation antenna
CN115882977A (en) Spectral noise processing method and device, storage medium and electronic equipment
CN102164014A (en) Radio frequency index testing method and system
CN114598375B (en) Non-signal-level satellite anti-interference simulation system supporting interference source access
CN117728884B (en) Method, device and storage medium for detecting voltage standing wave ratio of multi-system access platform
Shen et al. Design of MultiChannel Parallel Dynamic Testing Virtual Instrument Based on Dynamic Measurement
CN110377971A (en) A kind of chip drives EQ value optimal value determines method and device
CN111310311B (en) Accurate shaping design method and system for base station antenna
CN1297146A (en) Data processing method for evaluating size and quality in making car body

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