CN111076656A - Part tolerance distribution method and device based on four-parameter compensation - Google Patents
Part tolerance distribution method and device based on four-parameter compensation Download PDFInfo
- Publication number
- CN111076656A CN111076656A CN201811231186.6A CN201811231186A CN111076656A CN 111076656 A CN111076656 A CN 111076656A CN 201811231186 A CN201811231186 A CN 201811231186A CN 111076656 A CN111076656 A CN 111076656A
- Authority
- CN
- China
- Prior art keywords
- static balance
- measuring
- parameter
- measuring table
- sensor
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/30—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring angles or tapers; for testing the alignment of axes
- G01B7/31—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring angles or tapers; for testing the alignment of axes for testing the alignment of axes
- G01B7/312—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring angles or tapers; for testing the alignment of axes for testing the alignment of axes for measuring eccentricity, i.e. lateral shift between two parallel axes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Engines (AREA)
Abstract
The invention provides a part tolerance distribution method and a part tolerance distribution device based on four-parameter compensation.A four-parameter circular profile measurement model is established by considering four parameter components of a rotor eccentric error, a sensor measuring head offset, a sensor measuring ball radius and an included angle between a rotary axis and a geometric axis in circular profile measurement; and generating random numbers according to a Monte Carlo method, drawing a distribution function to obtain a probability density function, obtaining a probability relation between contact surface jumping and eccentric errors, and realizing single-stage rotor tolerance distribution.
Description
Technical Field
The invention belongs to the technical field of mechanical tolerance distribution, and particularly relates to a part tolerance distribution method and device based on four-parameter compensation.
Background
The aero-engine is used as a pearl on the crown of the human industry, and intensively embodies the national science and technology, industry and national defense strength. With the vigorous development of aviation industry in China, the requirements on the performance of the aero-engine are continuously improved. When the development of an aero-engine with a higher thrust-weight ratio is pursued, the vibration of the engine is also emphasized to be reduced so as to improve the performance of the engine and ensure the flight safety.
Reducing engine vibration is an important goal in engine development. The assembly quality of the rotor has a great influence on the performance of the aircraft engine, and therefore in order to reduce vibration, reasonable assignment needs to be carried out on the design of the common distribution of the rotor of the aircraft engine. Therefore, the machining precision of the engine rotor can be greatly improved, and the machining cost is reduced.
The current common method is to firstly select an effective part assembling sequence; secondly, reasonably selecting positioning parts for machining the matching parts, determining the number and the size of the matching parts, and positioning and fastening the parts relative to each other to form the assembling part; finally, the statistical dimensional tolerance of each individual component part is determined as the manufacturing requirement of the component part, and the tolerance is used for reducing the processing cost and enabling the assembly part to meet the dimensional tolerance requirement. In the prior art, the transmission relation of the circle center coordinates of n-level equipment after assembly can be determined by analyzing the positioning of the radial and axial measuring surfaces of the rotary equipment and the transmission process of the orientation tolerance in the assembly, so as to obtain the relation between the eccentricity of the assembled equipment and the positioning, orientation tolerance and rotation angle of each level of equipment; and obtaining the probability density of the coaxiality tolerance of the n-level equipment according to the target function of the coaxiality tolerance, and finally obtaining the probability relation between the radial eccentricity and axial perpendicularity tolerance of each level of rotary equipment and the final coaxiality tolerance of the multi-level equipment, thereby realizing the distribution of the tolerance of the large-scale high-speed rotary equipment.
The existing methods all have the following problems: only the tolerance transmission process and the tolerance distribution are considered, and the measurement error of the single-stage rotor circular profile is not considered. The precise measurement of the circular profile is the premise and the basis for realizing the accurate estimation of the concentricity error of the single-stage rotor, the measurement error is amplified in the error transmission process, and the tolerance distribution is greatly influenced.
Disclosure of Invention
The invention aims to solve the problems of low measurement precision, low assembly coaxiality and poor assembly quality of a circular profile of an aircraft engine rotor, improve the performance of the aircraft engine, and provide a part tolerance distribution method and device based on four-parameter compensation.
The purpose of the invention is realized by the following technical scheme: a part tolerance distribution method based on four-parameter compensation establishes a part four-parameter circle profile measurement model, and the model comprises four parameters of a measured part eccentric error, a sensor measuring head offset, a sensor measuring ball radius and an included angle between a rotary axis and a geometric axis;
the measurement equation of the four-parameter circle profile measurement model is as follows:
in the formula, ρiThe distance from the measuring head of the sensor to the measuring rotation center, e is the eccentric error of the measured part, and thetaiα is the eccentric angle of the measured part, r is the radius of the sphere measured by the sensor, n is the number of sampling points, r is the sampling angle relative to the rotation centeriTo fit the distance of the ith sample point of the ellipse to the geometric center, Δ riProcessing errors of the surfaces of the parts, and d is the offset of a sensor measuring head;
fitting the distance r from the ith sampling point of the ellipse to the geometric centeriExpressed as:
in the formula, r0In order to fit the minor axis of the ellipse,
β is an included angle between the projection of the geometric axis on the measuring plane and the initial measuring direction, and gamma is an included angle between the revolution axis and the geometric axis;
when the eccentric error of the measured part meets the requirement of e/r relative to the minor axis of the fitting ellipse of the part0＜10-3And then, the measurement equation is expanded through power series, and a simplified four-parameter circular profile measurement model is obtained as follows:
estimating the eccentric error of the measured part according to a simplified four-parameter circle profile measurement model, obtaining a target function of the eccentric error of the measured part, generating 10000 groups of circle profile data of the measured part of the aero-engine according to a Monte Carlo method, substituting a random number into the target function of the eccentric error of the measured part, further obtaining 10000 groups of eccentric error parameters of the measured part, solving a probability density function according to a drawn distribution function, further obtaining a probability relation between contact surface jumping information and the eccentric error of the aero-engine part, and realizing the distribution of the tolerance of the aero-engine part.
Further, the part is a single-stage rotor.
The invention also provides a device for realizing the part tolerance distribution method based on four-parameter compensation, which comprises a base 1, an air floatation shaft system 2, a static balance measuring table 3, an aligning and inclination adjusting workbench 4, an upright post 5, a lower transverse measuring rod 6, a lower telescopic inductive sensor 7, an upper transverse measuring rod 8 and an upper lever inductive sensor 9; the air floatation shaft system 2 is nested on the central position of the base 1, the static balance measuring table 3 is arranged on the central position of the air floatation shaft system 2, wherein the static balance measuring table 3 comprises a static balance measuring table lower plate 3a, a static balance measuring table upper plate 3b, two guide rods 3c and 3d, four driving rods 3e, 3f, 3g and 3h and three precise force sensors 3i, 3j and 3 k; the two guide rods 3c and 3d are uniformly and equidistantly arranged on the lower plate 3a of the static balance measuring table; four driving rods 3e, 3f, 3g and 3h are uniformly and equidistantly arranged on a lower plate 3a of the static balance measuring table, an upper plate 3b of the static balance measuring table is sleeved on two guide rods 3c and 3d, and three precise force sensors 3i, 3j and 3k are uniformly and equidistantly arranged on an upper plate 3b of the static balance measuring table; the aligning and inclination adjusting workbench 4 is arranged at the center of the static balance measuring table 3, and the upright columns 5 are distributed at the left side of the air floatation shaft system 2 and are fixedly arranged on the base 1; an upper transverse measuring rod 8 and a lower transverse measuring rod 6 are sleeved on the upright post 5 in a moving and adjusting manner from top to bottom, and an upper lever type inductive sensor 9 is fixedly connected with the upper transverse measuring rod 8; the lower telescopic inductive sensor 7 is fixedly connected with the lower transverse measuring rod 6.
The method analyzes the measurement error in the process of measuring the circular profile of the rotor of the aircraft engine, and establishes a four-parameter circular profile measurement model by considering four parameter components of the rotor eccentric error, the sensor measuring head offset, the sensor measuring ball radius and the included angle between the revolution axis and the geometric axis in the circular profile measurement; and generating random numbers according to a Monte Carlo method, drawing a distribution function to obtain a probability density function, obtaining a probability relation between contact surface jumping and eccentric errors, and realizing single-stage rotor tolerance distribution.
Drawings
FIG. 1 is a flow chart of a part tolerance assignment method based on four-parameter compensation according to the present invention;
FIG. 2 is a schematic diagram of an apparatus for implementing a part tolerance assignment method based on four-parameter compensation according to the present invention;
fig. 3 is a schematic view of a static balance measuring table of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
With reference to fig. 1, the invention provides a part tolerance allocation method based on four-parameter compensation, which establishes a part four-parameter circular profile measurement model, wherein the model comprises four parameters of a measured part eccentric error, a sensor measuring head offset, a sensor measuring ball radius and an included angle between a rotary axis and a geometric axis;
the measurement equation of the four-parameter circle profile measurement model is as follows:
in the formula, ρiThe distance from the measuring head of the sensor to the measuring rotation center, e is the eccentric error of the measured part, and thetaiα is the eccentric angle of the measured part, r is the radius of the sphere measured by the sensor, n is the number of sampling points, r is the sampling angle relative to the rotation centeriTo fit the distance of the ith sample point of the ellipse to the geometric center, Δ riProcessing errors of the surfaces of the parts, and d is the offset of a sensor measuring head;
fitting the distance r from the ith sampling point of the ellipse to the geometric centeriExpressed as:
in the formula, r0In order to fit the minor axis of the ellipse,
β is an included angle between the projection of the geometric axis on the measuring plane and the initial measuring direction, and gamma is an included angle between the revolution axis and the geometric axis;
when the eccentric error of the measured part meets the requirement of e/r relative to the minor axis of the fitting ellipse of the part0＜10-3And then, the measurement equation is expanded through power series, and a simplified four-parameter circular profile measurement model is obtained as follows:
estimating the eccentric error of the measured part according to a simplified four-parameter circle profile measurement model, obtaining a target function of the eccentric error of the measured part, generating 10000 groups of circle profile data of the measured part of the aero-engine according to a Monte Carlo method, substituting a random number into the target function of the eccentric error of the measured part, further obtaining 10000 groups of eccentric error parameters of the measured part, solving a probability density function according to a drawn distribution function, further obtaining a probability relation between contact surface jumping information and the eccentric error of the aero-engine part, and realizing the distribution of the tolerance of the aero-engine part.
The parts are single-stage rotors.
With reference to fig. 2 and 3, the present invention further provides a device for implementing a part tolerance allocation method based on four-parameter compensation, where the device includes a base 1, an air-floating shafting 2, a static balance measuring table 3, a centering and tilt-adjusting worktable 4, an upright column 5, a lower transverse measuring rod 6, a lower telescopic inductive sensor 7, an upper transverse measuring rod 8, and an upper lever inductive sensor 9; the air floatation shaft system 2 is nested on the central position of the base 1, the static balance measuring table 3 is arranged on the central position of the air floatation shaft system 2, wherein the static balance measuring table 3 comprises a static balance measuring table lower plate 3a, a static balance measuring table upper plate 3b, two guide rods 3c and 3d, four driving rods 3e, 3f, 3g and 3h and three precise force sensors 3i, 3j and 3 k; the two guide rods 3c and 3d are uniformly and equidistantly arranged on the lower plate 3a of the static balance measuring table; four driving rods 3e, 3f, 3g and 3h are uniformly and equidistantly arranged on a lower plate 3a of the static balance measuring table, an upper plate 3b of the static balance measuring table is sleeved on two guide rods 3c and 3d, and three precise force sensors 3i, 3j and 3k are uniformly and equidistantly arranged on an upper plate 3b of the static balance measuring table; the aligning and inclination adjusting workbench 4 is arranged at the center of the static balance measuring table 3, and the upright columns 5 are distributed at the left side of the air floatation shaft system 2 and are fixedly arranged on the base 1; an upper transverse measuring rod 8 and a lower transverse measuring rod 6 are sleeved on the upright post 5 in a moving and adjusting manner from top to bottom, and an upper lever type inductive sensor 9 is fixedly connected with the upper transverse measuring rod 8; the lower telescopic inductive sensor 7 is fixedly connected with the lower transverse measuring rod 6.
The method and the device for distributing the part tolerance based on the four-parameter compensation are described in detail, and a specific example is applied to illustrate the principle and the implementation mode of the invention, and the description of the embodiment is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
Claims (3)
1. A part tolerance distribution method based on four-parameter compensation is characterized in that,
establishing a four-parameter circular profile measurement model of the part, wherein the model comprises four parameters of the eccentric error of the part to be measured, the offset of a sensor measuring head, the radius of a sensor measuring ball and the included angle between the rotary axis and the geometric axis;
the measurement equation of the four-parameter circle profile measurement model is as follows:
in the formula, ρiThe distance from the measuring head of the sensor to the measuring rotation center, e is the eccentric error of the measured part, and thetaiα is the eccentric angle of the measured part, r is the radius of the sphere measured by the sensor, n is the number of sampling points, r is the sampling angle relative to the rotation centeriTo fit the distance of the ith sample point of the ellipse to the geometric center, Δ riProcessing errors of the surfaces of the parts, and d is the offset of a sensor measuring head;
fitting the distance r from the ith sampling point of the ellipse to the geometric centeriExpressed as:
in the formula, r0In order to fit the minor axis of the ellipse,
β is an included angle between the projection of the geometric axis on the measuring plane and the initial measuring direction, and gamma is an included angle between the revolution axis and the geometric axis;
when the eccentric error of the measured part meets the requirement of e/r relative to the minor axis of the fitting ellipse of the part0＜10-3And then, the measurement equation is expanded through power series, and a simplified four-parameter circular profile measurement model is obtained as follows:
estimating the eccentric error of the measured part according to a simplified four-parameter circle profile measurement model, obtaining a target function of the eccentric error of the measured part, generating 10000 groups of circle profile data of the measured part of the aero-engine according to a Monte Carlo method, substituting a random number into the target function of the eccentric error of the measured part, further obtaining 10000 groups of eccentric error parameters of the measured part, solving a probability density function according to a drawn distribution function, further obtaining a probability relation between contact surface jumping information and the eccentric error of the aero-engine part, and realizing the distribution of the tolerance of the aero-engine part.
2. The method of claim 1, wherein the component is a single stage rotor.
3. A device for realizing a part tolerance distribution method based on four-parameter compensation is characterized by comprising a base (1), an air floatation shaft system (2), a static balance measuring table (3), a centering and inclination adjusting workbench (4), a stand column (5), a lower transverse measuring rod (6), a lower telescopic inductive sensor (7), an upper transverse measuring rod (8) and an upper lever inductive sensor (9); the air floatation shaft system (2) is nested in the center of the base (1), the static balance measuring table (3) is arranged in the center of the air floatation shaft system (2), and the static balance measuring table (3) comprises a static balance measuring table lower plate (3a), a static balance measuring table upper plate (3b), two guide rods (3c and 3d), four driving rods (3e and 3f and 3g and 3h) and three precise force sensors (3i and 3j and 3 k); the two guide rods (3c, 3d) are uniformly and equidistantly arranged on the lower plate (3a) of the static balance measuring table; four driving rods (3e, 3f, 3g, 3h) are uniformly and equidistantly arranged on a lower plate (3a) of a static balance measuring table, an upper plate (3b) of the static balance measuring table is sleeved on two guide rods (3c, 3d), and three precise force sensors (3i, 3j, 3k) are uniformly and equidistantly arranged on an upper plate (3b) of the static balance measuring table; the aligning and inclination adjusting workbench (4) is arranged at the center of the static balance measuring table (3), and the upright columns (5) are distributed on the left side of the air floatation shaft system (2) and are fixedly arranged on the base (1); an upper transverse measuring rod (8) and a lower transverse measuring rod (6) are sleeved on the upright post (5) from top to bottom in a sequentially moving and adjusting manner, and an upper lever type inductive sensor (9) is fixedly connected with the upper transverse measuring rod (8); the lower telescopic inductive sensor (7) is fixedly connected with the lower transverse measuring rod (6).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811231186.6A CN111076656A (en) | 2018-10-22 | 2018-10-22 | Part tolerance distribution method and device based on four-parameter compensation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811231186.6A CN111076656A (en) | 2018-10-22 | 2018-10-22 | Part tolerance distribution method and device based on four-parameter compensation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111076656A true CN111076656A (en) | 2020-04-28 |
Family
ID=70309886
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811231186.6A Pending CN111076656A (en) | 2018-10-22 | 2018-10-22 | Part tolerance distribution method and device based on four-parameter compensation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111076656A (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004072600A2 (en) * | 2003-02-11 | 2004-08-26 | Durr Production Systems, Inc. | Apparatus for measuring radial displacement of a wheel |
| CN101339021A (en) * | 2008-08-13 | 2009-01-07 | 哈尔滨工业大学 | Double offset parameter circle contour measurement model and biased error separation method |
| CN103970965A (en) * | 2014-05-26 | 2014-08-06 | 中国航空动力机械研究所 | Test run method for accelerated life test of gas turbine engine |
| CN105426566A (en) * | 2015-10-15 | 2016-03-23 | 哈尔滨工业大学 | Monte Carlo technique based geometrical tolerance distribution method for large high-speed rotary equipment |
| CN105443169A (en) * | 2015-11-11 | 2016-03-30 | 上海交通大学 | Device and method for assembling multi-stage drum type rotor of aero-engine |
| CN105608263A (en) * | 2015-12-17 | 2016-05-25 | 北京航空航天大学 | Adaptive processing method oriented to service life probability analysis of turbine leaf disc structure |
| CN205272429U (en) * | 2015-12-24 | 2016-06-01 | 浙江澳尔法机械制造有限公司 | Static balance mould |
| CN206891656U (en) * | 2017-05-26 | 2018-01-16 | 广州市番禺远航螺旋桨制造有限公司 | Adjustable pitch propeller static balance measurement apparatus |
-
2018
- 2018-10-22 CN CN201811231186.6A patent/CN111076656A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004072600A2 (en) * | 2003-02-11 | 2004-08-26 | Durr Production Systems, Inc. | Apparatus for measuring radial displacement of a wheel |
| CN101339021A (en) * | 2008-08-13 | 2009-01-07 | 哈尔滨工业大学 | Double offset parameter circle contour measurement model and biased error separation method |
| CN103970965A (en) * | 2014-05-26 | 2014-08-06 | 中国航空动力机械研究所 | Test run method for accelerated life test of gas turbine engine |
| CN105426566A (en) * | 2015-10-15 | 2016-03-23 | 哈尔滨工业大学 | Monte Carlo technique based geometrical tolerance distribution method for large high-speed rotary equipment |
| CN105443169A (en) * | 2015-11-11 | 2016-03-30 | 上海交通大学 | Device and method for assembling multi-stage drum type rotor of aero-engine |
| CN105608263A (en) * | 2015-12-17 | 2016-05-25 | 北京航空航天大学 | Adaptive processing method oriented to service life probability analysis of turbine leaf disc structure |
| CN205272429U (en) * | 2015-12-24 | 2016-06-01 | 浙江澳尔法机械制造有限公司 | Static balance mould |
| CN206891656U (en) * | 2017-05-26 | 2018-01-16 | 广州市番禺远航螺旋桨制造有限公司 | Adjustable pitch propeller static balance measurement apparatus |
Non-Patent Citations (1)
| Title |
|---|
| 孙传智: "基于矢量投影的多级转子同轴度测量方法研究", 《中国博士学位论文全文数据库》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9695707B2 (en) | Five-degree-of-freedom adjustment and positioning method and apparatus for assembly/measurement of rotor and stator of aircraft engine | |
| 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 | |
| WO2015120744A1 (en) | Aircraft engine rotor assembly method and device based on concentricity and perpendicularity measurement | |
| CN111413031B (en) | Deep learning regulation and assembly method and device for large-scale high-speed rotation equipment based on dynamic vibration response characteristics | |
| CN103790647B (en) | Hydraulic pressure based on inductance sensing grabs cassette aeroengine rotor assembly apparatus | |
| CN103899367B (en) | The stacking assembly method of aeroengine rotor and device | |
| CN102997846A (en) | Full-airplane horizontal measurement method based on work space measurement location system | |
| CN103806958B (en) | Hydraulic pressure based on inductosyn grabs cassette aeroengine rotor assembly apparatus | |
| CN108340210A (en) | A kind of gang tool geometric error discrimination method measured based on ball bar | |
| CN105426565A (en) | Probability density technology based geometrical tolerance distribution method for large high-speed rotary equipment | |
| CN110598229B (en) | Large-scale high-speed rotation equipment multi-stage part optimization method based on Monte Carlo deviation evaluation | |
| CN103790653B (en) | Based on the planer-type aeroengine rotor assembly apparatus of gas magnetic composite support | |
| CN109871947A (en) | Multistage components initial unbalance, stacking method and device step by step are equipped in large high-speed revolution based on convolutional neural networks | |
| CN110906862B (en) | Geometric morphology and quality characteristic integrated measuring device for large-scale high-speed rotation equipment | |
| CN103790651A (en) | Aircraft engine rotor assembling method and device based on air/magnetism mixing guiding | |
| CN111891176A (en) | Track measurement data integral adjustment method based on track geometric state measuring instrument | |
| CN111076656A (en) | Part tolerance distribution method and device based on four-parameter compensation | |
| CN110595689B (en) | Large-scale high-speed rotation equipment multistage part unbalance amount prediction method and part assembly device | |
| CN110595413B (en) | Part tolerance distribution method and device based on five-parameter compensation | |
| CN110608666B (en) | Aero-engine rotor assembly measuring device based on four-point weighing and three-target optimization method | |
| CN110608667B (en) | Three-point weighing-based aeroengine rotor assembly measuring device and three-target optimization method | |
| CN103790648B (en) | The aeroengine rotor assembly method optimized based on multi-part concentricity and device | |
| CN203031075U (en) | Linear cutting machine of gas turbine outer ring dovetail groove | |
| CN103776368B (en) | Gas and magnetism composite supporting type aero-engine rotor assembling method and device based on concentricity optimization |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200428 |