楊 曉，王雪明，何 巍，顧文華

(中國特種飛行器研究所，湖北 荊門 448035)

1 引言

系留氣球利用囊體內(nèi)輕于空氣的氣體提供的靜升力長時間駐空，主要用于低空偵查、通信中繼以及環(huán)境監(jiān)測等領(lǐng)域。近年來，隨著薄膜結(jié)構(gòu)和新型復合材料等技術(shù)的快速發(fā)展，針對系留氣球的研究也受到越來越多人的關(guān)注。蝴蝶結(jié)也被稱為抓手或載荷簾，通常熱合在系留氣球的囊體表面，起著連接外部索具和傳遞索具荷載的作用，主要由襯布、承力條和拉環(huán)組成。襯布采用可以與囊體材料相融合的新型復合材料，承力條具有多種形式(織物類、織帶類)，熱合或者膠結(jié)在襯布表面，拉環(huán)用于連結(jié)承力條和外部索具。實際工程中根據(jù)承載力的不同，可對承力條的數(shù)量、材質(zhì)及扇形面積進行相應調(diào)整。

由于蝴蝶結(jié)具有良好的荷載擴展性，且制造工藝簡便，所以廣泛應用于系留氣球，飛艇等大型柔性結(jié)構(gòu)中，其優(yōu)良的力學性能對整體結(jié)構(gòu)的安全性具有非常重要的影響。蝴蝶結(jié)通常采用編織的新型復合材料薄膜，而纖維織物具有復雜的微觀結(jié)構(gòu)且編織結(jié)構(gòu)具有方向性，故在材料層次和結(jié)構(gòu)層次表現(xiàn)出明顯的非線性和各向異性。在實際中，蝴蝶結(jié)基本處于雙向張拉狀態(tài)，所以研究在相同荷載條件下的膜材力學響應特性非常關(guān)鍵。有關(guān)雙軸拉伸力學特性的研究表明，膜材的雙向受力性能具有明顯的非線性，且彈性參數(shù)受到不同應力比和應力組合的影響。魯國富等對蒙皮材料進行7種應力比下的雙軸拉伸試驗，建立應力彈性常數(shù)響應面，并通過有限元仿真驗證其非線性模型的正確性。陳建穩(wěn)等針對某芳綸纖維層合膜材進行了一系列雙軸拉伸試驗，得到了應力空間中材料彈性參數(shù)的響應曲面，并總結(jié)應力比和應力組合對膜材雙向拉伸力學性能的作用規(guī)律。目前國內(nèi)外對蒙皮材料的力學性能研究較多，但對蝴蝶結(jié)的研究相對較少，所以有必要深入研究它的力學響應特征。

首先設(shè)計了蝴蝶結(jié)試件的力學響應試驗。研究試件在雙向1.6kN張拉荷載下，在蝴蝶結(jié)拉環(huán)處分級施加面內(nèi)拉力2kN、4kN、6kN、8kN、10kN時膜面產(chǎn)生的變形、應力與應變響應；然后利用有限元軟件對試驗過程進行數(shù)值仿真分析；最后提出一種將試驗與仿真結(jié)果定量對比的方法，對試驗與仿真結(jié)果進行定量對比和誤差分析。驗證了仿真方案的可行性以及實驗數(shù)據(jù)的準確性。

2 蝴蝶結(jié)試驗

2.1 試件材料及結(jié)構(gòu)

蝴蝶結(jié)試件的蒙皮材料與承力條分別采用具有高強度，高比強，高模量的芳綸纖維層合膜材。

將蝴蝶結(jié)布置在蒙皮試件中部位置，如圖1所示。蝴蝶結(jié)承力條的彎折后成一定夾角，如圖2所示，每根承力條由多層膜材焊接而成。

圖1 蒙皮試件示意 圖2 蝴蝶結(jié)示意

2.2 試件預處理

首先將預制好的散斑貼在如圖3所示的矩形紅框內(nèi)，作為試驗測試初始區(qū)域(最終區(qū)域取決于試驗設(shè)備)。由于柔性膜材具有很強的非線性且加載過程容易產(chǎn)生褶皺，為了方便與數(shù)值結(jié)果對比且提高對比結(jié)果的準確性，提出一種定量對比的方法。在試件表面的重點關(guān)注區(qū)域，以50mm為間距，水平做編號為～的矩形區(qū)域，測量區(qū)域內(nèi)所有結(jié)點的水平應變(垂直加載方向)、剪切應變。豎直做編號為～的矩形區(qū)域，測量區(qū)域內(nèi)所有結(jié)點的豎直應變。沿承力帶方向做編號為～的矩形區(qū)域，測量區(qū)域內(nèi)所有結(jié)點的最大主應變。在每個工況下，取每個矩形區(qū)域中所有測試點的平均值作為最終結(jié)果。對比傳統(tǒng)的單節(jié)點數(shù)值對比方法，該方法顯著的提高了結(jié)果的準確性和可靠性。其中，矩形區(qū)域的邊長均為15mm，且應保證各區(qū)域測試數(shù)據(jù)點不少于10個。散斑區(qū)域分布位置如圖3中小方塊所示。

圖3 矩形區(qū)域分布位置

2.3 試驗設(shè)備

試驗環(huán)境溫度應在-5℃-40℃之間，相對濕度應小于85%RH，采用典型連接結(jié)構(gòu)強度試驗平臺，平臺采用自平衡框架的設(shè)計方式，無需額外約束，應用PLC程序控制，能夠?qū)崿F(xiàn)試驗件的自動張緊、加載、采集等功能。采用三維數(shù)字散斑應變測量及分析系統(tǒng)(DIC-3D)獲取試件加載后的位移、應力和應變數(shù)據(jù)。

2.4 試驗方法

首先將預處理后的試驗件安裝在平臺上，完全約束試件的四條邊，且施加雙向1.6kN的張緊力，作為試件的初始狀態(tài)；然后架設(shè)照明系統(tǒng)，使光圈照射在散斑區(qū)域，以提高圖像的分辨率；接著架設(shè)圖像采集系統(tǒng)，對其進行標定并鎖定焦距，鏈接遙控拍照控制器，設(shè)定相機在每循環(huán)中以固定間隔時間拍攝相同數(shù)量的照片，調(diào)試照相機，使其能清晰拍攝到測試區(qū)域的散斑情況；然后啟動試驗機，以60mm/min速度對蝴蝶結(jié)拉環(huán)施加面內(nèi)載荷至2kN，待穩(wěn)定后觀察散斑變形情況并拍照記錄后繼續(xù)分級加載至4kN，6kN，8kN，10kN，重復上一步。最后利用DIC-3D對每級的照片進行計算處理，得到散斑區(qū)域內(nèi)各工況下的應力、應變和位移場結(jié)果。采集標定區(qū)域內(nèi)所有測量點的數(shù)據(jù)后進行記錄。

2.5 試驗結(jié)果

由于試驗設(shè)備的局限性，只能獲取部分散斑測試區(qū)域的結(jié)果，最終選取10kN的工況結(jié)果進行分析。

2.5.1 h區(qū)域結(jié)果分析

圖4為試件的部分散斑區(qū)域相對于初始狀態(tài)的水平應變云圖。可以看出最大水平應變處于承力帶兩側(cè)。由于承力條由5層膜材粘合而成，剛度明顯增加，牽引兩側(cè)膜材發(fā)生變形。最小水平應變出現(xiàn)在中間位置的承力條區(qū)域，隨著荷載的增加，應變值由2kN荷載下的-0.4%逐漸增大10kN荷載下的-7.6%，與此同時，膜面局部發(fā)生受壓變形與褶皺，但褶皺占主要成分。

圖4 水平方向應變云圖

由于試驗設(shè)備原因，只能獲取～區(qū)域結(jié)果，各個工況下的數(shù)值結(jié)果如圖5所示。對比水平應變圖，隨著荷載的增加，各數(shù)據(jù)點的應變值呈現(xiàn)增加的趨勢。并且中間點的水平應變絕對值大于其它數(shù)據(jù)點，可見蝴蝶結(jié)荷載變形過程中，因承力條的剛度大，區(qū)域因局部剛度小，呈現(xiàn)水平負應變，表現(xiàn)出水平收縮。

圖5 h數(shù)據(jù)點水平應變結(jié)果

圖6為試件的剪切云圖，可見區(qū)域的剪切應變絕對值最小，向兩側(cè)擴展過程逐漸增大，在拉伸過程中，兩側(cè)表現(xiàn)出較大的剪切變形，區(qū)域以拉應變?yōu)橹?。圖7為各個工況下的數(shù)值結(jié)果，隨著荷載的增加，各區(qū)域沿兩側(cè)擴展時的應變絕對值逐漸增大，與云圖結(jié)果一致。

圖6 剪切方向應變云圖

圖7 h數(shù)據(jù)點剪切應變圖

252區(qū)域結(jié)果分析

圖8為試件的部分散斑區(qū)域相對于初始狀態(tài)的豎直應變云圖，可以看出承力帶區(qū)域的豎直應變相對較小。最大豎直應變出現(xiàn)在承力帶下側(cè)均勻區(qū)域，隨著荷載的增加，應變值由2荷載下的092增大至10荷載下的265。

圖8 豎直方向應變云圖

試驗獲取的～區(qū)域各工況下的數(shù)值結(jié)果如圖9所示。對比豎直應變云圖，隨著荷載的增加，各數(shù)據(jù)點的應變值呈現(xiàn)增大的趨勢。并且中間點、的豎直應變出現(xiàn)較大躍升，主要是承力條剛度較大，點以內(nèi)屬于承力條區(qū)域，因剛度大導致應變較小。點在承力條區(qū)域之外，因單層膜材剛度較小，導致應變值出現(xiàn)較大提高。

圖9 v區(qū)域豎直應變圖

253區(qū)域結(jié)果分析

圖10為試件的部分散斑區(qū)域相對于初始狀態(tài)的最大主應變云圖，試驗獲取的～區(qū)域的數(shù)值結(jié)果如圖11所示。隨著荷載的增加，各數(shù)據(jù)點的應變值呈現(xiàn)增大的趨勢。其中，、區(qū)域的豎直應變出現(xiàn)較大躍升，是因為承力條剛度較大，點以內(nèi)屬于承力條區(qū)域，因剛度大而應變較小。點在承力條區(qū)域之外，因單層膜材剛度較小，其應變值出現(xiàn)較大提高。

圖10 最大主應力應變云圖

圖11 f數(shù)據(jù)點最大主應變圖

3 數(shù)值模擬

3.1 有限元建模

利用三維建模軟件Catia對蝴蝶結(jié)試件建模后導入中預處理，根據(jù)圖3試驗件中標定的矩形區(qū)域位置和散斑位置，在模型中設(shè)置相應的節(jié)點集合，以方便在后處理中提取對應的數(shù)值結(jié)果。由于膜材為柔性結(jié)構(gòu)，實際中不可抗壓及抗彎，所以為模型劃分精細化網(wǎng)格時，采用能考慮幾何大變形且不抗彎的膜單元M3D3和M3D4。為模擬試件不同區(qū)域的鋪層效應，分區(qū)域設(shè)置材料厚度如圖12所示。

在膜材試驗過程中出現(xiàn)了較大的幾何變形，且為緩慢擬靜力加載，為準確模擬實驗過程，數(shù)值仿真模型采用考慮幾何非線性的通用靜力分析。為了提高模型的收斂性，分為兩個分析步，第一分析步通過設(shè)置位移邊界條件來施加與試驗相同的膜面預張力，以此作為模型的初始狀態(tài)。第二分析步在蝴蝶結(jié)拉環(huán)處逐漸施加拉力，均打開幾何非線性開關(guān)。

圖12 10kN工況下的應變云圖

依據(jù)試驗對模型施加邊界條件，左側(cè)邊界為固定約束，其余三邊均施加適當?shù)姆ㄏ蛭灰埔阅M試驗中1.6kN的預張力。實驗平臺每邊對膜材約束較弱，模擬自由雙向應力狀態(tài)，計算模型同時釋放橫向約束，如圖13所示。

圖13 數(shù)值模型邊界約束

為模擬集中加載條件，將5層厚度區(qū)域的邊結(jié)點與加載參考點(蝴蝶結(jié)拉環(huán))耦合，采用分布耦合方式連結(jié)，如圖14所示。分布耦合可將荷載傳遞到膜材的耦合點上，而不會約束膜材耦合點間的相對自由度，膜材耦合點依然保留獨立變形的能量，以準確實現(xiàn)加載過程。

圖14 分布耦合方式加載

3.2 數(shù)值模擬結(jié)果

選取10kN工況下的結(jié)果進行分析。圖15-圖22分別為模型相對于初始狀態(tài)下的水平，豎直，剪切方向的應變，應力與位移云圖。

圖15 水平方向應變云圖

圖16 水平方向應力云圖

圖15-17為水平方向的數(shù)值結(jié)果?？煽闯鰬υ茍D與應變云圖表現(xiàn)一致，應變、應力集中出現(xiàn)在承力條的頂端部位，在承力條兩側(cè)出現(xiàn)較大的應力分布，并持續(xù)延續(xù)到夾持邊界。可見在加載過程中，試件上下兩側(cè)邊界處對變形提供一定的約束。在承力條受拉荷載的相反方向，出現(xiàn)水平受壓的應變、應力區(qū)間。蝴蝶結(jié)試件在水平方向表現(xiàn)出對稱的位移云圖。

圖17 水平方向位移云圖

圖18-20為豎直方向的數(shù)值結(jié)果?？煽闯鰬υ茍D與應變云圖表現(xiàn)一致。承力條受荷載方向的另一側(cè)出現(xiàn)較大的應力分布，并持續(xù)延續(xù)到夾持邊界。在承力條外側(cè)，單層膜材區(qū)域出現(xiàn)較大應力，該區(qū)域是結(jié)點受力較重要的位置。蝴蝶結(jié)試件在豎直方向表現(xiàn)出對稱的位移云圖。在試件加載位置處具有最大的變形，變形自承力條區(qū)域向邊界逐漸遞減。

圖18 豎直方向應變云圖

圖19 豎直方向應力云圖

圖20 豎直方向位移云圖

圖21-22為試件相對于初始狀態(tài)的剪切應變云圖?？梢钥闯鰬兂式茖ΨQ分布，下端的應變最小，上端最大。受力過程中變形主要發(fā)生在承力帶兩側(cè)的單層材料部分，隨著荷載的增加，承力帶兩側(cè)的褶皺區(qū)域越來越明顯。

圖21 剪切方向應變云圖

圖22 剪切方向位移云圖

4 試驗與仿真結(jié)果對比分析

根據(jù)圖3標定的矩形區(qū)域位置，分別提取10kN工況下試驗與仿真的數(shù)值結(jié)果，做定量對比和誤差分析。

4.1 h區(qū)域結(jié)果對比

因水平應變絕對值接近零，對比計算值與試驗值的意義不大，僅對比了剪切應變的誤差情況，如表1所示?？梢钥闯鲈囼灁?shù)據(jù)與仿真結(jié)果對比良好。雖然膜單元依然具有抗壓的能力，但在集中荷載拉伸下，中間受壓單元的剪切應變接近于零。微小的測量、計算偏差就能引入巨大的相對計算誤差，因此沒有計算相對誤差。最終平均誤差為5954。

表1 h區(qū)域剪切應變數(shù)據(jù)對比表(%)

4.2 v區(qū)域結(jié)果對比

表2中列出了區(qū)域的豎直應變誤差情況，因點在測量區(qū)域邊緣，由于測試誤差較大將其除去。組區(qū)域與試驗結(jié)果吻合較好，在承力條區(qū)域，因剛度較大縱向應變較小。自點開始進入單層膜材區(qū)域，因單層膜材剛度較小，其應變值出現(xiàn)較大提高。最終平均誤差為8019。

表2 V區(qū)域豎直應變數(shù)據(jù)對比表(%)

4.3 f區(qū)域結(jié)果對比

表3中列出了區(qū)域的最大主應變誤差情況，因f點在測量區(qū)域邊緣，測試誤差較大，所以除去。f組數(shù)據(jù)與試驗結(jié)果趨勢基本吻合，在承力條區(qū)域，因剛度較大，應變值較小。f是承力區(qū)域到單層膜面區(qū)域的過渡位置，單層膜材剛度較小，其應變值出現(xiàn)較大提高，存在一定的褶皺變形。計算誤差與試驗誤差較大，最終平均誤差計算值6.757%。

表3 f區(qū)域最大主應變數(shù)據(jù)誤差計算表(%)

5 結(jié)論

通過數(shù)值仿真與試驗對蝴蝶結(jié)處于不同荷載下的力學響應進行分析。研究了蝴蝶結(jié)膜面在試驗工況下的變形、應力與應變響應及分布規(guī)律。提出一種試驗與數(shù)值模擬結(jié)果定量對比的方法，發(fā)現(xiàn)對比結(jié)果的誤差在10以內(nèi)，驗證了試驗數(shù)據(jù)的有效性與數(shù)值模擬方法的正確性。針對其它的同類試驗可以利用該仿真方法進行分析，對蝴蝶結(jié)的力學設(shè)計與試驗研究有著指導借鑒的作用。