Ω形可展開彎桿屈曲載荷的經(jīng)驗公式與試驗驗證

韓大為，朱曉海，靳宗向，楊 青，曹燕燕

Ω形可展開彎桿屈曲載荷的經(jīng)驗公式與試驗驗證

韓大為1，2，朱曉海2，靳宗向1，2，楊 青1，2，曹燕燕1，2

（1．上海市空間飛行器機構(gòu)重點實驗室，上海201109；2．上海宇航系統(tǒng)工程研究所，上海201109）

針對Ω形可展開彎桿多參數(shù)耦合、設計困難的問題，對Ω形可展開彎桿的屈曲載荷、兩側(cè)剛度差異特點等進行了仿真分析，擬合了計算公式；并研制了樣機，進行了試驗測試，仿真結(jié)果和試驗結(jié)果匹配性良好。研究內(nèi)容為多參數(shù)Ω形可展開彎桿的快速設計和推廣使用積累了技術(shù)。

Ω形可展開彎桿；可展開結(jié)構(gòu)；屈曲載荷；數(shù)據(jù)試驗

1 引言

近些年，隨著空間利用的快速發(fā)展，空間碎片數(shù)量急劇增加［1］。針對空間碎片清理需求，相關新型飛行器設計正引起國內(nèi)外專家的重視［2?3］。本文提出一種非破壞、可重復、具有大展收比、不需要目標飛行器提供特定捕獲接口的、半剛性捕獲機構(gòu)，如圖1所示，該捕獲機構(gòu)的核心部件是一種橫截面為Ω形的空間彎曲半剛性桿。這種Ω形桿的兩側(cè)剛度具有明顯差異（剛度差異特點），工作前可以卷曲收攏，工作時釋放展開，具有大展收比的特征。剛度大的內(nèi)側(cè)可以用作對太空碎片的柔性捕獲工作面，剛度小的外側(cè)可以用作卷曲收攏和展開釋放的接觸面。彎曲Ω形桿除了可以用作半剛性捕獲爪，高精度的多根彎曲Ω形桿還可用作新型拋物面天線傘狀支撐結(jié)構(gòu)，可以克服傳統(tǒng)拼接式拋物面天線整體弧度較差的問題。

傳統(tǒng)的直Ω形桿多拼接成豆莢桿進行應用，是一種一維直線伸展的結(jié)構(gòu)，國內(nèi)外已經(jīng)進行了較多的研究［4?6］；但彎曲Ω形桿國內(nèi)外還沒有相關的研究，是一種平面內(nèi)二維伸展結(jié)構(gòu)，如圖2所示。它具有內(nèi)外側(cè)剛度不同的特點，其內(nèi)側(cè)剛度更大，外側(cè)剛度更小，更適合于單側(cè)剛度要求高、收攏包絡要求嚴苛的應用場合。

Ω形可展開彎桿具有展收比大、內(nèi)側(cè)剛度大等優(yōu)點，是一種有廣泛應用前景的新型結(jié)構(gòu)，但是其橫截面曲率和縱向曲率具有耦合關系，如何計算其屈曲載荷和兩側(cè)剛度差異特點，對于Ω形彎桿的承載能力計算和參數(shù)優(yōu)化等具有重要影響，本文針對Ω形彎桿的屈曲載荷和剛度差異特點進行仿真分析和試驗驗證。

2 仿真分析

Ω形可展開彎桿的整體半徑R、橫截面半徑r、整體弦長L、橫截面弦長l、材料厚度t等指標參數(shù)示意如圖3、圖4所示，對Ω形彎桿的力學特性具有影響，但是何種參數(shù)組合情況下，Ω形彎桿的力學性能能夠滿足需求，通過試驗進行逐一驗證是不現(xiàn)實的，所以通過仿真掌握參數(shù)變化對于彎桿性能的影響具有重要意義。

2? 1 仿真數(shù)模幾何參數(shù)

Ω形可展開彎桿的內(nèi)側(cè)剛度較大，作為捕獲工作面；外側(cè)剛度較小，作為收攏工作面，用于彎桿發(fā)射入軌前的收攏，具體如圖5所示。在工作中，被捕獲的物體會對半剛性桿產(chǎn)生接觸力F，如果該力過大，桿件就會發(fā)生屈曲。因此需要對不同幾何參數(shù)下的桿件進行屈曲分析，尋找桿幾何參數(shù)與臨界屈曲載荷的關系。

針對多尺寸參數(shù)（橫截面曲率、縱向曲率、厚度、寬度、弦長等參數(shù)）的Ω形可展開彎桿進行分析。Ω形可展開彎桿的物理特性參數(shù)統(tǒng)計如表1所示，幾何參數(shù)組合統(tǒng)計如表2所示。

表1 半剛性桿物性參數(shù)統(tǒng)計表Table 1 Physical parameters of semi?rigid bars

通過對多參數(shù)組合的Ω形可展開彎桿進行仿真歸納，總結(jié)出關于Ω形彎桿屈曲載荷和剛度差異特點的經(jīng)驗公式或選型表格。其中，屈曲載荷定義為Ω形可展開彎桿內(nèi)側(cè)承載出現(xiàn)屈曲時候的最大承載；剛度差異特點定義為Ω形可展開彎桿內(nèi)外側(cè)同樣承載條件下，內(nèi)側(cè)抗彎剛度／外側(cè)抗彎剛度的比值。

表2 半剛性桿幾何參數(shù)組合統(tǒng)計表Table 2 Geometry parameters of semi?rigid bars／mm

2? 2 仿真前處理和參數(shù)設置

2? 2? 1 網(wǎng)格劃分

根據(jù)以上所述的試驗形式，在仿真分析中采用MSC．Patran 2012進行幾何建模。按照表2的某一組參數(shù)繪制三維曲面，使其橫截面及整體都呈圓弧形，與實際桿件完全一致［7］。在曲面的端部取4 mm分割出來，用來施加載荷，如圖6所示。將曲面劃分成殼單元，彎桿橫截面上定義6或12個單元，彎桿整體圓弧方向上，除用來施加載荷的4 mm段，定義100～200個單元，各組幾何參數(shù)下殼單元的長寬比都保持在1? 3～1? 7之間。

2? 2? 2 邊界條件及載荷

在屈曲試驗中，半剛性桿的固支段整體被固支，因此有限元分析中可以僅僅取桿的工作段進行建模，并約束其根部的所有自由度，對應實際桿件的橫截面。在分割出來的4 mm段曲面上，對所有節(jié)點施加外載荷，力矢量方向試驗工裝加載方向保持一致，對應實際試驗時砝碼重力方向不變的特點，仿真中不做跟隨力處理。

2? 2? 3 求解器參數(shù)設置

根據(jù)該分析任務大變形、大應變，材料非線性的特點，采用MSC．Nastran 2012的sol600求解器進行隱式非線性準靜態(tài)分析。半剛性桿幾何參數(shù)改變后，臨界屈曲力會發(fā)生較大改變，因此不同工況下所施加的外載荷大小也不盡相同，需要相應地略加調(diào)整。采用固定步長的方式施加載荷增量，但如果一個增量步內(nèi)求解不收斂，允許將載荷步自動二分。采用位移收斂判據(jù)，其余所有參數(shù)采用默認設置［8］。各組幾何參數(shù)情況所用的外載荷及增量步的數(shù)目也一并在后續(xù)表格中列出。

2? 3 仿真分析結(jié)果

2? 3? 1 Ω形可展開彎桿屈曲臨界載荷分析結(jié)果

Ω形可展開彎桿（R＝500 mm，r＝20 mm，L＝700 mm，l＝32 mm，t＝0? 4 mm，參數(shù)彎桿示例）受載荷作用后發(fā)生屈曲的變形情況及應力云圖如圖7所示。彎桿在載荷逐漸增大過程中，會不斷變形，彎桿各點位移也隨著載荷的增大逐漸變化。將該關系繪制成位移?載荷曲線，曲線坡度陡變處對應的即是臨界屈曲載荷點。

Ω形可展開彎桿多參數(shù)組合下屈曲載荷統(tǒng)計如表3所示。

如果能在設計Ω形彎桿時，可以通過公式估算出臨界屈曲載荷，將對Ω形彎桿的參數(shù)設計帶來極大的便利。接下來本文嘗試將典型工況下的仿真數(shù)據(jù)進行擬合，得到這樣的經(jīng)驗公式。如果采用非線性回歸需要事先明確模型的解析形式，再用實驗數(shù)據(jù)來擬合未知參數(shù)，但本課題還沒有對Ω形彎桿臨界屈曲載荷的解析形式進行研究，文獻調(diào)查的結(jié)果也沒有相關記錄，因此目前適合采用響應面方法。作為線性回歸的一種，該方法用多項式模型來描述響應，多項式模型可以有任意多個變量xi（i＝1，…，N）［9］。下面以完全二次函數(shù)為例（見式（1））、結(jié)合本次實例對響應面模型進行簡要介紹。

式（1）的模型包含橫截項、線性項、二次交互項、平方項。線性項產(chǎn)生超平面，交互項則使超平面發(fā)生翹曲，平方項生成使響應面取極值的最簡模型。半剛性桿結(jié)構(gòu)的弦長與半徑呈比例關系L＝1? 4R，l＝1? 6r，因此模型的變量僅取三個：整體半徑R，橫截面半徑r，材料厚度t。將其中兩個變量固定，臨界屈曲力隨另外一個變量變化的響應曲線如圖8所示，紅色虛線代表90%置信區(qū)間?？梢钥闯?，在徑向R＝700 mm，橫截面r＝27 mm，材料厚度t＝0? 4附近，臨界屈曲力的誤差最小，為±7 N左右。原因是上述尺寸接近各個工況的統(tǒng)計中心。當徑向R小于400 mm或大于800 mm時，響應的置信度降低；橫截面r小于20 mm時，響應的置信度降低；材料厚度t小于0? 3 mm時，響應的置信度降低。

表3 Ω形可展開彎桿多參數(shù)組合下屈曲載荷統(tǒng)計Table 3 Critical buckling loads of deployable Ω?shape bars with different combination of parameters

公式擬合后需要檢驗其誤差，可以采用調(diào)整多項式的項和階次的方法優(yōu)化結(jié)果。經(jīng)過試驗后發(fā)現(xiàn)，完全三次模型相比完全二次模型誤差更小，但三次以上的模型則改進甚微，有些檢驗點還會增大。完全三次模型下得到的臨界屈曲載荷經(jīng)驗公式如式（2）所示：

用各個試驗點來檢驗臨界屈曲載荷經(jīng)驗公式，以仿真值與公式計算值的偏差占仿真值的百分比為縱坐標，來繪制點圖，如圖9所示，可知在表3所用的幾何參數(shù)范圍內(nèi)，相對誤差小于15%。

采用完全三次多項式擬合，固定R變量時，響應面云圖如圖10所示。

彎桿整體半徑為200 mm時，橫截面半徑小，則材料越厚、臨界屈曲載荷越大；橫截面半徑大于50 mm時，改變材料厚度幾乎沒有影響。

彎桿整體半徑為400 mm時，橫截面半徑小，則材料越厚、臨界屈曲載荷越大；橫截面半徑大于40 mm時，改變材料厚度幾乎沒有影響，但橫截面半徑達到50 mm以上時，材料越厚，臨界屈曲載荷反而會下降。

2? 3? 2 Ω形可展開彎桿兩側(cè)抗彎剛度差異特點分析結(jié)果

結(jié)構(gòu)的剛度是指彈性體抵抗變形的能力。這里沿用屈曲分析中的約束及加載策略測試外側(cè)抗彎剛度，剛度定義為k＝F／δ，F(xiàn)為屈曲前的載荷大小，δ為對應的載荷方向位移。再保持端部載荷的大小不變，將其反向加載，用同樣方法計算內(nèi)側(cè)抗彎剛度，最后計算得出兩側(cè)抗彎剛度差異特點（抗彎剛度差異特點＝內(nèi)側(cè)抗彎剛度／外側(cè)抗彎剛度）。仍然采用MSC．Nastran 2012的sol600求解器進行隱式非線性準靜態(tài)分析。載荷以固定步長增量的方式施加，所用工況均統(tǒng)一為10步，在一個增量步內(nèi)求解不收斂的情況下允許將載荷步自動二分，采用位移收斂判據(jù)，其余所有參數(shù)采用默認設置。

在同一工況條件下（選取R＝250 mm，r＝10 mm，t＝0? 2 mm典型工況為例），在Ω形可展開彎桿內(nèi)、外側(cè)施加同樣大小的端部載荷的變形情況分別如圖11、圖12所示。

各組工況所用的外載荷及對應的抗彎剛度統(tǒng)計如表4所示。

表4 剛度差異特點仿真結(jié)果Table 4 Simulation results of the rigidity difference

由Ω形可展開彎桿剛度差異特點的仿真結(jié)果可知，Ω形彎桿的兩側(cè)抗彎剛度差異特點≥1? 5，內(nèi)外側(cè)剛度差異明顯。

3 試驗驗證

選取R250彎桿（R＝250 mm，r＝10 mm，t＝0? 2 mm）和R500彎桿（R＝500 mm，r＝20 mm，t＝0? 5 mm）兩種典型參數(shù)的Ω形可展開彎桿進行加工制備，并設計試驗加工試驗工裝，對Ω形彎桿的屈曲力和抗彎剛度進行試驗驗證。試驗工裝中把兩個Ω形彎桿中間通過剛性桿連接，在剛性桿中間進行加載，通過高度尺測量變形數(shù)據(jù)。剛性桿的彈性模量足夠大，忽略它自身的變形，則試驗中兩個Ω形彎桿并聯(lián)實測值變形除以2可以認為是單根Ω形彎桿的變形值，可用于與仿真結(jié)果對照。

3? 1 彎桿屈曲臨界載荷

分別對R500彎桿和R250彎桿進行正面懸吊砝碼逐步加載測試，直至彎桿出現(xiàn)屈曲，記錄其屈曲承載載荷，加載過程如圖13所示。

R500和R250彎桿的內(nèi)側(cè)屈曲承載能力測試結(jié)果如表5所示。

表5 R500和R250屈曲載荷實測結(jié)果Table 5 Measured buckling load of R500 and R250

3? 2 兩側(cè)抗彎剛度差異特點

分別對R250彎桿和R500彎桿進行正反向加載，測試二者在表4對應載荷下的剛度差異特點，正反向試驗加載和測試如圖14、圖15所示。

R250彎桿和R500彎桿正反向變形試驗結(jié)果如表6所示。

表6 剛度差異特點實測結(jié)果Table 6 Measured rigidity difference

3? 3 試驗結(jié)果分析

通過試驗結(jié)果和仿真結(jié)果的對比，R500彎桿和R250彎桿實測屈曲載荷，與仿真經(jīng)驗公式計算值之間的偏差不大于9? 2%；R500彎桿和R250彎桿實測剛度差異特點，與仿真結(jié)果之間的偏差不大于9? 5%，可以用于指導Ω形可展開彎桿的快速參數(shù)設計。

4 結(jié)論

本文提出一種可用于空間碎片清理捕獲機構(gòu)和拋物面天線機構(gòu)的Ω形可展開彎桿，并對其內(nèi)側(cè)承載屈曲載荷進行了仿真分析，歸納了經(jīng)驗公式，還仿真了Ω形彎桿的兩側(cè)剛度差異特點。針對具有代表性的R500彎桿和R250彎桿設計生產(chǎn)了測試工裝，并進行了試驗驗證。屈曲載荷和剛度差異特點實測結(jié)果與仿真結(jié)果偏差分別不大于9? 2%和9? 5%，仿真結(jié)果可行度較高。本文的研究內(nèi)容為Ω形可展開彎桿的多參數(shù)快速設計，尤其是為分析不同組合參數(shù)Ω形可展開彎桿的極限承載能力（屈曲載荷）、兩側(cè)剛度差異特點計算積累了經(jīng)驗公式和仿真計算方法，可作為Ω形彎桿推廣使用的基礎。

（References）

［1］ 林來興．空間碎片現(xiàn)狀與清理［J］．航天器工程，2012，21（3）：1?10．Lin Lai?xing．Status and removal of space debris［J］．Space?craft Engineering，2012，21（3）：1?10．（in Chinese）

［2］ 陳小前，袁建平，姚雯，等．航天器在軌服務技術(shù)［M］．北京：中國宇航出版社，2009：89?101．Chen Xiao?qian，Yuan Jian?ping，Yao Wen，et al．The Tech?nology of Spacecraft in Orbit Service［M］．Beijing：China As?tronautic Publishing House，2009：89?101．（in Chinese）

［3］ 蔡洪亮，高永明，邴啟軍，等．國外空間非合作目標抓捕系統(tǒng)研究現(xiàn)狀與關鍵技術(shù)分析［J］．裝備指揮技術(shù)學院學報，2010，21（6）：71?77．Cai Hong?liang，Gao Yong?ming，Bing Qi?jun，et al．The re?search status and key technology analysis of foreign non?coop?erative target in space capture system［J］．Journal of the A?cademy of Equipment Command＆Technology，2010，21（6）：71?77．（in Chinese）

［4］ 鞠博文，齊俊偉，李濤，等．針對半豆莢桿的先進拉擠成型壓力研究［J］．航空材料學報，2015，35（6）：20?26．Ju Bo?wen，Qi Jun?wei，Li Tao，et al．Research on molding pressure of advanced pultrusion for half lenticular tube［J］．Journal of Aeronautical Materials，2015，35（6）：20?26．（in Chinese）

［5］ 丁峻宏，咸奎成，韓軒，等．空間豆莢桿機構(gòu)收展過程的并行仿真計算［J］．宇航學報，2011，32（3）：676?682．Ding Jun?hong，Xian Kui?cheng，Han Xuan，et al．Parallel computing for simulation of stowing and deployment process of space lenticular boom mechanism［J］．Journal of Astronau?tics，2011，32（3）：676?682．（in Chinese）

［6］ 白江波，熊峻江，高軍鵬，等．可折疊復合材料豆莢桿的制備與驗證［J］．航空學報，2011，32（7）：1217?1223．Bai Jiang?bo，Xiong Jun?jiang，Gao Jun?peng，et al．Fabrica?tion and validation of collapsible composite lenticular tubes［J］．Acta Aeronautica ET Astronautica Sinica，2011，32（7）：1217?1223．（in Chinese）

［7］ 牛慶芳，尹永青．薄壁桿件的有限元分析法［J］．科技情報開發(fā)與經(jīng)濟，2006，16（13）：140?143．Niu Qing?fang，Yin Yong?qing．The finite element method for the thin?walled bar［J］．Sci?Tech Information Development＆Economy，2006，16（13）：140?143．（in Chinese）

［8］ 蔚春燕，陳軍．薄壁桿件結(jié)構(gòu)力學分析的方法探討［J］．黑龍江科技信息，2010（21）：251?251．Wei Chunyan，Chen Jun．Discussion on the method of me?chanical analysis of thin?walled bar［J］．Heilongjiang Science and Technology Information，2010（21）：251?251．（in Chi?nese）

［9］ 陳光南，黃文彬．開口薄壁桿件的彈塑性有限元分析［J］．北京農(nóng)業(yè)工程大學學報，1987（3）：104?112．Chen Guang?nan，Huang Wen?bing．Elastic?plastic finite ele?ment analysis of open thin?walled members［J］．Journal of Be?ing Agricultuer Engineering University，1987（3）：104?112．（in Chinese）

（責任編輯：康金蘭）

Empirical Formula and Verification Test of Buckling Load in
Deployable Ω?shape Curved Bars

HAN Dawei1，2，ZHU Xiaohai2，JIN Zongxiang1，2，YANG Qing1，2，CAO Yanyan1，2
（1．Shanghai Key Laboratory of Spacecraft and Machine，Shanghai 201109，China；2．Aerospace System Engineering Shanghai，Shanghai 201109，China）

In this paper，the problem of multi?parameter coupling and design difficulty of deployable Ω?shape curved bar was studied．The buckling loads and rigidity difference of the deployable Ω?shape curved bar were simulated and the fitting formulas were obtained．The prototype testing showed that the experimental results and the simulation results matched well．The research lays a good foundation for the rapid design and promotion of deployable Ω?shape curved bars．

deployable Ω?shape curved bar；deployable structure；buckling load；data test

TH123? 1

A

1674?5825（2017）04?0462?07

2017?01?10；

2017?03?02

國家高技術(shù)發(fā)展計劃（863計劃）（2013AA7043021）

韓大為，男，碩士，工程師，研究方向為空間載荷機構(gòu)設計。E?mail：handawei021＠163．com