CN112364448A - Method for realizing linear assembly of aero-engine rotor - Google Patents

Method for realizing linear assembly of aero-engine rotor Download PDF

Info

Publication number
CN112364448A
CN112364448A CN202011093679.5A CN202011093679A CN112364448A CN 112364448 A CN112364448 A CN 112364448A CN 202011093679 A CN202011093679 A CN 202011093679A CN 112364448 A CN112364448 A CN 112364448A
Authority
CN
China
Prior art keywords
rotor
assembly
stage
transformation matrix
stage 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.)
Pending
Application number
CN202011093679.5A
Other languages
Chinese (zh)
Inventor
孙清超
李明华
汪云龙
刘亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202011093679.5A priority Critical patent/CN112364448A/en
Publication of CN112364448A publication Critical patent/CN112364448A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/18Manufacturability analysis or optimisation for manufacturability

Abstract

The invention provides a method for realizing linear assembly of an aircraft engine rotor, and belongs to the technical field of precision assembly. The invention provides a feasible method aiming at optimizing rotor assembly phase as a target aiming at the 'bending' phenomenon after rotor assembly, which takes the actual rotation axis of a multi-stage rotor as a reference and adjusts the assembly phase of each stage of rotor step by step so as to ensure that the center of each stage of rotor is on the actual rotation axis as much as possible, thereby fundamentally lightening the 'bending' phenomenon after rotor assembly and improving the rotor assembly quality and the one-time assembly success rate.

Description

Method for realizing linear assembly of aero-engine rotor
Technical Field
The invention relates to a method for realizing linear assembly of an aircraft engine rotor, and belongs to the technical field of precision assembly.
Background
The core technology of the aircraft engine, as the bright bead on the industrial crown, is a typical high-precision technology product, is monopolized by a few developed countries for a long time, and strict technical blockade is implemented in China, so that a large amount of engineering experience and long-term technical accumulation far lead China. The rotor is taken as a core component of an aircraft engine, the assembly technology of the rotor is always the key point of the research field of the aircraft engine, although a large number of published documents simulating the rotor assembly technology can be found, due to the reasons of core technology blocking protection, military confidentiality, commercial confidentiality and the like, an engineering case and related result data aiming at the assembly technology of a real rotor system of the aircraft engine are almost blank, and the technical barrier is difficult to break.
The rotors are generally manufactured in a grading way and assembled in a stacking way, because the current domestic rotor assembly also depends on experience and trial and error methods in a large quantity, the offset of the centers of all levels of rotors relative to the actual rotation axis is often overlarge, and then a bent rotor is generated, and the bent rotor is easy to generate faults of vibration, friction and the like in a high-speed rotation state, so that the performance of an aeroengine is seriously influenced. Therefore, the assembly quality of the rotor is not high at present and only meets the requirement basically. With the increasing overall performance of the aircraft engine, the requirement on the core component rotor is also increased at a very high speed, and the assembly precision of the rotor is difficult to improve by purely depending on manual experience and trial and error.
Aiming at the problem of the bent rotor, the invention provides a linear assembly method based on the actual rotary axis, which seeks the optimal assembly phase between each stage of rotor through an assembly joint surface rotation matrix, further minimizes the distance from the center of each stage of rotor to the actual rotary axis, realizes the real linear assembly and improves the one-time assembly success rate of the rotor.
Disclosure of Invention
The method for realizing linear assembly of the rotor of the aircraft engine is characterized in that an actual rotation axis which passes through the center of the bottom surface of a bottom-stage rotor and is perpendicular to the bottom surface is used as a measurement standard of coaxiality during assembly, coordinate changes of the centers of all stages of rotors are expressed by means of a transformation matrix of assembly characteristics of two adjacent stages of rotors, the centers of all stages of rotors are enabled to be on the actual rotation axis as much as possible by the condition of minimizing the coaxiality, linear assembly of the rotors is really realized, and the success rate of primary assembly of the rotors is improved.
The technical scheme of the invention is as follows:
a method for realizing linear assembly of an aircraft engine rotor comprises the following steps:
step (1), a precision rotary table system is used as a basis, a monomer transformation matrix of each stage of rotor assembly characteristics is established through test data, the assembly process of two adjacent stages of rotors is represented by a rotation matrix of assembly joint surface characteristics, and a comprehensive transformation matrix from a bottom stage rotor to a top stage rotor is established by a rotor stacking theory;
four groups of data, namely the end face runout and the radial runout of the upper assembling face and the lower assembling face are measured according to the assembling characteristics of the single-stage rotor, and a transformation matrix from the lower assembling face to the upper assembling face of the single-stage rotor is established by using the groups of data as follows:
Figure BDA0002722958100000021
due to the existence of the processing error of the rotor, improper adjustment in the assembly test process and other factors, the joint end face of the spigot can incline when the two-stage rotor is assembled, the good centering and restraining effect of the spigot is considered, the assembly optimization of the two-stage rotor can be realized only by rotating the rotor around the z axis, and the assembly joint face rotating matrix T is introduced1-1’As follows:
Figure BDA0002722958100000022
thus, the transformation matrix T of the assembly characteristics 0-2 at the time of assembly of the two-stage rotor0-2Write as:
T0-2=T0-1*T1-1’*T1‘-2
similarly, when the n-level rotor is assembled, the comprehensive transformation matrix is as follows:
T0-n=T0-1*T1-1’*T1‘-2*……*Tm-m’*Tm’-n
wherein m is n-1;
T0-nthe structure of (A) is as follows:
Figure BDA0002722958100000031
wherein R isxyzRepresenting rotation of the bottom-stage rotor to the top-stage rotor about x, y, z, pxyz=[dxn dyn zn]TRepresenting the translation of the bottom-stage rotor to the top-stage rotor in the x, y, z directions;
step (2), an optimization method for seeking an optimal assembly phase by minimizing coaxiality is provided; the translation vector p in the comprehensive transformation matrix during the assembly of the n-level rotorxyzCenter O characterizing Assembly feature 00How to change to the center O of the assembly feature nnWherein the sum of the x and y components represents OnDistance to the actual axis of revolution; the coaxiality is defined as: 2 [ (d)xn)2+(dyn)2]1/2A function f (θ) as followszi) From the above analysis, f (θ) is known as an optimization targetzi) About thetaziIn min { f (θ) }zi) Solving to obtain the optimal phase when assembling the multistage rotor as an optimization target;
f(θzi)=[(dxn)2+(dyn)2];n=2,3,4……;i=1,2,3……n-1;
take two-stage rotor assembly as an example, i.e. when n is 2 and i is 1, θz1The direction of the adjustment phase of the lower-level rotor relative to the upper-level rotor is determined by the right-hand rule, and if the upper-level rotor needs to be adjusted during assembly, the rotating direction is the reverse direction; on the basis, 3 and 4 stages are optimized to n stages.
Step (3), rotor assembly visualization and comparison effect
And drawing an equivalent schematic diagram of the assembled rotor to check the actual pose of the rotor and compare the actual pose with the assembling effect during random assembly.
The invention has the beneficial effects that: the invention provides a feasible method aiming at optimizing rotor assembly phase as a target aiming at the 'bending' phenomenon after rotor assembly, which takes the actual rotation axis of a multi-stage rotor as a reference and adjusts the assembly phase of each stage of rotor step by step so as to ensure that the center of each stage of rotor is on the actual rotation axis as much as possible, thereby fundamentally lightening the 'bending' phenomenon after rotor assembly and improving the rotor assembly quality and the one-time assembly success rate.
Drawings
FIG. 1 is a flow chart of the present invention.
Fig. 2 is a single stage rotor transformation matrix definition.
Fig. 3a is a schematic diagram of random assembly of a two-stage rotor.
Fig. 3b is a schematic diagram of the optimized assembly of the two-stage rotor.
Fig. 4 shows the coaxiality contrast in the two assembly modes.
Detailed Description
The invention is further described below with reference to the accompanying drawings, taking a 3-stage rotor assembly as an example.
Step one, single-stage rotor transformation matrix definition
And fitting the radial run-out and end run-out data based on a least square method By considering the eccentric error of rotor installation to obtain a fitting circle center (dx, dy) and a fitting plane Ax + By + Cz + D which is 0. (dx, dy) is an eccentric coordinate, d θx-B/C is the angle of rotation of the fitting plane about the x-axis relative to the reference plane, d θyAnd A/C is the rotation angle of the fitting plane relative to the reference plane around the y axis, and the rotation direction follows the right-hand spiral rule. Wherein, OEIs the origin of a reference coordinate system, O0Is the center of the lower bottom surface of the rotor, O1The center of the upper bottom surface of the rotor, E represents a reference surface, 0 represents a lower bottom surface, and 1 represents an upper bottom surface, as shown in FIG. 2. The transformation matrix from E to 0 is TE-0Transformation matrix from E to 1 as TE-1So that the transformation matrix T is from 0 to 10-1=T-1 E-0*TE-1. The same method obtains the monomer transformation matrix of 2 and 3-stage rotors, as shown in table 1.
TABLE 1
Figure BDA0002722958100000051
Step two, establishing a stacking model
An optimized schematic diagram of a 2-stage rotor assembly is established considering the inclination of the joint surface of the two-stage rotor assembly, as shown in fig. 3. And simultaneously establishing a comprehensive change matrix during the assembly of the 2-stage and 3-stage rotors, thus obtaining the translation change of the center of the bottom-stage rotor to the centers of the 2 nd-stage and 3 rd-stage rotors about the x and y directions, namely f (theta)zi). Using matlab to solve f (theta)z1),f(θz2) Minimum assembly phase θz1=230°,θz2=215°。
And step three, establishing coaxiality contrast of optimized assembly and random assembly of the three-stage rotor.
Substituting the assembly phase solved in the second step into a comprehensive transformation matrix to obtain the coaxiality of the 2-level and 3-level rotor during optimization assembly; substituting a group of random assembly phases into the comprehensive transformation matrix to obtain the coaxiality of 2 and 3-level rotors during random assembly; line graphs were plotted for comparison, as in FIG. 4. Analysis shows that the coaxiality of the 2-level rotor and the 3-level rotor after optimized assembly is obviously reduced compared with random assembly, namely the rotors well realize 'linear assembly', and the assembly precision and quality are improved.

Claims (1)

1. A method for realizing linear assembly of an aircraft engine rotor is characterized by comprising the following steps:
step (1), a precision rotary table system is used as a basis, a monomer transformation matrix of each stage of rotor assembly characteristics is established through test data, the assembly process of two adjacent stages of rotors is represented by a rotation matrix of assembly joint surface characteristics, and a comprehensive transformation matrix from a bottom stage rotor to a top stage rotor is established by a rotor stacking theory;
four groups of data, namely the end face runout and the radial runout of the upper assembling face and the lower assembling face are measured according to the assembling characteristics of the single-stage rotor, and a transformation matrix from the lower assembling face to the upper assembling face of the single-stage rotor is established by using the groups of data as follows:
Figure FDA0002722958090000011
when the two-stage rotor is assembled, the joint end face of the spigot is inclined, the good centering and restraining effects of the spigot are considered, the assembly optimization of the two-stage rotor can be realized only by rotating the rotor around the z axis, and therefore the rotating matrix T of the joint surface is introduced1-1’As follows:
Figure FDA0002722958090000012
thus, the transformation matrix T of the assembly characteristics 0-2 at the time of assembly of the two-stage rotor0-2Write as:
T0-2=T0-1*T1-1’*T1‘-2
similarly, when the n-level rotor is assembled, the comprehensive transformation matrix is as follows:
T0-n=T0-1*T1-1’*T1‘-2*……*Tm-m’*Tm’-n
wherein m is n-1;
T0-nthe structure of (A) is as follows:
Figure FDA0002722958090000013
wherein R isxyzRepresenting rotation of the bottom-stage rotor to the top-stage rotor about x, y, z, pxyz=[dxn dyn zn]TRepresenting the translation of the bottom-stage rotor to the top-stage rotor in the x, y, z directions;
step (2), an optimization method for seeking an optimal assembly phase by minimizing coaxiality is provided; the translation vector p in the comprehensive transformation matrix during the assembly of the n-level rotorxyzCenter O characterizing Assembly feature 00How to change to the center O of the assembly feature nnWherein the sum of the x and y components represents OnDistance to the actual axis of revolution; the coaxiality is defined as: 2 [ (d)xn)2+(dyn)2]1/2A function f (θ) as followszi) From the above analysis, f (θ) is known as an optimization targetzi) About thetaziIn min { f (θ) }zi) Solving to obtain the optimal phase when assembling the multistage rotor as an optimization target;
f(θzi)=[(dxn)2+(dyn)2];n=2,3,4……;i=1,2,3……n-1;
take two-stage rotor assembly as an example, i.e. when n is 2 and i is 1, θz1The direction of the adjustment phase of the lower-level rotor relative to the upper-level rotor is determined by the right-hand rule, and if the upper-level rotor needs to be adjusted during assembly, the rotating direction is the reverse direction; on the basis, optimizing 3 and 4 to n levels;
step (3), rotor assembly visualization and comparison effect
And drawing an equivalent schematic diagram of the assembled rotor to check the actual pose of the rotor and compare the actual pose with the assembling effect during random assembly.
CN202011093679.5A 2020-10-14 2020-10-14 Method for realizing linear assembly of aero-engine rotor Pending CN112364448A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011093679.5A CN112364448A (en) 2020-10-14 2020-10-14 Method for realizing linear assembly of aero-engine rotor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011093679.5A CN112364448A (en) 2020-10-14 2020-10-14 Method for realizing linear assembly of aero-engine rotor

Publications (1)

Publication Number Publication Date
CN112364448A true CN112364448A (en) 2021-02-12

Family

ID=74507165

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011093679.5A Pending CN112364448A (en) 2020-10-14 2020-10-14 Method for realizing linear assembly of aero-engine rotor

Country Status (1)

Country Link
CN (1) CN112364448A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109117460A (en) * 2018-09-12 2019-01-01 大连理工大学 A method of survey calculation rotor is jumped based on end and assembles axis deviation
CN110119553A (en) * 2019-04-29 2019-08-13 西北工业大学 A kind of aeroengine rotor parts selection matching optimization method of seam allowance connection
US20200217211A1 (en) * 2019-01-07 2020-07-09 Harbin Institute Of Technology Method for Optimizing Multi-Stage Components of Large-Scale High-Speed Rotary Equipment Based on Monte Carlo Bias Evaluation
US20200309502A1 (en) * 2018-09-12 2020-10-01 Dalian University Of Technology Method for calculating eccentricity of rotor assembly axis based on radial runout measurement

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109117460A (en) * 2018-09-12 2019-01-01 大连理工大学 A method of survey calculation rotor is jumped based on end and assembles axis deviation
US20200309502A1 (en) * 2018-09-12 2020-10-01 Dalian University Of Technology Method for calculating eccentricity of rotor assembly axis based on radial runout measurement
US20200217211A1 (en) * 2019-01-07 2020-07-09 Harbin Institute Of Technology Method for Optimizing Multi-Stage Components of Large-Scale High-Speed Rotary Equipment Based on Monte Carlo Bias Evaluation
CN110119553A (en) * 2019-04-29 2019-08-13 西北工业大学 A kind of aeroengine rotor parts selection matching optimization method of seam allowance connection

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孙汕民;李永;: "航空发动机转子安装边制造误差对装配性能的影响研究", 现代制造技术与装备, no. 07, 15 July 2020 (2020-07-15), pages 126 - 131 *

Similar Documents

Publication Publication Date Title
CN110595690B (en) Large-scale high-speed rotation equipment measurement and intelligent learning assembly method and device based on centroid, gravity center and inertia center vector minimization
CN105426566B (en) Large high-speed revolution equipment form and position tolerance distribution method based on Monte Carlo technique
CN110598229B (en) Large-scale high-speed rotation equipment multi-stage part optimization method based on Monte Carlo deviation evaluation
CN111460677B (en) Method for establishing rotor stacking precision prediction model based on geometric algebraic theory
CN109960867B (en) Large-scale high-speed rotation equipment multi-stage part matching method based on inertial rotation center regulation optimization and intelligent learning
CN110119553B (en) Matching optimization method for rotor parts of aero-engine connected by spigot
CN108694294B (en) Method for predicting assembly eccentricity of high-pressure rotor of aircraft engine
CN105426565A (en) Probability density technology based geometrical tolerance distribution method for large high-speed rotary equipment
CN111384828A (en) Method for optimizing rotor assembly
CN112364448A (en) Method for realizing linear assembly of aero-engine rotor
CN110877751B (en) Large-scale high-speed rotation equipment rotational inertia stacking method based on vector projection
CN101972935A (en) Processing method and auxiliary processing device thereof of gear box body
CN1812347A (en) Protocol validity verifying and testing method based on mode conversion
CN110717293B (en) Rotor spigot bolt connection combination surface deformation rule fitting method
CN111076867A (en) Large-scale high-speed rotation equipment multistage part unbalance amount distribution method based on synchronous measurement and adjustment of mass center and inertia center
CN110595689B (en) Large-scale high-speed rotation equipment multistage part unbalance amount prediction method and part assembly device
CN115828446B (en) Method for improving rotor stacking optimization quality and calculation speed
CN116341038A (en) Method, system, equipment and storage medium for machining crankshaft center hole
CN116541980A (en) Rotor punching size parameter optimization method
CN107766647B (en) Numerical calculation method for profile of formed milling cutter for machining elliptical spiral rotor
CN111475890A (en) Method for distributing relative clearance parameters of parts of large-scale high-speed rotation equipment based on eccentric vector following measurement and adjustment
CN110929397B (en) Low-rotation-speed characteristic line extrapolation method of gas turbine based on zero-rotation-speed characteristic line
CN111666639B (en) Optimization method for eliminating assembly of crankshaft type rotor
CN102539071A (en) Fan rotor on-line dynamic balance method based on force parallel decomposition principle
CN105719310A (en) Collision detecting method and device

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