CN106096091A - A kind of airplane motion analogy method - Google Patents
A kind of airplane motion analogy method Download PDFInfo
- Publication number
- CN106096091A CN106096091A CN201610374361.1A CN201610374361A CN106096091A CN 106096091 A CN106096091 A CN 106096091A CN 201610374361 A CN201610374361 A CN 201610374361A CN 106096091 A CN106096091 A CN 106096091A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- axis system
- earth
- cos
- sin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Navigation (AREA)
Abstract
The invention discloses a kind of airplane motion analogy method, comprise the steps of S1, the original state parameter of input aircraft;S2, inputs aircraft current state parameter;S3, calculates the derivative of earth's axis system kinematic parameter;S4, using the original state parameter in described step S1 as initial value for integral;Using the derivative in described step S3 as integration rate, it is integrated, solves current four element values of aircraft, fly axis system current angular velocity, aircraft current location, aircraft earth's axis system present speed, aircraft relative to the current angle of attack of earth's axis system speed, aircraft relative to the current yaw angle of earth's axis system speed;S5, calculates aircraft true empty-running speed, the actual angle of attack and actual yaw angle;S6, delivers to aerodynamics evaluation module, the aerodynamic force suffered by accurate simulated aircraft and aerodynamic moment by described aircraft true empty-running speed, the actual angle of attack, actual yaw angle.It is an advantage of the current invention that: to realize the atmospheric perturbation fine analog to aircraft flight Parameters variation process, make airplane motion response the most true to nature.
Description
Technical field
The present invention relates to airplane design technical field, be specifically related to a kind of airplane motion analogy method.
Background technology
In field of airplane design and flight simulation field, it usually needs the assessment atmospheric perturbation impact on airplane motion, especially
It is the atmospheric perturbation that turbulent flow and wind shear etc. are frequently encountered by practical flight, carries out atmospheric perturbation in flight simulation environment
Journey simulation is requisite function in Project R&D and pilot training.Such as, aircraft enters experience wind shear of nearly stage, needs
Fine analog aircraft angle of attack and the progressive formation of yaw angle, so that airplane aerodynamic and the truest aircraft of motion are suitable.At present, greatly
Gas disturbance is generally directly injected in airplane motion equation with the form of terrestrial coordinate system three axle wind speed component, directly changes and flies
Machine ground velocity increment, have ignored ground velocity continually varying dynamic process in Live Flying, reduces people's fidelity at loop test,
Have impact on control law and the credibility of flight quality assessment result.
Summary of the invention
It is an object of the invention to provide a kind of airplane motion analogy method, to solve or at least to alleviate in background technology to be deposited
The problem at least one place.
The technical solution used in the present invention is: provides a kind of airplane motion analogy method, comprises the steps of
S1, the original state parameter of input aircraft, described original state parameter comprises the initial Eulerian angles of aircraft, aircraft body
Axle system initial angular velocity, aircraft initial position, aircraft earth's axis system initial velocity, aircraft are relative to the initial angle of attack of earth's axis system speed
With aircraft relative to the initial side-slip angle of earth's axis system speed;
S2, inputs aircraft current state parameter, and described aircraft current state parameter comprises the body axle of bonding force suffered by aircraft
It it is the body axle system three axle component of three axle components and bonding force square;Aircraft gross mass, the moment of inertia of each axle of axis system that is diversion and used
Property long-pending;The wind speed of current aircraft local environment three axle components in earth's axis system;
S3, calculates the derivative of earth's axis system kinematic parameter, calculates current four element values of aircraft correction,
q1=q1/norm
q2=q2/norm
q3=q3/norm
q4=q4/norm
In formula, norm is the modulus value of current four elements, q1、q2、q3And q4For revised four element currencys;
Current direction cosine matrix is calculated according to revised four element currencys;Calculate according to described direction cosine matrix
The aircraft earth's axis system present speed three axle components in earth's axis system, described component is also the derivative of aircraft current location;Calculating flies
The derivative of machine earth's axis system present speed;Aircraft is relative to the derivative of the earth's axis system current angle of attack of speed;Aircraft is relative to earth's axis system speed
Spend the derivative of current yaw angle;Calculating the derivative flying axis system current angular velocity, described current angular velocity comprises current rolling
Angular velocity, rate of pitch and yaw rate;Calculate the derivative of four elements;
Initial for aircraft in described step S1 Eulerian angles are converted to four element initial values by S4, and in addition to initial Eulerian angles
Other original state parameters are together as initial value for integral;By the derivative of four elements in described step S3, fly axis system roll angle
The derivative of speed, the derivative of rate of pitch, the derivative of yaw rate, the derivative of aircraft current location, aircraft earth's axis system work as
The derivative of front speed, aircraft currently break away relative to earth's axis system speed relative to derivative and the aircraft of the earth's axis system current angle of attack of speed
The derivative at angle, as integration rate, is integrated, and solves current four element values of aircraft, flies axis system current angular velocity, aircraft
Current location, aircraft earth's axis system present speed, aircraft relative to the current angle of attack of earth's axis system speed, aircraft relative to earth's axis system speed
The current yaw angle of degree;
S5, calculates aircraft airspeed at three axle components of earth's axis system, according to described aircraft airspeed at three axle components of earth's axis system
Calculate the aircraft airspeed three axle components in body axle system, calculate aircrafts according to described aircraft airspeed at three axle components of body axle system actual
Air speed, the actual angle of attack and actual yaw angle;
S6, delivers to aerodynamics evaluation module by described aircraft true empty-running speed, the actual angle of attack, actual yaw angle, accurately simulates
Aerodynamic force suffered by aircraft and aerodynamic moment.
Preferably, in described step S4, initial for aircraft Eulerian angles being converted to four element initial values, concrete conversion method is,
q10=cos (θ0/2)cos(φ0/2)cos(ψ0/2)+sin(θ0/2)sin(φ0/2)sin(ψ0/2)
q20=cos (θ0/2)sin(φ0/2)cos(ψ0/2)-sin(θ0/2)cos(φ0/2)sin(ψ0/2)
q30=sin (θ0/2)cos(φ0/2)cos(ψ0/2)+cos(θ0/2)sin(φ0/2)sin(ψ0/2)
q40=cos (θ0/2)cos(φ0/2)sin(ψ0/2)-sin(θ0/2)sin(φ0/2)cos(ψ0/2)
In formula, q10、q20、q30And q40It is four element initial values, φ0For roll angle, θ0For the angle of pitch, ψ0For course angle.
Preferably, in described step S3,
The specific algorithm of current direction cosine matrix is,
Ay=2 (q2q3-q1q4)
Az=2 (q2q4+q1q3)
Bx=2 (q2q3+q1q4)
Bz=2 (q3q4-q1q2)
Dx=2 (q2q4-q1q3)
Dy=2 (q3q4+q1q2)
In formula, q1、q2、q3And q4For revised four element currencys;
The aircraft earth's axis system present speed three axle component specific algorithms in earth's axis system are,
In formula, Vxg、VygAnd VzgFor the aircraft earth's axis system speed three axle components in earth's axis system, V is that aircraft earth's axis system is current
Speed, α is the aircraft current angle of attack relative to earth's axis system speed, and β is the aircraft current yaw angle relative to earth's axis system speed;
The derivative of aircraft earth's axis system present speed specifically,
Aircraft relative to the earth's axis system current angle of attack of speed derivative specifically,
Aircraft relative to the earth's axis system current yaw angle of speed derivative specifically,
P for flying the current angular velocity in roll of axis system, q is rate of pitch, r is yaw rate, and g is aircraft present bit
The gravity acceleration value put, Fxt、Fyt、FztBody axle system three axle component for aircraft currently suffered bonding force;
Fly the derivative of axis system current angular velocity specifically,
B1=L+ (Iy-Iz)qr+Izxpq
B2=M+ (Iz-Ix)rp-Izx(p2-r2)
B3=N+ (Ix-Iy)pq-Izxqr
For fly the current angular velocity in roll of axis system derivative,For the derivative of rate of pitch,For yaw rate
Derivative, B1、B2、B3For intermediate variable, the body axle system three axle component of L, M, N respectively aircraft currently suffered bonding force square, Ix、
Iy、Iz、IzxThe moment of inertia of axis system x-axis, y-axis and z-axis of being respectively diversion, IzxInertia for be diversion axis system z-axis and y-axis
Long-pending;
The derivative of four elements is,
Preferably, in described step S5, aircraft airspeed at three axle component specific algorithms of earth's axis system is,
In formula, Vx、Vy、VzIt is respectively the aircraft airspeed three axle components in earth's axis system, Vxg、Vyg、VzgIt is respectively aircraft ground
The speed three axle components in earth's axis system, Wxg、Wyg、WzgThe respectively wind speed of aircraft local environment three axle components in earth's axis system;
Aircraft airspeed at three axle components of body axle system is,
In formula, Vxt、Vyt、VztIt is respectively the aircraft airspeed three axle components in body axle system;
The specific algorithm of aircraft airspeed, the elevation angle and yaw angle is,
αV=arctan (Vzt/Vxt)
βV=arcsin (Vyt/VV)
In formula, VVFor aircraft airspeed, αVFor aircraft angle of attack, βVFor aircraft yaw angle.
The beneficial effects of the present invention is:
The method of the present invention can realize atmospheric perturbation to flight parameters such as aircraft angle of attack, yaw angle, true air speed and ground velocity
The fine analog of change procedure, makes airplane motion response the most true to nature.
Use this algorithm improvement full digital trigger technique platform, engineering simulator, comprehensive " iron bird " testing stand and flight training mould
Intend the airplane motion algorithm in the flight simulation such as device, it is possible to increase Design of Flight Control, test and the accuracy trained and
Effectiveness, increases people's credibility in loop test assessment result.
Accompanying drawing explanation
Fig. 1 is the flow chart of the airplane motion analogy method of one embodiment of the 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, 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 creative 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 effect of this innovatory algorithm for convenience of description, establishes a set of airplane motion simulated program.Emulation was set at 0 second
Starting, introduce the earth's axis system north orientation step wind when 0.2 second, when 0.4 second, step terminates, and introduces the earth's axis system east orientation step wind when 0.6 second,
When 0.8 second, step terminates, and when 1 second, emulation terminates.Wherein, step monsoon intensity is 5m/s.
As it is shown in figure 1, the airplane motion analogy method in the present embodiment, comprise the steps of
S1, the original state parameter of input aircraft, described original state parameter comprises the initial Eulerian angles of aircraft, aircraft body
Axle system initial angular velocity, aircraft initial position, aircraft earth's axis system initial velocity, aircraft are relative to the initial angle of attack of earth's axis system speed
With aircraft relative to the initial side-slip angle of earth's axis system speed;The initial Eulerian angles of described aircraft comprise aircraft pitch angle, aircraft rolling
Angle and vector angle;The described axis system initial angular velocity that flies comprises and flies the initial angular velocity in roll of axis system, flies axis
It is initial pitch angle speed, flies the initial yaw rate of axis system;Described aircraft initial position comprises aircraft initial north orientation position
Put, aircraft initial east orientation position, the initial sky of aircraft is to position.
In the present embodiment,
Aircraft pitch angle θ0=2deg;
Aircraft roll angle φ0=0deg;
Vector angle ψ0=45deg;
Fly axis system initial angular velocity in roll p0=0deg/s;
Fly axis system initial pitch angle speed q0=0deg/s;
Fly initial yaw rate r of axis system0=0deg/s;
Aircraft initial north orientation position x0=0m;
Aircraft initial east orientation position y0=0m;
The initial sky of aircraft is to position z0=500m;
Aircraft earth's axis system initial velocity V0=120m/s;
Aircraft is relative to the initial angle of attack α of earth's axis system speed0=2deg;
Aircraft is relative to the initial side-slip angle beta of earth's axis system speed0=0deg.
S2, inputs aircraft current state parameter, and described aircraft current state parameter comprises the body axle of bonding force suffered by aircraft
It it is the body axle system three axle component of three axle components and bonding force square;Aircraft gross mass, the moment of inertia of each axle of axis system that is diversion and used
Property long-pending;The wind speed of current aircraft local environment three axle components in earth's axis system.
In the present embodiment, input aircraft current state, be all set to fixed value
The body axle system three axle component of bonding force suffered by aircraft (without gravity) is respectively
Fxt=34.224191441307575;Fyt=0;Fzt=0.024677670778336;
The body axle system three axle component of bonding force square suffered by aircraft
L=0;M=0;N=0;
Aircraft Quality performance data
MassT=100kg;Ix=100000;Iy=100000;Iy=100000;Izx=-10;
The wind speed of aircraft local environment three axle components in earth's axis system
Earth's axis system north orientation step wind 5m/s when 0.2 second, when 0.4 second, step terminates;
Earth's axis system east orientation step wind 5m/s when 0.6 second, when 0.8 second, step terminates.
S3, calculates the derivative of earth's axis system kinematic parameter, calculates current four element values of aircraft correction,
q1=q1/norm
q2=q2/norm
q3=q3/norm
q4=q4/norm
In formula, norm is the modulus value of current four elements, q1、q2、q3And q4For revised four element currencys;
Current direction cosine matrix is calculated according to revised four element currencys;Calculate according to described direction cosine matrix
The aircraft earth's axis system present speed three axle components in earth's axis system, described component is also the derivative of aircraft current location;Calculating flies
The derivative of machine earth's axis system present speed;Aircraft is relative to the derivative of the earth's axis system current angle of attack of speed;Aircraft is relative to earth's axis system speed
Spend the derivative of current yaw angle;Calculating the derivative flying axis system current angular velocity, described current angular velocity comprises current rolling
Angular velocity, rate of pitch and yaw rate;Calculate the derivative of four elements;
Initial for aircraft in described step S1 Eulerian angles are converted to four element initial values by S4, and in addition to initial Eulerian angles
Other original state parameters are together as initial value for integral;By the derivative of four elements in described step S3, fly axis system roll angle
The derivative of speed, the derivative of rate of pitch, the derivative of yaw rate, the derivative of aircraft current location, aircraft earth's axis system work as
The derivative of front speed, aircraft currently break away relative to earth's axis system speed relative to derivative and the aircraft of the earth's axis system current angle of attack of speed
The derivative at angle, as integration rate, is integrated, and solves current four element values of aircraft, flies axis system current angular velocity, aircraft
Current location, aircraft earth's axis system present speed, aircraft relative to the current angle of attack of earth's axis system speed, aircraft relative to earth's axis system speed
The current yaw angle of degree;
S5, calculates aircraft airspeed at three axle components of earth's axis system, according to described aircraft airspeed at three axle components of earth's axis system
Calculate the aircraft airspeed three axle components in body axle system, calculate aircrafts according to described aircraft airspeed at three axle components of body axle system actual
Air speed, the actual angle of attack and actual yaw angle;
S6, delivers to aerodynamics evaluation module by described aircraft true empty-running speed, the actual angle of attack, actual yaw angle, accurately simulates
Aerodynamic force suffered by aircraft and aerodynamic moment.
In the present embodiment, initial for aircraft Eulerian angles are converted to four element initial value, the specifically sides of conversion by described step S4
Method is,
q10=cos (θ0/2)cos(φ0/2)cos(ψ0/2)+sin(θ0/2)sin(φ0/2)sin(ψ0/2)
q20=cos (θ0/2)sin(φ0/2)cos(ψ0/2)-sin(θ0/2)cos(φ0/2)sin(ψ0/2)
q30=sin (θ0/2)cos(φ0/2)cos(ψ0/2)+cos(θ0/2)sin(φ0/2)sin(ψ0/2)
q40=cos (θ0/2)cos(φ0/2)sin(ψ0/2)-sin(θ0/2)sin(φ0/2)cos(ψ0/2)
In formula, q10、q20、q30And q40It is four element initial values, φ0For roll angle, θ0For the angle of pitch, ψ0For course angle.
In the present embodiment, in described step S3,
The specific algorithm of current direction cosine matrix is,
Ay=2 (q2q3-q1q4)
Az=2 (q2q4+q1q3)
Bx=2 (q2q3+q1q4)
Bz=2 (q3q4-q1q2)
Dx=2 (q2q4-q1q3)
Dy=2 (q3q4+q1q2)
In formula, q1、q2、q3And q4For revised four element currencys;
The aircraft earth's axis system present speed three axle component specific algorithms in earth's axis system are,
In formula, Vxg、VygAnd VzgFor the aircraft earth's axis system speed three axle components in earth's axis system, V is that aircraft earth's axis system is current
Speed, α is the aircraft current angle of attack relative to earth's axis system speed, and β is the aircraft current yaw angle relative to earth's axis system speed;
The derivative of aircraft earth's axis system present speed specifically,
Aircraft relative to the earth's axis system current angle of attack of speed derivative specifically,
Aircraft relative to the earth's axis system current yaw angle of speed derivative specifically,
P for flying the current angular velocity in roll of axis system, q is rate of pitch, r is yaw rate, and g is aircraft present bit
The gravity acceleration value put, Fxt、Fyt、FztBody axle system three axle component for aircraft currently suffered bonding force;
Fly the derivative of axis system current angular velocity specifically,
B1=L+ (Iy-Iz)qr+Izxpq
B2=M+ (Iz-Ix)rp-Izx(p2-r2)
B3=N+ (Ix-Iy)pq-Izxqr
For fly the current angular velocity in roll of axis system derivative,For the derivative of rate of pitch,For yaw rate
Derivative, B1、B2、B3For intermediate variable, the body axle system three axle component of L, M, N respectively aircraft currently suffered bonding force square, Ix、
Iy、Iz、IzxThe moment of inertia of axis system x-axis, y-axis and z-axis of being respectively diversion, IzxInertia for be diversion axis system z-axis and y-axis
Long-pending;
The derivative of four elements is,
In the present embodiment, in described step S5, aircraft airspeed at three axle component specific algorithms of earth's axis system is,
In formula, Vx、Vy、VzIt is respectively the aircraft airspeed three axle components in earth's axis system, Vxg、Vyg、VzgIt is respectively aircraft ground
The speed three axle components in earth's axis system, Wxg、Wyg、WzgThe respectively wind speed of aircraft local environment three axle components in earth's axis system;
Aircraft airspeed at three axle components of body axle system is,
In formula, Vxt、Vyt、VztIt is respectively the aircraft airspeed three axle components in body axle system;
The specific algorithm of aircraft airspeed, the elevation angle and yaw angle is,
αV=arctan (Vzt/Vxt)
βV=arcsin (Vyt/VV)
In formula, VVFor aircraft airspeed, αVFor aircraft angle of attack, βVFor aircraft yaw angle.
Last it is noted that above example is only in order to illustrate technical scheme, it is not intended to limit.To the greatest extent
The present invention has been described in detail by pipe with reference to previous embodiment, it will be understood by those within the art that: it is still
Technical scheme described in foregoing embodiments can be modified, or wherein portion of techniques feature is carried out equivalent replace
Change;And these amendments or replacement, do not make the essence of appropriate technical solution depart from the essence of various embodiments of the present invention technical scheme
God and scope.
Claims (4)
1. an airplane motion analogy method, it is characterised in that comprise the steps of
S1, the original state parameter of input aircraft, described original state parameter comprises the initial Eulerian angles of aircraft, flies axis system
Initial angular velocity, aircraft initial position, aircraft earth's axis system initial velocity, aircraft are relative to the initial angle of attack of earth's axis system speed and fly
Machine is relative to the initial side-slip angle of earth's axis system speed;
S2, inputs aircraft current state parameter, and described aircraft current state parameter comprises the body axle system three of bonding force suffered by aircraft
The body axle system three axle component of axle component and bonding force square;Aircraft gross mass, the moment of inertia of each axle of axis system that is diversion and the product of inertia;
The wind speed of current aircraft local environment three axle components in earth's axis system;
S3, calculates the derivative of earth's axis system kinematic parameter, calculates current four element values of aircraft correction,
q1=q1/norm
q2=q2/norm
q3=q3/norm
q4=q4/norm
In formula, norm is the modulus value of current four elements, q1、q2、q3And q4For revised four element currencys;
Current direction cosine matrix is calculated according to revised four element currencys;Aircraft is calculated according to described direction cosine matrix
The earth's axis system present speed three axle components in earth's axis system, described component is also the derivative of aircraft current location;Calculate aircraft ground
The derivative of axle system present speed;Aircraft is relative to the derivative of the earth's axis system current angle of attack of speed;Aircraft is worked as relative to earth's axis system speed
The derivative of front yaw angle;Calculating the derivative flying axis system current angular velocity, described current angular velocity comprises current roll angle speed
Degree, rate of pitch and yaw rate;Calculate the derivative of four elements;
Initial for aircraft in described step S1 Eulerian angles are converted to four element initial values by S4, with other in addition to initial Eulerian angles
Original state parameter is together as initial value for integral;By the derivative of four elements in described step S3, fly axis system angular velocity in roll
Derivative, the derivative of rate of pitch, the derivative of yaw rate, the derivative of aircraft current location, aircraft earth's axis system the fastest
The derivative of degree, aircraft relative to the derivative of the earth's axis system current angle of attack of speed and aircraft relative to the earth's axis system current yaw angle of speed
Derivative, as integration rate, is integrated, and solves current four element values of aircraft, to fly axis system current angular velocity, aircraft current
Position, aircraft earth's axis system present speed, aircraft relative to the current angle of attack of earth's axis system speed, aircraft relative to earth's axis system speed
Current yaw angle;
S5, calculates the aircraft airspeed three axle components in earth's axis system, calculates at three axle components of earth's axis system according to described aircraft airspeed
Aircraft airspeed at three axle components of body axle system, calculates aircraft reality according to described aircraft airspeed at three axle components of body axle system empty
The angle of attack fast, actual and actual yaw angle;
S6, delivers to aerodynamics evaluation module, accurate simulated aircraft by described aircraft true empty-running speed, the actual angle of attack, actual yaw angle
Suffered aerodynamic force and aerodynamic moment.
2. airplane motion analogy method as claimed in claim 1, it is characterised in that: by initial for aircraft Euler in described step S4
Angle is converted to four element initial values, and concrete conversion method is,
q10=cos (θ0/2)cos(φ0/2)cos(ψ0/2)+sin(θ0/2)sin(φ0/2)sin(ψ0/2)
q20=cos (θ0/2)sin(φ0/2)cos(ψ0/2)-sin(θ0/2)cos(φ0/2)sin(ψ0/2)
q30=sin (θ0/2)cos(φ0/2)cos(ψ0/2)+cos(θ0/2)sin(φ0/2)sin(ψ0/2)
q40=cos (θ0/2)cos(φ0/2)sin(ψ0/2)-sin(θ0/2)sin(φ0/2)cos(ψ0/2)
In formula, q10、q20、q30And q40It is four element initial values, φ0For roll angle, θ0For the angle of pitch, ψ0For course angle.
3. airplane motion analogy method as claimed in claim 1, it is characterised in that: in described step S3,
The specific algorithm of current direction cosine matrix is,
Ay=2 (q2q3-q1q4)
Az=2 (q2q4+q1q3)
Bx=2 (q2q3+q1q4)
Bz=2 (q3q4-q1q2)
Dx=2 (q2q4-q1q3)
Dy=2 (q3q4+q1q2)
In formula, q1、q2、q3And q4For revised four element currencys;
The aircraft earth's axis system present speed three axle component specific algorithms in earth's axis system are,
In formula, Vxg、VygAnd VzgFor the aircraft earth's axis system speed three axle components in earth's axis system, V is that aircraft earth's axis system is the fastest
Degree, α is the aircraft current angle of attack relative to earth's axis system speed, and β is the aircraft current yaw angle relative to earth's axis system speed;
The derivative of aircraft earth's axis system present speed specifically,
Aircraft relative to the earth's axis system current angle of attack of speed derivative specifically,
Aircraft relative to the earth's axis system current yaw angle of speed derivative specifically,
P for flying the current angular velocity in roll of axis system, q is rate of pitch, r is yaw rate, and g is aircraft current location
Gravity acceleration value, Fxt、Fyt、FztBody axle system three axle component for aircraft currently suffered bonding force;
Fly the derivative of axis system current angular velocity specifically,
B1=L+ (Iy-Iz)qr+Izxpq
B2=M+ (Iz-Ix)rp-Izx(p2-r2)
B3=N+ (Ix-Iy)pq-Izxqr
For fly the current angular velocity in roll of axis system derivative,For the derivative of rate of pitch,Leading for yaw rate
Number, B1、B2、B3For intermediate variable, the body axle system three axle component of L, M, N respectively aircraft currently suffered bonding force square, Ix、Iy、Iz、
IzxThe moment of inertia of axis system x-axis, y-axis and z-axis of being respectively diversion, IzxThe product of inertia for be diversion axis system z-axis and y-axis;
The derivative of four elements is,
4. airplane motion analogy method as claimed in claim 1, it is characterised in that: in described step S5, aircraft airspeed is in the earth's axis
Three axle component specific algorithms of system are,
In formula, Vx、Vy、VzIt is respectively the aircraft airspeed three axle components in earth's axis system, Vxg、Vyg、VzgIt is respectively aircraft ground velocity to exist
Three axle components in earth's axis system, Wxg、Wyg、WzgThe respectively wind speed of aircraft local environment three axle components in earth's axis system;
Aircraft airspeed at three axle components of body axle system is,
In formula, Vxt、Vyt、VztIt is respectively the aircraft airspeed three axle components in body axle system;
The specific algorithm of aircraft airspeed, the elevation angle and yaw angle is,
αV=arctan (Vzt/Vxt)
βV=arcsin (Vyt/VV)
In formula, VVFor aircraft airspeed, αVFor aircraft angle of attack, βVFor aircraft yaw angle.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610374361.1A CN106096091B (en) | 2016-05-31 | 2016-05-31 | A kind of airplane motion analogy method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610374361.1A CN106096091B (en) | 2016-05-31 | 2016-05-31 | A kind of airplane motion analogy method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106096091A true CN106096091A (en) | 2016-11-09 |
CN106096091B CN106096091B (en) | 2019-02-12 |
Family
ID=57230527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610374361.1A Active CN106096091B (en) | 2016-05-31 | 2016-05-31 | A kind of airplane motion analogy method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106096091B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590353A (en) * | 2017-10-19 | 2018-01-16 | 郑州云海信息技术有限公司 | The cluster of the server of air turbulence field stimulation method and use KNL processors |
CN110585703A (en) * | 2019-09-18 | 2019-12-20 | 腾讯科技(深圳)有限公司 | Method, device, equipment and storage medium for controlling flight of virtual aircraft |
CN110709914A (en) * | 2017-03-31 | 2020-01-17 | Cae有限公司 | Continuous monitoring of models in interactive computer simulation stations |
CN111222082A (en) * | 2020-01-15 | 2020-06-02 | 南京航空航天大学 | Method for identifying position and speed of ascending airflow suitable for unmanned aerial vehicle |
CN111240297A (en) * | 2018-11-28 | 2020-06-05 | 中国科学院沈阳自动化研究所 | Spacecraft attitude control system fault diagnosis method based on Kalman filter |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698101A (en) * | 2013-12-15 | 2014-04-02 | 中国航空工业集团公司哈尔滨空气动力研究所 | High-load and high-precision pneumatic force measurement device and measurement method |
US20140343907A1 (en) * | 2013-05-15 | 2014-11-20 | Thales | Method and system for building at least one aircraft guideline in an airport navigation network |
KR101513099B1 (en) * | 2014-01-29 | 2015-04-17 | 국방과학연구소 | Vertification system for weapon aiming algorithm |
CN105243173A (en) * | 2015-08-25 | 2016-01-13 | 中国民航科学技术研究院 | Computer virtual environment simulation and check system for performance based navigation flight program |
-
2016
- 2016-05-31 CN CN201610374361.1A patent/CN106096091B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140343907A1 (en) * | 2013-05-15 | 2014-11-20 | Thales | Method and system for building at least one aircraft guideline in an airport navigation network |
CN103698101A (en) * | 2013-12-15 | 2014-04-02 | 中国航空工业集团公司哈尔滨空气动力研究所 | High-load and high-precision pneumatic force measurement device and measurement method |
KR101513099B1 (en) * | 2014-01-29 | 2015-04-17 | 국방과학연구소 | Vertification system for weapon aiming algorithm |
CN105243173A (en) * | 2015-08-25 | 2016-01-13 | 中国民航科学技术研究院 | Computer virtual environment simulation and check system for performance based navigation flight program |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110709914A (en) * | 2017-03-31 | 2020-01-17 | Cae有限公司 | Continuous monitoring of models in interactive computer simulation stations |
CN110709914B (en) * | 2017-03-31 | 2021-02-02 | Cae有限公司 | Continuous monitoring of models in interactive computer simulation stations |
CN107590353A (en) * | 2017-10-19 | 2018-01-16 | 郑州云海信息技术有限公司 | The cluster of the server of air turbulence field stimulation method and use KNL processors |
CN111240297A (en) * | 2018-11-28 | 2020-06-05 | 中国科学院沈阳自动化研究所 | Spacecraft attitude control system fault diagnosis method based on Kalman filter |
CN110585703A (en) * | 2019-09-18 | 2019-12-20 | 腾讯科技(深圳)有限公司 | Method, device, equipment and storage medium for controlling flight of virtual aircraft |
CN111222082A (en) * | 2020-01-15 | 2020-06-02 | 南京航空航天大学 | Method for identifying position and speed of ascending airflow suitable for unmanned aerial vehicle |
CN111222082B (en) * | 2020-01-15 | 2022-01-04 | 南京航空航天大学 | Method for identifying position and speed of ascending airflow suitable for unmanned aerial vehicle |
Also Published As
Publication number | Publication date |
---|---|
CN106096091B (en) | 2019-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106096091B (en) | A kind of airplane motion analogy method | |
CN105182989B (en) | A kind of aspect control method under the influence of wind field | |
Lu et al. | Real-time simulation system for UAV based on Matlab/Simulink | |
CN109612676B (en) | Pneumatic parameter back calculation method based on flight test data | |
CN106527122B (en) | Stratospheric airship height-lock control nonlinear PID controller method | |
CN111125935B (en) | Simulation system for near space aircraft | |
CN110471313B (en) | Flight simulation subsystem of simulation aircraft | |
CN106444836B (en) | It is a kind of without control sounding rocket Anti-interference Design method | |
CN108153330A (en) | Unmanned aerial vehicle three-dimensional track self-adaptive tracking method based on feasible region constraint | |
CN102425980A (en) | Control method for realizing overload autopilot by using accelerometer | |
CN107367941A (en) | Hypersonic aircraft angle of attack observation procedure based on non-linear gain | |
Tao et al. | Modeling and control of a powered parafoil in wind and rain environments | |
CN105759609A (en) | Carrier-based aircraft autonomous landing method based on explicitly nonlinear model predictive control | |
Wolowicz | Considerations in the determination of stability and control derivatives and dynamic characteristics from flight data | |
CN106570242B (en) | The big dynamic pressure monoblock type radome fairing high-speed separation fluid structurecoupling emulation mode in low latitude | |
CN108008645A (en) | Six-degree-of-freedom simulation modeling method | |
CN106774385A (en) | A kind of dirigible spot hover control method of use adaptive variable structure | |
Li et al. | Simulation method for wind tunnel based virtual flight testing | |
Shen et al. | A 6DOF mathematical model of parachute in Mars EDL | |
Smith et al. | Disturbance observer based control for gust alleviation of a small fixed-wing UAS | |
Caverly et al. | Dynamic modeling, trajectory optimization, and control of a flexible kiteplane | |
Wang et al. | Dynamics modeling and simulation of multi-rotor UAV based on the composite wind field model | |
Chen et al. | Modular aircraft simulation platform based on Simulink | |
Gonyea et al. | Aerodynamic stability and performance of next-generation parachutes for Mars descent | |
Ma et al. | Helicopter nonlinear dynamic inversion flight control model design |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |