CN106055788A - Combined modeling method of terrain and wind field - Google Patents
Combined modeling method of terrain and wind field Download PDFInfo
- Publication number
- CN106055788A CN106055788A CN201610374228.6A CN201610374228A CN106055788A CN 106055788 A CN106055788 A CN 106055788A CN 201610374228 A CN201610374228 A CN 201610374228A CN 106055788 A CN106055788 A CN 106055788A
- Authority
- CN
- China
- Prior art keywords
- landform
- wind field
- terrain
- point
- simulation
- 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.)
- Pending
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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Processing Or Creating Images (AREA)
Abstract
The invention relates to a combined modeling method of a terrain and a wind field. The combined modeling method comprises the following steps: 1) carrying out data acquisition on a terrain contour map; 2) according to acquired data, carrying out terrain modeling; and 3) inputting a wind field condition into a terrain model to carry out the combined calculation of the terrain and the wind field. The combined modeling method of the terrain and the wind field can cause the terrain to truly reappear, and finishes wind field calculation on the basis of terrain simulation. Only terrain contour data needs to be input and average wind speed in an area is simulated to realize the synchronous simulation of the terrain and the wind field; complex steps that the terrain is simulated, then, the terrain model is simplified and finally the wind field is calculated are not required; and the combined modeling method is suitable for the simulation of simple/ complex terrains and wind fields, and has the advantages of being simple in construction, small in calculated amounts, high in precision and the like.
Description
Technical field
The invention belongs to technical field of aerospace, particularly relate to a kind of landform and wind field combines modeling method.
Background technology
Existing landform is with simulation of wind method, and terrain simulation resolves with wind field and uses separate analysis model and meter
Calculation process.
The problem that terrain modeling mainly solves is according to given terrain data, sets up the curved surface of expression topographical surface true to nature
Model.At present, during Quantitative Terrain analysis is studied with three-dimensional modeling, the simplest and conventional method is data fitting method.
Wind field modeling is can restore wind field really as goal in research, according to the methods of the different simulations of research direction also
Have any different.At present wind field modeling and simulating method mainly has 3 classes: based on the theoretical emulation mode of atmospheric dynamics, actual measurement wind field number
According to method, and Engineering Simulation Method based on potential flow theory.The emulation mode theoretical based on atmospheric dynamics can be examined in detail
Considering the temperature of air, energy, the isoparametric change of air pressure, and provide the numerical solution that three-dimensional Flow over mountain is complete, the method is multiplex
In meteorological research, have the advantages that amount of calculation is huge, physical parameter is complicated, be difficult to directly apply to aeronautical engineering field.
Actual measurement wind field data method closing to reality situation the most, but need huge data space, measure work and be also required to expend greatly
Amount manpower and materials, and obtain wind field state be extremely limited, reproduction wind field difficulty the biggest.Engineering simulation method is with letter
The impact of single fluid mechanic model simulation mountain relief, the main feature of description Flow over mountain that can be macroscopical, but in landform
During modeling, parameter chooses complexity, workload relatively greatly, during analyzing, for making landform as close possible to reality, needs
Landform characteristic parameter is adjusted repeatedly, landform simplify the direct analog result affecting wind field.
Summary of the invention
It is an object of the invention to provide a kind of landform and wind field combines modeling method, solves current modeling method or builds
The problem that mould is inaccurate or workload is big.
For reaching above-mentioned purpose, the technical solution used in the present invention is: a kind of landform and the modeling method of combining of wind field, wraps
Include
Step one: landform contour map is carried out data acquisition；
Step 2: carry out terrain modeling according to the data collected；
Step 3: by wind field condition entry to relief model, carry out the Combined Calculation of landform and wind field.
Further, in step one, contour map to landform carries out data collection steps and includes
1) coordinate system is set up, with the contour map region of described landform for artificially generated terrain region；
2) by predetermined way, contour map being carried out data acquisition, the data of collection include plane coordinates and altitude data,
Described altitude data is height coordinate, and described data are referred to as known point.
Further, in described coordinate system with sea level as X/Y plane, with height above sea level Z axis coordinate.
Further, described predetermined way includes constant duration data acquisition and equidistantly height collection.
Further, in step 2, terrain modeling step includes
1) in artificially generated terrain region, the point being parallel to X/Y plane being carried out random acquisition, the point collected is as simulation
Point；
2) according to the plane coordinates of simulation point, triangular mesh is used to generate method, the network of triangle of constructing analog landform
Lattice；
3) according to plane coordinates and the altitude data of the known point collected under coordinate system, and the panel data of simulation point,
Utilize interpolation method, obtain the modality relation between known point and simulation point, and combine the height coordinate of known point, obtain mould
Intend the height coordinate of point, complete the simulation of dimensional topography, i.e. complete the modeling of landform.
Further, described interpolation method uses spline method.
Further, in step 3, described wind field condition entry to relief model will carry out combining of wind field and landform
Resolving, its concrete steps include
(1) the even direct current flow velocity that artificially generated terrain area inner measuring obtains is V∞, in the triangular mesh of artificially generated terrain, edge
Rectilinear vortex sections are arranged on three limits of each landform triangular mesh, form SPATIAL VORTEX ring element, can not penetrate bar in conjunction with object plane
Part, can obtain the relevant information arranging collar vortex on each triangle gridding；
(2) utilize the induced velocity formula of collar vortex, i.e. can get induced velocity V of space any pointPi, then in wind field
The wind speed of any point is this induced velocity and the linear superposition of even direct current speed, i.e.
WP=VPi+V∞
The most just the air speed value of space any point in artificially generated terrain region can be obtained.
The combining modeling method and can truly reappear landform, and the base of at this point shape simulation of a kind of landform of the present invention and wind field
Complete the calculating of wind field on plinth, the present invention only needs the mean wind speed inputting in topographic contour data and simulated domain to get final product real
Existing landform and the synchronously simulating of wind field；Need not move through first artificially generated terrain, again to landform model simplification, the most again calculate meet answering of field
Miscellaneous step, the simulation for simple/complicated landform and wind field is the most applicable, has simple structure, that amount of calculation is little, precision is high etc. is excellent
Point.
Accompanying drawing explanation
Accompanying drawing herein is merged in description and constitutes the part of this specification, it is shown that meet the enforcement of the present invention
Example, and for explaining the principle of the present invention together with description.
Fig. 1 is contour map and the given data extraction point schematic diagram of the present invention；
Fig. 2 is the artificially generated terrain stress and strain model schematic diagram of the present invention；
Fig. 3 is the dimensional topography numerical model schematic diagram obtained of the present invention；
Fig. 4 is the landform motion pattern with wind field of the present invention.
Detailed description of the invention
Clearer for the purpose making the present invention implement, technical scheme and advantage, below in conjunction with in the embodiment of the present invention
Accompanying drawing, the technical scheme in the embodiment of the present invention is further described in more detail.In the accompanying drawings, the most identical or class
As label represent same or similar element or there is the element of same or like function.Described embodiment is the present invention
A part of embodiment rather than whole embodiments.The embodiment described below with reference to accompanying drawing is exemplary type, it is intended to use
In explaining the present invention, and it is not considered as limiting the invention.Based on the embodiment in the present invention, ordinary skill people
The every other embodiment that member is obtained under not making creation type work premise, broadly falls into the scope of protection of the invention.Under
Face combines accompanying drawing and is described in detail embodiments of the invention.
In describing the invention, it is to be understood that term " " center ", " longitudinally ", " laterally ", "front", "rear",
The orientation of the instruction such as "left", "right", " vertically ", " level ", " top ", " end ", " interior ", " outward " or position relationship are for based on accompanying drawing institute
The orientation shown or position relationship, be for only for ease of and describe the present invention and simplify description rather than instruction or the dress of hint indication
Put or element must have specific orientation, with specific azimuth configuration and operation, therefore it is not intended that protect the present invention
The restriction of scope.
The dimensional topography in this region and corresponding wind field is simulated as a example by the contour map of Fig. 1 is as simulated target.
First, set up the coordinate system that artificially generated terrain region is suitable, in the present embodiment in coordinate system with sea level as X/Y plane,
With height above sea level Z axis coordinate, coordinate system is facilitated to match with the coordinate system of contour map.Special to ground based on equal pitch contour characteristic type
Levy and extract, the information provided according to contour map, if choosing n equal pitch contour node (point in Fig. 1) altogether, and will be all
The coordinate position of node two-dimensional matrix XkRepresenting, height uses column vector HkRepresent.Point in Fig. 1 is known point, obtains
Know plane coordinates a little and height coordinate.
Then, taking the 5000m × 2300m corresponding with actual landform is artificially generated terrain region, obtains artificially generated terrain region
Planar point and coordinate, for terrain simulation analog data or simulation point；And use triangular mesh to generate method, it is configured to
The triangular mesh of artificially generated terrain.
It is analog data or the simulation point constructing triangular mesh with the point in Fig. 2, trianglees for three simulation point compositions
Grid is the triangle element mesh of terrain simulation.
Afterwards, it is known that n equal pitch contour node coordinate XkCoordinate X with m terrain mesh node (or claiming simulation point)u, utilize
Morphological interpolation algorithm (concrete for spline method) determines interpolating matrix G between the twouk, i.e. analog data or simulation point with
Modality relation between given data or known point.Thus can obtain, with interpolation, the height that all terrain mesh nodes are corresponding
Data Hu, it is embodied as
Hu=GukHk
Terrain mesh coordinate data X obtaineduWith altitude information Hu.By visual software to actual landform modeling result
Carry out three dimensional display, as shown in Figure 3.
Finally, the even direct current flow velocity that this artificially generated terrain area measure obtains is V∞=[u0 v0 w0]T.With completed
Based on shape simulation lattice, rectilinear vortex section is arranged on three limits along each landform triangular mesh, forms SPATIAL VORTEX ring element, knot
Compound face can not penetrate condition, the relevant information arranging collar vortex on each triangle gridding can be obtained, relevant information is whirlpool
The position coordinates of ring and the intensity of collar vortex.
Wherein, Г is the intensity of collar vortex, AAICFor normal direction induced velocity coefficient matrix, A0It is and unit normal direction and remote front
Flow relevant constant coefficient array.
Utilize the induced velocity formula of collar vortex, i.e. can get induced velocity V of space any pointPi, then appoint in wind field
Anticipate any wind speed is this induced velocity and the linear superposition of even direct current speed, i.e.
WP=VPi+V∞
The most just the air speed value of space any point in this simulated domain can be obtained.Just can be shown by visual software
The motion pattern (as shown in Figure 4) of this region wind field is shown.
Terrain simulation is blended by landform and the associating modeling method of wind field that the present invention provides first with simulation of wind, opens
The new approaches simultaneously completing terrain simulation with wind field reproduction are warded off；The method of the present invention is applicable to the simulation of simple or complicated landform
And wind field reproduction, solve engineer applied mesorelief and the separate difficult problem of simulation of wind, not only increase terrain simulation
Verity, improves the accuracy of wind field reproduction simultaneously；It addition, the Method Modeling parameter acquiring of the present invention is simple, it is easy to solve,
And wind field distribution the most accurately can be given, meet engineering actual demand.By the phantom of the method foundation of the present invention not
It is only the analysis model that flight dynamics research provides, can be also loading analysis, the offer phantoms such as steady analysis are provided；Set up
Phantom can be used in flight dynamics research, and aircraft is the research of kinetic characteristic, the most expansible use in the complexity wind field of low latitude
In research fields such as dynamic response analysis, gust alleviation, prominent wind response analyses；The present invention is the research of follow-up complicated wind field modeling,
Low altitude flight performances for research aircraft etc. provide reliable environmental simulation basis, also extend to similar aviation simulator simultaneously
In analogue simulation field.
The above, the only optimum detailed description of the invention of the present invention, but protection scope of the present invention is not limited thereto,
Any those familiar with the art in the technical scope that the invention discloses, the change that can readily occur in or replacement,
All should contain within protection scope of the present invention.Therefore, protection scope of the present invention should be with the protection model of described claim
Enclose and be as the criterion.
Claims (7)
1. a landform and wind field combine modeling method, it is characterised in that including:
Step one: landform contour map is carried out data acquisition；
Step 2: carry out terrain modeling according to the data collected；
Step 3: by wind field condition entry to relief model, carry out the Combined Calculation of landform and wind field.
Landform the most according to claim 1 combines modeling method with wind field, it is characterised in that in step one, to landform
Contour map carry out data collection steps and include:
1) coordinate system is set up, with the contour map region of described landform for artificially generated terrain region；
2) by predetermined way, contour map being carried out data acquisition, the data of collection include plane coordinates and altitude data, described
Altitude data is height coordinate, and described data are referred to as known point.
Landform the most according to claim 2 combines modeling method with wind field, it is characterised in that with sea in described coordinate system
Plane is X/Y plane, with height above sea level for Z axis coordinate.
Landform the most according to claim 3 combines modeling method with wind field, it is characterised in that described predetermined way includes
Constant duration data acquisition and equidistantly height gather.
Landform the most according to claim 4 combines modeling method with wind field, it is characterised in that in step 2, landform is built
Mould step includes:
1) in artificially generated terrain region, the point being parallel to X/Y plane being carried out random acquisition, the point collected is as simulation point；
2) according to the plane coordinates of simulation point, triangular mesh is used to generate method, the triangular mesh of constructing analog landform；
3) according to plane coordinates and the altitude data of the known point collected under coordinate system, and the panel data of simulation point, utilize
Interpolation method, obtains the modality relation between known point and simulation point, and combines the height coordinate of known point, obtains simulating point
Height coordinate, complete the simulation of dimensional topography, i.e. complete the modeling of landform.
Landform the most according to claim 5 combines modeling method with wind field, it is characterised in that described interpolation method uses
It it is spline method.
Landform the most according to claim 6 combines modeling method with wind field, it is characterised in that in step 3, described general
Wind field condition entry is to carrying out the Combined Calculation of wind field and landform in relief model, and its concrete steps include:
(1) the even direct current flow velocity that artificially generated terrain area inner measuring obtains is V∞, in the triangular mesh of artificially generated terrain, along each
Rectilinear vortex section is arranged on three limits of landform triangular mesh, forms SPATIAL VORTEX ring element, in conjunction with the condition that can not penetrate of object plane, i.e.
The relevant information arranging collar vortex on each triangle gridding can be obtained；
(2) utilize the induced velocity formula of collar vortex, i.e. can get induced velocity V of space any pointPi, then in wind field arbitrarily
The wind speed of any is this induced velocity and the linear superposition of even direct current speed, i.e. WP=VPi+V∞
The most just the air speed value of space any point in artificially generated terrain region can be obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610374228.6A CN106055788A (en) | 2016-05-31 | 2016-05-31 | Combined modeling method of terrain and wind field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610374228.6A CN106055788A (en) | 2016-05-31 | 2016-05-31 | Combined modeling method of terrain and wind field |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106055788A true CN106055788A (en) | 2016-10-26 |
Family
ID=57172216
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610374228.6A Pending CN106055788A (en) | 2016-05-31 | 2016-05-31 | Combined modeling method of terrain and wind field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106055788A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107170044A (en) * | 2017-05-09 | 2017-09-15 | 福州大学 | A kind of dynamic and visual method of the wind based on dimensional topography |
| CN108664705A (en) * | 2018-04-13 | 2018-10-16 | 华中科技大学 | A method of the simulation complicated landform roughness of ground surface based on OpenFOAM |
| CN109033181A (en) * | 2018-06-26 | 2018-12-18 | 河海大学 | A kind of complicated landform area wind field geography method for numerical simulation |
| CN110111420A (en) * | 2019-04-12 | 2019-08-09 | 东北林业大学 | Skifield wind field modeling method |
| CN110136255A (en) * | 2019-04-12 | 2019-08-16 | 东北林业大学 | A kind of skifield sportsman windage energy consumption appraisal procedure |
| CN113160409A (en) * | 2021-02-26 | 2021-07-23 | 中国华能集团清洁能源技术研究院有限公司 | Complex terrain modeling method and system suitable for refined wind field simulation |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104361616A (en) * | 2014-11-05 | 2015-02-18 | 南车株洲电力机车研究所有限公司 | Terrain and landform document acquisition method for wind power plant wind resource assessment |
-
2016
- 2016-05-31 CN CN201610374228.6A patent/CN106055788A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104361616A (en) * | 2014-11-05 | 2015-02-18 | 南车株洲电力机车研究所有限公司 | Terrain and landform document acquisition method for wind power plant wind resource assessment |
Non-Patent Citations (3)
| Title |
|---|
| 吉柏锋: "风场数值模拟中复杂地形的ＮＵＲＢＳ建模", 《武汉理工大学学报》 * |
| 周志勇等: "大范围区域复杂地形风场数值模拟研究", 《力学季刊》 * |
| 唐矗等: "基于地形高程数据的复杂地形风场建模方法", 《北京航空航天大学学报》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107170044A (en) * | 2017-05-09 | 2017-09-15 | 福州大学 | A kind of dynamic and visual method of the wind based on dimensional topography |
| CN107170044B (en) * | 2017-05-09 | 2019-09-13 | 福州大学 | A kind of dynamic and visual method of the wind based on dimensional topography |
| CN108664705A (en) * | 2018-04-13 | 2018-10-16 | 华中科技大学 | A method of the simulation complicated landform roughness of ground surface based on OpenFOAM |
| CN109033181A (en) * | 2018-06-26 | 2018-12-18 | 河海大学 | A kind of complicated landform area wind field geography method for numerical simulation |
| CN109033181B (en) * | 2018-06-26 | 2020-06-30 | 河海大学 | Wind field geographic numerical simulation method for complex terrain area |
| CN110111420A (en) * | 2019-04-12 | 2019-08-09 | 东北林业大学 | Skifield wind field modeling method |
| CN110136255A (en) * | 2019-04-12 | 2019-08-16 | 东北林业大学 | A kind of skifield sportsman windage energy consumption appraisal procedure |
| CN110111420B (en) * | 2019-04-12 | 2023-08-01 | 东北林业大学 | Ski field wind field modeling method |
| CN110136255B (en) * | 2019-04-12 | 2023-08-01 | 东北林业大学 | Wind resistance and energy consumption evaluation method for skiing field athlete |
| CN113160409A (en) * | 2021-02-26 | 2021-07-23 | 中国华能集团清洁能源技术研究院有限公司 | Complex terrain modeling method and system suitable for refined wind field simulation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106055788A (en) | Combined modeling method of terrain and wind field | |
| Uccellini et al. | The coupling of upper and lower tropospheric jet streaks and implications for the development of severe convective storms | |
| CN106441311B (en) | A kind of non-cooperative Spacecraft relative pose measurement method based on laser imaging radar | |
| KR101332559B1 (en) | Numerical simulation system and method for atmospheric wind flow by computational fluid dynamics | |
| CN105513133A (en) | Method for making and displaying urban wind environment digital map | |
| CN113868971B (en) | Airport area three-dimensional refined wind field reconstruction method based on numerical simulation model and historical wind field characteristics | |
| CN113687447B (en) | Local area wind field monitoring method based on various wind measuring devices | |
| KR20120101949A (en) | System and method for wind field realizing of weather diagnostic model and diffusion model | |
| CN104778331A (en) | Spatial interpolation method for long-span bridge monitoring data | |
| CN114580310A (en) | Method for realizing scale reduction processing of WRF (hand-wrenching simulation) wind field based on PALM (PALM fiber laser) | |
| Willenborg et al. | Integration of semantic 3D city models and 3D mesh models for accuracy improvements of solar potential analyses | |
| CN107102339A (en) | A kind of automatic generation method of above-ground route data | |
| CN115994496A (en) | Urban park high resolution atmospheric CO 2 Numerical simulation method of concentration three-dimensional field | |
| CN105184667B (en) | The method of dual nested simulation wind farm wind velocity distribution | |
| CN103335608A (en) | Airborne LiDAR three-dimensional data acquisition method for establishing three-dimensional digital power transmission and transformation grid | |
| CN109033181B (en) | Wind field geographic numerical simulation method for complex terrain area | |
| CN101750616A (en) | Method for measuring wind resistance to vegetation and the system thereof | |
| CN115795402B (en) | Variational method-based multi-source precipitation data fusion method and system | |
| CN104573264A (en) | Method for simulating low-attitude wind shear area crossing of aircraft | |
| Zhang et al. | Voxel-based urban vegetation volume analysis with LiDAR point cloud | |
| Wei et al. | Elevation data fitting and precision analysis of Google Earth in road survey | |
| Haber | Eco-environmental surface modelling requires integration of both extrinsic and intrinsic informations | |
| Sun et al. | Mosaic research with 3D lidar point cloud of civil aircraft | |
| KR101089256B1 (en) | Geographic Information analysis method using of programing Invisible Depth Analysis | |
| CN116910691B (en) | Earthwork efficiency evaluation method based on influence conditions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161026 |