CN103576554B - The pneumatic error model component of aircraft based on demand for control, grading design method - Google Patents

The pneumatic error model component of aircraft based on demand for control, grading design method Download PDF

Info

Publication number
CN103576554B
CN103576554B CN201310549486.XA CN201310549486A CN103576554B CN 103576554 B CN103576554 B CN 103576554B CN 201310549486 A CN201310549486 A CN 201310549486A CN 103576554 B CN103576554 B CN 103576554B
Authority
CN
China
Prior art keywords
deviation
rudder
rolling
pitching
moment coefficient
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
Application number
CN201310549486.XA
Other languages
Chinese (zh)
Other versions
CN103576554A (en
Inventor
王颖
闵昌万
王毓栋
黄兴李
刘秀明
吴小华
刘全军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Launch Vehicle Technology CALT
Beijing Institute of Near Space Vehicles System Engineering
Original Assignee
China Academy of Launch Vehicle Technology CALT
Beijing Institute of Near Space Vehicles System Engineering
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Academy of Launch Vehicle Technology CALT, Beijing Institute of Near Space Vehicles System Engineering filed Critical China Academy of Launch Vehicle Technology CALT
Priority to CN201310549486.XA priority Critical patent/CN103576554B/en
Publication of CN103576554A publication Critical patent/CN103576554A/en
Application granted granted Critical
Publication of CN103576554B publication Critical patent/CN103576554B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

The present invention relates to the pneumatic error model component of a kind of aircraft based on demand for control, grading design method, comprise successively: one, Real-time Collection parameter; Two, obtain the raw pitching moment coefficient increment of pitching rudder oblique presentation; Obtain raw pitching moment coefficient increment, yawing moment coefficient increment, the rolling moment coefficient increment of rolling rudder oblique presentation; Three, obtain pitching and be often worth a deviation; Pitching steerage item deviation; The pitching moment coefficient deviation that rolling rudder brings out partially; The normal value deviation of driftage; The stable deviation of driftage; The yawing moment coefficient deviation that rolling rudder brings out partially; Roll and be often worth a deviation; Roll stabilization item deviation; Rolling steerage item deviation; Four, set up pitching moment coefficient deviation; Yawing moment coefficient deviation; Rolling moment coefficient deviation; Five, divide three grades of pneumatic error models. The present invention can balance master-plan in the design difficulty of each subsystem, accelerate overall plan closure, fully indicate flight test risk point, dwindle reference be provided for pneumatic error magnitude after flight test.

Description

The pneumatic error model component of aircraft based on demand for control, grading design method
Technical field
The present invention relates to the pneumatic error model component of a kind of aircraft based on demand for control, grading design sideMethod, particularly relates to one and realizes the severe pattern of control characteristic and draw inclined to one side combination, strengthens ground indication and controlsThe pneumatic error model component of aircraft based on demand for control of ability, grading design method.
Background technology
Pneumatic error model is most important group that in control system analysis, design and emulation, error state is chosenBecome part, framework, value and the occupation mode of pneumatic error model determined Control System Design resultReliability, therefore need to carry out based on demand for control the research of pneumatic error model design and use.
For hypersonic class glide vehicle, mainly rely on aerodynamic force control to realize gliding flight, needPneumatic error model is carried out to meticulous depiction. Traditional pneumatic error model form is " normal value item+proportional "Form, as an example of normal force coefficient deviation example explanation, be generally | Δ CN|≤(0.05+15%·|CN|) shown in form.There is certain problem in this kind of form, cannot meet Control System Design to stability, steerage, intersectionThe pneumatic increment item such as derivative draw inclined to one side demand, cannot realize control system key parameter limit of error stateExamination. For fully exciting the more severe pattern of hypersonic glide vehicle maneuverability and stability, withStrengthen the ability of ground indication and state control. According to aircraft ontological property novelty, a kind of component is proposedThe pneumatic error model of formula, this kind of model form departure is divided into " normal value+stablize item+steerage item+Control coupling terms " component form explain, can realize drawing partially of each component opposed polarity, form and dislikeThe assembled state of bad pattern.
On the basis of having set up at pneumatic error model, rationally effectively use and Assessment, on the one handThe design difficulty of each subsystem in can balance master-plan, accelerates overall plan closure, in addition, and can be abundantIndication flight test risk point, and provide reference for dwindling of pneumatic error magnitude after flight test. ForThe use of pneumatic error model and Assessment, general way is for one-level error band, carries out the limitDraw inclined to one side combination to carry out emulation, require all can realize under total dispersion assembled state stabilized flight, this kind of way warpThe problem often facing is to realize error source design and the balance of controlling design, and scheme is difficult to closure. ForThe closure of implementation, pneumatic error model need to be extremely dwindling of institute's energy to the greatest extent, control ability need to extremely be use up instituteThe enhancing of energy, its consequence is to have strengthened pneumatic error ground to indicate not enough risk, and for main dependenceAirvane is realized the aircraft of gliding flight, exists airvane control ability to excavate inadequate situation (thisIn be espespecially furnished with the aircraft of appearance control dynamical system, be the pneumatic error that adaptation is larger, in airvane controlIn the situation that ability processed not yet plays one's part to the full, adopt RCS to assist control).
Based on above pneumatic error model control application in exist two aspect problem, therefore need badly and provideThe pneumatic error model component of a kind of novel aircraft based on demand for control, grading design method.
Summary of the invention
The technical problem to be solved in the present invention is to provide one and realizes the severe pattern of control characteristic and draw inclined to one side combination,The indication of reinforcement ground and control ability is that tackling key problem order is set up in the convergence of pneumatic error size after flight testMark, the pneumatic error model of balance is realized difficulty and appearance control design difficulty, indicates that in advance pneumatic error is uncertainThe pneumatic error model component of aircraft based on demand for control, the classification of the stabilized flight risk that property is brought are establishedMeter method.
For solving the problems of the technologies described above, the pneumatic error model of a kind of aircraft based on demand for control of the present inventionComponent, grading design method, comprise the following steps successively:
The first step, Real-time Collection flying height Hi, flight Mach number Ma, flying drilling angle α, flight is breakked awayAngle β, aircraft pitching rudder is inclined to one sideThe inclined to one side δ of vehicle yaw rudderψ, the inclined to one side δ of Vehicle Roll rudderγ
Determine normal value coefficient a, b, k according to error band envelope;
Second step, the raw moment coefficient of acquisition pitching rudder oblique presentationObtain 0 °The moment coefficient that pitching rudder is on the lower sideSubtract each other and obtain the raw pitching of pitching rudder oblique presentationMoment coefficient increment
Obtain the raw moment coefficient of rolling rudder oblique presentationObtain 0 ° of rolling rudder inclined to one sideUnder moment coefficientSubtract each other and obtain rolling the raw pitching moment coefficient of rudder oblique presentationIncrementThe yawing moment coefficient increment that the oblique presentation of rolling rudder is rawThe rolling power that the oblique presentation of rolling rudder is rawMoment coefficient increment
The 3rd step, acquisition pitching are often worth a deviation delta Cmz0=±(a*α+b);
Pitching steerage item deviation
The pitching moment coefficient deviation that rolling rudder brings out partially
The normal value of a driftage deviation delta Cmy0=±(a*α+b);
The stable deviation of driftage ΔC my β = ± ( a * α + b ) ;
The yawing moment coefficient deviation that rolling rudder brings out partially
The normal value of rolling a deviation delta Cmx0=±(a*α+b);
Roll stabilization item deviation ΔC mx β = ± ( a * α + b ) ;
Rolling steerage item deviation
The 4th step, set up pitching moment coefficient deviation
Set up yawing moment coefficient deviation ΔC my = ΔC my 0 + ΔC my β · β + ΔC my _ δ γ ;
Set up rolling moment coefficient deviation ΔC mx = ΔC mx 0 + ΔC mx β · β + ΔC mx _ δ γ ;
The 5th step, after airvane starts to control, rely on the control of airvane independence can realize total dispersion state steadyFixed, adaptable in the case pneumatic error model is the pneumatic error model of the first order;
On first order error model basis, each component terms is carried out to amplification in various degree, when airvane risesAfter control, rely on RCS and the compound control of airvane to realize total dispersion in stable condition, adaptable in the casePneumatic error model is the pneumatic error model in the second level;
On error model basis, the second level, each component terms is carried out to amplification in various degree, when airvane risesAfter control, rely on RCS and the compound control of airvane, total dispersion assembled state is carried out to the probability of stability and know the real situation.
Adopt the present invention can excavate hypersonic glide vehicle evaluation of layout or control design and more dislikeBad state, realizes control characteristic and badly draws inclined to one side combination, fully excavated earth design risk and control ability;
The present invention breaks through the thinking that in traditional appearance control design, deviation is used, and makes design too not guard and advance rashly,For the design of pneumatic error and appearance control design have all proposed research objective, and balance pneumatic error model realExisting difficulty and appearance control design difficulty, on the one hand can balance master-plan in the design difficulty of each subsystem, addFast overall plan closure, in addition, can fully indicate flight test risk point, and be pneumatic after flight testDwindling of error magnitude provides reference.
Detailed description of the invention
For pitch channel, heart indication deviation is longitudinally pressed in the reaction of normal value. For driftage, roll channel,The reaction structure manufacture of normal value, installation, the asymmetric deviation of profile that the factor such as ablation produces in-flight. RightIn pitch channel, stability item reaction pitching static stability deviation. For driftage, roll channel, stableProperty item deviation is reacted respectively static yawing stability deviation, rolling static stability deviation. Inclined to one side for steerage itemPoor, react respectively the pneumatic rudder ability of triple channel deviation. For controlling coupling terms deviation, react respectively a certainThe deviation of other passage aerodynamic moment size that the pneumatic rudder of passage brings out.
For pitch channel, generally need to consider longitudinally to press heart deviation, steerage deviation, because of pitching rudder partially with rollIt is inclined to one side that moving rudder shares elevator partially, and the pitching moment coefficient that rolling rudder brings out partially changes can not ignore, so needThe special pitching moment coefficient bias term of considering that rolling is brought out.
For jaw channel, generally need to consider a normal value deviation, a stable deviation, steerage deviation, forHave the aircraft of larger trailing edge rudder layout, the yawing moment coefficient partially bringing out because of rolling rudder changes and can not neglectSlightly, so need the special yawing moment coefficient bias term of considering that rolling is brought out.
Specifically, the present invention comprises the following steps successively:
The first step, Real-time Collection flying height Hi, flight Mach number Ma, flying drilling angle α, flight is breakked awayAngle β, aircraft pitching rudder is inclined to one sideThe inclined to one side δ of vehicle yaw rudderψ, the inclined to one side δ of Vehicle Roll rudderγ
Determine normal value coefficient a, b, k according to error band envelope;
Second step, the raw moment coefficient of acquisition pitching rudder oblique presentationObtain 0 °The moment coefficient that pitching rudder is on the lower sideSubtract each other and obtain the raw pitching of pitching rudder oblique presentationMoment coefficient increment
Obtain the raw moment coefficient of rolling rudder oblique presentationObtain 0 ° of rolling rudder inclined to one sideUnder moment coefficientSubtract each other and obtain rolling the raw pitching moment coefficient of rudder oblique presentationIncrementThe yawing moment coefficient increment that the oblique presentation of rolling rudder is rawThe rolling power that the oblique presentation of rolling rudder is rawMoment coefficient increment
The 3rd step, for a normal value deviation delta Cmz0, the variation characteristic of its deviation shows as the primary line of the angle of attackProperty function, in the situation that pressing the heart amount of moving forward and backward consistent, obtain pitching and be often worth a deviationΔCmz0=±(a*α+b);
For steerage item deviation, the form of available rudder partial increment percentage, obtains pitching steerage item deviation
The pitching moment coefficient bias term of partially bringing out for rolling rudder, according to recognizing aircraft ontological propertyKnow the pitching moment coefficient increment size that deviation size is put with pitching rudder off normal and the rudder that rolls produces that it brings outRelevant, obtain the pitching moment coefficient deviation that rolling rudder brings out partially
The normal value of a driftage deviation delta Cmy0=±(a*α+b);
The stable deviation of driftage ΔC my β = ± ( a * α + b ) ;
The yawing moment coefficient deviation that rolling rudder brings out partially
The normal value of rolling a deviation delta Cmx0=±(a*α+b);
Roll stabilization item deviation ΔC mx β = ± ( a * α + b ) ;
Rolling steerage item deviation
The 4th step, set up pitching moment coefficient deviation
The in the situation that of the directionless rudder control of jaw channel, set up yawing moment coefficient deviation ΔC my = ΔC my 0 + ΔC my β · β + ΔC my _ δ γ ;
Set up rolling moment coefficient deviation ΔC mx = ΔC mx 0 + ΔC mx β · β + Δ C mx _ δ γ ;
The 5th step, on the definite basis of error model framework, the adaptable pneumatic mistake of Comprehensive Control abilityResidual quantity value and emulation mode, determine and error band be divided into the error band value of appropriate level to three grades and establishMeter and examination.
One: first determine the pneumatic error model of the first order. Consider pneumatic deviation and other design deviation,Carry out the total dispersion limit and draw inclined to one side combination to carry out emulation, after airvane starts to control, rely on airvane independently to controlCan realize total dispersion in stable condition, adaptable in the case pneumatic error model is the pneumatic mistake of the first orderDifferential mode type;
Two: on the pneumatic error model of first order basis, determine the pneumatic error model in the second level, this level mistakeThe amplification that differential mode type carries out in various degree each component terms on first order error model basis obtains. ComprehensivelyConsider this grade of pneumatic deviation and other design deviation, carry out the total dispersion limit and draw inclined to one side combination to carry out emulation, whenAfter airvane starts to control, rely on RCS and the compound control of airvane to realize total dispersion in stable condition, in the caseAdaptable pneumatic error model is the pneumatic error model in the second level;
Three: on pneumatic error model basis, the second level, determine the pneumatic error model of the third level, this level mistakeThe amplification that differential mode type carries out in various degree each component terms on error model basis, the second level obtains. ComprehensivelyConsider this grade of pneumatic deviation and other design deviation, pneumatic deviation is carried out to Sampling, to other deviationCarry out limit value, after airvane starts to control, rely on RCS and the compound control of airvane, to total dispersion groupThe state of closing carries out the probability of stability and knows the real situation.

Claims (1)

1. the pneumatic error model component of the aircraft based on demand for control, a grading design method, successivelyComprise the following steps:
The first step, Real-time Collection flying height Hi, flight Mach number Ma, flying drilling angle α, flight is breakked awayAngle β, aircraft pitching rudder is inclined to one sideThe inclined to one side δ of vehicle yaw rudderψ, the inclined to one side δ of Vehicle Roll rudderγ
Determine normal value coefficient a, b, k according to error band envelope;
Second step, the raw moment coefficient of acquisition pitching rudder oblique presentationObtain 0 °The moment coefficient that pitching rudder is on the lower sideSubtract each other and obtain the raw pitching of pitching rudder oblique presentationMoment coefficient increment
Obtain the raw moment coefficient of rolling rudder oblique presentationObtain 0 ° of rolling rudder inclined to one sideUnder moment coefficientSubtract each other and obtain rolling the raw pitching moment coefficient of rudder oblique presentationIncrementThe yawing moment coefficient increment that the oblique presentation of rolling rudder is rawThe rolling power that the oblique presentation of rolling rudder is rawMoment coefficient increment
The 3rd step, acquisition pitching are often worth a deviation delta Cmz0=±(a*α+b);
Pitching steerage item deviation
The pitching moment coefficient deviation that rolling rudder brings out partially
The normal value of a driftage deviation delta Cmy0=±(a*α+b);
The stable deviation of driftage ΔC m y β = ± ( a * α + b ) ;
The yawing moment coefficient deviation that rolling rudder brings out partially
The normal value of rolling a deviation delta Cmx0=±(a*α+b);
Roll stabilization item deviation ΔC m x β = ± ( a * α + b ) ;
Rolling steerage item deviation
The 4th step, set up pitching moment coefficient deviation
Set up yawing moment coefficient deviation ΔC m y = ΔC m y 0 + ΔC m y β · β + ΔC m y _ δ γ ;
Set up rolling moment coefficient deviation ΔC m x = ΔC m x 0 + ΔC m x β · β + ΔC m x _ δ γ ;
The 5th step, after airvane starts to control, rely on the control of airvane independence can realize total dispersion state steadyFixed, adaptable in the case pneumatic error model is the pneumatic error model of the first order;
On first order error model basis, each component terms is carried out to amplification in various degree, when airvane risesAfter control, rely on RCS and the compound control of airvane to realize total dispersion in stable condition, adaptable in the casePneumatic error model is the pneumatic error model in the second level;
On error model basis, the second level, each component terms is carried out to amplification in various degree, when airvane risesAfter control, rely on RCS and the compound control of airvane, total dispersion assembled state is carried out to the probability of stability and know the real situation.
CN201310549486.XA 2013-11-07 2013-11-07 The pneumatic error model component of aircraft based on demand for control, grading design method Active CN103576554B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310549486.XA CN103576554B (en) 2013-11-07 2013-11-07 The pneumatic error model component of aircraft based on demand for control, grading design method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310549486.XA CN103576554B (en) 2013-11-07 2013-11-07 The pneumatic error model component of aircraft based on demand for control, grading design method

Publications (2)

Publication Number Publication Date
CN103576554A CN103576554A (en) 2014-02-12
CN103576554B true CN103576554B (en) 2016-05-18

Family

ID=50048571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310549486.XA Active CN103576554B (en) 2013-11-07 2013-11-07 The pneumatic error model component of aircraft based on demand for control, grading design method

Country Status (1)

Country Link
CN (1) CN103576554B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104155983B (en) * 2014-08-08 2015-05-20 北京航天自动控制研究所 Crosslinking impact assessment method for aerodynamic coupling property between aircraft attitude movement channels
CN104155985B (en) * 2014-08-08 2015-05-20 北京航天自动控制研究所 Crosslinking impact assessment method for inertia coupling property between aircraft attitude movement channels
CN105628333B (en) * 2014-10-29 2018-06-22 北京临近空间飞行器系统工程研究所 A kind of pneumatic error under the conditions of high-altitude High Mach number determines method
CN104914736A (en) * 2015-05-07 2015-09-16 南京航空航天大学 Method for building general simulation model for hypersonic flight vehicles
CN105136422B (en) * 2015-09-10 2017-10-13 中国航天空气动力技术研究院 The method that dummy vehicle sideslip angular flexibility is corrected in wind tunnel test
CN111985169B (en) * 2020-08-26 2024-04-02 北京中科宇航技术有限公司 Paraxial symmetric aircraft aerodynamic characteristic modeling method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4970638A (en) * 1988-10-31 1990-11-13 The United States Of America As Represented By The Secretary Of The Air Force Control of unknown systems via deconvolution
CN101660914A (en) * 2009-08-19 2010-03-03 南京航空航天大学 Airborne starlight of coupling inertial position error and independent navigation method of inertial composition
CN103063216A (en) * 2013-01-06 2013-04-24 南京航空航天大学 Inertial and celestial combined navigation method based on star coordinate modeling
CN103245359A (en) * 2013-04-23 2013-08-14 南京航空航天大学 Method for calibrating fixed errors of inertial sensor in inertial navigation system in real time
CN103267531A (en) * 2013-04-23 2013-08-28 上海卫星工程研究所 Method for high-precision compensation of fiber-optic gyroscope random error

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2423377B (en) * 2002-12-09 2007-04-18 Georgia Tech Res Inst Adaptive Output Feedback Apparatuses And Methods Capable Of Controlling A Non-Minimum Phase System

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4970638A (en) * 1988-10-31 1990-11-13 The United States Of America As Represented By The Secretary Of The Air Force Control of unknown systems via deconvolution
CN101660914A (en) * 2009-08-19 2010-03-03 南京航空航天大学 Airborne starlight of coupling inertial position error and independent navigation method of inertial composition
CN103063216A (en) * 2013-01-06 2013-04-24 南京航空航天大学 Inertial and celestial combined navigation method based on star coordinate modeling
CN103245359A (en) * 2013-04-23 2013-08-14 南京航空航天大学 Method for calibrating fixed errors of inertial sensor in inertial navigation system in real time
CN103267531A (en) * 2013-04-23 2013-08-28 上海卫星工程研究所 Method for high-precision compensation of fiber-optic gyroscope random error

Also Published As

Publication number Publication date
CN103576554A (en) 2014-02-12

Similar Documents

Publication Publication Date Title
CN103576554B (en) The pneumatic error model component of aircraft based on demand for control, grading design method
CN105752314B (en) A kind of high-altitude low speed Natural Laminar Flow high-lift airfoil
CN102707624B (en) Design method of longitudinal controller region based on conventional aircraft model
CN102749181B (en) Wind tunnel test method on basis of momentum theory
CN103587680B (en) A kind of aircraft side-slipping turn control method
CN103587681A (en) Hypersonic speed aircraft control method capable of suppressing constant deviation influence of sideslip angle signal
CN103625637B (en) A kind of large aircraft side direction fitful wind mitigation process
CN102117362B (en) Light airplane horizontal tail design load determination method under slipstream influence
CN104881035A (en) Method and system for aircraft control coupling compensation, and method and system for attitude motion control
CN107463746A (en) Fuselage bulkhead circumference stress computational methods under a kind of airtight load
CN102052266B (en) After-load blunt trailing edge wing profile designed based on sharp trailing edge wing profile
CN104401503B (en) A kind of trim rudder drift angle defining method being applicable to many rudders face ablated configuration device
Hufford et al. Viscous flow around marine propellers using boundary-layer strip theory
Shindo Simplified tunnel correction method
Lamar Prediction of vortex flow characteristics of wings at subsonic and supersonic speeds
CN104155985B (en) Crosslinking impact assessment method for inertia coupling property between aircraft attitude movement channels
CN106168529B (en) A kind of wind tunnel test methods for correcting aircraft rigid body lift coefficient
CN106950982A (en) Reentry vehicle attitude control dynamical system high-altitude torque factor discrimination method
Rein Measurements of aerodynamic forces and moments on the DLR-F17E model in low-and high-speed flows
CN110160737A (en) The accurate acquisition methods of elastic pneumatic data based on engineering panel method
Huffman et al. Subsonic longitudinal and lateral-directional static aerodynamic characteristics for a model with swept back and swept forward wings
Nichols Analysis and compilation of missile aerodynamic data. Volume 1: Data presentation and analysis
Peristy et al. Wake Modelling with Embedded Lateral and Directional Stability Sensitivity Analysis for Aerial Refuelling or Formation Flight
Ferris Static aerodynamic characteristics of a model with a 17 percent thick supercritical wing
CN117669052A (en) Method for determining maximum hinge moment of spoiler of transport aircraft

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant