(中國空間技術研究院,北京 100094)
衛(wèi)星發(fā)動機是衛(wèi)星變軌的主要設備,發(fā)動機安裝在與衛(wèi)星中心承力筒相連的發(fā)動機支架上,其理論推力方向與衛(wèi)星機械坐標系的Z軸相同。由于發(fā)動機推力矢量偏差(即發(fā)動機實際推力方向和作用點與理論值存在偏差)及衛(wèi)星質心與機械坐標系原點不重合等因素的影響,導致發(fā)動機推力偏斜,在發(fā)動機工作時會存在干擾力矩。為減小或避免干擾力矩,發(fā)動機在安裝時,會有調整環(huán)節(jié),即通過確定安裝參數來調整發(fā)動機的實際推力方向和作用點。
目前,衛(wèi)星僅根據發(fā)動機推力矢量實際參數來確定發(fā)動機安裝參數,未考慮衛(wèi)星質心對干擾力矩的影響,也未考慮調整過程中的幾何限制問題[1]。這樣就存在實際干擾力矩不是最小的情況,甚至可能出現(xiàn)發(fā)動機調整后比調整前的干擾力矩更大的極端情況;同時,在某些情況下,設定的發(fā)動機安裝參數可能會因為發(fā)動機安裝時的物理限制而無法實現(xiàn)。
本文首先介紹了發(fā)動機安裝調整模型,其次分析了發(fā)動機安裝優(yōu)化命題的數學表述,最后給出了一個遺傳算法應用的優(yōu)化示例。在確定發(fā)動機安裝參數時,綜合考慮了發(fā)動機推力矢量參數、衛(wèi)星質心和發(fā)動機安裝相關的空間幾何約束等3方面因素。這些約束多為非線性約束,一般的優(yōu)化算法較難解決此類問題,所以本文選擇遺傳算法來解決這一優(yōu)化問題。
發(fā)動機與衛(wèi)星之間通過發(fā)動機支架連接,如圖1所示。發(fā)動機通過發(fā)動機安裝法蘭上的3個安裝孔與發(fā)動機支架相連。發(fā)動機安裝法蘭上的3個安裝孔的孔徑(即圖1中的D)大于安裝螺釘的外徑(即圖1中的d),這就使得發(fā)動機的安裝有了調整能力,可以通過圖1所示的9來調節(jié)發(fā)動機的安裝參數。
圖1 發(fā)動機安裝示意圖Fig.1 Sketch map of orbit-control engine installation
2.2.1 發(fā)動機調整相關坐標系
發(fā)動機的調整,牽涉到的坐標系有衛(wèi)星機械坐標系OscXscYscZsc和發(fā)動機安裝坐標系OnXnYnZn以及發(fā)動機坐標系OXYZ,3個坐標系之間的關系如圖2所示。其中,發(fā)動機安裝坐標系相對衛(wèi)星機械坐標系固定,兩者之間僅在縱軸上有h的相對位移;發(fā)動機坐標系相對發(fā)動機安裝坐標系在空間上有旋轉和平移(即圖2中的M1、M2和M3)。為簡化計算過程,本文后面的平移、旋轉等操作均相對衛(wèi)星機械坐標系;在描述發(fā)動機推力矢量參數時,使用發(fā)動機坐標系。
圖2 發(fā)動機安裝過程中各坐標系關系示意圖Fig.2 Relationship between several coordinates
2.2.2 發(fā)動機調整模型
發(fā)動機的理論推力方向與發(fā)動機的幾何軸線是重合的,即發(fā)動機未調整時,其理論推力主軸與衛(wèi)星機械坐標系的OscZsc軸重合。實際上,發(fā)動機的推力方向與發(fā)動機的幾何軸線成一定角度,作用點也有偏移。這可以通過推力矢量參數來表示。
發(fā)動機推力矢量參數的示意如圖3所示。
圖3中,OXYZ為發(fā)動機坐標系,其它參數定義如下:
F—發(fā)動機推力,單位為N;
FZ—發(fā)動機推力F在OZ軸上的投影;
FXOY—發(fā)動機推力F在XOY面內的投影;
P—發(fā)動機推力作用點;
α—推力矢量偏斜角(簡稱偏斜角),發(fā)動機推力F與推力理論軸線OZ軸夾角,單位為(°);
β—推力矢量橫移位置角(簡稱橫移位置角),OP與OX軸之間夾角,由發(fā)動機頂視方向逆時針為正,單位為(°);
γ—推力矢量偏斜位置角(簡稱偏斜位置角),側向推力FXY與OX軸之間夾角,由發(fā)動機頂視方向逆時針為正,單位為(°);
δ—推力矢量橫移量(簡稱橫移量),發(fā)動機推力作用點P到發(fā)動機坐標系原點O的距離,單位為mm。
圖3 發(fā)動機推力矢量參數在發(fā)動機坐標系中的示意Fig.3 Orbit-control engine’s thrust vector in engine coordinate system
為盡量減小干擾力矩,通??赏ㄟ^在發(fā)動機和發(fā)動機支架之間安裝調整墊片的方式(即圖1中的9),來調整發(fā)動機推力作用點和方向,使推力方向逼近衛(wèi)星機械坐標系的OscZsc。具體調整的原則如下:
不考慮橫移位置角和橫移量,僅考慮偏斜角和偏斜位置角,為保證調整后,發(fā)動機的推力方向與衛(wèi)星機械坐標系的OscZsc軸的夾角小于某值,按如下方法設置發(fā)動機向衛(wèi)星機械坐標系OscZsc軸調整的角度:當V1<α≤V2時,調整角度為0.5α;當α>V2時,調整角度為α-V2/2(α為發(fā)動機推力矢量偏斜角,V1和V2分別為調整角度參考值的下限和上限)。
上述調整方法,存在以下問題:
(1)未考慮橫移量等其它推力矢量參數,調整考慮的綜合因素不全;
(2)調整原則是將發(fā)動機推力方向盡量逼近衛(wèi)星機械坐標系的OscZsc軸,不是逼近衛(wèi)星質心所在的縱軸,這樣調整后,推力可能更偏離質心,干擾力矩反而加大;
(3)調整未考慮發(fā)動機和發(fā)動機支架之間機械接口的約束,可能設定的調整參數因為機械干涉而無法實現(xiàn)。
為解決上述問題,本文提出一種新的發(fā)動機安裝調整方法,計算過程如下:調整后的發(fā)動機推力作用點位置矢量為
式中:h為發(fā)動機安裝坐標系相對衛(wèi)星機械坐標系在縱軸上的相對位移;M1、M2和M3為發(fā)動機坐標系相對發(fā)動機安裝坐標系在空間上的平移量;[p1p2p3]為發(fā)動機推力作用點位置矢量P在衛(wèi)星機械坐標系下的空間位置;R為發(fā)動機調整后,發(fā)動機坐標系相對衛(wèi)星機械坐標系的旋轉矩陣,如式(2)所示。
式中:α1,α2,α3分別為發(fā)動機坐標系繞衛(wèi)星機械坐標系3個坐標軸的旋轉角度。
發(fā)動機推力為
式中:f1,f2,f3分別為發(fā)動機推力矢量在衛(wèi)星機械坐標系3個坐標軸上的分量值;F0為發(fā)動機理論推力。
發(fā)動機干擾力矩為
式中:t1、t2、t3分別為繞衛(wèi)星機械坐標系的OscXsc軸、OscYsc軸、OscZsc軸的干擾力矩分量,[xCyCzC]為衛(wèi)星質心在衛(wèi)星機械坐標系下的值。推力作用點和質心的單位為mm,推力的單位為N,干擾力矩單位為N·m。
式(1)~(4)給出了一種新的發(fā)動機安裝調整參數與干擾力矩的關系。發(fā)動機安裝調整的目的是如何確定α1、α2、α3以及M1、M2、M3,使干擾力矩T最小。從式(1)~(4)看,該問題屬于非線性優(yōu)化問題,可以利用遺傳算法求解。
遺傳算法(Genetic Algorithm,GA)是以自然選擇和遺傳機制為理論基礎,將生物進化過程中適者生存規(guī)則與種群內部染色體的隨機信息交換機制相結合的一種高效的全局尋優(yōu)搜索算法。遺傳算法主要應用在函數優(yōu)化[2]、組合優(yōu)化[3-4]、生產調度[5]等領域。遺傳算法除能處理線性問題的優(yōu)化外,還能處理非線性問題的優(yōu)化[6-7]。
按照遺傳算法解決問題的特點,確定優(yōu)化問題的數學表達為
式中:minT(xi)即求取干擾力矩最小值,為目標函數;xi為變量值,對應α1,α2,α3和M1,M2,M36個變量;X對應6個變量的取值空間,即約束條件。
發(fā)動機安裝優(yōu)化的目的是使得變軌期間的干擾力矩最小。一般情況下,衛(wèi)星最后一次變軌時的質心情況最為惡劣,干擾力矩也最大。所以,本文以最后一次變軌時的衛(wèi)星質心作為計算的輸入,以最后一次變軌時的綜合干擾力矩,即繞衛(wèi)星機械坐標系3個坐標軸的干擾力矩(式(6)中的t1,t2,t3)的平方根值最小作為優(yōu)化的目標函數:
發(fā)動機安裝的約束主要為幾何約束,下面結合發(fā)動機安裝情況進行分析。
結合圖1的發(fā)動機安裝方式和圖4所示的發(fā)動機安裝法蘭,可看出發(fā)動機安裝時,受到的幾何約束有:
(1)由于發(fā)動機支架相對衛(wèi)星是固定的,而發(fā)動機與支架之間的安裝螺釘(即圖1中的7)也是固定的,所以發(fā)動機安裝法蘭上的安裝孔在旋轉、平移時不得超過安裝螺釘的約束;
(2)發(fā)動機安裝法蘭的上表面邊緣(即圖1中的4)不得高于發(fā)動機安裝支架下表面,發(fā)動機安裝法蘭上3個安裝孔的下表面圓周上任意一點在衛(wèi)星機械坐標系下的Z坐標值不得低于安裝用墊圈(即圖1中的8)的上表面;
(3)發(fā)動機繞衛(wèi)星機械坐標系的OscZsc軸旋轉角度不得過大(考慮發(fā)動機燃口與管路的空間位置關系)。
圖4為發(fā)動機安裝法蘭示意圖(圖4中,H1、H2、H3處指示的即為發(fā)動機安裝法蘭上的3個安裝孔)。
將發(fā)動機安裝受到的上述3個約束轉化為表達式見下文。
圖4 發(fā)動機安裝法蘭俯視圖Fig.4 Top view of orbit-control engine’s fixing flange
3.3.1 第1個約束條件的數學表述
第1個約束條件可表述為
式中:C′ix_up和C′iy_up,C′ix_down和C′iy_down分別為發(fā)動機調整后,發(fā)動機安裝法蘭上3個安裝孔的上表面、下表面的圓周上任意一點在衛(wèi)星機械坐標系下的X、Y坐標值;Hix和Hiy分別為發(fā)動機調整前,發(fā)動機安裝法蘭上3個連接孔的中心在衛(wèi)星機械坐標系下的X、Y坐標值;i=1,2,3。
C′ix_up、C′iy_up、C′iz_up和C′ix_down、C′iy_down、C′iz_down的計算見式(8)。
式中:l為發(fā)動機安裝法蘭的厚度;C′iz_up、C′iz_down為發(fā)動機調整后,發(fā)動機安裝法蘭上3個安裝孔的上表面、下表面的圓周上任意一點在衛(wèi)星機械坐標系下的Z坐標值;i=1,2,3;κ∈[0°,360°)。
Hix和Hiy的計算見式(9)。
3.3.2 第2個約束條件的數學表述
第2個約束條件可表述為
式中:Zmax為發(fā)動機支架法蘭盤(即圖1中的6)下表面在衛(wèi)星機械坐標系下的Z向坐標值;Z′min為安裝用墊圈(即圖1中的8)上表面在衛(wèi)星機械坐標系下的Z向坐標值;Piz為發(fā)動機調整后,發(fā)動機安裝法蘭上6個凸耳的外緣在衛(wèi)星機械坐標系下的Z坐標值,i=1,2,…,6,其計算見式(11)。
式中:Pix和Piy分別為發(fā)動機調整后,發(fā)動機安裝法蘭上6個凸耳的外緣在衛(wèi)星機械坐標系下的X、Y坐標值;ψi為發(fā)動機安裝法蘭上6個凸耳(見圖4)上任意一點相對坐標軸OX的張角;i=1,2,…,6;j=4,5,…,9。
3.3.3 第3個約束條件的數學表述
第3個約束條件可表述為
式中:φ為限制發(fā)動機不能繞衛(wèi)星機械坐標系OscZsc軸旋轉過大的參數。
第3節(jié)給出了約束優(yōu)化問題的目標函數和約束條件,分別將目標函數(式(1)~(6))轉換為代碼,保存為Myfit函數;將約束條件(式(7)~(12))轉換為代碼,保存為Myconst函數(內部包含等式約束和不等式約束)。然后,利用Matlab軟件中的Gaoptimset函數設置種群規(guī)模和遺傳代數,終止條件可取默認值,使用GA 函數[8]運行程序,輸出干擾力矩最小時的發(fā)動機旋轉角α1,α2,α3和平移量M1,M2,M3。
在本問題中使用GA 函數進行尋優(yōu)時,群規(guī)模和遺傳代數及終止條件對結果并不敏感。
設定某臺發(fā)動機的偏斜角為0.25°,橫移位置角為282°,偏斜位置角為96°,橫移量為0.70 mm,衛(wèi)星質心為[5mm,-5mm,1550mm]。若發(fā)動機不調整,干擾力矩平方和的均方根值為2.73N·m。若按照目前的調整方法(即2.2.2節(jié)所述),經過調整后得到的干擾力矩平方和的均方根值為3.74N·m,該值超過了調整前的結果,其原因在于未考慮衛(wèi)星質心,所以會越調越偏。
其它相關參數不變,采用本文提出的優(yōu)化方法,在Gaoptimset函數中,設置種群規(guī)模為20,遺傳代數為30,程序代碼為
options=gaoptimset('populationsize',20,'generation',30,'plotfcns',{@gaplotbestf,@gaplotbestindiv});
[x,Tmin,flag,output,population,scores]=ga(@(x)myfit(x),6,[],[],[],[],[],[],@(x)myconstr(x),options)
x為1×6的數組,即α1,α2,α3以及M1,M2,M3,Tmin即為發(fā)動機安裝優(yōu)化的目標函數。
程序運行后,得到圖5所示結果。其中,圖5(a)為目標函數的優(yōu)化過程,橫坐標為優(yōu)化代數,縱坐標為目標函數值的優(yōu)化結果,得到的最小干擾力矩值為5.87×10-4N·m,此時的發(fā)動機旋轉角和平移量見圖5(b)(其中變量1~6即為目標函數最小時的α1,α2,α3以及M1,M2,M3的值)。
圖5 發(fā)動機安裝參數優(yōu)化Fig.5 Result of orbit-control engine installation parameters after being optimized
本文針對發(fā)動機安裝調整要求中考慮因素不全面的問題,提出一種綜合考慮全部因素的發(fā)動機調整方法,并將該方法利用遺傳算法進行了求解。相對目前的發(fā)動機調整方法來說,本文所提的方法有以下3點不同:
(1)考慮橫移位置角和橫移量后,增加平移調整環(huán)節(jié);
(2)確定發(fā)動機調整角度時,不僅考慮了發(fā)動機推力矢量參數,而且綜合了衛(wèi)星質心對干擾力矩的影響;
(3)發(fā)動機安裝參數的確定,考慮了發(fā)動機調整過程中的幾何約束問題。
通過對設定發(fā)動機參數的仿真,證明該方法比目前的方法更為科學,可更進一步減小衛(wèi)星變軌期間的干擾力矩,從而節(jié)省衛(wèi)星燃料,提高衛(wèi)星在軌服務壽命。
(References)
[1]胡偉,金煌煌,王慶華.衛(wèi)星490N 發(fā)動機安裝位置優(yōu)化研究[J].航天器工程,2013,22(6):44-47 Hu Wei,Jin Huanghuang,Wang Qinghua.Research on optimizing mounting position of satellite 490N thurster[J].Spacecraft Engineering,2013,22(6):44-47(in Chinese)
[2]Deng J J.A hybrid genetic algorithm for function optimization[J].Journal of Software.1999,10(8):819-823
[3]Zhang L,Wang L.Genetic ordinal optimization for stochastic traveling salesman problem[C]//Proceedings of the World Congress on Intelligent Control and Automa-tion(WCICA).Beijing:Dept.of Autom.,Tsinghua University,2004(3):2086-2090
[4]趙光,范陸海,馮碩.星載TWTA 備份環(huán)開關切換的遺傳算法研究[J].航天器工程,2011,20(5):93-98 Zhao Guang,F(xiàn)an Luhai,F(xiàn)eng Shuo.Redundancy switching solution of spaceborne TWTA using Genetic Algorithm [J].Spacecraft Engineering,2011,20 (5):93-98(in Chinese)
[5]Abe M,Matsumoto H,Kuroda C.An artificial neural network optimized by agenetic algorithm for real-time flow-shop rescheduling[J].International Journal of Knowledge-Based Intelligent Engineering Systems,2002,6(2):96-103
[6]蘇三買,廉小純.遺傳算法在航空發(fā)動機非線性數學模型中的應用[J].推進技術,2004,25(3):237-240 Su Sanmai,Lian Xiaochun.Application of genetic algorithm in aeroengine nonlinear mathematical models[J].Journal of Propulsioon Technology,2004,25(3):237-240(in Chinese)
[7]王玉瑩,李運澤,劉東曉.納衛(wèi)星隔熱層厚度與散熱面面積優(yōu)化設計[J].航天器工程,2010,19(2):46-51 Wang Yuying,Li Yunze,Liu Dongxiao.Optimized design of heat insulation layers'thickness and radiator surfaces for nano-satellite[J].Spacecraft Engineering,2010,19(2):46-51(in Chinese)
[8]杜東,馬震,孫曉明.MATLAB 遺傳算法工具箱(GAOT)在水資源優(yōu)化計算中的應用[J].水利科技與經濟,2007,13(2):73-75,78 Du Dong,Ma Zhen,Sun Xiaoming.Application of MATLAB genetic algorithms toolbox to optimization calculation of water resource[J].Water Conservancy Science and Technology and Economy,2007,13(2):73-75,78(in Chinese)