CN112364448A - Method for realizing linear assembly of aero-engine rotor - Google Patents
Method for realizing linear assembly of aero-engine rotor Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims description 27
- 230000009466 transformation Effects 0.000 claims description 20
- 238000005457 optimization Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 5
- 238000013519 translation Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 239000000178 monomer Substances 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000000452 restraining effect Effects 0.000 claims description 2
- 238000012800 visualization Methods 0.000 claims description 2
- 238000005452 bending Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 6
- 239000008358 core component Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/18—Manufacturability 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
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:
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:
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:
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
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:
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:
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:
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.
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)
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 |
-
2020
- 2020-10-14 CN CN202011093679.5A patent/CN112364448A/en active Pending
Patent Citations (4)
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)
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 |