CN103077315A - Method for rebuilding front and back oval edges of aerial blade based on error control - Google Patents
Method for rebuilding front and back oval edges of aerial blade based on error control Download PDFInfo
- Publication number
- CN103077315A CN103077315A CN2013100108308A CN201310010830A CN103077315A CN 103077315 A CN103077315 A CN 103077315A CN 2013100108308 A CN2013100108308 A CN 2013100108308A CN 201310010830 A CN201310010830 A CN 201310010830A CN 103077315 A CN103077315 A CN 103077315A
- Authority
- CN
- China
- Prior art keywords
- fitting
- point
- ellipse
- oval
- error
- 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.)
- Granted
Links
Images
Landscapes
- Feedback Control In General (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses a method for rebuilding front and back oval edges of an aerial blade based on error control, which is used for solving the technical problem of poor anti-noise performance existing in the conventional high-precision rebuilding method for a front edge of a blade. According to the technical scheme, the method comprises the following steps of: selecting region points covering the front edge or the back edge of the blade; sequencing the region points; performing spline fitting on the sequenced region points for three times; performing equal-arc-length encryption discretion on fitted spline curves, and marking beginning and end discrete points as P0i and P1i respectively; fitting discrete points at a P0iP1i section by using a least square method, and marking beginning and end discrete points on the intersection of the P0iP1i section and a fitting oval as P0(i+1) and P1(i+1) respectively; calculating an average error epsiloni ranging from the discrete points to an ith fitting oval; and setting an iteration control error as epsilon till the epsiloni is less than or equal to epsilon to finish iteration. According to the method, the front and back oval edges of the blade are rebuilt by adopting an error control-based iterative fitting method on the basis of fitting through the least square method, so that points on a non-oval circular arc in the selected region are prevented from participating in fitting, and the anti-noise performance is improved.
Description
Technical field
The present invention relates to the oval front and rear edge method for reconstructing of a kind of aerial blade, be specifically related to the oval front and rear edge method for reconstructing of a kind of aerial blade based on error control.
Background technology
The shape of blade front and rear edge plays key effect to the aeroperformance of whole blade.Traditional front and rear edge generally is designed to circular arc, and along with going deep into that blade is studied, numerous experiments and numerical value studies show that, adopt non-circular arc front edge, can obviously improve the aeroperformance of blade such as oval leading edge in recent years.
In the process that oval front and rear edge is rebuild based on the blade of measurement data, the selection area point for the treatment of match may comprise the point on leaf basin, the blade back, if the point on this non-elliptic arc is participated in match, can rebuild front and rear edge and have a huge impact.Document " high-precision reconstruction of blade leading edge research " aviation power journal " (2006; 21 (4): 722-726) " discloses a kind of high-precision reconstruction of blade leading edge, the method proposes based on the least square method of belt restraining data point to be carried out match, then use the result of match as initial value, determine oval center by the limit of ellipse-polar curve character, and then obtain oval characteristic parameter; Simultaneously, utilize the characteristics that lattice point increases error of fitting greatly that go out on the non-elliptic arc in the literary composition, determine lattice point: centered by leading edge point, carry out ellipse fitting after the enhancing data point, calculate the maximum relative error of oval feature parameter; Wherein, relative error can reduce along with increasing of normal point, if error of fitting increases manyfold (can set threshold values) before and after increasing certain point, and continues to increase data point, and error does not obviously reduce, and this point namely is lattice point; Point out in the literary composition that the error that random noise produces does not possess above-mentioned Variation Features; Therefore, for the measurement data of leaf packet Noise, utilize the method to carry out oval front and rear edge reconstruction and have certain difficulty.
Summary of the invention
In order to overcome the existing poor deficiency of high-precision reconstruction of blade leading edge noiseproof feature, the invention provides the oval front and rear edge method for reconstructing of a kind of aerial blade based on error control.The method is carried out in least square method on the basis of match for the measurement data of aerial blade, adopts the iterative fitting method based on error control, and the oval front and rear edge of blade is rebuild.By specification error the iterative fitting process is controlled, avoided as much as possible the participation of the point on non-elliptic arc match in the selection area; Under the prerequisite that guarantees reconstruction precision, have preferably noiseproof feature, can improve the reconstruction quality of the oval front and rear edge of blade.
The technical solution adopted for the present invention to solve the technical problems is: the oval front and rear edge method for reconstructing of a kind of aerial blade based on error control is characterized in may further comprise the steps:
(1) reads in the discrete measurement data of molded line of aerial blade, the selected region point that covers blade inlet edge or trailing edge;
(2) region point is sorted;
(3) region point after the ordering is carried out Cubic Spline Fitting;
(4) to the SPL of institute's match carry out etc. arc length encrypt discrete, and with the whole story discrete point be designated as respectively P
0i, P
1i
(5) with least square method to P
0iP
1iThe section discrete point carries out match, and with P
0iP
1iSection and fitted ellipse intersection the whole story discrete point be designated as respectively P
0 (i+1), P
1 (i+1)Ellipse fitting method is as follows:
1. as quafric curve f (x, y)=ax
2+ 2bxy+cy
2+ 2dx+2fy+g=0 the b that satisfies condition
2-ac<0 o'clock, quafric curve f (x, y) is expressed as ellipse.In the formula, a, b, c, d, f, g are the coefficient of quadratic polynomial equation.
Be provided with n data point P
i(x
i, y
i) participate in ellipse fitting, will be converted into based on the ellipse fitting of algebraic distance the least square problem of following belt restraining, the objective function of setting up ellipse fitting is as follows:
2. adopt least square method to try to achieve the coefficient a of quadratic polynomial equation, b, c, d, f behind the g, obtains fitted ellipse, and the elliptic geometry calculation of parameter is as follows:
The oval center of circle (x
0, y
0):
Transverse r
a:
Ellipse short shaft r
b:
In the formula: j=1,2,3 ..., N, d
jBe the distance that j discrete point arrives gained fitted ellipse in the step (5), N is for participating in counting out of match;
(7) setting the iteration control error is ε: to ε
iCompare with ε, if ε
i≤ ε calculates and finishes; Otherwise, work as ε
iDuring ε, be positioned at P after casting out the i time iteration
0 (i+1)P
0iAnd P
1 (i+1)P
1iBetween point on the non-elliptic arc; Make i=i+1, to P
0iP
1iBetween the match of discrete point repeating step (5), until satisfy ε
iDuring≤ε, iteration finishes.
The invention has the beneficial effects as follows: because the method is carried out in least square method on the basis of match for the measurement data of aerial blade, adopt the iterative fitting method based on error control, the oval front and rear edge of blade is rebuild.By specification error the iterative fitting process is controlled, avoided as much as possible the participation of the point on non-elliptic arc match in the selection area; Under the prerequisite that guarantees reconstruction precision, have preferably noiseproof feature, improved the reconstruction quality of the oval front and rear edge of blade.
Below in conjunction with drawings and Examples the present invention is elaborated.
Description of drawings
Fig. 1 is the oval front and rear edge method for reconstructing of the aerial blade process flow diagram that the present invention is based on error control.
Fig. 2 the present invention is based on the oval front and rear edge method for reconstructing of the aerial blade embodiment of error control through the reconstruction design sketch of seven matches.
Fig. 3 is the synoptic diagram that the present invention is based on the oval front and rear edge method for reconstructing of the aerial blade first fit of error control.
Fig. 4 is the synoptic diagram that the present invention is based on the i time match of the oval front and rear edge method for reconstructing of aerial blade of error control.
Fig. 5 is the synoptic diagram that the present invention is based on the final match of the oval front and rear edge method for reconstructing of aerial blade of error control.
Embodiment
With reference to Fig. 1~5.The oval front and rear edge method for reconstructing of the aerial blade concrete steps that the present invention is based on error control are as follows:
(1) reads in the discrete measurement data of molded line of blade, the selected region point that covers blade inlet edge or trailing edge;
(2) region point is sorted;
(3) region point after the ordering is carried out Cubic Spline Fitting;
(4) to the SPL of institute's match carry out etc. arc length encrypt discrete, and with the whole story discrete point be designated as respectively P
0i, P
1i
(5) with least square method to P
0iP
1iThe section discrete point carries out match, and with P
0iP
1iSection and fitted ellipse intersection the whole story discrete point be designated as respectively P
0 (i+1), P
1 (i+1)Ellipse fitting method is as follows:
1. set up the objective function of ellipse fitting:
In the plane, by the quadratic polynomial equation
F (x, y)=ax
2+ 2bxy+cy
2The curve that+2dx+2fy+g=0 is represented is called quafric curve.
In the formula, a,, c, d, f, g are the coefficient of quadratic polynomial equation.In above-mentioned equation, satisfy condition: b
2-ac<0 o'clock, quafric curve f (x, y) is expressed as ellipse.
Be provided with n data point P
i(x
i, y
i) participate in ellipse fitting, will be converted into based on the ellipse fitting of algebraic distance the least square problem of following belt restraining, the objective function of setting up ellipse fitting is as follows:
2. adopt least square method to try to achieve the coefficient a of quadratic polynomial equation,, c, d, f behind the g, obtains fitted ellipse, and the elliptic geometry calculation of parameter is as follows:
The oval center of circle (x
0, y
0):
Transverse r
a:
Ellipse short shaft r
b:
In the formula: j=1,2,3 ..., N, d
jBe the distance that j discrete point arrives gained fitted ellipse in (5), N is for participating in counting out of match;
(7) setting the iteration control error is ε: to ε
iCompare with ε, if ε
i≤ ε calculates and finishes; Otherwise, work as ε
iDuring ε, be positioned at P after casting out the i time iteration
0 (i+1)P
0iAnd P
1 (i+1)P
1iBetween point on the non-elliptic arc; Make i=i+1, to P
0iP
1iBetween the match of discrete point repeating step (5), until satisfy ε
iDuring≤ε, iteration finishes.
Preferred embodiment.Oval O
7For satisfying the fitted ellipse of specification error requirement, i.e. O
7Be its center, R
a, R
bBe respectively its long and short axle radius.Arbitrarily selected certain blade measurement data, concrete steps are as follows:
(1) reads in the discrete measurement data of molded line of blade, choose the region point that covers leading edge;
(2) region point is sorted, data point is from 0,1,2 ..., N;
(3) adopt cubic spline to carry out curve fitting, get curve L
0
(4) to the SPL L of institute's match
0The arc-length methods such as employing are encrypted discrete, get 1000 points, i=1, and the whole story, discrete point was designated as respectively P
0i, P
1iThe arc length such as the uniform sampling method of curved surface has, etc. parameter, etc. the action method, in same sampling precision situation, the image data amount is little by arc length sampling methods such as simulating, verifyings, it is discrete that the sample mode of the arc length such as present embodiment employing carries out the encryption of data.
(5) to comprising P
0i, P
1iAt interior P
0iP
1iBetween discrete point adopt least square method to carry out ellipse fitting; And with P
0iP
1iSection is designated as respectively P with the discrete point at the whole story of fitted ellipse intersection
0 (i+1), P
1 (i+1)
The average error of (6) establishing the i time match is ε
i, and
D in the formula
jBe the distance of a j discrete point to the (5) step gained fitted ellipse, N is for participating in counting out of match;
(7) preset error ε=0.005, to error ε
iCompare with predefined error ε; If ε
i≤ ε calculates and finishes.Otherwise, work as ε
iDuring ε, be positioned at P after casting out the i time iteration
0 (i+1)P
0iAnd P
1 (i+1)P
1iBetween point on the non-elliptic arc; Make i=i+1, to P
0iP
1iBetween the match of discrete point repeating step (5), until satisfy ε
iDuring≤ε, iteration finishes.
Table 1 is the oval leading edge process of reconstruction that this example is specifically controlled based on error, and iteration seven times satisfies the error requirements of setting, the average error of match this moment: ε
7=0.005mm, oval long and short semiaxis is respectively: R
a=1.011mm, R
b=0.604mm.
Table 1
According to embodiment as can be known, the present invention carries out iterative fitting by adopting the method for error control to discrete point, and casts out one by one the point on the non-elliptic arc in the selection area.In actual applications, the user only needs by setup control error ε process of reconstruction to be controlled, and can obtain required reconstructed results.
Claims (2)
1. oval front and rear edge method for reconstructing of the aerial blade based on error control is characterized in that may further comprise the steps:
(1) reads in the discrete measurement data of molded line of aerial blade, the selected region point that covers blade inlet edge or trailing edge;
(2) region point is sorted;
(3) region point after the ordering is carried out Cubic Spline Fitting;
(4) to the SPL of institute's match carry out etc. arc length encrypt discrete, and with the whole story discrete point be designated as respectively P
0i, P
1i
(5) with least square method to P
0iP
1iThe section discrete point carries out match, and with P
0iP
1iSection and fitted ellipse intersection the whole story discrete point be designated as respectively P
0 (i+1), P
1 (i+1)Ellipse fitting method is as follows:
1. as quafric curve f (x, y)=ax
2+ 2bxy+cy
2+ 2dx+2fy+g=0 the b that satisfies condition
2-ac<0 o'clock, quafric curve f (x, y) is expressed as ellipse; In the formula, a, b, c, d, f, g are the coefficient of quadratic polynomial equation;
Be provided with n data point P
i(x
i, y
i) participate in ellipse fitting, will be converted into based on the ellipse fitting of algebraic distance the least square problem of following belt restraining, the objective function of setting up ellipse fitting is as follows:
2. adopt least square method to try to achieve the coefficient a of quadratic polynomial equation, b, c, d, f behind the g, obtains fitted ellipse, and the elliptic geometry calculation of parameter is as follows:
The oval center of circle (x0, y0):
Transverse ra:
Ellipse short shaft r
b:
In the formula: j=1,2,3 ..., N, d
jBe the distance that j discrete point arrives gained fitted ellipse in the step (5), N is for participating in counting out of match;
(7) setting the iteration control error is ε: to ε
iCompare with ε, if ε
i≤ ε calculates and finishes; Otherwise, work as ε
iDuring ε, be positioned at P after casting out the i time iteration
0 (i+1)P
0iAnd P
1 (i+1)P
1iBetween point on the non-elliptic arc; Make i=i+1, to P
0iP
1iBetween the match of discrete point repeating step (5), until satisfy ε
iDuring≤ε, iteration finishes.
2. the oval front and rear edge method for reconstructing of the aerial blade based on error control according to claim 1 is characterized in that: described setting iteration control error ε=0.005.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310010830.8A CN103077315B (en) | 2013-01-11 | 2013-01-11 | Method for rebuilding front and back oval edges of aerial blade based on error control |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310010830.8A CN103077315B (en) | 2013-01-11 | 2013-01-11 | Method for rebuilding front and back oval edges of aerial blade based on error control |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103077315A true CN103077315A (en) | 2013-05-01 |
CN103077315B CN103077315B (en) | 2015-07-01 |
Family
ID=48153844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310010830.8A Expired - Fee Related CN103077315B (en) | 2013-01-11 | 2013-01-11 | Method for rebuilding front and back oval edges of aerial blade based on error control |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103077315B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103252536A (en) * | 2013-05-20 | 2013-08-21 | 西北工业大学 | Rolled blade flange head smooth transferring machining method |
CN105261065A (en) * | 2015-11-02 | 2016-01-20 | 南京航空航天大学 | Front and rear edge profile redesign method for adaptive processing |
CN107357959A (en) * | 2017-06-14 | 2017-11-17 | 西安交通大学 | A kind of high-precision blade profile characteristic parameter extraction method |
CN107526875A (en) * | 2017-07-31 | 2017-12-29 | 电子科技大学 | A kind of aerial blade type face mismachining tolerance method for visualizing |
CN109916534A (en) * | 2019-05-05 | 2019-06-21 | 国网黑龙江省电力有限公司哈尔滨供电公司 | A kind of temp measuring system and its temp measuring method of distribution type fiber-optic |
CN110909398A (en) * | 2019-11-18 | 2020-03-24 | 广东三维家信息科技有限公司 | Graph fitting method and device, electronic equipment and computer readable storage medium |
CN112967256A (en) * | 2021-03-09 | 2021-06-15 | 扬州大学 | Tunnel ovalization detection method based on spatial distribution |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102794488A (en) * | 2012-07-10 | 2012-11-28 | 上海交通大学 | Side milling processing method of resembled ruled surface integral wheel curved surfaces |
-
2013
- 2013-01-11 CN CN201310010830.8A patent/CN103077315B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102794488A (en) * | 2012-07-10 | 2012-11-28 | 上海交通大学 | Side milling processing method of resembled ruled surface integral wheel curved surfaces |
Non-Patent Citations (4)
Title |
---|
ANDREW FITZGIBBON 等: "《Direct Least Square Fitting of Ellipses》", 《TERN ANALYSIS AND MACHINE INTELLIGENCE》 * |
张力宁 等: "《叶片前缘高精度重建方法研究》", 《航空动力学报》 * |
杨密: "《复杂曲线轮廓度的误差评定》", 《中国优秀博硕士学位论文全文数据库 工程科技I辑》 * |
陆宏志 等: "《压气机叶片前缘形状的改进设计》", 《航空动力学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103252536B (en) * | 2013-05-20 | 2015-01-28 | 西北工业大学 | Rolled blade flange head smooth transferring machining method |
CN103252536A (en) * | 2013-05-20 | 2013-08-21 | 西北工业大学 | Rolled blade flange head smooth transferring machining method |
CN105261065B (en) * | 2015-11-02 | 2018-02-02 | 南京航空航天大学 | A kind of blade part front and rear edge profile towards adaptability processing redesigns method |
CN105261065A (en) * | 2015-11-02 | 2016-01-20 | 南京航空航天大学 | Front and rear edge profile redesign method for adaptive processing |
CN107357959B (en) * | 2017-06-14 | 2020-11-10 | 西安交通大学 | High-precision blade section characteristic parameter extraction method |
CN107357959A (en) * | 2017-06-14 | 2017-11-17 | 西安交通大学 | A kind of high-precision blade profile characteristic parameter extraction method |
CN107526875A (en) * | 2017-07-31 | 2017-12-29 | 电子科技大学 | A kind of aerial blade type face mismachining tolerance method for visualizing |
CN107526875B (en) * | 2017-07-31 | 2020-09-01 | 电子科技大学 | Visualization method for machining errors of aviation blade profile |
CN109916534A (en) * | 2019-05-05 | 2019-06-21 | 国网黑龙江省电力有限公司哈尔滨供电公司 | A kind of temp measuring system and its temp measuring method of distribution type fiber-optic |
CN109916534B (en) * | 2019-05-05 | 2021-08-13 | 国网黑龙江省电力有限公司哈尔滨供电公司 | Distributed optical fiber temperature measurement system and temperature measurement method thereof |
CN110909398A (en) * | 2019-11-18 | 2020-03-24 | 广东三维家信息科技有限公司 | Graph fitting method and device, electronic equipment and computer readable storage medium |
CN112967256A (en) * | 2021-03-09 | 2021-06-15 | 扬州大学 | Tunnel ovalization detection method based on spatial distribution |
CN112967256B (en) * | 2021-03-09 | 2023-11-24 | 扬州大学 | Tunnel ovalization detection method based on spatial distribution |
Also Published As
Publication number | Publication date |
---|---|
CN103077315B (en) | 2015-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103077315A (en) | Method for rebuilding front and back oval edges of aerial blade based on error control | |
CN101644915B (en) | Computer-control polishing method based on removal function prediction model | |
CN104493305B (en) | A kind of skiving processing method adjusted based on gear cutting machine | |
Boumali et al. | Comment on ‘Energy profile of the one-dimensional Klein–Gordon oscillator’ | |
CN103577635A (en) | Fitting method of blade surface data | |
CN103093091A (en) | Circular arc front-and-back edge of aviation blade reconstruction method based on error control | |
CN103676786A (en) | Curve fairing method based on accelerated speed principle | |
CN106991241A (en) | A kind of cutter chip pocket sharpening interference Forecasting Methodology | |
CN204095376U (en) | Pencil sharpener | |
CN206216002U (en) | A kind of aluminium alloy minute surface milling cutter | |
CN103317187A (en) | Groove cutting broach | |
Liu | A variant of Chebyshev inequality with applications | |
Michalski et al. | Surface topography of cylindrical gear wheels after smoothing in abrasive mass, honing and shot peening | |
CN205310086U (en) | Domestic cutter special abrasive disc of polishing | |
CN103464796A (en) | Complex molding comb tooth cutter | |
CN204107015U (en) | A kind of chessman | |
JP6376364B2 (en) | How to set the energy consumption of two tools when machining the pipe end | |
CN103692298A (en) | Angle adjusting device for grinding wheel on tool grinder | |
RU2663371C2 (en) | Method for processing blade blank nozzle apparatus of gas turbine engine | |
CN205254809U (en) | A governing valve for sand blasting equipment | |
CN103778652A (en) | Generation method for three-dimensional monster models | |
Höfener et al. | Deep learning nuclei detection | |
TH69208B (en) | Internal gear grinding method | |
TH161507A (en) | Internal gear grinding method | |
Carneiro et al. | A comparative thermophysiological study in tops for running |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150701 Termination date: 20160111 |
|
EXPY | Termination of patent right or utility model |