CN105550424B - A kind of screening technique based on RBF distortion of the mesh interpolating sequences - Google Patents
A kind of screening technique based on RBF distortion of the mesh interpolating sequences Download PDFInfo
- Publication number
- CN105550424B CN105550424B CN201510902301.8A CN201510902301A CN105550424B CN 105550424 B CN105550424 B CN 105550424B CN 201510902301 A CN201510902301 A CN 201510902301A CN 105550424 B CN105550424 B CN 105550424B
- Authority
- CN
- China
- Prior art keywords
- mesh
- distortion
- rbf
- initial
- interpolating
- 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.)
- Active
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
Abstract
The invention belongs to Computational fluid mechanics numerical simulations to calculate grid generation technique field, and in particular to a kind of screening technique based on RBF distortion of the mesh interpolating sequences.This approach includes the following steps：1, RBF distortion of the mesh modules are built；2, initial interpolating sequence is screened；3, interpolation error is calculated；4, using extreme value algorithm search interpolation error δiExtreme value；5, RBF interpolating sequences are screened；This method improves interpolating sequence screening efficiency using the method that extreme value algorithm and greedy algorithm combine, in screening iteration, interpolation error extreme value is found, screening error is more than the extreme point of threshold values and is added in sequence, realizes quickly simplifying for Interpolation Property of Radial Basis Function sequence.
Description
Technical field
The invention belongs to Computational fluid mechanics numerical simulations to calculate grid generation technique field, and in particular to one kind is based on
The screening technique of RBF distortion of the mesh interpolating sequences.
Background technology
It is the numerical simulations pre-treatments such as flight vehicle aerodynamic configuration design optimization, aeroelasticity to calculate grid automation and generate
Importance, currently, realizing that optimization configuration and gas bullet deform body fitted grids frequently with mesh deformation technique in engineer application
It automatically generates.Mesh deformation technique by calculate grid translated, rotate and torsional deflection operate, it is consistent in network topology
Under the premise of property, calculating grid is enable to accurately reflect the change in shape of aircraft.From grid applicability, deformability, deformation
Several aspect comparative analyses such as efficiency, texturing quality, algorithm complexity, it is integration capability to carry out distortion of the mesh based on RBF interpolation
Preferable grid deforming method obtains fast development in recent years, and has been successfully applied to elastomer Pneumatic Calculation distortion of the mesh analysis
In.
Calculation amount based on RBF distortions of the mesh and NspAt duplicate ratio relationship, memory consumption and NspAt cube than relationship, base
Middle NspIt is into the basic point number of the Interpolation Property of Radial Basis Function sequence of row interpolation, magnitude is generally 104~106.For grid node
It is directly that RBF interpolating sequences carry out in terms of distortion of the mesh by mesh point on all deformation boundaries faces for the more three-dimensional grid of number
Calculation amount is very big.In order to improve distortion of the mesh computational efficiency, the memory consumption in Reducing distortion, it is necessary to simplify RBF interpolating sequences.Mesh
Preceding to screen RBF interpolating sequences frequently with greedy algorithm, algorithm idea is that will occur the point of worst error in screening iteration
It brings into interpolating sequence, to which screening obtains N as small as possible under the premise of meeting interpolation precisionsp.The screening of greedy algorithm
Process an iteration only adds a node to interpolating sequence, and error convergence speed is slow, and screening process may need to sieve many times
Select iteration.
Invention content
In order to solve computationally intensive, the error receipts of greedy algorithm screening RBF interpolating sequences of distortion of the mesh in background technology
Slow-footed problem is held back, the present invention, which proposes one kind, can accelerate interpolation error convergence rate, improve interpolating sequence screening efficiency
The screening technique based on RBF distortion of the mesh interpolating sequences.
The specific technical solution of the present invention is：
A kind of screening technique based on RBF distortion of the mesh interpolating sequences, which is characterized in that include the following steps：
1) RBF distortion of the mesh modules are built；
2) initial interpolating sequence is screened；
Initial mesh deformation boundaries face S is set, arbitrary mess node is chosen on initial mesh deformation boundaries face, as
The basic point P of the initial radial base interpolating sequence of RBF distortion of the mesh modules0={ p1, p2..., pm}；M represents initial radial base interpolation
The number of the basic point of sequence；
3) RBF distortion of the mesh modules are run, the second distortion of the mesh boundary face S is obtained1, calculate the second distortion of the mesh boundary face
With the interpolation error δ on initial mesh deformation boundariesi；
4) extreme value algorithm search interpolation error δ is usediExtreme value；
4.1) initial mesh deformation boundaries face is divided into N number of region of search Ωj, j=1,2,3 ..., N；
4.2) each region of search Ω is searched forjWorst error；
Wherein k indicates search ΩjThe number of interior interpolation error；Complete N number of region of search ΩjExtremum search, form one
The extreme value of series
5) RBF interpolating sequences are screened；
5.1) step-up error threshold values σ, σ are the 10 of initial mesh deformation boundaries face maximum deformation quantity-3；
5.2) selecting step 4.2) formed extreme valueIn any oneIf
Then remember appearanceRegion of search on node be pm+j；
If
Then filter out appearanceRegion of search；
If 5.3) pm+jWith initial radial base interpolating sequence P0Basic point repeat, then filter out；Otherwise step 5.4) is carried out；
If 5.4) pm+jIt is the basic point P of initial radial base interpolating sequence0Adjacent node then filters out；Otherwise retain pm+j；
5.5) step 5.2) is repeated to 5.4) filtering out extreme value point range pm+1,…,pm+s, according to the principle of greedy algorithm, shape
The radial base interpolation sequence P of Cheng Xin1={ p1,p2,…,pm, pm+1,…,pm+s, behalf step 5.2) is inserted into initially to 5.4) middle
The node number of radial base interpolation sequence；
6) the radial base interpolation sequence formed according to step 5.5) repeats step 3) to step 5), until all extreme valuesNo more than given error threshold values σ, then screening process terminates.
Above-mentioned RBF distortions of the mesh module is for the input of distortion of the mesh control parameter, the input of initial calculation grid, establishment deformation
Boundary face creates distortion of the mesh region, creates Initial R BF interpolating sequences.
The advantage of the invention is that：
For method using the present invention in the iteration of screening radial basis function interpolating sequence, error convergence speed is fast, screening
It is efficient, realize quickly simplifying for Interpolation Property of Radial Basis Function sequence.The radial base interpolation sequence that this method is screened, Ke Yibao
The interpolation precision of distortion of the mesh is demonstrate,proved, while substantially having reduced the number of the calculation amount of distortion of the mesh, to shorten Flight Vehicle Design optimization and gas
The dynamic elasticity numerical simulation period provides technical support.
Description of the drawings
Fig. 1 is the Interpolation Property of Radial Basis Function sequence screening flow chart of the present invention
Specific implementation mode
Under the requirement for meeting interpolation precision, it is the most heavy of reduction RBF distortion of the mesh calculation amounts to simplify RBF interpolating sequences
Want technological means.On the basis of initial interpolating sequence, by screening the big section of RBF distortion of the mesh interpolation errors on deformation boundaries face
Basic point of the point as sequence, continuously decreases interpolation error amount, the error precision until reaching setting.The method of the present invention is intended to add
Fast error convergence speed, improves RBF distortion of the mesh iteration screening efficiencies, establishes the interpolation sequence of extreme value algorithm combination greedy algorithm
Column filter method.
1 the method for the present invention is described in detail below below in conjunction with the accompanying drawings：
1) RBF distortion of the mesh modules are built；
Wherein, RBF distortions of the mesh module grid becomes for the input of Deformation control parameter, the input of initial calculation grid, establishment
Shape boundary face, create distortion of the mesh region, create Initial R BF interpolating sequences etc., provide necessary pre- place for distortion of the mesh calculating
Manage information.
Specifically：
Distortion of the mesh control parameter main definitions are：
◆ define the type of distortion of the mesh, including Grid Mobile, rotation, torsional deflection and combinations thereof motion deformation.
◆ define grid rotation, moving movement form.
◆ define deflection, the variant of the data point of distortion of the mesh.
◆ define distortion of the mesh restricted area.
Initial calculation grid inputs：
◆ the grid point coordinates of input initial calculation grid.
◆ input initial mesh network topology and boundary information.
Creating deformation boundaries face is：
◆ the boundary face of initial calculation grid is created according to distortion of the mesh control parameter.
◆ the deformation boundaries face of distortion of the mesh is created according to distortion of the mesh control parameter.
◆ for grid torsional deflection, carry out geometric parameter modeling.
Creating distortion of the mesh region is：
◆ distortion of the mesh region is created in initial mesh according to distortion of the mesh control parameter, is simplified to be deformed in computational domain
Grid scale.
Creating Initial R BF interpolating sequences is：
◆ it is distributed according to network topology and mesh point, the node on selection deformation boundaries is as basic point.Initial R BF is created to insert
Value sequence.
◆ RBF interpolating sequence files are created, initial interpolating sequence is output to file.
2) initial interpolating sequence is screened；For the interpolating sequence of selection, the mesh point displacement of sequence basic point be it is accurate, it is non-
The mesh point displacement of sequence basic point is obtained by RBF method interpolation, and there is interpolation errors.By the basic point for increasing interpolating sequence
Number, is gradually reduced error amount, specific practice is：
Initial mesh deformation boundaries face S is set, arbitrary mess node is chosen on initial mesh deformation boundaries face, as
The basic point P of the initial radial base interpolating sequence of RBF distortion of the mesh modules0={ p1, p2..., pm}；M represents initial radial base interpolation
The number of the basic point of sequence
3) RBF distortion of the mesh modules are run, the second distortion of the mesh boundary face S is obtained1, calculate the second distortion of the mesh boundary face
With the interpolation error δ on initial mesh deformation boundariesi；
4) extreme value algorithm search interpolation error δ is usediExtreme value；
4.1) initial mesh deformation boundaries face is divided into N number of region of search Ωj；
4.2) each region of search Ω is searched forjWorst error；
Wherein k indicates search ΩjThe number of interior interpolation error；
Complete N number of region of search ΩjExtremum search, form a series of extreme value
5) RBF interpolating sequences are screened：The step is determined whether a upper program module is selected in flow extreme point
Whether it is added in RBF interpolating sequences.Screening mode includes threshold values screening, repeats point screening and isolatism screening.Then will expire
The node of the appearance extreme value of sufficient condition is added in existing RBF interpolating sequences and carries out error convergence judgement.
5.1) step-up error threshold values σ, σ are the 10 of initial mesh deformation boundaries face maximum deformation quantity-3；
5.2) selecting step 4.2) formed extreme valueIn any oneIf
Then remember appearanceRegion of search on node be pm+j；
If
Then filter out appearanceRegion of search；
If 5.3) pm+jWith initial radial base interpolating sequence P0Basic point repeat, then filter out；Otherwise step 5.4) is carried out；
If 5.4) pm+jIt is the basic point P of initial radial base interpolating sequence0Adjacent node then filters out；Otherwise retain pm+j；
5.5) step 5.2) is repeated to 5.4) filtering out extreme value point range pm+1,…,pm+s, according to the principle of greedy algorithm, shape
The radial base interpolation sequence P of Cheng Xin1={ p1,p2,…,pm,pm+1,…,pm+s, s indicates that step 5.2) is inserted into initially to 5.4) middle
The node number of radial base interpolation sequence；
6) the radial base interpolation sequence formed according to step 5.5) repeats step 3) to step 5), until all extreme valuesNo more than given error threshold values σ, then screening process terminates.
Claims (2)
1. a kind of screening technique based on RBF distortion of the mesh interpolating sequences, which is characterized in that include the following steps：
1) RBF distortion of the mesh modules are built；
2) initial interpolating sequence is screened；
Initial mesh deformation boundaries face S is set, arbitrary mess node is chosen on initial mesh deformation boundaries face, as RBF nets
The basic point P of the initial radial base interpolating sequence of lattice deformation module0={ p1, p2..., pm}；M represents initial radial base interpolating sequence
Basic point number；
3) RBF distortion of the mesh modules are run, the second distortion of the mesh boundary face S is obtained1, calculate the second distortion of the mesh boundary face with just
The borderline interpolation error δ of beginning distortion of the meshi；
4) extreme value algorithm search interpolation error δ is usediExtreme value；
4.1) initial mesh deformation boundaries face is divided into N number of region of search Ωl, j=1,2,3 ..., N；
4.2) each region of search Ω is searched forjWorst error；
Wherein k indicates search ΩjThe number of interior interpolation error；Complete N number of region of search ΩjExtremum search, form a series of
Extreme value
5) RBF interpolating sequences are screened；
5.1) step-up error threshold values σ, σ are the 10 of initial mesh deformation boundaries face maximum deformation quantity-3；
5.2) selecting step 4.2) formed extreme valueIn any oneIf
Then remember appearanceRegion of search on node be pm+j；
If
Then filter out appearanceRegion of search；
If 5.3) pm+jWith initial radial base interpolating sequence P0Basic point repeat, then filter out；Otherwise step 5.4) is carried out；
If 5.4) pm+jIt is the basic point P of initial radial base interpolating sequence0Adjacent node then filters out；Otherwise retain pm+j；
5.5) step 5.2) is repeated to 5.4) filtering out extreme value point range pm+1,…,pm+s, according to the principle of greedy algorithm, formed new
Radial base interpolation sequence
P1={ p1,p2,…,pm,pm+1,…,pm+s, behalf step 5.2) to 5.4) middle insertion initial radial base interpolating sequence
Node number；
6) the radial base interpolation sequence formed according to step 5.5) repeats step 3) to step 5), until all extreme valuesNo more than given error threshold values σ, then screening process terminates.
2. the screening technique according to claim 1 based on RBF distortion of the mesh interpolating sequences, it is characterised in that：The RBF
Distortion of the mesh module creates grid for the input of distortion of the mesh control parameter, the input of initial calculation grid, establishment deformation boundaries face
Deformed region creates Initial R BF interpolating sequences.
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 CN105550424A (en) | 2016-05-04 |
| CN105550424B true 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) |
Families Citing this family (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 |
| CN116227043B (en) * | 2023-05-10 | 2024-03-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 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 |
|---|
| Reduced surface point selection options for efficient mesh deformation using radial basis function;T.C.S. Rendall，et al;《Journal of Computational Physics》;20100104;全文 * |
| 一种基于径向基函数的非结构混合网格变形技术;王刚 等;《西北工业大学学报》;20111031;第29卷(第5期);全文 * |
| 基于径向基函数的高效网格变形算法研究;谢亮;《振动与冲击》;20131231;第32卷(第10期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105550424A (en) | 2016-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110705029A (en) | Flow field prediction method of oscillating flapping wing energy acquisition system based on transfer learning | |
| CN108153984B (en) | High-precision Galois gold artificial viscous shock wave capturing method based on flow field density step | |
| CN107273569B (en) | Aerodynamic configuration drag reduction optimization method based on mesh deformation technique | |
| Tang et al. | Generation of aerodynamic data using a design of experiment and data fusion approach | |
| CN107204040B (en) | Multipoint geostatistical modeling method and device and computer storage medium | |
| CN105550424B (en) | A kind of screening technique based on RBF distortion of the mesh interpolating sequences | |
| CN109388840B (en) | Special programming computer for numerical simulation of metal forming processes with predefined load paths and corresponding grid adjustment schemes | |
| CN104156546A (en) | Profile redesigning method of automobile covering part die based on T spline | |
| EP3486814A1 (en) | Dimensional reduction of finite element analysis solution for rapid emulation | |
| CN114996658B (en) | Projection-based hypersonic aircraft aerodynamic heat prediction method | |
| CN104992023A (en) | Aircraft parametric design method based on state type function | |
| CN111079326B (en) | Two-dimensional anisotropic grid cell measurement tensor field smoothing method | |
| CN104504227A (en) | Hydraulic forging press multi-dimensional virtual prototype collaborative simulation method and system | |
| CN114169100B (en) | Efficient design optimization method and system for super-large variable impeller machinery and application | |
| US9087165B2 (en) | Automatic extremum detection on a surface mesh of a component | |
| CN105243183A (en) | Product structure parameterized design optimization method and system | |
| CN111985137A (en) | Multivariable cutting level set optimization method for porous structure topology optimization | |
| CN103218493A (en) | Fast isogeometric analysis numerical simulation method based on multiple grids | |
| CN114341859A (en) | Object design process using coarse geometry cells and high resolution grid cells | |
| CN107391871A (en) | A kind of space lattice deformation method based on parallelization RBF | |
| CN114444012A (en) | Partial differential equation numerical solution method based on neural network architecture search | |
| CN110717271B (en) | Substance evolution simulation method based on exponential time difference format solution | |
| CN110555267A (en) | Parameterized horizontal set structure topology optimization method based on implicit B-spline | |
| Jiang et al. | Reduced multiscale computation on adapted grid for the convection-diffusion Robin problem | |
| JP2008276468A (en) | Lattice generation device for fluid analysis and lattice generation program for fluid analysis |
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 |