可展開雙Ω薄殼結(jié)構(gòu)動(dòng)態(tài)特性測(cè)試與分析

衛(wèi)劍征，丁海鑫，侯 雪，譚惠豐

式中，ω0為結(jié)構(gòu)的固有頻率，m為結(jié)構(gòu)軸向線密度，y為結(jié)構(gòu)自由端的振幅。彈性勢(shì)能形式儲(chǔ)存在結(jié)構(gòu)中，因此勢(shì)能的變化值可以通過彎矩做功的大小來表示。而彎矩做功可以表示為式（7）：

由于結(jié)構(gòu)自由端彎曲撓度數(shù)值很小，所以有式（8）：

可展開雙Ω薄殼結(jié)構(gòu)動(dòng)態(tài)特性測(cè)試與分析

衛(wèi)劍征，丁海鑫，侯 雪，譚惠豐

（哈爾濱工業(yè)大學(xué)復(fù)合材料與結(jié)構(gòu)研究所，哈爾濱150080）

針對(duì)一種碳纖維復(fù)合材料雙Ω薄殼結(jié)構(gòu)中折疊展開變形以及局部損傷對(duì)動(dòng)態(tài)特性的影響問題，基于Rayleigh方程，引入了纖維復(fù)合材料結(jié)構(gòu)中鋪層相關(guān)因子，討論卷曲展開程度及薄壁結(jié)構(gòu)局部損傷對(duì)結(jié)構(gòu)動(dòng)態(tài)特性的影響程度；建立了雙Ω薄殼結(jié)構(gòu)的動(dòng)態(tài)測(cè)試系統(tǒng)，采用了單點(diǎn)激勵(lì)多點(diǎn)響應(yīng)的動(dòng)態(tài)測(cè)試方法，得到了雙Ω薄殼結(jié)構(gòu)分別在展開、卷曲收攏和局部屈曲破壞后的一階、二階固有頻率。測(cè)試結(jié)果驗(yàn)證了數(shù)值仿真結(jié)果的正確性，并表明卷曲后雙Ω薄殼結(jié)構(gòu)基頻增大，局部屈曲破壞后的雙Ω薄殼結(jié)構(gòu)基頻減小。

雙Ω截面；薄殼結(jié)構(gòu)；動(dòng)態(tài)；卷曲；碳纖維復(fù)合材料

1 引言

本文引入纖維復(fù)合材料結(jié)構(gòu)中鋪層相關(guān)因子，討論卷曲展開程度及薄壁結(jié)構(gòu)局部損傷對(duì)動(dòng)態(tài)特性的影響程度，并建立了雙Ω薄殼管的動(dòng)態(tài)測(cè)試系統(tǒng)，采用了單點(diǎn)激勵(lì)多點(diǎn)響應(yīng)的動(dòng)態(tài)測(cè)試方法，研究了雙Ω薄殼結(jié)構(gòu)分別在展開、卷曲收攏和局部屈曲損傷后的動(dòng)態(tài)特性。

2 結(jié)構(gòu)動(dòng)態(tài)測(cè)試方法

為了測(cè)試雙Ω薄殼管不同展開狀態(tài)下的動(dòng)態(tài)特性，本文設(shè)計(jì)了結(jié)構(gòu)動(dòng)態(tài)特性測(cè)試方案：把雙Ω薄殼結(jié)構(gòu)水平放置，一端與支架固支安裝，使薄殼結(jié)構(gòu)為懸臂梁狀態(tài)；沿長度方向每間隔200 mm布置一個(gè)三軸加速度傳感器（PCB公司），共為5個(gè)；激勵(lì)方式為定點(diǎn)激勵(lì)，激勵(lì)位置為自由端，激勵(lì)方向?yàn)閥和z方向，測(cè)試方案如圖1所示。其中測(cè)試對(duì)象為外徑60 mm、長度100 cm的雙Ω薄殼結(jié)構(gòu)，由四層（45°／－45°／－45°／45°）碳纖維復(fù)合材料制備而成，其彈性模量為22? 4 GPa，剪切模量為8? 62 GPa，泊松比為0? 33，管壁線密度為62 g／m。振動(dòng)測(cè)試數(shù)據(jù)采集系統(tǒng)為DEWE?801。在管壁上均勻粘貼三軸傳感器，在自由端處利用端頭蓋保持截面的固定，防止開口出現(xiàn)局部振動(dòng)，利用帶有軟式錘頭的力錘對(duì)結(jié)構(gòu)自由端進(jìn)行激勵(lì)，以減小結(jié)構(gòu)出現(xiàn)局部變形，從而得到雙Ω薄殼結(jié)構(gòu)的整體動(dòng)態(tài)特性。

依據(jù)動(dòng)態(tài)特性測(cè)試方案，建立了動(dòng)態(tài)特性測(cè)試實(shí)物系統(tǒng)（圖2），包括DEWE?801振動(dòng)測(cè)試儀、薄殼結(jié)構(gòu)支架、傳感器等。測(cè)試過程中需要先設(shè)定傳感器靈敏系數(shù)，選擇測(cè)量的通道，再用軟頭力錘對(duì)自由端進(jìn)行激勵(lì)。采用模態(tài)測(cè)試方法選擇各階固有頻率和其相應(yīng)的阻尼，最后依據(jù)有限元計(jì)算得到振型，并選擇確定試驗(yàn)中對(duì)應(yīng)的各階基頻。

本文動(dòng)態(tài)試驗(yàn)中采用了單點(diǎn)激勵(lì)多點(diǎn)響應(yīng)的模態(tài)測(cè)試方法，由于試驗(yàn)測(cè)量的雙Ω薄殼結(jié)構(gòu)為一種低頻結(jié)構(gòu)，選擇采樣頻率為1 kHz／s，分辨率設(shè)定為1024，半帶寬度為500 Hz，得到試驗(yàn)頻率精度為0? 2918 Hz。利用力錘對(duì)雙Ω薄殼結(jié)構(gòu)進(jìn)行激勵(lì)，為降低人為激勵(lì)誤差導(dǎo)致出現(xiàn)非固有頻率峰值，采用三次激勵(lì)取平均值方法，對(duì)結(jié)構(gòu)端部進(jìn)行激勵(lì)，根據(jù)每個(gè)響應(yīng)點(diǎn)的振動(dòng)曲線，通過傅立葉變換擬合出雙Ω薄殼結(jié)構(gòu)的振型，并計(jì)算得到結(jié)構(gòu)固有頻率。本文采用三軸傳感器設(shè)置z方向沿著管長度方向，xy方向?yàn)殡pΩ薄殼結(jié)構(gòu)橫截面內(nèi)方向；固支端由四枚Φ4 mm螺絲加墊片緊密固定在金屬支架帶孔平板上；三軸傳感器分別對(duì)應(yīng)DEWE?801測(cè)試系統(tǒng)1、2、3、4和5測(cè)試通道，力錘與0通道連接，力錘錘頭選擇橡膠材料。

模態(tài)參數(shù)辨識(shí)方法采用頻域法，動(dòng)態(tài)信號(hào)測(cè)試儀得到頻響函數(shù)后，也得到導(dǎo)納圓的測(cè)試數(shù)據(jù)。其中導(dǎo)納圓法是利用動(dòng)態(tài)特性測(cè)試中得到的頻響函數(shù)，擬合成不同的固有頻率的模態(tài)圓。導(dǎo)納圓辨識(shí)方法是利用頻響函數(shù)的數(shù)據(jù)進(jìn)行辨識(shí)［13］。對(duì)于簡化所得單自由度結(jié)構(gòu)阻尼系統(tǒng)，其頻響函數(shù)如式（1）：

式中，（1＋jg）k為復(fù)剛度；g為結(jié)構(gòu)阻尼因子。將上式按實(shí)部與虛部分開寫得式（2）～（4）：

(4)從典型斷面流速變化來看，方案三(閘孔凈寬120 m)和方案四(閘孔凈寬240 m)實(shí)施后，相對(duì)于建閘前，流速增大的區(qū)域包括進(jìn)口單一段、新河段和匯流點(diǎn)以下單一段，流速減小的區(qū)域?yàn)槔虾拥缽澢?。且方案四?duì)水流的影響小于方案三。

由上兩式消去變量（1－w－2）后，得到圓方程如式（5）：

其圓心坐標(biāo)為（0，－1／2kg），半徑為1／2kg，此圓為導(dǎo)納圓或Nyquist圓。該圓是以HR（ω）為實(shí)部、HI（ω）為虛部在復(fù)平面上形成的圖形。

3 測(cè)試結(jié)果與分析

3? 1 展開狀態(tài)測(cè)試分析

用導(dǎo)納圓辨識(shí)法確定結(jié)構(gòu)的固有頻率，理論上試驗(yàn)結(jié)果應(yīng)精確地存在于導(dǎo)納圓上，但由于測(cè)試中存在誤差，數(shù)值結(jié)果不能精確地存在于導(dǎo)納圓上。本文通過曲線擬合的方法將測(cè)試結(jié)果數(shù)值擬合成一個(gè)圓，并得到其各階固有頻率。在測(cè)試中通過模態(tài)測(cè)試方法識(shí)別雙Ω薄殼結(jié)構(gòu)的各階頻率，主要包括頻率特性曲線、頻率隨振幅的響應(yīng)曲線和頻率隨相位角的響應(yīng)曲線。首先用雙重離散傅里葉變換（FFT）曲線（圖3）上相應(yīng)波峰作為固有頻率可能取值的點(diǎn)；同時(shí)頻率隨響應(yīng)變化的幅頻圖中（圖4），曲線波峰處作為薄殼結(jié)構(gòu)固有頻率可能取到的點(diǎn)，并結(jié)合導(dǎo)納圓辨識(shí)法（圖5），判斷出雙Ω薄殼結(jié)構(gòu)的固有頻率值。

測(cè)試得到懸臂狀的雙Ω薄殼結(jié)構(gòu)的yz面一階基頻為8? 5 Hz，xz面一階基頻為10? 2 Hz。由于結(jié)構(gòu)在x和y方向不對(duì)稱，其慣性矩Ixx為1? 6939×104mm4，Iyy為5? 5205×104mm4，所以得到的前兩階固有頻率不相等。根據(jù)Rayleigh方程，基于機(jī)械能守恒定律，結(jié)構(gòu)的固有頻率可由結(jié)構(gòu)振動(dòng)的動(dòng)能最大值和彈性勢(shì)能最大值相等得出。其動(dòng)能表示為式（6）：

式中，ω0為結(jié)構(gòu)的固有頻率，m為結(jié)構(gòu)軸向線密度，y為結(jié)構(gòu)自由端的振幅。彈性勢(shì)能形式儲(chǔ)存在結(jié)構(gòu)中，因此勢(shì)能的變化值可以通過彎矩做功的大小來表示。而彎矩做功可以表示為式（7）：

由于結(jié)構(gòu)自由端彎曲撓度數(shù)值很小，所以有式（8）：

其中，θ為結(jié)構(gòu)偏離初始位置角度，R為結(jié)構(gòu)的曲率半徑。又由材料力學(xué)中對(duì)梁彎曲曲率的描述有式（9）：

將式（8）、（9）代入式（7）中可以得到式（10）：

由能量守恒定律可知，當(dāng)形變等于零時(shí)能夠得到動(dòng)能最大值；當(dāng)形變達(dá)到最大值時(shí)，得到勢(shì)能最大值，即動(dòng)能最大值和彈性勢(shì)能最大值相等，所以有式（11）：

懸臂結(jié)構(gòu)上與固定端距離等于x處的撓度可以表示成式（12）：

根據(jù)Rayleigh方程有Dmax＝Umax，則得到結(jié)構(gòu)的固有頻率ω0如式（15）：

由于雙Ω薄殼結(jié)構(gòu)采用了復(fù)合材料層合結(jié)構(gòu)，與簡化的懸臂結(jié)構(gòu)計(jì)算得到的公式存在一定偏差，因此引入了與復(fù)合材料結(jié)構(gòu)鋪層相關(guān)因子如式（16）：

當(dāng)ζ＝5時(shí)，計(jì)算得到雙Ω薄殼結(jié)構(gòu)在yz面和xz面的基頻為9? 36 Hz、18? 74 Hz。

3? 2 收攏狀態(tài)測(cè)試分析

進(jìn)一步測(cè)試了雙Ω薄殼結(jié)構(gòu)在卷曲下的固有頻率，設(shè)計(jì)了一個(gè)直徑為150 mm的卷軸，其質(zhì)量與端蓋近似相同，質(zhì)量為58 g。自由端繞卷軸將雙Ω薄殼結(jié)構(gòu)卷起，并將其固定，如圖6所示。采用同樣的測(cè)試方法，通過導(dǎo)納圓、幅頻圖和FFT曲線得到一階基頻（圖7）。測(cè)得yz面一階基頻為9? 3 Hz，xy面一階基頻為17? 8 Hz（表1）。

表1 薄殼結(jié)構(gòu)卷曲前后測(cè)試結(jié)果對(duì)比Table 1 Comparison of test results of thin?shell struc?ture before and after curving

根據(jù)基頻計(jì)算公式（15），基頻與試件的長度和質(zhì)量有關(guān)，當(dāng)質(zhì)量或長度增大時(shí)，基頻減小。當(dāng)雙Ω薄殼結(jié)構(gòu)卷曲后，其長度減小，增加的卷軸質(zhì)量與卸掉的端頭蓋相當(dāng)，但卷曲部分的雙Ω薄殼結(jié)構(gòu)質(zhì)量也施加在了試件的自由端，限制了卷曲后的基頻變大。因此結(jié)構(gòu)的基頻在卷曲后仍然具有增大的趨勢(shì)，且二階基頻具有更明顯的變化。

3? 3 局部損傷測(cè)試分析

針對(duì)雙Ω薄殼結(jié)構(gòu)軸壓后發(fā)生的局部損傷程度進(jìn)行了動(dòng)態(tài)測(cè)試分析，局部破壞的雙Ω薄殼結(jié)構(gòu)如圖8所示。裂紋位置在距固支端80 mm的上半圓弧處，裂紋長度約為30 mm，四層鋪層均斷裂。

采用同樣無損傷的100 cm薄殼結(jié)構(gòu)測(cè)試方法和條件，得到相頻圖（圖9）和結(jié)構(gòu)前六階頻率（表2）。從表2中看出經(jīng)過軸壓試驗(yàn)后產(chǎn)生局部破壞的雙Ω薄殼結(jié)構(gòu)固有頻率有所降低，這是由于端部結(jié)構(gòu)開裂導(dǎo)致結(jié)構(gòu)的剛度退化，根據(jù)式（11）對(duì)基頻的計(jì)算，剛度的退化會(huì)導(dǎo)致結(jié)構(gòu)基頻下降。

表2 局部裂紋的雙Ω薄殼結(jié)構(gòu)固有頻率Table 2 The inherent frequency of locally cracking double Ω thin?shell boom

4 結(jié)構(gòu)動(dòng)態(tài)特性有限元分析

4? 1 有限元計(jì)算

本文建立了雙Ω薄殼結(jié)構(gòu)模型［14］并進(jìn)行了動(dòng)態(tài)特性有限元計(jì)算，在有限元計(jì)算中采用四節(jié)點(diǎn)減縮積分的薄殼單元，建立相同的一端固支、一端自由邊界，并在自由端部建立端頭蓋剛體，與雙Ω薄殼結(jié)構(gòu)通過Tie方式連接，兩端截面設(shè)置耦合約束以防止開口局部振動(dòng)，得到結(jié)構(gòu)的振動(dòng)頻率，計(jì)算得到y(tǒng)z面一階頻率為11? 42 Hz，xz面一階頻率為20? 85 Hz（圖10）。表3為雙Ω結(jié)構(gòu)試驗(yàn)得到的前兩階固有頻率與有限元計(jì)算基頻結(jié)果對(duì)比，結(jié)果表明：與展開狀態(tài)基頻對(duì)比，卷曲基頻比初始基頻增大，軸壓后局部裂紋破壞的雙Ω薄殼結(jié)構(gòu)基頻降低44? 8%。

4? 2 誤差分析

盡管雙Ω薄殼結(jié)構(gòu)的Iyy約為Ixx的4倍，由公式（11）中ω0和慣性矩的關(guān)系，xz面內(nèi)基頻應(yīng)是yz面基頻的2倍關(guān)系，證明結(jié)果的正確性，但是理論與測(cè)試和有限元仿真之間還是存在誤差。引起誤差的原因主要有三方面：首先，薄殼結(jié)構(gòu)截面慣性矩的理論計(jì)算，引入了結(jié)構(gòu)輔層相關(guān)因子，通過經(jīng)驗(yàn)確定的ζ，該參數(shù)需要多次、多種細(xì)長比結(jié)構(gòu)的測(cè)試進(jìn)行修正；其次，測(cè)試過程中存在試件尺寸、材料阻尼、激勵(lì)方向等耦合，以及系統(tǒng)誤差等因素；最后有限元模態(tài)計(jì)算為線性計(jì)算方法，沒有考慮實(shí)際測(cè)試對(duì)象碳纖維復(fù)合材料中樹脂的阻尼。

表3 有限元計(jì)算固有頻率結(jié)果與試驗(yàn)及理論值對(duì)比Table 3 Comparison of inherent frequencies calculated by finite element analysis，by experiment and by theory

5 結(jié)論

1）通過對(duì)懸臂結(jié)構(gòu)的基頻進(jìn)行推導(dǎo)分析，得到了基頻與慣性矩，長度和質(zhì)量之間的關(guān)系，驗(yàn)證了不對(duì)稱截面雙Ω薄殼結(jié)構(gòu)在展開和卷曲狀態(tài)下基頻的準(zhǔn)確性。

2）雙Ω薄殼結(jié)構(gòu)軸壓局部裂紋破壞后基頻降低，比無損傷的基頻降低了44? 8%；且結(jié)構(gòu)在卷曲狀態(tài)下的基頻與初始態(tài)基頻進(jìn)行了對(duì)比，卷曲收攏狀態(tài)的基頻比展開狀態(tài)的基頻大。

（References）

［1］ 衛(wèi)劍征，毛麗娜，杜星文．空間卷曲折疊管充氣控制展開動(dòng)力學(xué)研究［J］．工程力學(xué)，2009，26（1）：227?232．Wei J Z，Mao L N，Du X W．Study for inflatable control de?ployment dynamics of rolled booms［J］．Engineering Mechan?ics，2009，26（1）：227?232．（in Chinese）

［2］ Bartonb A，Puiga L，Randoa N．A review on large deployable structures for astrophysics missions［J］．Acta Astronautica，2010，67：12?26．

［3］ Blockn J，Straubel M，Wiedemann M．Ultralight deployable booms for solar sails and other large gossamer structures in space［J］．ACTA Astronautica，2011，68：984?992．

［4］ Blandino J，Duncan R，Nuckels M，et al．Three?dimensional shape sensing for inflatable booms［C］／／46th AIAA／ASME／ASCE／AHS／ASC Structures，Structural Dynamics and Mate?rials Conference，AIAA，2005：1807．

［5］ Roh J H，Kim H I，Lee S Y．Viscoelastic effect on unfolding behaviors of shape memory composite booms［J］．Composite Structures，2015，133：235?245．

［6］ Aguirre M M，Bowen D H，Davidson R，et al．Development of a continuous manufacturing method for a deployable satellite mast in CFRP［C］／／15th Reinforced Plastics Congress，Not?tingham，Engl，1986：107?110．

［7］ Sickinger C，Herbeck L．Deployment strategies，analyses and tests for the CFRP booms of a solar sail［J］．DLR?Mitteilung，2004，1：215?225．

［8］ Khoddam S，Hakkak F．On calculation of preliminary design parameters for lenticular booms［J］．Journal of Aerospace Engineering，2007，211：377?384．

［9］ 房光強(qiáng)，彭福軍．航天器可展開支撐桿的研制及其收攏展開特性研究［J］．材料工程，2009，（S2）：157?160．Fang G Q，Peng F J．Fabrication and retraction／deployment testings of space deployable booms［J］．Materials Engineer?ing，2009，（S2）：157?160．（in Chinese）

［10］ 白江波，熊峻江，高軍鵬，等．可折疊復(fù)合材料豆莢桿的制備與驗(yàn)證［J］．航空學(xué)報(bào)，2011，32（7）：1217?1223．Bai J B，Xiong J J，Gao J P，et al．Fabrication and validation of collapsible composite lenticular tubes［J］．Astronautica Sinica，2011，32（7）：1217?1223．（in Chinese）

［11］ 劉靖，周偉．CFRP薄壁管空間伸展臂模態(tài)與影響參數(shù)分析及試驗(yàn)［J］．東華大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，39（5）：577?582．Liu J，Zhou W．Experiments and parametric influence analy?sis on modal for CFRP thin?walled tube space boom［J］．Journal of Donghua University（Natural Science），2013，39（5）：577?582．（in Chinese）

［12］ 鄒濤．薄壁管狀空間伸展臂收展數(shù)值模擬與試驗(yàn)［D］．上海：上海交通大學(xué)，2013：35?60．Zou T．Numerical Simulation and Tests of Stowing and Deploying Process of Thin?walled Tubular Space Booms［D］．Shanghai：Shanghai Jiao Tong University，2013：35?60．（in Chinese）

［13］ 鄒經(jīng)湘，于開平．結(jié)構(gòu)動(dòng)力學(xué)（第二版）［M］．哈爾濱：哈爾濱工業(yè)大學(xué)出版社，2009：7?32．Zhou J X，Yu K P．Dynamics of Structures（Second Edition）［M］．Harbin：Haerbin Institute of Technology Press，2009：7?32．（in Chinese）

［14］ 侯雪．可展開雙Ω薄殼結(jié)構(gòu)力學(xué)性能研究［D］．哈爾濱：哈爾濱工業(yè)大學(xué)，2016：58?63．Hou X．Research of Mechanical Properties of Deployable double?ω Thin?shell Structure［D］．Harbin：Harbin Institute of Technology，2016：58?63．（in Chinese）

（責(zé)任編輯：龐迎春）

Dynamic Testing and Analysis of Deployable Double?Ω Thin?Shell Structures

WEI Jianzheng，DING Haixin，HOU Xue，TAN Huifeng
（Center for Composite Materials and Structures，Harbin Institute of Technology，Harbin 150080，China）

As far as the influences of the local damage and folded deformation on the dynamic char?acteristics are concerned，a double?Ω thin?shell structure was studied．A related layer factor was in?troduced in the fiber composite material structure based on Rayleigh equation．Its dynamic system was also built，and the method of excitating at a single?point and responding at multi points was used to test the first and the second order fundamental frequency of the double?Ω boom under folding，un?folding or local cracking conditions respectively．The test results are consistent with the numerical simulation．The results show that the double?Ω thin?shell boom had a higher frequency under folding and a lower frequency under cracking condition．

double?Ω boom；thin shell structure；dynamics；folded；carbon fibrous composite ma?terial

V414

A

1674?5825（2017）04?0551?05

2017?03?09；

2017?05?31

中央高校基本科研專項(xiàng)資金資助（HIT．MKSTISP．201609）

衛(wèi)劍征，男，博士，副教授，研究方向?yàn)榭臻g柔性可展開結(jié)構(gòu)動(dòng)力學(xué)。E?mail：weijz＠hit．edu．cn