CN105869212B - A kind of ice shape simplification method - Google Patents
A kind of ice shape simplification method Download PDFInfo
- Publication number
- CN105869212B CN105869212B CN201610153144.XA CN201610153144A CN105869212B CN 105869212 B CN105869212 B CN 105869212B CN 201610153144 A CN201610153144 A CN 201610153144A CN 105869212 B CN105869212 B CN 105869212B
- Authority
- CN
- China
- Prior art keywords
- ice
- point
- shape
- ice shape
- curve
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
A kind of ice shape simplifies method and belongs to aviation field, the complicated irregular two-dimentional ice shape obtained in icing tunnel is simplified, using B-spline curves method, for the ice of two kinds of fundamental types of horn shape type ice and streamlined ice, the parametric description method based on ice shape characteristic point position is provided, and is disclosed and a kind of optimal or quasi-optimal scheme ice shape is sought as object function using ice shape similarity simplified method；The present invention provides technical support for fields such as aircraft icing research, icing wind tunnel debugging and the processing of icing tunnel experimental result, has preferable engineering application value.
Description
Technical field
The present invention relates to aviation fields, provide a kind of ice shape simplification method, can be applied to aircraft icing research, icing wind tunnel
Debugging and the processing of icing tunnel experimental result.
Background technology
Aircraft may meet with icing phenomenon in the cloud layer containing super-cooling waterdrop when flight, icing is the main of flight safety
One of hidden danger is then gently that the safe flight range of aircraft reduces, and major accident that is heavy then leading to fatal crass, ice shape is icing pair
Aircraft safety influences the input condition of research and foundation and the basis of anti-deicing research.
Icing wind tunnel experiment is to obtain one of the main means of ice shape, and the icing that icing wind tunnel is tested is complex-shaped more
Sample, some shape irregulars, especially in the case where temperature is higher, Liquid water content is larger, there are liquid in freezing process
State water frozen surface overflow and it is secondary freeze, pinniform region can be formed near icing critical, the shape of this ice is extremely
Complexity will investigate influence of this ice to aerodynamic characteristic, it is necessary to make ice model in conventional wind-tunnel with the method for experiment.Though
So modern Computerized Numerical Control processing technology is very advanced, but still is unable to processing so complicated and being easily mounted on object table for entirely accurate
The ice model in face.It must be set up rational ice shape thus and simplify method so that treated ice graphic data can both keep ice
Main feature, and the making of ice model is facilitated to process.
In the world, icing tunnel is all had in the aviations developed country such as the U.S., Europe, and knot is yet there are no from the data that can be inquired
Icing tunnel ice shape simplifies method.The country starts late in aircraft icing area research, is used for the large-scale icing wind tunnel of productive experiment
And just build up, although having several small-sized icing wind tunnels, mating test on Mianyang, Shenyang, Wuhan, Harbin and other places
Equipment, experimental technique and theory show that icing wind tunnel ice shape simplifies grinding for method at present it is still necessary to perfect, from the open data consulted
Study carefully and still belongs to blank at home.
Invention content
The purpose of the invention provides a kind of letter of the complicated ice shape of experiment aiming at above-mentioned existing technical problem
Change method, the area researches such as research, icing wind tunnel debugging and the processing of icing tunnel experimental result that can freeze for aircraft provide technology branch
Support has preferable engineering application value.
What the purpose of the invention was realized in：
(1) type of ice shape is determined
Determine that ice shape is typical streamlined ice or horn shape type ice by observation；
(2) characteristic point of ice shape is determined
1. five characteristic points of horn shape type ice are upper ice angle point, lower ice angle point, stationary point, limes superiors point and limit inferior point are adopted
It obtains with the following method：
Two-dimentional ice shape is divide into upper part and lower part by the first step using aerofoil profile edge as line of demarcation；
Vertically ice shape is searched for from object plane farthest point in surface to second step respectively, this 2 points are upper and lower ice angle point, are searched for
To ice shape, each elementary section carries out interpolated value in journey, calculates separately each point with a distance from object plane；
Point between the upper and lower ice angle point of three-wave mixing from object plane most nearby, which is stationary point；
4th step is scanned for according to band ice aerofoil profile data and initial aerofoil profile data, and misaligned part is the areas Bing Xing
Domain, the upper and lower critical point in ice shape region are ice shape Upper-lower Limit point；
2. three characteristic points of streamlined ice are stationary point, limes superiors point and limit inferior point, obtain with the following method：
The first step searches for point of the ice shape from object plane farthest, which is stationary point；
Second step is scanned for according to band ice aerofoil profile data and initial aerofoil profile data, and misaligned part is the areas Bing Xing
Domain, the upper and lower critical point in ice shape region are ice shape Upper-lower Limit point；
(3) the ice shape of feature based point position simplifies process
It is in the case where maintaining original ice shape characteristic point position constant, using passing through each characteristic point that ice shape, which simplifies process,
B-spline curves replace the parametric description method of original ice shape curve, and process is：
1. ice shape is described using B-spline Curve, B-spline curves are made of sectional curve, 4 in B-spline Curve
A control point controls one section of curve, then the hop count of the B-spline curves at a N number of control point is (N-3) section；At Upper-lower Limit point
Control point take triple vertex, i.e., all corresponding X, Y coordinates for being taken as the limit point of X, Y coordinates value at the control point of three coincidences
Value；Stationary point, upper and lower apex control point point-blank should take the vertical direction on the vertex and aerofoil profile range line, three controls
The distance between system point is a small amount of for one, and the coordinate at intermediate control point is taken as corresponding apex coordinate；Between two vertex
Ice shape curve negotiating one or more snippets curve composition, corresponding control point can be adjusted according to ice graphic data；
Keep the ice shape curve after simplifying similar to original ice shape curve 2. adjusting the parameters in above-mentioned B-spline curves
Degree is optimal, and process seeks optimal or quasi-optimal scheme optimization as one using ice shape similarity as object function
Journey, finding method use Gene hepatitis B vaccine or other optimization algorithms；The ice shape similarity is by calculating ice shape variance rate
It obtains, the ice shape variance rate is two ice shape areas not the ratio between the sum of overlapping area and two ice shape average areas, ice shape
The sum of the variance rate of similarity and ice shape be 100%；
3. the B-spline curves after above-mentioned optimization are the ice shape curve after simplifying.
The present invention provides a kind of simplification method of ice shape, makes ice by the characteristic feature of the two-dimentional ice shape of research aircraft icing
Mould processing is more convenient, simple, is ground for fields such as aircraft icing research, icing wind tunnel debugging and the processing of icing tunnel experimental result
Offer technical support is provided, there is preferable engineering application value.
Description of the drawings
Fig. 1 typical case's horn shape type ice
The typical streamlined ice of Fig. 2
The original ice shapes of Fig. 3
Fig. 4 simplifies front and back ice shape and control point
Specific implementation mode
The embodiment of the invention is described in detail below in conjunction with the accompanying drawings：A kind of ice shape simplification method, it is special
Sign is to include the following steps：
(1) type of ice shape is determined
Determine that ice shape is typical streamlined ice or horn shape type ice by observation；
(2) characteristic point of ice shape is determined
1. five characteristic points of horn shape type ice are upper ice angle point, lower ice angle point, stationary point, limes superiors point and limit inferior point are adopted
It obtains with the following method：
Two-dimentional ice shape is divide into upper part and lower part by the first step using aerofoil profile edge as line of demarcation；
Vertically ice shape is searched for from object plane farthest point in surface to second step respectively, this 2 points are upper and lower ice angle point, are searched for
To ice shape, each elementary section carries out interpolated value in journey, calculates separately each point with a distance from object plane；
Point between the upper and lower ice angle point of three-wave mixing from object plane most nearby, which is stationary point；
4th step is scanned for according to band ice aerofoil profile data and initial aerofoil profile data, and misaligned part is the areas Bing Xing
Domain, the upper and lower critical point in ice shape region are ice shape Upper-lower Limit point；
2. three characteristic points of streamlined ice are stationary point, limes superiors point and limit inferior point, obtain with the following method：
The first step searches for point of the ice shape from object plane farthest, which is stationary point；
Second step is scanned for according to band ice aerofoil profile data and initial aerofoil profile data, and misaligned part is the areas Bing Xing
Domain, the upper and lower critical point in ice shape region are ice shape Upper-lower Limit point；
(3) the ice shape of feature based point position simplifies process
It is in the case where maintaining original ice shape characteristic point position constant, using passing through each characteristic point that ice shape, which simplifies process,
B-spline curves replace the parametric description method of original ice shape curve, and process is：
1. ice shape is described using B-spline Curve, B-spline curves are made of sectional curve, 4 in B-spline Curve
A control point controls one section of curve, then the hop count of the B-spline curves at a N number of control point is (N-3) section；At Upper-lower Limit point
Control point take triple vertex, i.e., all corresponding X, Y coordinates for being taken as the limit point of X, Y coordinates value at the control point of three coincidences
Value；Stationary point, upper and lower apex control point point-blank should take the vertical direction on the vertex and aerofoil profile range line, three controls
The distance between system point is a small amount of for one, and the coordinate at intermediate control point is taken as corresponding apex coordinate；Between two vertex
Ice shape curve negotiating one or more snippets curve composition, corresponding control point can be adjusted according to ice graphic data；
Keep the ice shape curve after simplifying similar to original ice shape curve 2. adjusting the parameters in above-mentioned B-spline curves
Degree is optimal, and process seeks optimal or quasi-optimal scheme optimization as one using ice shape similarity as object function
Journey, finding method use Gene hepatitis B vaccine or other optimization algorithms；The ice shape similarity is by calculating ice shape variance rate
It obtains, the ice shape variance rate is two ice shape areas not the ratio between the sum of overlapping area and two ice shape average areas, ice shape
The sum of the variance rate of similarity and ice shape be 100%；
3. the B-spline curves after above-mentioned optimization are the ice shape curve after simplifying.
Claims (1)
1. a kind of ice shape simplifies method, it is characterised in that include the following steps：
(1) type of ice shape is determined
Determine that ice shape is typical streamlined ice or horn shape type ice by observation；
(2) characteristic point of ice shape is determined
1. five characteristic points of horn shape type ice are upper ice angle point, lower ice angle point, stationary point, limes superiors point and limit inferior point, using such as
Lower method obtains：
Two-dimentional ice shape is divide into upper part and lower part by the first step using aerofoil profile edge as line of demarcation；
Vertically ice shape is searched for from object plane farthest point in surface to second step respectively, this 2 points are upper and lower ice angle point, in search process
To ice shape, each elementary section carries out interpolated value, calculates separately each point with a distance from object plane；
Point between the upper and lower ice angle point of three-wave mixing from object plane most nearby, which is stationary point；
4th step is scanned for according to band ice aerofoil profile data and initial aerofoil profile data, and misaligned part is ice shape region, ice
The upper and lower critical point in shape region is ice shape Upper-lower Limit point；
2. three characteristic points of streamlined ice are stationary point, limes superiors point and limit inferior point, obtain with the following method：
The first step searches for point of the ice shape from object plane farthest, which is stationary point；
Second step is scanned for according to band ice aerofoil profile data and initial aerofoil profile data, and misaligned part is ice shape region, ice
The upper and lower critical point in shape region is ice shape Upper-lower Limit point；
(3) the ice shape of feature based point position simplifies process
It is in the case where maintaining original ice shape characteristic point position constant, using the B samples by each characteristic point that ice shape, which simplifies process,
Curve replaces the parametric description method of original ice shape curve, process to be：
1. ice shape is described using B-spline Curve, B-spline curves are made of sectional curve, 4 controls in B-spline Curve
System point one section of curve of control, then the hop count of the B-spline curves at a N number of control point is (N-3) section；Control at Upper-lower Limit point
System point takes triple vertex, i.e., all corresponding X, Y coordinates value for being taken as the limit point of X, Y coordinates value at the control point of three coincidences；It stays
Point, limes superiors point, the control points limit inferior Dian Chu point-blank should take the vertical direction on the vertex and aerofoil profile range line, and three
The distance between a control point is one a small amount of, and the coordinate at intermediate control point is taken as corresponding apex coordinate；Two characteristic points
Between ice shape curve negotiating one or more snippets curve composition, corresponding control point can be adjusted according to ice graphic data；
The ice shape curve after simplifying and the similarity of original ice shape curve is set to reach 2. adjusting the parameters in above-mentioned B-spline curves
To optimal, process seeks optimal or quasi-optimal scheme optimization process as one using ice shape similarity as object function,
Finding method uses Gene hepatitis B vaccine or other optimization algorithms；The ice shape similarity is obtained by calculating ice shape variance rate
It arrives, the ice shape variance rate is two ice shape areas not the ratio between the sum of overlapping area and two ice shape average areas, ice shape
The sum of similarity and the variance rate of ice shape are 100%；
3. the B-spline curves after above-mentioned optimization are the ice shape curve after simplifying.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610153144.XA CN105869212B (en) | 2016-03-11 | 2016-03-11 | A kind of ice shape simplification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610153144.XA CN105869212B (en) | 2016-03-11 | 2016-03-11 | A kind of ice shape simplification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105869212A CN105869212A (en) | 2016-08-17 |
| CN105869212B true CN105869212B (en) | 2018-10-12 |
Family
ID=56625289
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610153144.XA Expired - Fee Related CN105869212B (en) | 2016-03-11 | 2016-03-11 | A kind of ice shape simplification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105869212B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109899248B (en) * | 2019-03-04 | 2020-11-03 | 天津工业大学 | Horizontal axis wind turbine blade frost ice shape forming method based on polynomial fitting |
| CN110889233B (en) * | 2019-12-04 | 2023-03-28 | 中国空气动力研究与发展中心 | Icing wind tunnel test ice-shaped three-dimensional entity generation method and three-dimensional entity ice |
| CN111368367B (en) * | 2020-02-25 | 2023-08-25 | 空气动力学国家重点实验室 | Parameterization method for icing shape of icing wind tunnel icing test |
| CN111323200B (en) * | 2020-05-11 | 2020-08-07 | 中国空气动力研究与发展中心低速空气动力研究所 | Icing area calculation method for icing wind tunnel test |
| CN111539072B (en) * | 2020-07-08 | 2020-09-29 | 中国空气动力研究与发展中心低速空气动力研究所 | Quantitative evaluation method and system for ice-shaped ice angle in icing wind tunnel test |
| CN111739048B (en) * | 2020-07-27 | 2020-11-24 | 中国空气动力研究与发展中心低速空气动力研究所 | Icing wind tunnel test ice-shaped geometric contour line digitalization method |
| CN112985753B (en) * | 2021-05-06 | 2021-10-19 | 中国空气动力研究与发展中心计算空气动力研究所 | Ice shape parameterization method for wind tunnel test |
| CN113112568B (en) * | 2021-06-08 | 2021-09-07 | 中国空气动力研究与发展中心低速空气动力研究所 | Two-dimensional ice shape simplification method based on arc fitting |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101865668A (en) * | 2010-04-29 | 2010-10-20 | 北京航空航天大学 | Three-dimensional ice form detection instrument |
| CN101894370A (en) * | 2010-07-14 | 2010-11-24 | 苏州大学 | Automatic generation method of shape parameter-adaptive oracle-bone inscription contour glyphs |
| CN102494635A (en) * | 2011-12-13 | 2012-06-13 | 中国空气动力研究与发展中心设备设计及测试技术研究所 | Wind tunnel model ice shape measuring method and device |
| CN103057691A (en) * | 2011-09-06 | 2013-04-24 | 空中客车西班牙运营有限责任公司 | Aircraft tail surface with leading edge section of undulated shape |
| CN103387055A (en) * | 2012-05-08 | 2013-11-13 | 波音公司 | Ice protection for aircraft using electroactive polymer surfaces |
| CN103569366A (en) * | 2012-08-02 | 2014-02-12 | 罗斯蒙特航天公司 | Rotor ice protection systems and methods |
| CN104802997A (en) * | 2013-10-10 | 2015-07-29 | 波音公司 | Methods and apparatus for detecting ice formation on aircraft |
-
2016
- 2016-03-11 CN CN201610153144.XA patent/CN105869212B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101865668A (en) * | 2010-04-29 | 2010-10-20 | 北京航空航天大学 | Three-dimensional ice form detection instrument |
| CN101894370A (en) * | 2010-07-14 | 2010-11-24 | 苏州大学 | Automatic generation method of shape parameter-adaptive oracle-bone inscription contour glyphs |
| CN103057691A (en) * | 2011-09-06 | 2013-04-24 | 空中客车西班牙运营有限责任公司 | Aircraft tail surface with leading edge section of undulated shape |
| CN102494635A (en) * | 2011-12-13 | 2012-06-13 | 中国空气动力研究与发展中心设备设计及测试技术研究所 | Wind tunnel model ice shape measuring method and device |
| CN103387055A (en) * | 2012-05-08 | 2013-11-13 | 波音公司 | Ice protection for aircraft using electroactive polymer surfaces |
| CN103569366A (en) * | 2012-08-02 | 2014-02-12 | 罗斯蒙特航天公司 | Rotor ice protection systems and methods |
| CN104802997A (en) * | 2013-10-10 | 2015-07-29 | 波音公司 | Methods and apparatus for detecting ice formation on aircraft |
Non-Patent Citations (4)
| Title |
|---|
| 民用飞机翼面冰形的工程修正;赵宾宾 等;《民用飞机设计与研究》;20131231(第2013年第4期);第52-55页 * |
| 运12F飞机结冰风洞试验的参数模拟;王强;《探索 创新 交流——第五届中国航空学会青年科技论坛文集(第5集)》;20121008;第102-108页 * |
| 飞机机翼除冰过程的数值研究;马成樑;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20080115(第2008年第01期);第C031-39页 * |
| 飞机结冰冰形预测的建模与仿真;潘环 等;《系统仿真学报》;20140131;第26卷(第1期);第221-229页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105869212A (en) | 2016-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105869212B (en) | A kind of ice shape simplification method | |
| CN106772387A (en) | A kind of wind shear recognition methods | |
| CN104636750B (en) | A kind of pavement crack recognizer and system based on double scale clustering algorithms | |
| CN107942411B (en) | Atmospheric visibility prediction method | |
| CN106364697B (en) | The multistage variant gliding Waverider aircraft design method in the fast domain of width that theory is led based on cone | |
| CN108460217A (en) | A kind of unstable state three-dimensional icing method for numerical simulation | |
| CN106650095B (en) | Correction method of unmanned aerial vehicle control matrix based on wind tunnel test data and CFD calculation | |
| CN105957055B (en) | A kind of ice shape similarity quantitative estimation method | |
| CN109448087A (en) | A kind of shield tunnel surface point cloud data segmentation method | |
| CN103675923B (en) | Foreign matter detecting method and equipment for detecting foreign matter based on static clutter map | |
| CN103927537A (en) | Extraction algorithm of self-adapting intersection regions of strokes of Chinese characters | |
| CN109408941B (en) | Aircraft pneumatic optimization method based on data mining and genetic algorithm | |
| CN103839274B (en) | Extended target tracking method based on geometric proportion relation | |
| CN110188378B (en) | Pneumatic data fusion method based on neural network | |
| CN109583131A (en) | A kind of optimum design method of the surface micro-structure dimensional parameters with drag-reduction effect | |
| Bland | AGARD two-dimensional aeroelastic configurations | |
| CN112926132A (en) | Fixed wing airfoil aerodynamic shape design method considering influence of three-dimensional effect | |
| CN112985753B (en) | Ice shape parameterization method for wind tunnel test | |
| Woodward | USSAERO computer program development, versions B and C | |
| CN106381831A (en) | Method for determining aeration facility vent hole sectional area | |
| CN109000682A (en) | The method of aerial coarse alignment based on rectilinear path | |
| CN109325257B (en) | Modeling method for rapidly generating flight load analysis model | |
| CN106023064A (en) | Assessment method for field of vision of airplane cabin | |
| CN111666366A (en) | Flight monitoring data display method and system | |
| Bisset et al. | Aspects of turbulent/non-turbulent interfaces |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181012 Termination date: 20190311 |