CN108195558B - A kind of ascending air measurement method for weather modification aircraft - Google Patents
A kind of ascending air measurement method for weather modification aircraft Download PDFInfo
- Publication number
- CN108195558B CN108195558B CN201711280576.8A CN201711280576A CN108195558B CN 108195558 B CN108195558 B CN 108195558B CN 201711280576 A CN201711280576 A CN 201711280576A CN 108195558 B CN108195558 B CN 108195558B
- Authority
- CN
- China
- Prior art keywords
- cos
- sin
- cosβ
- ascending air
- vtz
- 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
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention belongs to aeronautical meteorology field, specifically a kind of ascending air measurement method for weather modification aircraft.Solve the problems, such as that conventional method ascending air can not be settled accounts, include the following steps, 101~obtain true air speed Vt, angle of attack and yaw angle β that fixed wing aircraft flies, 102~by vector VtSb=(VtX, VtY, VtZ) earth axes are transformed to by body coordinate system, 103~obtain the three-dimensional ground velocity VgXg of fixed wing aircraft flight, VgYg, VgZg, pitching angle theta, roll angle Φ and course angle Ψ；104~in ground coordinate system, three-dimensional ground velocity subtracts three-dimensional air speed and obtains ascending air.The present invention can accurately measure ascending air, and the measurement of ascending air can be used to judge cloud interior air-flow direction and size in the field of increasing rain artificially, and for the formulation of artificial rain operation scheme, improve catalytic efficiency.
Description
Technical field
The invention belongs to aeronautical meteorology field, specifically a kind of ascending air measurement side for weather modification aircraft
Method.
Background technique
Meteorological field mainly passes through the measurement that release balloon carries out upper-level winds at present, and the balloon that will ascend to heaven is as mobile with air-flow
Particle, with ground installation tracking balloon track of ascending to heaven, read the elevation angle, azimuth, oblique distance of its time interval, determine its sky
Between position coordinate values, the mean wind speed wind direction in the passed through height of balloon can be found out.This method be only capable of measuring horizontal wind speed and
Horizontal wind direction cannot detect ascending air, and a detection time is not substantially less than 2h, and space-time is representative poor.
Some carry out upper air wind finding using fixed wing aircraft, when aircraft keeps horizontal flight posture, the speed of a ship or plane, air speed and
The vector triangle relationship that wind speed is constituted calculates upper-level winds, is also only capable of resolving and obtains horizontal wind, is unable to get ascending air.
Summary of the invention
The present invention provides a kind of for weather modification to solve the problems, such as that conventional method ascending air can not be settled accounts
The ascending air measurement method of aircraft.
The present invention takes following technical scheme: a kind of ascending air measurement method for weather modification aircraft, packet
Include following steps,
101~true air speed Vt, angle of attack and yaw angle β that fixed wing aircraft flies are obtained, according to angle of attack, yaw angle β meter
Three-dimensional component VtX, VtY, the VtZ of true air speed Vt in the body coordinate system are calculated, vector VtSb=(VtX, VtY, VtZ) is denoted as,
Wherein: VtX=Vt*cos β * cos α, VtY=Vt*sin β, VtZ=Vt*cos β * sin α；
102~be denoted as vector VtSb=(VtX, VtY, VtZ) by the vector that body coordinate system transforms to earth axes
VtSg=(VtXg, VtYg, VtZg)；Transition matrix of the body coordinate system Sb to earth axes Sg are as follows:
VtSg=SθψФ× VtSb, wherein
VtYg=VtX*cos θ * sin ψ+VtY* (sin θ * sin ψ * sin φ+cos ψ * cos φ)+VtZ* (sin θ * sin ψ *
cosφ-cosψ*sinφ)
VtZg=VtY*cos θ * sin φ+VtZ*cos θ * cos φ-VtX*sin θ；
VtX, VtY, VtZ are substituted into joint account in above-mentioned three formula to obtain:
VtXg=Vt* (cos α * cos β * cos θ * cos ψ+sin β * sin θ * cos ψ * sin φ-sin β * sin ψ * cos φ+
sinα*cosβ*sinθ*cosψ*cosφ+sinα*cosβ*sinψ*sinφ)
VtYg=Vt* (cos α * cos β * cos θ * sin ψ+sin β * sin θ * sin ψ * sin φ+sin β * cos ψ * cos φ+
sinα*cosβ*sinθ*sinψ*cosφ-sinα*cosβ*cosψ*sinφ)
103~obtain three-dimensional ground velocity VgXg, VgYg that fixed wing aircraft flies, VgZg, pitching angle theta, roll angle Φ and boat
To angle Ψ；
104~in earth axes, three-dimensional ground velocity subtracts three-dimensional air speed and obtains ascending air.By three-dimensional ground velocity and three-dimensional space
Speed vector subtraction in earth axes obtains three-dimensional velocity, relative to earth axes, i.e.,
VwXg=VgXg-VtXg
VwYg=VgYg-VtYg
VwZg=VgZg-VtZg
VwXg=VgXg-Vt (cos α * cos β * cos θ * cos ψ+sin β * sin θ * cos ψ * sin φ-sin β * sin ψ * cos
φ+sinα*cosβ*sinθ*cosψ*cosφ+sinα*cosβ*sinψ*sinφ)
VwYg=VgYg-Vt (cos α * cos β * cos θ * sin ψ+sin β * sin θ * sin ψ * sin φ+sin β * cos ψ * cos
φ+sinα*cosβ*sinθ*sinψ*cosφ-sinα*cosβ*cosψ*sinφ)
VwZg=VgZg-Vt (sin β * cos θ * sin φ+sin α * cos β * cos θ * cos φ-cos α * cos β * sin θ)
Wherein: VwXg is the east orientation wind that north orientation wind speed, the VwYg relative to earth axes are relative to earth axes
Speed, VwZg are the ascending air relative to earth axes.
The measurement of ascending air is the important indicator increased rain artificially, and one side Yun Zhongyou ascending air proves that cloud is in development
Stage, catalysis potentiality with higher prove that cloud is in the mature and extinction stage if Yun Zhongwei down current, catalysis potentiality compared with
It is low.On the other hand there is the place of ascending air that can sow catalyst in cloud base, brought into entire cloud by ascending air, so that
Catalyst distribution is more uniform.
Compared with prior art, the present invention can calculate rising by the atmospheric parameter and flight parameter of fixed wing aircraft
Air-flow, the deficiency of ascending air cannot accurately be measured at present by solving meteorological field.The present invention can accurately measure ascending air, rise
The measurement of air-flow can be used to judge cloud interior air-flow direction and size in the field of increasing rain artificially, for the formulation of artificial rain operation scheme,
Improve catalytic efficiency.
Detailed description of the invention
Fig. 1 is flow chart of the present invention.
Specific embodiment
To keep the technical problem to be solved in the present invention, technical solution and advantage clearer, below in conjunction with attached drawing and tool
Body embodiment is described in detail.
A kind of ascending air measurement method for weather modification aircraft, includes the following steps,
101~true air speed Vt, angle of attack and yaw angle β that fixed wing aircraft flies are obtained, according to angle of attack, yaw angle β meter
Three-dimensional component VtX, VtY, the VtZ of true air speed Vt in the body coordinate system are calculated, vector VtSb=(VtX, VtY, VtZ) is denoted as,
Wherein: VtX=Vt*cos β * cos α, VtY=Vt*sin β, VtZ=Vt*cos β * sin α；
102~be denoted as vector VtSb=(VtX, VtY, VtZ) by the vector that body coordinate system transforms to earth axes
VtSg=(VtXg, VtYg, VtZg)；Transition matrix of the body coordinate system Sb to earth axes Sg are as follows:
VtSg=sθψФ× VtSb, wherein
VtYg=VtX*cos θ * sin ψ+VtY* (sin θ * sin ψ * sin φ+cos ψ * cos φ)+VtZ* (sin θ * sin ψ *
cosφ-cosψ*sinφ)
VtZg=VtY*cos θ * sin φ+VtZ*cos θ * cos φ-VtX*sin θ；
VtX, VtY, VtZ are substituted into joint account in above-mentioned three formula to obtain:
VtXg=Vt* (cos α * cos β * cos θ * cos ψ+sin β * sin θ * cos ψ * sin φ-sin β * sin ψ * cos φ+
sinα*cosβ*sinθ*cosψ*cosφ+sinα*cosβ*sinψ*sinφ)
VtYg=Vt* (cos α * cos β * cos θ * sin ψ+sin β * sin θ * sin ψ * sin φ+sin β * cos ψ * cos φ+
sinα*cosβ*sinθ*sinψ*cosφ-sinα*cosβ*cosψ*sinφ)
103~obtain three-dimensional ground velocity VgXg, VgYg that fixed wing aircraft flies, VgZg, pitching angle theta, roll angle Φ and boat
To angle Ψ；
104~in earth axes, three-dimensional ground velocity subtracts three-dimensional air speed and obtains ascending air.By three-dimensional ground velocity and three-dimensional space
Speed vector subtraction in ground coordinate system obtains three-dimensional velocity, relative to earth axes, i.e.,
VwXg=VgXg-VtXg
VwYg=VgYg-VtYg
VwZg=VgZg-VtZg
VwXg=VgXg-Vt (cos α * cos β * cos θ * cos ψ+sin β * sin θ * cos ψ * sin φ-sin β * sin ψ * cos
φ+sinα*cosβ*sinθ*cosψ*cosφ+sinα*cosβ*sinψ*sinφ)
VwYg=VgYg-Vt (cos α * cos β * cos θ * sin ψ+sin β * sin θ * sin ψ * sin φ+sin β * cos ψ * cos
φ+sinα*cosβ*sinθ*sinψ*cosφ-sinα*cosβ*cosψ*sinφ)
VwZg=VgZg-Vt (sin β * cos θ * sin φ+sin α * cos β * cos θ * cos φ-cos α * cos β * sin θ)
Wherein: VwXg is the east orientation wind that north orientation wind speed, the VwYg relative to earth axes are relative to earth axes
Speed, VwZg are the ascending air relative to earth axes.
The above is a preferred embodiment of the present invention, it is noted that for those skilled in the art
For, without departing from the principles of the present invention, several improvements and modifications can also be made, these improvements and modifications
It should be regarded as protection scope of the present invention.
Claims (3)
1. a kind of ascending air measurement method for weather modification aircraft, it is characterised in that:
Include the following steps,
101~true air speed Vt, angle of attack and yaw angle β that fixed wing aircraft flies are obtained, according to angle of attack, yaw angle β is calculated
Three-dimensional component VtX, VtY, the VtZ of true air speed Vt in the body coordinate system, are denoted as vector VtSb=(VtX, VtY, VtZ), in which:
VtX=Vt*cos β * cos α, VtY=Vt*sin β, VtZ=Vt*cos β * sin α；
102~vector VtSb=(VtX, VtY, VtZ) is denoted as VtSg by the vector that body coordinate system transforms to earth axes
=(VtXg, VtYg, VtZg)；
103~obtain three-dimensional ground velocity VgXg, VgYg that fixed wing aircraft flies, VgZg, pitching angle theta, roll angle Φ and course angle
Ψ；
104~in earth axes, three-dimensional ground velocity subtracts three-dimensional air speed and obtains ascending air.
2. the ascending air measurement method according to claim 1 for weather modification aircraft, it is characterised in that: institute
It states in step 102, transition matrix of the body coordinate system Sb to earth axes Sg are as follows:
VtSg=SθψФ× VtSb, wherein
VtYg=VtX*cos θ * sin ψ+VtY* (sin θ * sin ψ * sin φ+cos ψ * cos φ)+VtZ* (sin θ * sin ψ * cos
φ-cosψ*sinφ)
VtZg=VtY*cos θ * sin φ+VtZ*cos θ * cos φ-VtX*sin θ；
VtX, VtY, VtZ are substituted into joint account in above-mentioned three formula to obtain:
VtXg=Vt* (cos α * cos β * cos θ * cos ψ+sin β * sin θ * cos ψ * sin φ-sin β * sin ψ * cos φ+sin α *
cosβ*sinθ*cosψ*cosφ+sinα*cosβ*sinψ*sinφ)
VtYg=Vt* (cos α * cos β * cos θ * sin ψ+sin β * sin θ * sin ψ * sin φ+sin β * cos ψ * cos φ+sin α *
cosβ*sinθ*sinψ*cosφ-sinα*cosβ*cosψ*sinφ)
3. the ascending air measurement method according to claim 2 for weather modification aircraft, it is characterised in that: institute
It states in step 104, three-dimensional ground velocity and three-dimensional air speed vector subtraction in earth axes is obtained into three-dimensional velocity, relative to ground
Coordinate system, i.e.,
VwXg=VgXg-VtXg
VwYg=VgYg-VtYg
VwZg=VgZg-VtZg
VwXg=VgXg-Vt (cos α * cos β * cos θ * cos ψ+sin β * sin θ * cos ψ * sin φ-sin β * sin ψ * cos φ+
sinα*cosβ*sinθ*cosψ*cosφ+sinα*cosβ*sinψ*sinφ)
VwYg=VgYg-Vt (cos α * cos β * cos θ * sin ψ+sin β * sin θ * sin ψ * sin φ+sin β * cos ψ * cos φ+
sinα*cosβ*sinθ*sinψ*cosφ-sinα*cosβ*cosψ*sinφ)
VwZg=VgZg-Vt (sin β * cos θ * sin φ+sin α * cos β * cos θ * cos φ-cos α * cos β * sin θ)
Wherein: VwXg is the east orientation wind speed that north orientation wind speed, the VwYg relative to earth axes are relative to earth axes,
VwZg is the ascending air relative to earth axes.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711280576.8A CN108195558B (en) | 2017-12-06 | 2017-12-06 | A kind of ascending air measurement method for weather modification aircraft |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711280576.8A CN108195558B (en) | 2017-12-06 | 2017-12-06 | A kind of ascending air measurement method for weather modification aircraft |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108195558A CN108195558A (en) | 2018-06-22 |
| CN108195558B true CN108195558B (en) | 2019-09-13 |
Family
ID=62573913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711280576.8A Active CN108195558B (en) | 2017-12-06 | 2017-12-06 | A kind of ascending air measurement method for weather modification aircraft |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108195558B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109710961B (en) * | 2018-10-26 | 2023-01-13 | 中国飞行试验研究院 | High-altitude unmanned aerial vehicle limit rising data processing method based on GPS data |
| CN110414110B (en) * | 2019-07-19 | 2023-01-06 | 中仿智能科技（上海）股份有限公司 | Airplane stress simulation method used in flight stall state |
| CN111222082B (en) * | 2020-01-15 | 2022-01-04 | 南京航空航天大学 | Method for identifying position and speed of ascending airflow suitable for unmanned aerial vehicle |
| CN115438293B (en) * | 2022-09-09 | 2023-06-06 | 青海省人工影响天气办公室 | Ground smoke furnace artificial precipitation operation condition identification method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2048508A2 (en) * | 2007-10-12 | 2009-04-15 | Interactive Flow Studies LLC | Fluid flow computation, visualization, and analysis |
| CN104459193A (en) * | 2014-12-05 | 2015-03-25 | 中国航天空气动力技术研究院 | Crosswind information estimation method based on unmanned aerial vehicle crabbing method |
| CN105075759A (en) * | 2015-09-22 | 2015-11-25 | 许琬晨 | Cloud water reduction and rain enhancement scale prediction method |
| CN105718727A (en) * | 2016-01-19 | 2016-06-29 | 中国科学院光电研究院 | Stratospheric airship flight performance parameter estimation method and system |
-
2017
- 2017-12-06 CN CN201711280576.8A patent/CN108195558B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2048508A2 (en) * | 2007-10-12 | 2009-04-15 | Interactive Flow Studies LLC | Fluid flow computation, visualization, and analysis |
| CN104459193A (en) * | 2014-12-05 | 2015-03-25 | 中国航天空气动力技术研究院 | Crosswind information estimation method based on unmanned aerial vehicle crabbing method |
| CN105075759A (en) * | 2015-09-22 | 2015-11-25 | 许琬晨 | Cloud water reduction and rain enhancement scale prediction method |
| CN105718727A (en) * | 2016-01-19 | 2016-06-29 | 中国科学院光电研究院 | Stratospheric airship flight performance parameter estimation method and system |
Non-Patent Citations (1)
| Title |
|---|
| 层状暖云降雨及其催化的数值模拟;胡志晋 等;《气象学报》;19830228;第41卷(第1期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108195558A (en) | 2018-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108195558B (en) | A kind of ascending air measurement method for weather modification aircraft | |
| CN205139117U (en) | Perpendicular monitoring system of atmosphere pollution based on on -vehicle captive balloon | |
| Landwehr et al. | Direct flux measurements from mobile platforms at sea: Motion and airflow distortion corrections revisited | |
| CN110647171B (en) | Method for plant protection unmanned aerial vehicle to adjust flight line independently | |
| Wildmann et al. | Measuring the local wind field at an escarpment using small remotely-piloted aircraft | |
| CN109059916A (en) | A kind of aerostatics occultation forecasting procedure based on inertial navigation | |
| CN109579841B (en) | High-precision positioning method for vehicle-mounted fire-fighting high-load rotor unmanned aerial vehicle under GPS (Global positioning System) rejection condition | |
| US20150260879A1 (en) | Tethersonde system and observation method thereby | |
| CN110441767A (en) | The measurement method and system of power transmission line corridor screen of trees head room | |
| CN1036185A (en) | Wind shear detection system | |
| CN111239788A (en) | Unmanned aerial vehicle autonomous positioning system and method for bridge detection | |
| CN107064545A (en) | A kind of detection method of less than 200 meters shallow-layer wind | |
| Schelbergen et al. | Validation of the quasi-steady performance model for pumping airborne wind energy systems | |
| CN208559732U (en) | A kind of high-altitude mapping balloon robot | |
| CN110011032B (en) | Retraction control device and method for airborne emergency communication system antenna | |
| CN113325880B (en) | Sun and moon tracking device and method suitable for motion state | |
| CN202631567U (en) | GPS (Global Positioning System) based mobile anemometry device | |
| CN111222082B (en) | Method for identifying position and speed of ascending airflow suitable for unmanned aerial vehicle | |
| JP3571268B2 (en) | Fog observation radar device | |
| Kocer et al. | Wake structure of a 2mw wind turbine measured using an instrumented uav | |
| CN104698510A (en) | Atmosphere thermal inversion layer detection method based on captive balloon sonde | |
| CN102407934A (en) | Hitching type movable radar | |
| CN219295700U (en) | Survey and drawing unmanned aerial vehicle that wind-resistant performance is strong | |
| CN211206243U (en) | Unmanned aerial vehicle with accurate height measurement device | |
| CN109178300A (en) | A kind of wind measuring device based on multi-rotor unmanned aerial vehicle platform |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |