白昀博，許 鋒

(南京航空航天大學(xué) 飛行器先進(jìn)設(shè)計技術(shù)國防重點學(xué)科實驗室，江蘇 南京 210016)

基于LMS Virtual.Lab的一種折疊機(jī)翼動力學(xué)仿真

白昀博，許 鋒

(南京航空航天大學(xué) 飛行器先進(jìn)設(shè)計技術(shù)國防重點學(xué)科實驗室，江蘇 南京 210016)

對艦載機(jī)折疊機(jī)翼的折疊及展開運動過程進(jìn)行動力學(xué)仿真研究，提出了一套基于LMS Virtual.Lab的折疊機(jī)翼動力學(xué)仿真建模方法，實現(xiàn)了機(jī)構(gòu)的動力學(xué)仿真計算，得到了折疊機(jī)翼作動筒載荷、翼稍運動包絡(luò)面、結(jié)構(gòu)應(yīng)力水平等動力學(xué)指標(biāo)。

折疊翼；艦載機(jī)；動力學(xué)仿真；LMS Virtual.Lab

艦載機(jī)是航空母艦的主要作戰(zhàn)武器，艦載機(jī)的數(shù)量標(biāo)志了其戰(zhàn)斗能力的強(qiáng)弱，這使得折疊機(jī)翼的研究成為了艦載機(jī)研發(fā)的重要組成部分。對于艦載機(jī)折疊翼的研究，國外起步較早，有著較為成熟的基礎(chǔ)，而國內(nèi)對于艦載機(jī)折疊翼的研究剛剛起步，具有較大的研究意義。本文將艦載機(jī)折疊翼的動力學(xué)仿真作為研究內(nèi)容，基于仿真平臺LMS Virtual.Lab建立折疊翼的動力學(xué)模型，著重對機(jī)構(gòu)作動筒載荷、翼稍運動包絡(luò)面及結(jié)構(gòu)應(yīng)力水平等指標(biāo)進(jìn)行了仿真計算，得出相關(guān)結(jié)果，為折疊翼設(shè)計提供依據(jù)。

1 總體方法

1.1仿真系統(tǒng)

LMS Virtual.Lab Motion基于計算多體動力學(xué)建模理論及計算方法，是專門為仿真多體機(jī)械系統(tǒng)的真實運動和載荷所設(shè)計的仿真平臺。它為快速創(chuàng)建和改進(jìn)多體模型提供了有效的方法，集成了應(yīng)用廣泛的CAD建模工具CATIA，可以實現(xiàn)建模分析一體化，可重復(fù)使用CAD和有限元模型，并快速反復(fù)模擬不同設(shè)計方案。這使得工程師可以在早期開發(fā)階段靈活地進(jìn)行動力學(xué)及運動學(xué)研究，并在后期階段結(jié)合實驗數(shù)據(jù)進(jìn)行綜合設(shè)計評估[1]。用戶可以直接通過軟件交互界面為機(jī)械系統(tǒng)輸入基本參數(shù)、建立復(fù)雜運動學(xué)和動力學(xué)模型，并利用軟件內(nèi)嵌的數(shù)學(xué)模型計算方法及樹枝計分方法進(jìn)行程式化求解，得到運動學(xué)規(guī)律及動力學(xué)響應(yīng)。LMS Virtual.Lab Motion同時還提供了可視化的后處理方法，采用動畫顯示、圖表、曲線等方式使得仿真數(shù)據(jù)更加直觀。此外，LMS Virtual.Lab Motion可以幫助用戶分析復(fù)雜機(jī)械系統(tǒng)的真實性能，為耐久性和振動噪聲分析及結(jié)構(gòu)分析提供精確的載荷[2]。

1.2仿真流程

在機(jī)翼折疊過程中，除結(jié)構(gòu)自重等內(nèi)部載荷外,最主要的外部載荷來自于甲板風(fēng)造成的氣動載荷，本文將來流與機(jī)身軸線夾角30°、由機(jī)頭指向機(jī)尾、流速31m/s(約0.1馬赫)的甲板風(fēng)設(shè)定為典型工況進(jìn)行仿真分析。

根據(jù)仿真過程，折疊機(jī)翼動力學(xué)剛?cè)狁詈辖＜胺抡媪鞒倘鐖D1所示。折疊機(jī)翼剛?cè)狁詈夏Ｐ偷慕⒅饕譃閯傮w模型的建立及有限元柔性模型的建立。剛體動力學(xué)模型是整個仿真研究的基礎(chǔ)，將折疊翼的CAD模型導(dǎo)入仿真軟件后對其進(jìn)行簡化、生成仿真分析實體，再添加運動副、約束及載荷等生成剛體分析模型。然后將需要柔性化的部件進(jìn)行有限元建模，之后在LMS Virtual.Lab Motion仿真平臺中將需要柔性化的部件用柔性體替換得到剛?cè)狁詈夏Ｐ?。最后進(jìn)行整個折疊機(jī)翼系統(tǒng)的動力學(xué)仿真，從而得到其動力學(xué)及運動學(xué)特征。需要注意的是，在建立有限元柔性體模型的過程中，需要在運動副及載荷的作用點添加MPC作為邊界條件，用以模擬系統(tǒng)實際約束關(guān)系，否則在剛?cè)狁詈戏抡鏁r將無法計算[3-4]。

2 仿真建模

2.1外部載荷計算

氣動載荷作為最主要的外部載荷是本文主要的研究對象，機(jī)翼在折疊運動過程中外翼受到的氣動載荷大小及方向都在不斷變化，故本文采用離散化的方法將外翼繞鉸鏈軸旋轉(zhuǎn)過程離散為多個穩(wěn)定的獨立狀態(tài)進(jìn)行氣動載荷的求解，然后再使用最小二乘法將其擬合生成外翼氣動載荷曲線。離散點的選取原則為載荷變化率越大取樣越密集，載荷變化率較小的區(qū)域取樣可相對減少。在專業(yè)流體分析軟件Fluent中計算得到的折疊過程中外翼阻、升力系數(shù)變化規(guī)律見表1。

然后將離散數(shù)據(jù)通過MATLAB進(jìn)行擬合，并根據(jù)機(jī)翼參數(shù)忽略翼面切向力,最終得到外翼氣動載荷關(guān)于旋轉(zhuǎn)角度的多項式函數(shù)：

F(x)=1.716×10-6×x5-0.000 893 3×x4+0.167 2×x3-14.32×x2+563.9×x-2 554.0式中：x為機(jī)翼轉(zhuǎn)角，°；F(x)為機(jī)翼壓力，N。氣動載荷正方向為翼面正法向。

2.2剛體建模

在LMSVirtual.Lab動力學(xué)仿真軟件Motion模塊中建立某型折疊機(jī)翼的剛體模型，作為其后各項分析的基礎(chǔ)。飛機(jī)折疊機(jī)翼結(jié)構(gòu)部件很多，各種結(jié)構(gòu)之間運動關(guān)系以及受力關(guān)系較為復(fù)雜。因此，在利用LMSVirtual.Lab進(jìn)行機(jī)翼折疊運動仿真分析時有必要對模型進(jìn)行一系列簡化，這樣既能夠降低多體系統(tǒng)仿真時的出錯概率，減少模型調(diào)試的時間，又能夠在保證分析精度的前提下提高分析計算的效率。

在完成以上模型簡化和準(zhǔn)備的基礎(chǔ)之上，獲得了簡化后的折疊機(jī)構(gòu)CAD模型，在LMSVirtual.Lab中進(jìn)行分析文件的建立。依次將各部件的CATIA模型通過“BodyFromExistingPart”命令導(dǎo)入到軟件Motion模塊中，生成分析實體，建立分析文件并生成多體系統(tǒng)中的分析體。然后根據(jù)折疊機(jī)構(gòu)各個部件之間的實際運動情況添加相應(yīng)的運動副，如定義圓柱副(CylindricalJoint)、固定副(BracketJoint)、轉(zhuǎn)動副(RevoluteJoint)及萬向節(jié)副(UniversalJoint)。創(chuàng)建運動副后，軟件會默認(rèn)在分析體中添加相應(yīng)的裝配約束，自動將機(jī)翼各個部件裝配成完整的折疊機(jī)構(gòu)整體，通過運動副以及裝配關(guān)系的添加，得到折疊機(jī)翼系統(tǒng)的剛體模型如圖2所示。同時為盡可能準(zhǔn)確地模擬折疊機(jī)構(gòu)實際運動情況，在拉桿與作動筒之間的圓柱副上施加一個運動副速度驅(qū)動“JointVelocityDriver”，使作動筒相對于拉桿在公共軸線的方向上以規(guī)定的相對速度運動[5]。

由于軟件限制及為了兼顧仿真計算精度和計算效率，本文將外機(jī)翼氣動載荷的分布力等效為作用在若干個翼肋剛心處的集中力。通過“ScalarForceExpression”命令在翼肋剛心處添加隨外翼運動并保持垂直翼面方向的函數(shù)作用力，設(shè)置為“ACTION_ONLY”,即只有作用力作用于外翼。

2.3剛?cè)狁詈辖?/p>

由于LMSVirtual.Lab軟件自身的網(wǎng)格劃分功能并不強(qiáng)，對于由蒙皮、翼肋等薄片結(jié)構(gòu)的部件無法進(jìn)行可用的有限元網(wǎng)格劃分，本文采用有限元前處理軟件Patran，利用其強(qiáng)大的有限元建模功能，使用二維網(wǎng)格對外機(jī)翼進(jìn)行有限元建模，通過對二維網(wǎng)格屬性的添加使有限元網(wǎng)格匹配零件實際形狀，同時使用HyperMesh更為優(yōu)秀的三維網(wǎng)格功能對鏈接耳片等形狀較為規(guī)則的零件進(jìn)行有限元網(wǎng)格的建立[6]。按照折疊機(jī)翼各部件實際運動關(guān)系及約束，在有限元前處理軟件中添加REB2類型的MPC(多點約束)點來模擬折疊機(jī)構(gòu)各個部件位置、運動關(guān)系以作為邊界條件，從而建立起折疊機(jī)翼的有限元模型。

將折疊機(jī)翼各部件的有限元模型導(dǎo)入至仿真平臺LMSVirtual.Lab中，利用軟件中的“MakeFlexiblewithExistingData”命令將剛體模型中的對應(yīng)分析實體替換為柔性體。柔性體生成后，LMSVirtual.Lab軟件界面中外機(jī)翼模型將由柔性體分析文件中的網(wǎng)格模型替換。在外機(jī)翼柔性體模型與剛體模型成功關(guān)聯(lián)后，上文所述的MPC點會被自動檢測到，通過“RigidSpiderbyPicking”命令建立“蜘蛛網(wǎng)格”將節(jié)點與網(wǎng)格的實際受力部分相連作為柔性體的約束條件。然后通過“NastranCraig-BamptonSolution”使用修正的Craig-Bampton模態(tài)綜合法對各柔性體逐一進(jìn)行模態(tài)分析生成中性文件，并激活相應(yīng)的模態(tài)，從而完成折疊機(jī)翼的剛?cè)狁詈辖7-8]，如圖3所示。

3 仿真與分析

3.1仿真參數(shù)

仿真過程利用建立的折疊機(jī)翼多體動力學(xué)剛?cè)狁詈夏Ｐ?，采用作動筒與拉桿公共軸線方向相對速度驅(qū)動，并在開始與結(jié)束階段設(shè)定一定的緩沖。折疊運動耗時19.4s，作動筒/拉桿相對速度驅(qū)動曲線如圖4所示。在分析之前必須對求解器進(jìn)行設(shè)置，在整個系統(tǒng)運動的動態(tài)分析計算過程中，為加快計算速度采用顯式算法計算力雅可比矩陣，在時域內(nèi)優(yōu)先采用BDF(向后差分法)進(jìn)行積分，如不收斂則采用PECE(預(yù)估校正法)進(jìn)行積分求解，計算的步長設(shè)置為1E-9s。由于柔性體的存在，整個折疊機(jī)翼系統(tǒng)在計算初始時刻存在落振，故在機(jī)構(gòu)開始折疊及展開前各設(shè)置足夠長的緩沖時間(70s)，使折疊翼機(jī)構(gòu)進(jìn)入穩(wěn)定狀態(tài)，各參數(shù)平穩(wěn)后機(jī)構(gòu)開始運動，以便更好地觀測和分析，設(shè)置仿真時間180s。

3.2仿真結(jié)果分析

經(jīng)過仿真分析，得到折疊機(jī)構(gòu)作動筒在折疊與展開運動過程中載荷譜線，如圖5所示。圖中載荷正方向為拉力，負(fù)方向為壓力。氣動載荷對外翼的作用由不斷減小的升力連續(xù)變化為增大的壓力，然后基本穩(wěn)定，這個過程對折疊運動的作用由阻力變?yōu)閯恿?。由于結(jié)構(gòu)自重與氣動載荷的耦合作用，作動筒軸向載荷譜線的變化規(guī)律為拉力載荷先增大后減小，外翼存在翼展長的特點，氣動載荷相對轉(zhuǎn)軸會產(chǎn)生較大的力矩，這使得氣動載荷在折疊運動的后段提供了較大的動力使外翼具有加速運動的趨勢，與此同時重力的阻力力矩也在增大，但無法抵消氣動載荷的作用，為使外翼平衡，作動筒在整個折疊運動過程中均提供拉力，即受到拉載荷。展開過程作動筒載荷譜線與折疊過程基本對稱。折疊過程最大拉載荷49 421.625N，展開過程最大拉載荷49 380.246N。

以外翼翼稍下側(cè)邊緣中點作為監(jiān)測點，仿真得到外翼在折疊及展開運動中翼稍形變包絡(luò)面，如圖6所示，X方向為機(jī)翼弦向向后、Y方向為翼面法向向上。由于存在風(fēng)載的升力作用，折疊初始位置及展開終止位置的Y方向形變位移量均約為-3.8mm。此工況下翼稍Y方向形變位移量在翼稍經(jīng)過最低點時出現(xiàn)最大值，約-14.5mm，折疊終止位置X，Y方向翼稍偏移量分別為-0.466mm、-8.754mm。

對于機(jī)構(gòu)運動過程中承載最大的雙耳鉸鏈軸部位零件仿真分析后得到了其應(yīng)力分布情況，固連于內(nèi)翼的鉸鏈雙耳接頭最大應(yīng)力為450MPa，固連于外翼的上、下耳片最大應(yīng)力分別為391MPa、273MPa。

經(jīng)過仿真計算，得到了折疊機(jī)翼機(jī)構(gòu)的作動筒載荷曲線、外翼翼稍運動包絡(luò)面及鉸鏈軸位置部件的結(jié)構(gòu)應(yīng)力水平，與機(jī)構(gòu)實際運動情況對比分析，其結(jié)果真實可靠。

4 結(jié)束語

本文基于多體動力學(xué)仿真平臺LMSVirtual.LabMotion建立了一套折疊機(jī)翼的多剛體動力學(xué)建模方法，對該折疊翼機(jī)構(gòu)在折疊及展開過程中的動力學(xué)特性進(jìn)行了仿真分析，填補了國內(nèi)對艦載機(jī)折疊翼動力學(xué)仿真研究的空白。同時，由于條件限制,使用運動副來模擬了機(jī)構(gòu)的運動規(guī)律，忽略了機(jī)械結(jié)構(gòu)在運動過程中必然存在的結(jié)構(gòu)間隙所產(chǎn)生的非線性影響，這將是后續(xù)研究的重點，以便更精確地對折疊機(jī)翼的動力學(xué)特性進(jìn)行仿真分析。

[1] 時峰．LMSVirtual.LabMotion新一代多體動力學(xué)軟件[J]．CAD/CAM與制造業(yè)信息化，2008(12)：52-56.

[2]LMS公司.LMSVirtual.Lab第11版新功能[J].CAD/CAM與制造業(yè)信息化, 2012 (12)：98-99.

[3] 陳晏，余躍慶，蘇麗穎，等．基于LMSVirtual.Lab和ANSYS的剛?cè)狁詈蠙C(jī)器人仿真系統(tǒng)[J]．輕工機(jī)械，2007，25(1):46-49.

[4] 崔飛，馬東立．基于LMSVirtual.Lab的起落架動態(tài)性能仿真分析[J]．計算機(jī)輔助工程，2012，21(2): 25-29.

[5] 郭小良，裴錦華，楊忠清，等．無人機(jī)折疊機(jī)翼展開運動特性研究[J]．南京航空航天大學(xué)學(xué)報，2006(4):438-441.

[6] 徐孝武，張煒．折疊機(jī)翼變體飛機(jī)的動力學(xué)建模與分析[J]．西北工業(yè)大學(xué)學(xué)報，2012(5):681-688.

[7] 鄧峰巖，和興鎖，張娟，等．修正的Craig-Bampton方法在多體系統(tǒng)動力學(xué)建模中的應(yīng)用[J]．機(jī)械設(shè)計，2004，21(3):41-43.

[8] 高星斗，畢世華，陳陣．基于改進(jìn)Craig-Bampton法的導(dǎo)彈發(fā)射過程多柔體動力學(xué)研究[J]．固體火箭技術(shù)，2011，34(5)：559-563.

The dynamic simulation of a folding wing based on LMS Virtual.lab

BAI Yunbo, XU Feng

(Nanjing University of Aeronautics and Astronautics, Jiangsu Nanjing, 210016, China)

It introduces the dynamic simulation for folding wings of aircraft in folding and unfolding process, proposes an effective dynamic simulation modeling method of folding wings based on LMS Virtual.lab, realizes the dynamic simulation, obtains a series of dynamic parameters including the cylinder load and the wingtips' moving envelope surface with the structure stress.

folding wing; ship based aircraft; dynamic simulation; LMS Virtual.lab

10.3969/j.issn.2095-509X.2015.01.007

2014-12-16

白昀博(1989—)，男，遼寧阜新人，南京航空航天大學(xué)碩士研究生，主要研究方向為飛機(jī)結(jié)構(gòu)動力學(xué)仿真。

V224

A

2095-509X(2015)01-0029-04