Method for calculating eccentricity of rotor assembly axis based on radial runout measurement

Patent metadata

Field	Value
Publication number	US-11307015-B2
Application number	US-201816604037-A
Country	US
Kind code	B2
Filing date	Sep 12, 2018
Priority date	Sep 12, 2018
Publication date	Apr 19, 2022
Grant date	Apr 19, 2022

Abstract

A method for calculating eccentricity of rotor assembly axis based on radial runout measurement comprises matrix characterization of data and calculation of relative runout value at each point, establishment of a spring equivalent model and calculation of contact force, eccentric direction and magnitude; calculation of relative runout value; establishment of a spring equivalent model to analyze the relationship between force and displacement in each phase of a contact process, and then an uneven contact force at each point is obtained; and determination of eccentricity is to determine the magnitude of eccentricity. Based on the measured radial runout data in production practice, this method realizes the prediction of eccentricity of axis before assembly, improves the coaxiality of rotors after assembly, and has important practical guiding significance for axis prediction as well as assembly phase adjustment and optimization in the assembly process of aero-engine rotor pieces.

First claim

The invention claimed is:

1. A method for calculating eccentricity of rotor assembly axis based on radial runout measurement, wherein comprising the following steps:
step A, for two adjacent stage disks, positioning the joint of a lower stage rotor A and an upper stage rotor B by a spigot interference fit, then the centers of surfaces from the bottom up are respectively A O1 , A O2 , B O1 , B O2 ; accessing the measured radial runout values of the spigot joint surfaces of the two rotors and characterizing the values by a matrix, with the data in the form of a circle, using a polar coordinate representation method for characterization, then the radial runout data of a lower end surface of an upper end piece B is expressed as: B O1 (α,z b1 ), and the radial runout data of an upper end surface of a lower end piece A is expressed as: A O2 (α,z a2 ); α refers to the measurement angle, Z b1 and Z a2 refer to the radial runout value when the measurement angle is α; a coordinate of a reference circle center is O(0,0), and the fit spigot radius R is given, R is the radius of the reference circle;
step B, calculating the relative runout value at each point of the upper end surface and the lower end surface are in an interference fit state; the formed “burr” means that in a part manufacturing process, due to the existence of a machining error, the actual diameter of a rotor at a fit spigot is randomly generated within the range of the machining error, wherein, the radius value of each point relative to center of the reference circle is also randomly generated, and the burr dimension is the distance of each point relative to the reference circle;
step C, due to the existence of spigots interference amount, uneven elastic forces directing to the center of the reference circle are generated around both spigots of the rotor A and the rotor B when the two rotors are assembled; in order to calculate the elastic force at each measuring point, the burr dimension at each measuring point needs to be calculated; the rigidity of the burr is negligible compared with the reference circle, and the relative position of the center of the reference circle is temporarily assumed to be unchanged in the assembly process; a force caused the burr portion is equivalent to a elastic force of a spring, assuming that the burr rigidity function is k 1 (n), where n is the burr width; the reference circle rigidity is k 2 (y), where y is the distance from the surface of the reference circle; h 1 is the burr dimension of B 1 at one point, h 2 is the burr dimension of A 2 at another point; if h 2 >h 1 , x represents a displacement distance, and {circle around (1)} represents the first phase of two springs being compressed at the same time; {circle around (2)} represents the second phase of h 1 spring being compressed; {circle around (3)} represents the third phase of the reference circle compression phase; F represents the suffered elastic force; the relationship between force and displacement in each phase is as follows:

in the first phase:

F
=

k
1

⁡

(
n
)

⁢

x
2

(
1
)

in the second phase:
F=k 1 ( n ) h 2 +k 1 ( n )( x− 2 h 2 )  (2)

in the third phase:
F=k 2 ( y )( x−h 1 −h 2 )+ k 1 ( n ) h 2 +k 1 ( n )( h 1 −h 2 )  (3)

if h 2 <h 1 , it is obtained in a similar way that:
in the first phase:

F
=

k
1

⁡

(
n
)

⁢

x
2

(
4
)

in the second phase:
F=k 1 ( n ) h 1 +k 1 ( n )( x− 2 h 1 )  (5)

in the third phase:
F=k 2 ( y )( x−h 1 −h 2 )+ k 1 ( n ) h 1 +k 1 ( n )( h 2 −h 1 )  (6)

step D, calculating the contact force at each point in an ideal state, includes: considering the contact surfaces of spigot of the rotor A and spigot of the rotor B as springs, the uneven radial runout values of the two spigots are equivalent to uneven elastic forces, and the solid portion of a base body is regarded as a rigid body without deformation; calculating the elastic force at each point in the case where the centers of the two joint surfaces coincides with each other in the ideal state first according to formulas (1) to (6);
the formulas are the same as formulas (1) to (6);
step E, calculating resultant force vector, includes: combining the spring force of each measuring point at the center of the reference circle, and obtaining a resultant force F n ;
F n =Σ 1 n F i   (7)

where, i represents the number of measuring points;
step F, calculating eccentricity e: the offset direction of the actual centroid relative to the center of the reference circle, wherein, the direction of eccentricity is the direction of the resultant force F n ; when the relative positions of the centers of the reference circle of the two rotors are moved in F n direction, the contact portions of each pair of measuring points are changed, and the elastic force of each measuring point is changed until the force is balanced and an equilibrium state is achieved; the centroid offset amount and eccentric angle are calculated according to this principle;
since the center of the reference circle is moved by a certain distance, the displacement distance of each point is different; e is the eccentricity, O is the center of the reference circle, O 2 is the center of an eccentric circle, r is the radius of the reference circle of the rotors, d is the displacement of a measuring point, θ is the eccentric angle, θ 2 is the angle corresponding to the measuring point, and γ is the angle corresponding to a triangle;
it is obtained according to the geometrical relationship that:

{

r
2

=

e
2

+

Assignees

• Univ Dalian Tech

Inventors

• Sun Qingchao
• Liu Xin
• Gao Yichao
• Wang Yunlong

Classifications

• G01B5/201 Primary

  for measuring roundness · CPC title

• G06F17/16

  Matrix or vector computation {, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization (matrix transposition G06F7/78)} · CPC title

• G01B5/252 Primary

  for measuring eccentricity, i.e. lateral shift between two parallel axes · CPC title