CN109703770A - Based on the carrier-borne machine aided of anemometry laser radar and CFD database, method drops - Google Patents
Based on the carrier-borne machine aided of anemometry laser radar and CFD database, method drops Download PDFInfo
- Publication number
- CN109703770A CN109703770A CN201811519169.2A CN201811519169A CN109703770A CN 109703770 A CN109703770 A CN 109703770A CN 201811519169 A CN201811519169 A CN 201811519169A CN 109703770 A CN109703770 A CN 109703770A
- Authority
- CN
- China
- Prior art keywords
- velocity
- cfd
- laser radar
- carrier
- anemometry laser
- 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
Landscapes
- Traffic Control Systems (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
Method and system are dropped based on the carrier-borne machine aided of anemometry laser radar and CFD database the invention discloses a kind of.The present invention is based on the CFD models of anemometry laser radar amendment target naval vessels, when establishing CFD database, only retain the wind speed information of carrier-borne aircraft glissade near zone, the optimum data table of CFD database is matched according to the radial velocity of anemometry laser radar real-time measurement, and then the wind speed information of carrier-borne aircraft drop zone is obtained, auxiliary carrier-borne aircraft landing.Result can not be slowly provided in real time the present invention overcomes existing CFD calculating and the shortcomings that anemometry laser radar can only measure radial velocity, with precision of prediction height, stability is good, predicted time complexity is low and the good characteristics such as spatial resolution height, it can be applied to carrier landing decision and control system, improve carrier-borne the accuracy and safety of warship.
Description
Technical field
The present invention relates to computational fluid dynamics field (CFD), aerodynamic database and anemometry laser radar fields, especially
Be related to combine anemometry laser radar measured data and CFD database estimation carrier-borne aircraft glissade wind speed carrier-borne machine aided drop method and
System.
Background technique
Some large-scale water surface naval vessels can for carrier-borne aircraft takeoff and landing, primary challenge strength of the carrier-borne aircraft as naval vessels,
Its key technology be how to ensure the safe falling under very rugged environment, if not can solve safety and precise warship problem,
Naval vessels just lose due fighting capacity, wherein the interference of naval vessels wake flow be influence carrier-borne aircraft succeed warship principal element it
One, horizontal disturbance influences aircraft airspeed and pitching movement, vertical disturbing influence flying height, side direction perturbation are easy to cause aircraft side
Rolling, therefore obtain the crucial problem that each velocity component of wake flow is carrier landing.
Research for large-scale naval vessels wake flow at present relies primarily on wind tunnel test, numerical simulation and real ship and directly measures three kinds
Mode.
The present inventor has found after study: wind tunnel test can save a large amount of expense, and can control more
The test environment of change is the important means of Flow Field Distribution around research naval vessel, but the dimensional effect of scale model and wind-tunnel wall
Face will influence test result.And local flow field details can not be captured, therefore results of wind tunnel and real flow of ship field exist
Gap.Traditional reality ship measurement method, if hot line is popped one's head in, anemobiagraph somewhat expensive, safety is bad and measurement range is very limited.
Doppler anemometry laser radar is acknowledged as the most effectual way of atmospheric wind remote sensing, with the high and low empty nothing of spatial and temporal resolution
Blind area, measurement accuracy is high, Electro Magnetic Compatibility is good, round the clock continuous observation, can be achieved it is excellent from ground to 110km height all standing
Gesture, but it is directed to this kind of specific question of carrier landing, the distance resolution of anemometry laser radar is still not high enough at present, and
Can only real-time measurement radial velocity, the component velocity in three directions can not be directly obtained, not can be used directly carrier-borne machine control system.
The differential equation that computational fluid dynamics (CFD) is flowed by solving control fluid obtains the flow field of fluid flowing in continuum
Discrete distribution on domain has at low cost thus approximate simulation fluid mobility status, and information is complete, distance resolution is high and can mould
In quasi- real process the advantages of various states, with the development of high-performance computer technology (HPC), Pneumatic Calculation has also become to mention
For the main means of flow numerical simulation data, it is widely used in the research of ship deck airflow field, and its numerical value calculates
Abundant verifying has been obtained in result reliability, but since CFD is the thought based on numerical discretization, computationally intensive, Wu Fashi
When provide analog result, therefore not can be used directly yet in carrier-borne aircraft in real time warship process.Therefore, in order to which carrier-borne aircraft succeeds
Warship is badly in need of a set of auxiliary system that can be realized accurate acquisition glissade wind speed information in real time.
Summary of the invention
In order to overcome the above-mentioned deficiencies of the prior art, the object of the present invention is to provide one kind based on anemometry laser radar and
System drops in the carrier-borne machine aided of CFD database, the radial velocity that can be measured according to anemometry laser radar, real in conjunction with CFD database
When output carrier-borne aircraft glissade on horizontal velocity, side velocity and vertical speed, provide speed in time for carrier-borne machine control system
Information is spent, and provides reference by prediction flight path for pilot and signal official, auxiliary carrier-borne aircraft is capable of the drop of safety and precise
It falls.
The present invention is implemented as follows:
It is a kind of that method is dropped based on the carrier-borne machine aided of anemometry laser radar and CFD database, comprising:
The Computational Fluid Dynamics model of target naval vessels is established according to the geometric parameter of target naval vessels；
The CFD model is corrected according to the air speed data of anemometry laser radar measurement；The anemometry laser radar setting exists
On target naval vessels；
It selects scheduled parameter to carry out discretization to the boundary condition of revised CFD model, establishes and dropped for carrier-borne aircraft
Fall on the CFD database of target naval vessels；
Obtain the radial velocity in the carrier-borne aircraft glissade direction of anemometry laser radar real-time measurement；
The radial velocity is matched with the tables of data in the CFD database of foundation, obtains the optimum number being matched to
According to table；It include multiple component velocities in the optimum data table；The component velocity includes horizontal velocity, side velocity and vertical speed
Degree；
According to the ship's speed that GPS module built in anemometry laser radar exports, the earth is converted by the multiple component velocity and is referred to
The earth component velocity under system；
The earth component velocity is fed back to the flight control system of carrier-borne aircraft.
Further, described to include: according to the air speed data amendment CFD model of anemometry laser radar measurement
Radial wind speed and CFD model calculated result that anemometry laser radar measurement obtains are compared into verifying, according to right
Than computational domain, mesh quality, boundary condition and the turbulence model in verification result adjustment CFD model, revised CFD mould is obtained
Type.
Further, the scheduled parameter section of selection carries out discretization to the boundary condition of revised CFD model,
Include:
The size and Orientation of deck wind is selected to carry out discretization to the boundary condition of revised CFD model；The deck
Wind is the speed of front incoming flow relative target naval vessels, reduces and calculates cost.Deck wind stream field influence of spatial distribution is maximum；Front
The direction advanced for target naval vessels.
Further, the downslide in target area when establishing CFD database, near batch signatures carrier-borne aircraft glissade
Line wind speed information, and stored according to reference format, reduce memory space.
Glissade wind speed information includes 11 variables, is followed successively by the elevation angle, deflection angle, deck wind size, deck wind direction, diameter
It is sat to velocity magnitude, horizontal component velocity, lateral component velocity, vertical component velocity, horizontal direction coordinate, lateral coordinate and vertical direction
Mark.
It is further, described to match the radial velocity with the tables of data in the CFD database of foundation, comprising:
According to the number for corresponding to the elevation angle and pitch angle in the elevation angle of anemometry laser radar and deflection angle index CFD database
According to table, the radial velocity is matched with the tables of data indexed.
Further, the ship's speed exported according to GPS module built in anemometry laser radar, converts the multiple component velocity to
Before the earth component velocity under earth reference system, further includes: the radial velocity based on anemometry laser radar real-time measurement, according to pre-
If Optimized model multiple component velocities in optimum data table are optimized；
Optimized model are as follows:
x0=[ui,vi,wi]；
Ub=[ui+ub,vi+vb,wi+wb]；
Lb=[ui-ub,vi-vb,wi-wb]；
Wherein, f (x) is optimization object function, and f (x) is each measuring point radial velocity measurement value and the optimum data table diameter
To the sum of residual absolute value of the difference of speed, minf (x) is to minimize to f (x)；Vi eFor the radial direction of anemometry laser radar real-time measurement
Speed,For the elevation angle, θ is deflection angle, uiHorizontal velocity, v for measuring pointiFor the side velocity of measuring point, wiFor the vertical speed of measuring point
Degree, ub are the coboundary of optimized variable, and lb is the lower boundary of optimized variable, ubFor preset horizontal velocity maximum variation, vbFor
Preset side velocity maximum variation, wbFor preset vertical speed maximum variation.
Correspondingly, a kind of drop system based on the carrier-borne machine aided of anemometry laser radar and CFD database, comprising:
CFD model establishes module, for establishing the Fluid Mechanics Computation of target naval vessels according to the geometric parameter of target naval vessels
CFD model；
CFD model correction module, the air speed data for being measured according to anemometry laser radar correct the CFD model；Institute
Anemometry laser radar is stated to be arranged on target naval vessels；
CFD Database module, for select scheduled parameter to the boundary condition of revised CFD model carry out from
Dispersion establishes the CFD database that target naval vessels are dropped to for carrier-borne aircraft；
Radial velocity obtains module, the radial direction in the carrier-borne aircraft glissade direction for obtaining anemometry laser radar real-time measurement
Speed；
Optimum data table obtains module, for carrying out the tables of data in the CFD database of the radial velocity and foundation
Matching, obtains the optimum data table being matched to；It include multiple component velocities in the optimum data table；The component velocity includes level
Speed, side velocity and vertical speed；
The earth component velocity obtains module, the ship's speed for being exported according to GPS module built in anemometry laser radar, will be described more
A component velocity is converted into the earth component velocity under earth reference system.
Drop module is helped, for the earth component velocity to be fed back to the flight control system of carrier-borne aircraft.
Further, further includes:
Optimization module, for the radial velocity based on anemometry laser radar real-time measurement, according to preset Optimized model pair
Multiple component velocities in optimum data table optimize；
Optimized model are as follows:
x0=[ui,vi,wi]；
Ub=[ui+ub,vi+vb,wi+wb]；
Lb=[ui-ub,vi-vb,wi-wb]；
Wherein, f (x) is optimization object function, and f (x) is each measuring point radial velocity measurement value and the optimum data table diameter
To the sum of residual absolute value of the difference of speed, Vi eFor the radial velocity of anemometry laser radar real-time measurement,For the elevation angle, θ is deflection angle,
uiHorizontal velocity, v for measuring pointiFor the side velocity of measuring point, wiFor the vertical speed of measuring point, ub is the coboundary of optimized variable,
Lb is the lower boundary of optimized variable, ubFor preset horizontal velocity maximum variation, vbFor the variation of preset side velocity maximum
Amount, wbFor preset vertical speed maximum variation.
It is provided by the invention that method and system are dropped based on the carrier-borne machine aided of anemometry laser radar and CFD database, it overcomes
CFD calculating can not slowly provide result and in real time the shortcomings that anemometry laser radar can only measure radial velocity, have precision of prediction
Height, stability is good, and predicted time complexity is low and the good characteristics such as spatial resolution height, can be applied to carrier landing and determines
Plan and control system assist and improve carrier-borne the accuracy and safety of warship.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill in field, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is provided in an embodiment of the present invention based on the carrier-borne machine aided of anemometry laser radar and CFD database drop method
Flow chart；
Fig. 2 is provided in an embodiment of the present invention based on the carrier-borne machine aided of anemometry laser radar and CFD database drop system
Flow chart；
Fig. 3 is the method for building up and data structure schematic diagram of CFD air speed data provided in an embodiment of the present invention；
Fig. 4 is System Working Principle schematic diagram provided in an embodiment of the present invention.
Specific embodiment
With reference to the attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based on this
The embodiment of invention, every other implementation obtained by those of ordinary skill in the art without making creative efforts
Example, belongs to protection scope of the present invention.
Embodiment:
Fig. 1 is provided in an embodiment of the present invention based on the carrier-borne machine aided of anemometry laser radar and CFD database drop method
Flow chart；As shown in Figure 1, a kind of drop method based on the carrier-borne machine aided of anemometry laser radar and CFD database, it is characterised in that:
Include:
Step S1, the Computational Fluid Dynamics model of target naval vessels is established according to the geometric parameter of target naval vessels.
Wherein, target naval vessels are used for for carrier-borne aircraft takeoff and landing.
Step S1 is specifically included: the geometric parameter of target naval vessels is obtained, it is soft using professional grid dividing according to geometric parameter
Part obtains the grid file of target naval vessels, tentatively establishes the Computational Fluid Dynamics model of target naval vessels.
CFD is the abbreviation of computational fluid dynamics (Computational Fluid Dynamics), be hydrodynamics and
The emerging cross discipline that computer science mutually merges, it quickly calculates energy from calculation method, using computer
Power obtains the approximate solution of fluid governing equation.
By the model foundation target ship for choosing reasonable mesh quality, the setting of boundary condition type, computational domain and turbulent flow
The CFD model of ship.
Step S2, the CFD model is corrected according to the air speed data of anemometry laser radar measurement；The anemometry laser radar
It is arranged on target naval vessels.
Anemometry laser radar, also referred to as Doppler anemometry laser radar, LDV technique, it is distant to be acknowledged as atmospheric wind
The most effectual way of sense.Anemometry laser radar is to atmospheric emission laser pulse (ultraviolet to infrared), with atmospheric interaction, optics
Telescope inputs optical receiver after collecting atmospheric aerosol particle and atmospheric molecule backscatter signal, is swashed by analysis transmitting
The radial Doppler frequency shift of light carrys out Wind Speed Inversion.It is simultaneous with the high and low empty non-blind area of spatial and temporal resolution, measurement accuracy height, electromagnetism
Capacitive is good, round the clock continuous observation, can be achieved from ground to the advantage of 110km height all standing.
The scan pattern of anemometry laser radar includes: PPI: constant zenith angle mode (azimuthal variation)；RHI: constant side
Parallactic angle mode (elevation angle variation)；DBS: Vertical Profile；LOS: fixed position continuous observation.
In one embodiment of the invention, anemometry laser radar can be scanned with PPI, i.e., constant zenith angle Mode scans,
Obtain the radial wind speed information of measurement direction.
Specifically, including: according to the air speed data amendment CFD model of anemometry laser radar measurement described in step S2
Radial wind speed and CFD model calculated result that anemometry laser radar measurement obtains are compared into verifying, wind is surveyed and swashs
The radial air speed data that optical radar measurement obtains includes radial wind speed size and corresponding measurement angle, location information.According to right
Than computational domain, mesh quality, boundary condition and the turbulence model in verification result adjustment CFD model, revised CFD mould is obtained
Type.The consistency of revised CFD model and the radial wind speed that anemometry laser radar measurement obtains is preferable.
Using the real ship measurement data of anemometry laser radar, the CFD model tentatively established is modified and is verified, it can
Guarantee the reasonability and accuracy of revised CFD model.
Step S3, it selects scheduled parameter to carry out discretization to the boundary condition of revised CFD model, establishes and be directed to warship
Carrier aircraft drops to the CFD database of target naval vessels.
In one embodiment, all be likely to occur is calculated by discrete boundary condition using the CFD model in step S2
Target naval vessels tail flow field under inlet flow conditions extracts the coordinate of carrier-borne aircraft glissade near zone and wind speed letter in each operating condition
Breath constructs small-sized wind speed database according to standard data format.
In one embodiment, perimeter strip of the scheduled parameter section of selection to revised CFD model described in step S3
Part carries out discretization, comprising:
Select the size and Orientation of the maximum deck wind of stream field influence of spatial distribution to the boundary of revised CFD model
Condition carries out discretization；The deck wind is the speed of front incoming flow relative target naval vessels.
When establishing CFD database, when progress boundary condition is discrete, the maximum first of stream field influence of spatial distribution is selected
The size and Orientation progress of plate wind (wind speed of relative target naval vessels referential) is discrete, can greatly reduce calculating operating condition, save
Calculating cost.Secondly the wind speed information of carrier-borne aircraft near zone is only remained, the size of data of single operating condition is only in 10kB
Left and right, for the size of population of CFD database in 100,000,000 magnitudes, this greatly reduces greatly carrying cost and index time.
Glissade wind speed letter when establishing CFD database, in the target area near batch signatures carrier-borne aircraft glissade
Breath, and stored according to reference format；The glissade is the landing track that carrier-borne aircraft drops to target naval vessels.
Glissade wind speed information includes 11 variables, is followed successively by the elevation angle, deflection angle, deck wind size, deck wind direction, diameter
It is sat to velocity magnitude, horizontal component velocity, lateral component velocity, vertical component velocity, horizontal direction coordinate, lateral coordinate and vertical direction
Mark.
After having corrected CFD database, so that it may according to the radial velocity of laser radar real-time measurement and revised CFD number
Speed reference information is provided according to library for carrier-borne aircraft landing；Specifically include step S4-S7.
Step S4, the radial velocity in the carrier-borne aircraft glissade direction of anemometry laser radar real-time measurement is obtained.
The anemometry laser radar realizes pitching movement and beat movement by two servo motors, can export scanning in real time
The pitch angle and deflection angle of line, built-in GPS module is for recording and exporting real-time ship's speed；The ship's speed is target naval vessels
Movement velocity.
The anemometry laser radar mounted in target naval vessels larboard Fresnel lens optics help drop system from steady platform, both
It can guarantee that laser radar light beam is not influenced by warship body is wiggly, and can guarantee that laser scanning line most connects with carrier-borne aircraft glissade
Closely.
The measurement direction of anemometry laser radar includes carrier-borne aircraft glissade direction, i.e., carrier-borne aircraft warship direction.
Step S5, the radial velocity is matched with the tables of data in the CFD database of foundation, obtains and is matched to
Optimum data table；It include multiple component velocities in the optimum data table；The component velocity includes horizontal velocity, side velocity and erects
Straight speed.
In one embodiment, in step S5, the tables of data by the radial velocity and the CFD database of foundation
It is matched, comprising:
According to the number for corresponding to the elevation angle and pitch angle in the elevation angle of anemometry laser radar and deflection angle index CFD database
According to table, the radial velocity is matched with the tables of data indexed.
Preferably, radial velocity is carried out based on the anemometry laser radar alignment carrier-borne aircraft isogonism glissade in step S4 to sweep
It retouches, according to the radial velocity of the pitch angle of the output of the servo-system of anemometry laser radar, pivot angle and real-time measurement in miniature number
According to matching in library with the immediate tables of data of measured data as optimum data table, and extract corresponding optimum level speed, side
To speed and vertical speed.
Step S6, the ship's speed exported according to GPS module built in anemometry laser radar, converts the multiple component velocity to greatly
The earth component velocity under ground referential.
It in a preferred embodiment, further include step S7, the flight control that the earth component velocity is fed back to carrier-borne aircraft
System.
In a preferred embodiment, the ship's speed that step S6 is exported according to GPS module built in anemometry laser radar, by institute
Before stating the earth component velocity that multiple component velocities are converted under earth reference system, further includes:
Step S8, the radial velocity based on anemometry laser radar real-time measurement, according to preset Optimized model to optimum number
It is optimized according to multiple component velocities in table；
Optimized model are as follows:
x0=[ui,vi,wi]；
Ub=[ui+ub,vi+vb,wi+wb]；
Lb=[ui-ub,vi-vb,wi-wb]；
Wherein, f (x) is optimization object function, and f (x) is each measuring point radial velocity measurement value and the optimum data table diameter
To the sum of residual absolute value of the difference of speed, Vi eFor the radial velocity of anemometry laser radar real-time measurement,For the elevation angle, θ is deflection angle,
uiHorizontal velocity, v for measuring pointiFor the side velocity of measuring point, wiFor the vertical speed of measuring point, ub is the coboundary of optimized variable,
Lb is the lower boundary of optimized variable, ubFor preset horizontal velocity maximum variation, vbFor the variation of preset side velocity maximum
Amount, wbFor preset vertical speed maximum variation.
Based on matched in step S5 come optimum data table in best component velocity and radar real-time measurement radial speed
Degree, using optimization algorithm, further corrects component velocity, excellent by measuring point with the wind speed of external position except radar measuring point on glissade
Change result interpolation to obtain.
Based on the ship's speed information that glissade wind speed information revised in step S6 and GPS provide, earth reference system is obtained
Under glissade wind speed and feed back to carrier-borne machine control system, help carrier-borne aircraft correct state of flight in time, inhibit the shadow of wake flow
It rings.
In one embodiment, method further include: the earth component velocity on three directions is inputted into carrier-borne dynamics side
Cheng Zhong further predicts that the flight path of carrier-borne aircraft, real-time display on a user interface, provide for pilot and signal official's decision
With reference to.
As shown in Fig. 2, correspondingly, the present invention also provides a kind of carrier-borne based on anemometry laser radar and CFD database
System drops in machine aided, comprising:
CFD model establishes module, for establishing the Fluid Mechanics Computation of target naval vessels according to the geometric parameter of target naval vessels
CFD model；
CFD model correction module, the air speed data for being measured according to anemometry laser radar correct the CFD model；Institute
Anemometry laser radar is stated to be arranged on target naval vessels；
CFD Database module, for select scheduled parameter to the boundary condition of revised CFD model carry out from
Dispersion establishes the CFD database that target naval vessels are dropped to for carrier-borne aircraft；
Radial velocity obtains module, the radial direction in the carrier-borne aircraft glissade direction for obtaining anemometry laser radar real-time measurement
Speed；
Optimum data table obtains module, for carrying out the tables of data in the CFD database of the radial velocity and foundation
Matching, obtains the optimum data table being matched to；It include multiple component velocities in the optimum data table；The component velocity includes level
Speed, side velocity and vertical speed；
The earth component velocity obtains module, the ship's speed for being exported according to GPS module built in anemometry laser radar, will be described more
A component velocity is converted into the earth component velocity under earth reference system；
Drop module is helped, for the earth component velocity to be fed back to the flight control system of carrier-borne aircraft.
In one embodiment, the system also includes:
Optimization module, for the radial velocity based on anemometry laser radar real-time measurement, according to preset Optimized model pair
Multiple component velocities in optimum data table optimize；
Optimized model are as follows:
x0=[ui,vi,wi]；
Ub=[ui+ub,vi+vb,wi+wb]；
Lb=[ui-ub,vi-vb,wi-wb]；
Wherein, f (x) is optimization object function, and f (x) is each measuring point radial velocity measurement value and the optimum data table diameter
To the sum of residual absolute value of the difference of speed, Vi eFor the radial velocity of anemometry laser radar real-time measurement,For the elevation angle, θ is deflection angle,
uiHorizontal velocity, v for measuring pointiFor the side velocity of measuring point, wiFor the vertical speed of measuring point, ub is the coboundary of optimized variable,
Lb is the lower boundary of optimized variable, ubFor preset horizontal velocity maximum variation, vbFor the variation of preset side velocity maximum
Amount, wbFor preset vertical speed maximum variation.
In the following, method of the invention is described in detail with specific application scenarios.
According to target naval vessels contour structures, geometrical model is established, in order to reduce number of grid and computation complexity, in geometry
The main structure for retaining target naval vessels in modeling process, which ignores stream field, influences smaller local micro-structure.Target naval vessels tail flow field
Simulation belong to Flow Field outside simulation, the outer boundary of zoning theoretically should in the infinite point of target naval vessels periphery, but
It is practical to calculate, a limited computational domain can be used to substitute unlimited computational domain, when computational domain is more much bigger than model, example
Such as, the volume of computational domain is 500-1000 times of target ship hull product, can reduce meter while the accuracy for guaranteeing to calculate
Calculation amount.
Entire computational domain is handled using the hybrid grid mode that unstructured grid and structured grid combine, target ship
The warship body geometry of ship is complex, preferable using tetrahedron unstructured grid effect around, and to warship body surface surface grids
It is encrypted with accurate catch surface information of flow, remaining periphery is all made of hexahedron structure grid, for of interest
Carrier-borne aircraft glissade region use refined net, remaining is set with controlling grid in ten million magnitude in boundary condition using coarse grid
When setting, entrance is speed entrance (Velocity-inlet), and outlet and coboundary are set as pressure export (Pressure-
Outlet), ship surface and sea level are set as wall surface (Wall), zoning right boundary according to wind direction be set as speed entrance and
Pressure export, i.e. windward side are speed entrance, and leeward is set as pressure export.Ship air wake belongs to high reynolds number turbulent flow,
In numerical solution, turbulence model selection criteria k- ε model carries out stable state calculating, and pressure-velocity coupling terms use in the equation of momentum
SIMPLE algorithm, the equation of momentum and turbulent viscosity use Second-order Up-wind format, and pressure equation uses second order spatial discrete scheme.
When one timing of front inlet flow conditions, theoretically the spatial distribution of the pulsatile flow field of target naval vessels tail portion is unique,
Steady state values simulation is carried out to tail flow field using CFD model, the Flow Field Distribution under the inlet flow conditions can be obtained.
Before establishing CFD database, carries out Modifying model first with the ship trial data of anemometry laser radar and test
Card, concrete operations are to extract the radial velocity and test of the corresponding plane of radar PPI scanning (fixed elevation) in CFD analog result
Data comparison corrects CFD model by improving computational domain, mesh quality and boundary condition repeatedly, until CFD analog result with
Test result Integral-fit obtains preferably, wherein the radial velocity expression formula of CFD simulation are as follows:
Wherein Q is position coordinates of the laser radar under target naval vessels coordinate system, and H is radar measuring point in target naval vessels coordinate
The position of system, V are the conjunction speed of measuring point under target naval vessels coordinate system.
The method for building up of CFD air speed data of the present invention with data structure as shown in figure 3, first have to carry out inlet flow conditions it is discrete
Change.Speed, direction, temperature, viscosity and the pressure and other parameters of front incoming flow and the travel speed of target naval vessels can all influence warship
Wake flow distribution, wherein front speed of incoming flow size, direction and ship's speed are maximum to the influence of spatial distribution of tail flow field.Select deck wind
(speed of front incoming flow relative target naval vessels) size and wind direction this two parameters carry out discrete velocity entrance boundary condition, the earth
Referential lower tail flow velocity degree can be modified by the ship's speed that GPS is exported, and substantially reduce the complexity of speed entrance in this way,
Therefore it also saves memory space and calculates cost.
It, can be by the size of deck wind and side for all inlet flow conditions being likely to occur in coverage goal naval vessels driving process
To being divided into plurality of discrete operating condition according to scheduled interval and precision respectively.
When handling CFD calculated result, first obtain carrier-borne aircraft glissade angle, batch signatures carrier-borne aircraft glissade it is attached
The wind speed information in close-target region.In order to cover the region of all possible glissade, target area can be divided at equal intervals
It is more parts, generates multiple groups glissade.
Every glissade wind speed information includes 11 variables, be followed successively by the elevation angle, pivot angle, deck wind size, deck wind direction,
Radial velocity size, horizontal component velocity, lateral component velocity, vertical component velocity, horizontal direction coordinate, lateral coordinate and vertical direction
Coordinate.According to the computing capability of current high-performance calculation (HPC) server, the calculating of all operating conditions can be completed within one week,
CFD data base establishment is can be completed into according to reference format storage as shown in Figure 3 in all glissade wind speed informations.
Shown in carrier-borne machine aided drop System Working Principle schematic diagram 4 of the present invention, specific work process is as follows:
The Doppler anemometry laser radar realizes pitching and beat movement by 2 servo electricity, and built-in GPS module is used for
The motion conditions of record and output target naval vessels, wherein within 1 ', laser radar is not influencing the precision controlling of servo motor
Under the premise of Fresnel lens optics helps drop system to work, it is mounted on target naval vessels larboard Fresnel lens optics and helps drop system
From steady platform, it not only can guarantee that laser radar light beam was not influenced by warship body is wiggly, but also can guarantee laser scanning line and warship
Carrier aircraft glissade is closest.Laser radar at work, the elevation angle and beat by servo-control system control radar probe
Angle is to be directed at carrier-borne aircraft glissade.Laser signal exports radial velocity, the position letter of measurement after the processing of radar signal system
Breath.
In order to improve database index efficiency, measured data is watched according to laser radar first in matching database data
The tables of data for corresponding to the elevation angle and pitch angle in the elevation angle and deflection angle index data base of control system output is taken, is greatly reduced in this way
Data Matching range.
In Data Matching, with radial velocity, component velocity and the coordinate of tables of data D record current data table, calculate simultaneously
With the sum of residual absolute value of the measured data of record anemometry laser radar and radial velocity in current data table δ, when δ minimum pair
The tables of data answered is optimum data table, it may be assumed that
Wherein Vi eFor the radar measured data of i-th of measuring point, ViFor the database data of i-th of measuring point, n is that measuring point is total
Number.Then start to match next tables of data, when δ is less than the δ of current record, current data table information be transmitted to tables of data D,
Otherwise next tables of data is directly matched, is repeated the above steps, until matched all tables of data, the wind speed of tables of data D at this time
It is nearest with practical glissade wind speed.
In order to further increase the approximation ratio of component velocity Yu true wind field, optimization algorithm is called to repair component velocity
Just, wherein optimization object function f (x) be each measuring point radial velocity measurement value and tables of data D radial velocity residual absolute value it
It is each point component velocity with, optimized variable.The Optimized model of use is as follows:
x0=[ui,vi,wi]；
Ub=[ui+ub,vi+vb,wi+wb]；
Lb=[ui-ub,vi-vb,wi-wb]；
Wherein ui、viAnd wiThe respectively horizontal velocity of measuring point, side velocity and vertical speed, ub and lb are respectively to optimize
The up-and-down boundary of variable, ui、viAnd wiComponent velocity maximum variation, recommendation 3m/s, 1m/s and 1m/s.The optimization mould
Type is Multivariable Linear optimization problem, preferably simplex method, can be quickly obtained optimal solution.Further, laser radar GPS is utilized
The ship's speed information of output is modified the speed under target naval vessels coordinate system, thus the horizontal velocity under obtaining earth coordinates
(U), side velocity (V) and vertical speed (W).For the component velocity in three directions there are two types of effect, one is believe component velocity
Breath feeds back to carrier-borne aircraft flight control system, the corresponding control response for helping carrier-borne aircraft flight control system to make, to inhibit tail
The influence of stream.Another kind effect is that component velocity information is substituted into carrier-borne aircraft kinetics equation, predicts the flight path of carrier-borne aircraft, and
Real-time display provides reference on a user interface, for signal official and pilot's decision.
Result and anemometry laser radar can not slowly be provided in real time the present invention overcomes CFD calculating can only measure radial velocity
The shortcomings that, have precision of prediction high, stability is good, and predicted time complexity is low and the good characteristics such as spatial resolution height, can be with
Applied to carrier landing decision and control system, carrier-borne the accuracy and safety of warship is assisted and improved.
The present invention has following gain effect compared with prior art:
(1) according to the real ship measurement data of laser radar, target naval vessels CFD is corrected and is verified repeatedly, standard can be obtained
True reasonable CFD model.
(2) when establishing CFD database, when progress boundary condition is discrete, select stream field influence of spatial distribution maximum
The size and Orientation of deck wind (relative target naval vessels referential wind speed) carry out discrete, greatly reduce calculating operating condition, save
Calculating cost.Secondly the wind speed information of carrier-borne aircraft near zone is only remained, the size of data of single operating condition is only in 10kB
Effect, the size of population of CFD database greatly reduce carrying cost and index time in 100,000,000 magnitudes.
(3) anemometry laser radar be mounted on target naval vessels larboard Fresnel lens optics help drop system from steady platform,
Not only it can guarantee that laser radar light beam was not influenced by warship body is wiggly, but also can guarantee laser scanning line and carrier-borne aircraft glissade most
It is close.
(4) it is first defeated according to anemometry laser radar SERVO CONTROL when anemometry laser radar measured data is with CFD database matching
The elevation angle and deflection angle out indexes the data area that the elevation angle and pitch angle are corresponded in CFD database, reduces Data Matching range,
Improve recall precision.
(5) based on match come best component velocity and radar real-time measurement radial velocity, using optimization algorithm, to point
Speed is further corrected, and the accuracy of component velocity is further improved.
(6) component velocity is converted to using GPS the component velocity of earth reference system, can be applied to the posture system of flight control system
System, can also apply Flight Trajectory Prediction.
Here it must be noted that other unaccounted parts that the present invention provides all be for well known to those skilled in the art,
Title or function according to the present invention, those skilled in the art can find the document of related record, therefore not into one
Walk explanation.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
Within the technical scope of the present disclosure, any changes or substitutions that can be easily thought of by anyone skilled in the art,
It should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with the protection model of claims
Subject to enclosing.
Claims (10)
1. a kind of drop method based on the carrier-borne machine aided of anemometry laser radar and CFD database, it is characterised in that: include:
The Computational Fluid Dynamics model of target naval vessels is established according to the geometric parameter of target naval vessels；
The CFD model is corrected according to the air speed data of anemometry laser radar measurement；The anemometry laser radar is arranged in target
On naval vessels；
Scheduled parameter is selected to carry out discretization to the boundary condition of revised CFD model, foundation is dropped to for carrier-borne aircraft
The CFD database of target naval vessels；
Obtain the radial velocity in the carrier-borne aircraft glissade direction of anemometry laser radar real-time measurement；
The radial velocity is matched with the tables of data in the CFD database of foundation, obtains the optimum data table being matched to；
It include multiple component velocities in the optimum data table；The component velocity includes horizontal velocity, side velocity and vertical speed；
The earth component velocity the multiple component velocity converted under earth reference system.
2. the method according to claim 1, wherein the air speed data according to anemometry laser radar measurement is repaired
The just described CFD model includes:
Radial wind speed and CFD model calculated result that anemometry laser radar measurement obtains are compared into verifying, tested according to comparison
Computational domain, mesh quality, boundary condition and the turbulence model in result adjustment CFD model are demonstrate,proved, revised CFD model is obtained.
3. method according to claim 1 or 2, which is characterized in that the scheduled parameter section of selection is to revised
The boundary condition of CFD model carries out discretization, comprising:
The size and Orientation of deck wind is selected to carry out discretization to the boundary condition of revised CFD model；The deck wind is
The speed of front incoming flow relative target naval vessels.
4. according to the method in claim 2 or 3, which is characterized in that when establishing CFD database, batch signatures carrier-borne aircraft
Glissade wind speed information in target area near glissade, and stored according to reference format；
Glissade wind speed information includes 11 variables, is followed successively by the elevation angle, deflection angle, deck wind size, deck wind direction, radial speed
Spend size, horizontal component velocity, lateral component velocity, vertical component velocity, horizontal direction coordinate, lateral coordinate and vertical direction coordinate.
5. the method according to claim 1, wherein the CFD database by the radial velocity and foundation
In tables of data matched, comprising:
According to the data for corresponding to the elevation angle and pitch angle in the elevation angle of anemometry laser radar and deflection angle index CFD database
Table matches the radial velocity with the tables of data indexed.
6. the method according to claim 1, wherein converting the multiple component velocity under earth reference system
Before the earth component velocity, further includes: the radial velocity based on anemometry laser radar real-time measurement, according to preset Optimized model pair
Multiple component velocities in optimum data table optimize；
Optimized model are as follows:
x0=[ui,vi,wi]；
Ub=[ui+ub,vi+vb,wi+wb]；
Lb=[ui-ub,vi-vb,wi-wb]；
Wherein, f (x) is optimization object function, Vi eFor the radial velocity of anemometry laser radar real-time measurement,For the elevation angle, θ is inclined
Pivot angle, uiHorizontal velocity, v for measuring pointiFor the side velocity of measuring point, wiFor the vertical speed of measuring point, ub is the upper of optimized variable
Boundary, lb are the lower boundary of optimized variable, ubFor preset horizontal velocity maximum variation, vbIt is maximum for preset side velocity
Variation, wbFor preset vertical speed maximum variation.
7. the method according to claim 1, wherein converting the multiple component velocity under earth reference system
After the earth component velocity, further includes: the earth component velocity is fed back to the flight control system of carrier-borne aircraft.
8. a kind of drop system based on the carrier-borne machine aided of anemometry laser radar and CFD database, it is characterised in that: include:
CFD model establishes module, for establishing the computation fluid dynamics mould of target naval vessels according to the geometric parameter of target naval vessels
Type；
CFD model correction module, the air speed data for being measured according to anemometry laser radar correct the CFD model；The survey
Wind laser radar is arranged on target naval vessels；
CFD Database module, it is discrete for selecting scheduled parameter to carry out the boundary condition of revised CFD model
Change, establishes the CFD database for dropping to target naval vessels for carrier-borne aircraft；
Radial velocity obtains module, the radial speed in the carrier-borne aircraft glissade direction for obtaining anemometry laser radar real-time measurement
Degree；
Optimum data table obtains module, for the radial velocity to be matched with the tables of data in the CFD database of foundation,
Obtain the optimum data table being matched to；It include multiple component velocities in the optimum data table；The component velocity include horizontal velocity,
Side velocity and vertical speed；
The earth component velocity obtains module, the earth component velocity for converting the multiple component velocity under earth reference system.
9. system according to claim 8, which is characterized in that further include:
Optimization module, for the radial velocity based on anemometry laser radar real-time measurement, according to preset Optimized model to best
Multiple component velocities in tables of data optimize；
Optimized model are as follows:
x0=[ui,vi,wi]；
Ub=[ui+ub,vi+vb,wi+wb]；
Lb=[ui-ub,vi-vb,wi-wb]；
Wherein, f (x) is optimization object function, and f (x) is each measuring point radial velocity measurement value and the radial speed of the optimum data table
Spend the sum of residual absolute value of the difference, Vi eFor the radial velocity of anemometry laser radar real-time measurement,For the elevation angle, θ is deflection angle, uiFor
The horizontal velocity of measuring point, viFor the side velocity of measuring point, wiFor the vertical speed of measuring point, ub is the coboundary of optimized variable, and lb is
The lower boundary of optimized variable, ubFor preset horizontal velocity maximum variation, vbFor preset side velocity maximum variation, wb
For preset vertical speed maximum variation.
10. system according to claim 8, which is characterized in that further include:
Drop module is helped, for the earth component velocity to be fed back to the flight control system of carrier-borne aircraft.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811519169.2A CN109703770B (en) | 2018-12-12 | 2018-12-12 | Shipboard aircraft landing assisting method based on wind-finding laser radar and CFD database |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811519169.2A CN109703770B (en) | 2018-12-12 | 2018-12-12 | Shipboard aircraft landing assisting method based on wind-finding laser radar and CFD database |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109703770A true CN109703770A (en) | 2019-05-03 |
| CN109703770B CN109703770B (en) | 2021-08-06 |
Family
ID=66256413
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811519169.2A Active CN109703770B (en) | 2018-12-12 | 2018-12-12 | Shipboard aircraft landing assisting method based on wind-finding laser radar and CFD database |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109703770B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110263369A (en) * | 2019-05-10 | 2019-09-20 | 珠海市公共气象服务中心(珠海市防雷所)(珠海市突发事件预警信息发布中心) | Building surface wind resistance grade design method based on climatic analysis and numerical simulation |
| CN111898687A (en) * | 2020-08-03 | 2020-11-06 | 成都信息工程大学 | Radar reflectivity data fusion method based on Dilongnie triangulation |
| CN113961543A (en) * | 2021-10-25 | 2022-01-21 | 成都飞机工业（集团）有限责任公司 | Aerodynamic database generation method based on mgaero |
| CN115345091A (en) * | 2022-09-02 | 2022-11-15 | 北京瑞科同创能源科技有限公司 | CFD-based wind finding radar correction method and device, electronic equipment and storage medium |
| CN115408962A (en) * | 2022-11-02 | 2022-11-29 | 南京信息工程大学 | Wind field reconstruction method and system based on CFD simulation and wind lidar |
| CN117949972A (en) * | 2024-03-26 | 2024-04-30 | 中国人民解放军国防科技大学 | Laser radar-based method, system and device for measuring wind field of take-off and landing channel of airplane |
| CN117949972B (en) * | 2024-03-26 | 2024-06-11 | 中国人民解放军国防科技大学 | Laser radar-based method, system and device for measuring wind field of take-off and landing channel of airplane |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1124349A (en) * | 1995-11-08 | 1996-06-12 | 中国民用航空华东管理局 | Laser navigation aiding system and its using method |
| CN105260508A (en) * | 2015-09-16 | 2016-01-20 | 南京航空航天大学 | Method for predicting optimal release point of airdropped material |
| CN102610126B (en) * | 2010-12-20 | 2016-03-30 | 塞莱斯系统集成公司 | The rapid vertical trajectory predictions method of air traffic control and relevant ATM system |
| CN106405541A (en) * | 2016-11-14 | 2017-02-15 | 苏州途视电子科技有限公司 | Fully-coherent continuous-wave Doppler radar and distance measurement and velocity measurement method thereof |
| CN108974373A (en) * | 2018-07-19 | 2018-12-11 | 西安恒宇众科空间技术有限公司 | Based on binocular vision aircraft independent landing device |
-
2018
- 2018-12-12 CN CN201811519169.2A patent/CN109703770B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1124349A (en) * | 1995-11-08 | 1996-06-12 | 中国民用航空华东管理局 | Laser navigation aiding system and its using method |
| CN102610126B (en) * | 2010-12-20 | 2016-03-30 | 塞莱斯系统集成公司 | The rapid vertical trajectory predictions method of air traffic control and relevant ATM system |
| CN105260508A (en) * | 2015-09-16 | 2016-01-20 | 南京航空航天大学 | Method for predicting optimal release point of airdropped material |
| CN106405541A (en) * | 2016-11-14 | 2017-02-15 | 苏州途视电子科技有限公司 | Fully-coherent continuous-wave Doppler radar and distance measurement and velocity measurement method thereof |
| CN108974373A (en) * | 2018-07-19 | 2018-12-11 | 西安恒宇众科空间技术有限公司 | Based on binocular vision aircraft independent landing device |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110263369A (en) * | 2019-05-10 | 2019-09-20 | 珠海市公共气象服务中心(珠海市防雷所)(珠海市突发事件预警信息发布中心) | Building surface wind resistance grade design method based on climatic analysis and numerical simulation |
| CN110263369B (en) * | 2019-05-10 | 2023-12-29 | 珠海市公共气象服务中心(珠海市防雷所)(珠海市突发事件预警信息发布中心) | Building surface wind resistance grade design method based on climate analysis and numerical simulation |
| CN111898687A (en) * | 2020-08-03 | 2020-11-06 | 成都信息工程大学 | Radar reflectivity data fusion method based on Dilongnie triangulation |
| CN111898687B (en) * | 2020-08-03 | 2021-07-02 | 成都信息工程大学 | Radar reflectivity data fusion method based on Dilongnie triangulation |
| CN113961543A (en) * | 2021-10-25 | 2022-01-21 | 成都飞机工业（集团）有限责任公司 | Aerodynamic database generation method based on mgaero |
| CN115345091A (en) * | 2022-09-02 | 2022-11-15 | 北京瑞科同创能源科技有限公司 | CFD-based wind finding radar correction method and device, electronic equipment and storage medium |
| CN115408962A (en) * | 2022-11-02 | 2022-11-29 | 南京信息工程大学 | Wind field reconstruction method and system based on CFD simulation and wind lidar |
| CN117949972A (en) * | 2024-03-26 | 2024-04-30 | 中国人民解放军国防科技大学 | Laser radar-based method, system and device for measuring wind field of take-off and landing channel of airplane |
| CN117949972B (en) * | 2024-03-26 | 2024-06-11 | 中国人民解放军国防科技大学 | Laser radar-based method, system and device for measuring wind field of take-off and landing channel of airplane |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109703770B (en) | 2021-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109703770A (en) | Based on the carrier-borne machine aided of anemometry laser radar and CFD database, method drops | |
| CN103257342B (en) | Three-dimension laser sensor and two-dimension laser sensor combined calibration method | |
| CN116644608B (en) | Real sea area ship motion forecasting method and system based on marine environment data | |
| CN110009037B (en) | Short-term engineering wind speed prediction method and system based on physical information coupling | |
| CN113868971B (en) | Airport area three-dimensional refined wind field reconstruction method based on numerical simulation model and historical wind field characteristics | |
| CN102176003B (en) | Optimization design method for aerial survey parameter of airborne laser radar | |
| CN108061901A (en) | The method that 3D electric power line models are rebuild based on airborne laser radar point cloud data | |
| CN109978275B (en) | Extreme strong wind speed prediction method and system based on mixed CFD and deep learning | |
| CN113777623B (en) | Prediction and alarm method for airplane wake threat area | |
| CN101957317B (en) | Altitude distribution mode measurer of refractive index structural constants of atmospheric turbulence | |
| CN107145647B (en) | Method for correcting deviation of measured data of sea surface wind speed and wind direction of ship | |
| CN112965084B (en) | Airport wind field characteristic detection method, device and equipment based on laser radar | |
| CN107621628A (en) | One kind placement angle error calibration method | |
| CN106706133A (en) | Spot-like target attitude estimation method and system | |
| CN112163381B (en) | Lateral boundary condition setting method suitable for complex terrain wind field flow numerical simulation | |
| CN113687447A (en) | Local area wind field monitoring method based on multiple wind measuring devices | |
| CN104504255A (en) | Method for determining lifting force and resistance moment of spiral wing | |
| CN105652271A (en) | Super-resolution processing method for augmented Lagrangian real-beam radar angle | |
| CN111693999A (en) | Multi-sensor fusion wind speed and direction estimation method based on radar wind measurement combination strategy | |
| CN112446844A (en) | Point cloud feature extraction and registration fusion method | |
| CN111624623A (en) | Wind field inversion method based on laser radar non-uniform scanning | |
| CN108050995B (en) | Oblique photography non-image control point aerial photography measurement area merging method based on DEM | |
| CN113504543A (en) | Unmanned aerial vehicle LiDAR system positioning and attitude determination system and method | |
| CN115585814B (en) | Aircraft variable-scale terrain following method based on settlement mechanism | |
| Gao et al. | A hybrid method for fine-scale wind field retrieval based on machine learning and data assimilation |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20210514 Address after: 250101 Room 501, building 5, 747 Shunhua Road, high tech Zone, Jinan City, Shandong Province Applicant after: Shandong Guoyao quantum Radar Technology Co.,Ltd. Address before: No.99, xiupu Road, Pudong New Area, Shanghai, 201315 Applicant before: GUOYAO QUANTUM RADAR TECHNOLOGY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |