邱慧 李瀟 樊俊峰 林秋紅 檀傈錳

(北京空間飛行器總體設(shè)計(jì)部，北京 100094)

?

航天器平面薄膜結(jié)構(gòu)模態(tài)分析和試驗(yàn)

邱慧 李瀟 樊俊峰 林秋紅 檀傈錳

(北京空間飛行器總體設(shè)計(jì)部，北京 100094)

為分析航天器平面薄膜結(jié)構(gòu)(包含薄膜和繩索張拉系統(tǒng))的動(dòng)力學(xué)特性，對(duì)真空和空氣中的多花邊平面薄膜結(jié)構(gòu)的模態(tài)進(jìn)行了有限元仿真分析，且構(gòu)建了一套測(cè)試平面薄膜結(jié)構(gòu)模態(tài)的試驗(yàn)系統(tǒng)，分別利用攝影測(cè)量法和掃描式激光測(cè)振儀對(duì)平面薄膜結(jié)構(gòu)的平面度和模態(tài)進(jìn)行了測(cè)量，試驗(yàn)結(jié)果驗(yàn)證了模態(tài)仿真分析的正確性。進(jìn)一步的分析表明，平面薄膜結(jié)構(gòu)的基頻近似與繩索拉力的平方根呈正比，且基頻隨花邊半徑增大而減小,但減小幅度很小，說明增大繩索拉力對(duì)平面薄膜結(jié)構(gòu)基頻的提高比減小花邊半徑對(duì)基頻的提高要明顯得多。

航天器；平面薄膜結(jié)構(gòu)；模態(tài)分析；空氣附加質(zhì)量

1 引言

平面薄膜結(jié)構(gòu)因具有質(zhì)量輕、收納體積小等優(yōu)點(diǎn)而在大型航天器薄膜天線[1-2](面積約為幾百平方米或者更大)、太陽(yáng)帆[3-4]等方面具有良好的應(yīng)用前景。平面薄膜結(jié)構(gòu)的模態(tài)直接決定或影響著結(jié)構(gòu)型面精度保持、振動(dòng)控制等，而模態(tài)分析的正確性需通過試驗(yàn)來驗(yàn)證，所以開展平面薄膜結(jié)構(gòu)模態(tài)分析和試驗(yàn)研究具有重要意義。

航天器空間環(huán)境為高真空，在地面對(duì)平面薄膜結(jié)構(gòu)進(jìn)行模態(tài)測(cè)試時(shí)需用真空設(shè)備，而真空設(shè)備成本高且易受場(chǎng)地限制，為了降低成本和試驗(yàn)復(fù)雜度，地面試驗(yàn)大多在空氣中進(jìn)行。由于薄膜結(jié)構(gòu)為輕質(zhì)柔性結(jié)構(gòu)，其與周圍介質(zhì)空氣之間的耦合因素不可忽略，所以關(guān)于平面薄膜結(jié)構(gòu)模態(tài)分析和試驗(yàn)研究大多需考慮空氣因素。2002年，劍橋大學(xué)的S. Kukathasan等對(duì)三角平面薄膜結(jié)構(gòu)在空氣中的模態(tài)進(jìn)行了有限元仿真分析，并對(duì)仿真分析的正確性進(jìn)行了試驗(yàn)驗(yàn)證[5]，但分析的薄膜對(duì)象僅有1個(gè)花邊。2006年，美國(guó)噴氣推進(jìn)實(shí)驗(yàn)室(JPL)的Houfei Fang等利用雙參數(shù)模型和分布函數(shù)傳遞法，對(duì)平面薄膜天線結(jié)構(gòu)的模態(tài)進(jìn)行了數(shù)值計(jì)算，并通過試驗(yàn)對(duì)數(shù)值結(jié)果進(jìn)行了驗(yàn)證[6]，但數(shù)值計(jì)算中未涉及空氣對(duì)平面薄膜結(jié)構(gòu)模態(tài)的影響。2010年，上海交通大學(xué)肖薇薇等分析計(jì)算了不同花邊形狀下三角平面薄膜結(jié)構(gòu)的基頻，結(jié)果表明薄膜結(jié)構(gòu)的基頻隨花邊矢跨比(花邊曲線弓形高與弦長(zhǎng)之比)減小而降低[7]，但分析結(jié)果未得到試驗(yàn)驗(yàn)證。2011年，同濟(jì)大學(xué)王磊等利用真空箱測(cè)試了圓形平面薄膜在不同密度空氣中的模態(tài)，且基于試驗(yàn)結(jié)果提出了薄膜振動(dòng)的附加質(zhì)量分布的簡(jiǎn)化模型[8]。2015年，上海交通大學(xué)陳宇峰等基于流固耦合勢(shì)流體原理，對(duì)預(yù)應(yīng)力矩形平面薄膜在空氣中的模態(tài)進(jìn)行了仿真分析及試驗(yàn)驗(yàn)證，結(jié)果表明空氣對(duì)薄膜模態(tài)的影響與薄膜材料和空氣密度有關(guān)[9]。2016年，上海交通大學(xué)邱振宇等基于勢(shì)流理論推導(dǎo)了三角形膜單元的附加質(zhì)量，且在有限元軟件中完成了三角形附加質(zhì)量單元的開發(fā)，并通過圓形平面薄膜的模態(tài)試驗(yàn)驗(yàn)證了該方法[10]。但文獻(xiàn)[8-10]只針對(duì)純薄膜，未涉及繩索張拉系統(tǒng)。

綜上，目前對(duì)于平面薄膜結(jié)構(gòu)模態(tài)分析與試驗(yàn)研究大多局限于單花邊的純薄膜，且未考慮薄膜邊緣的繩索張拉系統(tǒng)。針對(duì)這一問題，本文對(duì)多花邊平面薄膜結(jié)構(gòu)(含薄膜和繩索張拉系統(tǒng))的模態(tài)進(jìn)行仿真分析，且對(duì)仿真分析結(jié)果的正確性進(jìn)行試驗(yàn)驗(yàn)證，在此基礎(chǔ)上獲得平面薄膜結(jié)構(gòu)參數(shù)對(duì)模態(tài)的影響規(guī)律，為大型平面薄膜結(jié)構(gòu)的設(shè)計(jì)提供借鑒。

2 平面薄膜結(jié)構(gòu)模態(tài)分析

平面薄膜結(jié)構(gòu)(含薄膜和張拉系統(tǒng))在自然狀態(tài)下處于無應(yīng)力的狀態(tài)，需通過繩索張拉系統(tǒng)對(duì)薄膜施加預(yù)應(yīng)力。張拉系統(tǒng)是指平面薄膜結(jié)構(gòu)和外框架之間的繩索，通過繩索實(shí)現(xiàn)對(duì)薄膜施加均勻預(yù)應(yīng)力及提供薄膜和外框架間的柔性連接。

在真空狀態(tài)下，薄膜結(jié)構(gòu)的模態(tài)一般稱為干模態(tài)[9]。此時(shí)薄膜結(jié)構(gòu)自由振動(dòng)方程為

(1)

當(dāng)考慮薄膜結(jié)構(gòu)與周圍空氣之間的耦合因素時(shí)，薄膜結(jié)構(gòu)的模態(tài)稱為濕模態(tài)。此時(shí)可把空氣的作用作為附加質(zhì)量來分析，根據(jù)流固耦合基本原理[11]，可以得到薄膜結(jié)構(gòu)的自由振動(dòng)方程為

(2)

式(1)和(2)中：Ms、Cs和Ks分別表示薄膜結(jié)構(gòu)自身質(zhì)量矩陣、阻尼矩陣和剛度矩陣；Ma表示薄膜結(jié)構(gòu)與空氣耦合作用產(chǎn)生的附加質(zhì)量矩陣；Ca為氣動(dòng)阻尼矩陣；r為結(jié)構(gòu)位移矢量。

對(duì)于四邊固定的矩形純薄膜結(jié)構(gòu)，不考慮空氣因素時(shí)，薄膜結(jié)構(gòu)自由振動(dòng)基頻的解析解為

(3)

式中：T0為單位長(zhǎng)度薄膜上張力；m0為單位面積薄膜質(zhì)量；A、B為矩形純薄膜的長(zhǎng)和寬。

根據(jù)文獻(xiàn)[12]，當(dāng)考慮薄膜結(jié)構(gòu)與周圍空氣之間的耦合因素時(shí)，可取薄膜結(jié)構(gòu)的附加質(zhì)量為

ma=0.68ρa(bǔ)L

(4)

式中：ρa(bǔ)為空氣密度；L為薄膜結(jié)構(gòu)的特征尺寸。

假設(shè)附加質(zhì)量均勻分布，根據(jù)式(3)可得到薄膜結(jié)構(gòu)濕模態(tài)的自振頻率為

(5)

本文的多花邊平面薄膜結(jié)構(gòu)(含薄膜和張拉系統(tǒng))構(gòu)型如圖1所示。

圖1 平面薄膜結(jié)構(gòu)構(gòu)型Fig.1 Configuration of flat membrane structure

薄膜邊緣被裁剪制作成曲線索套(花邊)，繩索穿在索套內(nèi)，索套穿出的繩索拉緊固連于外框架(圖中未顯示)上，每個(gè)花邊對(duì)應(yīng)兩個(gè)繩索張拉點(diǎn)，每個(gè)張拉點(diǎn)處有2根繩索。設(shè)薄膜長(zhǎng)為a、寬為b。每個(gè)花邊的跨距為l，花邊半徑為R，花邊角點(diǎn)處2根繩索的間距為d。

當(dāng)薄膜內(nèi)的應(yīng)力處于各向等拉應(yīng)力狀態(tài)時(shí)，薄膜花邊為圓弧形。繩索上拉力和薄膜內(nèi)拉應(yīng)力的關(guān)系[13]為

F=σtR

(6)

式中：F為繩索上的拉力；σ為薄膜內(nèi)應(yīng)力；t為薄膜厚度；R為花邊圓弧半徑。

薄膜材料采用Kapton，繩索材料采用Kelvar，薄膜和繩索的材料參數(shù)如表1所示。

表 1 平面薄膜結(jié)構(gòu)各部分材料參數(shù)

利用ABAQUS有限元軟件分別對(duì)上述平面薄膜結(jié)構(gòu)在真空和空氣中的模態(tài)進(jìn)行仿真分析。平面薄膜結(jié)構(gòu)參數(shù)取值分別為：a=1320 mm，b=790 mm，R=250 mm，l=300 mm，d=30 mm。其中薄膜邊緣一共有14個(gè)繩索張拉點(diǎn)，每個(gè)張拉點(diǎn)處有2根繩索。薄膜采用膜單元(M3D3)進(jìn)行模擬，繩索采用桿單元(T3D2)進(jìn)行模擬，薄膜花邊和繩索之間的連接采用ABAQUS里的綁定(Tie)約束進(jìn)行模擬。

在對(duì)平面薄膜結(jié)構(gòu)干模態(tài)和濕模態(tài)進(jìn)行有限元仿真分析時(shí)，若將濕模態(tài)仿真分析過程中的空氣模型的密度設(shè)為0，即此時(shí)Ma=0,Ca=0，根據(jù)式(1)和(2)，可以得到，此時(shí)濕模態(tài)的仿真分析結(jié)果就為平面薄膜結(jié)構(gòu)干模態(tài)的仿真分析結(jié)果。

首先對(duì)上述平面薄膜結(jié)構(gòu)的濕模態(tài)進(jìn)行仿真分析。利用3D聲單元AC3D20或AC3D15(AC3D20和AC3D15分別代表ABAQUS軟件中的20節(jié)點(diǎn)和15節(jié)點(diǎn)實(shí)體聲音單元)來模擬空氣，此時(shí)空氣被近似為聲音介質(zhì)，它與平面薄膜結(jié)構(gòu)的相互作用近似為作用在薄膜表面的聲壓力。定義一個(gè)長(zhǎng)方體空氣域，平面薄膜結(jié)構(gòu)和空氣域的分析模型如圖2所示，空氣密度[5]為1.23 kg/m3，體積壓縮模量為1.42×105N/m2。薄膜位于空氣域的中央，薄膜單元和聲單元之間的耦合通過面—面接觸實(shí)現(xiàn)，其中一個(gè)面為薄膜表面，另一個(gè)面為空氣域中與薄膜的接觸面，兩者之間采用綁定(Tie)進(jìn)行連接。濕模態(tài)仿真分析分為兩步：①靜力分析，對(duì)各個(gè)繩索的外部張拉點(diǎn)施加相同拉力，使得薄膜和繩索獲得拉應(yīng)力，最終可得到平面薄膜結(jié)構(gòu)的預(yù)應(yīng)力狀態(tài)；②模態(tài)分析，在上步分析得到平面薄膜結(jié)構(gòu)的預(yù)應(yīng)力狀態(tài)后，固定薄膜邊界處各條繩索的外部張拉點(diǎn)，對(duì)薄膜上某一點(diǎn)施加激振力，利用平面薄膜結(jié)構(gòu)對(duì)此激勵(lì)的響應(yīng)得到結(jié)構(gòu)自振頻率。圖3表示了仿真分析得到的空氣中平面薄膜結(jié)構(gòu)的前兩階振型。

將空氣密度設(shè)為0，對(duì)平面薄膜結(jié)構(gòu)的干模態(tài)進(jìn)行仿真分析。分析步驟與濕模態(tài)仿真分析一樣，仿真分析得到的平面薄膜結(jié)構(gòu)的前兩階振型與濕模態(tài)一樣，見圖3。(U并不是代表結(jié)構(gòu)的實(shí)際位移，而是振幅的比值并歸一化的結(jié)果，無量綱。)

圖2 平面薄膜結(jié)構(gòu)與空氣域的分析模型Fig.2 Analytic model of flat membrane structure in air

圖3 平面薄膜結(jié)構(gòu)振型Fig.3 Mode shapes of flat membrane structure

通過仿真分析，得到平面薄膜結(jié)構(gòu)在不同繩索拉力F下的干濕模態(tài)基頻fd和fw，如表2所示。另外，由式(6)，可根據(jù)繩索拉力F得到薄膜內(nèi)拉應(yīng)力σ，進(jìn)而可計(jì)算出薄膜單位長(zhǎng)度上的張力T0。L為薄膜結(jié)構(gòu)的特征尺寸,其值為矩形薄膜的寬度。最后由式(3)、(4)和(5)可計(jì)算出A=1320 mm、B=790 mm、L=790 mm時(shí)，四邊固支的矩形純薄膜的干濕模態(tài)基頻值，見表2。

表2 薄膜結(jié)構(gòu)基頻仿真分析結(jié)果和解析解比較

從表2中可以看出，相對(duì)于矩形純薄膜結(jié)構(gòu)，平面薄膜結(jié)構(gòu)的基頻略有降低，說明加入了繩索后，平面薄膜結(jié)構(gòu)的基頻會(huì)降低。平面薄膜結(jié)構(gòu)的干模態(tài)基頻和濕模態(tài)基頻之間相差很大。這是由于分析的薄膜面積較大，從而附加質(zhì)量較大，導(dǎo)致干濕模態(tài)基頻之間相差較大?？梢娍諝庖蛩貙?duì)平面薄膜結(jié)構(gòu)的模態(tài)影響較大，所以將平面薄膜結(jié)構(gòu)模態(tài)仿真分析結(jié)果和地面室內(nèi)環(huán)境下的模態(tài)試驗(yàn)結(jié)果進(jìn)行對(duì)比時(shí)，仿真分析中不可忽略空氣因素。

3 平面薄膜結(jié)構(gòu)模態(tài)的試驗(yàn)研究

3.1 平面薄膜結(jié)構(gòu)的試驗(yàn)狀態(tài)

為了驗(yàn)證上述仿真結(jié)果的正確性，開展平面薄膜結(jié)構(gòu)模態(tài)試驗(yàn)。試驗(yàn)裝置圖如圖4所示。

圖4 試驗(yàn)裝置圖Fig.4 Schematic of test device

薄膜由繩索懸掛于外框架上，通過繩索調(diào)力裝置，將不同的預(yù)應(yīng)力導(dǎo)入到平面薄膜結(jié)構(gòu)中。外框架尺寸長(zhǎng)2600 mm、高2200 mm，材料選用6060型鋁型材。外框架上分布14個(gè)繩支座，對(duì)應(yīng)薄膜邊緣的14個(gè)繩索張拉點(diǎn)，繩支座可在外框架上任意位置固定。每個(gè)繩支座上安裝2個(gè)繩索調(diào)力裝置，繩索調(diào)力裝置的細(xì)節(jié)圖如圖5所示，電子秤固定在繩支座上，電子秤下面掛著花籃螺栓，花籃螺栓下面掛著彈簧，彈簧與薄膜角點(diǎn)處的繩索相連。通過旋轉(zhuǎn)花籃螺栓調(diào)節(jié)彈簧長(zhǎng)度，從而改變繩索上拉力，繩索上拉力可通過電子秤讀出。

圖5 繩索調(diào)力裝置Fig.5 Device for cable tension adjustment

薄膜試件根據(jù)尺寸裁剪完成后，將薄膜邊緣翻邊制為索套，繩索穿在索套內(nèi)，以此實(shí)現(xiàn)繩索和薄膜的連接。之后將薄膜安裝在外框架的中心位置上。調(diào)節(jié)各個(gè)花籃螺栓的長(zhǎng)度，使各個(gè)繩支座上的電子秤讀數(shù)相同，從而使各根繩索上拉力相同。對(duì)平面薄膜結(jié)構(gòu)進(jìn)行模態(tài)測(cè)試前，需將其調(diào)整到基本無褶皺的狀態(tài)，如圖6所示。

圖6 平面薄膜結(jié)構(gòu)試驗(yàn)前的最終狀態(tài)Fig.6 Final state of flat membrane structure before test

3.2 試驗(yàn)測(cè)量步驟與工況

調(diào)整薄膜各個(gè)角點(diǎn)處繩索的拉力為預(yù)定值后，分別利用攝影測(cè)量法測(cè)量平面薄膜結(jié)構(gòu)的平面度和激光測(cè)振儀測(cè)量平面薄膜結(jié)構(gòu)的模態(tài)。試驗(yàn)所用儀器設(shè)備主要包括：一套攝影測(cè)量設(shè)備，Polytec OFV-505掃描式激光測(cè)振儀，激振器等。試驗(yàn)包括三個(gè)工況：繩索拉力為10 N、20 N和30 N。試驗(yàn)主要步驟如下：

(1)分布測(cè)點(diǎn)。在薄膜上均勻分布35個(gè)測(cè)點(diǎn)，橫向間距為200 mm，縱向間距為150 mm，測(cè)點(diǎn)的分布如圖7所示。

圖7 測(cè)點(diǎn)分布Fig.7 Distribution of measuring points

(2)測(cè)量平面度。不斷調(diào)整薄膜各個(gè)角點(diǎn)處的繩索拉力為10 N后，利用非接觸式的攝影測(cè)量法[14]測(cè)量薄膜的平面度，使薄膜的平面度滿足要求。

(3)測(cè)量模態(tài)。薄膜的平面度滿足要求后，利用Polytec OFV-505掃描式激光測(cè)振儀對(duì)薄膜的模態(tài)進(jìn)行測(cè)試。測(cè)點(diǎn)不變，固定外框架，只對(duì)平面薄膜結(jié)構(gòu)進(jìn)行模態(tài)測(cè)試。給薄膜表面一個(gè)初始激振力(本試驗(yàn)通過激振器給膜面某一點(diǎn)進(jìn)行激振來實(shí)現(xiàn))，利用激光測(cè)振儀，分別對(duì)每個(gè)測(cè)點(diǎn)進(jìn)行掃描，且記錄下各個(gè)測(cè)點(diǎn)的位移、速度、加速度隨時(shí)間的變化規(guī)律。對(duì)各點(diǎn)的掃描值進(jìn)行整合便可得到整個(gè)薄膜的幅值圖像，從而得到平面薄膜結(jié)構(gòu)的基頻和振型。

(4)改變工況。對(duì)薄膜進(jìn)行卸載，調(diào)節(jié)花籃螺栓長(zhǎng)度，使得各個(gè)電子秤讀數(shù)分別為20 N、30 N，在這兩種工況下重復(fù)步驟(2)、(3)的內(nèi)容，可得到這兩種工況下平面薄膜結(jié)構(gòu)的模態(tài)。完成工況為30 N的測(cè)試后，試驗(yàn)結(jié)束。

3.3 試驗(yàn)結(jié)果分析

通過攝影測(cè)量法得到薄膜上各測(cè)點(diǎn)的三維坐標(biāo)，從而得到膜面的平面度。通過不斷調(diào)整，3種工況下測(cè)量的平面度均小于1 mm，膜面基本無褶皺。根據(jù)各測(cè)點(diǎn)的坐標(biāo)在MATLAB軟件中畫出薄膜的三維圖像，如圖8所示。

平面度滿足要求后，便可測(cè)量平面薄膜結(jié)構(gòu)的模態(tài)。通過激光測(cè)振儀得到薄膜上各測(cè)點(diǎn)的振動(dòng)位移隨時(shí)間的變化曲線，可得出平面薄膜結(jié)構(gòu)的基頻，根據(jù)各測(cè)點(diǎn)的位移在MATLAB軟件中可擬合出薄膜的振型圖。圖9表示了當(dāng)繩索拉力為10 N時(shí)平面薄膜結(jié)構(gòu)中心矩形區(qū)域內(nèi)的前兩階振型的仿真分析結(jié)果圖和試驗(yàn)測(cè)量圖。

從圖9中看出，仿真分析和試驗(yàn)測(cè)量得到的振型圖基本一致。試驗(yàn)所測(cè)的振型并不完全對(duì)稱，其原因是試驗(yàn)過程中各個(gè)角點(diǎn)處繩索上的拉力與預(yù)期值稍有偏差，造成薄膜內(nèi)應(yīng)力分布不均勻，從而導(dǎo)致平面薄膜結(jié)構(gòu)的振型并非完全對(duì)稱。

圖 8 平面薄膜結(jié)構(gòu)平面度的測(cè)試結(jié)果Fig.8 Testing result of membrane structure flatness

圖9 平面薄膜結(jié)構(gòu)振型仿真與試驗(yàn)的對(duì)比Fig.9 Comparison of simulating and testing results for modal shapes

表3給出了繩索拉力分別為10 N、20 N、30 N時(shí)平面薄膜結(jié)構(gòu)基頻的仿真和試驗(yàn)結(jié)果。

表 3 平面薄膜結(jié)構(gòu)基頻的仿真與試驗(yàn)結(jié)果比較

從表3中可以看出，試驗(yàn)結(jié)果ft與仿真分析中的干模態(tài)基頻fd相差較大(70%～74%)，而與濕模態(tài)基頻fw相近,ft與fw的誤差均在10%以內(nèi)。因試驗(yàn)是在室內(nèi)環(huán)境下進(jìn)行，而濕模態(tài)仿真分析過程中考慮了空氣因素，所以試驗(yàn)結(jié)果與濕模態(tài)仿真分析結(jié)果相近，從而說明了平面薄膜結(jié)構(gòu)濕模態(tài)仿真分析的正確性。由式(1)和(2)可以得到，將濕模態(tài)仿真分析中空氣模型的密度設(shè)為0，就可以得到平面薄膜結(jié)構(gòu)的干模態(tài)。所以，濕模態(tài)仿真分析的正確性也進(jìn)一步說明了干模態(tài)仿真分析的正確性。

4 平面薄膜結(jié)構(gòu)參數(shù)對(duì)模態(tài)的影響分析

在上文仿真分析的基礎(chǔ)上，分析繩索上拉力和花邊半徑對(duì)平面薄膜結(jié)構(gòu)干模態(tài)基頻的影響規(guī)律。保持其他參數(shù)不變，改變繩索上拉力，平面薄膜結(jié)構(gòu)的基頻隨繩索拉力的變化如圖10(a)所示，從圖中可以看出，平面薄膜結(jié)構(gòu)的基頻近似與繩索拉力的平方根成正比。

保持繩索上拉力不變，改變花邊半徑大小，平面薄膜結(jié)構(gòu)基頻隨花邊半徑的變化如圖10(b)所示，從圖中可以看出，基頻隨花邊半徑的增大而減小，但減小的幅度很小。

圖10 平面薄膜結(jié)構(gòu)基頻隨結(jié)構(gòu)參數(shù)的變化規(guī)律Fig.10 Fundamental frequency of flat membrane structure versus structural parameters

通過對(duì)比圖10(a)和10(b)可以看出，當(dāng)繩索拉力從10 N增大至20～60 N時(shí)，平面薄膜結(jié)構(gòu)的基頻呈41.4%～144.7%幅度不等的提高；而當(dāng)花邊半徑從200 mm增大至250～450 mm時(shí)，平面薄膜結(jié)構(gòu)的基頻只呈1.2%～1.8%幅度不等的降低。說明通過增大繩索拉力來提高平面薄膜結(jié)構(gòu)基頻的效果比減小花邊半徑來提高基頻的效果要明顯得多。

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

(1)本文對(duì)多花邊平面薄膜結(jié)構(gòu)的干模態(tài)和濕模態(tài)進(jìn)行了仿真分析，且通過試驗(yàn)對(duì)仿真結(jié)果進(jìn)行了驗(yàn)證，結(jié)果對(duì)比見表3，試驗(yàn)與仿真結(jié)果吻合(相對(duì)誤差小于10%)，說明本文仿真分析的正確性。

(2)空氣因素對(duì)平面薄膜結(jié)構(gòu)的模態(tài)有較大影響，所以將模態(tài)仿真結(jié)果與地面室內(nèi)環(huán)境下的模態(tài)測(cè)試結(jié)果進(jìn)行比對(duì)時(shí)，仿真分析中不可忽略空氣因素。

(3)平面薄膜結(jié)構(gòu)(包括薄膜和繩索張拉系統(tǒng))的基頻近似與繩索上拉力的平方根呈正比；且結(jié)構(gòu)的基頻隨花邊半徑增大而減小，但減小幅度很小。兩者對(duì)比說明通過增大繩索拉力來提高平面薄膜結(jié)構(gòu)基頻的效果比減小花邊半徑來提高基頻的效果要明顯得多。

References)

[1]Delapierre M,Pellegrino S. Spin-stabilized membrane antenna structures[C]//2nd AIAA Spacecraft Structures Conference. Washington D.C.:AIAA,2015:1-15

[2]Potvin M J,Montminy S,Brunel S,et al.Testing of a deployable SAR membrane antenna mechanical prototype[C]//49th AIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics,and Materials Conference. Washington D.C.:AIAA,2008:1-12

[3]左華平，馮煜東，王虎，等．太陽(yáng)帆航天器研究進(jìn)展及其關(guān)鍵技術(shù)分析[J]. 真空科學(xué)與技術(shù)學(xué)報(bào)，2016，36(1):117-124

Zuo Huaping,Feng Yudong,Wang Hu,et al. Latest development of solar sail research and analysis of its key technology [J]. Chinese Journal of Vacuum Science and Technology,2016,36(1):117-124 (in Chinese)

[4]Okuizumi N,Satou Y,Mori O,et al.Investigation of the deformation and stiffness of spinning solar sail membrane of IKAROS[C]//4th AIAA Spacecraft Structures Conference. Washington D.C.:AIAA,2017:1-10

[5]Kukathasan S,Pellegrino S. Vibration of prestressed membrane structures in air [C]//43rd AIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics,and Materials Conference. Washington D.C.:AIAA,2002:1-11

[6]Fang H，Yang B，Ding H，et al．Dynamic analysis of large in-space deployable membrane antennas [C]//The 13th International Congress on Sound and Vibration. Vienna:ICSV,2006:1-6

[7]肖薇薇，陳務(wù)軍，付功義．空間薄膜陣面預(yù)應(yīng)力導(dǎo)入效應(yīng)及影響因素[J]．宇航學(xué)報(bào)，2010，31(3)：845-849

Xiao Weiwei,Chen Wujun,Fu Gongyi. Pre-stress introduction effects and influence factors investigation for the space planar film reflect-array [J]. Journal of Astronautics,2010,31(3):845-849 (in Chinese)

[8]王磊，李元齊，沈祖炎．薄膜振動(dòng)附加質(zhì)量試驗(yàn)研究[J]．振動(dòng)工程學(xué)報(bào)，2011，24(2)：125-132

Wang Lei,Li Yuanqi,Shen Zuyan. Experimental investigation on the added mass of membranes vibrating in air[J]. Journal of Vibration Engineering,2011,24(2):125-132 (in Chinese)

[9]陳宇峰，陳務(wù)軍，邱振宇，等．空氣對(duì)預(yù)應(yīng)力薄膜結(jié)構(gòu)模態(tài)的影響[J]．浙江大學(xué)學(xué)報(bào)(工學(xué)版)，2015，49(6)：1123-1127

Chen Yufeng,Chen Wujun,Qiu Zhenyu,et al. Effects of air on modal behavior of pre-stressed membrane structure [J]. Journal of Zhejiang University (Engineering Science),2015,49(6):1123-1127 (in Chinese)

[10]邱振宇，陳務(wù)軍，趙兵，等．充氣尾翼濕模態(tài)分析與試驗(yàn)研究[J]．振動(dòng)與沖擊，2016，35(15)：140-148

Qiu Zhenyu,Chen Wujun,Zhao Bing,et al. Wet modal analysis and tests for pneumatic empennages[J]. Journal of Vibration and Shock,2016,35(15):140-148 (in Chinese)

[11]李慶祥，孫炳楠，顧正維．考慮空氣耦合的薄膜結(jié)構(gòu)動(dòng)力分析[J]．計(jì)算力學(xué)學(xué)報(bào)，2006，23(4)：476-482

Li Qingxiang,Sun Bingnan,Gu Zhengwei.Dynamic analysis of membrane structures interacting with air [J]. Chinese Journal of Computational Mechanics,2006,23(4):476-482 (in Chinese)

[12]Minami H. Added mass of a membrane vibrating at finite amplitude[J]. Journal of Fluids and Structures,1998,12:919-932

[13]Greschik G,White C V,Salama M A. On the precisely uniform and uniaxial tensioning of a film sheet via integrated catenary[C]//44th AIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics,and Materials Conference. Washington D.C.: AIAA,2003: 1-11

[14]Shimoda Y,Watanabe K,Sakamoto N,et al. Development of stereo camera system for accurate observation of deployable membranes onboard cubesat[C]//55th AIAA Aerospace Sciences Meeting. Washington D.C.: AIAA,2017: 1-7

(編輯：張小琳)

Modal Analysis and Experiment of Spacecraft Flat Membrane Structures

QIU Hui LI Xiao FAN Junfeng LIN Qiuhong TAN Limeng

(Beijing Institute of Spacecraft System Engineering,Beijing 100094,China)

In order to analyze dynamic characteristics of flat membrane structures (including membrane and tensioning system),modal analysis for a flat membrane structure with multiple pockets in vacuum and air is carried out by finite element simulation and a modal testing system is built. Flatness of the flat membrane structure is measured by photogrammetric testing technique and the modal is tested by a scanning laser vibrometer. Validity of the modal simulation analysis is demonstrated by experimental results. Based on the verified modeling technique,further research results show that the fundamental frequency of flat membrane structures is approximately proportional to the square root of cable tension force and the fundamental frequency decreases with the in creasing of radius of pocket,but the decrease magnitude is small. Therefore,in order to increase the fundamental frequency of flat membrane structures,increasing cable tension force is much more efficient than decreasing radius of pocket.

spacecraft; flat membrane structure; modal analysis; air added mass

2017-03-30；

2017-05-16

邱慧，女，碩士研究生，研究方向?yàn)楹教炱鹘Y(jié)構(gòu)與機(jī)構(gòu)設(shè)計(jì)。Email:qh_sofa@126.com。

V414

A

10.3969/j.issn.1673-8748.2017.03.007