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 PDF

Info

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
Application number
CN202010116493.0A
Other languages
Chinese (zh)
Other versions
CN111368367B (en
Inventor
李伟斌
宋超
罗凤
赵凡
马洪林
王跃军
王应宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AERODYNAMICS NATIONAL KEY LABORATORY
Original Assignee
AERODYNAMICS NATIONAL KEY LABORATORY
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by AERODYNAMICS NATIONAL KEY LABORATORY filed Critical AERODYNAMICS NATIONAL KEY LABORATORY
Priority to CN202010116493.0A priority Critical patent/CN111368367B/en
Publication of CN111368367A publication Critical patent/CN111368367A/en
Application granted granted Critical
Publication of CN111368367B publication Critical patent/CN111368367B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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

Parameterization method for ice shape of icing wind tunnel icing test
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 pointsjj,η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 { lambdajj,ηj) 1, 2, …, m into a corresponding set of polar points kjj,εj)|j=1,2,…,m};
Each polar coordinate ρjSubstitution into interpolation model Y (X) yields a new set of polar coordinate points { κ'jj,ε′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 pointsjj,η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).
Specifically, the model function y (x) is:
Figure BDA0002391656640000071
wherein f isj(X) is a basis function, here chosen as constant 1, βjIs the coefficient corresponding to the basis function,
Figure BDA0002391656640000072
a mathematical expectation value representing y (x); z (X) is mean zero and variance
Figure BDA0002391656640000073
The static random process of (a); for any two positions X and X' in the design space, the covariance of the random quantity is:
Figure BDA0002391656640000074
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:
Figure BDA0002391656640000075
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 { lambdajj,ηj) 1, 2, …, m into a corresponding set of polar points { κ |jj,ε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 modeljj,ε'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 pointsjj,η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 { lambdajj,ηj) 1, 2, …, m into a corresponding set of polar points { κ |jj,εj)|j=1,2,…,m};
Each polar coordinate ρjSubstitution into interpolation model Y (X) yields a new set of polar coordinate points { κ'jj,ε′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.
CN202010116493.0A 2020-02-25 2020-02-25 Parameterization method for icing shape of icing wind tunnel icing test Active CN111368367B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (9)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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