范震 郭翔鷹 張偉

（北京工業(yè)大學機電學院，機械結構非線性振動與強度北京市重點實驗室，北京 100124）

大型空間環(huán)型天線結構動力學分析*

范震 郭翔鷹 張偉

（北京工業(yè)大學機電學院，機械結構非線性振動與強度北京市重點實驗室，北京 100124）

本文對大型空間環(huán)型天線結構振動的基本形式進行研究，選擇邊界條件為伸展臂與桁架結構連接面的固定約束.網(wǎng)格劃分采用四面體單元與六面體單元結合的方法，桿件體結構規(guī)則采用六面體單元格，剛性連接處由于結構復雜采用四面體單元格.前四階模態(tài)分析結果表明天線存在扭轉振動、彎曲振動、橫搖振動、呼吸振動等四種基本振動形式.研究了材料參數(shù)和邊界條件的選取對振動模態(tài)和固有頻率的影響.所得結果能為結構非線性動力學分析中模態(tài)函數(shù)選取和共振形式的分析提供重要參考.施加簡諧激勵研究大型空間環(huán)型天線穩(wěn)態(tài)受迫振動響應，該結果對于預測結構的動力學特性以及避免共振具有重要的指導意義.

大型空間環(huán)型天線結構，扭轉運動，模態(tài)分析，材料參數(shù)，諧響應

引言

隨著航天技術的迅速發(fā)展，對于未來新型衛(wèi)星天線的要求也在不斷提高.大口徑，高精度，低質量，高強度已經(jīng)成為大型衛(wèi)星天線設計的必然趨勢，同時衛(wèi)星天線又面臨航天器運載空間及運載能力的限制.解決這一矛盾的關鍵技術在于天線結構型式的合理設計，使大型衛(wèi)星天線能被現(xiàn)有的運載工具送入軌道，同時滿足衛(wèi)星多功能、大容積、長壽命等功能的需要.目前大型空間衛(wèi)星普遍設計成可展結構，完成展開后鎖定來保持工作狀態(tài)［1-5］.為了國防技術以及國民經(jīng)濟的迫切需求，尤其是航空航天技術的發(fā)展需要，空間大型網(wǎng)架式可展天線的研制在許多國家都非常受重視，我國研究現(xiàn)狀處于起步階段，因此對大型環(huán)狀天線振動特性進行研究對保證天線穩(wěn)定工作具有重要意義［6-7］.

空間可展開桁架結構形式多樣，并且具有收縮體積小、剛度高等特點，可展開桁架結構可作為桁架式索網(wǎng)拋物面天線反射器的周邊支撐結構.無論是理論、試驗研究方面還是實際應用方面，桁架式可展開天線都是國際宇航界較為關注的熱點［8］. 1990年，Astro航空航天公司研制了周邊環(huán)形析架可展開天線，該天線代表了目前可展開天線的技術水平.國內(nèi)該天線研究現(xiàn)狀處于起步階段［9-18］.

本文利用FEM的方法對實際可展天線結構進行仿真，通過四面體六面體相結合的方式進行網(wǎng)格劃分，并選取適當?shù)倪吔鐥l件，進行模態(tài)分析得到了天線結構的諧振頻率和振型.研究邊界條件和材料參數(shù)對結構振型和固有頻率的影響.在非線性動力學建模過程中選取模態(tài)振型時，很難初步了解天線可能的振動形式，現(xiàn)在進行模態(tài)分析計算，對可能的振動形式進行分析，為以后非線性動力學建模過程中模態(tài)振型選擇提供重要依據(jù)，對提高非線性動力學計算精度具有重要指導意義.同時研究結果為后續(xù)天線振動模態(tài)的選取以及參數(shù)共振的比例關系的研究提供重要參考［23］.

1 大型空間環(huán)型桁架結構建模

ANSYS軟件是融合了結構、熱、流體等于一體的通用有限元分析軟件，其優(yōu)點之一就是能夠分析復雜多體結構.

在結構動力學中，系統(tǒng)的振動特性可以用模態(tài)來描述，表征模態(tài)的各階參數(shù)是振動系統(tǒng)的各階固有頻率、固有振型等，建立用模態(tài)參數(shù)表示的振動系統(tǒng)的運動方程并確定其模態(tài)參數(shù)的過程就是模態(tài)分析［19-21］.用有限元法對網(wǎng)架式天線結構進行仿真共三個部分.第一部分是前處理，包括建模、選取單元類型、定義材料參數(shù)和劃分網(wǎng)格.第二部分是求解模塊，ANSYS中提供了強大的結構分析模塊，可以用于結構的動態(tài)性能分析，包括模態(tài)分析.第三部分是后處理，通過后處理模塊，得到模態(tài)分析［22］.

空間環(huán)形桁架是由n個結構完全相同的平行四邊形單元組成，這些單元在相交處通過固定夾角的鉸鏈相連接成為一個封閉的多邊環(huán)型結構.其中平行四邊形單元結構由豎桿、橫桿、斜桿三種桿件組合而成.桿件為空心管狀結構，桿件直徑為20mm和18mm兩種，桿件壁厚度均為1.5mm，空心桿件結構可以有效減輕天線整體的質量，從而提高天線機動性與安全性.由于環(huán)型桁架天線鉸鏈連接處包含復雜結構，這種復雜結構對網(wǎng)格劃分收斂性影響較大，本文將天線桿件間的鉸鏈連接簡化為與桿件壁厚相等的剛性連接，從而使得模型獲得較好的收斂性.本文簡化后的天線有限元結構模型如圖1所示.

圖1 大型環(huán)型桁架天線有限元模型Fig.1 The FEM model of the large deployable truss antenna

1.1 模型尺寸與材料參數(shù)

環(huán)型桁架天線的結構尺寸參數(shù)如表1所示.

表1 天線整體尺寸參數(shù)Table 1 Dimension parameter of antenna

本文主要研究了結構尺寸與材料參數(shù)對環(huán)型桁架天線結構振動形態(tài)的影響，選取的材料為某種超強碳纖維纏繞，材料參數(shù)如表2所示.

表2 天線材料參數(shù)Table 2 Material parameter of antenna

1.2 邊界條件的選取

環(huán)型桁架天線結構由各桿件通過三桿節(jié)點和五桿節(jié)點連接而成.選取桁架結構邊界條件為同一豎桿兩端一個三桿節(jié)點和一個五桿節(jié)點處的固定約束，如圖1所示，A為三桿節(jié)點B為五桿節(jié)點，其他位置為自由約束條件.如圖2所示.

圖2 大型環(huán)型桁架天線邊界條件Fig.2 The boundary conditions of the large deployable truss antenna

1.3 網(wǎng)格的劃分

網(wǎng)格常用的劃分方式有四面體、六面體網(wǎng)格.其中四面體常應變單元具有簡單，計算速度快，邊界適應性強，易于自動劃分且便于由表面向實體內(nèi)部生成網(wǎng)格等特點.六面體單元由于變形特性好、計算精度高、相同精度要求所需網(wǎng)格單元數(shù)量較四面體少等諸多優(yōu)點被廣泛應用于三維有限元仿真等工程領域.

本文中環(huán)型桁架的桿件結構規(guī)則，為了提高精度采用六面體單元格.環(huán)型桁架結構剛性連接處由于結構復雜，需要采用邊界適用性更強的四面體單元格.因此本文網(wǎng)格劃分采用四面體單元與六面體單元結合的方法，兼顧兩者優(yōu)勢.從劃分過程中發(fā)現(xiàn)隨著單元數(shù)目的增加，即單元尺寸的縮小，解的近似程度將不斷改進，如果單元是滿足收斂性要求的，其近似解最后將收斂于精確解.劃分結果顯示，共有212069個單元格.其劃分情況如表3和圖3所示.

表3 天線網(wǎng)格劃分情況Table 3 Mesh status of antenna

圖3 大型環(huán)型桁架天線網(wǎng)格劃分Fig.3 Themeshing of the large deployable truss antenna

2 大型空間環(huán)型桁架結構仿真計算結果

模態(tài)分析模塊主要分析結構的振動特性，利用模態(tài)分析得到大型空間環(huán)型天線結構振動的基本模式，計算結構前四階振動頻率和振動形式.并通過改變邊界條件，研究不同邊界條件對振動形式和振動頻率的影響.

2.1 四種振動形式計算結果

采用1.2節(jié)邊界條件.第一階振型主要是剛體運動的彎曲變形.第二階振型主要是整體的周向扭轉振動形式.第三階振型主要是徑向呼吸振動形式.第四階振型是橫搖振動形式.

通過分析結果可以發(fā)現(xiàn)，振動的低階模態(tài)主要以彎曲、扭轉、呼吸運動為主.如圖4所示，第一階模態(tài)為整體的大范圍彎曲運動，最大位移出現(xiàn)在固定約束的最遠端，最大位移方向主要為軸向.如圖5所示，第二階模態(tài)為周向的大范圍扭轉運動，最大位移出現(xiàn)在固定約束的最遠端，最大位移方向主要為周向.第三階模態(tài)為徑向的呼吸運動，最大位移出現(xiàn)在固定約束的最遠端，最大位移方向主要為徑向.第四階模態(tài)為整體的橫搖運動，最大位移出現(xiàn)在固定約束的兩側端，在整體橫搖的同時，出現(xiàn)局部柔性變形.隨著模態(tài)階數(shù)和固有頻率的升高，天線振動出現(xiàn)高階柔性體變形.

圖4 大型環(huán)型桁架天線第一階振型Fig.4 The firstmode shape of the large deployable truss antenna

圖5 大型環(huán)型桁架天線第二階振型Fig.5 The second mode shape of the large deployable truss antenna

圖6 大型環(huán)型桁架天線第三階振型Fig.6 The third mode shape of the large deployable truss antenna

圖7 大型環(huán)型桁架天線第四階振型Fig.7 The fourth mode shape of the large deployable truss antenna

2.2 四種振動頻率計算結果

天線的各階模態(tài)對應的頻率分析結果如表4所示.

表4 天線前4階固有頻率Table 4 First four order frequency of antenna

由頻率分析可以看出，前三階振動的頻率低且頻率接近，天線振動可能存在幾種振動形式的運動耦合的情況，可以進一步研究幾種振動形式耦合的復雜非線性運動.

2.3 呼吸振動形式計算結果

改變邊界條件，在1.2節(jié)設置邊界條件的基礎上，增加豎桿在天線結構軸向的約束，計算天線前六階振動形式和振動頻率.結果顯示，前六階呼吸運動為具有振動峰值個數(shù)與振動階數(shù)相同的呼吸運動，如圖8～13所示.

圖8 大型環(huán)型桁架天線呼吸第一階振型Fig.8 The first breathingmode shape of the large deployable truss antenna

圖9 大型環(huán)型桁架天線呼吸第二階振型Fig.9 The second breathingmode shape of the large deployable truss antenna

圖10 大型環(huán)型桁架天線呼吸第三階振型Fig.10 The third breathingmode shape of the large deployable truss antenna

圖11 大型環(huán)型桁架天線呼吸第四階振型Fig.11 The fourth breathingmode shape of the large deployable truss antenna

可見邊界條件不同，結構發(fā)生振動形式不同.在此種邊界條件下，振動均為徑向的呼吸運動.振動最大峰值點出現(xiàn)規(guī)律變化.一階呼吸運動共有一個振動峰值點，位置在距離固定邊界條件處周向180度處.隨著階數(shù)增加振動峰值點個數(shù)也增加，n階呼吸運動共有n個振動峰值點，且振動峰值點周向均勻分布.

圖12 大型環(huán)型桁架天線呼吸第五階振型Fig.12 The fifth breathingmode shape of the large deployable truss antenna

圖13 大型環(huán)型桁架天線呼吸第六階振型Fig.13 The sixth breathingmode shape of the large deployable truss antenna

2.4 呼吸振動形式頻率計算結果

天線的呼吸振動形式的前6階模態(tài)對應的頻率分析結果如表5所示.

表5 天線前6階固有頻率Table 5 The first six order frequency of antenna

由頻率分析可以看出，前二階振動的頻率低且頻率接近.此種邊界條件的前四階振動頻率分別為：1.8668Hz，1.9712Hz，3.5555Hz，4.9482Hz.1.2節(jié)所設置邊界條件下前四階振動頻率為：1.7301Hz，1.8520Hz，1.9855Hz，3.2223Hz.可見，改變天線邊界條件增加豎桿在天線結構軸向的約束，使得相同階數(shù)振動形式的固有頻率升高.因此為提高計算精度，應當重視對邊界條件的合理選擇.

3 材料參數(shù)對振動頻率的影響分析

天線發(fā)射過程中，需要考慮發(fā)射衛(wèi)星天線的質量，因此會采用較小密度，減輕重量.本文采用的材料是強度大的碳纖維材料，為了分析采用高強度材料對固有頻率的影響，下面針對普通工程鋼材料的天線進行動力學分析，所采用工程鋼參數(shù)如表6所示，頻率計算結果如表7所示.同時為了研究不同彈性模量對前4階固有頻率的影響，計算了不同強度碳纖維材料天線的固有頻率.

表6 天線采用工程鋼的材料參數(shù)Table 6 Material parameters of antenna steel

表7 線不同材料參數(shù)的前4階固有頻率Table 7 The first four order frequency of the antenna with differentmaterial

表7中，不同彈性模量的碳纖維材料最低階固有頻率分別為1.7301Hz，1.8207Hz，1.9016Hz明顯高于傳統(tǒng)工程鋼材料的固有頻率1.3204Hz.不同彈性模量的碳纖維材料，隨著材料彈性模量的增加天線各階固有頻率緩慢增加.應當選擇合適的彈性模量的碳纖維材料.高強度材料使得天線低階振動頻率明顯增高，需考慮結構阻尼對降低固有頻率的影響，使振動頻率保持在合理范圍，保證天線工作精度.

4 大型空間環(huán)型天線結構諧響應分析

采用基于模態(tài)疊加法的諧響應分析，分析大型空間環(huán)型桁架結構的穩(wěn)態(tài)受迫振動.施加隨時間變化的簡諧外激勵作用點在周向距離固定邊界條件180度處的豎桿下端點.余弦形式外激勵力為：F =F0cos（Ωt），其中F0=100N.

圖14 大型環(huán)型桁架天線激勵作用點與觀測點Fig.14 The loading point and observation point of the large deployable truss antenna

圖15中可以看到，設置掃頻范圍為0～4Hz時，在頻率為1.7Hz時，在觀測點A的出現(xiàn)應力最大值為1.7Mpa，另提取1.7Hz時的整體天線結構整體的應力云圖，如圖16所示，天線的最大應力位于施加固定邊界條件的豎桿件右側第一根豎桿的上端點上，大小為162.36 Mpa.

圖15 0～4Hz頻率范圍觀測點應力-頻率曲線Fig.15 The stress-frequency curve of the observation point in the frequency range of 0-4Hz

圖16 1.7Hz頻率的應力云圖Fig.16 The stress contour of the antenna at the frequency of 1.7Hz

圖17中可以看到，設置掃頻范圍為0～4Hz時，在頻率為1.7Hz時，在觀測點A的出現(xiàn)位移最大值為7.37cm，另提取1.7Hz時天線結構整體的位移云圖，如圖18所示，天線的最大位移位于周向距離固定邊界條件180度處的豎桿下端點處.

圖17 0～4Hz頻率范圍觀測點位移-頻率曲線Fig.17 The displacement-frequency curve of the observation point in the frequency range of 0～4Hz

圖18 1.7Hz頻率的為位移云圖Fig.18 The displacement contour of the antenna at the frequency of 1.7Hz

在諧響應分析結果可以看到，應力響應的峰值與固有頻率相對應.該結果對于預測結構的持續(xù)性動力特性和避免共振產(chǎn)生具有重要指導意義.

5 小結

模態(tài)分析與頻率計算結果表明，大型環(huán)型桁架天線結構的低階振動模態(tài)多以扭轉運動與呼吸運動為主，這兩種振動屬于低頻振動，振動形式集中產(chǎn)生的高能量對振動響應的影響大.故后續(xù)研究應以環(huán)形天線的低階扭轉振動與呼吸振動為主要研究對象.結構有限元計算結果能否更精確反應結構在實際振動形式和頻率，除材料參數(shù)的選取外，還在于邊界條件的確定合理與否.同時大型環(huán)型桁架天線結構的模態(tài)分析與頻率計算為結構非線性動力學分析中模態(tài)函數(shù)選取和共振形式的分析都能提供重要參考，諧響應分析對于避免共振產(chǎn)生具有重要指導意義.

1 胡海巖，田強，張偉等.大型網(wǎng)架式可展開空間結構的非線性動力學與控制.力學進展，2013，43（4）：390～414 （Hu H Y，Tian Q，et al.Nonlinear dynamics and control of large deployable space structures composed of trusses and meshes.Advances in Mechanics，2013，43（4）：390～414（in Chinese））

2 Nayfeh A H，Hefzy M S.Effective constitutive relations for large repetitive frame-like structures.International Journal of Solids and Structures，1982，18（11）：975～987

3 Burgardt B，Cartraud P.Continuum modeling of beam like lattice trusses using averagingmethods.Computers&structures，1999，73（1）：267～279

4 Stephen N G，Zhang Y.Coupled tension-torsion vibration of repetitive beam-like structures.Journal of Sound and Vibration，2006，293（1）：253～265

6 閔桂榮.衛(wèi)星熱控制技術.北京：宇航出版社，1991 （Min G R.The Satellite thermal control technology.Beijing：Aerospace Press，1991（in Chinese））

7 閔桂榮，郭舜.航天器熱控制.北京：科學出版社，1998 （Min G R，Guo S.Spacecraft thermal control.Beijing：Science Press，1998（in Chinese））

8 安翔，馮剛，張鐸.大型空間析架結構熱分析的有限元方法.強度與環(huán)境，2000（2）：18～23（An X，F(xiàn)eng G，Zhang D.Finite-elementmode for the thermal analysis of largemain truss of space station′solar-array.Structure& Environment Engineering，2000（2）：18～23（in Chinese））

9 Tibert G.Deployable tensegrity structures for space applications［PhD Thesis］.Royal Institute of Technology，2002

10 Thomson M W.The AstroMesh deployable reflector.IUTAM-IASS Symposium on Deployable Structures：Theory and Applications.Springer Netherlands，2000：435～446

11 Thomson M W.The Astromesh deployable reflector.Antennas and Propagation Society International Symposium，IEEE，1999，3：1516～1519

12 Thomson M W.AstroMesh deployable reflectors for Kuand Ka-band commercial satellites.In：20th AIAA International Communication Satellite Systems Conference and Exhibit，2002：12～15

13 劉明治，高桂芳.空間可展開天線結構研究進展.宇航學報，2003，24（1）：82～87（Liu M Z，Gao G F.Advances in the study on structure for space deployable antenna.Journal of Astronautics，2003，24（1）：82～87（in Chinese））

14 Liu X N，Huang G L，Hu G K.Chiral effect in plane isotropic micropolar elasticity and its application to chiral lattices.Journal of the Mechanics and Physics of Solids，2012，60（11）：1907～1921

15 Luo GW，Zhang Y L，Xie JH.1：3 resonance bifurcation of a three-degree-of-freedom vibratory system with symmetrical rigid stops.Physics Letters A，2008，372（12）：2026～2031

16 Mahapatra D R.Equivalent constitutive model-based de-sign ofwave-absorbingmaterial system and controller.Journal of sound and vibration，2006，289（3）：509～528

17 張偉，胡海巖.非線性動力學理論與應用的新進展.北京：科學出版社，2009（Zhang W，Hu H Y.Advances in theory and application of nonlinear dynamics.Beijing：Science Press，2009（in Chinese））

18 趙孟良，關富玲.考慮摩擦的周邊桁架式可展天線展開動力學分析.空間科學學報，2006（26）：220～227（Zhao M N，Guan F L.Deployment dynamic analysis of circular truss deployable antenna with friction.Chinese Journal of Space Science，2006（26）：220～227（in Chinese））

19 陳國榮.有限單元法原理及應用.北京：科學出版社，2009：1～3（Chen G R.The principle and application of finite elementmethod.Bejing：Science Press，2009：1～3 （in Chinese））

20 王助成，邵敏.有限單元法基本原理和數(shù)值方法.北京：清華大學出版社，1997：443～81（Wang Z C，Shao M. The principle and numericalmethod of finite elementmethod.Beijing：Tsinghua University Press 1997：443～81（in Chinese））

21 應隆安，陳傳森.有限元理論與方法.北京：科學出版社，2009：688～700（Ying L A，Chen CS.The finite element theory and method.Beijing：Science Press，2009：688～700（in Chinese））

22 張洪信，吳俊飛，沈孝芹.有限元基礎理論與ANSYS應用.北京：機械工業(yè)出版社，2006：1～3（Zhang H X，Wu JF，Shen X Q.The finite element theory and application of ANSYS.Beijing：Machinery Industry Press，2006：1～3 （in Chinese））

23 胡海良，張偉.兩自由度非線性系統(tǒng)1：3內(nèi)共振的漸進攝動分析.動力學與控制學報，2015，13：97～100（Hu H L，Zhang W.Two degree of freedom nonlinear system asymptotically 1：3 internal resonance perturbation analysis. Journal of Dynamics and Control，2015，13：97～100（in Chinese））

DYNAM ICS ANALYSISOF LARGE DEPLOYABLE MESH ANTENNA STRUCTURE*

Fan Zhen Guo Xiangying ZhangWei

（Beijing Key Laboratory of Nonlinear Vibrations and Strength of Mechanical Structures，College of Mechanical Engineering，Beijing University of Technology，Beijing 100124，China）

This paper focuses on the fundamental vibration mode of the large deployable antenna structure. Here，the boundary conditions are considered as fixed Supportaccording to the actual contact surface between the extended arm and truss.The gridmesh for the large deployable antenna structure is combined by both tetrahedral element and hexahedral element，where the hexahedral element is employed in the rod parts of the large deployable antenna structure，while tetrahedral element is in the joint parts of the large deployable antenna structure.The analysis results of the first four ordermodels show that the four basic vibration forms of the antenna contain the twistmotion，the bending motion，the roll motion and the breathing motion.These four vibration modes are important for the responses of the large deployable antenna，since they have the similar frequencies.It also analyzes the influence ofmaterial parameters and boundary conditions on themode shape and frequencies. The results also provide the important reference to the selection of modal function and the nonlinear dynamic analysis of resonance parameters for the large deployable antenna.The study of the forced vibration through the harmonic response analysis is alsomeaningful to predict the dynamic behaviour and avoid resonance.

large deployable mesh antenna structure，twisting motion，modal analysis，material parameters，harmonic response analysis Received 17 April 2015，revised 17 May 2015.

E-mail：sandyzhang0@yahoo.com

10.6052/1672-6553-2015-038

*The project supported by the National Natural Science Foundation of China（11202009，11290152）

2015-04-17收到第1稿，2015-05-17收到修改稿.

*國家自然科學基金資助項目（11202009，11290152）

E-mail：eagle2008guo@yeah.net