CN111767671B - Unbalance parameter identification method suitable for multi-face rotor - Google Patents
Unbalance parameter identification method suitable for multi-face rotor Download PDFInfo
- Publication number
- CN111767671B CN111767671B CN202010566363.7A CN202010566363A CN111767671B CN 111767671 B CN111767671 B CN 111767671B CN 202010566363 A CN202010566363 A CN 202010566363A CN 111767671 B CN111767671 B CN 111767671B
- Authority
- CN
- China
- Prior art keywords
- particle
- chaotic
- particle swarm
- unbalance
- rotor
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Testing Of Balance (AREA)
Abstract
An unbalance parameter identification method suitable for a multi-face rotor comprises the steps of firstly establishing a dynamic model of a rotor shaft system to obtain theoretical unbalance response, and then combining actual unbalance response to construct an optimization target; then determining each particle dimension of the particle swarm, initializing the particle swarm, including randomly initializing the position and the speed of the particle swarm, and determining the chaotic inertial weight factor of each dimension according to the chaotic attenuation principle; performing particle swarm optimization, updating the particle speed and the particle position, iterating until the total iteration frequency reaches 80%, and entering a chaotic particle swarm optimization stage; determining a tth generation optimal particle, performing ten chaotic mapping processing by taking the tth generation optimal particle as an initial value to obtain a chaotic particle swarm, performing particle swarm optimization on the chaotic particle swarm to obtain the chaotic optimal particle, and randomly replacing one particle in an original particle swarm; repeating the steps until the determined iteration times are reached, wherein each parameter of the obtained optimal particles is the unbalance parameter; the invention meets the requirements of high-precision and simple-operation dynamic balance occasions of the rotor.
Description
Technical Field
The invention relates to the technical field of rotor dynamic balance, in particular to an unbalance parameter identification method suitable for a multi-face rotor.
Background
The stability of the rotor operation is an important standard for measuring the quality of the rotary mechanical performance, and in an actual field, the rotor can be unstable due to the interference of a plurality of factors, so that the dynamic unbalance of the rotor is caused. In order to ensure the normal operation of the rotor, the performance of the rotor needs to be evaluated periodically, and the influence caused by the imbalance fault is reduced. The faults that cause the unbalance of the rotor can be divided into two categories according to the occurrence reasons: static faults and dynamic faults. The static faults refer to unbalance faults still existing when the equipment stops running, such as faults of uneven rotor mass distribution, mass center deviation and the like caused by machining errors, assembly errors and the like. Dynamic faults are unbalance faults caused by damage of parts of equipment under long-term operation, such as abrasion of a rotor, bearing faults and the like. For static faults, field balancing is generally performed after the rotor is assembled, and dynamic faults need to be eliminated by adopting a real-time compensation method. For the existing large-scale rotating machinery, such as large-scale steam turbine units, generator units and the like, the rotor imbalance correction is usually completed by adopting a multi-face balancing technology.
Accurate identification of imbalance parameters is a prerequisite for dynamic balancing of the rotor. In order to accurately identify the imbalance parameters, methods have been proposed, which can be roughly divided into two types according to different research objects, one is an influence coefficient method based on experiments and a series of improved influence coefficient methods, and the other is an imbalance identification method based on models. The two methods are essentially two complementary unbalance parameter identification methods, the former method pays attention to the influence coefficient of a calibration system through multiple tests, and the method is suitable for the field of balancing of large-batch small-class rotors due to simple and convenient operation. The latter is an unbalance identification method based on a mathematical model of a research object, so that a plurality of tests of an influence coefficient method are omitted, and the method is suitable for identifying unbalance parameters of a specific rotor. However, both of the two methods have the problem of imbalance parameter calculation accuracy, and the number of the measuring point surfaces for rotor balance detection is smaller than that of the surfaces to be balanced, so that the ill-conditioned problem of the calculation equation caused by the imbalance parameter calculation is generally solved by increasing the number of tests, but the method still has certain limitation. Therefore, a method for identifying imbalance parameters of a rotor without multiple tests and with high precision is urgently needed.
The intelligence is the trend of future manufacturing development, and the intelligent connotation is spontaneity, namely, the algorithm or the equipment has the capability of continuous self discovery and self updating, and the particle swarm optimization algorithm is just the intelligent algorithm. The particle swarm optimization algorithm is an optimization algorithm simulating the population of a living being, and the optimization algorithm is an optimal algorithm by continuously performing self-adjustment in the system by linking individuals and a collective. In the problem of rotor unbalance parameter identification, the accuracy of unbalance parameter solution directly determines the quality of a method, some students indirectly realize high identification accuracy by continuously improving the accuracy of a solution model, and some students continuously improve a solution algorithm to obtain a parameter value with higher accuracy, but the method is difficult to be practically applied due to the increased calculation complexity caused by excessive pursuit of accuracy. And the particle swarm optimization algorithm is simple and convenient to calculate, and the solving precision is high, so that the method is very suitable for identifying the unbalanced parameters of the rotor.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide an unbalance parameter identification method suitable for a multi-face rotor, which can meet the requirements of high-precision and simple-operation dynamic balance occasions of the rotor.
In order to achieve the above purpose, the invention adopts the following technical scheme:
an imbalance parameter identification method suitable for a multi-face rotor comprises the following steps:
step 1, establishing a dynamic model of a rotor shaft system to obtain a theoretical unbalance response equation:
Z=Ae iΩt
wherein, A = Ω 2 [-MΩ 2 +JΩ 2 +K+iΩC] -1 Q and omega are revolution speed, M is mass matrix of the system, J is inertia matrix of the system, K is rigidity matrix of the system, C is damping matrix of the system,wherein Q ξ Distributed in xi direction in dynamic coordinate system η The unbalanced distribution of the moving coordinate system eta;
wherein n is the number of the set measuring points; meanwhile, giving an optimization range phi;
step 3, determining the dimensions of each particle of the particle swarm:
p=[x 1 ,θ 1 ,x 2 ,θ 2 ,…,x N ,θ N ]
wherein N is the number of balance surfaces of the system;
step 4, initializing a particle swarm, including randomly initializing the position and the speed of the particle swarm, and determining the chaotic inertia weight factor of each dimension according to the chaotic attenuation principle:
wherein w i t Inertial weight factor, w, for the t-th generation of particle i imax And w imin For the set maximum and minimum inertial weight factors, T is the total number of iterations, r t Attenuation radius of the t-th generation, U t Mapping sequences for the tent of the t generation;
step 6, determining the t generation optimal particle p in the chaotic particle swarm optimization stage global_best And carrying out tent chaotic mapping treatment by taking the particles as initial values to obtain chaotic particle swarms:
p=p global_best ×U
wherein, U is a ten chaotic mapping sequence;
performing particle swarm optimization on the chaotic particle swarm to obtain chaotic optimal particles p chaos_global_best Randomly replacing one particle in the original particle swarm;
step 7, repeating the step 5 and the step 6 until reaching the set iteration times to obtain the optimal particle p global_best The parameters of (1) are the imbalance parameters.
Compared with the prior art, the invention has the following beneficial technical effects:
1. according to the method, based on a dynamic model of a rotor system, after a characteristic equation of unbalance response is obtained, a particle swarm optimization solution idea is introduced to replace the original unbalance identification direct solution process, so that the ill-conditioned problem existing in the multi-surface rotor unbalance identification is effectively solved;
2. combining the chaos optimization idea and the particle swarm optimization algorithm, adopting uniformly distributed tend mapping as a mapping function, chaotizing multi-dimensional variables regarded as particles in all dimensions, and entering the chaos particle swarm optimization process when the iteration times reach 80% of the total times, thereby greatly increasing the unbalance identification precision;
3. the invention does not need complex calculation process, has simple and convenient operation and is very suitable for the actual application occasions that the measuring surface is smaller than the balance surface.
Drawings
FIG. 1 is a schematic view of a rotor system model according to an embodiment of the present invention
Fig. 2 is a schematic diagram of a rotor system measurement according to an embodiment of the present invention.
FIG. 3 is a flow chart of an embodiment of the present invention.
FIG. 4 is a diagram showing the effect of the test point S1 before and after balancing.
FIG. 5 is a diagram showing the effect of the test point S2 before and after balancing.
Detailed Description
The invention is described in detail below with reference to the figures and examples.
As shown in fig. 1, the rotor model used in the dynamic balance experimental system of the embodiment is a single-span four-balance-surface double-measurement-point rotor model, the rotor has a full length of 560mm, each 8mm shaft segment is divided into one unit, and the unit is divided into 70 units from left to right; b1 and B2 are rolling bearing supports and are positioned at a node 8 and a node 65, P1, P2, P3 and P4 are weighted balance surfaces and are positioned at a node 18, a node 30, a node 42 and a node 56, S1 and S2 are vibration measuring sensors, and unbalance response of the system is measured and is positioned at a node 24 and a node 48.
As shown in fig. 2, the experimental platform of the rotor system of the embodiment is composed of a motor speed regulator, a multi-balance-plane rotor system, a sensor, a signal pre-processor, a signal acquisition instrument, an upper computer and the like; in an experiment, a multi-balance-surface rotor system adopts a Terry RK4 experiment table, a Germany and Mimo high-precision displacement sensor is used as the sensor and is used for measuring the runout of a rotor during operation, and meanwhile, the Terry RK4 experiment table is provided with a groove used for measuring the actual rotating speed, so that the signal acquisition instrument outputs 5 paths of signals to upper computer software for signal analysis.
As shown in fig. 3, a method for identifying imbalance parameters of a multi-faceted rotor includes the following steps:
step 1, establishing a dynamic model of a rotor shaft system to obtain a theoretical unbalance response equation:
Z=Ae iΩt
wherein, A = Ω 2 [-MΩ 2 +JΩ 2 +K+iΩC] -1 Q and omega are revolution speed, M is mass matrix of the system, J is inertia matrix of the system, K is rigidity matrix of the system, C is damping matrix of the system,wherein Q ξ Is distributed in the xi direction of the dynamic coordinate system in an unbalanced way, Q η The unbalanced distribution of the moving coordinate system eta is obtained;
wherein n is the number of the set measuring points; meanwhile, giving an optimization range phi;
step 3, determining each particle dimension of the particle swarm:
p=[x 1 ,θ 1 ,x 2 ,θ 2 ,…,x N ,θ N ]
wherein N is the number of balance surfaces of the system;
step 4, initializing a particle swarm, including randomly initializing the position and the speed of the particle swarm, and determining the chaos inertial weight factor of each dimension according to the chaos attenuation principle:
wherein, w i t Inertial weight factor, w, for the t-th generation of particle i imax And w imin For the set maximum and minimum inertial weight factors, T is the total number of iterations, r t Attenuation radius of the t-th generation, U t A tent mapping sequence for the t generation;
step 6, determining the t generation optimal particle p in the chaotic particle swarm optimization stage global_best And with the particles as initial values, performing tend chaotic mapping treatment to obtain chaotic particle swarms:
p=p global_best ×U
wherein, U is a tent chaotic mapping sequence;
performing particle swarm optimization on the chaotic particle swarm to obtain chaotic optimal particles p chaos_global_best Randomly replacing one particle in the original particle swarm;
step 7, repeating the steps 5 and 6 until reaching the set iteration number to obtain the optimal particle p global_best The parameters of (1) are the imbalance parameters.
And carrying out unbalance correction according to the calculated unbalance parameters, specifically comprising the following steps: if weighting correction is adopted, selecting a proper correction screw or other correction mass blocks according to the calculated mass parameters, and then carrying out weighting correction according to the calculated phase at the balance surface; if the deduplication correction is adopted, the deduplication correction is carried out on the position where the phase is rotated by 180 degrees in the reverse direction of the phase calculated on the balance plane according to the calculated quality parameters.
The dynamic balance implementation effect of the embodiment is as follows:
as shown in FIGS. 4 and 5, the weight-increasing correction experiment is carried out at the rotation speed of the rotor system from 240r/min to 2800r/min, and it can be seen from the figure that the amplitude of the rotor system is obviously reduced after the balance, and the reduction amplitude is maximum at the rotation speed of 2040r/min, the amplitude at the point of the sensor S1 is reduced from 83.12 μm to 4.93 μm, and the reduction amplitude reaches 94.07%, and the amplitude at the point of the sensor S2 is reduced from 64.65 μm to 2.89 μm, and the reduction amplitude reaches 95.93%. The experimental result verifies the effectiveness of the unbalance parameter identification method.
The above examples are merely illustrative of the present invention and should not be construed as limiting the scope of the present invention, and any design similar or identical to the present invention is within the scope of the present invention.
Claims (1)
1. An unbalance parameter identification method suitable for a multi-face rotor is characterized by comprising the following steps:
step 1, establishing a dynamic model of a rotor shaft system to obtain a theoretical unbalance response equation:
Z=Ae iΩt
wherein, A = Ω 2 [-MΩ 2 +JΩ 2 +K+iΩC] -1 Q and omega are rotation speed, M is mass matrix of the system, J is inertia matrix of the system, K is rigidity matrix of the system, C is damping matrix of the system,wherein Q ξ Is distributed in the xi direction of the dynamic coordinate system in an unbalanced way, Q η The unbalanced distribution of the moving coordinate system eta is obtained;
step 2, obtaining theoretical unbalance response Z according to the step 1, and obtaining actual unbalance response according to actual measurementConstructing an optimization target:
wherein n is the number of the set measuring points; meanwhile, giving an optimization range phi;
step 3, determining the dimensions of each particle of the particle swarm:
p=[x 1 ,θ 1 ,x 2 ,θ 2 ,…,x N ,θ N ]
wherein N is the number of balance surfaces of the system;
step 4, initializing a particle swarm, including randomly initializing the position and the speed of the particle swarm, and determining the chaos inertial weight factor of each dimension according to the chaos attenuation principle:
wherein w i t The inertial weight factor, w, of the t-th generation of particle i imax And w imin For the set maximum and minimum inertial weight factors, T is the total number of iterations, r t Attenuation radius of the t-th generation, U t A tent mapping sequence for the t generation;
step 5, performing particle swarm optimization, updating the particle speed and the particle position, continuously iterating until the total iteration frequency reaches 80%, and entering a chaotic particle swarm optimization phase;
step 6, determining the t generation optimal particle p in the chaotic particle swarm optimization stage global_best And carrying out tent chaotic mapping treatment by taking the particles as initial values to obtain chaotic particle swarms:
p=p global_best ×U
wherein, U is a tent chaotic mapping sequence;
performing particle swarm optimization on the chaotic particle swarm to obtain chaotic optimal particles p chaos_global_best Randomly replacing one particle in the original particle swarm;
step 7, repeating the step 5 and the step 6 until reaching the set iteration times to obtain the optimal particle p global_best The parameters of (a) are the imbalance parameters.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010566363.7A CN111767671B (en) | 2020-06-19 | 2020-06-19 | Unbalance parameter identification method suitable for multi-face rotor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010566363.7A CN111767671B (en) | 2020-06-19 | 2020-06-19 | Unbalance parameter identification method suitable for multi-face rotor |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111767671A CN111767671A (en) | 2020-10-13 |
CN111767671B true CN111767671B (en) | 2022-12-09 |
Family
ID=72721144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010566363.7A Active CN111767671B (en) | 2020-06-19 | 2020-06-19 | Unbalance parameter identification method suitable for multi-face rotor |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111767671B (en) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103972908A (en) * | 2014-05-23 | 2014-08-06 | 国家电网公司 | Multi-target reactive power optimization method based on adaptive chaos particle swarm algorithm |
CN108133258B (en) * | 2016-12-01 | 2020-11-13 | 南京理工大学 | Hybrid global optimization method |
CN107240923A (en) * | 2017-08-10 | 2017-10-10 | 广东工业大学 | A kind of Method for Reactive Power Optimization in Power and device |
CN108493951A (en) * | 2018-03-21 | 2018-09-04 | 中南大学 | A kind of multi-objective reactive optimization method based on Chaos particle swarm optimization algorithm |
-
2020
- 2020-06-19 CN CN202010566363.7A patent/CN111767671B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN111767671A (en) | 2020-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101639395B (en) | Improved holographic dynamic balancing method of high-speed main shaft | |
CN109115408B (en) | Dynamic balance test method of large hydroelectric generating set based on centrifugal force equation | |
CN111767671B (en) | Unbalance parameter identification method suitable for multi-face rotor | |
CN113190786B (en) | Vibration prediction method for large-scale rotating equipment by utilizing multidimensional assembly parameters | |
CN112525533A (en) | Online detection method for contact angle of ball bearing of aero-engine | |
CN110646138A (en) | Dynamic balance method and analysis device for rotary machine without key phase and trial weight | |
Skowron et al. | Permanent magnet synchronous motor fault detection system based on transfer learning method | |
CN113357101A (en) | Method for identifying electromagnetic space-time coupling load of wind driven generator | |
CN110702314B (en) | Non-trial-weight transient high-speed dynamic balancing method based on flexible rotor acceleration response information | |
CN116839889A (en) | Method for identifying faults of rotor system and supporting system of rotary machine | |
CN111783274A (en) | Bearing fault simulation method and device | |
CN116011112A (en) | Aeroengine bearing fault data simulation method based on transmission characteristics of casing | |
CN116754134A (en) | Rotor imbalance state identification method based on test and simulation data fusion | |
CN116070090A (en) | Rolling bearing fault sample generation method based on digital twin model | |
Srinivas et al. | Role of an artificial neural network and a wavelet transform for condition monitoring of the combined faults of unbalance and cracked rotors | |
CN1103443C (en) | Holographic dynamic balance method for asymmetric rotor | |
CN109847952B (en) | Dynamic balance method of double-shaft precision centrifuge turntable based on driving current | |
Yu et al. | Unbalanced variation after assembly and double-speed influence coefficient method in the threshing drum | |
CN106951654A (en) | Complicated coupling structural dynamic characteristics parameter recognition methods based on bounded-but-unknown uncertainty | |
CN102121862B (en) | Method for measuring amount of unbalance of tire dynamic balance and tire balancer | |
Konopiński et al. | Two planes balancing method of UAV motors using a single three-axis MEMS accelerometer | |
Bastakoti et al. | Model-based initial residual unbalance identification for rotating machines in one and two planes using an iterative inverse approach | |
RU2790897C1 (en) | Method of dynamic balancing of the rotor in two planes of correction | |
Kumhar et al. | Dynamic-balance monitoring scheme for industrial fans using neural-network based machine learning approach | |
zhan et al. | Experimental Study on Unbalanced Vibration Signal of High Speed Spindle Based on All Phase Fast Fourier Transform |
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 |