李 濤,朱春艷,陳必發(fā),張美艷,唐國安
(1. 復(fù)旦大學(xué) 航空航天系,上海 200433; 2. 上海宇航系統(tǒng)工程研究所,上海 201108)
航天器入軌后,太陽翼在驅(qū)動機(jī)構(gòu)的作用下轉(zhuǎn)動,實現(xiàn)對日定向,吸收太陽輻射能并轉(zhuǎn)化為電能,為航天器工作提供能量[1]。航天器大型化、輕質(zhì)化及低剛度的發(fā)展趨勢對其結(jié)構(gòu)設(shè)計及安全運(yùn)行提出了更高要求[2]。與航天器本體相比,太陽翼具有跨度大、剛度低的特點。由于太陽翼低階頻率ωn很小,衰減項e-ωnζt隨時間變化慢,使調(diào)姿、定向及變軌中形成的低頻振動很難自行衰減。這些振動不僅影響航天器的姿態(tài)穩(wěn)定、定位精度及機(jī)載設(shè)備的正常工作,嚴(yán)重時還會降低航天器使用壽命甚至導(dǎo)致結(jié)構(gòu)破壞,威脅航天器結(jié)構(gòu)的安全[3]。因此對太陽翼振動抑制開展研究尤為必要。
輸入整形技術(shù)是一種前饋控制方法,將原始輸入信號與依據(jù)系統(tǒng)動力學(xué)特性設(shè)計的特定脈沖整形序列卷積生成新的輸入信號,以減小結(jié)構(gòu)在特定頻率的振動[4-7]。文獻(xiàn)[8]提出了一種在太陽翼根部和表面加裝阻尼元件和彈性元件的被動控制方法,具有結(jié)構(gòu)簡單、易操作、成本低等優(yōu)點,是較為成熟的振動抑制方法。文獻(xiàn)[9]采用開關(guān)式噴射發(fā)動機(jī)、運(yùn)用不同的開關(guān)策略實現(xiàn)振動抑制,但作動需要消耗工質(zhì),不宜頻繁采用。近年來,以壓電材料為作動器的主動控制方法發(fā)展迅速,在航天器柔性結(jié)構(gòu)的振動抑制中得到廣泛應(yīng)用,控制效果顯著[10-15]。文獻(xiàn)[16]利用反作用飛輪抑制太陽翼振動,取得了不錯效果。文獻(xiàn)[17]提出了以太陽翼自身配備的驅(qū)動機(jī)構(gòu)為作動器的振動控制方法,但仿真模型基于帶配重的鋼尺,與實際太陽翼結(jié)構(gòu)相差較大,不能很好體現(xiàn)太陽翼動力學(xué)特性,且由于采用了較簡單的以驅(qū)動機(jī)構(gòu)轉(zhuǎn)速為輸出量的控制策略,控制結(jié)束后驅(qū)動機(jī)構(gòu)會發(fā)生角度偏移,使太陽翼無法復(fù)位且只開展了仿真分析,缺乏實驗支撐,振動抑制效果的驗證還顯不足。
在此基礎(chǔ)上,建立了太陽翼真實比例的有限元模型,提高仿真的真實性,并提出通過添加驅(qū)動機(jī)構(gòu)角度反饋的轉(zhuǎn)速項來解決驅(qū)動機(jī)構(gòu)角度的復(fù)位問題,開展了太陽翼縮比模型實驗,在控制算法中引入積分環(huán)節(jié),實現(xiàn)驅(qū)動機(jī)構(gòu)角度反饋轉(zhuǎn)速的等效作用,驗證了基于驅(qū)動機(jī)構(gòu)作動的振動控制方法的有效性。
哈勃望遠(yuǎn)鏡由望遠(yuǎn)鏡本體及展開的太陽翼組成,如圖1所示。實際中,航天器本體質(zhì)量遠(yuǎn)大于太陽翼質(zhì)量,為簡化考慮,在分析太陽翼振動時,將航天器本體視為靜止體。展開后太陽翼的簡化模型如圖2所示,建立隨驅(qū)動機(jī)構(gòu)運(yùn)動的非慣性坐標(biāo)系O-XYZ。
圖1 哈勃望遠(yuǎn)鏡Fig.1 Hubble space telescope
圖2 太陽翼簡化模型Fig.2 Simplified model of solar array
太陽翼在YOZ平面內(nèi)橫向彎曲振動主坐標(biāo)下的運(yùn)動方程[17]為
(1)
(2)
其中,Φ為根部固支的固定界面主模態(tài);M為太陽翼有限元模型的質(zhì)量矩陣;Ψ為繞Y軸轉(zhuǎn)動的剛體約束模態(tài)。
對于線性結(jié)構(gòu),太陽翼根部彎矩正比于相對位移,可表示為
Mb=κTΦξ
(3)
太陽翼振動主要由一階橫向彎曲振動構(gòu)成,式(1)可簡化為
(4)
式(4)右端為驅(qū)動機(jī)構(gòu)提供的輸入扭矩,控制律為
(5)
式中:ζa為通過反饋控制力矩添加的人工阻尼系數(shù),ζa>0。將式(5)代入式(4),進(jìn)一步可得
(6)
式中:ζa起到了加快振動衰減的作用,對小阻尼情況,ζa越大,太陽翼振動衰減速度越快。
(7)
式(7)為驅(qū)動機(jī)構(gòu)轉(zhuǎn)速與模態(tài)坐標(biāo)之間的關(guān)系,但在實際情況中,ξ1無法測得。根據(jù)式(3),在僅考慮一階模態(tài)的情況下,太陽翼根部的彎矩表示為
Mb=κTΦ1ξ1
(8)
根據(jù)式(7)和(8),驅(qū)動機(jī)構(gòu)轉(zhuǎn)速與太陽翼根部彎矩之間的反饋關(guān)系為
Ω=βMb
(9)
采用太陽翼有限元模型仿真驗證振動控制方法有效性,仿真模型如圖3所示。太陽翼長度為12 m,寬度為2.8 m,等效均質(zhì)板的厚度為3.5 mm,總質(zhì)量為317.5 kg,彈性模量E=7×1010Pa,泊松比ν=0.3。建立空間坐標(biāo)系O-XYZ,太陽翼在驅(qū)動機(jī)構(gòu)的作用下繞Y軸轉(zhuǎn)動。
圖3 太陽翼有限元模型Fig.3 Finite element model of solar array
太陽翼根部與航天器驅(qū)動機(jī)構(gòu)之間的連接簡化為1個扭簧,剛度ks=50 kN·m/rad。在進(jìn)行控制過程仿真時,可利用其計算太陽翼根部的彎矩
Mb=ks(θ1-θ0)
(10)
式中:θ0,θ1分別為驅(qū)動機(jī)構(gòu)和太陽翼根部的轉(zhuǎn)角。
太陽翼有限元模型采用Rayleigh比例阻尼,即
C=c1K+c2M
(11)
式中:c1=4.5×10-3;c2=2.7×10-3。太陽翼一、二階阻尼比為0.005。
仿真過程分為激振階段和抑振階段。在激振階段,模擬太陽翼調(diào)姿過程的驅(qū)動機(jī)構(gòu)轉(zhuǎn)動,采用梯形信號輸入,經(jīng)過勻加速、勻速、勻減速階段后停止。由于驅(qū)動機(jī)構(gòu)的轉(zhuǎn)速與太陽翼根部彎矩成正比,為保護(hù)驅(qū)動電機(jī),抑振階段設(shè)定在激振階段結(jié)束后、太陽翼根部彎矩接近于0時開始啟動。為此,在激振階段仿真時,有限元模型的節(jié)點位移和速度響應(yīng)計算結(jié)果被導(dǎo)出,作為振動抑制階段仿真的初始條件。
施加控制時驅(qū)動機(jī)構(gòu)的轉(zhuǎn)速為
(12)
式(12)可視為2個節(jié)點轉(zhuǎn)角之間的約束關(guān)系,用Nastran軟件仿真時可利用傳遞函數(shù)輸入功能來實現(xiàn),其標(biāo)準(zhǔn)形式為
(B0+B1p+B2p2)ud+
(13)
對式(12)進(jìn)行Laplace變換,并按式(13)的形式整理得
(1+η·s+0·s2)L(θ0)+
(-1+0·s+0·s2)L(θ1)=0
(14)
式中:η為控制的增益,η=(βks)-1。
在用Nastran軟件進(jìn)行振動控制仿真時,式(14)可用TF卡的形式表示,并添加到有限元模型[18]。但仿真結(jié)果表明,控制過程中驅(qū)動機(jī)構(gòu)會發(fā)生角度偏移,第226號節(jié)點位移如圖4所示。
圖4 不添加反饋轉(zhuǎn)速第226號節(jié)點位移Fig.4 Displacement of node 226 without feedback rotational speed
為使驅(qū)動機(jī)構(gòu)正確復(fù)位,對式(12)修正,整理得
(15)
式中:-λθ0為由驅(qū)動機(jī)構(gòu)轉(zhuǎn)動角度θ0反饋的轉(zhuǎn)速,當(dāng)θ0沒有復(fù)位時,提供反向轉(zhuǎn)速使其復(fù)位。式(15)也可用TF卡的形式添加至有限元模型,整理得
(λη+η·s+0·s2)L(θ0)+
(-1+0·s+0·s2)L(θ1)=0
(16)
結(jié)果表明:當(dāng)λ=0.1時,既能保證驅(qū)動機(jī)構(gòu)控制前后位置一致,又不影響控制效果。添加反饋轉(zhuǎn)速第226號節(jié)點位移如圖5所示。
圖5 添加反饋轉(zhuǎn)速第226號節(jié)點位移Fig.5 Displacement of node 226 with feedback rotational speed
圖6 不同增益下太陽翼根部彎矩衰減曲線Fig.6 Attenuation curve of bending moment at the root of solar array with different gains
以太陽翼縮比模型來驗證驅(qū)動機(jī)構(gòu)抑制振動方法的有效性。太陽翼縮比模型如圖7所示。由圖可見:太陽翼縮比模型由彈簧鋼片制成,中心用夾具夾持,模擬太陽翼調(diào)姿轉(zhuǎn)軸,底座通過螺釘與由步進(jìn)電機(jī)驅(qū)動的小轉(zhuǎn)臺連接,提供驅(qū)動力矩;太陽翼側(cè)向放置,彎曲振動與重力方向垂直;裝配4個配重,用于調(diào)節(jié)實驗?zāi)P偷墓逃蓄l率。經(jīng)實驗測量,太陽翼模型一階固有頻率為0.86 Hz。太陽翼縮比模型的參數(shù)見表1。
圖7 太陽翼縮比模型Fig.7 Scaled model of solar array
類型參數(shù)彈簧鋼片長度、高度、厚度/mm×mm×mm750×150×0.3內(nèi)測配重塊質(zhì)量m1/g220外側(cè)配重塊質(zhì)量m2/g122內(nèi)測配重塊至根部距離l1/mm137.5外測配重塊至根部距離l2/mm362.5
縮比模型實驗過程如圖8所示。由圖可知:應(yīng)變片將測得的根部應(yīng)變電壓輸入控制模塊,由設(shè)計的控制律進(jìn)行計算,輸出相應(yīng)的驅(qū)動電壓,以驅(qū)動步進(jìn)電機(jī)及轉(zhuǎn)臺,向太陽翼提供控制力矩,進(jìn)行振動控制。實驗中,應(yīng)變片貼在縮比模型根部,其電阻值R=350 Ω,采用全橋電路連接,橋路輸入電壓為5 V,應(yīng)變放大儀增益為1 000,濾波模式為低通3 Hz。在線性變形范圍內(nèi),太陽翼根部彎矩Mb正比于應(yīng)變放大儀測得的應(yīng)變電壓Uε,根部彎矩根據(jù)力及力臂參數(shù)計算,可標(biāo)定應(yīng)變電壓至根部彎矩的轉(zhuǎn)換系數(shù)為
圖8 縮比模型實驗流程圖Fig.8 Flow diagram of scaled model
(17)
實驗測得kV-M約為0.04 N·m/V。
根據(jù)式(9)和(17),轉(zhuǎn)臺轉(zhuǎn)速與應(yīng)變電壓之間的關(guān)系為
Ω=βkV-MUε
(18)
實驗中轉(zhuǎn)臺由步進(jìn)電機(jī)驅(qū)動。轉(zhuǎn)臺轉(zhuǎn)速與電機(jī)轉(zhuǎn)速的關(guān)系為
Ω=2πN/(60×i)
(19)
式中:N為步進(jìn)電機(jī)轉(zhuǎn)速,r/min;i為電機(jī)至轉(zhuǎn)臺的減速比,i=180。
實驗利用步進(jìn)電機(jī)細(xì)分驅(qū)動器,放大來自Arduino UNO單片機(jī)的脈沖信號,實現(xiàn)對步進(jìn)電機(jī)的精確平穩(wěn)控制。脈沖信號的頻率正比于驅(qū)動電壓,設(shè)定0~5 V的驅(qū)動電壓,電機(jī)便可以在0~285 r/min的轉(zhuǎn)速范圍內(nèi)轉(zhuǎn)動。電機(jī)的方向由映射至UNO單片機(jī)上數(shù)字信號的高低電平控制。驅(qū)動電壓Udri與電機(jī)轉(zhuǎn)速之間的關(guān)系為線性的,其表達(dá)式為
Ω=2π×285/(60×180)×(Vdri/5)
(20)
結(jié)合式(18)和(20),可得電機(jī)驅(qū)動電壓與應(yīng)變電壓的關(guān)系為
Udri=γUε
(21)
式中:增益為
γ=5×60×180/(2π×285)×βkV-M
(22)
根據(jù)理論分析可知,增益γ越大,振動抑制效果越好。增益的選取也應(yīng)考慮電機(jī)實際的轉(zhuǎn)速能力。當(dāng)太陽翼根部彎矩Mb最大時,轉(zhuǎn)臺轉(zhuǎn)速達(dá)到最大值,即應(yīng)變電壓最大時,驅(qū)動電壓達(dá)最大值。經(jīng)計算,增益最大值γmax=1.3。
應(yīng)變電壓采集過程中存在噪聲等干擾信號,因此需對輸入信號進(jìn)行濾波處理。選用巴特沃茲(Butterworth)帶通濾波器,該濾波器低通為0.75 Hz,高通為1.5 Hz。為了保護(hù)電機(jī)驅(qū)動器,需要對輸出的驅(qū)動電壓進(jìn)行飽和限制。完整的控制算法如圖9所示。與仿真一樣,在實驗中,驅(qū)動機(jī)構(gòu)存在角度偏移。為解決該問題,在控制算法中增加了積分環(huán)節(jié)。由經(jīng)驗可得,當(dāng)積分增益為-0.5時,轉(zhuǎn)臺位置在施加控制前后能保持一致,且不影響振動控制的效果。采用固定時間步長的算法,時間步長Δt=0.01 s。實驗激振采用梯形信號輸入,在1.14 s激振階段內(nèi)轉(zhuǎn)臺轉(zhuǎn)速Ω變化情況如圖10所示。
為驗證控制增益對振動抑制效果的影響,選取3組不同的增益(γ=0.3,0.7,1.3)進(jìn)行振動控制,實驗結(jié)果如圖12所示。
針對航天器太陽翼在調(diào)姿或變軌后低頻振動衰減速度緩慢的問題,本文以太陽翼自身驅(qū)動機(jī)構(gòu)為作動器,產(chǎn)生正比于系統(tǒng)主振動速度的反饋控制力矩,加快振動衰減速度。該方法無需增加額外機(jī)構(gòu),具有較強(qiáng)的可操作性。在有限元仿真時,通過添加驅(qū)動機(jī)構(gòu)角度反饋的轉(zhuǎn)速,并在實驗控制算法中引入積分環(huán)節(jié),解決了控制結(jié)束后驅(qū)動機(jī)構(gòu)發(fā)生角度偏移的問題。結(jié)果表明:基于驅(qū)動機(jī)構(gòu)作動的太陽翼振動控制方法能明顯減振;對于小阻尼情況,等效的人工阻尼系數(shù)越大,振動衰減的速率越快。本文方法能在保證太陽翼不偏離目標(biāo)位置的同時,有效抑制太陽翼的低頻振動。
圖9 控制程序框圖Fig.9 Block diagram of control program
圖10 轉(zhuǎn)臺轉(zhuǎn)速曲線Fig.10 Curve of rotational speed of spinning platform
圖11 太陽翼根部應(yīng)變電壓衰減曲線Fig.11 Attenuation curve of strain voltage at the root of solar array
圖12 不同增益下太陽翼根部應(yīng)變電壓衰減曲線Fig.12 Attenuation curve of strain voltage at the root of solar array with different gains
[1] 袁家軍. 衛(wèi)星結(jié)構(gòu)設(shè)計與分析[M]. 北京:中國宇航出版社, 2004.
[2] 黃文虎, 王心清, 張景繪, 等. 航天柔性結(jié)構(gòu)振動控制的若干新進(jìn)展[J]. 力學(xué)進(jìn)展, 1997, 27(1): 6-19.
[3] 李東旭. 撓性航天器結(jié)構(gòu)動力學(xué)[M]. 北京:科學(xué)出版社, 2010.
[4] 劉瑋, 郭其威, 朱春艷, 等. 柔性航天器艙體搬運(yùn)時輸入整形振動抑制方法和原理試驗[J]. 上海航天, 2015, 32(5): 31-34.
[5] 那帥, 朱春艷, 彭福軍, 等. 基于輸入整形技術(shù)的太陽翼調(diào)姿殘余振動抑制實驗研究[J]. 振動與沖擊, 2013, 32(7): 107-112.
[6] 朱春艷, 邵濟(jì)明, 那帥, 等. 太陽電池翼調(diào)姿后殘余振動抑制的整形器設(shè)計[J]. 振動與沖擊, 2012, 31(8): 176-180.
[7] 朱春艷, 紹濟(jì)明, 唐國安. 雙輸入整形器抑制兩軸轉(zhuǎn)動太陽翼調(diào)姿后的殘余振動[J]. 噪聲與振動控制, 2011, 31(6): 1-4.
[8] 許雲(yún)淞. 撓性衛(wèi)星姿態(tài)機(jī)動及太陽帆板振動抑制研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué), 2014.
[9] 任萱. 人造地球衛(wèi)星軌道力學(xué)[M]. 長沙: 國防科技大學(xué)出版社, 1988.
[10] MEYER J L, HARRINGTON W B, AGRAWAL B N, et al. Vibration suppression of a spacecraft flexible appendage using smart material [J]. Smart Materials and Structures, 1998, 7(1): 95-104.
[11] SONG G, AGRAWAL B N. Vibration suppression of flexible spacecraft during attitude control [J]. Acta Astronautica, 2001, 49(2): 73-83.
[12] HU Q L, MA G F. Spacecraft vibration suppression using variable structure output feedback control and smart materials [J]. Journal of Vibration and Acoustics, 2006,(2): 221.
[13] SALES T P, RADE D A, DE SOUZA L C G. Passive vibration control of flexible spacecraft using shunted piezoelectric transducers [J]. Aerospace Science and Technology, 2013, 29(1): 403-412.
[14] WU D, HUANG L, PAN B, et al. Experimental study and numerical simulation of active vibration control of a highly flexible beam using piezoelectric intelligent material [J]. Aerospace Science and Technology, 2014, 37: 10-19.
[15] HALIM D, LUO X, TRIVAILO P M. Decentralized vibration control of a multi-link flexible robotic manipulator using smart piezoelectric transducers[J]. Acta Astronautica, 2014, 104(1): 186-196.
[16] 陳必發(fā),呂亮亮,彭福軍,等.利用反作用飛輪抑制太陽電池陣振動的原理和方法研究[J].上海航天,2018,34(1):104-109.
[17] 李濤, 陳必發(fā), 呂亮亮,等. 基于驅(qū)動機(jī)構(gòu)作動的柔性附件振動控制原理與算法[J]. 載人航天, 2017, 23(4):482-486.