Method for Distributing Relative Gap Parameters of Large-Scale High-Speed Rotary Equipment Components Based on Eccentricity Vector Following Measurement and Adjustment

Abstract

The present invention provides a method for distributing relative gap parameters of large-scale high-speed rotary equipment components based on eccentricity vector following measurement and adjustment. According to the present invention, a propagation process of location and orientation errors of rotors and stators of an aero-engine during assembly are analyzed, a propagation relationship of eccentricity errors after n-stage rotor and stator assembly is determined, and a coaxiality prediction model after multi-stage rotor and stator assembly is obtained; and the relative concentricity and relative runout of the rotors and stators can be further obtained by using an offset of the rotors and stators, thereby implementing the calculation of a relative gap; thereafter, a dual-objective optimization model for multi-stage rotor and stator coaxiality and relative gap amount based on an angular orientation mounting position of all stages of rotors and stators is established, the angular orientation mounting position of all stages of rotors and stators is optimized by using a genetic algorithm, to obtain an optimal mounting phase of all stages of rotors and stators; and finally, relative gap parameters of the rotor and stator can be distributed by using a probability density method.

Claims

1. A method for distributing relative gap parameters of large-scale high-speed rotary equipment components based on eccentricity vector following measurement and adjustment, wherein during multi-stage rotor and stator assembly, rotor and stator location and orientation errors are propagated and accumulated during the assembly process, an eccentricity error propagation matrix T.sub.0-n caused by location and orientation errors of all stages of rotors and stators after n-stage rotor and stator assembly being: $T_{0-n}=\left[\begin{array}{cc}\underset{i=1}{\overset{n}{.Math.}}.Math.S_{\mathrm{ri}}.Math.S_{\mathrm{xi}}.Math.S_{\mathrm{yi}}& \underset{i=1}{\overset{n}{.Math.}}.Math.\left(\underset{j=2}{\overset{i}{.Math.}}.Math.S_{\mathrm{rj}-1}.Math.S_{\mathrm{xj}-1}.Math.S_{\mathrm{yj}-1}\right).Math.S_{\mathrm{ri}}\left(p_i+{\mathrm{dp}}_i\right)\\ 0^T& 1\end{array}\right]$ where p.sub.i is an ideal position vector of a center of a radial measurement plane of an ith stage of rotor or stator, dp.sub.i is a machining error vector of a center position of the radial measurement plane of the ith stage of rotor or stator, S.sub.ri is a rotation matrix of the ith stage of rotor or stator rotating around a Z axis by an angle .sub.ri, S.sub.r1 is a unit matrix, S.sub.xi is rotation matrix of the ith stage of rotor or stator reference plane rotating around an X axis by an angle .sub.xi, S.sub.yi is a rotation matrix of the ith stage of rotor or stator reference plane rotating around a Y axis by an angle .sub.yi, S.sub.xj-1 is a rotation matrix of a (j-1)th stage of rotor or stator reference plane rotating around an X axis by an angle .sub.xj-1, S.sub.yj-1 is a rotation matrix of the (j-1)th stage of rotor or stator reference plane rotating around a Y axis by an angle .sub.yj-1, and S.sub.rj-1 is a rotation matrix of the (j-1)th stage of rotor or stator reference plane rotating around a Z axis by an angle .sub.rj-1; wherein the accumulative offset of a kth stage rotor or stator after n-stage rotor and stator assembly is expressed as: $\left[\begin{array}{c}{\mathrm{dx}}_{0-k}\\ {\mathrm{dy}}_{0-k}\end{array}\right]=\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\end{array}\right].Math.\underset{i=1}{\overset{k}{.Math.}}.Math.\left(\underset{j=2}{\overset{i}{.Math.}}.Math.S_{\mathrm{rj}-1}.Math.S_{\mathrm{xj}-1}.Math.S_{\mathrm{yj}-1}\right).Math.S_{\mathrm{ri}}\left(p_i+{\mathrm{dp}}_i\right)$ where dx.sub.0-k is the accumulative offset of a center of a measurement plane of the kth stage of rotor or stator in an X-axis direction after n-stage rotor and stator assembly, and dy.sub.0-k is the accumulative offset of the center of the measurement plane of the kth stage of rotor or stator in a Y-axis direction after n-stage rotor and stator assembly; wherein, according to an ISO standard definition of coaxiality, an expression of coaxiality after n-stage rotor and stator assembly is:
coaxiality=max{2{square root over (dx.sup.2.sub.0-k+dy.sup.2.sub.0-k)}, k=1,2, . . . , n} a coaxiality prediction model after multi-stage rotor and stator assembly is established accordingly; wherein an offset after multi-stage rotor and stator assembly is analyzed, the relative concentricity and relative runout of a rotor and a stator can be obtained by calculating offsets of the rotor and the stator, and the calculation of a relative gap after multi-stage rotor and stator assembly is implemented; wherein a dual-objective optimization model for multi-stage rotor and stator coaxiality and relative gap amount based on an angular orientation mounting position of all stages of rotors and stators is established according to a relationship between multi-stage rotor and stator coaxiality, relative concentricity, relative runout and angular orientation mounting position, the angular orientation mounting position of all stages of rotors and stators is optimized by using a genetic algorithm, so that an optimal mounting phase of all stages of rotors and stators can be obtained; and wherein an objective function of a relative gap can be obtained by using a multi-stage rotor and stator relative gap measurement model, then the probability density of the relative gap is further obtained, and then a probability relationship between contact surface runout information and relative gaps is obtained, so that relative gap parameters of a multi-stage rotor and stator are distributed.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0015] FIG. 1 is a flowchart of a relative gap parameter distribution method according to the present invention.

DETAILED DESCRIPTION OF THE INVENTION

[0016] The technical solutions in the embodiments of the present invention are clearly and completely described in the following with reference to the drawings in the embodiments of the present invention. It is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present invention without creative efforts are within the scope of protection of the present invention.

[0017] Referring to FIG. 1, the present invention provides a method for distributing relative gap parameters of large-scale high-speed rotary equipment components based on eccentricity vector following measurement and adjustment.

[0018] During multi-stage rotor and stator assembly, rotor and stator location and orientation errors are propagated and accumulated during the assembly process, an eccentricity error propagation matrix T.sub.0-n caused by location and orientation errors of all stages of rotors and stators after n-stage rotor and stator assembly being:

[00003]$T_{0-n}=\left[\begin{array}{cc}\underset{i=1}{\overset{n}{.Math.}}.Math.S_{\mathrm{ri}}.Math.S_{\mathrm{xi}}.Math.S_{\mathrm{yi}}& \underset{i=1}{\overset{n}{.Math.}}.Math.\left(\underset{j=2}{\overset{i}{.Math.}}.Math.S_{\mathrm{rj}-1}.Math.S_{\mathrm{xj}-1}.Math.S_{\mathrm{yj}-1}\right).Math.S_{\mathrm{ri}}\left(p_i+{\mathrm{dp}}_i\right)\\ 0^T& 1\end{array}\right]$

[0019] where p.sub.i is an ideal position vector of a center of a radial measurement plane of the ith stage of rotor or stator, dp.sub.i is a machining error vector of a center position of the radial measurement plane of the ith stage of rotor or stator, S.sub.ri is a rotation matrix of the ith stage of rotor or stator rotating around a Z axis by an angle .sub.ri, S.sub.r1 is a unit matrix, S.sub.xi is a rotation matrix of the ith stage of rotor or stator reference plane rotating around an X axis by an angle .sub.xi, S.sub.yi is a rotation matrix of the ith stage of rotor or stator reference plane rotating around a Y axis by an angle .sub.yi, S.sub.xj-1 is a rotation matrix of a (j-1)th stage of rotor or stator reference plane rotating around a X axis by an angle .sub.xj-1, S.sub.yj-1 is a rotation matrix of the (j-1)th stage of rotor or stator reference plane rotating around a Y axis by an angle .sub.yj-1 and S.sub.rj-1 is a rotation matrix of the (j-1)th stage of rotor or stator reference plane rotating around a Z axis by an angle .sub.rj-1.

[0020] The kth stage of rotor or stator accumulative offset after n-stage rotor and stator assembly is expressed as:

[00004]$\left[\begin{array}{c}{\mathrm{dx}}_{0-k}\\ {\mathrm{dy}}_{0-k}\end{array}\right]=\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\end{array}\right].Math.\underset{i=1}{\overset{k}{.Math.}}.Math.\left(\underset{j=2}{\overset{i}{.Math.}}.Math.S_{\mathrm{rj}-1}.Math.S_{\mathrm{xj}-1}.Math.S_{\mathrm{yj}-1}\right).Math.S_{\mathrm{ri}}\left(p_i+{\mathrm{dp}}_i\right)$

[0021] where dx.sub.0-k is the accumulative offset of a center of a measurement plane of the kth stage of rotor or stator in an X-axis direction after n-stage rotor and stator assembly, and dy.sub.0-k is the accumulative offset of the center of the measurement plane of the kth stage of rotor or stator in a Y-axis direction after n-stage rotor and stator assembly.

[0022] According to an ISO standard definition of coaxiality, an expression of coaxiality after n-stage rotor and stator assembly is:

coaxiality=max{2{square root over (dx.sup.2.sub.0-k+dy.sup.2.sub.0-k)}, k=1,2, . . . , n}

[0023] A coaxiality prediction model after multi-stage rotor and stator assembly is established accordingly.

[0024] An offset after multi-stage rotor and stator assembly is analyzed, the relative concentricity and relative runout of a rotor and a stator can be obtained by calculating offsets of the rotor and the stator, so that the calculation of a relative gap between rotor and stator after multi-stage rotor and stator assembly is implemented.

[0025] A dual-objective optimization model for multi-stage rotor and stator coaxiality and relative gap amount based on an angular orientation mounting position of all stages of rotors and stators is established according to a relationship between multi-stage rotor and stator coaxiality, relative concentricity, relative runout and angular orientation mounting position, and the angular orientation mounting position of all stages of rotors and stators is optimized by using a genetic algorithm, so that an optimal mounting phase of all stages of rotors and stators can be obtained.

[0026] An objective function of a relative gap can be obtained by using a multi-stage rotor and stator relative gap measurement model, then the probability density of the relative gap is further obtained, and then a probability relationship between contact surface runout information and relative gaps is obtained, so that relative gap parameters of a multi-stage rotor and stator can be distributed.

[0027] According to the present invention, a propagation process of location and orientation errors of a rotor and stator of an aero-engine during assembly are analyzed, a propagation relationship of eccentricity errors after n-stage rotor and stator assembly is determined, and a coaxiality prediction model after multi-stage rotor and stator assembly is obtained; the relative concentricity and relative runout of the rotor and stator can be further obtained by using an offset of the rotor and stator, so that the calculation of a relative gap is implemented; thereafter, a dual-objective optimization model for multi-stage rotor and stator coaxiality and relative gap amount based on an angular orientation mounting position of all stages of rotors and stators is established, and the angular orientation mounting position of all stages of rotors and stators is optimized by using a genetic algorithm, so that an optimal mounting phase of all stages of rotors and stators is obtained; and finally, relative gap parameters of the rotor and stator can be distributed by using a probability density method.

[0028] The above is a detailed description of the method for distributing relative gap parameters of large-scale high-speed rotary equipment components based on eccentricity vector following measurement and adjustment provided by the present invention. The principle and implementation of the present invention are described herein by using specific examples. The foregoing descriptions for the embodiments are only used to help understand the method of the present invention and its core ideas; at the same time, for a person of ordinary skill in the art, according to the idea of the present invention, there will be changes in specific implementations and application scopes. To sum up, the description is not to be construed as limiting the present invention.