CN108388747B - The multichannel circumferential direction class Sine distribution sample implementation method of error of fixed angles blade - Google Patents
The multichannel circumferential direction class Sine distribution sample implementation method of error of fixed angles blade Download PDFInfo
- Publication number
- CN108388747B CN108388747B CN201810198451.9A CN201810198451A CN108388747B CN 108388747 B CN108388747 B CN 108388747B CN 201810198451 A CN201810198451 A CN 201810198451A CN 108388747 B CN108388747 B CN 108388747B
- Authority
- CN
- China
- Prior art keywords
- error
- blade
- fixed angles
- random
- multichannel
- 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.)
- Active
Links
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
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
A kind of multichannel circumferential direction class Sine distribution sample implementation method of error of fixed angles blade, pass through the setting blade angle margin of tolerance and Gaussian Profile probability density function, it is random to generate error of fixed angles and obtain multichannel analog sample, circumferential class Sine distribution sample is obtained after resetting.The present invention combines sine curve with the circumferential variation tendency of blade angle error, and to reduce the deterioration degree of Capability of Compressor caused by error of fixed angles, the processing cost of compressor blade is controlled while adaptation turbomachine high performance demand for development.
Description
Technical field
The present invention relates to a kind of technology in turbomachine field, the multichannel of specifically a kind of error of fixed angles blade
Circumferential class Sine distribution sample implementation method.
Background technique
In aviation engine shaft stream high-pressure compressor, blade is to realize that the part of energy conversion and quantity are most, add
The most complicated part of work, the compressor efficiency and aerodynamic stability that the processing quality of blade designs current high load capacity are played to pass
Important role.However due to the complexity of blade processing, the geometry that actual processing obtains easily deviates design parameter and generates mistake
The error of fixed angles of difference, especially compressor blade.A large amount of error of fixed angles blade will will lead to compressor aeroperformance
Sharply decline, especially the violent reduction of surge margin.
In actual processing, by reducing predetermined tolerance range as much as possible, blade processing error can be effectively reduced
Size, to reduce performance loss, but this mode will lead to dramatically increasing for processing cost.If can be in blade error range
Under the premise of constant, the performance loss of compressor is reduced by other methods, will bring weight to the manufacturing field of compressor blade
Big meaning.
Summary of the invention
The present invention In view of the above shortcomings of the prior art, proposes a kind of multichannel circumferential direction class of error of fixed angles blade
Sine distribution sample implementation method combines sine curve with the circumferential variation tendency of blade angle error, to reduce peace
The deterioration degree for filling Capability of Compressor caused by angle error, control pressure while adaptation turbomachine high performance demand for development
The processing cost of mechanism of qi blade.
The present invention is achieved by the following technical solutions:
The present invention is random to generate by the setting blade angle margin of tolerance and Gaussian Profile probability density function (PDF)
Error of fixed angles simultaneously obtains multichannel analog sample, obtains circumferential class Sine distribution sample after resetting.
The blade angle margin of tolerance refers to: the criterion at different levels with reference to as defined in navigation mark determine compressor blade
The established angle margin of tolerance uses but is not limited to [- 1.5 ° ,+1.5 °].
The Gaussian Profile probability density function is established under the selected margin of tolerance
Wherein: μ is the average value of geometric error, and being taken as 0, σ is standard deviation, depending on the margin of tolerance；Meet the distribution of geometrical deviation:
State of the probability density close to 0 when deviateing that prototype is remoter, and probability density is lower, and reaching tolerance boundary.
The error of fixed angles constructs multiple minizones, is based on probability by the way that the margin of tolerance is carried out classifying rationally
Cumulative distribution function corresponding to density functionCalculate the generation of each minizone
Probability；Total radix of vane manufacturing needed for determining, calculates separately each cell by the probability of happening value of obtained each minizone
The interior practical number of blade, then constructs one group of random number using random function in these minizones, and the quantity of random number is
The quantity of the section intra vane, these random numbers are the error of fixed angles of each blade.
The multichannel analog sample, by determining that the calculating cycle of compressor blade (is equivalent to and determines the logical of calculating
Road number).It by out-of-order random alignment after the vaned error of fixed angles mixing of institute and is grouped, every group of the number of blade is to calculate channel
Number, thus obtains the random combine of a large amount of error of fixed angles blade, as the multichannel analog sample of error of fixed angles blade.
The circumferential class Sine distribution mode refers to: the peace by redistributing each blade in multichannel analog sample
It fills angle error and carries out distribution rearrangement by reference object of sine curve, selecting sine curve equation is y=Asin [α (x-1)],
Middle x and y is respectively that blade is circumferentially arranged serial number and error of fixed angles value, and amplitude A is initially selected margin of tolerance boundary, and α is then
It is the related coefficient depending on compressor blade sum and calculating cycle, then this equation is based on the circumferential direction for updating error blade
Circumferential direction class Sine distribution sample obtained from distribution sequence.
Technical effect
Compared with prior art, the present invention matches with actual blade processing and installation situation, for the peace for studying blade
Dress angle error provides more true data result to the affecting laws of Capability of Compressor.Meanwhile this circumferentially class sine
Continuity error of fixed angles Leaf positional distribution mode, can effectively reduce the feelings of the big change of gradient of error of fixed angles between adjacent blades
Condition, reduction adjacency channel flow field unevenly changes and the probability of happening of adjacent blades load mismatch case, so that adjacent logical
The degree of flow field change is extremely subtle between road, more consistent with design point, to reduce evil of the aeroperformance relative to prototype
Change degree.
Detailed description of the invention
Fig. 1 is compressor rotor schematic diagram of the present invention；
Fig. 2 defines and its generates geometric representation when deviation for certain compressor rotor blade established angle of embodiment；
In figure: γ indicates that established angle, c indicate that chord length, number are error of fixed angles size example；
Fig. 3 is Gaussian Profile probability density function (PDF) curve of embodiment in a certain margin of tolerance；
Fig. 4 is a certain multichannel analog sample of the invention and its corresponding class Sine distribution pattern diagram.
Specific embodiment
As shown in Figure 1, the present embodiment specifically includes the following steps:
The selection of the step 1. blade angle margin of tolerance and PDF
1.1 firstly, for this compressor blade, and selecting its established angle margin of tolerance is [- 1.5 ° ,+1.5 °].
1.2 under the selected margin of tolerance, establish reasonable Gaussian Profile probability density function, meet geometry distribution:
State of the probability density close to 0 when deviateing that prototype is remoter, and probability density is lower, and reaching tolerance boundary, therefore this example
In Gaussian Profile probability density function, standard deviation sigma is taken as 0.5, then functional image is as shown in Figure 3.
The random generation of step 2. error of fixed angles
2.1, by [- 1.5 ° ,+1.5 °] progress classifying rationallies of the margin of tolerance, construct multiple minizones, divide in this example
For [- 1.5 °, -0.5 °], [- 0.5 ° ,+0.5 °], [+0.5 ° ,+1.5 °] three sections carry out subsequent calculating.
2.2 based on cumulative distribution function corresponding to the probability density function in step 1, these three minizones are calculated
Probability of happening be respectively 0.16,0.68,0.16.
2.3 simultaneously, total radix of vane manufacturing needed for determining, in this instance, for convenience of description, it is assumed that vane manufacturing sum
It is 100, then by above-mentioned probability value, the practical number of blade obtained in three sections is respectively 16,68,16.
2.4 then, using random function, in these three minizones, construct the random number of one group of respective numbers, i.e., respectively
16 random numbers are generated respectively in [- 1.5 °, -0.5 °] and [+0.5 ° ,+1.5 °] section, in [- 0.5 ° ,+0.5 °] section then
Generate 68 random numbers.These random numbers are to represent the error of fixed angles size of each blade.
The generation of step 3. multichannel analog sample
3.1 determine the calculating cycle (being equivalent to the port number for determining calculating) of compressor blade.In the present embodiment, rotor
Blade amt is 64, according to actual numerical computation, selects for 1/4 period for calculating cycle in this example, i.e. calculating sample
Port number be selected as 16 channels.
100 random error blades obtained in step 2 are mixed simultaneously random ordering by 3.2, are carried out random alignment and are simultaneously grouped, every group
The number of blade is calculating port number, therefore is 16, it is hereby achieved that at least 6 groups of error of fixed angles blade random combine, this
A little combinations are the multichannel analog sample of error of fixed angles blade.
3.3 if desired more geometry samples need to only increase the total radix of blade, and the step of repeating 2.3,2.4 and 3.2
?.
The realization of step 4. circumferential direction class Sine distribution mode
4.1 in the multichannel analog sample that step 3 obtains, and error of fixed angles blade circumferentially is random distribution,
The circumferential array sequence for now redistributing its installation, by the error of fixed angles size of each blade using sine curve as reference object
It is distributed.Selected sine curve equation is y=Asin [α (x-1)], and amplitude A is initially selected margin of tolerance boundary, therefore
It is then the related coefficient for depending on compressor blade total (64) and its calculating cycle (1/4 period) for 1.5, α in this example, because
It is π/8 in this this example, so that the distribution of blade error at least meets 1 week of SIN function in the simulation in this 16 channel
Phase (as shown in Figure 4).
4.2 with equationOn the basis of (see Fig. 4 black curve) update error blade circumferential direction
Multichannel sample obtained from distribution sequence is class Sine distribution sample.For a certain random geometry sample in selection 3.2
(see Fig. 4 red open column shape figure), the sequence of each of which blade is rearranged according to respective error size, obtains figure
4 solid black histograms, variation tendency circumferentially substantially conform to above-mentioned sine curve, i.e., the class of random sample is sinusoidal thus
It is distributed sample, therefore such location mode is known as " class Sine distribution method ".
Indexes of Evaluation Effect of the invention
Step 4 of the invention gives maximum contribution to final effect, so that the deterioration degree of overall performance substantially reduces.
In the present embodiment, the assessment of this method is carried out using isentropic efficiency, surge margin and opposite total pressure loss coefficient as performance indicator,
Class Sine distribution mode, the results of property of random distribution mode and prototype are compared, as shown in table 1.Seen from table 1, relative to original
Type, error of fixed angles cause the deterioration of Blade Properties, still, by being the distribution pattern of circumferential class sine by blade control,
Three kinds of indicators of overall performance have relative to the result of random distribution to be extremely obviously improved.
The comparison of each performance indicator of table 1
Compared with prior art, the present invention combines sine curve with the circumferential variation tendency of blade angle error,
Considerably reduce the deterioration degree of Capability of Compressor.In compressor process industry, if applying the present invention to practical blade
Processing installation process, that is, practical blade installation when, apply class sine circumferentially distributed method, then not reducing tolerance
In the case where the even slight increase margin of tolerance of range, error of fixed angles can be effectively reduced, the deterioration of Capability of Compressor is made
With reduction manufacturing cost also meets performance requirement.
Above-mentioned specific implementation can by those skilled in the art under the premise of without departing substantially from the principle of the invention and objective with difference
Mode carry out local directed complete set to it, protection scope of the present invention is subject to claims and not by above-mentioned specific implementation institute
Limit, each implementation within its scope is by the constraint of the present invention.
Claims (4)
1. a kind of multichannel circumferential direction class Sine distribution sample implementation method of error of fixed angles blade, which is characterized in that by setting
The blade angle margin of tolerance and Gaussian Profile probability density function are set, it is random to generate error of fixed angles and obtain multichannel analog
Sample obtains circumferential class Sine distribution sample after resetting；
The error of fixed angles constructs multiple minizones, is based on probability density by the way that the margin of tolerance is carried out classifying rationally
Cumulative distribution function corresponding to functionThe generation for calculating each minizone is general
Rate；Total radix of vane manufacturing needed for determining, calculates separately each minizone by the probability of happening value of obtained each minizone
Then the interior practical number of blade constructs one group of random number using random function in these minizones, the quantity of random number is should
The quantity of section intra vane, these random numbers are the error of fixed angles of each blade.
2. according to the method described in claim 1, it is characterized in that, the Gaussian Profile probability density function, in selected public affairs
Under poor range, establishWherein: μ is the average value of geometric error, and being taken as 0, σ is standard deviation.
3. according to the method described in claim 1, it is characterized in that, the multichannel analog sample, pass through determine gas compressor blade
The calculating cycle of piece；It by out-of-order random alignment after the vaned error of fixed angles mixing of institute and is grouped, every group of the number of blade is
Port number is calculated, the random combine of a large amount of error of fixed angles blade, the as multichannel of error of fixed angles blade are thus obtained
Analog sample.
4. according to the method described in claim 1, it is characterized in that, the circumferential class Sine distribution mode refers to: by again
The error of fixed angles of each blade carries out distribution rearrangement by reference object of sine curve in distribution multichannel analog sample, selectes
Sine curve equation is y=A sin [α (x-1)], and wherein x and y is respectively that blade is circumferentially arranged serial number and error of fixed angles value,
Amplitude A is initially selected margin of tolerance boundary, and α is then the related coefficient depending on compressor blade sum and calculating cycle,
Then this equation is circumferential direction class Sine distribution sample obtained from the circumferentially distributed sequence of datum renewal error blade.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810198451.9A CN108388747B (en) | 2018-03-12 | 2018-03-12 | The multichannel circumferential direction class Sine distribution sample implementation method of error of fixed angles blade |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810198451.9A CN108388747B (en) | 2018-03-12 | 2018-03-12 | The multichannel circumferential direction class Sine distribution sample implementation method of error of fixed angles blade |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108388747A CN108388747A (en) | 2018-08-10 |
| CN108388747B true CN108388747B (en) | 2019-08-02 |
Family
ID=63067376
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810198451.9A Active CN108388747B (en) | 2018-03-12 | 2018-03-12 | The multichannel circumferential direction class Sine distribution sample implementation method of error of fixed angles blade |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108388747B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111241653B (en) * | 2019-04-19 | 2024-03-19 | 上海交通大学 | Multi-dimensional blade error performance influence simulation implementation method based on Monte Carlo method |
| CN110059414B (en) * | 2019-04-22 | 2020-09-29 | 北京理工大学 | Two-dimensional blade modeling method for directly controlling channel |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0903404D0 (en) * | 2009-03-02 | 2009-04-08 | Rolls Royce Plc | Surface profile evaluation |
| UA99876C2 (en) * | 2011-05-19 | 2012-10-10 | Мита-Текник А/С | Method for control of orientation of wind turbine and wind turbine |
| CN107489651B (en) * | 2017-10-10 | 2019-05-07 | 北京航空航天大学 | A kind of blade profile optimization method that can inhibit fan shock wave noise based on quadratic function |
| CN108804853B (en) * | 2018-06-28 | 2020-09-15 | 东北大学 | Elastic support lower torsional shoulder blade dynamic modeling method based on variable cross-section beam |
| CN109340042A (en) * | 2018-11-19 | 2019-02-15 | 内蒙古工业大学 | Turbine blade and its design method and solar chimney electricity generation system |
-
2018
- 2018-03-12 CN CN201810198451.9A patent/CN108388747B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108388747A (en) | 2018-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108388747B (en) | The multichannel circumferential direction class Sine distribution sample implementation method of error of fixed angles blade | |
| CN106126925B (en) | A method of it improving reactor core three-dimensional netron-flux density and is finely distributed | |
| CN104866692A (en) | Aircraft multi-objective optimization method based on self-adaptive agent model | |
| CN103942599A (en) | Particle swarm optimization method based on survival of the fittest and step-by-step selection | |
| CN111859746B (en) | Method for predicting variable working condition performance of turbomachinery based on flow field reconstruction | |
| CN108108528B (en) | One-dimensional matching design method for power turbine of split-shaft type aeroderivative | |
| CN112270139B (en) | Pneumatic optimization design method for centrifugal compressor of fuel cell based on mother type library | |
| CN114254460B (en) | Turbomachine pneumatic robustness optimization method based on graph convolution neural network | |
| CN107357959A (en) | A kind of high-precision blade profile characteristic parameter extraction method | |
| CN110110349A (en) | A kind of rotary blade type mixed-flow pump spatial guide blade optimization design under multiple working conditions method | |
| CN110288196A (en) | Transmission line galloping grading forewarning system method and system | |
| Youngren | Analysis and design of transonic cascades with splitter vanes | |
| CN107679336A (en) | Reflector antenna surface random error analysis method based on Two-order approximation formula | |
| CN110532726A (en) | A kind of non local life expectance appraisal procedure of the most weak ring of the turbine disk based on Bayes's calibration | |
| CN111611646B (en) | Rapid design method for aerodynamic configuration of aircraft precursor | |
| CN112487570A (en) | Centrifugal compressor shunting blade shape optimization method based on free deformation technology | |
| Lyes et al. | Experimental evaluation of the high-to-low speed transformation process for a highly loaded core compressor stage | |
| CN105488521B (en) | A kind of dilatation screening sample method based on kernel function | |
| CN105916113B (en) | A kind of localization method and device | |
| CN109753695B (en) | Method for expanding characteristics of engine parts | |
| CN105760647A (en) | Method for calculating mass weighted mean | |
| CN109086245B (en) | VIC parameter multi-target method based on substitution model | |
| Garza et al. | Evaluation of two novel wake models in offshore wind farms | |
| KR101890494B1 (en) | Method for prediction the solar proton particle | |
| Mulloth et al. | A high fidelity quality assessment of high pressure turbine blades using surface curvature and gradient-based adjoint |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |