撓性衛(wèi)星快速機動光滑控制

饒衛(wèi)東

1. 北京控制工程研究所,北京 100190 2. 空間智能控制技術重點實驗室,北京 100190

?

撓性衛(wèi)星快速機動光滑控制

饒衛(wèi)東1,2

1. 北京控制工程研究所,北京 100190 2. 空間智能控制技術重點實驗室,北京 100190

著力求出一種光滑的正弦序列控制，該控制律執(zhí)行撓性衛(wèi)星快速機動任務的同時抑制振動。推導得到撓性振動的解析解，當正弦序列加權系數滿足一定條件時，中心剛體自動滿足靜止-靜止邊界條件，然后通過非線性規(guī)劃優(yōu)化加權系數，得到一組殘余振動指標次優(yōu)解。相關公式表明該優(yōu)化解能較大程度抑制撓性振動，且相應殘余振動對模態(tài)參數魯棒性好。數學仿真在5階模態(tài)撓性對象上開展，結果顯示該方法有較強的抑振能力，對頻率誤差不敏感，且控制方法是開環(huán)的，適合作為衛(wèi)星快速機動前饋輸入。 關鍵詞 撓性衛(wèi)星； 快速機動； 前饋輸入； 非線性規(guī)劃

航天任務中質量更輕意味著成本更低廉，能量消耗更少。航天器的輕型化和大型化使得大撓性問題日益突出，這使得航天器在姿態(tài)機動時，更容易產生撓性振動。而太空環(huán)境的微重力和微阻尼，又使得振動很難衰減。很多航天任務中姿態(tài)控制都遇到了困難，比如機械臂抓取、移動重物，敏捷衛(wèi)星姿態(tài)快速機動等。如何在這些任務中既快速又穩(wěn)定地控制撓性對象是一個重大挑戰(zhàn)。

研究人員經過幾十年積累，充分解答了撓性衛(wèi)星的時間最優(yōu)控制。Singh等[1]給出了純粹時間最優(yōu)控制的數值解法，得到了包含4階撓性模態(tài)對象的Bang-Bang控制最優(yōu)解，結果顯示撓性衛(wèi)星實現靜止到靜止的快速機動，需要的時間與純剛體衛(wèi)星的相差不多。張洪華等[2]在脈沖序列控制的基礎上，提出了基于非約束模態(tài)方程的方波序列控制方法，并給出了使得機動時間最小的方波規(guī)劃方法，其優(yōu)點是控制輸出是方波而非脈沖，因而特別適合衛(wèi)星上的噴氣推力器控制模式，也適合飛輪控制模式。

Bang-Bang控制給出快速機動時間的下界，其他方法的機動時間都要更長，更多的研究人員著手研究光滑的次優(yōu)控制。一部分學者直接研究Bang-Bang控制的光滑逼近，增大一點機動時間得到連續(xù)平滑的控制函數，使得激發(fā)的過程撓性振動大幅度減少。Albassam等[3]研究了一類撓性衛(wèi)星的近似時間最優(yōu)控制，考慮柔化Bang-Bang控制曲線，即考慮輸入受限并且一階導數受限的剛體衛(wèi)星時間最優(yōu)控制，仿真表明該柔化曲線比Bang-Bang控制激發(fā)的撓性振動小。另一部分學者從能量最優(yōu)角度研究光滑控制，Hindle[4]等研究了撓性衛(wèi)星加權最小能量沖擊控制，得到了最優(yōu)控制解析解，作者研究了模型參數不確定性對于殘余振動的影響，通過引入新狀態(tài)變量描述不確定性，得到了對撓性頻率有魯棒性的控制解析解。

其他研究人員還考慮了非最優(yōu)控制方法，如變結構方法。王曉磊等[5]研究撓性衛(wèi)星在控制受限情況下，采用拋物線型滑動模態(tài)設計一種過渡過程性能很好的變結構控制器，并且分析了這種非線性滑動模態(tài)的穩(wěn)定性。周連文等[6]采用等速趨近率的滑模變結構控制策略進行大角度機動控制，并通過最優(yōu)控制理論設計彈性穩(wěn)態(tài)器，抑制由于剛體運動而激發(fā)的彈性振動，實現旋轉機動的同時，有效抑制彈性振動。這些方法在干擾和不確定性影響下的表現比最優(yōu)控制好，但只關心了整個系統閉環(huán)穩(wěn)定性，對快速機動過程中撓性振動的抑制和規(guī)避分析很少。

Chen[7]等研究前饋和反饋相結合的控制方法，其中前饋輸入是光滑的正弦函數形狀，仿真表明平滑前饋控制有助于抑制撓性振動，而反饋控制有助于應對外界干擾。正弦函數平滑且意義明確易于設計，常用于撓性結構參考軌跡設計。文獻[8]和[9]利用單正弦或者正弦序列設計撓性對象快速機動任務控制器，前者能方便地分析執(zhí)行機構的機動能力，后者則對撓性振動進行優(yōu)化。

本文采用正弦序列控制方法，考慮了對象控制受限和撓性振動受限，并優(yōu)化了快速機動后的殘余振動。該方法滿足一定意義上的振動指標最小，設計簡單，且受撓性參數影響小。作為對比，同時研究了其它兩種開環(huán)控制方法，包括光滑剛體控制[3]和Bang-Bang控制。數學仿真表明：作為開環(huán)控制，本文方法在快速機動任務中，產生的過程振動和殘余振動減少，適合工程應用。

1 研究對象和問題描述

研究的對象為中心剛體帶撓性附件類，如圖1所示，對象僅能繞中心軸轉動。中心剛體中的動量輪和角度敏感器為執(zhí)行機構和測量部件，撓性板通過剛撓連接體固定于中心剛體。假設在快速機動時，撓性附件為歐拉-伯努利懸臂梁，在任何時候形變很小可近似線性化。

圖1 試驗平臺示意圖

采用混合坐標法建立撓性衛(wèi)星模型，給出動力學方程

(1)

式中x為n+1維模態(tài)向量；u為控制輸入；y為角度輸出；M為慣量陣；K為剛度陣；G為輸入系數陣。

(2)

式中，z為非約束剛性模態(tài)；q=(q1,q2,…,qn)T為非約束撓性模態(tài)；φ0為剛體增益；φ1=(φ1,φ2,…,φn)T為模態(tài)增益陣。

進一步，上述模型寫成狀態(tài)空間形式

(3)

非約束模態(tài)；A,B,C為系統陣、輸入輸出陣。

2 正弦序列優(yōu)化控制

2.1 正弦序列介紹

利用已知的一類平滑控制函數，經過加權組合優(yōu)化得到控制解滿足“較小程度”激發(fā)撓性振動。這類平滑控制函數如果本身性質良好，那么得到的優(yōu)化控制解對撓性參數也不敏感。正弦函數就是這樣一類函數，其任意階導數有界，平滑特性非常好。

定義θf為機動角度，tf為機動時間，φ0為剛體增益。給定正弦系列控制輸入u(t)

(4)

其中參數滿足

那么基于方程(2)的剛體姿態(tài)響應x0(t)為

(5)

其中，

(6)

很顯然，控制輸入式(4)的參數條件使得剛體姿態(tài)x0(t)滿足邊界條件

也就是說，在控制律式(4)的作用下，衛(wèi)星的剛體部分將實現靜止到靜止的姿態(tài)定向任務。為了得到激發(fā)撓性振動較小的控制，需要對參數c1,c2,…,cm進行優(yōu)化。一般而言，只要m不是特別大，輸入控制式(4)都很平滑。特殊的當c1=1,c2=…=cm=0時，曲線為單正弦控制輸入，常在工程中使用。

2.2 殘余振動的優(yōu)化

下面討論正弦序列作用下的殘余振動。正弦序列由一組正弦函數構成，系數之和等于1保證剛體姿態(tài)到達目標。加權系數仍然還有自由度，利用這一點對殘余振動能量優(yōu)化。對于第i階撓性模態(tài)，在控制律(4)作用下，系統(3)在tf時刻殘余振動的表達式為

式中，λi只與第i階模態(tài)有關，與正弦序列無關。

殘余振動確定項，與機動時間無關。

2)ei=ωitf

殘余振動隨機項，與機動時間有關。

為降低殘余振動隨機項ei大小，只有調節(jié)機動時間tf。若此時衛(wèi)星撓性耦合很強，需要同時優(yōu)化多個模態(tài)殘余振動時，tf顯然不足以優(yōu)化所有模態(tài)。另外，殘余振動隨機項ei對撓性頻率ωi相當敏感，此時優(yōu)化ei的實用性不大。圖2可以驗證上述分析：對于快速機動tf=15s的撓性對象，若撓性頻率ωi為3rad/s左右，則隨機項ei誤差±20%置信頻率區(qū)間為2.99～3.08rad/s，即ωi容許頻率誤差僅為±1.5%左右。

圖2 隨機項敏感曲線

圖3 確定項敏感曲線

為了優(yōu)化殘余振動確定項fi，計算頻率ωi偏離時fi的變化

(8)

2.3 控制算法求解

由于控制輸入被限制為正弦序列，得到的將是次優(yōu)解，構造殘余振動指標，并對其中的確定項優(yōu)化，得到下面的確定時間次優(yōu)控制

次優(yōu)指標：

(9)

邊界條件：自動滿足

由式(2)及控制律(4)可求撓性模態(tài)

(10)

非線性規(guī)劃

minJ(C)

s.tU(C)≤um,Q(C)≤qmn

c1+c2+…+cm=1

(11)

其中:

3 數學仿真

3.1 正弦序列優(yōu)化控制

為了驗證方法的有效性，本節(jié)利用相關動力學模型[10]進行數學仿真。表1給出了該裝置的詳細參數，總慣量為I=20kg·m2，最大飽和力矩um=1N·m。保留了前5階撓性模態(tài)信息用于數學仿真，本文的控制方法優(yōu)化前2階主要模態(tài)。

表1 對象標稱參數

圖4 控制輸入曲線(15s機動90°)

仿真90°的快速機動，本文方法的機動時間是19s，而光滑剛體控制的機動時間是15s。2種方法最大控制力矩相同時，本文方法的機動時間稍長。本文2.3節(jié)中的參數取值為：Q1=Q2=I(單位陣)，正弦序列包括2個頻率分量，即m=2。圖5給出本組仿真的控制力矩曲線，2種控制方法始終都光滑連續(xù)，在初始時刻、末端時刻和開關時刻，兩種方法的曲線一階導數和二階導數都連續(xù)有界。尤其是本文的方法在整個機動過程都沒有切換存在，任意階導數有界，一致性非常好。

圖5 姿態(tài)輸出曲線(15s機動90°)

圖6和7給出了前兩階模態(tài)的振動情況。本文采用的控制方法整體性能很好，整個機動過程各階撓性模態(tài)振動都得到抑制。尤其是在最大撓性振幅處，本文采用的方法要明顯優(yōu)于光滑剛體控制。唯一不足的是，延長了機動時間，犧牲了執(zhí)行機構快速機動能力。

圖6 一階模態(tài)(15s機動90°)

圖7 二階模態(tài)(15s機動90°)

接下來的仿真是針對非標稱系統的。由于光滑剛體控制律的設計和撓性參數無關，所以該方法對撓性參數的魯棒性較好。2.2節(jié)分析了本文采用的控制方法受參數不確定性的影響，因此可以預見本文的控制方法對撓性參數的魯棒性也較好。

同3.1節(jié)的仿真(15s機動90°)，此時撓性對象頻率偏離標稱值20%，前兩階撓性頻率為ω1=2.3rad/s,ω2=3.9rad/s。圖8為本文控制方法的姿態(tài)曲線，從圖中可以看出，本文方法的魯棒性非常好，在撓性頻率偏差情況下，抑制殘余振動的能力受影響很小。

圖8 本文方法控制的姿態(tài)角(頻率偏差20%)

3.2 單頻率正弦控制

工程中為了設計簡單，有時采用單頻率控制(m=1)，這樣就失去了對殘余振動的優(yōu)化。本節(jié)的仿真我們將能看到，這種簡單控制方法抑制撓性振動不盡如人意，但是和一般的非光滑控制相比，激發(fā)的撓性振動還是要小很多。

圖9 控制力矩(15s機動90°)

圖10 一階模態(tài)(15s機動90°)

但作為光滑的控制方法，單頻率正弦控制激發(fā)的撓性振動要遠小于非光滑控制，為此本節(jié)還額外進行了多組仿真，這其中包括了Bang-Bang控制方法的仿真。并且引入過程振動能量這個概念，定義過程振動能量E：

表2給出了各種方法用于姿態(tài)快速機動90°時的過程振動能量E。單頻率正弦控制(一)和光滑剛體控制(二)的指標E很接近，這說明方法(一)和(二)激發(fā)撓性振動較小，是溫和、光滑的控制方法。最后還給出Bang-Bang控制(三)的過程能量E，要遠大于其他兩種控制方法。

4 結論

正弦序列優(yōu)化控制設計簡單，由幾個倍頻的正弦信號組成，整個機動過程中一致性好，激發(fā)的過程振動并不大。更重要的是，該方法能夠通過控制加權系數，減少殘余振動能量并且對模態(tài)信息不敏感，所以不失為一種較好的開環(huán)輸入控制。正弦序列優(yōu)化控制的缺點主要有2點：1)正弦序列個數受到對象撓性頻率的限制，在進行殘余振動優(yōu)化時，只能考慮少數幾個撓性模態(tài)，如果快速機動時間很短，正弦序列個數大大受限，優(yōu)化方法幾乎就不能使用；2)為了限制殘余振動，該方法需要更大的輸入力矩，會犧牲一定程度的快速機動能力。

表2 機動90°過程振動能量

把該方法和反饋控制相結合，為應對快速機動時外界的干擾力矩將是下一步的工作重點。

[1] Singh G, Kabamba P T, McClamroch N H. Planar Time-Optimal Rest-to-Rest Slewing Maneuvers of Flexible Spacecraft[J]. Journal of Guidance, Control and Dynamics, 1989, 12(1):71-81.

[2] 張洪華，黎康，趙宇.撓性衛(wèi)星姿態(tài)快速機動控制[J].中國空間科學技術，2005，25(1):53～59. (Zhang Honghua， Li Kang， Zhao Yu . Fast Maneuver Control of Flexible Satellite System[J].Chinese Space Science And Technology , 2005，25(1):53～59.)

[3] Albassam B A. Optimal Near-Minimum-Time Control Design for Flexible Structures[J]. Journal of Guidance, Control and Dynamics, 2002, 25(4): 618-625.

[4] Hindle T A, Singh T. Robust Minimum Power/Jerk Control of Maneuvering Structures[J]. Journal of Guidance, Control and Dynamics, 2001, 24(4):816-826.

[5] 王曉磊，吳宏鑫，李智斌.撓性航天器非線性滑動模態(tài)控制和試驗研究[J].宇航學報，2001，22(6):50-56. (Wang Xiaolei,Wu Hongxin,Li Zhibin. The Nonlinear Sliding Mode Control and Experiment Research for Flexible Spacecraft[J]. Journal of Astronautics，2001，22(6):50-56.)

[6] 周連文，周軍.撓性航天器大角度機動的變結構控制[J]. 航天控制，2003,21(3):48-52.(ZHOU Lian wen, ZHOU Jun. Variable Structure Control for Large Angle Maneuvering of Flexible Spacecraft[J].Aerospace Control, 2003, 21(3): 48-52.)

[7] Chen H J, Agrawal B N. Method of Slewing the Spacecraft to Minimize Settling Time[C]. AIAA Guidance, Navigation, and Control Conference and Exhibit. Monterey, California, 5-8 August, 2002.

[8] 宗紅，雷擁軍.基于控制力矩陀螺的衛(wèi)星姿態(tài)快速機動控制方法研究[J]．空間控制技術與應用，2007,33(5):11-16.(Zong Hong, Lei Yongjun. Research on Control of Satellite Attitude Rapidly Maneuvering with Control Moment Gyros[J].Aerospace Control and Application,2007,33(5):11-16.)

[9] Zhao H C, Chen D. Optimal Motion Planning for Flexible Space Robots[C]. International Conference on Robotics and Automation, Minneapolis, Minnesota, April 1996.

[10] 饒衛(wèi)東，徐李佳．單軸撓性衛(wèi)星快速機動試驗平臺[J]．空間控制技術與應用，2011，37(5):59-62.(Rao Weidong,Xu Lijia.A Simulation TestBed for Rapid Maneuver Experiment of Flexible Satellite[J].Aerospace Control and Application, 2011，37(5):59-62.)

Smooth Control for Quick Maneuver of Elastic Satellites

RAO Weidong1,2

1. Beijing Institute of Control Engineering, Beijing 100190, China 2. Science and Technology on Space Intelligent Control Laboratory, Beijing 100190, China

Thisarticleconcentratesonakindofsmoothcontrolwhichconsistsofseriesofsinusoids.Themethodcarriesoutquickslewingforelasticsatellitesandmeanwhilerestrainsvibrationsuppression.Analyticsolutionsofelasticvibrationareobtainedbyderivation.Rigidbodywillberesttorestautomaticwhencoefficientsofsinusoidssatisfysomeequations,andthenoptimizationisdonetowardsthecoefficientswhichleadtoaresidualvibrationsuboptimalsolutionbynonlinearprogramming.Relevantformulashowstheoptimizedsolutionsuppressesmuchvibrationandpossessesgreatrobustnesstowardselasticparameters.Numericsimulationiscarriedoutonaplantwithfiveelasticmodes.Theresultsindicatethemethod’sstrongabilityinrestrainingresidualvibration.Themethodhereisanopen-looponewhichcanbeusedasfeedforwardcontrol.

Elasticsatellites;Quickmaneuver;Feedforwardcontrol;Nonlinearprogramming

2012-06-13

饒衛(wèi)東(1985-)，男，江西人，博士研究生，主要研究方向為撓性衛(wèi)星控制。

V448.22

A

1006-3242(2013)05-0069-07