折疊翼緩沖裝置參數(shù)確定分析

李玉亮，付 凡

（1．中國航空工業(yè)洪都，江西 南昌330024；2．空軍駐江西地區(qū)軍事代表室，江西 南昌330024）

折疊翼緩沖裝置參數(shù)確定分析

李玉亮1，付 凡2

（1．中國航空工業(yè)洪都，江西 南昌330024；2．空軍駐江西地區(qū)軍事代表室，江西 南昌330024）

對某折疊翼所采用的金屬橡膠緩沖器的設(shè)計參數(shù)進行了研究。本文采用準(zhǔn)靜態(tài)接觸有限元分析與ADAMS動力學(xué)仿真結(jié)合的方法，將剛體動力學(xué)接觸碰撞剛度與有限元接觸剛度相融合，確定了等效的緩沖裝置線彈性材料參數(shù)，最后得到了不同驅(qū)動力與緩沖裝置組合下折疊機構(gòu)展開到位的沖擊載荷，可為展開機構(gòu)設(shè)計及緩沖裝置工程研制提供參考。

折疊翼；金屬橡膠；緩沖；接觸分析；動力學(xué)仿真

0 引言

在某縱向折疊翼展開到位過程中，翼面沖擊面碰撞壓縮金屬橡膠緩沖片，通過緩沖片變形將折疊翼動能轉(zhuǎn)化為緩沖片內(nèi)能，從而起到緩沖作用。

為實現(xiàn)折疊翼展開之后的緩沖要求，需要確定合理的金屬橡膠的屬性參數(shù)。在折疊翼剛性體動力學(xué)仿真過程中，翼面展開到位時的減速緩沖是通過剛性體翼面與緩沖片之間定義的接觸碰撞來實現(xiàn)的，緩沖效果主要受到ADAMS接觸參數(shù)的影響，不能夠得到緩沖裝置的材料屬性，因此需要一種方法，將ADAMS的接觸參數(shù)（碰撞剛度k）與緩沖材料屬性參數(shù)（彈性模量E）建立一種關(guān)系，以確定合適的材料參數(shù)。

金屬橡膠緩沖片宏觀彈性模量表現(xiàn)為非線性，其參數(shù)更多需通過試驗測得。本文采用準(zhǔn)靜態(tài)接觸有限元分析與ADAMS仿真結(jié)合的方法，首先進行線性假設(shè)，將線性彈性模量設(shè)為優(yōu)化變量，確定線彈性參數(shù)；得到合理的金屬橡膠緩沖裝置屬性參數(shù)彈性模量E。再根據(jù)能量等效方法，結(jié)合確定的線彈性參數(shù)、金屬橡膠性質(zhì)和試驗，設(shè)計出符合指標(biāo)的金屬橡膠緩沖裝置。因此，金屬橡膠宏觀彈性模量在動力學(xué)仿真時按照線彈性Elinear計算，緩沖片設(shè)計時按照緩沖片從0壓縮至Δlstandard時能量等效原則進行非線性彈性模量Enon設(shè)計。

1 準(zhǔn)靜態(tài)接觸非線性有限元分析

1.1 接觸碰撞定義

在ADAMS/View中有兩種計算接觸力的方法，一種是恢復(fù)系數(shù)法；另一種是沖擊函數(shù)法。翼面與緩沖裝置的接觸碰撞屬沖擊接觸，在折疊翼虛擬樣機建模中的接觸定義中，采用沖擊函數(shù)法。沖擊函數(shù)法是根據(jù)IMPACT函數(shù)來計算兩個構(gòu)件之間的接觸力，接觸力由兩部分組成：一個是由于兩個構(gòu)件之間的相互切入而產(chǎn)生的彈性力；另一個是由相對速度產(chǎn)生的阻尼力。

在ADAMS的接觸碰撞函數(shù)中，定義的IMPACT碰撞力函數(shù)形式為：

式中，q0—兩個要碰撞物體的初始距離；

q—兩物體碰撞過程中的實際距離；

dq/dt—兩個物體間距離隨時間的變化率，即速度；

e—剛度貢獻指數(shù)；

cmax—最大阻尼系數(shù)；

d—全阻尼穿透值，決定阻尼力何時達到最大。為了防止碰撞過程中阻尼力的不連續(xù)，式中采用了step函數(shù)，其形式為step(x，x0，h0，x，h)

3）材料屬性定義

由于金屬橡膠材料的泊松比υ≤0.02，對緩沖片性能的影響很小，這里取0.02。

4）定義接觸

接觸設(shè)置里定義緩沖裝置為變形體，定義接觸曲面為剛性面。對緩沖裝置施加約束：圓孔底部約束6個自由度；約束內(nèi)側(cè)表面法向上的自由度，如圖1所示。對剛性面定義運動方式為準(zhǔn)靜態(tài)的繞軸轉(zhuǎn)動，直至緩沖裝置壓縮2mm，轉(zhuǎn)軸相對位置與實際情況相同，如圖2所示。

式中：a=h1-h0;Δ=(x-x0)/(x1-x0)。

1.2 接觸有限元分析

基于Patran進行建模并使用Marc求解，對緩沖裝置進行接觸分析，具體步驟如下：

1）創(chuàng)建幾何模型

在Patran中創(chuàng)建新數(shù)據(jù)庫文件，為保證緩沖裝置接觸分析的精度與幾何相似度，直接將緩沖裝置與翼面CATIA幾何模型導(dǎo)入。通過對翼面幾何模型的編輯，建立與翼面接觸面相同的曲面，然后將翼面實體模型刪除。

2）劃分有限元網(wǎng)格

采用手動劃分網(wǎng)格的方式，對緩沖裝置采用六面體網(wǎng)格劃分網(wǎng)格。

5）分析求解

分別對彈性模量E=150，300，400MPa時的模型進行接觸分析，輸出壓縮2mm過程中緩沖裝置的各節(jié)點上的位移變形和接觸力。400MPa情況下的位移變形和接觸受力分別見圖3和圖4。

2 緩沖裝置參數(shù)等效計算

2.1 彈性模量E與接觸力的關(guān)系

根據(jù)有限元分析計算，設(shè)接觸面從接觸緩沖裝置到壓縮2mm過程中轉(zhuǎn)動角度為φ，將φ九等分為9個角度，每個角度為φ，計算不同彈性模量下接觸面每轉(zhuǎn)過角度φ時的接觸力，記為角度1,2,…,9對應(yīng)的接觸力，如表1所示。

表1 不同彈性模量下接觸過程中的接觸力(N)

對于線彈性材料，其接觸力與彈性模量成正比。對300MPa情況下轉(zhuǎn)動角度與接觸力關(guān)系曲線進行擬合，得到力與角度之間的關(guān)系：

因此對于任意彈性模量E，對應(yīng)的接觸力：

圖5為采用Patran/Marc計算300Mpa情況下接觸力和角度關(guān)系后，擬合得到的曲線。

2.2 碰撞剛度k與接觸力的關(guān)系

在ADAMS中對翼面和緩沖裝置之間的準(zhǔn)靜態(tài)接觸過程進行分析，定義緩沖裝置固定，翼面繞軸旋轉(zhuǎn)，如圖6。剛度貢獻指數(shù)取e=1，最大阻尼cmax取k=0，全阻尼穿透值取d=2mm。分別取k=10000，15000，20000，同樣計算翼面轉(zhuǎn)動9個角度φ達到壓縮2mm過程中的接觸力，結(jié)果如表2所示。

由表中的數(shù)據(jù)得出接觸力與碰撞剛度k也成正比。k=20000時，對接觸力與轉(zhuǎn)動角度關(guān)系曲線進行擬合，得到接觸力與角度關(guān)系式為：

表2 不同參數(shù)k下接觸過程中的接觸力(N)

同樣對于任意的系數(shù)k有：

圖7為Adams計算k=20000情況下接觸力和角度關(guān)系，并擬合得到的曲線。

2.3 緩沖裝置材料彈性模量的確定

通過Patran&Marc分析和ADAMS仿真，分別得到了彈性模量E與接觸力、碰撞剛度k與接觸力之間的關(guān)系，以角度為9時的接觸力為標(biāo)準(zhǔn)，即壓縮量達到2mm時，令Fe=Fk，若確定k的值，可以根據(jù)式（4）和式（6）求出相應(yīng)的緩沖裝置等效彈性模量E。例如假設(shè)根據(jù)要求確定出k=20000，根據(jù)式（6）可以得到Fk=53981N，再利用式（4）即可反解出E= 235MPa，此時Fk與Fe對比如圖8。

3 緩沖裝置參數(shù)計算

分別對作動筒峰值壓強為8MPa、9MPa、10MPa三種工況改變不同的k值進行仿真計算，通過分析翼面展開到位時的緩沖效果：緩沖片最大壓縮量、速度下降比等因素得出，取緩沖片的最大壓縮量在2.7mm到2.8mm之間為最優(yōu)，確定碰撞參數(shù)k后，再根據(jù)第二節(jié)計算緩沖裝置參數(shù)的理論方法，求出三種工況下對應(yīng)的緩沖裝置彈性模量E。表3給出了三種工況下的計算結(jié)果。

表3 三種工況下仿真計算結(jié)果

4 結(jié)論

根據(jù)各工況下的計算結(jié)果與處理，可得：

1）通過有限元非線性接觸計算壓縮位移-載荷曲線與剛體動力學(xué)接觸計算位移-載荷曲線比對和迭代修正剛性體接觸參數(shù)，獲得的剛性體動力學(xué)接觸參數(shù)可用于模擬緩沖裝置壓縮過程。

2）利用該緩沖參數(shù)確定方法選取的緩沖材料參數(shù)，能使緩沖片最大壓縮量在2.2mm到3.05mm之間，緩沖碰撞壓縮2mm后，動能衰減超過50%，且最大壓縮量和最終壓縮量也均能滿足研制要求。

[1]白爭鋒,趙陽，等.含間隙太陽翼展開過程碰撞動力學(xué)研究.強度與環(huán)境，2008，35（4）.

[2]夏宇宏,姜洪源，等.金屬橡膠隔振器抗沖擊性能研究.振動與沖擊，2009,28(1).

[3]陳萌,張新訪，等.對接機構(gòu)隔框鎖的接觸碰撞動力學(xué)研究及仿真.宇航學(xué)報，2002，23(6).

[4]郭寶亭,朱梓根，等.金屬橡膠材料的理論模型研究[J].航空動力學(xué)報,2004,19(3):314-319.

[5]李玉亮，田偉，等.折疊翼展開到位碰撞特性分析.南昌大學(xué)學(xué)報工科版，2012,34(3).

[6]李莉,任茶仙，等.折疊翼機構(gòu)展開動力學(xué)仿真及優(yōu)化.強度與環(huán)境，2007,34(1).

>>>作者簡介

李玉亮，男，1981年3月出生，2003年畢業(yè)于西北工業(yè)大學(xué)，高級工程師，現(xiàn)從事飛行器設(shè)計工作。

付凡，男，1974年10月出生，1999年畢業(yè)于空軍工程學(xué)院，碩士，工程師，現(xiàn)從事飛行器設(shè)計及科研管理工作。

Analysis on the parameters of buffer device for folded-wing

Li Yuliang1，F(xiàn)u Fan2
(1.AVIC-HONGDU，Nanchang,Jiangxi，330024；2.Air Force Military Representative Office in Jiangxi,Nanchang，Jiangxi，330024)

After researching for overall characteristics of metal rubber buffer analysis,this thesis defines the equivalent elastic material parameters of the metal rubber buffer by combining the quasi-static contact finite element analysis to ADAMS dynamics simulation and by merging the rigid body dynamics contact collision rigidity with finite element contact rigidity.Ultimately,the impact load under different driving forces and metal rubbers combination is obtained when the folded-wing deploys to expected position,which provides reference for design of deploy mechanism and buffer engineering.

Folded-wing;Metal rubber;Cushion;Contact analysis;Dynamics simulation

2017-02-18）