柴嘉薪,王 聰,施國(guó)興,王 輝
(北京航天自動(dòng)控制研究所,北京,100854)
在運(yùn)載火箭執(zhí)行任務(wù)過(guò)程中,慣性器件作為控制系統(tǒng)的關(guān)鍵單機(jī),可能因自身質(zhì)量問(wèn)題或受到外部環(huán)境的影響發(fā)生故障,直接影響任務(wù)成敗。
提升慣組的可靠性,可從提升單一元部件的可靠性入手,但受生產(chǎn)工藝水平等因素限制,此項(xiàng)研究投入大、周期長(zhǎng)、見(jiàn)效慢,提升空間有限[1],且無(wú)法完全避免未知的外部環(huán)境因素對(duì)慣性器件性能造成的不良影響。傳統(tǒng)的運(yùn)載火箭設(shè)計(jì)中通常采用慣組冗余技術(shù),以大幅度提高運(yùn)載火箭控制系統(tǒng)的可靠性[2]。
對(duì)于重大工程任務(wù),須遵循“一度故障正常飛行,二度故障保證安全”的設(shè)計(jì)原則,中國(guó)工程實(shí)踐中現(xiàn)有可滿足此項(xiàng)原則的方案僅有雙七表?yè)闲詰T組冗余方案[3]。隨著光學(xué)陀螺技術(shù)的發(fā)展,成熟度逐步提高,光學(xué)陀螺因其測(cè)量精度高、靈敏度高、動(dòng)態(tài)范圍大、體積小、重量輕、易于集成的優(yōu)點(diǎn),已成為航天應(yīng)用的主流,目前航天主流運(yùn)載火箭除中國(guó)的載人運(yùn)載火箭外,均無(wú)機(jī)械陀螺構(gòu)成的慣組,已全部更新?lián)Q代為光學(xué)陀螺慣組[4-5],因此不再考慮沿用此成熟方案。
十表光學(xué)慣組可組成成本較低、節(jié)約空間、質(zhì)量較輕且冗余度更高的慣組冗余配置?,F(xiàn)已有成熟的單十表光學(xué)慣組冗余方案[6]經(jīng)過(guò)了多型運(yùn)載火箭飛行試驗(yàn)驗(yàn)證,該方案可實(shí)現(xiàn)慣組在一度故障下正常工作,但會(huì)損失一定的導(dǎo)航精度,對(duì)于重型運(yùn)載火箭的長(zhǎng)時(shí)間飛行任務(wù),無(wú)法保證一度故障下的高精度入軌。部分文獻(xiàn)中研究了斜置表與正交軸相對(duì)角度的最優(yōu)化配置設(shè)計(jì),以提高單十表慣組配置下斜表導(dǎo)航的精度[7-8]。但受方案冗余度和斜置表性能的限制,單十表光學(xué)慣組冗余方案仍有較大局限性。
針對(duì)以上問(wèn)題,本文研究了一種基于雙十表光學(xué)捷聯(lián)慣組的主從冗余方法及信息重構(gòu)技術(shù)。依據(jù)在主/從單套慣組內(nèi)分別對(duì)五表陀螺儀和五表加速度計(jì)進(jìn)行故障診斷和進(jìn)一步通過(guò)主從慣組之間一致性判別策略實(shí)現(xiàn)故障定位的結(jié)果,完成異常表的故障計(jì)數(shù)和故障表信息的隔離。通過(guò)冗余信息重構(gòu)技術(shù),重型運(yùn)載火箭在長(zhǎng)時(shí)間飛行過(guò)程中,慣組的同類表發(fā)生三度以下故障后,控制系統(tǒng)仍能通過(guò)系統(tǒng)級(jí)冗余得到可靠有效的慣性導(dǎo)航信息。
雙十表光學(xué)捷聯(lián)慣組配置了雙五表陀螺儀和雙五表加速度計(jì),其組成如圖1 所示。在箭體坐標(biāo)系OX1Y1Z1下,定義Ax、Ay、Az為3個(gè)正交安裝的石英加速度計(jì);AS和AT為兩個(gè)斜置安裝的石英加速度計(jì);GX、GY、GZ為3個(gè)正交安裝的光學(xué)陀螺儀,GS和GT為兩個(gè)斜置安裝的光學(xué)陀螺儀,圖中箭頭指向即為對(duì)應(yīng)表的敏感軸向。兩套慣組中斜表軸向與正交軸之間的定向夾角一致,其中OS軸與三正交軸之間的夾角分別為αS、βS、γS,OT軸與三正交軸之間的夾角分別為αT、βT、γT。
圖1 雙十表光學(xué)捷聯(lián)慣組的組成及安裝定向示意Fig.1 Composition and installation orientation of double tensensor optical strapdown inertial unit
在不考慮兩套慣組之間的安裝誤差及擺放不水平度差異的理想情況下,主慣組與從慣組對(duì)應(yīng)軸向陀螺儀/加速度計(jì)的測(cè)量信息是一致的。除了主從慣組之間的測(cè)量信息比對(duì)外,每套慣組同類表中兩個(gè)斜置表與3 個(gè)正交表輸出的測(cè)量增量值另有投影關(guān)系如下(此處以陀螺儀輸出為例,加速度計(jì)輸出同理):
式 中δθαi為 各 軸 陀 螺 儀 輸 出 的 角 速 度,α=x、y、z、s、t;Esi_x,Esi_y,Esi_z分別為GS表在各正交軸的投影系數(shù);Eti_x、Eti_y、Eti_z分別為GT表在各正交軸的投影系數(shù);i= 1、2,分別為主慣組和從慣組。
雙十表光學(xué)慣組冗余信息處理采用主從冗余方案。首先,在主慣組和從慣組內(nèi)分別進(jìn)行單套慣組內(nèi)五表陀螺儀/五表加速度計(jì)的故障診斷與故障隔離。對(duì)于單套慣組內(nèi)出現(xiàn)同類表二度及二度以上一致性故障的情況,再進(jìn)行主從慣組之間的一致性故障檢測(cè),以判定故障慣組各正交表是否損壞。利用隔離故障表后的冗余測(cè)量信息進(jìn)行冗余信息重構(gòu)與導(dǎo)航計(jì)算,分別通過(guò)主、從慣組五表陀螺儀/五表加速度計(jì)測(cè)量信息得到角增量信息和視速度增量信息。最終根據(jù)故障定位結(jié)果和冗余信息使用真值表,從主慣組或從慣組中擇優(yōu)選取角增量/視速度增量信息用于導(dǎo)航控制。
慣性器件的故障模式多且機(jī)理復(fù)雜,但從慣組冗余信息使用的角度來(lái)說(shuō),主要關(guān)注慣性器件的輸出特性是否正常,雙十表慣組冗余信息處理流程見(jiàn)圖2。
圖2 雙十表光學(xué)慣組冗余信息處理流程Fig.2 Redundant information processing flow of double tensensor optical strapdown inertial unit
單套慣組內(nèi)故障診斷措施主要包括陀螺儀/加速度計(jì)脈沖常零值輸出判別、角速度/視加速度極大值判別、視加速度/角速度/角度一致性判別。其中,常零值輸出判別與極大值判別的方法簡(jiǎn)單、明確且成熟,可直接沿用現(xiàn)有方法[9],但需根據(jù)任務(wù)特點(diǎn)和慣性器件技術(shù)指標(biāo)謹(jǐn)慎設(shè)計(jì)故障判別和故障隔離的門限。
針對(duì)陀螺儀測(cè)量信息和加速度計(jì)測(cè)量信息的非零值故障和非極大值故障,采用正交軸與斜表在正交軸投影共判、短周期與長(zhǎng)周期相結(jié)合的方式,將一度故障直接定位到故障單表或診斷出二度故障,并根據(jù)故障判別和定位結(jié)果分情況進(jìn)行后續(xù)的故障隔離和冗余信息使用。
以陀螺儀輸出信息為例,介紹一種雙五表配置下的一致性故障判別方法。
針對(duì)單套慣組內(nèi)同類表的一致性故障判別,以陀螺儀輸出的角速度信息為例,建立如下一致性判別方程:
式中δθαi0為陀螺儀輸出的角速度;εGiα為故障門限值;α=x、y、z、s、t;i= 1、2,分別為主慣組和從慣組;Δt為滾動(dòng)累加周期,為降低高頻干擾信號(hào)對(duì)故障判別的影響,對(duì)角速度一致性判別進(jìn)行滾動(dòng)累加,本文取滾動(dòng)累加周期Δt為控制周期的25倍。
在單套慣組內(nèi),利用建立的一致性判別方程,對(duì)五表陀螺儀分別進(jìn)行故障表一致性檢測(cè),并執(zhí)行故障表隔離操作。分別判斷主慣組和從慣組內(nèi)是否存在五表陀螺儀出現(xiàn)一度故障或二度故障的情況。主、從慣表故障定位見(jiàn)表1、表2。
表1 主慣組五表陀螺儀輸出信息一致性故障定位Tab.1 Consistency fault location of gyroscope output information from main strapdown inertial unit
表2 從慣組五表陀螺儀輸出信息一致故障定位Tab.2 Consistency fault location of gyroscope output information from slave strapdown inertial unit
設(shè)單套慣組的一致性故障標(biāo)志字N主/N從,用以表示單套慣組內(nèi)部一致性故障判別的結(jié)果:
結(jié)合表1 和表2 單套慣組內(nèi)五表陀螺儀輸出信息一致性故障的定位結(jié)果,可得到主、從慣組的故障標(biāo)志字見(jiàn)表3。
表3 單套慣組內(nèi)一致性故障判別結(jié)果Tab.3 Consistent fault discrimination results in a single strapdown inertial unit
單套慣組內(nèi)五表加速度計(jì)一致性故障判別方法與五表陀螺儀類似,不同之處在于,對(duì)五表陀螺儀的一致性故障判別采用了短周期累加和長(zhǎng)周期累加相結(jié)合的方式。但對(duì)于五表加速度計(jì)一致性故障的判別,考慮到在飛行狀態(tài)下火箭軸向過(guò)載與橫法向過(guò)載差異很大,軸向過(guò)載在橫法向投影,不同慣組之間的安裝誤差會(huì)導(dǎo)致橫法向表輸出有差異,長(zhǎng)時(shí)間累計(jì)下來(lái)對(duì)單表輸出影響很大。若隨之放大故障判別門限將使判別失去意義,因此五表加速度計(jì)采用短周期累加和中等周期累加相結(jié)合的一致性故障判別方式。
針對(duì)第2 節(jié)中通過(guò)式(3)~(12)判斷出主慣組和從慣組中有一個(gè)慣組出現(xiàn)同類表二度故障,且另一套慣組無(wú)故障的情況,需要進(jìn)一步通過(guò)主從慣組之間一致性故障判別確定發(fā)生一度以上的慣組的正交表是否出現(xiàn)故障。
針對(duì)主從慣組之間的一致性故障判別,以陀螺儀輸出的角速度信息為例,建立如下一致性判別方程:
式中δθαi0為陀螺儀輸出的角速度;εGiα為故障門限值;α=x,y,z,s,t;i= 1、2,分別為主慣組和從慣組;Δt為滾動(dòng)累加周期。
針對(duì)第2.2 節(jié)中判斷出主慣組和從慣組中有一個(gè)慣組出現(xiàn)同類表二度故障,且另一套慣組無(wú)故障的情況,利用式(14)~(16)進(jìn)一步判別出現(xiàn)二度及二度以上故障的慣組中3個(gè)正交表是否為故障表,診斷結(jié)果用于后續(xù)的冗余信息重構(gòu)。
針對(duì)第2.2 節(jié)中判斷出一套慣組為二度及二度以上故障,另一套慣組為一度故障或二度及二度以上故障的情況,屬于系統(tǒng)三度或三度以上故障,考慮控制系統(tǒng)的工程需求與復(fù)雜性,不再針對(duì)此類情況作進(jìn)一步故障表判別處理。
同類表二度及二度以上慣組正交表故障判別真值見(jiàn)表4。
表4 故障慣組中正交表故障判別Tab.4 Fault discrimination of orthogonal axis sensors in fault strapdown inertial unit
根據(jù)表4單套慣組之間正交表故障判別結(jié)果,得到主、從慣組故障標(biāo)志字見(jiàn)表5。
表5 兩套捷聯(lián)慣組之間一致性故障判別結(jié)果Tab.5 Consistent fault discrimination results between two strapdown inertial units
與雙五表陀螺儀類似,雙五表加速度計(jì)經(jīng)過(guò)故障判別與故障隔離,在主慣組和從慣組內(nèi)分別利用各自的無(wú)故障三正交加速度計(jì)測(cè)量信息或冗余信息重構(gòu)結(jié)果進(jìn)行導(dǎo)航解算得到視速度增量,待系統(tǒng)進(jìn)行冗余信息的選用。
通過(guò)冗余信息故障診斷,得到當(dāng)前五表陀螺儀/加速度計(jì)輸出的測(cè)量信息的健康狀態(tài),基于此,開(kāi)展單慣組內(nèi)同類表冗余信息的重構(gòu)與導(dǎo)航解算。具體流程如圖3所示。
圖3 冗余信息重構(gòu)流程Fig.3 Redundant information reconstruction process
如圖3所示,對(duì)于在當(dāng)前控制周期內(nèi)未達(dá)到三度故障切除的五表陀螺儀/加速度計(jì),依次進(jìn)行是否有故障計(jì)數(shù)、是否有正交表異常計(jì)數(shù)或正交表故障、是否有兩個(gè)正交表異常計(jì)數(shù)或故障、是否有斜置表異常計(jì)數(shù)或故障等判斷。針對(duì)各步驟判斷結(jié)果,依據(jù)冗余信息重構(gòu)流程,利用相應(yīng)的測(cè)量信息進(jìn)行導(dǎo)航解算或信息重構(gòu)。
慣組選型時(shí),受成本控制、安裝空間和總質(zhì)量要求等方面的約束,通常配置在正交軸方向上的陀螺儀和加速度計(jì)比斜置軸方向上的精度高,為了實(shí)現(xiàn)更高的導(dǎo)航制導(dǎo)精度,在故障狀態(tài)下系統(tǒng)重構(gòu)時(shí),優(yōu)先選用無(wú)故障的正交表提供的測(cè)量信息進(jìn)行導(dǎo)航計(jì)算。
隔離故障慣性儀表并完成冗余信息重構(gòu)后,系統(tǒng)的故障適應(yīng)能力明顯提高,可以應(yīng)對(duì)二度及二度以上的同類單表故障。
當(dāng)五表陀螺儀/加速度計(jì)的正交軸出現(xiàn)故障時(shí),通過(guò)信息重構(gòu)可由其他斜置軸的冗余信息準(zhǔn)確估計(jì)出故障正交軸向的導(dǎo)航信息,具體方法如下:
五表配置下的陀螺儀/加速度計(jì)的輸出方程為[10]
式中Z∈R5為五軸向慣性儀表的測(cè)量信息;ω∈R3為系統(tǒng)三正交軸的角增量/視速度增量信息;H為慣組的測(cè)量矩陣;ε為均值為零的高斯白噪聲。
將兩斜置軸OS和OT分別投影在XOZ平面上,得到斜置軸投影線OS'和OT',如圖4 所示,αS、αT分別為從OX1軸順時(shí)針旋轉(zhuǎn)至OS和OT軸的角度,βS'、βT'分別為從OY1軸逆時(shí)針轉(zhuǎn)至OS'和OT'軸的角度。
圖4 十表光學(xué)慣組斜置軸投影示意Fig.4 Schematic of oblique axis projection in ten-sensor optical strap down inertial unit
假設(shè)陀螺儀在五表配置下,經(jīng)故障檢測(cè)到安裝Y軸的陀螺儀GY發(fā)生故障,其輸出信息被系統(tǒng)進(jìn)行故障隔離,需要使用其他無(wú)故障陀螺儀進(jìn)行信息重構(gòu)。
a)使用無(wú)故障表GX、GZ、GS進(jìn)行信息重構(gòu)。
由此,得到:
在當(dāng)前假設(shè)條件下,wx=ZGX,wz=ZGZ,可以進(jìn)一步得到GY的正確輸出為
b)使用無(wú)故障表GX、GZ、GT進(jìn)行信息重構(gòu)??傻玫剑?/p>
在當(dāng)前假設(shè)條件下,wx=ZGX,wz=ZGZ,可以進(jìn)一步得到GY的正確輸出為:
c)使用無(wú)故障表GX、GS、GT進(jìn)行信息重構(gòu)。可得到:
求解可得到:
d)使用無(wú)故障表GZ、GS、GT進(jìn)行信息重構(gòu)??傻玫剑?/p>
求解可得到:
基于第2.2節(jié)和第2.3節(jié)中以陀螺儀為例的各單套慣組內(nèi)部同類表一致性故障判別和兩慣組之間一致性故障判別和定位結(jié)果,隔離慣組中故障陀螺儀/加速度計(jì)的測(cè)量輸出信息。再利用冗余信息重構(gòu)技術(shù)分別得到由兩慣組測(cè)量輸出導(dǎo)航解算得到的角增量/視速度增量。最后,遵從優(yōu)先使用主慣組信息以及盡量選用無(wú)故障的正交表慣組測(cè)量信息進(jìn)行導(dǎo)航計(jì)算這兩項(xiàng)原則,完成冗余信息的選用,依據(jù)冗余信息使用真值表確定用于導(dǎo)航控制的角增量/視速度增量信息,使得在慣組發(fā)生故障的情況下,火箭經(jīng)過(guò)長(zhǎng)時(shí)間飛行后仍能盡可能減少精度損失,確保在本文提出的慣組冗余方案冗余度可覆蓋的故障模式下實(shí)現(xiàn)精確入軌。冗余信息使用真值見(jiàn)表6。
表6 冗余信息使用真值Tab.6 Truth for redundant information
本文針對(duì)重型運(yùn)載火箭的任務(wù)特點(diǎn)及可靠性需求,提出了一種基于雙十表光學(xué)捷聯(lián)慣組的主從冗余方法與信息重構(gòu)技術(shù)。首先,給出了冗余信息處理方案,并設(shè)計(jì)了雙五表陀螺儀/雙五表加速度計(jì)一致性故障診斷算法。再基于故障診斷結(jié)果,設(shè)計(jì)了冗余信息重構(gòu)流程以及十表慣組正交軸慣性儀表故障下的信息重構(gòu)算法。最后,遵從優(yōu)先使用主慣組信息以及盡量選用無(wú)故障的正交表慣組測(cè)量信息進(jìn)行導(dǎo)航計(jì)算這兩項(xiàng)原則,完成冗余信息的選用。
本文提出的技術(shù)方法實(shí)現(xiàn)了在雙十表光學(xué)慣組配置下最大限度地對(duì)故障進(jìn)行檢出與隔離,滿足長(zhǎng)時(shí)間飛行任務(wù)中,在慣組三度以下同類表故障的情況下,控制系統(tǒng)仍可靠工作的需求,并確保運(yùn)載火箭完成發(fā)射任務(wù)。