Patent application title: TYPICAL ROTATIONAL PART CHARACTERIZATION METHOD BASED ON ACTUALLY MEASURED RUN-OUT DATA
Inventors:
IPC8 Class: AG01B2110FI
USPC Class:
1 1
Class name:
Publication date: 2021-08-12
Patent application number: 20210247183
Abstract:
The present invention provides a typical rotational part characterization
method based on actually measured run-out data. Aiming at the
characterization of rotational parts containing morphology data, the
present invention proposes a matrix form characterization method in which
microscopic run-out data and macroscopic axial size are comprehensively
considered. In addition, the method can be applied to an assembly
accuracy calculation process, and can characterize a single part
containing morphology feature quantities by using only one matrix M. The
calculation process of accuracy transfer is simplified, and a
high-efficiency calculation model is provided for the prediction of
assembly accuracy.Claims:
1. A typical rotational part characterization method based on actually
measured run-out data, comprising the following steps: 1) measuring a
mating face of a rotational part by a cylindricity measuring instrument
to obtain run-out data D.sub.bot of a bottom end face, radial run-out
data dR.sub.bot of a bottom spigot, run-out data D.sub.top of a top end
face, and radial run-out data dR.sub.top of a top spigot; 2) processing
the original run-out data obtained in step 1): since the data measured by
the cylindricity measuring instrument is a vector matrix with n row(s)
and 1 column, i.e., the data of each end face is an axial one-dimensional
run-out value, the data at each spigot is a radial one-dimensional
run-out value; and a corresponding method is used to process and obtain
three-dimensional coordinate data of the mating face according to
actually measured radius values r.sub.bot and r.sub.top at the circular
end faces and measured radius values R.sub.bot and R.sub.top at the
spigots in combination with actually measured run-out data; the
processing method is as follows: for the data of the bottom end face,
letting .theta. = [ 2 .times. .pi. n , 4 .times. .pi. n ,
6 .times. .pi. n , .times. , 2 .times. .pi. ] T ,
##EQU00007## then the X and Y coordinates at a bottom end face measuring
point are: X.sub.Dbot(i)=r.sub.bot.times.cos .theta..sub.(i)i=1,2 . . .
n-1,n Y.sub.Dbot(i)=r.sub.bot.times.sin .theta..sub.(i)i=1,2 . . .
n-1,n integrating the X and Y coordinates X.sub.Dbot and Y.sub.Dbot at
the bottom end face measuring point and the run-out data D.sub.bot of the
bottom end face to obtain a processed bottom end face spatial coordinate
matrix D.sub.bot', and a top end face spatial coordinate matrix
D.sub.top' can be obtained in the same way; for the radial run-out data
of the bottom spigot, according to the radial run-out data dR.sub.bot and
the measured radius value R.sub.bot, the X and Y coordinates at a bottom
spigot measuring point are:
X.sub.Rbot(i)=(R.sub.bot+dR.sub.Rbot(i)).times.cos .theta..sub.(i)
Y.sub.Rbot(i)=(R.sub.bot+dR.sub.Rbot(i)).times.cos .theta..sub.(i) due
to the spigot plays a centering role in assembly, the main concern is
about the position of a circle center, so letting
Z.sub.Rbot=0.sub.n.times.1; integrating the X, and Z coordinates
X.sub.Rbot, Y.sub.Rbot, and Z.sub.Rbot the bottom spigot run-out
measuring point to obtain the processed bottom spigot face spatial
coordinate matrix dR.sub.bot', and the top spigot face spatial coordinate
matrix dR.sub.top' can be obtained in the same way; 3) performing least
square fitting on the data obtained in step 2), and extracting the
corresponding feature quantities; the extracting method is as follows:
fitting the processed end face data D'.sub.bot and D'.sub.top by a least
square plane, and the equation of the fitted plane is: Ax+By+Cz+D=0 this
plane can be regarded as an ideal plane rotated by a certain angle around
X axis and Y axis respectively, and the corresponding deflection angles
are respectively: d .times. .theta. x = - B C ; d .times.
.theta. y = A C ##EQU00008## for a typical rotational part, four
deflection feature quantities can be extracted from the processed end
face data, which are respectively: d.theta..sub.x bot,
d.theta..sub.y_bot, d.theta..sub.x_top and d.theta..sub.y_top; fitting
the processed spigot face data dR'.sub.bot and dR'.sub.top a least square
circle, and the equation of the fitted circle is:
R.sup.2=(x-dX).sup.2+(y-dY).sup.2 for a typical rotational part, four
eccentricity feature quantities can be extracted from the processed
spigot face data, which are respectively: dX.sub.bot, dY.sub.bot,
dX.sub.top and dY.sub.top; therefore, for any rotational part with
spigots, the deflection feature quantities d.theta..sub.x bot,
d.theta..sub.x_bot, d.theta..sub.x_top and d.theta..sub.y_top of the top
and bottom end faces and the eccentricity feature quantities dX.sub.botd,
dY.sub.bot, dX.sub.top and dY.sub.top of the top and bottom spigots of
the rotational part can be obtained by performing corresponding data
processing on the actually measured run-out data of the mating faces; 4)
expressing the feature quantities of the part extracted in step 3) in a
matrix form: since most spigots adopt a connection form of short spigot
connection, and the spigot measuring point is very close to an adjacent
end face, compared with the axial height Z of the part, the axial
distance between the spigot measuring point and the adjacent end face can
be ignored; therefore, the end surface morphology feature quantity and
the spigot morphology feature quantity are coupled into a spatial
circular plane; any rotational part with spigots will include a bottom
spatial circular plane and a top spatial circular plane, and the
corresponding bottom circular plane and top circular plane are
respectively expressed as: P b .times. o .times. t = [ 1 0
d .times. .theta. y .times. .times. _ .times. .times. bot d
.times. X b .times. o .times. t 0 1 - d .times. .theta.
x .times. .times. _ .times. .times. bot d .times. Y b
.times. o .times. t - d .times. .theta. y .times. .times.
_ .times. .times. bot d .times. .theta. x .times. .times. _
.times. .times. bot 1 0 0 0 0 1 ] ##EQU00009##
P t .times. o .times. p = [ 1 0 d .times. .theta. y
.times. .times. _ .times. .times. top d .times. X top
0 1 - d .times. .theta. x .times. .times. _ .times.
.times. top d .times. Y top - d .times. .theta. y
.times. .times. _ .times. .times. top d .times. .theta. x
.times. .times. _ .times. .times. top 1 Z 0 0 0 1
] ##EQU00009.2## 5) performing spatial eccentricity and deflection
adjustment on the bottom circular plane of the part obtained in step 4):
first adjusting the circle center of the bottom circular plane to the
origin of absolute coordinates, and then adjusting the spatial deflection
amount of the bottom circular plane to 0; i.e., the bottom plane is
transformed from a spatial circular plane with a certain eccentricity
amount and deflection amount into an ideal circular plane of which the
circle center is located at the origin of absolute coordinates; the ideal
circular plane can be expressed by a fourth-order unit matrix E, and the
whole transformation process is: i . e . : P bot .times. T
.times. E T .times. P bot = E ##EQU00010## then the top
circular plane undergoes the same transformation, and the transformation
process is: i . e . : P top .times. T .times. P top '
T .times. P top = P top ' ##EQU00011## at this moment, the
bottom circular plane has been transformed into an ideal circular plane,
which no longer contains morphology feature quantities, and the
morphology of the bottom circular plane is coupled to the top circular
plane; letting M=F.sub.top', using a matrix M to characterize a
rotational part that contains microscopic morphology features and
macroscopic axial height, and using the matrix in sub sequent assembly
accuracy calculation process.Description:
TECHNICAL FIELD
[0001] The present invention relates to a characterization method for a rotational part containing interference spigots, in particular to a rotational part characterization method based on actually measured mating face run-out data.
BACKGROUND
[0002] Rotational parts are typical parts in rotating machines such as wind power equipment and engine rotors. Since the mating face of a rotational part is not an ideal plane, but a surface with certain morphology features, when assembling the rotational part, if the influence of the morphology features of the mating faces is not considered, a certain deviation will be caused between the predicted value and the true value of assembly accuracy, making it impossible to guarantee the product assembly quality, and even leading to product failure. Therefore, in an assembly prediction process, it is particularly important to realize accurate characterization of feature quantities of a part.
[0003] For a long time, scholars at home and abroad have carried out a lot of research on the characterization of parts with surface morphology. At present, the most common method is to use Small Displacement Torsor (SDT) to characterize each mating face of a part respectively. However, one part often has multiple mating faces, so it is necessary to use multiple SDT matrices to completely characterize the micro-morphology of all mating faces of the part. In addition, this characterization method can only realize the characterization of the micro-topography of the mating faces, and the macro-size of the part is ignored; whereas in an assembly process, the coupling between the macro-size of the part and the morphology of the mating faces will have a certain influence on the prediction of the assembly accuracy. Therefore, in order to realize accurate and high-efficiency prediction of the assembly accuracy, a characterization model of a typical rotational part that contains both micro-morphology features of mating faces and macroscopic key size is urgently needed.
[0004] To solve the above problems, the present invention starts with point cloud fitting technology and homogeneous coordinate transformation technology, and establishes a calculation model for characterizing a typical rotational part with spigots based on the actually measured run-out data of mating faces of the typical rotational part by using corresponding fitting methods to fit the actually measured run-out data of end faces and radial run-out data, extracting corresponding micro-morphology feature quantities, and combining with the nominal macro-size of the part itself.
SUMMARY
[0005] The purpose of the present invention is to provide a typical rotational part characterization method based on actually measured run-out data, thus to realize accurate and high-efficiency prediction of the assembly accuracy subsequently.
[0006] The technical solution of the present invention is as follows:
[0007] A typical rotational part characterization method based on actually measured run-out data, comprising the following steps:
[0008] 1) Measuring a mating face of a rotational part by a cylindricity measuring instrument to obtain run-out data D.sub.bot of a bottom end face, radial run-out data dR.sub.bot of a bottom spigot, run-out data D.sub.top of a top end face, and radial run-out data dR.sub.top of a top spigot;
[0009] 2) Processing the original run-out data obtained in step 1): since the data measured by the cylindricity measuring instrument is a vector matrix with n row(s) and 1 column, i.e., the data of each end face is an axial one-dimensional run-out value, the data at each spigot is a radial one-dimensional run-out value; and a corresponding method is used to process and obtain three-dimensional coordinate data of the mating face according to actually measured radius values r.sub.bot and r.sub.top at the circular end faces and measured radius values R.sub.bot and R.sub.top at the spigots in combination with actually measured run-out data;
[0010] The processing method is as follows:
[0011] For the data of the bottom end face, letting
.theta. = [ 2 .times. .pi. n , 4 .times. .pi. n , 6 .times. .pi. n , .times. , 2 .times. .pi. ] T , ##EQU00001##
then the X and Y coordinates at a bottom end face measuring point are:
X.sub.Dbot(i)=r.sub.bot.times.cos .theta..sub.(i),i=1,2 . . . n-1,n
Y.sub.Dbot(i)=r.sub.bot.times.sin .theta..sub.(i),i=1,2 . . . n-1,n
Integrating the X and Y coordinates X.sub.Dbot and Y.sub.Dbot at the bottom end face measuring point and the run-out data D.sub.bot of the bottom end face to obtain a processed bottom end face spatial coordinate matrix N.sub.bot', and a top end face spatial coordinate matrix D.sub.top' can be obtained in the same way.
[0012] For the radial run-out data of the bottom spigot, according to the radial run-out data dR.sub.bot and the measured radius value R.sub.bot, the X and Y coordinates at a bottom spigot measuring point are:
X.sub.Rbot(i)=(R.sub.bot+dR.sub.bot(i).times.cos .theta..sub.(i)
Y.sub.Rbot(i)=(R.sub.bot+dR.sub.bot(i).times.sin .theta..sub.(i)
[0013] Due to the spigot plays a centering role in assembly, the main concern is about the position of a circle center, so letting Z.sub.Rbot=0.sub.n.times.1; integrating the X, Y and Z coordinates X.sub.Rbot, Y.sub.Rbot and Z.sub.Rbot at the bottom spigot run-out measuring point to obtain the processed bottom spigot face spatial coordinate matrix dR.sub.bot', and the top spigot face spatial coordinate matrix dR.sub.top' can be obtained in the same way.
[0014] 3) Performing least square fitting on the data obtained in step 2), and extracting the corresponding feature quantities;
[0015] The extracting method is as follows:
[0016] Fitting the processed end face data D'.sub.bot and D'.sub.top by a least square plane, and the equation of the fitted plane is:
Ax+By+Cz+D=0
[0017] This plane can be regarded as an ideal plane rotated by a certain angle around X axis and Y axis respectively, and the corresponding deflection angles are respectively:
d .times. .theta. x = - B C ; d .times. .theta. y = A C ##EQU00002##
[0018] For a typical rotational part, four deflection feature quantities can be extracted from the processed end face data, which are respectively: d.theta..sub.x_bot, d.theta..sub.y_bot, d.theta..sub.x_top and d.theta..sub.y_top.
[0019] Fitting the processed spigot face data dR'.sub.bot and dR'.sub.top by a least square circle, and the equation of the fitted circle is:
R.sup.2=(x-dX.sub.bot).sup.2+(y-dY.sub.bot).sup.2
[0020] For a typical rotational part, four eccentricity feature quantities can be extracted from the processed spigot face data, which are respectively: dX.sub.bot, dY.sub.bot, dX.sub.top and dY.sub.top.
[0021] Therefore, for any rotational part with spigots, the deflection feature quantities d.theta..sub.x_bot, d.theta..sub.y_bot, d.theta..sub.x_top and d.theta..sub.y_top of the top and bottom end faces and the eccentricity feature quantities dX.sub.bot, dY.sub.bot, dX.sub.top and dY.sub.top of the top and bottom spigots of the rotational part can be obtained by performing corresponding data processing on the actually measured run-out data of the mating faces;
[0022] 4) Expressing the feature quantities of the part extracted in step 3) in a matrix form: since most spigots adopt a connection form of short spigot connection, and the spigot measuring point is very close to an adjacent end face, compared with the axial height Z of the part, the axial distance between the spigot measuring point and the adjacent end face can be ignored; therefore, the end surface morphology feature quantity and the spigot morphology feature quantity are coupled into a spatial circular plane; any rotational part with spigots will include a bottom spatial circular plane and a top spatial circular plane, as shown in FIG. 1, and the corresponding bottom circular plane and top circular plane are respectively expressed as:
P b .times. o .times. t = [ 1 0 d .times. .theta. y .times. .times. _ .times. .times. bot d .times. X b .times. o .times. t 0 1 - d .times. .theta. x .times. .times. _ .times. .times. bot d .times. Y b .times. o .times. t - d .times. .theta. y .times. .times. _ .times. .times. bot d .times. .theta. x .times. .times. _ .times. .times. bot 1 0 0 0 0 1 ] ##EQU00003## P t .times. o .times. p = [ 1 0 d .times. .theta. y .times. .times. _ .times. .times. top d .times. X top 0 1 - d .times. .theta. x .times. .times. _ .times. .times. top d .times. Y top - d .times. .theta. y .times. .times. _ .times. .times. top d .times. .theta. x .times. .times. _ .times. .times. top 1 Z 0 0 0 1 ] ##EQU00003.2##
[0023] 5) Performing spatial eccentricity and deflection adjustment on the bottom circular plane of the part obtained in step 4), and the whole process is shown in FIG. 2: first adjusting the circle center of the bottom circular plane to the origin of absolute coordinates, and then adjusting the spatial deflection amount of the bottom circular plane to 0; i.e., the bottom plane is transformed from a spatial circular plane with a certain eccentricity amount and deflection amount into an ideal circular plane of which the circle center is located at the origin of absolute coordinates; the ideal circular plane can be expressed by a fourth-order unit matrix E, and the whole transformation process is:
i . e . : P bot .times. T .times. E T .times. P bot = E ##EQU00004##
Then the top circular plane undergoes the same transformation, and the transformation process is:
i . e . : P top .times. T .times. P top ' T .times. P top = P top ' ##EQU00005##
[0024] At this moment, the bottom circular plane has been transformed into an ideal circular plane, which no longer contains morphology feature quantities, and the morphology of the bottom circular plane is coupled to the top circular plane;
[0025] Letting M=P.sub.top', using a matrix M to characterize a typical rotational part that contains microscopic morphology features and macroscopic axial height, and using the matrix in subsequent assembly accuracy calculation process.
[0026] The present invention has the following beneficial effects aiming at the characterization of rotational parts containing morphology data, the present invention proposes a matrix form characterization method in which microscopic run-out data and macroscopic axial size are comprehensively considered. In addition, the method can be applied to the assembly accuracy calculation process, and can characterize a single part containing morphology feature quantities by using only one matrix M. The calculation process of accuracy transfer is simplified, and a high-efficiency calculation model is provided for the prediction of assembly accuracy.
DESCRIPTION OF DRAWINGS
[0027] FIG. 1 is a schematic diagram of a spatial circular plane of a typical rotational part.
[0028] FIG. 2 shows a transformation process of a bottom circular plane of a part.
[0029] FIG. 3a shows actual measured run-out data of end faces of an aeroengine compressor rotor part.
[0030] FIG. 3b shows actual measured radial run-out data of spigots of an aeroengine compressor rotor part.
[0031] FIG. 4a is a least square fitting plane.
[0032] FIG. 4b is a least square fitting circle.
DETAILED DESCRIPTION
[0033] To make the purpose, the technical solution and the advantages of the present invention more clear, the technical solution in the present invention will be fully described below by taking a typical rotational part (a certain type of engine rotor part) as an example in combination with the drawings of the present invention.
[0034] An existing iMap4 integrated measurement and assembly platform from a company is used to measure the rotor, wherein a set of inner ring run-out data and a set of outer ring run-out data are measured for each end face, and only one set of radial run-out data is measured for each spigot. The measured data is shown in FIG. 3.
[0035] The measuring point positions in a test process are: r.sub.bot1=123, r.sub.bot2=133, r.sub.top1=168, r.sub.top2=178, R.sub.bot=120 and R.sub.top=165. The axial height of the part is Z=120. The spatial point cloud data D.sub.bot_n.times.3 and D.sub.top_n.times.3 at the end faces and the spatial point cloud data dR.sub.bot_n.times.3 and dR.sub.top_n.times.3 at the spigots can be obtained by using the method in step 2 to process the actually measured data.
[0036] Least square plane fitting is performed on the processed spatial point cloud data at the end faces, and least square circle fitting is performed on the processed spatial point cloud data at the spigots. The fitting effect is shown in FIG. 4.
[0037] The corresponding morphology feature quantities of the mating faces can be extracted by fitting, as shown in Table 1:
TABLE-US-00001 Morphology feature quantities Value d.theta..sub.x.sub.--.sub.bot (10.sup.-5 rad) 0.5129 d.theta..sub.y.sub.--.sub.bot (10.sup.-5 rad) 2.0986 d.theta..sub.x.sub.--.sub.top (10.sup.-5 rad) -1.0230 d.theta..sub.y.sub.--.sub.top (10.sup.-5 rad) 1.5391 dX.sub.bot (10.sup.-6 m) -4.1 dY.sub.bot (10.sup.-6 m) -2.7 dX.sub.top (10.sup.-6 m) -2.2 dY.sub.top (10.sup.-6 m) 3.5
[0038] The morphology feature quantities in Table 1 and the axial height Z of the part are respectively substituted into the matrices P.sub.bot and P.sub.top to obtain the bottom circular plane matrix and the top circular plane matrix of the part.
[0039] Eccentricity and deflection adjustment is performed on the bottom circular plane by the method of matrix transformation, and a transformation matrix T=P.sub.not.sup.-1 of the bottom circular plane can be obtained by the method described in step 5. Multiplying the top circular plane matrix P.sub.top by the matrix T, the transformed top circular plane matrix P.sub.top' can be obtained.
P top ' = [ 1 0 - 0 . 5 .times. 5 .times. 9 .times. 5 .times. E - 5 0 . 6 .times. E - 6 0 1 1 . 5 .times. 3 .times. 5 .times. 9 .times. E - 5 6 . 8 .times. E - 6 0 . 5 .times. 5 .times. 9 .times. 5 .times. E - 5 - 1 . 5 .times. 3 .times. 5 .times. 9 .times. E - 5 1 1 .times. 2 .times. 0 0 0 0 1 ] ##EQU00006##
[0040] Letting M=P.sub.top', the matrix M is used to characterize the engine rotor part with the morphology features of the mating faces and the macroscopic axial size taken into consideration.
User Contributions:
Comment about this patent or add new information about this topic: