陳振中，張 弛，王璐璐，王 震

(1.沈陽航空航天大學(xué) 民用航空學(xué)院，沈陽 110136；2.中國南方航空股份有限公司 沈陽維修基地，沈陽 110000)

0 引言

供油管道是飛機輔助動力裝置(APU)中的重要零部件，工作狀態(tài)下管道中燃油的流動會產(chǎn)生振動，進而引發(fā)管道螺栓連接件振動，造成螺母與螺栓之間的螺紋連接產(chǎn)生松動，導(dǎo)致螺母中螺紋的疲勞損傷，降低管道螺母的使用壽命。

祖衛(wèi)華等人運用ANSYS軟件對管壁板結(jié)構(gòu)進行特性分析，通過增加管道約束的方式，提高整體結(jié)構(gòu)各階自振頻率，有效地避開了工作頻率范圍。李少靜等人在流固耦合作用下對管道進行模態(tài)分析，得到管道固有頻率和振型圖。康力等人對航空發(fā)動機多個位置的外管路進行振動響應(yīng)分析，得出各個部位管道各異振動特征。呂彪等人研究航空發(fā)動機供油管道裂紋，認為管道裂紋為高周疲勞裂紋。國內(nèi)外學(xué)者對航空發(fā)動機上的管道振動特性研究較多，對于飛機APU供油管道研究較少。APU在實際使用中經(jīng)常發(fā)生供油管路振動，進而引起供油管路接頭螺紋損傷，因此對APU供油管道進行模態(tài)分析，給出供油管道的固有頻率及振型，對研究管道連接接口疲勞損傷具有重要價值。

1 三維建模及數(shù)值模擬

1.1 三維建模

本文利用SolidWorks軟件對APU供油管道進行1：1等比例建模，導(dǎo)管三維模型如圖1所示。管道的材料為鈦合金Ti-3Al-2.5V，密度為4 620 kg·m。在ANSYS DM模塊中建立流體模型，選中管道內(nèi)壁，在Fluent模塊中利用“fill”命令生成流體，設(shè)定材質(zhì)為液體燃油Gasoil Liquid，密度為770 kg·m，流體模型如圖2所示。

圖1 導(dǎo)管三維模型

圖2 流體模型

APU中供油管道種類多種多樣，本文利用SolidWorks軟件對其中一種型號進行等比例三維實體建模，然后運用Workbench軟件對管道進行仿真分析，最終得出管道的振型特點。

1.2 流固耦合及模態(tài)分析

流固耦合力學(xué)是研究變形固體在流場作用下的各種行為以及固體位形對流場影響這二者相互作用的一門科學(xué)。為了計算燃油與管道之間的相互作用，本文采用流固耦合進行供油管道的模態(tài)分析，因為管道變形很小，所以使用Fluent模塊與靜力學(xué)模塊Static Structural連結(jié)組合成的單向流固耦合模塊，在Fluent模塊中計算流體模型在流經(jīng)管道產(chǎn)生的壓力，然后在靜力學(xué)模塊Static Structural中計算管道受流體影響產(chǎn)生的應(yīng)力載荷。

1.2.1 流體模型數(shù)據(jù)計算

采用掃掠法sweep method對流體模型進行網(wǎng)格劃分，設(shè)置流體進出口，對進口面設(shè)定8層邊界層，網(wǎng)格類型設(shè)置為三棱柱，并對進口面設(shè)定網(wǎng)格大小為1 mm，共生成214 901個節(jié)點，277 140個單元，流體網(wǎng)格劃分如圖3所示。根據(jù)APU供油管道實際維修數(shù)據(jù)，在Fluent模塊中設(shè)置管道內(nèi)壓強為1.379 MPa，燃油流量0.025 2 kg·s，流速0.014 2 m/s。由于管道出口唯一，根據(jù)進口流速和壓強設(shè)置出口邊界條件為自由流出，流出量權(quán)重設(shè)置為1。計算得到殘差圖如圖4所示，從殘差圖可以看出各參數(shù)收斂曲線隨計算步數(shù)的增加而降低，最后趨于平緩，說明流體計算收斂，可導(dǎo)入管道進行下一步計算。

圖3 流體網(wǎng)格劃分

圖4 殘差圖

通過CFD-Post軟件得到流體流速矢量圖如圖5所示，從矢量圖可以看出流體流速在管道中間會更快，并且在管道彎曲處流速增加。

圖5 流體流速矢量圖

1.2.2 靜力學(xué)模塊計算及模態(tài)分析

在模型中對供油管道采用六面體網(wǎng)格劃分，設(shè)置網(wǎng)格單元大小為1 mm，對于管道彎曲部分進行面映射網(wǎng)格劃分，共生成144 900個節(jié)點，26 310個單元，管道網(wǎng)格劃分如圖6所示。將Fluent中流體對管道產(chǎn)生的載荷設(shè)置壓力導(dǎo)出，設(shè)置邊界條件：由于管道較短，管道之間只有在接頭處的螺母固定，所以在接頭處設(shè)置為固定支點，然后進行靜力學(xué)分析，得出管道在流固耦合狀態(tài)下的載荷分布。

圖6 管道網(wǎng)格劃分

在靜力學(xué)模塊中進行模態(tài)分析，得到管道固有頻率如表1所示，管道前六階振型圖如圖7所示。從表1中可以看出前三階固有頻率平穩(wěn)增長，第四階到第六階固有頻率出現(xiàn)跳躍式增長。最大變形主要在管道中間彎曲處，其中第三階和第五階變形最大；第一階振型中管道主要沿著軸方向變形，第二階振型中管道主要沿著軸方向變形，第三階振型中管道主要沿著軸和軸方向變形，四到六級階振型管道變形并無明顯規(guī)律。

圖7 前六階振型圖

表1 固有頻率

2 結(jié)論

本文通過SolidWorks軟件對某型APU供油導(dǎo)管進行三維建模，在ANSYS Workbench軟件中進行流固耦合分析，得到供油狀態(tài)下管道的固有頻率和前六階振型。研究結(jié)果表明第一階振型中管道主要沿著軸方向變形，第二階振型中管道主要沿著軸方向變形，第三階振型中管道主要沿著軸和軸方向變形，四到六級階振型管道變形并無明顯規(guī)律；管道變形主要發(fā)生在管道中間彎曲處，其中第三階和第五階振型變形最大。