牛寶華 姚瀚宇 艾永強 王鋒斌 徐曉龍 賀爾銘
?
星載天線地面動平衡試驗仿真研究
牛寶華1姚瀚宇2艾永強1王鋒斌1徐曉龍1賀爾銘2
(1.西安空間無線電技術研究所,西安 710100;2.西北工業(yè)大學航空學院,西安 710072)
為了實現某型衛(wèi)星的高精度姿態(tài)控制,需要對其大型回轉天線探測頭部結構進行嚴格的動平衡設計仿真與試驗。結合星載天線研制需求,首先針對該型回轉天線設計了動平衡配平方案,通過試驗仿真及結果比較,甄選出動平衡試驗最佳配平方案并給出指導性建議;探索了在地面動平衡試驗時補償重力影響的幾種措施,驗證了加支撐桿消除重力影響的可行性。最終與動平衡試驗相結合,以很小的配重成本實現了衛(wèi)星天線的配平要求,保證了衛(wèi)星的按期發(fā)射。
大型回轉載荷;動平衡設計;補償重力
某型號衛(wèi)星的天線探測頭部為大型回轉結構,回轉部分質量在整星質量中占比很大。衛(wèi)星入軌后會以工作轉速持續(xù)旋轉進行掃描觀測,由于衛(wèi)星質量分布不均,衛(wèi)星旋轉會引入靜不平衡與動不平衡,產生不平衡力,使衛(wèi)星反射器發(fā)生形變,影響衛(wèi)星的飛行姿態(tài)控制。為實現衛(wèi)星的高精度姿態(tài)控制,需要對天線探測頭為主體的大型回轉載荷進行精確的動平衡設計仿真與動平衡試驗與配平。
動平衡的概念最早出現在轉子動力學領域,各種動平衡理論與實踐方法在航空發(fā)動機[1]、直升機旋翼[2]等領域得到了廣泛應用。對于衛(wèi)星,在太空微重力環(huán)境下,回轉載荷支承受力和潤滑狀態(tài)與地面大不相同,需要對其動平衡特性進行更精確的計算與控制[3~5],因此必須通過建立仿真衛(wèi)星的精確有限元模型[6,7],基于有限元模型進行動平衡仿真與試驗,提出配平方案,減小衛(wèi)星的動靜不平衡量,并盡量減小反射器最大形變量與支撐形變量。某型天線設計指標要求配平后靜不平衡量≤5kg·mm,動不平衡量≤2000kg·mm2。此外,在太空中衛(wèi)星處于失重狀態(tài),而在地面試驗中重力的影響不可避免,這會導致衛(wèi)星結構的形變,因此需要設法消除重力對配平結果的影響,使試驗結果與實際工作狀態(tài)一致。
通過有限元分析軟件MSC/PATRAN,建立了天線探測頭部結構展開兩種狀態(tài)下的有限元模型。其中高頻箱、反射器及展開機構、冷空反射器及支架、鎖緊釋放裝置均采用shell單元;高頻箱內部設備、饋源以及定標源采用質量相同的shell單元進行等效;高頻箱內波導管采用矩形截面的梁單元;各結構之間的螺栓連接采用MPC模擬;其他膠接面則通過合并重節(jié)點實現。最終建立了如圖1所示的有限元模型,有限元模型與原始物理模型各部分質量誤差很小,可以認為有限元模型十分準確。
圖1 有限元模型及位移云圖
在試驗中直接反應衛(wèi)星動平衡情況的指標就是衛(wèi)星反射面位移與支撐變形量,因此我們重點關注工作轉速下的衛(wèi)星形變結果。由有限元分析結果中反射器位移云圖可以獲得反射器的反射面最大位移。
由于結構有限元模型中,每個軸承支撐處的所有節(jié)點自由度是耦合在一起的,圓心處節(jié)點的合位移就是支撐變形量。由有限元分析軟件PATRAN導出上軸承支撐處支柱圓心處節(jié)點的合位移,得到25r/min轉速下的支撐形變量結果:反射器位移4.09mm,支撐變形量0.149 mm。
采用有限元網格法計算配平前載荷艙本身的靜平衡殘余量以及動平衡殘余量。有限元網格法計算公式如下所示:
剩余靜不平衡量
剩余動不平衡量
其中:m——轉子第個單元質量; (x,y,z)——轉子第個單元質心坐標。
根據初始動、靜不平衡量的計算公式,編寫MATLAB程序,利用PATRAN中的Mass Properties功能,將展開狀態(tài)模型各個單元的質量質心信息導入MATLAB中,計算得到天線展開模型的初始靜、動不平衡量分別為2.38238E+04kg·mm,1.55228194E+07kg·mm2,嚴重超過了設計指標要求。
總結SMAP衛(wèi)星以及XOVWM衛(wèi)星的配平經驗,主要采用整體配平法。整體配平中遵循雙面動平衡配平原理:對于每一個不平衡力,都可以將其分解到兩個配平面上,得到配平面上的等效不平衡力,在這兩個配平面分別加上大小相等、方向相反的平衡質量,即可消除該不平衡力。
除整體配平外,還可以嘗試分部件配平法,即將衛(wèi)星分成兩部分分別配平。而根據SMAP衛(wèi)星配平的經驗,也可嘗試僅對主反射器配平的方案。
根據配平原理與展開狀態(tài)有限元模型,在不考慮重力的前提下,選取不同配平面,可計算得到多種配平方案的配平結果,如表1所示。
表1 各方案動平衡配平詳細結果
對各個方案配平后的模型進行應力分析,得到工作轉速下的衛(wèi)星形變情況,如表2所示。
表2 25r/min各方案形變結果
對比表1、表2:
a. 由表1可知方案1、2、3、4的動靜不平衡量都滿足配平標準,且分部件配平與整體配平相比,只需更小的配平質量即得到更小的動靜不平衡量。
b. 由表1、表2可知,僅配平反射器的方案雖然能以很小的質量顯著減小反射面位移與支撐形變量,但配平后的動靜不平衡量過大,不符合配平指標。
c. 由表2可知,雙面法可最大程度地減小軸承的支反力,分部件方法雖然能更好地減小動靜不平衡量,但也會增大反射面位移與支撐形變量
綜合考慮配平質量、反射面最大位移、支撐變形量等方面的因素,在不考慮重力的前提下,最終選定方案2為最優(yōu)方案,方案4為次優(yōu)方案。建議在實際配平中采取整體雙平面配平,在反射器和高頻箱上以盡可能大的間距選取兩個配平面,配平效果最佳。
前文用于動平衡分析的有限元模型是理想模型,未考慮重力的影響。而在實際試驗中,重力必然會導致衛(wèi)星結構的形變,因此需要研究重力對動平衡的影響。
對按最優(yōu)方案2配平后的有限元模型研究,添加重力載荷,靜力分析得到模型各個單元的形變情況,重新進行動靜不平衡量分析。不考慮重力時,按方案2配平后的動不平衡量為1272.60kg·mm2,靜不平衡量為1.21kg·mm,都已經滿足配平標準。但加載重力后,動不平衡量變?yōu)?2342.12kg·mm2,靜不平衡量變?yōu)?8.47kg·mm,動靜不平衡量將嚴重惡化,不符合配平要求。
在地面試驗中必然存在重力,但實際工作狀態(tài)中處于失重狀態(tài),因此必須要設法消除重力的影響,使試驗結果與實際工作狀態(tài)一致。主要有兩種思路:
a. 先在考慮重力的前提下進行配平,并在地面試驗中驗證,在實際工作狀態(tài)中設法補償重力消失的影響,使試驗結果與實際工作狀態(tài)一致。
b. 在地面試驗中先設法消除重力的影響,在消除重力影響后進行配平與試驗,使試驗結果與實際工作狀態(tài)情況一致。
根據第一種思路,需先對衛(wèi)星模型在加載重力后進行配平,再設法補償重力消失的影響。加載重力后的配平方案與4.1中類似,在此不再贅述,重點關注如何補償重力消失影響,本文主要研究通過提高轉速方式補償重力消失的影響。
因為有限元模型近似關于平面=0對稱,所以重力和離心力作用下幾乎不產生方向位移。離心力和重力的最終作用效果都是使有限元模型反射器部分產生向和向位移。所以在太空中無重力狀態(tài)下,可以通過提高轉速來補償重力消失反射器的反射面部分向和向位移減小的影響。至于轉速提高多少才能最大程度地補償重力消失對反射面向和向位移的影響,并沒有直接計算公式,需要手動提高轉速一步步去試算。計算結果見表3。
由表3可知,在地面配平后,太空中為補償重力消失的影響,25r/min的工作轉速需要提高到76r/min。在工作狀態(tài)中將轉速提高數倍,明顯不符合實際,因此提高轉速消除重力影響只在理論上可行,實際空間應用十分困難。
表3 提高轉速結果
在地面動平衡試驗中,通過兩根碳纖維支撐桿將反射器支起,模擬太空中重力消失的情況,盡可能使反射器達到太空中重力消失后狀態(tài),經過靜力分析與地面試驗驗證后得到支撐桿數據如表4所示,支撐桿示意圖如2所示。
表4 支撐桿數據
圖2 支撐桿示意圖
4.3.1 帶支撐桿配平
因為加支撐桿后的重力模型與無重力模型相似,對加支撐桿后的重力模型配平時,使用無重力模型數據。每個支撐桿劃分100個單元,共200個單元。
添加配重位置選擇反射器給定位置、高頻箱上表面和高頻箱下表面。根據類似4.1的方法配平,共添加配重14.929kg,結果如表5所示。
由表5可知:高頻箱上表面,添加配重的質量很?。?.35E-16kg),可以忽略,實際添加配重的位置僅為反射器根部給定位置和高頻箱下表面。
具體配平位置如圖3所示。
表5 帶支撐桿配平結果
圖3 添加配重位置
4.3.2 補償太空中去掉支撐桿影響
在地面上對帶支撐桿模型配平可以達到平衡狀態(tài),但在太空中去掉支撐桿后,反射器又進入新的不平衡狀態(tài)。為了保證太空中去掉支撐桿后,整個轉子結構依然處于平衡狀態(tài),需要修正帶撐桿配平方案,在方案中給定的ABCD配平位置添加或減少配平質量,補償支撐桿去掉的影響,使去掉撐桿配平后的動靜不平衡量與帶支撐桿配平的動靜不平衡量誤差最小。
求解原理與過程與上文類似,結果如表6所示。
表6 帶支撐桿配平結果
由此可見,通過在給定位置添加少量配重,幾乎可以完全補償太空中去掉支撐桿對剩余動、靜不平衡量的影響,誤差大小數量級可以達到10E-3。
本文結合某型號衛(wèi)星研制需求,開展了針對星載天線的動平衡設計仿真與試驗,通過相對較小的配重質量極大地改善了產品的動平衡特性,最終衛(wèi)星天線的靜不平衡量和動不平衡量達到設計指標要求。通過研究重力對天線動平衡的影響,總結出一套補償重力對星載天線動不平衡量影響的試驗方法,為后續(xù)星載天線地面平衡試驗積累了寶貴經驗。該型號衛(wèi)星已成功發(fā)射,衛(wèi)星入軌后運行良好,其回轉天線指向精度滿足設計指標要求。
1 姜廣義,王德友,焦業(yè)秋. 大型發(fā)動機轉子本機平衡技術試驗研究[J]. 航空發(fā)動機,2008(1):19~22
2 周秋峰. 直升機旋翼動平衡優(yōu)化研究[D]. 南京:南京航空航天大學,2013
3 王朋朋,牛寶華,艾永強,等. 海洋二號衛(wèi)星微波輻射計的動平衡設計仿真與試驗[J]. 振動與沖擊,2016,35(9):1~6
4 薛孝補,朱海江,呂旺,等. 衛(wèi)星大型旋轉載荷動平衡控制應用研究[J]. 宇航學報,2012,33(8):1001~1005
5 王秋曉,汪飛,孔波. 一種改進的衛(wèi)星旋轉載荷動平衡的振動結構[J]. 振動與沖擊,2016,35(3):204~208
6 He Erming, Yin Geliang, Yaqi Hu. Random vibration analysis and structural modification of satellite-borne equipment. Applied Mechanics and Materials, 2014(455):310~313
7 Mobrem M, Kuehn S, Spier C, et al. Design and performance of Astromesh reflector onboard soil moisture active passive spacecraft. Aerospace Conference, 2012 IEEE, 1~10
Dynamic Balancing Design Simulation and Ground Test for Satellite Antenna
Niu Baohua1Yao Hanyu2Ai Yongqiang1Wang Fengbin1Xu Xiaolong1He Erming2
(1. China Academy of Space Technology(Xi’an), Xi’an 710100;2. School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072)
In order to realize high precision attitude control of the satellite, the dynamic balancing design simulation and test for large rotating payload of microwave radiometer should be conducted. At first, combining with the research and development of the satellite antenna, the dynamic balancing schemes of satellite antenna were designed, and balancing simulation was performed to provide the best balancing scheme and suggestions for ground balancing test. Then, several measures to compensate the influence of gravity during the ground dynamic balancing test were simulated, and the feasibility of adding support rods to eliminate the influence of gravity was verified. Finally, combined with the ground balancing test, antenna’s balance target was achieved with very small added weights, and this study ensured the satellite’s successful launch on schedule.
large rotating payload;dynamic balancing design;compensating gravity
牛寶華(1976),高級工程師,機械工程專業(yè);研究方向:航天技術領域星載天線研究設計。
2019-04-08