Method for Optimizing Multi-Stage Components of Large-Scale High-Speed Rotary Equipment Based on Monte Carlo Bias Evaluation

Abstract

The present invention provides a method for optimizing multi-stage components of large-scale high-speed rotary equipment based on Monte Carlo bias evaluation. The method comprises: obtaining an offset of a contact surface between all stages of rotors according to a multi-stage rotor propagation relationship, and calculating coaxiality according to a coaxiality formula; calculating a cross sectional moment of inertia of the contact surface, and obtaining a bending stiffness according to a bending stiffness formula; obtaining the amount of unbalance of a rotor according to a rotor error propagation relationship; and obtaining a probability relationship between the assembly surface runout of all stages of aero-engine rotors and the final geometric concentricity, the amount of unbalance and stiffness of multi-stage rotors by using a Monte Carlo method, and optimizing the tolerance distribution and bending stiffness of the aero-engine multi-stage rotors.

Claims

1. A method for optimizing multi-stage components of a large-scale high-speed rotary equipment based on Monte Carlo bias evaluation, comprising: during n rotors assembly, single-stage rotor location and orientation errors are propagated and accumulated to affect an accumulative offset of a single-stage rotor for n rotors assembly, wherein a kth-stage rotor accumulative offset after n-stage rotor assembly may be expressed as: $\left[\begin{array}{c}d.Math.x_{0.Math..Math.k}\\ d.Math.y_{0.Math..Math.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(\stackrel{i}{\underset{j=2}{.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)$ $k=1,2,\text{...}.Math.,n,$ where dx.sub.0-k is the accumulative offset of a center of a measurement plane of a kth-stage rotor in an X-axis direction after n-stage rotor assembly, dy.sub.0-k is the accumulative offset of the center of the measurement plane of the kth-stage rotor in a Y-axis direction after n-stage rotor assembly, p.sub.i is an ideal position vector of a center of a radial measurement plane of an ith-stage rotor, dp.sub.i is a machining error vector of a center position of the radial measurement plane of the ith-stage rotor, S.sub.ri is a rotation matrix of the ith-stage rotor rotating around a Z axis for an angle .sub.ri, S.sub.r1 is unit matrix, S.sub.xj-1 is the rotation matrix of a j-1th-stage rotation stator reference plane rotating around an X axis for an angle .sub.xj-1, S.sub.yj-1 is the rotation matrix of the j-1th-stage rotation stator reference plane rotating around a Y axis for an angle .sub.yj-1, and S.sub.r j-1 is the rotation matrix of the j-1th-stage rotation stator reference plane rotating around a Z axis for an angle .sub.rj-1; according to an ISO standard definition of coaxiality, an expression of coaxiality after n-stage rotor 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 cross-sectional moment of inertia I of an inter-rotor assembly contact surface after assembly is:
I=*(R.sup.4r.sup.4)/642*.sub.0.sup.de.sub.0.sup.d*(R.sup.4r.sup.4)/64ded where R is an outer diameter of the contact surface, r is an inner diameter of the contact surface, the eccentricity is de={square root over ((dx.sub.0-k).sup.2 +(dy.sub.0-k).sup.2)}, the eccentricity angle is d=arctan(dy.sub.0-k/dx.sub.0-k) and the bending stiffness of a rotor is EI, where E is the elasticity modulus of a material, and a bending stiffness objective function is obtained; during n rotors assembly, single-stage rotation stator location and orientation errors are propagated and accumulated to affect an amount of unbalance for n rotors assembly, wherein the amount of unbalance of an nth-stage rotor caused by location and orientation errors of all stages of rotors is expressed as: $\left[\begin{array}{c}U.Math.x_{0-n}\\ U.Math.y_{0-n}\end{array}\right]=\left[\begin{array}{ccc}{m}_{0-n}& 0& 0\\ 0& m_{0-n}& 0\end{array}\right].Math.\underset{i=1}{\overset{n}{.Math.}}.Math.\left(\stackrel{i}{\underset{j=2}{.Math.}}.Math.S_{\mathrm{rj}-1}.Math.S_{\mathrm{xj}-1}.Math.S_{\mathrm{yj}-1}\right).Math.S_{\mathrm{ri}}\left(p_i.Math.+{\mathrm{dp}}_i\right)$ where Ux.sub.0-n is the amount of unbalance of a measurement plane of an assembled nth-stage rotor in an X-axis direction, Uy.sub.0-n is the amount of unbalance of the measurement plane of the assembled nth-stage rotor in a Y-axis direction, and m.sub.0-n is the mass of the assembled nth-stage rotor; performing vector addition on the amount of unbalance of a single-stage rotor and the amount of unbalance introduced by location and orientation errors during an assembly process to obtain the amount of unbalance of any stage of rotor for n rotors assembly, projecting unbalances of all stages of rotors to two correction planes respectively, combining the amount of unbalance according to a dynamic balance formula, and establishing a prediction model for the amount of unbalance of multi-stage rotors; and generating, according to a Monte Carlo method, 10,000 sets of assembly surface runout data of multi-stage rotors, bringing a random number into an objective function of multi-stage rotor coaxiality, bending stiffness and unbalance, rotating a rotation angle of each stage of aero-engine to obtain 10,000 sets of coaxiality, bending stiffness and parameters of the amount of unbalance of the multi-stage rotors, solving a probability density function according to a drawn distribution function to obtain a probability relationship between the assembly surface runout of all stages of aero-engine rotors and the final coaxiality, bending stiffness and the amount of unbalance of the multi-stage rotors, and optimizing the tolerance distribution and bending stiffness of the aero-engine multi-stage rotors.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0015] FIG. 1 is a flowchart of a method for optimizing multi-stage components of large-scale high-speed rotary equipment based on Monte Carlo bias evaluation according to the present invention.

DETAILED DESCRIPTION OF THE INVENTION

[0016] The technical solutions in the examples of the present invention are clearly and completely described in the following with reference to the drawings in the examples of the present invention. It is obvious that the described examples are only a part of the examples of the present invention, and not all of the examples. All other examples obtained by those skilled in the art based on the examples 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 optimizing multi-stage components of large-scale high-speed rotary equipment based on Monte Carlo bias evaluation.

[0018] During n rotors assembly, single-stage rotor location and orientation errors are propagated and accumulated to affect an accumulative offset of a single-stage rotor for n rotors assembly, wherein a kth-stage rotor accumulative offset after n-stage rotor assembly may be expressed as:

[00003]$\left[\begin{array}{c}d.Math.x_{0.Math..Math.k}\\ d.Math.y_{0.Math..Math.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(\stackrel{i}{\underset{j=2}{.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)$$k=1,2,\text{...}.Math.,n,$

[0019] where dx.sub.0-k is the accumulative offset of a center of a measurement plane of a kth-stage rotor in an X-axis direction after n-stage rotor assembly, dy.sub.0-k is the accumulative offset of the center of the measurement plane of the kth-stage rotor in a Y-axis direction after n-stage rotor assembly, p.sub.i is the ideal position vector of a center of a radial measurement plane of an ith-stage rotor, dp.sub.i is the machining error vector of a center position of the radial measurement plane of the ith-stage rotor, S.sub.ri is the rotation matrix of the ith-stage rotor rotating around a Z axis for an angle .sub.ri, S.sub.r1 is unit matrix, S.sub.xj-1 is the rotation matrix of a j-1th-stage rotation stator reference plane rotating around an X axis for an angle .sub.xj-1, S.sub.yj-1 is the rotation matrix of the j-1th-stage rotation stator reference plane rotating around a Y axis for an angle .sub.yj-1, and S.sub.r j-1 is the rotation matrix of the j-1th-stage rotation stator reference plane rotating around a Z axis for an angle .sub.rj-1.

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

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

[0021] A cross sectional moment of inertia I of an inter-rotor assembly contact surface after assembly is:

I=*(R.sup.4r.sup.4)/642*.sub.0.sup.de.sub.0.sup.d*(R.sup.4r.sup.4)/64ded0

[0022] where R is the outer diameter of the contact surface, r is the inner diameter of the contact surface, the eccentricity is de={square root over ((dx.sub.0-k).sup.2 (dy.sub.0-k).sup.2)}, the eccentricity angle is d=arctan(dy.sub.0-k/dx.sub.0-k), and the bending stiffness of a rotor is EI, where E is the elasticity modulus of a material, and a bending stiffness objective function is obtained.

[0023] During n rotors assembly, single-stage rotation stator location and orientation errors are propagated and accumulated to affect the amount of unbalance for n rotors assembly, wherein the amount of unbalance of an nth-stage rotor caused by location and orientation errors of all stages of rotors is expressed as:

[00004]$\left[\begin{array}{c}U.Math.x_{0-n}\\ U.Math.y_{0-n}\end{array}\right]=\left[\begin{array}{ccc}{m}_{0-n}& 0& 0\\ 0& m_{0-n}& 0\end{array}\right].Math.\underset{i=1}{\overset{n}{.Math.}}.Math.\left(\stackrel{i}{\underset{j=2}{.Math.}}.Math.S_{\mathrm{rj}-1}.Math.S_{\mathrm{xj}-1}.Math.S_{\mathrm{yj}-1}\right).Math.S_{\mathrm{ri}}\left(p_i.Math.+{\mathrm{dp}}_i\right)$

[0024] where Ux.sub.0-n is the amount of unbalance of a measurement plane of an assembled nth-stage rotor in an X-axis direction, Uy.sub.0-n is the amount of unbalance of the measurement plane of the assembled nth-stage rotor in a Y-axis direction, and m.sub.0-n is the mass of the assembled nth-stage rotor.

[0025] Vector addition is performed on the amount of unbalance of a single-stage rotor and the amount of unbalance introduced by location and orientation errors during the assembly process to obtain the amount of unbalance of any stage of rotor for n rotors assembly, unbalances of all stages of rotors are projected to two correction planes respectively, the amount of unbalance is combined according to a dynamic balance formula, and a prediction model for the amount of unbalance of the multi-stage rotors can be established.

[0026] According to a Monte Carlo method, 10,000 sets of assembly surface runout data of multi-stage rotors are generated, a random number is brought into an objective function of multi-stage rotor coaxiality, bending stiffness and the amount of unbalance, a rotation angle of each stage of aero-engine is rotated to obtain 10,000 sets of coaxiality, bending stiffness and parameters of the amount of unbalance of the multi-stage rotors, a probability density function is solved according to a drawn distribution function to obtain a probability relationship between the assembly surface runout of all stages of aero-engine rotors and the final coaxiality, bending stiffness and the amount of unbalance of the multi-stage rotors, and the tolerance distribution and bending stiffness of the aero-engine multi-stage rotors can be optimized.

[0027] The method of the present invention includes: obtaining an offset of a contact surface between all stages of rotors according to a multi-stage rotor propagation relationship, and calculating coaxiality according to a coaxiality formula; calculating a cross sectional moment of inertia of the contact surface, and obtaining a bending stiffness according to a bending stiffness formula; obtaining the amount of unbalance of a rotor according to a rotor error propagation relationship; and obtaining a probability relationship between the assembly surface runout of all stages of aero-engine rotors and the final geometric concentricity, the amount of unbalance and stiffness of multi-stage rotors by using a Monte Carlo method, and optimizing the tolerance distribution and bending stiffness of the aero-engine multi-stage rotors.

[0028] The above is a detailed description of the method for optimizing multi-stage components of large-scale high-speed rotary equipment based on Monte Carlo bias evaluation provided by the present invention. The principle and implementation manner of the present invention are described herein by using specific examples. The foregoing descriptions for the examples are only used to help understand the method of the present invention and its core ideas; at the same time, for those of ordinary skill in the art, according to the idea of the present invention, there will be changes in specific implementation manners and application scopes. To sum up, the description is not to be construed as limiting the present invention.