CN105550424A - Screening method of interpolating sequences based on RBF grid deformation - Google Patents
Screening method of interpolating sequences based on RBF grid deformation Download PDFInfo
- Publication number
- CN105550424A CN105550424A CN201510902301.8A CN201510902301A CN105550424A CN 105550424 A CN105550424 A CN 105550424A CN 201510902301 A CN201510902301 A CN 201510902301A CN 105550424 A CN105550424 A CN 105550424A
- Authority
- CN
- China
- Prior art keywords
- mesh
- distortion
- initial
- rbf
- sequence
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Image Processing (AREA)
- Complex Calculations (AREA)
Abstract
The invention belongs to the technical field of numerical simulation of computational grid generation of computational fluid mechanics, and specifically relates to a screening method of interpolating sequences based on RBF grid deformation. The method comprises the following steps of: 1, building a RBF grid deformation module; 2, screening initial interpolating sequences; 3, calculating an interpolation error; 4, searching an extremum value of the interpolation error delta i through an extremum algorithm; 5, screening a RBF interpolating sequence. In the method provided by the invention, efficiency of screening the interpolating sequences is improved by means of combination of the extremum algorithm and a greedy algorithm, the extremum value of the interpolation error is searched in screening iteration, an extremum value point of which the error is greater than a threshold value is screened and added into the sequence, and thus, fast simplification of the radial basis function interpolating sequences is realized.
Description
Technical field
The invention belongs to Computational fluid mechanics numerical simulation computing grid generation technique field, be specifically related to a kind of screening technique based on RBF distortion of the mesh interpolating sequence.
Background technology
Computing grid robotization generation is the importance of the numerical simulation pre-treatments such as flight vehicle aerodynamic configuration design optimization, aeroelasticity, at present, and the normal automatic generation adopting mesh deformation technique to realize optimizing configuration and gas bullet distortion body fitted grids in engineer applied.Mesh deformation technique, by carrying out translation, rotation and torsional deflection operation to computing grid, under the conforming prerequisite of network topology, enables computing grid accurately reflect the change of shape of aircraft.From several aspects comparative analyses such as grid applicability, deformability, deformation efficiency, texturing quality, algorithm complexes, carrying out distortion of the mesh based on RBF interpolation is the good grid deforming method of integration capability, in recent years obtain fast development, and be successfully applied in elastic body Pneumatic Calculation distortion of the mesh analysis.
Based on calculated amount and the N of RBF distortion of the mesh
spbecome duplicate ratio relation, memory consumption and N
spbecome cube than relation, N in base
spbe the basic point number of the Interpolation Property of Radial Basis Function sequence of carrying out interpolation, its magnitude is generally 10
4~ 10
6.For the 3D grid that grid node number is more, directly very large for RBF interpolating sequence carries out distortion of the mesh calculated amount with net point on all deformation boundaries faces.In order to improve distortion of the mesh counting yield, the memory consumption in Reducing distortion, must simplify RBF interpolating sequence.Often adopt greedy algorithm to the screening of RBF interpolating sequence at present, its algorithm idea will occur that the point of maximum error is brought in interpolating sequence in screening iteration, thus screening obtains little as far as possible N under the prerequisite meeting interpolation precision
sp.Screening process iteration of greedy algorithm only adds a node to interpolating sequence, and error convergence speed is slow, and screening process may need to screen iteration many times.
Summary of the invention
Large in order to solve the calculated amount of distortion of the mesh in background technology, the slow-footed problem of error convergence of greedy algorithm screening RBF interpolating sequence, the present invention proposes one and can accelerate interpolation error speed of convergence, improve the screening technique based on RBF distortion of the mesh interpolating sequence of interpolating sequence screening effeciency.
Concrete technical scheme of the present invention is:
Based on a screening technique for RBF distortion of the mesh interpolating sequence, it is characterized in that, comprise the following steps:
1) RBF distortion of the mesh module is built:
2) initial interpolating sequence is screened
Setting initial mesh deformation boundaries face S, arbitrary mess node is chosen in initial mesh deformation boundaries face, as the basic point P of the initial radial base interpolating sequence of RBF distortion of the mesh module
0={ p
1, p
2..., p
m; M represents the number of the basic point of initial radial base interpolating sequence
3) run RBF distortion of the mesh module, obtain the second distortion of the mesh boundary surface S
1, calculate the interpolation error S on the second distortion of the mesh boundary surface and initial mesh deformation boundaries
i;
4) extreme value algorithm search interpolation error δ is adopted
iextreme value
4.1) initial mesh deformation boundaries face is divided into N number of region of search Ω
i
4.2) each region of search Ω is searched for
imaximum error
Wherein k represents search Ω
ithe number of interior interpolation error;
Complete N number of region of search Ω
iextremum search, define a series of extreme value
5) RBF interpolating sequence is screened
5.1) step-up error threshold values σ, σ are 10 of initial mesh deformation boundaries face maximum deformation quantity
-3
5.2) selecting step 4.2) extreme value that formed
in any one
if
Then remember appearance
region of search on node be p
m+i;
If
Then filtering occurs
region of search;
5.3) if p
m+iwith initial radial base interpolating sequence P
0basic point repeat, then filtering; Otherwise carry out step 5.4)
5.4) if p
m+ithe basic point P of initial radial base interpolating sequence
0adjacent node, then filtering; Otherwise retain p
m+i
5.5) step 5.2 is repeated) to 5.4) filter out extreme value point range p
m+1..., p
m+s, according to the principle of greedy algorithm, form new radial base interpolation sequence
P
1={ p
1, p
2..., p
m, p
m+1..., p
m+s, behalf step 5.2) to 5.4) the middle node number inserting initial radial base interpolating sequence;
6) according to step 5.5) the radial base interpolation sequence that formed, repeat step 3) to step 5), until all extreme values
all be not more than given error threshold values σ, then screening process terminates.
Above-mentioned RBF distortion of the mesh module is used for the input of distortion of the mesh controling parameters, the input of initial calculation grid, creates deformation boundaries face, creates distortion of the mesh region, creates Initial R BF interpolating sequence.
The invention has the advantages that:
Adopt method of the present invention in the iteration of screening radial basis function interpolating sequence, error convergence speed is fast, and screening effeciency is high, realizes simplifying fast of Interpolation Property of Radial Basis Function sequence.This method screens the radial base interpolation sequence obtained, and can ensure the interpolation precision of distortion of the mesh, has significantly reduced the number of the calculated amount of distortion of the mesh simultaneously, for shortening Flight Vehicle Design optimization and aeroelasticity numerical simulation cycle provide technical support.
Accompanying drawing explanation
Fig. 1 is Interpolation Property of Radial Basis Function sequence screening process flow diagram of the present invention
Embodiment
Under the requirement meeting interpolation precision, simplify RBF interpolating sequence be reduce RBF distortion of the mesh calculated amount most important technological means.On initial interpolating sequence basis, by deformation boundaries face being screened the basic point of the large node of RBF distortion of the mesh interpolation error as sequence, reduce interpolation error amount gradually, until reach the error precision of setting.Method of the present invention is intended to accelerate error convergence speed, improves RBF distortion of the mesh iteration screening effeciency, establishes the interpolating sequence screening technique of extreme value algorithm in conjunction with greedy algorithm.
Specific descriptions below in conjunction with accompanying drawing 1 pair of the inventive method are as follows:
1) RBF distortion of the mesh module is built;
Wherein, RBF distortion of the mesh module grid is used for Deformation control parameters input, the input of initial calculation grid, creates deformation boundaries face, creates distortion of the mesh region, creates Initial R BF interpolating sequence etc., provides necessary pretreatment information for distortion of the mesh calculates.
Specifically:
Distortion of the mesh controling parameters main definitions is:
◆ the type of definition distortion of the mesh, comprises Grid Mobile, rotation, torsional deflection and aggregate motion distortion thereof.
◆ definition grid rotates, moving movement form.
◆ deflection, the form of distortion of the data point of definition distortion of the mesh.
◆ definition distortion of the mesh restricted area.
Initial calculation grid is input as:
◆ the net point coordinate of input initial calculation grid.
◆ input initial mesh network topology and boundary information.
Creating deformation boundaries face is:
◆ the boundary surface of initial calculation grid is created according to distortion of the mesh controling parameters.
◆ the deformation boundaries face of distortion of the mesh is created according to distortion of the mesh controling parameters.
◆ for grid torsional deflection, carry out geometric parameter modeling.
Creating distortion of the mesh region is:
◆ in initial mesh, create distortion of the mesh region according to distortion of the mesh controling parameters, simplify grid scale to be deformed in computational fields.
Creating Initial R BF interpolating sequence is:
◆ according to network topology and net point distribution, choose node on deformation boundaries as basic point.Create Initial R BF interpolating sequence.
◆ create RBF interpolating sequence file, initial interpolating sequence is outputted to file.
2) initial interpolating sequence is screened; For the interpolating sequence chosen, the net point displacement of sequence basic point is accurate, and the net point displacement of non-sequence basic point is obtained by RBF method interpolation, there is interpolation error.By increasing the basic point number of interpolating sequence, reduce error amount gradually, specific practice is:
Setting initial mesh deformation boundaries face S, arbitrary mess node is chosen in initial mesh deformation boundaries face, as the basic point P of the initial radial base interpolating sequence of RBF distortion of the mesh module
0={ p
1, p
2..., p
m; M represents the number of the basic point of initial radial base interpolating sequence
3) run RBF distortion of the mesh module, obtain the second distortion of the mesh boundary surface S
1, calculate the interpolation error δ on the second distortion of the mesh boundary surface and initial mesh deformation boundaries
i;
4) extreme value algorithm search interpolation error δ is adopted
iextreme value
4.1) initial mesh deformation boundaries face is divided into N number of region of search Ω
i
4.2) each region of search Ω is searched for
imaximum error
Wherein k represents search Ω
ithe number of interior interpolation error;
Complete N number of region of search Ω
iextremum search, define a series of extreme value
5) RBF interpolating sequence is screened: this step determines whether whether added in RBF interpolating sequence by the extreme point that a program module upper in flow process is selected.Screening mode comprises threshold values screening, repeats a screening and isolatism screening.Then the node of the appearance extreme value satisfied condition to be added in existing RBF interpolating sequence and to carry out error convergence judgement.
5.1) step-up error threshold values σ, σ are 10 of initial mesh deformation boundaries face maximum deformation quantity
-3
5.2) selecting step 4.2) extreme value that formed
in any one
if
Then remember appearance
region of search on node be p
m+i;
If
Then filtering occurs
region of search;
5.3) if p
m+iwith initial radial base interpolating sequence P
0basic point repeat, then filtering; Otherwise carry out step 5.4)
5.4) if p
m+ithe basic point P of initial radial base interpolating sequence
0adjacent node, then filtering; Otherwise retain p
m+i
5.5) step 5.2 is repeated) to 5.4) filter out extreme value point range p
m+1..., p
m+s, according to the principle of greedy algorithm, form new radial base interpolation sequence
P
1={ p
1, p
2..., p
m, p
m+1..., p
m+s, behalf step 5.2) to 5.4) the middle node number inserting initial radial base interpolating sequence
6) according to step 5.5) the radial base interpolation sequence that formed, repeat step 3) to step 5), until all extreme values
all be not more than given error threshold values σ, then screening process terminates.
Claims (2)
1. based on a screening technique for RBF distortion of the mesh interpolating sequence, it is characterized in that, comprise the following steps:
1) RBF distortion of the mesh module is built:
2) initial interpolating sequence is screened
Setting initial mesh deformation boundaries face S, arbitrary mess node is chosen in initial mesh deformation boundaries face, as the basic point P of the initial radial base interpolating sequence of RBF distortion of the mesh module
0={ p
1, p
2..., p
m; M represents the number of the basic point of initial radial base interpolating sequence
3) run RBF distortion of the mesh module, obtain the second distortion of the mesh boundary surface S
1, calculate the interpolation error δ on the second distortion of the mesh boundary surface and initial mesh deformation boundaries
i;
4) extreme value algorithm search interpolation error δ is adopted
iextreme value
4.1) initial mesh deformation boundaries face is divided into N number of region of search Ω
i
4.2) each region of search Ω is searched for
imaximum error
Wherein k represents search Ω
ithe number of interior interpolation error; Complete N number of region of search Ω
iextremum search, define a series of extreme value
5) RBF interpolating sequence is screened
5.1) step-up error threshold values σ, σ are 10 of initial mesh deformation boundaries face maximum deformation quantity
-3
5.2) selecting step 4.2) extreme value that formed
in any one
If
Then remember appearance
region of search on node be p
m+i;
If
Then filtering occurs
region of search;
5.3) if p
m+iwith initial radial base interpolating sequence P
0basic point repeat, then filtering; Otherwise carry out step 5.4)
5.4) if p
m+ithe basic point P of initial radial base interpolating sequence
0adjacent node, then filtering; Otherwise retain p
m+i;
5.5) step 5.2 is repeated) to 5.4) filter out extreme value point range p
m+1..., p
m+s, according to the principle of greedy algorithm, form new radial base interpolation sequence
P
1={ p
1, p
2..., p
m, p
m+1..., p
m+s, behalf step 5.2) to 5.4) the middle node number inserting initial radial base interpolating sequence;
6) according to step 5.5) the radial base interpolation sequence that formed, repeat step 3) to step 5), until all extreme values
all be not more than given error threshold values σ, then screening process terminates.
2. the screening technique based on RBF distortion of the mesh interpolating sequence according to claim 1, it is characterized in that: described RBF distortion of the mesh module is used for the input of distortion of the mesh controling parameters, the input of initial calculation grid, creates deformation boundaries face, creates distortion of the mesh region, creates Initial R BF interpolating sequence.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510902301.8A CN105550424B (en) | 2015-12-08 | 2015-12-08 | A kind of screening technique based on RBF distortion of the mesh interpolating sequences |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510902301.8A CN105550424B (en) | 2015-12-08 | 2015-12-08 | A kind of screening technique based on RBF distortion of the mesh interpolating sequences |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105550424A true CN105550424A (en) | 2016-05-04 |
CN105550424B CN105550424B (en) | 2018-10-26 |
Family
ID=55829613
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510902301.8A Active CN105550424B (en) | 2015-12-08 | 2015-12-08 | A kind of screening technique based on RBF distortion of the mesh interpolating sequences |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105550424B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107403466A (en) * | 2017-08-03 | 2017-11-28 | 中国空气动力研究与发展中心计算空气动力研究所 | Ultra-large unstrctured grid generation method based on overall situation encryption |
CN116227043A (en) * | 2023-05-10 | 2023-06-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Aircraft numerical simulation method, system, equipment and computer storage medium |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049932A (en) * | 2012-12-20 | 2013-04-17 | 北京农业信息技术研究中心 | Radial basis function-based plant three-dimensional configuration virtual modeling method |
CN103530472A (en) * | 2013-10-24 | 2014-01-22 | 山东师范大学 | Three-dimensional model automation simplification method based on importance sampling |
-
2015
- 2015-12-08 CN CN201510902301.8A patent/CN105550424B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049932A (en) * | 2012-12-20 | 2013-04-17 | 北京农业信息技术研究中心 | Radial basis function-based plant three-dimensional configuration virtual modeling method |
CN103530472A (en) * | 2013-10-24 | 2014-01-22 | 山东师范大学 | Three-dimensional model automation simplification method based on importance sampling |
Non-Patent Citations (3)
Title |
---|
T.C.S. RENDALL,ET AL: "Reduced surface point selection options for efficient mesh deformation using radial basis function", 《JOURNAL OF COMPUTATIONAL PHYSICS》 * |
王刚 等: "一种基于径向基函数的非结构混合网格变形技术", 《西北工业大学学报》 * |
谢亮: "基于径向基函数的高效网格变形算法研究", 《振动与冲击》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107403466A (en) * | 2017-08-03 | 2017-11-28 | 中国空气动力研究与发展中心计算空气动力研究所 | Ultra-large unstrctured grid generation method based on overall situation encryption |
CN116227043A (en) * | 2023-05-10 | 2023-06-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Aircraft numerical simulation method, system, equipment and computer storage medium |
CN116227043B (en) * | 2023-05-10 | 2024-03-12 | 中国空气动力研究与发展中心计算空气动力研究所 | Aircraft numerical simulation method, system, equipment and computer storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN105550424B (en) | 2018-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7191161B1 (en) | Method for constructing composite response surfaces by combining neural networks with polynominal interpolation or estimation techniques | |
CN107273569B (en) | Aerodynamic configuration drag reduction optimization method based on mesh deformation technique | |
CN105631125A (en) | Aerodynamic-thermal-structural coupling analysis method based on reduced-order model | |
CN111597631B (en) | Automobile wind resistance coefficient optimization method based on self-adaptive agent model | |
CN114996658B (en) | Projection-based hypersonic aircraft aerodynamic heat prediction method | |
CN109726465A (en) | The three-dimensional method for numerical simulation streamed without viscous low speed based on non-structural curl grid | |
CN111898212A (en) | Impeller mechanical profile design optimization method based on BezierGAN and Bayesian optimization | |
CN104750948A (en) | Optimization method for processing multiple extreme values and multiple restricted problems in flight vehicle design | |
CN115688276A (en) | Aircraft appearance automatic optimization method, system, equipment and medium based on discrete companion method | |
Tang et al. | Generation of aerodynamic data using a design of experiment and data fusion approach | |
CN105550424A (en) | Screening method of interpolating sequences based on RBF grid deformation | |
CN115544815A (en) | Method and device for generating fan model | |
CN114936428A (en) | ADAMS-based steering mechanism multi-body dynamics modeling and parameter confirmation method | |
US9087165B2 (en) | Automatic extremum detection on a surface mesh of a component | |
CN108197353A (en) | A kind of solid propellant rocket Fixture Design method of the APDL language based on ANSYS | |
CN105243183A (en) | Product structure parameterized design optimization method and system | |
CN114169100B (en) | Efficient design optimization method and system for super-large variable impeller machinery and application | |
CN103218493A (en) | Fast isogeometric analysis numerical simulation method based on multiple grids | |
US20200082035A1 (en) | Simulation method, physical quantity calculation program, and physical quantity calculation apparatus | |
CN105512352A (en) | Analysis method aiming at pneumatic loads | |
CN109241585B (en) | High-low pressure turbine transition flow passage profile inverse problem design method | |
CN110555267A (en) | Parameterized horizontal set structure topology optimization method based on implicit B-spline | |
CN110633509A (en) | Cartesian robot cantilever optimization simulation method based on MOGA algorithm | |
CN106407600B (en) | It is a kind of that water simulating analysis is entered based on fluid-wall interaction and agent model | |
CN115048730B (en) | Axial symmetry supersonic velocity spray pipe optimization design method and device based on grid displacement |
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 |