1.一種針對三旋翼無人機舵機堵塞故障的容錯控制方法,其特征是,步驟如下:首先定義慣性坐標(biāo)系{I}、機體坐標(biāo)系{B}和目標(biāo)坐標(biāo)系{Bd},通過分析舵機對三旋翼無人機的作用原理,并考慮外部擾動對其動力學(xué)特性的影響,得到三旋翼無人機執(zhí)行器發(fā)生故障時的非線性動力學(xué)模型:
式(1)中各變量定義如下:ω=[ω1 ω2 ω3]T∈R3×1表示機體坐標(biāo)系{B}相對于慣性坐標(biāo)系{I}的姿態(tài)角速度,ω1,ω2,ω3分別表示滾轉(zhuǎn)角速度、俯仰角速度和偏航角速度,[·]T表示矩陣的轉(zhuǎn)置,∈表示集合間的“屬于”關(guān)系,R3×1表示3行1列的實數(shù)向量,表示求取ω的一階時間導(dǎo)數(shù);J=diag{[J1 J2 J3]T}∈R3×3為轉(zhuǎn)動慣量矩陣,diag{[J1 J2 J3]T}表示向量[J1 J2 J3]張成的對角矩陣,J1,J2,J3分別表示繞各個坐標(biāo)軸的轉(zhuǎn)動慣量;S(ω)表示求取ω對應(yīng)的反對稱矩陣;為一系數(shù)矩陣,其中l(wèi)表示前面某一電機中心到無人機軸心的距離,l3表示舵機中心到無人機軸心的距離,α表示前面兩個電機連線與某一電機和無人機軸心連線之間的夾角,k為電機的升力系數(shù),δf為舵機發(fā)生堵塞時的偏轉(zhuǎn)角度,l,l3,α,k均為已知常數(shù),δf為未知常數(shù),sin(·)和cos(·)分別表示正弦和余弦函數(shù);fδ=[fδ1 fδ2 fδ3]T∈R3×1表示故障發(fā)生后的升力向量,fδ1,fδ2,fδ3分別表示故障發(fā)生后三個電機產(chǎn)生的升力,D=diag{[d1 d2 d3]T}∈R3×3為外部擾動矩陣,diag{[d1 d2 d3]T}表示向量[d1 d2 d3]張成的對角矩陣,d1,d2,d3分別表示作用于各個通道的外部擾動;
定義變量λ1=-l3 cosδf,λ2=k cosδf+l3 sinδf,則λ1,λ2為未知常數(shù),A(δf)可寫為式(1)可表示為
為避免姿態(tài)表示奇異性問題,采用基于單位四元數(shù)的姿態(tài)表示方法,機體坐標(biāo)系{B}在慣性坐標(biāo)系{I}下的表達用“等效軸角坐標(biāo)系”方法,將{B}和{I}重合,將{B}繞矢量k0∈R3×1按右手定則旋轉(zhuǎn)角,得到當(dāng)前姿態(tài)單位四元數(shù)其中且滿足k0∈R3×1為定義在坐標(biāo)系{I}中的任意單位矢量,為坐標(biāo)系{B}繞矢量k旋轉(zhuǎn)的任意角度;由機體坐標(biāo)系{B}到慣性坐標(biāo)系{I}的坐標(biāo)變換矩陣用四元數(shù)表示為I3為3×3的單位矩陣,S(qv)表示求取qv對應(yīng)的反對稱矩陣,同理,目標(biāo)坐標(biāo)系{Bd}在慣性坐標(biāo)系{I}下的表達也用“等效軸角坐標(biāo)系”方法,將{Bd}和{I}重合,將{Bd}繞矢量kd∈R3×1按右手定則旋轉(zhuǎn)角,得到目標(biāo)姿態(tài)單位四元數(shù)其中且滿足kd∈R3×1同樣為定義在坐標(biāo)系{I}中的任意單位矢量,為坐標(biāo)系{Bd}繞矢量kd旋轉(zhuǎn)的任意角度;由目標(biāo)坐標(biāo)系{Bd}到慣性坐標(biāo)系{I}的坐標(biāo)變換矩陣用四元數(shù)表示為S(qvd)表示求取qvd對應(yīng)的反對稱矩陣,為了描述三旋翼無人機當(dāng)前姿態(tài)與目標(biāo)姿態(tài)之間的差異,定義姿態(tài)誤差四元數(shù)
其中e0和ev同樣滿足由目標(biāo)坐標(biāo)系{Bd}到機體坐標(biāo)系{B}的坐標(biāo)變換矩陣示為S(ev)表示求取ev對應(yīng)的反對稱矩陣,定義角速度誤差其中ωd∈R3×1表示目標(biāo)坐標(biāo)系{Bd}相對于慣性坐標(biāo)系{I}的姿態(tài)角速度;
為了對三旋翼無人機舵機堵塞故障進行更有針對性的容錯控制,采用基于自適應(yīng)滑模方法的觀測器技術(shù)對故障進行觀測,設(shè)計觀測器為:
其中表示對ω的估計值,表示求取的一階時間導(dǎo)數(shù),v=[v1 v2 v3]T∈R3×1,表示求取v的一階時間導(dǎo)數(shù),定義ω的估計誤差為分別表示對λ1,λ2的估計值,SIG2=[k21sign(eω1) k22sign(eω2) k23sign(eω3)]T,SIG1=[k11|eω1|1/2 sign(eω1) k12|eω2|1/2 sign(eω2) k13|eω3|1/2 sign(eω3)]T,其中k11,k12,k13,k21,k22,k23均為正常數(shù),eω1,eω2,eω3為eω的三個元素,sign(·)表示符號函數(shù),|·|1/2表示求取絕對值的次方;
定義變量z1i=|eωi|1/2 sign(eωi),z2i=vi,i=1,2,3,設(shè)計控制器為:
其中表示求取矩陣的逆矩陣,表示求取ωd的一階時間導(dǎo)數(shù),G=diag{[g1 g2 g3]T}∈R3×3為正常數(shù)增益矩陣,diag{[g1 g2 g3]T}表示向量[g1 g2 g3]張成的對角矩陣,g1,g2,g3均為正常數(shù)。為外部擾動估計矩陣,表示向量張成的對角矩陣,是對d1,d2,d3的估計值;
若設(shè)計和的自適應(yīng)律滿足:
其中ε1,ε2,m,n,γ1,γ2,γ3均為正常數(shù),分別表示滾轉(zhuǎn)角速度誤差、俯仰角速度誤差和偏航角速度誤差,則姿態(tài)誤差四元數(shù)和角速度誤差漸近穩(wěn)定。
驗證步驟具體是,采用基于Lyapunov的分析方法可以證明當(dāng)時間趨于無窮時,和ev分別漸近收斂到[0 0 0]T。