CN111368367A - Parameterization method for ice shape of icing wind tunnel icing test - Google Patents
Parameterization method for ice shape of icing wind tunnel icing test Download PDFInfo
- Publication number
- CN111368367A CN111368367A CN202010116493.0A CN202010116493A CN111368367A CN 111368367 A CN111368367 A CN 111368367A CN 202010116493 A CN202010116493 A CN 202010116493A CN 111368367 A CN111368367 A CN 111368367A
- Authority
- CN
- China
- Prior art keywords
- curve
- icing
- numerical
- deicing
- ice
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention discloses a parameterization method for an ice shape of an icing wind tunnel icing test, which comprises the following steps of: selecting characteristic points: calculating the curvature change of the numerical ice-shaped curve, and selecting a characteristic point on the ice-shaped curve according to a curvature change rule; a curve parameterization step: and carrying out regression interpolation on the characteristic points on the numerical deicing curve by using a regression model to obtain a smooth parametric deicing curve. Aiming at irregular complex ice-shaped curves, the mechanism of filtering noise by using a regression model is utilized, and characteristic points are specified by artificially or setting a selection rule, so that the aim of quickly parameterizing the ice-shaped curves can be fulfilled under the condition of keeping the original ice-shaped characteristics.
Description
Technical Field
The invention relates to a parameterization method, in particular to a parameterization method for an icing test ice shape of an icing wind tunnel.
Background
During the flight of an airplane, water droplets in a cloud layer can be condensed at a certain height, and when the temperature is zero degrees centigrade and below, ice can be formed on the surface of the airplane and can further form accumulated ice. Icing of an aircraft is one of the main hidden dangers threatening the flight safety of the aircraft, and has received wide attention in recent years.
At present, the research means of airplane icing includes flight test, wind tunnel test and numerical simulation, wherein the wind tunnel test becomes an important means for the research of icing and deicing by virtue of the icing reproducing capability with high reliability. The icing shape is one of the important attention results of the icing wind tunnel test, and can be used for evaluating the aerodynamic characteristics of the frozen airplane. During the simulation evaluation of the pneumatic characteristic numerical value, the complexity of the ice shape directly influences the grid generation quality, thereby influencing the difficulty and precision of calculation; the complexity of the ice shape can cause great difficulty in the formation of the ice shape when the aerodynamic characteristic test is evaluated. Most of the time, the ice shape obtained by the ice wind tunnel is extremely irregular, and more burrs exist, so that the ice shape is parameterized in a later period, the parameterized ice shape can keep the characteristics of the original ice shape, and the original ice shape can be simplified to a certain extent. Therefore, a parameterization method of the ice shape obtained by an icing wind tunnel icing test needs to be developed to enhance the operability and scientificity of the later use of the ice shape.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a parameterization method for an icing wind tunnel icing test ice shape, aiming at a two-dimensional ice shape curve, a regression Kriging model is applied to filter noise on the ice shape curve to obtain a final parameterized ice shape curve, so that the aim of correcting and optimizing the ice shape is fulfilled.
The purpose of the invention is realized by the following technical scheme: a parameterization method for an icing shape of an icing wind tunnel icing test comprises the following steps:
selecting characteristic points: calculating the curvature change of the numerical ice-shaped curve, and selecting a characteristic point on the ice-shaped curve according to a curvature change rule;
a curve parameterization step: and carrying out regression interpolation on the characteristic points on the numerical deicing curve by using a regression model to obtain a smooth parametric deicing curve.
Before the step of selecting the characteristic points, a step of acquiring a numerical deicing curve is required to be completed.
The step of acquiring the numerical icing curve comprises the step of acquiring a two-dimensional numerical data point set { (x) of the test icing in the wind tunnel icing test by adopting a hot knife methodi,yi) I |, 1, 2, …, N }, and an ice-shaped curve is obtained after two-dimensional numerical representation.
The calculating the curvature change of the numerical ice curve and selecting the characteristic points on the numerical ice curve according to the curvature change rule specifically comprise:
calculating the curvatures of all data points on the numerical deicing curve and the curvature change between adjacent data points;
and m characteristic points are selected on the numerical deicing curve according to the magnitude relation of curvature change to obtain a coordinate point set { lambda } of all the characteristic pointsj(ξj,ηj)|j=1,2,…,m}。
The method for selecting m characteristic points on the numerical deicing curve according to the size relation of curvature change comprises the following steps:
selecting a characteristic points at the positions of which the curvature change is smaller than a given threshold value;
b characteristic points are selected at all positions with curvature change larger than a given threshold value.
And the relationship between a, b and m is a + b.
The method for carrying out regression interpolation on the characteristic points on the numerical deicing curve by using the regression model to obtain the smooth parametric deicing curve comprises the following steps:
taking the numerical deicing curve as an input object, and performing regression interpolation by using a regression Kriging model to obtain an interpolation model;
and outputting interpolation results of the positions of all the characteristic points to obtain a smooth parametric deicing curve with noise points on the numerical deicing curve filtered.
Taking the numerical deicing curve as an input object, and performing regression interpolation by using a regression Kriging model to obtain an interpolation model, which specifically comprises the following steps:
the data point set of the numerical icing curve { (x)i,yi) 1, 2, …, N into a polar coordinate data point set { (r)i,θi)|i=1,2,…,N};
Using a set of polar data points { (r)i,θi) And (3) performing regression Kriging model on the | i ═ 1, 2, …, N to obtain an interpolation model Y (X).
Outputting interpolation results of positions of all the characteristic points to obtain a smooth parametric deicing curve with noise points on the numerical deicing curve filtered, wherein the smooth parametric deicing curve specifically comprises the following steps:
set of coordinate points of m feature points { lambdaj(ξj,ηj) 1, 2, …, m into a corresponding set of polar points kj(ρj,εj)|j=1,2,…,m};
Each polar coordinate ρjSubstitution into interpolation model Y (X) yields a new set of polar coordinate points { κ'j(ρj,ε′j)|j=1,2,…,m};
Converting the coordinate points into a coordinate point set { (x) 'in a Cartesian coordinate system'j,y′j) And j equals to 1, 2, …, m, and a parametric deicing curve for filtering noise is obtained.
The invention has the beneficial effects that: a parameterization method for an icing test ice shape of an icing wind tunnel aims at an irregular complex ice shape curve, utilizes a mechanism of noise filtering of a regression model, and can achieve the purpose of rapid parameterization of the ice shape curve under the condition of keeping original ice shape characteristics by specifying characteristic points through artificial or set selection rules.
Drawings
FIG. 1 is a schematic illustration of the ice parameterization process of the present invention;
FIG. 2 is a schematic diagram of the transformation of coordinates of an ice-shaped curve;
FIG. 3 is a graph showing the comparison between the ice-shaped curve and the parameterized curve
In the figure, 1-an ice-shaped profile curve, 1-1-an ice-shaped profile curve under a polar coordinate system, 2-icing of the leading edge of the blade, 3-icing of the leading edge of the blade, 4-characteristic points and 5-parametric icing of the ice-shaped curve.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. 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.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
In the description of the present invention, it should be noted that the terms "upper", "inner", "outer", etc. indicate orientations or positional relationships based on those shown in the drawings or orientations or positional relationships that the products of the present invention conventionally use, which are merely for convenience of description and simplification of description, but do not indicate or imply that the devices or elements referred to must have a specific orientation, be constructed in a specific orientation, and be operated, and thus, should not be construed as limiting the present invention.
In the description of the present invention, it should also be noted that, unless otherwise explicitly specified or limited, the terms "disposed," "mounted," and "connected" are to be construed broadly, e.g., as meaning fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
The technical solutions of the present invention are further described in detail below with reference to the accompanying drawings, but the scope of the present invention is not limited to the following.
As shown in fig. 1, a parameterization method for an icing test ice shape of an icing wind tunnel includes the following steps:
s1, acquiring a numerical deicing curve: blade leading edge icing 2 is formed on the blade leading edge 3, and a two-dimensional numerical data point set { (x) of a test ice shape in a wind tunnel icing test is obtained by adopting a hot knife methodi,yi) 1, 2, …, N }, and obtaining a two-dimensional numerical ice-shaped profile curve 1;
s2, selecting characteristic points: calculating the curvature change of the numerical ice-shaped profile curve 1, and selecting a characteristic point 4 on the ice-shaped profile curve 1 according to the curvature change rule;
s3, curve parameterization step: and (3) carrying out regression interpolation on the characteristic points on the numerical deicing contour curve 1 by using a regression model to obtain a smooth parametric deicing curve 5.
Further, the calculating the curvature change of the numerical icing shape profile curve 1 and selecting the feature point 4 on the numerical icing shape profile curve 1 according to the curvature change rule specifically include:
s21, calculating the curvatures of all data points on the numerical deicing contour curve 1 and the curvature change between adjacent data points;
s22, selecting m characteristic points on the numerical deicing contour curve 1 according to the magnitude relation of curvature change to obtain a coordinate point set { lambda } of all the characteristic pointsj(ξj,ηj)|j=1,2,…,m}。
Further, the selection of m feature points on the numerical icing profile curve 1 according to the magnitude relation of the curvature change includes the following steps:
selecting a characteristic points at the positions of which the curvature change is smaller than a given threshold value;
b characteristic points are selected at all positions with curvature change larger than a given threshold value.
And the relationship between a, b and m is a + b.
The position where the curvature change is smaller than the given threshold value represents that the curve is relatively smooth and has fewer burrs, and the position where the curvature change is larger than the given threshold value represents that the curve has larger curvature or is more prominent and has more burrs; more characteristic points are selected at positions with larger curvature change, so that the geometric characteristics of the icing curve can be retained to the greatest extent; in fig. 1, a total of 15 feature points are selected on the ice-shaped profile curve 1, 8 feature points are selected on the middle projection portion (the position where the curvature change is large) of the ice-shaped profile curve 1, and the remaining feature points are selected on both sides (the position where the curvature change is small).
Further, the obtaining of the smooth parametric icing curve by performing regression interpolation on the characteristic points on the numerical icing curve by using the regression model includes the following steps:
s31, taking the numerical deicing curve as an input object, and performing regression interpolation by using a regression Kriging model to obtain an interpolation model;
and S32, outputting interpolation results of all the feature point positions to obtain a smooth parametric deicing curve with noise points on the numerical deicing curve filtered.
Further, taking the numerical deicing curve as an input object, and performing regression interpolation by using a regression Kriging model to obtain an interpolation model specifically as follows:
the data point set of the numerical icing curve { (x)i,yi) 1, 2, …, N into a polar coordinate data point set { (r)i,θi)|i=1,2,…,N};
Using a set of polar data points { (r)i,θi) Carrying out regression Kriging model interpolation on | i ═ 1, 2, …, N } to obtainInterpolation model y (x).
wherein f isj(X) is a basis function, here chosen as constant 1, βjIs the coefficient corresponding to the basis function,a mathematical expectation value representing y (x); z (X) is mean zero and varianceThe static random process of (a); for any two positions X and X' in the design space, the covariance of the random quantity is:
wherein, the correlation function R (X, X') is only related to the space distance and represents the correlation between random variables at different positions, when the distance between two positions in the correlation function is infinite, R is 0, and when the distance between two positions is zero, R is 1; r decreases with increasing distance. Wherein, the correlation matrix and the correlation vector are defined as:
R(R(X(i),X(j))i,j
rx=[R(X(1),X),R(X(2),X),…,R(X(n),X)]T
the correlation model satisfies the following conditions: 1. when the sample point X(i)And X(j)When the distance tends to 0, the corresponding correlation function value tends to 1; 2. when the sample point distance tends to infinity, the correlation function value tends to 0; 3. as the sample point distance increases, the correlation function value decreases smoothly; 4. at least one order of conductability is satisfied. Selecting a gaussian function as a correlation model:
wherein theta is a Kriging hyper-parameter and can be obtained by utilizing a maximum likelihood estimation method; for the regression Kriging model needing to filter noise, a constant needs to be added to the diagonal line of the correlation matrix to improve the positive nature of the matrix:
R′=R+γI
wherein, I is a unit diagonal matrix, γ is a constant, the larger the value thereof is, the stronger the noise filtering effect is, further, the value range is 0-10, preferably, the value is 10, and can be adjusted according to the characteristics of the ice-shaped contour curve 1.
Outputting interpolation results of positions of all the characteristic points to obtain a smooth parametric deicing curve with noise points on the numerical deicing curve filtered, wherein the smooth parametric deicing curve specifically comprises the following steps:
set of coordinate points of m feature points { lambdaj(ξj,ηj) 1, 2, …, m into a corresponding set of polar points { κ |j(ρj,εj)|j=1,2,…,m};
Each polar coordinate ρjSubstitution into interpolation model Y (X) yields a new set of polar points, { kappa', in the interpolation modelj(ρj,ε'j)|j=1,2,…,m};
Converting the coordinate points into a coordinate point set { (x) 'in a Cartesian coordinate system'j,y′j) And j equals to 1, 2, …, m, and a parametric deicing curve for filtering noise is obtained.
Further, as shown in fig. 1 and fig. 2, the two-dimensional ice-shaped profile curve 1 obtained by one test is described by 123 discrete data points, and one x coordinate of the data point set in the rectangular coordinate system of the two-dimensional ice-shaped profile curve 1 corresponds to at least two y coordinates, which is inconvenient for performing regression interpolation, so that the data point set is converted into the polar coordinate system to ensure the one-to-one correspondence of the coordinates, and the ice-shaped profile curve 1-1 in the polar coordinate system is obtained.
The ice-shaped profile curve 1-1 under the polar coordinate is also described by 123 discrete data points, regression Kriging model interpolation is carried out by using the 123 discrete data points, a positive fixed parameter gamma is adjusted, and the best effect is obtained when the gamma is 10 through test verification, so that the parameterized ice-shaped curve 5 which can reflect the characteristics of the ice-shaped profile curve 1 and has smooth curvature is obtained.
As shown in fig. 3, a lot of burrs exist on the icing test ice-shaped profile curve 1, that is, "noise" exists in the curve data, and through regression Kriging interpolation, the outputted parameterized ice-shaped curve 5 well retains the main characteristics of the icing test ice-shaped profile curve 1, successfully filters the noise of the curve data, and obtains the parameterized ice-shaped curve 5 with smooth curvature.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes performed by the present specification and drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.
Claims (8)
1. A parameterization method for an icing shape of an icing wind tunnel icing test is characterized by comprising the following steps: the method comprises the following steps:
selecting characteristic points: calculating the curvature change of the numerical ice-shaped curve, and selecting a characteristic point on the ice-shaped curve according to a curvature change rule;
a curve parameterization step: and carrying out regression interpolation on the characteristic points on the numerical deicing curve by using a regression model to obtain a smooth parametric deicing curve.
2. The parameterization method for the ice shape of the icing wind tunnel icing test according to claim 1, characterized by comprising the following steps of: before the step of selecting the characteristic points, a step of acquiring a numerical deicing curve is required to be completed.
3. The parameterization method for the ice shape of the icing wind tunnel icing test according to claim 2, characterized by comprising the following steps of: the step of obtaining the numerical icing curve comprises the step of obtaining a two-dimensional numerical data point set { (x) of test icing in a wind tunnel icing test by adopting a hot knife methodi,yi) I |, 1, 2, …, N }, and an ice-shaped curve is obtained after two-dimensional numerical representation.
4. The parameterization method for the ice shape of the icing wind tunnel icing test according to claim 1, characterized by comprising the following steps of: the calculating the curvature change of the numerical ice curve and selecting the characteristic points on the numerical ice curve according to the curvature change rule specifically comprise:
calculating the curvatures of all data points on the numerical deicing curve and the curvature change between adjacent data points;
and m characteristic points are selected on the numerical deicing curve according to the magnitude relation of curvature change to obtain a coordinate point set { lambda } of all the characteristic pointsj(ξj,ηj)|j=1,2,…,m}。
5. The parameterization method for the ice shape of the icing wind tunnel icing test according to claim 4, characterized by comprising the following steps of: the method for selecting m characteristic points on the numerical deicing curve according to the size relation of curvature change comprises the following steps:
selecting a characteristic points at the positions of which the curvature change is smaller than a given threshold value;
b characteristic points are selected at all positions with curvature change larger than a given threshold value.
6. The parameterization method for the ice shape of the icing wind tunnel icing test according to claim 1, characterized by comprising the following steps of: the method for carrying out regression interpolation on the characteristic points on the numerical deicing curve by using the regression model to obtain the smooth parametric deicing curve comprises the following steps:
taking the numerical deicing curve as an input object, and performing regression interpolation by using a regression Kriging model to obtain an interpolation model;
and outputting interpolation results of the positions of all the characteristic points to obtain a smooth parametric deicing curve with noise points on the numerical deicing curve filtered.
7. The parameterization method for the ice shape of the icing wind tunnel icing test according to claim 6, characterized by comprising the following steps of: taking the numerical deicing curve as an input object, and performing regression interpolation by using a regression Kriging model to obtain an interpolation model, wherein the interpolation model specifically comprises the following steps:
the data point set of the numerical icing curve { (x)i,yi) 1, 2, …, N into a polar coordinate data point set { (r)i,θi)|i=1,2,…,N};
Using a set of polar data points { (r)i,θi) And (3) performing regression Kriging model on the | i ═ 1, 2, …, N to obtain an interpolation model Y (X).
8. The parameterization method for the ice shape of the icing wind tunnel icing test according to claim 7, characterized by comprising the following steps of: outputting interpolation results of positions of all the characteristic points to obtain a smooth parametric deicing curve with noise points on the numerical deicing curve filtered, wherein the smooth parametric deicing curve specifically comprises the following steps:
set of coordinate points of m feature points { lambdaj(ξj,ηj) 1, 2, …, m into a corresponding set of polar points { κ |j(ρj,εj)|j=1,2,…,m};
Each polar coordinate ρjSubstitution into interpolation model Y (X) yields a new set of polar coordinate points { κ'j(ρj,ε′j)|j=1,2,…,m};
Converting the coordinate points into a coordinate point set { (x) 'in a Cartesian coordinate system'j,y′j) And j equals to 1, 2, …, m, and a parametric deicing curve for filtering noise is obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010116493.0A CN111368367B (en) | 2020-02-25 | 2020-02-25 | Parameterization method for icing shape of icing wind tunnel icing test |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010116493.0A CN111368367B (en) | 2020-02-25 | 2020-02-25 | Parameterization method for icing shape of icing wind tunnel icing test |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111368367A true CN111368367A (en) | 2020-07-03 |
CN111368367B CN111368367B (en) | 2023-08-25 |
Family
ID=71204404
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010116493.0A Active CN111368367B (en) | 2020-02-25 | 2020-02-25 | Parameterization method for icing shape of icing wind tunnel icing test |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111368367B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112985753A (en) * | 2021-05-06 | 2021-06-18 | 中国空气动力研究与发展中心计算空气动力研究所 | Ice shape parameterization method for wind tunnel test |
CN114169256A (en) * | 2022-02-14 | 2022-03-11 | 中国空气动力研究与发展中心低速空气动力研究所 | Method for obtaining optimal wind speed of icing wind tunnel in scaling test working condition |
CN114516403A (en) * | 2022-03-07 | 2022-05-20 | 中国空气动力研究与发展中心低速空气动力研究所 | Power control method of electric heating ice prevention and removal system based on optical fiber icing detection |
CN115795696A (en) * | 2023-01-31 | 2023-03-14 | 中国空气动力研究与发展中心计算空气动力研究所 | 5363 method, device, equipment and medium for generating ice shape in process of designing Rong Bing airfoil profile |
CN118351336A (en) * | 2024-06-17 | 2024-07-16 | 上海交通大学四川研究院 | Icing wind tunnel test similarity conversion method based on ice complexity |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014006787A (en) * | 2012-06-26 | 2014-01-16 | Honda Motor Co Ltd | Feature point determination device, feature point determination method and program |
CN105869212A (en) * | 2016-03-11 | 2016-08-17 | 空气动力学国家重点实验室 | Ice shape simplification method |
JP2016218936A (en) * | 2015-05-26 | 2016-12-22 | アイシン・エィ・ダブリュ株式会社 | Curve section identification system, curve section identification method, and curve section identification program |
US20170039761A1 (en) * | 2014-05-14 | 2017-02-09 | Huawei Technologies Co., Ltd. | Image Processing Method And Apparatus |
JP2017045105A (en) * | 2015-08-24 | 2017-03-02 | 横浜ゴム株式会社 | Method and device for designing die shape of tire, and program |
CN106951657A (en) * | 2017-03-31 | 2017-07-14 | 西安石油大学 | One kind abrasion steel rail grinding target profile Fast design method |
CN107945266A (en) * | 2017-11-29 | 2018-04-20 | 中国航空工业集团公司沈阳飞机设计研究所 | A kind of generation method of wing three-dimensional icing model |
CN109766578A (en) * | 2018-12-10 | 2019-05-17 | 南京航空航天大学 | A kind of three-dimensional wing icing rear surface reconstructing method |
CN109899248A (en) * | 2019-03-04 | 2019-06-18 | 天津工业大学 | A kind of Blades For Horizontal Axis Wind rime ice ice shape forming method based on fitting of a polynomial |
-
2020
- 2020-02-25 CN CN202010116493.0A patent/CN111368367B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014006787A (en) * | 2012-06-26 | 2014-01-16 | Honda Motor Co Ltd | Feature point determination device, feature point determination method and program |
US20170039761A1 (en) * | 2014-05-14 | 2017-02-09 | Huawei Technologies Co., Ltd. | Image Processing Method And Apparatus |
JP2016218936A (en) * | 2015-05-26 | 2016-12-22 | アイシン・エィ・ダブリュ株式会社 | Curve section identification system, curve section identification method, and curve section identification program |
JP2017045105A (en) * | 2015-08-24 | 2017-03-02 | 横浜ゴム株式会社 | Method and device for designing die shape of tire, and program |
CN105869212A (en) * | 2016-03-11 | 2016-08-17 | 空气动力学国家重点实验室 | Ice shape simplification method |
CN106951657A (en) * | 2017-03-31 | 2017-07-14 | 西安石油大学 | One kind abrasion steel rail grinding target profile Fast design method |
CN107945266A (en) * | 2017-11-29 | 2018-04-20 | 中国航空工业集团公司沈阳飞机设计研究所 | A kind of generation method of wing three-dimensional icing model |
CN109766578A (en) * | 2018-12-10 | 2019-05-17 | 南京航空航天大学 | A kind of three-dimensional wing icing rear surface reconstructing method |
CN109899248A (en) * | 2019-03-04 | 2019-06-18 | 天津工业大学 | A kind of Blades For Horizontal Axis Wind rime ice ice shape forming method based on fitting of a polynomial |
Non-Patent Citations (6)
Title |
---|
刘藤 等: "基于降阶模型的翼型结冰冰形预测方法", 北京航空航天大学学报, vol. 45, no. 05, pages 1033 - 1041 * |
周志宏;易贤;郭龙;桂业伟;郭民;: "一种结冰外形相似度评估方法", no. 05, pages 14 - 19 * |
容江磊;谷正气;杨易;江涛;尹郁琦;: "基于Kriging模型的跑车尾翼断面形状的气动优化", 中国机械工程, no. 02, pages 121 - 125 * |
易贤 等: "飞机结冰冰形测量方法研究进展", vol. 38, no. 02, pages 18 - 29 * |
李伟斌 等: "基于变分分割模型的结冰冰形测量方法", 航空学报, vol. 38, no. 01, pages 100 - 107 * |
李永红,黄勇,陶洋: "基于Kriging模型的栅格翼外形参数优化研究", pages 631 - 636 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112985753A (en) * | 2021-05-06 | 2021-06-18 | 中国空气动力研究与发展中心计算空气动力研究所 | Ice shape parameterization method for wind tunnel test |
CN112985753B (en) * | 2021-05-06 | 2021-10-19 | 中国空气动力研究与发展中心计算空气动力研究所 | Ice shape parameterization method for wind tunnel test |
CN114169256A (en) * | 2022-02-14 | 2022-03-11 | 中国空气动力研究与发展中心低速空气动力研究所 | Method for obtaining optimal wind speed of icing wind tunnel in scaling test working condition |
CN114169256B (en) * | 2022-02-14 | 2022-04-19 | 中国空气动力研究与发展中心低速空气动力研究所 | Method for obtaining optimal wind speed of icing wind tunnel in scaling test working condition |
CN114516403A (en) * | 2022-03-07 | 2022-05-20 | 中国空气动力研究与发展中心低速空气动力研究所 | Power control method of electric heating ice prevention and removal system based on optical fiber icing detection |
CN114516403B (en) * | 2022-03-07 | 2023-05-05 | 中国空气动力研究与发展中心低速空气动力研究所 | Power control method of electrothermal ice prevention and removal system based on optical fiber icing detection |
CN115795696A (en) * | 2023-01-31 | 2023-03-14 | 中国空气动力研究与发展中心计算空气动力研究所 | 5363 method, device, equipment and medium for generating ice shape in process of designing Rong Bing airfoil profile |
CN118351336A (en) * | 2024-06-17 | 2024-07-16 | 上海交通大学四川研究院 | Icing wind tunnel test similarity conversion method based on ice complexity |
CN118351336B (en) * | 2024-06-17 | 2024-08-20 | 上海交通大学四川研究院 | Icing wind tunnel test similarity conversion method based on ice complexity |
Also Published As
Publication number | Publication date |
---|---|
CN111368367B (en) | 2023-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111368367A (en) | Parameterization method for ice shape of icing wind tunnel icing test | |
CN109284677B (en) | Bayesian filtering target tracking algorithm | |
CN106202804B (en) | Complex appearance aircraft distributed heat environmental parameter prediction technique based on database | |
CN113538689B (en) | Three-dimensional model grid simplifying method based on feature fusion of neural network | |
CN101639934B (en) | SAR image denoising method based on contour wave domain block hidden Markov model | |
CN101853531B (en) | Helicopter flight state identification method based on presort technology and RBF (Radial Basis Function) neural network | |
CN110930387A (en) | Fabric defect detection method based on depth separable convolutional neural network | |
CN107238500B (en) | Vehicle handling stability tests RES(rapid evaluation system) method for building up | |
CN106126791A (en) | A kind of hypersonic wing aerodynamic force/heat analysis method considering geometrical uncertainty | |
CN111754447B (en) | Infrared and visible light image fusion method based on multi-state up-down Wen Yin Markov model | |
CN113777931B (en) | Icing wing type pneumatic model construction method, device, equipment and medium | |
CN111144581A (en) | Machine learning hyper-parameter adjusting method and system | |
CN105868853A (en) | Method for predicting short-term wind power combination probability | |
CN102706573A (en) | Fault classification diagnosis method of equipment | |
CN108983800A (en) | A kind of aspect control method based on deep learning | |
CN103684349B (en) | A kind of kalman filter method based on recursion Estimates on Covariance Matrix | |
CN115311535A (en) | Icing wind tunnel cloud and mist field parameter acquisition method and storage medium | |
CN118133430A (en) | Aerodynamic electromagnetic coupling design method for aircraft based on generated diffusion model | |
CN118134984A (en) | Three-dimensional reconstruction algorithm and three-dimensional reconstruction system based on characteristic boundary line | |
CN105978733A (en) | Network flow modelling method and system based on Weibull distribution | |
CN117195402A (en) | Flying wing type unmanned aerial vehicle aircraft nose wave speed prediction model modeling method based on convolutional neural network | |
CN112085272A (en) | Wind power prediction method | |
CN116885711A (en) | Wind power prediction method, device, equipment and readable storage medium | |
CN106934124B (en) | Adaptive variable window method based on measurement change detection | |
CN115828421A (en) | Helicopter noise early warning method |
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 |