CN102320378A - 多操纵面飞机的一种均衡操纵分配方法 - Google Patents
多操纵面飞机的一种均衡操纵分配方法 Download PDFInfo
- Publication number
- CN102320378A CN102320378A CN201110165732A CN201110165732A CN102320378A CN 102320378 A CN102320378 A CN 102320378A CN 201110165732 A CN201110165732 A CN 201110165732A CN 201110165732 A CN201110165732 A CN 201110165732A CN 102320378 A CN102320378 A CN 102320378A
- Authority
- CN
- China
- Prior art keywords
- msub
- mover
- msubsup
- mrow
- mtd
- 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
- 238000000034 method Methods 0.000 title claims abstract description 88
- 239000011159 matrix material Substances 0.000 claims abstract description 131
- 238000010606 normalization Methods 0.000 claims abstract description 47
- 238000012163 sequencing technique Methods 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 314
- 238000005096 rolling process Methods 0.000 claims description 98
- 230000007704 transition Effects 0.000 claims description 29
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000032683 aging Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000001965 increasing effect Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 241001398959 Onychorhynchus coronatus Species 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
Images
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明提出的多操纵面飞机的一种均衡操纵分配方法,主要包括:将多操纵面的控制分配问题转化为虚拟操纵面与真实操纵面的等效控制问题、确定对真实操纵面进行标陈化及控制效率均衡分配的基本方法、将真实操纵面分类、进行标陈化并按照标陈化控制效率的绝对值排序、确定虚拟操纵面的主导控制效率、将标陈化操纵面分级、计算分配系数等几个步骤。本发明将多操纵面的控制等效为特定控制效率及偏转范围限制的虚拟操纵面的控制,简化了控制律设计,避免了操纵面角速率超限的问题,加权均衡分配矩阵计算简便,综合考虑了操纵面的偏转范围、控制效率、控制级别、权限系数,充分利用了有效操纵面,有效避免了某些操纵面负担过多的控制效率而引起的饱和。
Description
技术领域
本发明属于飞行控制领域,具体涉及多操纵面飞机的一种均衡操纵分配方法。
背景技术
传统的飞机往往采用副翼产生滚转运动,方向舵产生偏航运动,升降舵产生俯仰运动。为了提高机动性能及可靠性,当今的飞机往往采用多操纵面布局。例如,美国先进战斗机YF-22有12个独立的操纵面(平尾、副翼、方向舵、襟副翼、前缘襟翼、俯仰推力矢量各两个操纵面)。由于产生滚转、偏航及俯仰力矩的操纵面不唯一,多操纵面布局提高了飞机操纵的灵活性及飞机的控制余度,从而增强了飞行的控制性能及安全性能。
多操纵面的冗余及不同操纵面间的耦合给飞行员的驾驶和飞机自动驾驶仪的控制律设计带来很大的难度,需要设计一种分配方法将控制指令分配到各个操纵面。当前许多操纵面的控制分配方法是将多操纵面的控制分配问题归结为等效期望力矩的概念,即在操纵面的偏转范围及角速率限制下,根据控制效率矩阵来求解所有操纵面的偏转角,使得操纵面在该偏转角下产生的力矩等于期望力矩。工程上对于多操纵面的控制分配往往采用串接链式分配方法,该方法是将操纵面进行分组,不同组的操纵面对应不同的优先级,优先级高的操纵面首先分担期望力矩,当该组中某些操纵面的偏转角达到门限值时依次启用下一优先级的操纵面。串接链式分配方法物理概念清晰,但是对于操纵面的分组及优先级的划分往往只是根据工程经验或者控制效率矩阵给出,由于不同操纵面的偏转范围不同,控制效率矩阵不能体现操纵面的饱和度,优先级高的操纵面往往导致较大的偏转甚至饱和,且没有充分利用优先级低的操纵面。加权伪逆方法基于操纵面的偏转范围及转动角速率限制设计加权矩阵,采用伪逆的方法由期望力矩解算所有操纵面的偏转角。加权矩阵选取要求设计者有比较丰富的经验,某些参数可能需要试凑与调整,增加了设计的复杂性。线性规划、直接几何分配及二次规划的方法是根据操纵面的偏转范围、角速率限制及期望力矩建立指标函数,利用寻优的方法通过复杂的计算得到满足期望指标的所有操纵面的偏转角。这些方法数学性强,实时性难以保证,物理概念不够清晰,在实际工程应用中应用有一定的限制。
可见,当前许多控制分配的方法是根据控制效率矩阵由期望力矩求解所有操纵面的偏转角,计算量大,需要考虑可能出现的偏转角速率超限的情况,驾驶员的驾驶指令或者飞机自动驾驶仪的控制指令需要转换成期望力矩指令,增加了控制律设计的复杂性,工程应用受到限制。
发明内容
针对现有技术中存在的问题,本发明提出多操纵面飞机的一种均衡操纵分配方法,从工程应用出发,将多操纵面的控制分配问题转化为虚拟操纵面与真实操纵面的等效控制问题,确定标称化操纵面及控制效率均衡分配的基本方法。根据操纵面的滚转、偏航、俯仰控制效率的大小,将真实操纵面划分为滚转操纵面、偏航操纵面与俯仰操纵面。分别对于滚转操纵面、偏航操纵面、俯仰操纵面进行标称化,并根据标称化控制效率设计虚拟操纵面的主导的控制效率,且将陈化操纵面划分为主、副、辅三级,根据设定的主、副权限系数确定主、副、辅各级中每个标称化操纵面期望分担的主导控制效率,进而求解虚拟操纵面向真实操纵面的分配系数,建立加权均衡分配矩阵,并得到虚拟操纵面矢量的非主导控制效率。本发明的加权均衡分配矩阵计算简便,综合考虑了操纵面的偏转范围、控制效率、控制级别、权限系数,充分利用了有效操纵面,有效避免了某些操纵面负担过多的控制效率而引起的饱和。
本发明提出的多操纵面飞机的一种均衡操纵分配方法,具体包括以下几个步骤:
步骤一:将多操纵面的控制分配问题转化为虚拟操纵面与真实操纵面的等效控制问题:
飞机的n个真实操纵面组成真实操纵面矢量U=[u1 u2 …ui …un-1 un]′,其中ui表示第i个真实操纵面,真实操纵面偏转角矢量为δ=[δ1 δ2 …δi …δn-1 δn]′,其中δi表示第i个真实操纵面的偏转角,且控制效率矩阵为 其中,(i=1,2,…,n-1,n)表示第i个真实操纵面ui的控制效率矢量,且 分别表示第i个真实操纵面ui的滚转控制效率、偏航控制效率、俯仰控制效率,且第i个真实操纵面ui(i=1,2,…,n-1,n)的偏转角δi的偏转范围为[-Ri,Ri],角速率限制为[-ρi,ρi];
建立虚拟操纵面矢量 其中,表示虚拟副翼,虚拟副翼偏转范围为[-Rk Rk];表示虚拟方向舵,虚拟方向舵偏转范围为[-Rl Rl];表示虚拟升降舵,虚拟升降舵偏转范围为[-Rm Rm];虚拟操纵面偏转角矢量 其中分别表示虚拟副翼偏转角、虚拟方向舵偏转角、虚拟升降舵偏转角;虚拟操纵面矢量的控制效率矩阵为:
其中,为虚拟副翼的滚转控制效率,为虚拟副翼的偏航控制效率,为虚拟副翼的俯仰控制效率,为虚拟方向舵的滚转控制效率,为虚拟方向舵的偏航控制效率,为虚拟方向舵的俯仰控制效率,为虚拟升降舵的滚转控制效率,为虚拟升降舵的偏航控制效率,为虚拟升降舵的俯仰控制效率;
真实操纵面的偏转角矢量δ与虚拟操纵面偏转角矢量满足关系K为加权均衡分配矩阵,且虚拟操纵面矢量控制的效率矩阵满足控制器根据驾驶员的驾驶指令或自动驾驶仪期望的轨迹、姿态指令Cmd及飞机当前状态矢量X计算虚拟操纵面偏转角矢量 根据加权均衡矩阵K计算真实操纵面偏转角矢量真实操纵面偏转角矢量根据真实操纵面矢量的控制效率矩阵B,计算得到期望的力矩矢量其中 虚拟操纵面矢量根据虚拟操纵面矢量的控制效率矩阵计算得到力矩矢量为
步骤二:确定对真实操纵面进行标称化及控制效率加权均衡分配的基本方法:
定义真实操纵面ui(i=1,2,…,n-1,n)的偏转角δi时对应的真实操纵面的饱和度Si为Si=fabs(δi/Ri),其中函数fabs(x)表示小数x的绝对值;
将控制效率矢量为偏转范围为[-Ri Ri]的真实操纵面ui(i=1,2,…,n-1,n)对于偏转范围[-Rj Rj]进行标称化,得到标称化滚转控制效率的绝对值、标称化偏航控制效率的绝对值、标称化俯仰控制效率的绝对值分别为 其中函数fabs(x)表示小数x的绝对值;
设n个真实操纵面中有m个真实操纵面是属于A转动方向的操纵面,其中tj∈[1,2,…,n-1,n],且真实操纵面(j=1,2,…,m-1,m)偏转角偏转范围分别为 控制效率矢量为且在A转动方向的控制效率为设某虚拟操纵面与m个真实操纵面在A转动方向的控制等效,虚拟操纵面的偏转角为虚拟操纵面在A转动方向的控制效率为偏转角的偏转范围为[-Rh Rh];真实操纵面(j=1,2,…,m-1,m)对于偏转范围[-Rh Rh]进行标称化,标称化操纵面的标称化控制效率矢量:
在A转动方向的标称化控制效率:
标称化偏转角:
步骤三:根据真实操纵面的滚转、偏航及俯仰控制效率的大小,将所有真实操纵面划分为滚转操纵面、偏航操纵面与俯仰操纵面:
滚转操纵面的数目计为n1,偏航操纵面的数目计为n2,俯仰操纵面的数目计为n3,n1+n2+n3=n;
步骤四:分别将滚转操纵面、偏航操纵面、俯仰操纵面对于所对应的虚拟操纵面的偏转范围进行标称化,并按照标称化控制效率的绝对值排序:
确定虚拟副翼的滚转控制效率与所有滚转操纵面中标称化滚转控制效率的绝对值最大的操纵面的滚转控制效率同号。将所有滚转操纵面对于虚拟副翼的偏转范围[-Rk Rk]进行标称化,并按照标称化滚转控制效率由大到小排序,排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面为(pj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的滚转操纵面的控制效率矢量及滚转控制效率,表示虚拟副翼的滚转控制效率,滚转操纵面的偏转范围为 且 分别表示标称化滚转操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率,建立标称化滚转操纵面矢量 其中(pj∈[1,2,…,n-1,n])表示按照标称化滚转控制效率的绝对值排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面,标称化滚转操纵面矢量的控制效率矩阵 其中表示按照标称化滚转控制效率的绝对值排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面(pj∈[1,2,…,n-1,n])的标称化控制效率矢量;
确定虚拟方向舵的偏航控制效率与所有偏航操纵面中标称化偏航控制效率的绝对值最大的操纵面的偏航控制效率同号。将所有偏航操纵面对于虚拟方向舵的偏转范围[-Rl Rl]进行标称化,并按照标称化偏航控制效率的绝对值由大到小排序,排序中第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面为(qj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的偏航操纵面的控制效率矢量及偏航控制效率,表示虚拟方向舵的偏航控制效率,偏航操纵面的偏转范围为 且 分别表示标称化偏航操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率,建立标称化偏航操纵面矢量 其中(qj∈[1,2,…,n-1,n])表示按照标称化偏航控制效率的绝对值排序中第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面,标称化偏航操纵面矢量的控制效率矩阵 其中表示按照标称化偏航控制效率的绝对值排序中第j(j=1,2,…n2-1,n2)个标称化偏航操纵面(qj∈[1,2,…,n-1,n])的标称化控制效率矢量;
确定虚拟升降舵的俯仰控制效率与所有俯仰操纵面中标称化俯仰控制效率的绝对值最大的操纵面的俯仰控制效率同号。将所有俯仰操纵面对于虚拟升降舵的偏转范围[-Rm Rm]进行标称化,并按照标称化俯仰控制效率的绝对值由大到小排序,排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面为(rj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的俯仰操纵面的控制效率矢量及俯仰控制效率,表示虚拟升降舵的俯仰控制效率,俯仰操纵面的偏转范围为 且 分别表示标称化俯仰操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率;则标称化俯仰操纵面矢量 其中(rj∈[1,2,…,n-1,n])表示按照标称化俯仰控制效率的绝对值排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面,标称化俯仰操纵面矢量的控制效率矩阵为 其中,表示按照标称化俯仰控制效率的绝对值排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面(rj∈[1,2,…,n-1,n])的标称化控制效率矢量;
步骤五:确定虚拟操纵面的主导控制效率:
步骤六:根据标称化控制效率将标称化操纵面划分成主、副、辅三级,设定标称化主操纵面与标称化副操纵面的权限系数,并计算主、副、辅各级标称化操纵面期望分担的主导控制效率:
(1)分别将标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面的划分为主、副、辅三级:
确定标称化滚转副操纵面序号临界值nx2满足nx1<nx2≤n1,nx2是整数且满足:
将位于标称化滚转操纵面矢量中前nx1个位置的操纵面划分为标称化滚转主操纵面;将位于标称化滚转操纵面矢量中从(nx1+1)到nx2个位置的操纵面划分为标称化滚转副操纵面;将位于标称化滚转操纵面矢量中从(nx2+1)到n1个位置的操纵面划分为标称化滚转辅操纵面;
则标称化偏航副操纵面序号临界值ny2满足ny1<ny2≤n2,ny2是整数且满足:
将位于标称化偏航操纵面矢量中前ny1个位置的操纵面划分为标称化偏航主操纵面;将位于标称化偏航操纵面矢量中从(ny1+1)到ny2个位置的操纵面划分为标称化偏航副操纵面;将位于标称化偏航操纵面矢量中从(ny2+1)到n2个位置的操纵面划分为标称化偏航辅操纵面;
则标称化俯仰副操纵面序号临界值nz2满足nz1<nz2≤n3,nz2是整数且满足:
将位于标称化俯仰操纵面矢量中前nz1个位置的操纵面划分为标称化俯仰主操纵面;将位于标称化俯仰操纵面矢量中从(nz1+1)到nz2个位置的操纵面划分为标称化俯仰副操纵面;将位于标称化俯仰操纵面矢量中从(nz2+1)到n3个位置的操纵面划分为标称化俯仰辅操纵面;
(2)分别对标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面设定主、副权限系数,并计算标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面的主、副、辅各级标称化操纵面期望分担的主导控制效率:
按照控制效率均衡分配的方法得到所有的标称化滚转主操纵面所分担的虚拟副翼的滚转控制效率的比率Cxm为:
所有的标称化滚转副操纵面所分担的虚拟副翼的滚转控制效率的比率Cxs为:
所有的标称化滚转辅操纵面所分担的虚拟副翼的滚转控制效率的比率Cxf为:
Cxf=1-Cxm-Cxs
标称化滚转主操纵面的权限系数为Qxm及标称化滚转副操纵面的权限系数为Qxs,则所有的标称化滚转主操纵面期望分担的滚转控制效率则所有的标称化滚转副操纵面期望分担的滚转控制效率所有的标称化滚转辅操纵面期望分担的滚转控制效率
按照控制效率均衡分配的方法,得到所有的标称化偏航主操纵面所分担的虚拟方向舵的偏航控制效率的比率Cym为:
所有的标称化偏航副操纵面所分担的虚拟方向舵的偏航控制效率的比率Cys为:
所有的标称化偏航辅操纵面所分担的虚拟方向舵的偏航控制效率的比率Cyf为:
Cyf=1-Cym-Cys
标称化偏航主操纵面的权限系数为Qym及标称化偏航副操纵面的权限系数为Qys,则所有的标称化偏航主操纵面期望分担的偏航控制效率所有的标称化偏航副操纵面期望分担的偏航控制效率所有的标称化偏航辅操纵面期望分担的偏航控制效率
按照控制效率均衡分配的方法,得到所有的标称化俯仰主操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czm为:
所有的标称化俯仰副操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czs为:
所有的标称化俯仰辅操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czf为:
Czf=1-Czm-Czs
标称化俯仰主操纵面的权限系数为Qzm,标称化俯仰副操纵面的权限系数为Qzs,所有的标称化俯仰主操纵面期望分担的俯仰控制效率所有的标称化俯仰副操纵面期望分担的俯仰控制效率所有的标称化俯仰辅操纵面期望分担的俯仰控制效率
步骤七:采用控制效率均衡分配的方法计算每个标称化操纵面期望分担的主导控制效率:
(1)计算每个标称化滚转操纵面期望分担的滚转控制效率:
采用控制效率均衡分配的方法将所有的标称化滚转副操纵面期望分担的滚转控制效率分配给每个标称化滚转副操纵面,则第j(j=nx1+1,nx1+2,…,nx2-1,nx2)个标称化滚转副操纵面期望分担的滚转控制效率为:
采用控制效率均衡分配的方法将所有的标称化滚转辅操纵面期望分担的滚转控制效率分配每个标称化滚转辅操纵面,则第j(j=nx2+1,nx2+2,…,n1-1,n1)个标称化滚转辅操纵面期望分担的滚转控制效率为:
(2)计算每个标称化偏航操纵面期望分担的偏航控制效率:
采用控制效率均衡分配的方法将所有的标称化偏航副操纵面期望分担的偏航控制效率分配给每个标称化偏航副操纵面,则第j(j=ny1+1,ny1+2,…,ny2-1,ny2)个标称化偏航副操纵面期望分担的偏航控制效率为:
采用控制效率均衡分配的方法将所有的陈化偏航辅操纵面期望分担的偏航控制效率分配给每个标称化偏航辅操纵面,则第j(j=ny2+1,ny2+2,…,n2-1,n2)个标称化偏航辅操纵面期望分担的偏航控制效率为:
(3)计算每个标称化俯仰操纵面期望分担的俯仰控制效率:
采用控制效率均衡分配的方法将所有的标称化俯仰副操纵面期望分担的俯仰控制效率分配给每个标称化俯仰副操纵面,则第j(j=nz1+1,nz1+2,…,nz2-1,nz2)个标称化俯仰副操纵面期望分担的俯仰控制效率为:
采用控制效率均衡分配的方法将所有的标称化俯仰辅操纵面期望分担的俯仰控制效率分配给每个标称化俯仰辅操纵面,则第j(j=nz2+1,nz2+2,…,n3-1,n3)个标称化俯仰辅操纵面期望分担的俯仰控制效率为:
步骤八:计算所有真实操纵面的分配系数,建立加权均衡分配矩阵,求解虚拟操纵面矢量的控制效率矩阵中的非主导控制效率:
(1)计算所有真实操纵面的分配系数:
①计算每个滚转操纵面的分配系数:
由于第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面与标称化之前的滚转操纵面对应,标称化滚转操纵面期望分担的滚转控制效率为滚转操纵面的偏转范围为滚转控制效率为则虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数为:
②计算每个偏航操纵面的分配系数:
由于第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面与标称化之前的偏航操纵面对应,标称化偏航操纵面期望分担的偏航控制效率为偏航操纵面的偏转范围为 偏航控制效率为则虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数为:
③计算每个俯仰操纵面的分配系数:
由于第j(j=1,2,…,n3-1,n3)个标称化俯仰操纵面与标称化之前的俯仰操纵面对应,标称化俯仰操纵面期望分担的俯仰控制效率为俯仰操纵面的偏转范围为 俯仰控制效率为则虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数为:
(2)建立加权均衡分配矩阵,求解虚拟操纵面矢量的控制效率矩阵中的非主导控制效率,从而确定虚拟操纵面矢量的控制效率矩阵:
①虚拟操纵面矢量到真实操纵面矢量的加权均衡分配矩阵:
建立过渡操纵面矢量 过渡操纵面的偏转角矢量 过渡操纵面矢量的控制效率矩阵
其中,(i=1,2,…,n1-1,n1)分别表示滚转操纵面的偏转角及控制效率矢量,(j=1,2,…,n2-1,n2)分别表示偏航操纵面的偏转角及控制效率矢量,(t=1,2,…,n3-1,n3)分别表示滚转操纵面的偏转角及控制效率矢量,虚拟操纵面矢量 对应的虚拟操纵面偏转角矢量 虚拟操纵面矢量的控制效率矩阵 其中,表示虚拟副翼的偏转角及控制效率矢量, 表示虚拟方向舵的偏转角及控制效率矢量,分别表示虚拟升降舵的偏转角及控制效率矢量;则虚拟操纵面偏转角矢量到过渡操纵面偏转角矢量的传递关系为其中,表示虚拟操纵面矢量到过渡操纵面矢量的加权均衡分配矩阵:
对应的真实操纵面的偏转角矢量δ=[δ1 δ2 …δi…δn-1 δn]′,真实操纵面的控制效率矩阵 其中δi、分别表示真实操纵面ui的偏转角、控制效率矢量,过渡操纵面矢量的偏转角矢量通过n行n列的转换矩阵T转化为真实操纵面矢量U的偏转角矢量δ,即且过渡操纵面矢量控制效率矩阵转换矩阵T的第a行b列表示为T(a,b),对于所有的j=1,2,…,n1-1,n1,T(pj,j)=1;对于所有的j=1,2,…,n2-1,n2,T(qj,j+n1)=1;对于所有的j=1,2,…,n3-1,n3,T(rj,j+n1+n2)=1;矩阵T的所有其他元素都为0;
②由加权均衡分配矩阵计算虚拟操纵面的非主导控制效率,得到虚拟操纵面矢量的控制效率矩阵;
其中,分别表示滚转操纵面(i=1,2,…,n1-1,n1)的滚转控制效率、偏航控制效率及俯仰控制效率,分别表示偏航操纵面(j=1,2,…,n2-1,n2)的滚转控制效率、偏航控制效率及俯仰控制效率,分别表示俯仰操纵面(t=1,2,…,n3-1,n3)的滚转控制效率、偏航控制效率及俯仰控制效率;
虚拟操纵面矢量的控制效率矩阵为:
其中,虚拟副翼的滚转控制效率虚拟方向舵的偏航控制效率虚拟升降舵的俯仰控制效率为主导控制效率,均由步骤五得出。进而确定虚拟操纵面的非主导控制效率:虚拟副翼的偏航控制效率 (i=1,2,…,n1-1,n1)表示虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数;虚拟副翼的俯仰控制效率 (i=1,2,…,n1-1,n1)表示虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数;虚拟方向舵的滚转控制效率 (i=1,2,…,n2-1,n2)表示虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数;虚拟方向舵的俯仰控制效率 (i=1,2,…,n2-1,n2)表示虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数;虚拟升降舵的滚转控制效率 (i=1,2,…,n3-1,n3)表示虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数;虚拟升降舵的偏航控制效率 (i=1,2,…,n3-1,n3)表示虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数。
本发明具有的优点在于:
1、本发明提出的多操纵面飞机的一种均衡操纵分配方法,从操纵面提供控制效率的物理机理出发,将多操纵面的控制等效为特定控制效率及偏转范围限制的虚拟副翼、虚拟方向舵、虚拟升降舵的控制,物理概念清晰,简化了控制律设计。
2、本发明提出的多操纵面飞机的一种均衡操纵分配方法,由控制律计算得到虚拟操纵面偏转角,并通过加权均衡分配矩阵运算得到真实操纵面的偏转角,计算简便,便于工程实现,也避免了操纵面角速率超限的问题。
3、本发明提出的多操纵面飞机的一种均衡操纵分配方法,对操纵面引入标称化的概念,将操纵面的控制效率与操纵面的饱和度联系起来,从而为控制效率的均衡分配提供指导。
4、本发明提出的多操纵面飞机的一种均衡操纵分配方法,根据标称化控制效率及所设计的主、副权限系数对于主、副、辅三级分配控制效率,综合考虑了操纵面的偏转范围、控制效率、控制级别、权限系数,并充分利用了有效操纵面。
5、本发明提出的多操纵面飞机的一种均衡操纵分配方法,采用均衡分配的方法将陈化主、副、辅操纵面期望标称化主导控制效率分配给各个标称化操纵面,防止因某些操纵面分担过多的控制效率而引起操纵面饱和。
6、本发明提出的多操纵面飞机的一种均衡操纵分配方法,加权均衡分配矩阵的计算简便,概念清晰。
附图说明
图1本发明中虚拟操纵面与真实操纵面的等效控制图;
图2本发明中真实操纵面的分类、标称化及分级过程图;
图3本发明中滚转操纵面的分配系数矢量的计算过程图。
具体实施方式
下面结合附图对本发明进一步详细说明:
本发明提出的多操纵面飞机的一种均衡操纵分配方法,具体包括以下几个步骤:
步骤一:将多操纵面的控制分配问题转化为虚拟操纵面与真实操纵面的等效控制问题。
设某飞机的n个真实操纵面组成真实操纵面矢量U=[u1 u2 …ui …un-1 un]′,其中ui表示第i个真实操纵面,真实操纵面偏转角矢量为δ=[δ1 δ2 …δi …δn-1 δn]′,其中δi表示第i个真实操纵面的偏转角,且控制效率矩阵为 其中,(i=1,2,…,n-1,n)表示第i个真实操纵面ui的控制效率矢量,且 分别表示第i个真实操纵面ui的滚转控制效率、偏航控制效率及俯仰控制效率,且第i个真实操纵面ui(i=1,2,…,n-1,n)的偏转角δi的偏转范围为[-Ri,Ri],角速率限制为[-ρi,ρi]。
传统的控制分配问题往往是给出期望力矩矢量 其中分别表示期望的滚转力矩、期望的偏航力矩、期望的俯仰力矩,求解真实操纵面偏转角矢量δ,使得满足并满足约束-Ri≤δi≤Ri,(i=1,2,…,n-1,n)。要求控制律输出期望力矩矢量指令,增加了控制律设计的复杂性,并且由期望力矩解算真实操纵面偏转角矢量往往采用复杂的优化过程,计算量大,且需要考虑可能出现的偏转角速率过大的问题。
本发明将多操纵面的控制分配问题转化为虚拟操纵面与真实操纵面的等效控制问题。建立虚拟操纵面矢量 其中,表示虚拟副翼,且虚拟副翼偏转范围为[-Rk Rk];表示虚拟方向舵,且虚拟方向舵偏转范围为[-Rl Rl];表示虚拟升降舵,且虚拟升降舵偏转范围为[-Rm Rm]。虚拟操纵面偏转角矢量 其中 分别表示虚拟副翼偏转角、虚拟方向舵偏转角、虚拟升降舵偏转角;虚拟操纵面矢量的控制效率矩阵定义为:
其中,为虚拟副翼的滚转控制效率,为虚拟副翼的偏航控制效率,为虚拟副翼的俯仰控制效率,为虚拟方向舵的滚转控制效率,为虚拟方向舵的偏航控制效率,为虚拟方向舵的俯仰控制效率,为虚拟升降舵的滚转控制效率,为虚拟升降舵的偏航控制效率,为虚拟升降舵的俯仰控制效率。
控制分配目的是综合考虑操纵面偏转范围、控制效率、控制级别、权限系数来设计加权均衡分配矩阵K,使得真实操纵面的偏转角矢量δ与虚拟操纵面偏转角矢量满足关系且虚拟操纵面矢量控制的效率矩阵满足图1表明了虚拟操纵面与真实操纵面等效控制的物理意义,控制器根据驾驶员的驾驶指令或者是自动驾驶仪期望的轨迹、姿态指令Cmd及飞机当前状态矢量X计算虚拟操纵面偏转角矢量 从真实操纵面矢量U来看,根据加权均衡矩阵K计算真实操纵面偏转角矢量真实操纵面偏转角矢量根据真实操纵面矢量的控制效率矩阵B,计算得到期望的力矩矢量其中 分别表示期望的滚转力矩、期望的偏航力矩及期望的俯仰力矩。从虚拟操纵面矢量来看,图1中的虚线框内的加权均衡矩阵K与真实操纵面矢量的控制效率矩阵B组合成虚拟操纵面矢量的控制效率矩阵因而,虚拟操纵面矢量根据虚线框的虚拟操纵面矢量的控制效率矩阵计算得到力矩矢量为即虚拟操纵面矢量产生与真实操纵面矢量相同的力矩矢量 从而将控制分配问题转化为虚拟操纵面与真实操纵面的等效控制问题。
步骤二:确定对真实操纵面进行标称化及控制效率加权均衡分配的基本方法。
定义真实操纵面ui(i=1,2,…,n-1,n)的偏转角δi时对应的真实操纵面的饱和度为Si=fabs(δi/Ri)。其中,函数fabs(x)表示小数x的绝对值。真实操纵面的控制效率矢量体现了单位偏转角所产生的滚转控制效率、偏航控制效率及俯仰控制效率的大小,而不同真实操纵面的偏转范围不同,因此控制效率不能体现该真实操纵面提供控制效率所产生的饱和度,需要将真实操纵面的控制效率按照其偏转范围进行统一。
将控制效率矢量为偏转范围为[-Ri Ri]的真实操纵面ui(i=1,2,…,n-1,n)对于偏转范围[-Rj Rj]进行标称化,得到标称化滚转控制效率的绝对值、标称化偏航控制效率的绝对值、标称化俯仰控制效率的绝对值分别为 其中函数fabs(x)表示小数x的绝对值。
设n个真实操纵面中有m个真实操纵面是属于A转动方向(代表滚转转动方向、偏航转动方向或者俯仰转动方向)的操纵面,其中tj∈[1,2,…,n-1,n],且j=1,2,…,m-1,m。真实操纵面(j=1,2,…,m-1,m)偏转角偏转范围分别为控制效率矢量为且在A转动方向的控制效率为设某虚拟操纵面与m个真实操纵面在A转动方向的控制等效,虚拟操纵面的偏转角为虚拟操纵面在A转动方向的控制效率为偏转角的偏转范围为[-Rh Rh]。由于不同的真实操纵面在A转动方向的控制效率可能与虚拟操纵面在A转动方向的控制效率同号或者反号,因此,真实操纵面(j=1,2,…,m-1,m)对于偏转范围[-Rh Rh]进行标称化,需要考虑真实操纵面与虚拟操纵面在A转动方向的控制效率符号的一致性,即标称化操纵面的标称化控制效率矢量:
在A转动方向的标称化控制效率:
标称化偏转角:
控制效率均衡分配是将真实操纵面(j=1,2,…,m-1,m)对于偏转范围[-Rh Rh]进行标称化得到标称化操纵面对应标称化控制效率为要求标称化操纵面按照标称化操纵面的控制效率来分担虚拟操纵面在A转动方向的控制效率即标称化操纵面(j=1,2,…,m-1,m)在A转动方向期望分担的控制效率:
可见,采用均衡分配的方法保证了所有真实操纵面按照标称化操纵面的控制效率来分担虚拟操纵面的控制效率,真实操纵面在负担期望的控制效率时对应的标称化偏转角相同,避免了某些真实操纵面分担过大的控制效率而出现饱和。
步骤三:根据真实操纵面的滚转、偏航及俯仰控制效率的大小,将所有真实操纵面划分为滚转操纵面、偏航操纵面与俯仰操纵面三类。
其中,滚转操纵面的数目计为n1,偏航操纵面的数目计为n2,俯仰操纵面的数目计为n3。则有n1+n2+n3=n。图2给出了将真实操纵面划分为滚转操纵面、偏航操纵面及俯仰操纵面的分类结构。
步骤四:分别将滚转操纵面、偏航操纵面、俯仰操纵面对于所对应的虚拟操纵面的偏转范围进行标称化,并按照标称化控制效率的绝对值排序。
确定虚拟副翼的滚转控制效率与所有滚转操纵面中标称化滚转控制效率的绝对值最大的操纵面的滚转控制效率同号。将所有滚转操纵面对于虚拟副翼的偏转范围[-Rk Rk]进行标称化,并按照标称化滚转控制效率的绝对值由大到小排序,假设排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面为(pj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的滚转操纵面的控制效率矢量及滚转控制效率,表示虚拟副翼的滚转控制效率,滚转操纵面的偏转范围为 且 分别表示标称化滚转操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率。其中,标称化滚转操纵面与真实操纵面对应。则可以建立标称化滚转操纵面矢量 其中(pj∈[1,2,…,n-1,n])表示按照标称化滚转控制效率的绝对值排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面,标称化滚转操纵面矢量的控制效率矩阵 其中表示按照标称化滚转控制效率的绝对值排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面(pj∈[1,2,…,n-1,n])的标称化控制效率矢量。
确定虚拟方向舵的偏航控制效率与所有偏航操纵面中标称化偏航控制效率的绝对值最大的操纵面的偏航控制效率同号。将所有偏航操纵面对于虚拟方向舵的偏转范围[-Rl Rl]进行标称化,并按照标称化偏航控制效率的绝对值由大到小排序,假设排序中第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面为(qj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的偏航操纵面的控制效率矢量及偏航控制效率,表示虚拟方向舵的偏航控制效率,偏航操纵面的偏转范围为 且 分别表示标称化偏航操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率。其中,标称化偏航操纵面与真实操纵面对应。则可以建立标称化偏航操纵面矢量 其中(qj∈[1,2,…,n-1,n])表示按照标称化偏航控制效率的绝对值排序中第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面,标称化偏航操纵面矢量的控制效率矩阵 其中表示按照标称化偏航控制效率的绝对值排序中第j(j=1,2,…n2-1,n2)个标称化偏航操纵面(qj∈[1,2,…,n-1,n])的标称化控制效率矢量。
确定虚拟升降舵的俯仰控制效率与所有俯仰操纵面中标称化俯仰控制效率的绝对值最大的操纵面的俯仰控制效率同号。将所有俯仰操纵面对于虚拟升降舵的偏转范围[-Rm Rm]进行标称化,并按照标称化俯仰控制效率的绝对值由大到小排序,假设排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面为(rj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的俯仰操纵面的控制效率矢量及俯仰控制效率,表示虚拟升降舵的俯仰控制效率,俯仰操纵面的偏转范围为 且 分别表示标称化俯仰操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率。其中,标称化俯仰操纵面与真实操纵面对应。则标称化俯仰操纵面矢量
表示按照标称化俯仰控制效率的绝对值排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面,标称化俯仰操纵面矢量的控制效率矩阵为 其中,表示按照标称化俯仰控制效率的绝对值排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面(rj∈[1,2,…,n-1,n])的标称化控制效率矢量。
步骤五:确定虚拟操纵面的主导控制效率。
虚拟副翼的主导控制效率为滚转控制效率且 为虚拟副翼的滚转控制效率比例系数,考虑到飞行控制系统的稳定性,一般选取在0.8~2.0之间,表示标称化滚转操纵面矢量中第j(j=1,2,…,n1-1,n1)的标称化滚转操纵面的标称化滚转控制效率。虚拟副翼的滚转控制效率等于所有的标称化滚转操纵面期望分担的滚转控制效率。虚拟方向舵的主导控制效率为偏航控制效率且 为虚拟方向舵的偏航控制效率比例系数,考虑到飞行控制系统的稳定性,一般选取在0.8~2.0之间,表示标称化偏航操纵面矢量中第j(j=1,2,…,n2-1,n2)的标称化偏航操纵面的标称化偏航控制效率。虚拟方向舵的偏航控制效率等于所有的标称化偏航操纵面期望分担的偏航控制效率。
虚拟升降舵的主导控制效率为俯仰控制效率且 为虚拟升降舵的俯仰控制效率比例系数,考虑到飞行控制系统的稳定性,一般选取在0.8~2.0之间,表示标称化俯仰操纵面矢量中第j(j=1,2,…,n3-1,n3)的标称化俯仰操纵面的标称化俯仰控制效率。虚拟升降舵的俯仰控制效率等于所有的标称化俯仰操纵面期望分担的俯仰控制效率。
步骤六:根据标称化控制效率将标称化操纵面划分成主、副、辅三级,设定标称化主操纵面与标称化副操纵面的权限系数,并计算主、副、辅各级标称化操纵面期望分担的主导控制效率。
(1)分别将标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面的划分为主、副、辅三级。
确定标称化滚转副操纵面序号临界值nx2满足nx1<nx2≤n1,nx2是整数且满足:
将位于标称化滚转操纵面矢量中前nx1个位置的操纵面划分为标称化滚转主操纵面;将位于标称化滚转操纵面矢量中从(nx1+1)到nx2个位置的操纵面划分为标称化滚转副操纵面;将位于标称化滚转操纵面矢量中从(nx2+1)到n1个位置的操纵面划分为标称化滚转辅操纵面。
则标称化偏航副操纵面序号临界值ny2满足ny1<ny2≤n2,ny2是整数且满足:
将位于标称化偏航操纵面矢量中前ny1个位置的操纵面划分为标称化偏航主操纵面;将位于标称化偏航操纵面矢量中从(ny1+1)到ny2个位置的操纵面划分为标称化偏航副操纵面;将位于标称化偏航操纵面矢量中从(ny2+1)到n2个位置的操纵面划分为标称化偏航辅操纵面。
则标称化俯仰副操纵面序号临界值nz2满足nz1<nz2≤n3,nz2是整数且满足:
将位于标称化俯仰操纵面矢量中前nz1个位置的操纵面划分为标称化俯仰主操纵面;将位于标称化俯仰操纵面矢量中从(nz1+1)到nz2个位置的操纵面划分为标称化俯仰副操纵面;将位于标称化俯仰操纵面矢量中从(nz2+1)到n3个位置的操纵面划分为标称化俯仰辅操纵面。
图2给出了将滚转操纵面、偏航操纵面及俯仰操纵面进行标称化得到标称化滚转操纵面、标称化偏航操纵面及标称化俯仰操纵面,并分别将标称化滚转操纵面、标称化偏航操纵面及标称化俯仰操纵面划分为主、副、辅三级的划分结构图。
(2)分别对标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面设定主、副权限系数,并计算标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面的主、副、辅各级标称化操纵面期望分担的主导控制效率。
虚拟副翼的滚转控制效率由所有的标称化滚转操纵面来分担。按照控制效率均衡分配的方法得到所有的标称化滚转主操纵面所分担的虚拟副翼的滚转控制效率的比率Cxm为:
所有的标称化滚转副操纵面所分担的虚拟副翼的滚转控制效率的比率Cxs为:
所有的标称化滚转辅操纵面所分担的虚拟副翼的滚转控制效率的比率Cxf为:
Cxf=1-Cxm-Cxs
图3给出了虚拟副翼的滚转控制效率的分配过程,根据飞行控制的要求及飞机的配置情况,设定标称化滚转主操纵面的权限系数为Qxm及标称化滚转副操纵面的权限系数为Qxs。一般选取Qxm在0.8~1.5之间,Qxs在0.5~2.0之间。则所有的标称化滚转主操纵面期望分担的滚转控制效率则所有的标称化滚转副操纵面期望分担的滚转控制效率所有的标称化滚转辅操纵面期望分担的滚转控制效率
虚拟方向舵的偏航控制效率由所有的标称化偏航操纵面来分担。按照控制效率均衡分配的方法,得到所有的标称化偏航主操纵面所分担的虚拟方向舵的偏航控制效率的比率Cym为:
所有的标称化偏航副操纵面所分担的虚拟方向舵的偏航控制效率的比率Cys为:
所有的标称化偏航辅操纵面所分担的虚拟方向舵的偏航控制效率的比率Cyf为:
Cyf=1-Cym-Cys
根据飞行控制的要求及飞机的配置情况,设计标称化偏航主操纵面的权限系数为Qym及标称化偏航副操纵面的权限系数为Qys,一般选取Qym在0.8~1.5之间,Qys在0.5~2.0之间。则所有的标称化偏航主操纵面期望分担的偏航控制效率所有的标称化偏航副操纵面期望分担的偏航控制效率所有的标称化偏航辅操纵面期望分担的偏航控制效率
虚拟升降舵的俯仰控制效率由所有的标称化俯仰操纵面来分担。按照控制效率均衡分配的方法,得到所有的标称化俯仰主操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czm为:
所有的标称化俯仰副操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czs为:
所有的标称化俯仰辅操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czf为:
Czf=1-Czm-Czs
根据飞行控制的要求及飞机的配置情况,设计标称化俯仰主操纵面的权限系数为Qzm及标称化俯仰副操纵面的权限系数Qzs,一般选取Qzm在0.8~1.5之间,Qzs在0.5~2.0之间。则所有的标称化俯仰主操纵面期望分担的俯仰控制效率所有的标称化俯仰副操纵面期望分担的俯仰控制效率所有的标称化俯仰辅操纵面期望分担的俯仰控制效率
步骤七:采用控制效率均衡分配的方法计算每个标称化操纵面期望分担的主导控制效率。
(1)计算每个标称化滚转操纵面期望分担的滚转控制效率:
如图3所示,采用控制效率均衡分配的方法将所有的标称化滚转主操纵面期望分担的滚转控制效率分配给每个标称化滚转主操纵面,则第j(j=1,2,…,nx1-1,nx1)个标称化滚转主操纵面期望分担的滚转控制效率为:
如图3所示,采用控制效率均衡分配的方法将所有的标称化滚转副操纵面期望分担的滚转控制效率分配给每个标称化滚转副操纵面,则第j(j=nx1+1,nx1+2,…,nx2-1,nx2)个标称化滚转副操纵面期望分担的滚转控制效率为:
如图3所示,采用控制效率均衡分配的方法将所有的标称化滚转辅操纵面期望分担的滚转控制效率分配每个标称化滚转辅操纵面,则第j(j=nx2+1,nx2+2,…,n1-1,n1)个标称化滚转辅操纵面期望分担的滚转控制效率为:
(2)计算每个标称化偏航操纵面期望分担的偏航控制效率:
采用控制效率均衡分配的方法将所有的标称化偏航副操纵面期望分担的偏航控制效率分配给每个标称化偏航副操纵面,则第j(j=ny1+1,ny1+2,…,ny2-1,ny2)个标称化偏航副操纵面期望分担的偏航控制效率为:
采用控制效率均衡分配的方法将所有的陈化偏航辅操纵面期望分担的偏航控制效率分配给每个标称化偏航辅操纵面,则第j(j=ny2+1,ny2+2,…,n2-1,n2)个标称化偏航辅操纵面期望分担的偏航控制效率为:
(3)计算每个标称化俯仰操纵面期望分担的俯仰控制效率:
采用控制效率均衡分配的方法将所有的标称化俯仰副操纵面期望分担的俯仰控制效率分配给每个标称化俯仰副操纵面,则第j(j=nz1+1,nz1+2,…,nz2-1,nz2)个标称化俯仰副操纵面期望分担的俯仰控制效率为:
采用控制效率均衡分配的方法将所有的标称化俯仰辅操纵面期望分担的俯仰控制效率分配给每个标称化俯仰辅操纵面,则第j(j=nz2+1,nz2+2,…,n3-1,n3)个标称化俯仰辅操纵面期望分担的俯仰控制效率为:
步骤八:计算所有真实操纵面的分配系数,建立加权均衡分配矩阵,求解虚拟操纵面矢量的控制效率矩阵中的非主导控制效率。
(1)计算所有真实操纵面的分配系数:
①计算每个滚转操纵面的分配系数:
由于第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面与标称化之前的滚转操纵面对应,标称化滚转操纵面期望分担的滚转控制效率为滚转操纵面的偏转范围为滚转控制效率为则虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数为:
从而得到滚转操纵面的分配系数矢量 其中(j=1,2,…,n1-1,n1)表示虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数。图3给出了根据虚拟副翼的滚转控制效率比例系数按照控制效率均衡分配的方法得到所有的标称化滚转主操纵面所分担的虚拟副翼的滚转控制效率的比率Cxm、按照控制效率均衡分配的方法得到的所有的标称化滚转副操纵面所分担的虚拟副翼的滚转控制效率的比率Cxs、标称化滚转主操纵面的权限系数Qxm、标称化滚转副操纵面的权限系数Qxs计算滚转操纵面的分配系数矢量的过程。
②计算每个偏航操纵面的分配系数:
由于第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面与标称化之前的偏航操纵面对应,标称化偏航操纵面期望分担的偏航控制效率为偏航操纵面的偏转范围为 偏航控制效率为则虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数为:
③计算每个俯仰操纵面的分配系数:
由于第j(j=1,2,…,n3-1,n3)个标称化俯仰操纵面与标称化之前的俯仰操纵面对应,标称化俯仰操纵面期望分担的俯仰控制效率为俯仰操纵面的偏转范围为俯仰控制效率为则虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数为:
(2)建立加权均衡分配矩阵,求解虚拟操纵面矢量的控制效率矩阵中的非主导控制效率,从而确定虚拟操纵面矢量的控制效率矩阵。
①虚拟操纵面矢量到真实操纵面矢量的加权均衡分配矩阵:
建立过渡操纵面矢量 过渡操纵面的偏转角矢量 过渡操纵面矢量的控制效率矩阵
其中,(i=1,2,…,n1-1,n1)分别表示滚转操纵面的偏转角及控制效率矢量,(j=1,2,…,n2-1,n2)分别表示偏航操纵面的偏转角及控制效率矢量,(t=1,2,…,n3-1,n3)分别表示滚转操纵面的偏转角及控制效率矢量。虚拟操纵面矢量 对应的虚拟操纵面偏转角矢量 虚拟操纵面矢量的控制效率矩阵 其中,表示虚拟副翼的偏转角及控制效率矢量, 表示虚拟方向舵的偏转角及控制效率矢量,分别表示虚拟升降舵的偏转角及控制效率矢量。则虚拟操纵面偏转角矢量到过渡操纵面偏转角矢量的传递关系为其中,表示虚拟操纵面矢量到过渡操纵面矢量的加权均衡分配矩阵:
对应的真实操纵面的偏转角矢量δ=[δ1 δ2 …δi …δn-1 δn]′,真实操纵面的控制效率矩阵 其中δi、分别表示真实操纵面ui的偏转角、控制效率矢量。过渡操纵面矢量的偏转角矢量可以通过n行n列的转换矩阵T转化为真实操纵面矢量U的偏转角矢量δ,即且过渡操纵面矢量控制效率矩阵转换矩阵T的第a行b列表示为T(a,b),对于所有的j=1,2,…,n1-1,n1,T(pj,j)=1;对于所有的j=1,2,…,n2-1,n2,T(qj,j+n1)=1;对于所有的j=1,2,…,n3-1,n3,T(rj,j+n1+n2)=1;矩阵T的所有其他元素都为0。
可见,通过求解加权均衡分配矩阵K实现了多操纵面的控制分配,该计算过程简单,物理概念清晰,便于工程实现。
②由加权均衡分配矩阵计算虚拟操纵面的非主导控制效率,得到虚拟操纵面矢量的控制效率矩阵。
其中,分别表示滚转操纵面(i=1,2,…,n1-1,n1)的滚转控制效率、偏航控制效率及俯仰控制效率,分别表示偏航操纵面(j=1,2,…,n2-1,n2)的滚转控制效率、偏航控制效率及俯仰控制效率,分别表示俯仰操纵面(t=1,2,…,n3-1,n3)的滚转控制效率、偏航控制效率及俯仰控制效率。
虚拟操纵面矢量的控制效率矩阵为:
其中,虚拟副翼的滚转控制效率虚拟方向舵的偏航控制效率虚拟升降舵的俯仰控制效率为主导控制效率,均由步骤五确定。可以确定虚拟操纵面的非主导控制效率:虚拟副翼的偏航控制效率 (i=1,2,…,n1-1,n1)表示虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数;虚拟副翼的俯仰控制效率 (i=1,2,…,n1-1,n1)表示虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数;虚拟方向舵的滚转控制效率 (i=1,2,…,n2-1,n2)表示虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数;虚拟方向舵的俯仰控制效率 (i=1,2,…,n2-1,n2)表示虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数;虚拟升降舵的滚转控制效率 (i=1,2,…,n3-1,n3)表示虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数;虚拟升降舵的偏航控制效率 (i=1,2,…,n3-1,n3)表示虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数。
根据虚拟操纵面矢量的控制效率矩阵来设计控制律,且控制律输出平稳变化的虚拟操纵面偏转角指令,避免了控制律输出期望力矩、复杂优化及偏转角速率过大的问题。
Claims (7)
1.多操纵面飞机的一种均衡操纵分配方法,其特征在于:包括以下几个步骤:
步骤一:将多操纵面的控制分配问题转化为虚拟操纵面与真实操纵面的等效控制问题:
飞机的n个真实操纵面组成真实操纵面矢量U=[u1 u2 …ui …un-1 un]′,其中ui表示第i个真实操纵面,真实操纵面偏转角矢量为δ=[δ1 δ2 …δi …δn-1 δn]′,其中δi表示第i个真实操纵面的偏转角,且控制效率矩阵为 其中,(i=1,2,…,n-1,n)表示第i个真实操纵面ui的控制效率矢量,且 分别表示第i个真实操纵面ui的滚转控制效率、偏航控制效率、俯仰控制效率,且第i个真实操纵面ui(i=1,2,…,n-1,n)的偏转角δi的偏转范围为[-Ri,Ri],角速率限制为[-ρi,ρi];
建立虚拟操纵面矢量 其中,表示虚拟副翼,虚拟副翼偏转范围为[-Rk Rk];表示虚拟方向舵,虚拟方向舵偏转范围为[-Rl Rl];表示虚拟升降舵,虚拟升降舵偏转范围为[-Rm Rm];虚拟操纵面偏转角矢量 其中分别表示虚拟副翼偏转角、虚拟方向舵偏转角、虚拟升降舵偏转角;虚拟操纵面矢量的控制效率矩阵为:
其中,为虚拟副翼的滚转控制效率,为虚拟副翼的偏航控制效率,为虚拟副翼的俯仰控制效率,为虚拟方向舵的滚转控制效率,为虚拟方向舵的偏航控制效率,为虚拟方向舵的俯仰控制效率,为虚拟升降舵的滚转控制效率,为虚拟升降舵的偏航控制效率,为虚拟升降舵的俯仰控制效率;
真实操纵面的偏转角矢量δ与虚拟操纵面偏转角矢量满足关系K为加权均衡分配矩阵,且虚拟操纵面矢量控制的效率矩阵满足控制器根据驾驶员的驾驶指令或自动驾驶仪期望的轨迹、姿态指令Cmd及飞机当前状态矢量X计算虚拟操纵面偏转角矢量 根据加权均衡矩阵K计算真实操纵面偏转角矢量真实操纵面偏转角矢量根据真实操纵面矢量的控制效率矩阵B,计算得到期望的力矩矢量其中 虚拟操纵面矢量根据虚拟操纵面矢量的控制效率矩阵计算得到力矩矢量为
虚拟操纵面矢量产生与真实操纵面矢量相同的力矩矢量 其中 分别表示期望的滚转力矩、期望的偏航力矩、期望的俯仰力矩;
步骤二:确定对真实操纵面进行标称化及控制效率加权均衡分配的基本方法:
定义真实操纵面ui(i=1,2,…,n-1,n)的偏转角δi时对应的真实操纵面的饱和度Si为Si=fabs(δi/Ri),其中函数fabs(x)表示小数x的绝对值;
将控制效率矢量为偏转范围为[-Ri Ri]的真实操纵面ui(i=1,2,…,n-1,n)对于偏转范围[-Rj Rj]进行标称化,得到标称化滚转控制效率的绝对值、标称化偏航控制效率的绝对值、标称化俯仰控制效率矢量为的绝对值分别为 其中函数fabs(x)表示小数x的绝对值;
设n个真实操纵面中有m个真实操纵面是属于A转动方向的操纵面,其中tj∈[1,2,…,n-1,n],且真实操纵面(j=1,2,…,n-1,m)偏转角偏转范围分别为 控制效率矢量为且在A转动方向的控制效率为设某虚拟操纵面与m个真实操纵面在A转动方向的控制等效,虚拟操纵面的偏转角为虚拟操纵面在A转动方向的控制效率为偏转角的偏转范围为[-Rh Rh];真实操纵面(j=1,2,…,m-1,m)对于偏转范围[-Rh Rh]进行标称化,标称化操纵面的标称化控制效率矢量:
在A转动方向的标称化控制效率:
标称化偏转角:
步骤三:根据真实操纵面的滚转、偏航及俯仰控制效率的大小,将所有真实操纵面划分为滚转操纵面、偏航操纵面与俯仰操纵面:
滚转操纵面的数目计为n1,偏航操纵面的数目计为n2,俯仰操纵面的数目计为n3,n1+n2+n3=n;
步骤四:分别将滚转操纵面、偏航操纵面、俯仰操纵面对于所对应的虚拟操纵面的偏转范围进行标称化,并按照标称化控制效率的绝对值排序:
确定虚拟副翼的滚转控制效率与所有滚转操纵面中标称化滚转控制效率的绝对值最大的操纵面的滚转控制效率同号,将所有滚转操纵面对于虚拟副翼的偏转范围[-Rk Rk]进行标称化,并按照标称化滚转控制效率由大到小排序,排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面为(pj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的滚转操纵面的控制效率矢量及滚转控制效率,表示虚拟副翼的滚转控制效率,滚转操纵面的偏转范围为 且 分别表示标称化滚转操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率,建立标称化滚转操纵面矢量 其中(pj∈[1,2,…,n-1,n])表示按照标称化滚转控制效率的绝对值排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面,标称化滚转操纵面矢量的控制效率矩阵 其中表示按照标称化滚转控制效率的绝对值排序中第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面(pj∈[1,2,…,n-1,n])的标称化控制效率矢量;
确定虚拟方向舵的偏航控制效率与所有偏航操纵面中标称化偏航控制效率的绝对值最大的操纵面的偏航控制效率同号,将所有偏航操纵面对于虚拟方向舵的偏转范围[-Rl Rl]进行标称化,并按照标称化偏航控制效率的绝对值由大到小排序,排序中第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面为(qj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的偏航操纵面的控制效率矢量及偏航控制效率,表示虚拟方向舵的偏航控制效率,偏航操纵面的偏转范围为 且 分别表示标称化偏航操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率,建立标称化偏航操纵面矢量 其中(qj∈[1,2,…,n-1,n])表示按照标称化偏航控制效率的绝对值排序中第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面,标称化偏航操纵面矢量的控制效率矩阵 其中表示按照标称化偏航控制效率的绝对值排序中第j(j=1,2,…n2-1,n2)个标称化偏航操纵面(qj∈[1,2,…,n-1,n])的标称化控制效率矢量;
确定虚拟升降舵的俯仰控制效率与所有俯仰操纵面中标称化俯仰控制效率的绝对值最大的操纵面的俯仰控制效率同号,将所有俯仰操纵面对于虚拟升降舵的偏转范围[-Rm Rm]进行标称化,并按照标称化俯仰控制效率的绝对值由大到小排序,排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面为(rj∈[1,2,…,n-1,n]),其标称化控制效率矢量为:
分别表示标称化之前的俯仰操纵面的控制效率矢量及俯仰控制效率,表示虚拟升降舵的俯仰控制效率,俯仰操纵面的偏转范围为 且 分别表示标称化俯仰操纵面的标称化滚转控制效率、标称化偏航控制效率及标称化俯仰控制效率;则标称化俯仰操纵面矢量 其中(rj∈[1,2,…,n-1,n])表示按照标称化俯仰控制效率的绝对值排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面,标称化俯仰操纵面矢量的控制效率矩阵为 其中,表示按照标称化俯仰控制效率的绝对值排序中第j(j=1,2,…n3-1,n3)个标称化俯仰操纵面(rj∈[1,2,…,n-1,n])的标称化控制效率矢量;
步骤五:确定虚拟操纵面的主导控制效率:
步骤六:根据标称化控制效率将标称化操纵面划分成主、副、辅三级,设定标称化主操纵面与标称化副操纵面的权限系数,并计算主、副、辅各级标称化操纵面期望分担的主导控制效率:
(1)分别将标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面的划分为主、副、辅三级:
确定标称化滚转副操纵面序号临界值nx2满足nx1<nx2≤n1,nx2是整数且满足:
将位于标称化滚转操纵面矢量中前nx1个位置的操纵面划分为标称化滚转主操纵面;将位于标称化滚转操纵面矢量中从(nx1+1)到nx2个位置的操纵面划分为标称化滚转副操纵面;将位于标称化滚转操纵面矢量中从(nx2+1)到n1个位置的操纵面划分为标称化滚转辅操纵面;
则标称化偏航副操纵面序号临界值ny2满足ny1<ny2≤n2,ny2是整数且满足:
将位于标称化偏航操纵面矢量中前ny1个位置的操纵面划分为标称化偏航主操纵面;将位于标称化偏航操纵面矢量中从(ny1+1)到ny2个位置的操纵面划分为标称化偏航副操纵面;将位于标称化偏航操纵面矢量中从(ny2+1)到n2个位置的操纵面划分为标称化偏航辅操纵面;
则标称化俯仰副操纵面序号临界值nz2满足nz1<nz2≤n3,nz2是整数且满足:
将位于标称化俯仰操纵面矢量中前nz1个位置的操纵面划分为标称化俯仰主操纵面;将位于标称化俯仰操纵面矢量中从(nz1+1)到nz2个位置的操纵面划分为标称化俯仰副操纵面;将位于标称化俯仰操纵面矢量中从(nz2+1)到n3个位置的操纵面划分为标称化俯仰辅操纵面;
(2)分别对标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面设定主、副权限系数,并计算标称化滚转操纵面、标称化偏航操纵面、标称化俯仰操纵面的主、副、辅各级标称化操纵面期望分担的主导控制效率:
按照控制效率均衡分配的方法得到所有的标称化滚转主操纵面所分担的虚拟副翼的滚转控制效率的比率Cxm为:
所有的标称化滚转副操纵面所分担的虚拟副翼的滚转控制效率的比率Cxs为:
所有的标称化滚转辅操纵面所分担的虚拟副翼的滚转控制效率的比率Cxf为:
Cxf=1-Cxm-Cxs
标称化滚转主操纵面的权限系数为Qxm及标称化滚转副操纵面的权限系数为Qxs,则所有的标称化滚转主操纵面期望分担的滚转控制效率则所有的标称化滚转副操纵面期望分担的滚转控制效率所有的标称化滚转辅操纵面期望分担的滚转控制效率
按照控制效率均衡分配的方法,得到所有的标称化偏航主操纵面所分担的虚拟方向舵的偏航控制效率的比率Cym为:
所有的标称化偏航副操纵面所分担的虚拟方向舵的偏航控制效率的比率Cys为:
所有的标称化偏航辅操纵面所分担的虚拟方向舵的偏航控制效率的比率Cyf为:
Cyf=1-Cym-Cys
标称化偏航主操纵面的权限系数为Qym及标称化偏航副操纵面的权限系数为Qys,则所有的标称化偏航主操纵面期望分担的偏航控制效率所有的标称化偏航副操纵面期望分担的偏航控制效率所有的标称化偏航辅操纵面期望分担的偏航控制效率
按照控制效率均衡分配的方法,得到所有的标称化俯仰主操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czm为:
所有的标称化俯仰副操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czs为:
所有的标称化俯仰辅操纵面所分担的虚拟升降舵的俯仰控制效率的比率Czf为:
Czf=1-Czm-Czs
标称化俯仰主操纵面的权限系数为Qzm,标称化俯仰副操纵面的权限系数为Qzs,所有的标称化俯仰主操纵面期望分担的俯仰控制效率所有的标称化俯仰副操纵面期望分担的俯仰控制效率所有的标称化俯仰辅操纵面期望分担的俯仰控制效率
步骤七:采用控制效率均衡分配的方法计算每个标称化操纵面期望分担的主导控制效率:
(1)计算每个标称化滚转操纵面期望分担的滚转控制效率:
采用控制效率均衡分配的方法将所有的标称化滚转副操纵面期望分担的滚转控制效率分配给每个标称化滚转副操纵面,则第j(j=nx1+1,nx1+2,…,nx2-1,nx2)个标称化滚转副操纵面期望分担的滚转控制效率为:
采用控制效率均衡分配的方法将所有的标称化滚转辅操纵面期望分担的滚转控制效率分配每个标称化滚转辅操纵面,则第j(j=nx2+1,nx2+2,…,n1-1,n1)个标称化滚转辅操纵面期望分担的滚转控制效率为:
(2)计算每个标称化偏航操纵面期望分担的偏航控制效率:
采用控制效率均衡分配的方法将所有的标称化偏航副操纵面期望分担的偏航控制效率分配给每个标称化偏航副操纵面,则第j(j=ny1+1,ny1+2,…,ny2-1,ny2)个标称化偏航副操纵面期望分担的偏航控制效率为:
采用控制效率均衡分配的方法将所有的陈化偏航辅操纵面期望分担的偏航控制效率分配给每个标称化偏航辅操纵面,则第j(j=ny2+1,ny2+2,…,n2-1,n2)个标称化偏航辅操纵面期望分担的偏航控制效率为:
(3)计算每个标称化俯仰操纵面期望分担的俯仰控制效率:
采用控制效率均衡分配的方法将所有的标称化俯仰副操纵面期望分担的俯仰控制效率分配给每个标称化俯仰副操纵面,则第j(j=nz1+1,nz1+2,…,nz2-1,nz2)个标称化俯仰副操纵面期望分担的俯仰控制效率为:
采用控制效率均衡分配的方法将所有的标称化俯仰辅操纵面期望分担的俯仰控制效率分配给每个标称化俯仰辅操纵面,则第j(j=nz2+1,nz2+2,…,n3-1,n3)个标称化俯仰辅操纵面期望分担的俯仰控制效率为:
步骤八:计算所有真实操纵面的分配系数,建立加权均衡分配矩阵,求解虚拟操纵面矢量的控制效率矩阵中的非主导控制效率:
(1)计算所有真实操纵面的分配系数:
①计算每个滚转操纵面的分配系数:
由于第j(j=1,2,…,n1-1,n1)个标称化滚转操纵面与标称化之前的滚转操纵面对应,标称化滚转操纵面期望分担的滚转控制效率为滚转操纵面的偏转范围为滚转控制效率为则虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数为:
②计算每个偏航操纵面的分配系数:
由于第j(j=1,2,…,n2-1,n2)个标称化偏航操纵面与标称化之前的偏航操纵面对应,标称化偏航操纵面期望分担的偏航控制效率为偏航操纵面的偏转范围为 偏航控制效率为则虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数为:
③计算每个俯仰操纵面的分配系数:
由于第j(j=1,2,…,n3-1,n3)个标称化俯仰操纵面与标称化之前的俯仰操纵面对应,标称化俯仰操纵面期望分担的俯仰控制效率为俯仰操纵面的偏转范围为 俯仰控制效率为则虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数为:
(2)建立加权均衡分配矩阵,求解虚拟操纵面矢量的控制效率矩阵中的非主导控制效率,从而确定虚拟操纵面矢量的控制效率矩阵:
①虚拟操纵面矢量到真实操纵面矢量的加权均衡分配矩阵:
建立过渡操纵面矢量 过渡操纵面的偏转角矢量 过渡操纵面矢量的控制效率矩阵
其中,(i=1,2,…,n1-1,n1)分别表示滚转操纵面的偏转角及控制效率矢量,(j=1,2,…,n2-1,n2)分别表示偏航操纵面的偏转角及控制效率矢量,(t=1,2,…,n3-1,n3)分别表示滚转操纵面的偏转角及控制效率矢量,虚拟操纵面矢量 对应的虚拟操纵面偏转角矢量 虚拟操纵面矢量的控制效率矩阵 其中,表示虚拟副翼的偏转角及控制效率矢量, 表示虚拟方向舵的偏转角及控制效率矢量,分别表示虚拟升降舵的偏转角及控制效率矢量;则虚拟操纵面偏转角矢量到过渡操纵面偏转角矢量的传递关系为其中,表示虚拟操纵面矢量到过渡操纵面矢量的加权均衡分配矩阵:
虚拟操纵面矢量的控制效率矩阵与过渡操纵面矢量的控制效率矩阵满足真实操纵面矢量U=[u1 u2 …ui …un-1 un]′,对应的真实操纵面的偏转角矢量δ=[δ1 δ2 …δi …δn-1 δn]′,真实操纵面的控制效率矩阵 其中δi、分别表示真实操纵面ui的偏转角、控制效率矢量,过渡操纵面矢量的偏转角矢量通过n行n列的转换矩阵T转化为真实操纵面矢量U的偏转角矢量δ,即且过渡操纵面矢量控制效率矩阵转换矩阵T的第a行b列表示为T(a,b),对于所有的j=1,2,…,n1-1,n1,T(pj,j)=1;对于所有的j=1,2,…,n2-1,n2,T(qj,j+n1)=1;对于所有的j=1,2,…,n3-1,n3,T(rj,j+n1+n2)=1;矩阵T的所有其他元素都为0;
②由加权均衡分配矩阵计算虚拟操纵面的非主导控制效率,得到虚拟操纵面矢量的控制效率矩阵;
其中,分别表示滚转操纵面(i=1,2,…,n1-1,n1)的滚转控制效率、偏航控制效率及俯仰控制效率,分别表示偏航操纵面(j=1,2,…,n2-1,n2)的滚转控制效率、偏航控制效率及俯仰控制效率,分别表示俯仰操纵面(t=1,2,…,n3-1,n3)的滚转控制效率、偏航控制效率及俯仰控制效率;
虚拟操纵面矢量的控制效率矩阵为:
其中,虚拟副翼的滚转控制效率虚拟方向舵的偏航控制效率虚拟升降舵的俯仰控制效率为主导控制效率,均由步骤五得出,进而确定虚拟操纵面的非主导控制效率:虚拟副翼的偏航控制效率 (i=1,2,…,n1-1,n1)表示虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数;虚拟副翼的俯仰控制效率 (i=1,2,…,n1-1,n1)表示虚拟副翼的偏转角到滚转操纵面的偏转角的分配系数;虚拟方向舵的滚转控制效率 (i=1,2,…,n2-1,n2)表示虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数;虚拟方向舵的俯仰控制效率 (i=1,2,…,n2-1,n2)表示虚拟方向舵的偏转角到偏航操纵面的偏转角的分配系数;虚拟升降舵的滚转控制效率 (i=1,2,…,n3-1,n3)表示虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数;虚拟升降舵的偏航控制效率 (i=1,2,…,n3-1,n3)表示虚拟升降舵的偏转角到俯仰操纵面的偏转角的分配系数。
5.根据权利要求1所述的多操纵面飞机的一种均衡操纵分配方法,其特征在于:步骤六中的标称化滚转主操纵面的权限系数为Qxm的取值为0.8~1.5,标称化滚转副操纵面的权限系数Qxs的取值为0.5~2.0之间。
6.根据权利要求1所述的多操纵面飞机的一种均衡操纵分配方法,其特征在于:步骤六中的标称化偏航主操纵面的权限系数为Qym的取值为0.8~1.5,标称化偏航副操纵面的权限系数Qys的取值为0.5~2.0。
7.根据权利要求1所述的多操纵面飞机的一种均衡操纵分配方法,其特征在于:步骤六中的标称化俯仰主操纵面的权限系数Qzm的取值为0.8~1.5,标称化俯仰副操纵面的权限系数Qzs的取值为0.5~2.0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110165732 CN102320378B (zh) | 2011-06-20 | 2011-06-20 | 多操纵面飞机的一种均衡操纵分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110165732 CN102320378B (zh) | 2011-06-20 | 2011-06-20 | 多操纵面飞机的一种均衡操纵分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102320378A true CN102320378A (zh) | 2012-01-18 |
CN102320378B CN102320378B (zh) | 2013-07-24 |
Family
ID=45448264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110165732 Expired - Fee Related CN102320378B (zh) | 2011-06-20 | 2011-06-20 | 多操纵面飞机的一种均衡操纵分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102320378B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104216417A (zh) * | 2014-09-22 | 2014-12-17 | 北京航空航天大学 | 一种基于四矢量螺旋桨的平流层飞艇的动态控制分配方法 |
CN104331084B (zh) * | 2014-09-30 | 2017-05-03 | 中国运载火箭技术研究院 | 一种基于方向舵控滚转策略的气动舵偏范围计算方法 |
CN109270947A (zh) * | 2018-12-13 | 2019-01-25 | 北京航空航天大学 | 倾转旋翼无人机飞行控制系统 |
CN110187634A (zh) * | 2018-02-23 | 2019-08-30 | 北京京东尚科信息技术有限公司 | 飞行器的控制方法、装置和计算机可读存储介质 |
CN110398961A (zh) * | 2019-07-09 | 2019-11-01 | 哈尔滨工程大学 | 一种气垫船航迹保持舵桨协同控制与分配方法 |
CN110411289A (zh) * | 2019-06-13 | 2019-11-05 | 上海航天控制技术研究所 | 一种抑制强机弹干扰的分离稳定控制方法 |
CN112198817A (zh) * | 2020-09-23 | 2021-01-08 | 深圳市领峰电动智能科技有限公司 | 无人机控制方法、装置、设备、无人机和介质 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106598056B (zh) * | 2016-11-23 | 2019-05-17 | 中国人民解放军空军工程大学 | 一种提升固定翼飞机隐身性能的舵面优先级调节方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5791598A (en) * | 1996-01-16 | 1998-08-11 | Globalstar L.P. and Daimler-Benz Aerospace AG | Dynamic bias for orbital yaw steering |
CN101804862A (zh) * | 2010-04-07 | 2010-08-18 | 南京航空航天大学 | 一种无人机推力变向装置及其控制方法 |
CN101833761A (zh) * | 2010-04-20 | 2010-09-15 | 南京航空航天大学 | 基于合作目标特征线的无人机位姿估计方法 |
CN102073755A (zh) * | 2010-11-10 | 2011-05-25 | 南京航空航天大学 | 近空间高超声速飞行器运动控制仿真方法 |
-
2011
- 2011-06-20 CN CN 201110165732 patent/CN102320378B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5791598A (en) * | 1996-01-16 | 1998-08-11 | Globalstar L.P. and Daimler-Benz Aerospace AG | Dynamic bias for orbital yaw steering |
CN101804862A (zh) * | 2010-04-07 | 2010-08-18 | 南京航空航天大学 | 一种无人机推力变向装置及其控制方法 |
CN101833761A (zh) * | 2010-04-20 | 2010-09-15 | 南京航空航天大学 | 基于合作目标特征线的无人机位姿估计方法 |
CN102073755A (zh) * | 2010-11-10 | 2011-05-25 | 南京航空航天大学 | 近空间高超声速飞行器运动控制仿真方法 |
Non-Patent Citations (1)
Title |
---|
余晓艇等: "无人机轮式着陆横侧向控制", 《系统仿真学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104216417A (zh) * | 2014-09-22 | 2014-12-17 | 北京航空航天大学 | 一种基于四矢量螺旋桨的平流层飞艇的动态控制分配方法 |
CN104216417B (zh) * | 2014-09-22 | 2017-09-12 | 北京航空航天大学 | 一种基于四矢量螺旋桨的平流层飞艇的动态控制分配方法 |
CN104331084B (zh) * | 2014-09-30 | 2017-05-03 | 中国运载火箭技术研究院 | 一种基于方向舵控滚转策略的气动舵偏范围计算方法 |
CN110187634A (zh) * | 2018-02-23 | 2019-08-30 | 北京京东尚科信息技术有限公司 | 飞行器的控制方法、装置和计算机可读存储介质 |
CN109270947A (zh) * | 2018-12-13 | 2019-01-25 | 北京航空航天大学 | 倾转旋翼无人机飞行控制系统 |
CN109270947B (zh) * | 2018-12-13 | 2020-07-10 | 北京航空航天大学 | 倾转旋翼无人机飞行控制系统 |
CN110411289A (zh) * | 2019-06-13 | 2019-11-05 | 上海航天控制技术研究所 | 一种抑制强机弹干扰的分离稳定控制方法 |
CN110411289B (zh) * | 2019-06-13 | 2021-10-15 | 上海航天控制技术研究所 | 一种抑制强机弹干扰的分离稳定控制方法 |
CN110398961A (zh) * | 2019-07-09 | 2019-11-01 | 哈尔滨工程大学 | 一种气垫船航迹保持舵桨协同控制与分配方法 |
CN112198817A (zh) * | 2020-09-23 | 2021-01-08 | 深圳市领峰电动智能科技有限公司 | 无人机控制方法、装置、设备、无人机和介质 |
CN112198817B (zh) * | 2020-09-23 | 2022-07-12 | 峰飞航空科技(昆山)有限公司 | 无人机控制方法、装置、设备、无人机和介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102320378B (zh) | 2013-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102320378B (zh) | 多操纵面飞机的一种均衡操纵分配方法 | |
CN103488814B (zh) | 一种适用于再入飞行器姿态控制的闭环仿真系统 | |
Usta et al. | Aircraft roll control system using LQR and fuzzy logic controller | |
CN104635741B (zh) | 可重复使用运载器再入姿态控制方法 | |
CN106874617B (zh) | 一种高效的直升机机动飞行品质等级评估方法 | |
CN104216417A (zh) | 一种基于四矢量螺旋桨的平流层飞艇的动态控制分配方法 | |
CN104331084B (zh) | 一种基于方向舵控滚转策略的气动舵偏范围计算方法 | |
CN110347170A (zh) | 可重复使用运载器再入段鲁棒容错制导系统及工作方法 | |
CN114942649B (zh) | 一种基于反步法的飞机俯仰姿态与航迹角解耦控制方法 | |
Akyazı et al. | A self-tuning fuzzy logic controller for aircraft roll control system | |
CN112487551A (zh) | 一种飞翼无人机多舵面链式直接控制分配与重构方法 | |
CN105173051A (zh) | 一种平流层飞艇的制导控制一体化及控制分配方法 | |
Acheson et al. | Examination of unified control incorporating generalized control allocation | |
CN108919827A (zh) | 一种推力矢量飞行器控制的双次优化快速分配方法 | |
CN102692928B (zh) | 基于飞行器四元数模型的控制器区域设计方法 | |
CN111897219A (zh) | 基于在线逼近器的倾转四旋翼无人机过渡飞行模式最优鲁棒控制方法 | |
Sadien et al. | A detailed comparison of control allocation techniques on a realistic on-ground aircraft benchmark | |
Yu et al. | Predefined-time fault-tolerant attitude control for tailless aircraft considering actuator input saturation | |
Movahhed et al. | Output tracking of a 6-DOF flying wing UAV in longitudinal motion using LQR optimized sliding mode control with integral action | |
CN112009669A (zh) | 一种基于空气舵的飞行器的减速方法及装置 | |
Reinhardt et al. | Nonlinear model predictive attitude control for fixed-wing unmanned aerial vehicle based on a wind frame formulation | |
McIntosh et al. | A Switching-Free Control Architecture for Transition Maneuvers of a Quadrotor Biplane Tailsitter | |
CN114771807A (zh) | 一种飞机多操纵面协调控制管理系统 | |
Zhang et al. | Fault tolerant tracking control for quad-rotor helicopter via robust adaptive technique | |
Timchenko et al. | Robust-optimal control system of quadrocopter for maritime traffic's monitoring |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130724 Termination date: 20140620 |
|
EXPY | Termination of patent right or utility model |