基于Kriging代理模型的液壓閥塊管道結(jié)構(gòu)的優(yōu)化設(shè)計

蘇猛猛， 張紅哲， 梁延德， 何福本

(大連理工大學(xué),遼寧大連 116024)

0 引言

液壓閥塊在液壓系統(tǒng)中應(yīng)用范圍十分廣泛，不僅可以簡化液壓系統(tǒng)的設(shè)計和安裝，而且便于實現(xiàn)液壓系統(tǒng)的集成化和標(biāo)準(zhǔn)化。但液壓閥塊在使用過程中存在一定的壓力損失，導(dǎo)致系統(tǒng)的能量利用率降低。針對該問題，國內(nèi)外學(xué)者做了大量研究。Drago Matko[1]介紹了不同的管道模型及其仿真，并根據(jù)實際管道上的測量結(jié)果對導(dǎo)出的模型進(jìn)行了比較驗證；Baghalian[2]通過使用人工智能方法，在廣泛查閱文獻(xiàn)綜述的基礎(chǔ)上，對90°彎道中的速度場進(jìn)行了分析、實驗和數(shù)值方法的研究；田樹軍[3]運用現(xiàn)代智能計算方法及其混合優(yōu)化策略，實現(xiàn)了滿足性能約束的液壓集成塊布局布孔集成方案的自動優(yōu)化設(shè)計；王永安[4]針對液壓集成塊內(nèi)部孔道中液流的復(fù)雜性流動，探討了造成液流壓力損失的原因；趙宏林[5]分析了不同出油口直徑等不同工藝參數(shù)對流道壓力產(chǎn)生的影響，得到了使流道壓力損失最小的流道尺寸配比；胡建軍[6]分析了交連形式等參數(shù)對集成塊內(nèi)流道壓力損失的影響，探討了圓弧過渡直角轉(zhuǎn)彎流道的壓力損失特性。以上國內(nèi)外學(xué)者的研究均是針對傳統(tǒng)加工方法，在提升液壓閥塊能量利用率的研究過程中，仍會受到傳統(tǒng)加工工藝方法的限制。

本文結(jié)合當(dāng)下增材制造的快速發(fā)展趨勢[7-9]，對液壓閥塊內(nèi)部的正交管道結(jié)構(gòu)進(jìn)行了參數(shù)化設(shè)計，并使用Kriging代理模型優(yōu)化方法對設(shè)計模型進(jìn)行了優(yōu)化，并結(jié)合COMSOL軟件進(jìn)行了流體仿真試驗，得到了壓降值最小的結(jié)構(gòu)模型。

1 參數(shù)化設(shè)計

在液壓閥塊內(nèi)部，正交管道是常見結(jié)構(gòu)形式[10-13]。流體在正交管道轉(zhuǎn)彎處，靠近內(nèi)側(cè)壁的位置，會產(chǎn)生一定大小的渦流，使液壓閥塊產(chǎn)生壓力損失，導(dǎo)致液壓閥塊的能量利用率降低[14-15]。針對該問題，本文對液壓閥塊內(nèi)部常見的正交管道結(jié)構(gòu)進(jìn)行了參數(shù)化設(shè)計，圖1是正交管道結(jié)構(gòu)，圖2是對正交管道結(jié)構(gòu)進(jìn)行參數(shù)化設(shè)計后的模型。

參數(shù)化設(shè)計管道的直徑r1和r2均與原管道直徑相同。彎管與兩端連接處相切，切點b坐標(biāo)為（x1,y1），切點c坐標(biāo)為（x2,y2）。點o為彎管中心線對應(yīng)的曲率中心，曲率半徑為R。點a坐標(biāo)為（x0,y0），a點到b點的距離為l2，a點到o點的距離為l1。由式（1）可用x0、y0、R、x2、y2表示出坐標(biāo)點（x1,y1）的坐標(biāo)值。公式如下所示:

圖1 正交管道

圖2 參數(shù)化設(shè)計管道

該模型以曲率半徑R和c點坐標(biāo)(x2,y2)為自變量，變量范圍 R∈［6,x2］，x2∈［6,50］。

2 優(yōu)化設(shè)計

2.1 代理模型

Kriging代理模型是指在不降低計算精度情況下構(gòu)造的一個計算量小、計算周期短，但計算結(jié)果與數(shù)值分析或物理實驗結(jié)果相近的數(shù)學(xué)模型。應(yīng)用于優(yōu)化領(lǐng)域時，適用于結(jié)構(gòu)功能函數(shù)非線性程度高甚至沒有具體表達(dá)式的情況，對于參數(shù)化模型快速通用。本文使用Kriging代理模型優(yōu)化方法，以參數(shù)化設(shè)計管道入口處平均壓力值p(x)為目標(biāo)函數(shù)，設(shè)計變量為彎管中心線的曲率半徑R和c點坐標(biāo)值，其優(yōu)化列式為

式中：p(x)是流體仿真過程管道入口處的平均壓力；x1和x2為設(shè)計變量。

2.2 流體仿真

運用COMSOL軟件對正交管道和優(yōu)化后參數(shù)化設(shè)計管道結(jié)構(gòu)進(jìn)行流體仿真。對于建好的模型進(jìn)行網(wǎng)格劃分，其中網(wǎng)格的細(xì)分方法選用Delaunay劃分方法，優(yōu)化級別選擇高水平，使其生成的網(wǎng)格更符合該幾何結(jié)構(gòu)，對于生成的網(wǎng)格單元，其尺寸和形狀也盡可能達(dá)到最優(yōu)狀態(tài)。物理場類型選擇標(biāo)準(zhǔn)湍流模型。流體仿真過程采用的相關(guān)邊界條件如下：管道內(nèi)選用的流體為水，溫度為25℃，動力黏度為0.8949×10-3Pa·s，密度為997.044 kg/m3。入口處邊界條件為質(zhì)量流，質(zhì)量流率為0.45 kg/s，湍流強度選擇0.02，湍流長度為0.09 m，出口處壓力設(shè)置為14 MPa。在管道入口處添加邊界探針，用來檢測入口處的平均壓力值。

由仿真結(jié)果可得，正交管道模型入口處的平均壓力值為14.166 MPa，進(jìn)出口壓力差為0.166 MPa。優(yōu)化后參數(shù)化設(shè)計管道入口處的平均壓力值為14.027 MPa，進(jìn)出口壓力差為0.027 MPa，較正交管道模型壓降減少了83.73%，如表1所示。

表1 優(yōu)化結(jié)果

圖3和圖4為正交管道的壓力等值線圖和速度云圖，圖5和圖6為優(yōu)化后參數(shù)化設(shè)計管道的壓力等值線圖與速度云圖。

圖3 正交管道壓力等值線圖

圖4 正交管道速度云圖

圖5 優(yōu)化后設(shè)計管道壓力等值線圖

圖6 優(yōu)化后設(shè)計管道速度云圖

2.3 結(jié)果分析

由伯努利方程[16]式（3）可知，管道內(nèi)流速大時，壓力值則小。

式中：p為壓強；ρ為密度；v為速度。

對比分析圖3和圖5可以發(fā)現(xiàn)，優(yōu)化后參數(shù)化設(shè)計管道內(nèi)部流體在轉(zhuǎn)彎處對管道外側(cè)壁產(chǎn)生的壓力較正交管道有所減小，在轉(zhuǎn)彎處的流速將增大，從圖4和圖6可以得到，優(yōu)化后管道在轉(zhuǎn)彎處流速確實較原管道有所增大。流體經(jīng)過轉(zhuǎn)彎后，在正交管道內(nèi)側(cè)壁存在渦流，并向外側(cè)擴展，導(dǎo)致主流截面減小，流速又會快速增加。

在最大和最小壓力值點的分布上，對比圖3和圖5，可以發(fā)現(xiàn)正交管道的最大和最小壓力值點均分布在轉(zhuǎn)彎處，形成了較大的壓力差，主要是因為流體在轉(zhuǎn)彎處產(chǎn)生的離心力對外壁造成了一定擠壓，內(nèi)側(cè)壁產(chǎn)生的渦流較大，因渦流不參與主流液體流動，會加速液體之間摩擦，造成能量損失。優(yōu)化后參數(shù)化設(shè)計管道結(jié)構(gòu)減小了管道內(nèi)側(cè)壁附近渦流強度，流速也更加平穩(wěn)，從而減小了轉(zhuǎn)彎處的壓力損失，進(jìn)而達(dá)到提升液壓閥塊能量利用率的目的。

3 結(jié)論

本文對液壓閥塊中常見正交管道結(jié)構(gòu)進(jìn)行了參數(shù)化設(shè)計，運用Kriging代理模型優(yōu)化方法對參數(shù)化設(shè)計管道結(jié)構(gòu)進(jìn)行了優(yōu)化，并結(jié)合COMSOL進(jìn)行了流體仿真計算。結(jié)果顯示，優(yōu)化后參數(shù)化設(shè)計管道出入口壓降值較正交管道出入口壓降值降低了83.73%。結(jié)合目前增材制造快速成型的方法，該優(yōu)化設(shè)計有效降低了液壓閥塊的壓力損失，提升了能量利用率，為液壓閥塊整體管道結(jié)構(gòu)的優(yōu)化設(shè)計提供了參考。