賀迪華，羅衛(wèi)東,2*，黃琳森，劉忠

?

某軌姿控管路系統(tǒng)模擬分析及優(yōu)化

賀迪華1，羅衛(wèi)東1,2*，黃琳森1，劉忠3

（1.貴州大學(xué)機(jī)械工程學(xué)院，貴州 貴陽 550025；2.貴州大學(xué)明德學(xué)院，貴州 貴陽 550025；3.貴州航天朝陽科技有限責(zé)任公司，貴州 遵義 563000）

根據(jù)某軌姿控系統(tǒng)中關(guān)鍵組件的固定位置進(jìn)行管路布置，并通過ANSYS Workbench對管路和管路接頭進(jìn)行了仿真分析，針對分析結(jié)果、結(jié)合管路布置原則對原管路進(jìn)行改進(jìn)。通過與原管路對比發(fā)現(xiàn)：改進(jìn)管路模型的固有頻率避開了系統(tǒng)工作頻率910 Hz；在隨機(jī)振動(dòng)中方向上的等效應(yīng)力和位移分別減小了66.62 MPa和0.72878 mm且在、方向上的等效應(yīng)力和位移增加幅度較小；四通管路接頭使得管路中流量分配更加均勻且降低了將近10%的水擊壓力，有效地改善了軌姿控系統(tǒng)中管路的布置形式。

軌姿控管路系統(tǒng)；四通管路接頭；管路模型；仿真分析

軌姿控管路系統(tǒng)是管子、管子連接件、閥門等連接成用于輸送液體推進(jìn)劑的裝置，是典型的動(dòng)力支撐和維持系統(tǒng)。管路系統(tǒng)的結(jié)構(gòu)、布置形式以及工作品質(zhì)直接影響到結(jié)構(gòu)系統(tǒng)的正常工作[1]。管路系統(tǒng)設(shè)計(jì)和布置不合理導(dǎo)致工作系統(tǒng)失效的事故鮮有發(fā)生，嚴(yán)重的甚至造成管路破裂，給系統(tǒng)帶來巨大的經(jīng)濟(jì)損失。因此，對軌姿控管路系統(tǒng)的優(yōu)化改進(jìn)和分析研究有著重要的意義。

目前，管路系統(tǒng)的研究存在諸多方面。文獻(xiàn)[2]采用了一維有限元方法對管路中的液體瞬變流進(jìn)行計(jì)算。文獻(xiàn)[3-4]對液體火箭發(fā)動(dòng)機(jī)系統(tǒng)的充填過程動(dòng)態(tài)特性進(jìn)行了研究。文獻(xiàn)[5]采用有限差分格式的特征線方法，針對某姿控推進(jìn)系統(tǒng)發(fā)動(dòng)機(jī)真實(shí)氧化劑管路的關(guān)機(jī)瞬變特性進(jìn)行了數(shù)值研究，結(jié)果發(fā)現(xiàn)噴注前的最大瞬變壓力達(dá)到貯箱壓力的三倍。文獻(xiàn)[6]對液壓管路流固耦合線性化動(dòng)力學(xué)模型和非線性振動(dòng)模型的研究成果進(jìn)行分析，并對液壓管路被動(dòng)、主動(dòng)及半主動(dòng)振動(dòng)控制的研究進(jìn)展及研究成果進(jìn)行評述，提出了管路系統(tǒng)流固耦合振動(dòng)機(jī)理及振動(dòng)控制研究的發(fā)展趨勢。文獻(xiàn)[7]針對船舶管路系統(tǒng)抗沖擊性能仿真方法精度不足的問題，研究用于艦船管路系統(tǒng)的高精度抗沖擊仿真方法。文獻(xiàn)[8]針對原四通接頭在工作過程中存在的缺陷，對四通接頭進(jìn)行改進(jìn)設(shè)計(jì)，結(jié)果表明改進(jìn)后的四通接頭效果良好。

本文通過對四通接頭和原管路布置形式進(jìn)行了優(yōu)化改進(jìn)，將優(yōu)化的四通接頭運(yùn)用到改進(jìn)的管路中，最后對改進(jìn)的管路模型進(jìn)行模擬分析。

1 計(jì)算理論

1.1 湍流計(jì)算理論

標(biāo)準(zhǔn)－模型是目前應(yīng)用范圍最廣的兩方程湍流模型[9]，其方程為：

其中：

1.2 模態(tài)分析理論

模態(tài)分析是最為基礎(chǔ)的內(nèi)容，也是動(dòng)力學(xué)分析中必不可少的環(huán)節(jié)。工程中進(jìn)行模態(tài)分析的主要目的在于找出系統(tǒng)的固有頻率，從而避免可能引起的共振現(xiàn)象。

自由度系統(tǒng)自由振動(dòng)的方程為[10]：

式中：[]、[]和[]分別為系統(tǒng)質(zhì)量、阻尼和剛度矩陣。

因?yàn)槭剑?）為耦合方程，計(jì)算時(shí)需要對其進(jìn)行解耦處理。解耦方式如下：

令：

即將式（3）變化為：

將質(zhì)量、阻尼和剛度矩陣同時(shí)對角化，依據(jù)振型矩陣正交性關(guān)系，有：

將自由度系統(tǒng)自由振動(dòng)方程解耦成模態(tài)坐標(biāo)下個(gè)相互獨(dú)立的自由振動(dòng)方程組，第個(gè)方程為：

2 四通接頭計(jì)算與分析

2.1 四通接頭建模

運(yùn)用繪圖軟件Creo對四通接頭進(jìn)行建模，如圖1所示3種四通接頭模型。圖（a）為直通型流道空槽，圖（b）為半徑＝13 mm的球形流道空槽，圖（c）為半徑＝23 mm的球形流道空槽，接頭壁厚均為2 mm，管路壁厚均為1 mm，管路內(nèi)徑為10 mm。

2.2 流域模型及邊界條件

將三種四通接頭導(dǎo)入ICEM CFD軟件中抽取流域并劃分網(wǎng)格，劃分結(jié)果如圖2所示。

圖1 三種四通接頭模型

圖2 流域網(wǎng)格模型

根據(jù)軌姿控發(fā)動(dòng)機(jī)實(shí)際工作要求，設(shè)置三個(gè)出口邊界條件均為0.4 MPa；選用液體推進(jìn)劑MMH作為工質(zhì)，其密度為874.40 kg/m3，動(dòng)力粘度為0.775×10-3Pa·s[11]；計(jì)算3.0579 m/s、 6.1157 m/s、9.1736 m/s和12.2315 m/s四種工況下推進(jìn)劑入口速度對三個(gè)出口質(zhì)量流量分配的影響。

2.3 計(jì)算結(jié)果分析

圖3為三種四通接頭在入口速度6.1157 m/s下的速度矢量分布情況。

圖3 四通接頭速度矢量圖

在直通流域中，Outlet1和Outlet3附近存在明顯的渦流現(xiàn)象，入口速度增大時(shí)，渦流現(xiàn)象越加明顯，對四通接頭產(chǎn)生不利影響。在＝13 mm球形流域中，避免了上述現(xiàn)象的發(fā)生，但由于流體慣性作用，使得Outlet2的出口速度依然顯著高于Outlet1和Outlet3。在＝23 mm球形流域中，由于球形容積較大，有效地減緩了流體進(jìn)入球形域的這種慣性，因而使得各出口的速度基本相近。

通過對三種四通接頭在不同入口速度下的仿真分析，繪制出各四通接頭出口質(zhì)量流量百分率隨速度的變化情況，如圖4所示。

圖4 不同出口流量分配率與速度關(guān)系

從圖4中可以看出，＝23 mm球形流域各出口的質(zhì)量流量分配率基本相等，即說明三個(gè)出口同時(shí)工作時(shí)不會(huì)出現(xiàn)某個(gè)出口處的質(zhì)量流量滿足不了給定需求；此外，各出口質(zhì)量流量分配率與入口速度基本無關(guān)。

根據(jù)文獻(xiàn)[12]減小水擊的措施可知，當(dāng)管路串聯(lián)一個(gè)容積50 ml的集液腔后，系統(tǒng)管路的水擊壓力會(huì)降低10%左右。由于＝23 mm球形流域的體積為50.965 ml，所以在布置空間允許的條件下，盡量采用＝23 mm球形空槽四通接頭能夠顯著的降低水擊壓力。

3 原始管路計(jì)算與分析

3.1 管路系統(tǒng)原始模型

根據(jù)某軌姿控系統(tǒng)中關(guān)鍵組件的固定位置和發(fā)動(dòng)機(jī)的工作要求對管路系統(tǒng)進(jìn)行布置，如圖5所示。

圖5 原始管路結(jié)構(gòu)

3.2 模型網(wǎng)格及邊界條件

將管路系統(tǒng)模型導(dǎo)入網(wǎng)格劃分軟件，管路采用掃掠方式而接頭采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行網(wǎng)格劃分。管路系統(tǒng)的邊界條件根據(jù)空間結(jié)構(gòu)形式進(jìn)行合理的布置，選取5個(gè)固定支點(diǎn)。管路的彈性模量＝204 GPa，密度＝7930 kg/m3，泊松比＝0.285，抗拉強(qiáng)度σ＝520 MPa，屈服強(qiáng)度σ＝206 MPa。網(wǎng)格模型和邊界條件如圖6所示。

圖6 模型網(wǎng)格及邊界條件

3.3 模型計(jì)算及結(jié)果分析

采用ANSYS Workbench對原始管路模型進(jìn)行仿真模擬，得出原始模型的前20階固有頻率值如表1所示。

表1 原始管路模態(tài)

根據(jù)軌姿控系統(tǒng)工作要求，系統(tǒng)的工作頻率在910 Hz附近。從表1中可以看出，原始管路第10階和第11階模態(tài)剛好在系統(tǒng)的工作頻率范圍附近，即容易產(chǎn)生共振現(xiàn)象。

在模態(tài)分析的基礎(chǔ)上，對各個(gè)方向施加載加速度功率譜密度曲線對其進(jìn)行隨機(jī)振動(dòng)分析，其曲線如圖7所示。

圖7 加速度功率譜密度曲線

根據(jù)上述曲線，計(jì)算隨機(jī)振動(dòng)下各個(gè)方向的最大應(yīng)力、位移值，結(jié)果如表2所示。由表2可知，、和方向上的最大應(yīng)力值均小于材料的屈服強(qiáng)度；但在方向上的位移值明顯偏大，并且在方向上的其他部件布置相對較緊密，容易與其他部件產(chǎn)生干涉，進(jìn)而破壞其結(jié)構(gòu)。

表2 各方向隨機(jī)振動(dòng)載荷下的最大應(yīng)力和位移（一）

4 改進(jìn)管路計(jì)算與分析

4.1 管路系統(tǒng)改進(jìn)模型

由于整個(gè)軌姿控系統(tǒng)中外圍和中間部位無組件占用空間位置且不影響管路的布置；遂針對原始管路系統(tǒng)存在的缺陷和產(chǎn)生的不利影響，并結(jié)合管路布置的原則對原始管路系統(tǒng)進(jìn)行如下改進(jìn)：

（1）考慮管路的補(bǔ)償，原始管路中的直管改成彎管布置；

（2）根據(jù)“在考慮管路補(bǔ)償?shù)那疤嵯?，管路彎得越少越好”的原則，原始管路中直管改成彎管不應(yīng)彎得太大；

（3）采用上述計(jì)算得到的＝23 mm球形空槽四通接頭替換直通槽接頭進(jìn)行管路連接；

（4）考慮三通接頭流量均勻分配原則，選用垂直端為入口、水平兩端為出口。

上述對原始管路系統(tǒng)具體改進(jìn)如圖8所示。

圖8 改進(jìn)管路結(jié)構(gòu)

4.2 模型網(wǎng)格及邊界條件

采用與原始模型相同方式對其進(jìn)行網(wǎng)格劃分；管路系統(tǒng)的邊界條件根據(jù)空間結(jié)構(gòu)形式進(jìn)行合理布置，選取7個(gè)固定支點(diǎn)，材料選用相同的材料，其網(wǎng)格模型和邊界條件如圖9所示。

圖9 模型網(wǎng)格及邊界條件

4.3 模型計(jì)算及結(jié)果分析

對改進(jìn)管路模型進(jìn)行模擬分析，得出改進(jìn)模型的前20階固有頻率值如表3所示。

表3 改進(jìn)管路模態(tài)

從表3中可以看出，改進(jìn)管路的固有頻率都遠(yuǎn)離了系統(tǒng)工作的頻率910 Hz，即很好的避免了共振現(xiàn)象發(fā)生。

同樣在模態(tài)分析的基礎(chǔ)上，采用上述加速度功率譜密度曲線對改進(jìn)模型進(jìn)行隨機(jī)振動(dòng)計(jì)算，得到各個(gè)方向上的最大應(yīng)力、位移值，具體如表4所示。

由表4可知，、和方向上的最大應(yīng)力值均小于材料的屈服強(qiáng)度，并且在方向上的最大位移為0.16502 mm，相比原始管路減小了0.72878 mm，效果顯著。

表4 各方向隨機(jī)振動(dòng)載荷下的最大應(yīng)力和位移（二）

5 分析結(jié)果對比

通過模態(tài)結(jié)果對比發(fā)現(xiàn)：原始管路模型的固有頻率剛好位于軌姿控系統(tǒng)工作頻率910 Hz附近，即容易引發(fā)共振現(xiàn)象；而改進(jìn)管路模型的固有頻率最接近軌姿控系統(tǒng)工作頻率的也只有974.33 Hz，距離系統(tǒng)工作頻率較遠(yuǎn)，很難引發(fā)共振現(xiàn)象。

通過隨機(jī)振動(dòng)分析對比發(fā)現(xiàn)：原始管路模型在各個(gè)方向上的最大應(yīng)力值均小于材料的屈服強(qiáng)度，但在方向上的位移值明顯偏大，很容易與其他部件發(fā)生干涉；而改進(jìn)管路模型在方向上的應(yīng)力和位移值相對原始模型分別減小了66.62 MPa和0.72878 mm，降低效果顯著，繼而避免了與其他部件發(fā)生干涉的現(xiàn)象，并且在和方向上的應(yīng)力和位移值相比原始管路模型增加幅度較小。

6 結(jié)論

通過上述3種四通接頭和2款管路模型研究分析，得到以下結(jié)論：

（1）在3種四通接頭中，＝23 mm的球形四通接頭各個(gè)出口的質(zhì)量流量基本相等，且能夠顯著的降低水擊壓力，但其受到布置空間的限制。

（2）改進(jìn)管路模型很好的避開了系統(tǒng)工作的頻率范圍，即避免了共振現(xiàn)象的發(fā)生。

（3）改進(jìn)管路模型減小了在方向上隨機(jī)振動(dòng)的應(yīng)力和位移值，有效的避免了管路系統(tǒng)與其他部件之間發(fā)生干涉。

[1]蘇爾皇. 液壓流體力學(xué)[M]. 北京：國防工業(yè)出版社，1982.

[2]程謀森，劉昆，張育林. 推進(jìn)劑供應(yīng)管路內(nèi)液體瞬變流一維有限元計(jì)算[J]. 推進(jìn)技術(shù)，2000，21（4）：12-15.

[3]程謀森，張育林. 航天器推進(jìn)系統(tǒng)管路充填過程動(dòng)態(tài)特征(I)理論模型與仿真結(jié)果[J]. 推進(jìn)技術(shù)，2000，21（2）：25-28.

[4]程謀森，張育林，姜春林，等.航天器推進(jìn)系統(tǒng)管路充填過程動(dòng)態(tài)特性(II)實(shí)驗(yàn)?zāi)M與結(jié)果評估[J]. 推進(jìn)技術(shù)，2000，21（3）：1-5.

[5]聶萬勝，陳新華，戴德海，夏鵬，莊逢辰. 姿控推進(jìn)系統(tǒng)發(fā)動(dòng)機(jī)關(guān)機(jī)的管路瞬變特性[J]. 推進(jìn)技術(shù)，2003，24（1）：6-8.

[6]權(quán)凌霄，孔祥東，俞濱，白歡歡. 液壓管路流固耦合振動(dòng)機(jī)理及控制研究現(xiàn)狀與發(fā)展[J]. 機(jī)械工程學(xué)報(bào)，2015，51（18）：175-183.

[7]于雯，嵇春艷，李曙生. 基于 ANSYS 的高精度管路系統(tǒng)抗沖擊仿真方法及試驗(yàn)驗(yàn)證[J]. 船舶工程，2017，39（1）：73-85.

[8]賈凱軍，馮光明，鄭艷棟，范瑋剛，李鵬舉. 混合管與充填軟管間[J]. 礦山機(jī)械，2010，38（20）：56-58.

[9]王福軍. 計(jì)算流體動(dòng)力學(xué)分析[M]. 北京：清華大學(xué)出版社，2004.

[10]杜平，胡夏夏. 管路系統(tǒng)的模態(tài)分析[J]. 輕工機(jī)械，2012，30（3）：29-31.

[11]張煒，鮑桐，周星. 火箭推進(jìn)劑[M]. 北京：國防工業(yè)出版社，2014.

[12]張崢岳，康乃全. 軌姿控液體火箭發(fā)動(dòng)機(jī)水擊仿真模擬[J]. 火箭推進(jìn)，2012，38（3）：12-16.

Simulation Analysis and Optimization of Pipeline System in an Orbit and Attitude Control System

HE Dihua1，LUO Weidong1,2，HUANG Linsen1，LIU Zhong3

( 1.School of Mechanical Engineering, Guizhou University, Guiyang 550025, China; 2.Mingde College of Guizhou University, Guiyang 550025, China; 3.Limited Liability Company of Guizhou Space Flight Sun Technology, Zunyi 563000, China )

According to the fixed position of key components in an orbit and attitude control system, using by ANSYS Workbench, the pipeline and pipe joints are analyzed and the original pipeline is enhanced by the following analysis results and the principle of pipeline arrangement. Compared with the original pipeline, this paper finds that the natural frequency of the improved pipeline model avoids the working frequency 910 Hz of the operating system. In the random vibration, the equivalent stress and displacement in thedirection are significantly reduced to 66.62 MPa and 0.72878 mm and the equivalent stress and increment of displacement in theanddirections are small. The flow distribution of the four-way pipe joint is more reasonable and it can reduce the nearly 10% water hammer pressure which can effectively improve the orbit and attitude control system layout of the pipeline.

the orbit and attitude control system；the four-way pipe joint；pipe model；simulation analysis

TH138.8+4

A

10.3969/j.issn.1006-0316.2018.01.003

1006－0316 (2018) 01－0011－07

2017－06－09

貴州省自主創(chuàng)新項(xiàng)目（KYKZ－001）；貴州省科技支撐計(jì)劃項(xiàng)目（黔科合支撐[2016]2319）

賀迪華（1992－），男，江西吉安人，碩士研究生，主要研究方向?yàn)樽塑壙貏?dòng)力系統(tǒng)。

通訊作者：羅衛(wèi)東（1962－），男，貴州習(xí)水人，教授，碩士生導(dǎo)師，主要研究方向?yàn)樽塑壙貏?dòng)力系統(tǒng)、汽車及其零部件設(shè)計(jì)理論與方法。