劉 陽(yáng)，李團(tuán)結(jié)，陳聰聰，唐雅瓊

(西安電子科技大學(xué) 機(jī)電工程學(xué)院，陜西 西安 710071)

在衛(wèi)星通信、對(duì)地觀測(cè)、深空探測(cè)等航天領(lǐng)域，空間可展開天線起著重要的作用，其應(yīng)用十分廣泛。隨著空間技術(shù)的發(fā)展，許多科學(xué)應(yīng)用領(lǐng)域?qū)Υ罂趶胶统罂趶娇臻g可展天線的需求越來(lái)越高。由于受火箭載荷及空間的限制，口徑大、質(zhì)量輕、收縮比高的可展開天線成為各航天大國(guó)的研究重點(diǎn)[1-3]，為此各國(guó)學(xué)者提出了形式眾多的空間可展開天線結(jié)構(gòu)概念，環(huán)柱天線是其中最具潛力的結(jié)構(gòu)形式之一[4]。

可展開天線的結(jié)構(gòu)形式多樣，對(duì)于大口徑天線，最常用的反射面結(jié)構(gòu)為網(wǎng)狀反射面。網(wǎng)狀反射面具有曲面成型精度高、收攏體積小、天線的展收比大等優(yōu)點(diǎn)，在衛(wèi)星通信及遙感等領(lǐng)域得到了廣泛應(yīng)用，受到各航天大國(guó)普遍關(guān)注，具有制作大口徑和超大口徑天線的前景[5-7]。文中所研究的環(huán)柱天線是可展開索網(wǎng)天線的一種重要形式，美國(guó)國(guó)家航空航天局(NASA)的蘭利研究中心(Langley Reseacherch Center)與Harris公司在1980年前后共同提出環(huán)柱天線的概念[8]，它是一種預(yù)拉壓結(jié)構(gòu)，主要分為周邊圓環(huán)、中心伸展立柱、索網(wǎng)系統(tǒng)和金屬絲網(wǎng)反射面幾部分，當(dāng)天線收攏時(shí)，周邊圓環(huán)在電機(jī)的驅(qū)動(dòng)下向中心收攏，收攏體積可以很小[9-10]。為了驗(yàn)證這種方案的可行性，美國(guó)國(guó)家航空航天局制作了口徑為15 m的環(huán)柱天線樣機(jī)[11-12]，其收攏直徑為0.9 m，收攏高度為2.7 m，天線總重為291 kg，并對(duì)其性能進(jìn)行了測(cè)試[13]，這類天線的展開結(jié)構(gòu)較為簡(jiǎn)單，但索網(wǎng)系統(tǒng)復(fù)雜，網(wǎng)面設(shè)計(jì)和調(diào)整難度比較大[14]。

為此，筆者根據(jù)環(huán)柱天線結(jié)構(gòu)組成，確定了環(huán)柱天線索網(wǎng)的幾何構(gòu)型，采用分步求解的方式計(jì)算其預(yù)應(yīng)力分布，并修正節(jié)點(diǎn)坐標(biāo)，得到環(huán)柱天線的最終節(jié)點(diǎn)坐標(biāo)和預(yù)應(yīng)力分布。最后通過(guò)有限元仿真驗(yàn)證了所述索網(wǎng)結(jié)構(gòu)設(shè)計(jì)與分析方法的有效性。

1 索網(wǎng)結(jié)構(gòu)幾何設(shè)計(jì)

環(huán)柱可展開天線的主要結(jié)構(gòu)如圖1和圖2所示。主要包含周邊圓環(huán)、中心伸展立柱、索網(wǎng)系統(tǒng)和金屬絲網(wǎng)反射面。

圖1 環(huán)柱天線三維圖

圖2 環(huán)柱天線結(jié)構(gòu)圖

天線的中心立柱由可伸縮的套管組成，分為中心轂和上下立柱。周邊圓環(huán)是由48節(jié)桿件鉸接而成，通過(guò)上下各48根環(huán)支撐索分別與上下立柱端點(diǎn)相連，給周邊圓環(huán)提供支撐，天線網(wǎng)面固定在周邊圓環(huán)與下立柱之間，起著支撐和定位周邊圓環(huán)的作用，前網(wǎng)面采用三角形拓?fù)?，通過(guò)豎向和斜向張緊索與背索相連，使前網(wǎng)面的節(jié)點(diǎn)落在理想拋物面上，網(wǎng)面的背后連接著調(diào)整索，用以形面精度調(diào)整。這些索是由高剛度低熱膨脹系數(shù)的石墨材料組成的，給金屬絲網(wǎng)提供了一個(gè)非常穩(wěn)定的支撐結(jié)構(gòu)。饋源支架是一個(gè)獨(dú)立的結(jié)構(gòu)，用以支撐饋源陣列。收攏時(shí)，周邊環(huán)可以折疊起來(lái)收縮到中心柱周圍一個(gè)很小的體積內(nèi)。

圖3 1/2片索網(wǎng)結(jié)構(gòu)示意圖

文中的重點(diǎn)是研究環(huán)柱天線的索網(wǎng)系統(tǒng)設(shè)計(jì)，對(duì)于這種正饋旋轉(zhuǎn)拋物面可展開天線的網(wǎng)面設(shè)計(jì)，可以從單片索網(wǎng)入手，首先確定單片索網(wǎng)節(jié)點(diǎn)坐標(biāo)與拓?fù)?，然后將單片索網(wǎng)進(jìn)行坐標(biāo)旋轉(zhuǎn)得到完整索網(wǎng)，完成其幾何設(shè)計(jì)。整個(gè)索網(wǎng)分為48片。1/2片索網(wǎng)結(jié)構(gòu)如圖3所示。

1.1 前網(wǎng)面設(shè)計(jì)

前網(wǎng)面是一焦距為60 m的拋物面，以拋物面最低點(diǎn)為坐標(biāo)原點(diǎn)，前網(wǎng)面的節(jié)點(diǎn)位置與編號(hào)如圖4所示，前網(wǎng)面內(nèi)圈和外圈半徑分別為8 m和48 m，單片網(wǎng)面對(duì)應(yīng)圓心角為7.5°。為了使索網(wǎng)在靜力平衡時(shí)張力分布均勻，應(yīng)使前網(wǎng)面節(jié)點(diǎn)盡量呈對(duì)稱分布。將邊界拉索AB沿弧長(zhǎng)方向均分為27份得到AB間節(jié)點(diǎn)坐標(biāo)，其他節(jié)點(diǎn)坐標(biāo)由邊界拉索AB上節(jié)點(diǎn)旋轉(zhuǎn)相應(yīng)角度得到。引入垂跨比調(diào)整環(huán)連接索節(jié)點(diǎn)位置，用以降低環(huán)連接索對(duì)周邊圓環(huán)的拉力，減小周邊圓環(huán)變形。

圖4 單片索網(wǎng)結(jié)構(gòu)前網(wǎng)面示意圖

1.2 背網(wǎng)面節(jié)點(diǎn)坐標(biāo)設(shè)計(jì)

圖5 背網(wǎng)面結(jié)構(gòu)示意圖

背網(wǎng)面節(jié)點(diǎn)分布如圖5所示，調(diào)整索連接懸鏈線索與下立柱端點(diǎn)，調(diào)整索與x軸的夾角分別為68°、46°、36.2°和32°，豎向張緊索連接邊界索與懸鏈線索，使前網(wǎng)面位于理想拋物面上，方向沿拋物線在該點(diǎn)處的法線方向。

綜上所述，可以得到單片索網(wǎng)的節(jié)點(diǎn)坐標(biāo)與拓?fù)洌炀€反射面為正饋拋物面，由對(duì)稱特性將單片索網(wǎng)旋轉(zhuǎn)陣列即可得到索網(wǎng)整體坐標(biāo)與拓?fù)洹?/p>

2 索網(wǎng)結(jié)構(gòu)預(yù)應(yīng)力分析

索網(wǎng)天線的形面精度由前網(wǎng)面節(jié)點(diǎn)位置決定，而索網(wǎng)結(jié)構(gòu)只有在施加預(yù)應(yīng)力后才具有剛度，天線在這組預(yù)應(yīng)力的作用下保持平衡且前索網(wǎng)節(jié)點(diǎn)均位于理想拋物面上，因此索網(wǎng)結(jié)構(gòu)的預(yù)應(yīng)力設(shè)計(jì)是否合理直接影響著天線的形面精度。索網(wǎng)預(yù)應(yīng)力求解問(wèn)題中，常用的方法有力密度法、有限元法和動(dòng)力松弛法等，文中采用力密度法求解環(huán)柱天線索網(wǎng)預(yù)應(yīng)力。

2.1 力密度法基本原理

力密度法是一種常用的索網(wǎng)結(jié)構(gòu)找形方法，最初由Schek等提出[15]，力密度定義為索段的張力和索段長(zhǎng)度的比值。通過(guò)力密度的概念對(duì)索網(wǎng)結(jié)構(gòu)力平衡方程進(jìn)行變形，解決了原來(lái)求解非線性方程組較復(fù)雜的問(wèn)題，極大地降低了求解的難度，更加有利于運(yùn)算[16]。在力密度法中，網(wǎng)面的平衡位置與力密度值是相對(duì)應(yīng)的，只要已知一組力密度，對(duì)應(yīng)的網(wǎng)面節(jié)點(diǎn)坐標(biāo)就能夠直接求出。

對(duì)索網(wǎng)結(jié)構(gòu)列節(jié)點(diǎn)力平衡方程，以x方向?yàn)槔?，?jié)點(diǎn)i的平衡方程如下所示：

(1)

其中ci為與節(jié)點(diǎn)i相連的所有點(diǎn)，節(jié)點(diǎn)i、j構(gòu)成的索段稱為索段ij，Tij、Lij分別為索段ij的張力和長(zhǎng)度，Pix為節(jié)點(diǎn)i在x方向上所受的合力。力密度qij=Tij/Lij，則式(1)變換成

(2)

對(duì)環(huán)柱天線所有節(jié)點(diǎn)列平衡方程，如下所示：

CTQCx=Px,

(3)

其中C為索網(wǎng)結(jié)構(gòu)的拓?fù)渚仃?，維數(shù)為m×n，m、n分別為索段和節(jié)點(diǎn)的數(shù)目；Q為索段力密度矩陣，維數(shù)為m×m，它是一個(gè)對(duì)角陣，對(duì)角線上每一個(gè)數(shù)值分別對(duì)應(yīng)相應(yīng)索段的力密度。索網(wǎng)節(jié)點(diǎn)可分為自由節(jié)點(diǎn)和邊界固定節(jié)點(diǎn)兩類，數(shù)量分別為nf和ng，因此可以將拓?fù)渚仃嘋拆分為Cf和Cg，分別對(duì)應(yīng)自由節(jié)點(diǎn)和邊界固定節(jié)點(diǎn)；同理，節(jié)點(diǎn)的x坐標(biāo)也可以分為xf和xg，分別對(duì)應(yīng)自由節(jié)點(diǎn)和邊界固定節(jié)點(diǎn)的x坐標(biāo)。索網(wǎng)平衡時(shí)，自由節(jié)點(diǎn)所受的合力為零，邊界固定節(jié)點(diǎn)所受的合力為支撐框架給索網(wǎng)的支反力，設(shè)其大小為Pxg，則式(3)可變換為

(4)

(5)

在索網(wǎng)力密度已知的情況下，求解式(5)可以解出所有自由節(jié)點(diǎn)的x坐標(biāo)。運(yùn)用同樣的方法，可以計(jì)算出它們y、z方向的坐標(biāo)。通過(guò)此方法，給定一組力密度，無(wú)需迭代，可以快速求解與之對(duì)應(yīng)的自由節(jié)點(diǎn)坐標(biāo)。文中索網(wǎng)結(jié)構(gòu)的幾何拓?fù)湟阎?，需要尋找與之對(duì)應(yīng)的預(yù)應(yīng)力，所以需要對(duì)力密度法進(jìn)行改進(jìn)。首先假定一組初始預(yù)應(yīng)力，通過(guò)力密度法求解出與之對(duì)應(yīng)的平衡狀態(tài)下的自由節(jié)點(diǎn)坐標(biāo)，通過(guò)節(jié)點(diǎn)坐標(biāo)計(jì)算得到當(dāng)前索段長(zhǎng)度。為了得到新的力密度，可以通過(guò)下式進(jìn)行迭代[17]：

qi+1=qili/l0，

(6)

式中,li為通過(guò)平衡后坐標(biāo)計(jì)算得到的當(dāng)前索段長(zhǎng)度，l0為給定的初始索段長(zhǎng)度。

2.2 預(yù)應(yīng)力求解

進(jìn)行索網(wǎng)預(yù)應(yīng)力求解時(shí)，將周邊圓環(huán)結(jié)構(gòu)和中心立柱作為剛性結(jié)構(gòu)，不產(chǎn)生變形，相關(guān)節(jié)點(diǎn)固定。整個(gè)索網(wǎng)結(jié)構(gòu)共有346×48=16 608根索，每一根索對(duì)應(yīng)一個(gè)力密度值，如果對(duì)索網(wǎng)整體進(jìn)行力密度迭代，則變量太多，求解較為復(fù)雜，甚至有可能找不到解。為此，將索網(wǎng)拆分為前、后索網(wǎng)兩部分分步求解。前索網(wǎng)包含前網(wǎng)面、張緊索和懸鏈線索，其余索為后索網(wǎng)。首先由力密度法計(jì)算前索網(wǎng)的預(yù)應(yīng)力，然后根據(jù)求得的預(yù)應(yīng)力由節(jié)點(diǎn)力平衡計(jì)算剩余索預(yù)應(yīng)力?？紤]到旋轉(zhuǎn)拋物面索網(wǎng)結(jié)構(gòu)的對(duì)稱性，此處僅需對(duì)其2/48分塊的索網(wǎng)進(jìn)行預(yù)應(yīng)力求解，圖6為索網(wǎng)結(jié)構(gòu)2/48分塊的示意圖。

圖6 前索網(wǎng)結(jié)構(gòu)示意圖

2.2.1 前索網(wǎng)預(yù)應(yīng)力求解

圖7 修正節(jié)點(diǎn)示意圖

前索網(wǎng)的結(jié)構(gòu)如圖6所示，對(duì)索網(wǎng)結(jié)構(gòu)進(jìn)行預(yù)應(yīng)力優(yōu)化設(shè)計(jì)時(shí)，將前索網(wǎng)的邊界節(jié)點(diǎn)作為固定節(jié)點(diǎn)處理，固定節(jié)點(diǎn)的位置如圖6中圓點(diǎn)標(biāo)記所示。上下兩個(gè)分塊關(guān)于豎直平面P3對(duì)稱，上分塊與下分塊自身分別關(guān)于豎直平面P1和P2對(duì)稱。為了進(jìn)一步減少變量，同時(shí)保證由分塊旋轉(zhuǎn)得到整體時(shí)索力能夠平衡，選取P1和P2之間的索力密度為設(shè)計(jì)變量，其余索力密度由對(duì)稱特性直接賦值。

圖7所示為前網(wǎng)面中需要修正的節(jié)點(diǎn)位置，受到結(jié)構(gòu)的限制，這些節(jié)點(diǎn)背后并沒(méi)有連接張緊索，故結(jié)構(gòu)平衡時(shí)會(huì)偏離拋物面。由于這些節(jié)點(diǎn)所在的三角形面片尺寸明顯小于其他面片，所以節(jié)點(diǎn)的偏離帶來(lái)的誤差可以歸于面片擬合誤差。在預(yù)應(yīng)力設(shè)計(jì)階段，為了使索網(wǎng)結(jié)構(gòu)更容易達(dá)到穩(wěn)定平衡狀態(tài)，在迭代過(guò)程中，根據(jù)有限元軟件計(jì)算結(jié)果更新索網(wǎng)中需要修正節(jié)點(diǎn)的坐標(biāo)。

為了索網(wǎng)結(jié)構(gòu)的預(yù)應(yīng)力分布更加合理，懸鏈線索的節(jié)點(diǎn)坐標(biāo)也需要進(jìn)行調(diào)整，采用的方法與前網(wǎng)面中節(jié)點(diǎn)修正方法一致。

綜上前索網(wǎng)預(yù)應(yīng)力求解流程如圖8所示，求解程序編制的步驟如下：

(1)給定所有索段初始力密度q0，單步迭代次數(shù)n；

(2)求解式(5)得到第i次節(jié)點(diǎn)平衡坐標(biāo)(初次計(jì)算時(shí)i=0)；

(3)根據(jù)平衡后節(jié)點(diǎn)坐標(biāo)更新索段力密度qi+1=ti/l0，并給對(duì)稱索力密度賦值；

(4)解式(5)得到第i+1次節(jié)點(diǎn)平衡坐標(biāo)；

(5)判斷迭代次數(shù)是否達(dá)到單步迭代次數(shù)n。若達(dá)到，則轉(zhuǎn)入(6)，若未達(dá)到，則i=i+1，重復(fù)(2)～(4)；

(6)將所求得的索力代入ANSYS軟件中進(jìn)行靜力分析；

(7)判斷當(dāng)前索網(wǎng)節(jié)點(diǎn)是否滿足精度要求。若滿足，則輸出此時(shí)索段力密度，退出計(jì)算；若不滿足，則根據(jù)ANSYS軟件計(jì)算結(jié)果修正無(wú)張緊索和懸鏈線索節(jié)點(diǎn)坐標(biāo)，重復(fù)(2)～(6)。

2.2.2 后索網(wǎng)的預(yù)應(yīng)力求解

圖8 前索網(wǎng)預(yù)應(yīng)力求解流程圖

如圖5所示，前索網(wǎng)的預(yù)應(yīng)力已經(jīng)計(jì)算得到，針對(duì)剩余索段列節(jié)點(diǎn)力平衡方程，只需給定其中一根索的預(yù)應(yīng)力，即可計(jì)算出其余索的預(yù)應(yīng)力。文中指定下輪轂索力密度值，求解平衡方程得到剩余索段力密度。

2.3 仿真驗(yàn)證

利用上述方法計(jì)算得到所有索段預(yù)應(yīng)力。為了驗(yàn)證文中方法的有效性，將分析所得的預(yù)應(yīng)力作為初始應(yīng)變導(dǎo)入索網(wǎng)結(jié)構(gòu)有限元模型中，周邊圓環(huán)和中心立柱作為剛性結(jié)構(gòu)處理，不考慮重力的影響，進(jìn)行非線性有限元分析，分析結(jié)果如表1所示。

表1 仿真分析結(jié)果

圖9和圖10分別為位移云圖和索力分布云圖?？梢钥吹阶畲蟮墓?jié)點(diǎn)位移為1.93×10-7m。該節(jié)點(diǎn)位移已經(jīng)很小，可以認(rèn)為有限元仿真中出現(xiàn)的節(jié)點(diǎn)位移來(lái)自于預(yù)應(yīng)力計(jì)算過(guò)程中的截?cái)嗾`差。也就是說(shuō)，通過(guò)該方法得到的預(yù)應(yīng)力是正確的，且具有較高的精度。

圖9 位移云圖

圖10 索力分布云圖

3 結(jié)束語(yǔ)

環(huán)柱天線具有構(gòu)建大口徑和超大口徑天線的潛力。文中根據(jù)環(huán)柱天線結(jié)構(gòu)特點(diǎn)確定了環(huán)柱天線索網(wǎng)結(jié)構(gòu)的成形方法，將索網(wǎng)拆分為前后索網(wǎng)兩部分分別進(jìn)行預(yù)應(yīng)力求解。對(duì)于前索網(wǎng)，根據(jù)對(duì)稱特性進(jìn)行變量歸并，采用力密度法求解預(yù)應(yīng)力，并修正節(jié)點(diǎn)坐標(biāo)；由節(jié)點(diǎn)力平衡求解后索網(wǎng)預(yù)應(yīng)力。最后通過(guò)仿真分析，驗(yàn)證了文中所提方法的可行性與有效性，為今后復(fù)雜索網(wǎng)的找形提供了一種思路。