張順鋒，張 勝，王 強，何曉暉，徐 磊

(解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院， 南京 210007)

【機械制造與檢測技術(shù)】

流固耦合作用下的液控換向閥動態(tài)特性仿真

張順鋒，張 勝，王 強，何曉暉，徐 磊

(解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院， 南京 210007)

基于流固耦合理論和動網(wǎng)格技術(shù)，運用CFD-ACE+軟件建立了計算流體力學(xué)三維模型。對在彈簧力、驅(qū)動力及流體力的作用下液控換向閥的開啟過程進行了仿真，仿真結(jié)果與理論值、AMESim數(shù)值計算基本一致。對3種節(jié)流結(jié)構(gòu)的液控換向閥開啟過程瞬時的流場、閥芯位移、液動力、von Mises應(yīng)力分布和形變進行了對比。

液控換向閥；流固耦合；動網(wǎng)格；動態(tài)仿真

simulation

在各類液壓系統(tǒng)中液控換向閥應(yīng)用廣泛，其性能影響著整個液壓系統(tǒng)的工作?；y式換向閥的主要性能是平穩(wěn)換向，迅速反應(yīng)?；y內(nèi)部流場、液動力以及動態(tài)響應(yīng)特性是目前研究的熱點。趙蕾等[1-3]利用計算流體動力學(xué)方法對液壓滑閥開啟過程或關(guān)閉時內(nèi)部的流場分布及液動力進行了仿真分析，其中動態(tài)模擬時閥芯速度是預(yù)先給定，實際上滑閥開啟中閥芯速度是不斷變化著的。姜濤[4]對換向閥P-B通道(進行了動態(tài)仿真研究，換向閥的動態(tài)開啟過程是通過假定閥芯不動來實現(xiàn)的，但在實際運行中閥芯和流體相互影響，需要考慮流固耦合的影響。張蕊華[5]對閥芯及其周圍流場數(shù)值模擬時考慮了流固耦合，先對流場分析進行了計算，然后將數(shù)值結(jié)果加載到固體結(jié)構(gòu)上，顯然這種流固耦合的數(shù)值計算是單向的。在換向閥的實際工作中，換向閥的開合是一個動態(tài)過程，閥芯的運動及形變相反又會對流體運動產(chǎn)生重要影響[6]。

針對上述問題，將流固耦合理論和動網(wǎng)格技術(shù)相結(jié)合，以流體力學(xué)軟件CFD-ACE+為平臺，對3種不同節(jié)流結(jié)構(gòu)的液控換向閥在彈簧力、驅(qū)動力及流體力作用下的開啟過程進行動態(tài)模擬分析，為改善液壓滑閥內(nèi)部流體流動狀態(tài)、動態(tài)響應(yīng)特性和閥芯受力狀況提供依據(jù)。

1 三維幾何模型和網(wǎng)格劃分

液控換向滑閥由多個閥體內(nèi)腔組成，本文主要對其中的一個閥體內(nèi)腔進行研究[7]。以P-B口開啟過程為例，其結(jié)構(gòu)示意圖如圖1所示，在出口節(jié)流處閥芯臺肩上設(shè)有3種不同的節(jié)流結(jié)構(gòu)。

借鑒朱鈺[8]的簡化方法，將液控換向滑閥閥芯兩端部分及其附近流域略去，彈簧模型加載到閥芯一端。由于閥體內(nèi)腔、流域和閥的內(nèi)芯均以閥中央對稱布置，為提高計算效率取流體區(qū)域和閥芯的對稱作為研究對象。

在劃分網(wǎng)格時采用高精度六面體結(jié)構(gòu)網(wǎng)格，考慮流固耦合導(dǎo)致網(wǎng)格形變較大，三維瞬態(tài)流場解析誤差太大，因此采用單一的網(wǎng)格變形方式，劃分網(wǎng)格的變形能力和變形品質(zhì)就很難控制。本研究將彈簧法和超限插值法聯(lián)系起來處理結(jié)構(gòu)網(wǎng)格的形變。基于分塊網(wǎng)格的觀點，將計算區(qū)域分成16個網(wǎng)格區(qū)。采取彈簧聯(lián)系不同網(wǎng)格區(qū)的角點，通過彈簧啟動各角點的形變[9]。3種不同節(jié)流結(jié)構(gòu)的液控換向滑閥網(wǎng)格劃分情況如圖2所示，在閥芯出口臺肩處為A、B、C 3種不同的節(jié)流結(jié)構(gòu)。對節(jié)流處進行細化網(wǎng)格可以大大提高數(shù)值模擬結(jié)果的準(zhǔn)確性。

圖1 液控換向滑閥結(jié)構(gòu)示意圖

圖2 網(wǎng)格劃分

2 計算條件及參數(shù)設(shè)置

入口條件：壓力入口，選取6 MPa；

出口條件：壓力出口，設(shè)為大氣壓；

閥芯驅(qū)動壓差：1 MPa；

彈簧剛度：80 N/mm；

液壓油：型號VG46，密度ρ=890 kg/m3；運動粘度μ=46 mm2/s；設(shè)為不可壓縮流體；

湍流模型：k-ε湍流模型。

仿真中不考慮溫度、泄漏的影響，閥芯和閥體之間配合良好，無徑向間隙，閥芯只考慮軸向位移。

3 仿真結(jié)果分析

3.1 閥芯運動情況分析

閥芯的運動影響因素主要是驅(qū)動力、液動力以及彈簧力。為驗證數(shù)值計算結(jié)果的正確性，利用AMESim軟件建立如圖3所示的液控換向閥模型，與CFD-ACE+所建立的無節(jié)流結(jié)構(gòu)理想滑閥模型仿真結(jié)果進行對比，兩種軟件仿真得到的閥芯位移、速度和加速度對比情況如圖4所示。

圖3 液控換向閥AMESim模型

閥芯穩(wěn)態(tài)位移的計算公式為：

(3)

式中：D為流道直徑；xv為閥芯穩(wěn)態(tài)位移；W為閥口的面積梯度；k為彈簧剛度；F為滑閥閥芯所受的驅(qū)動力；Fs為閥芯所受的穩(wěn)態(tài)液動力；Δp為滑閥進出口壓差；Δp1為驅(qū)動壓差。

從式(3)，解得閥芯的穩(wěn)態(tài)位移為0.71 mm，而CFD-ACE+和AMESim仿真的結(jié)果與理論計算值誤差分別為3.1%和11.3%，CFD-ACE+的數(shù)值計算結(jié)果因考慮了流固耦合作用和流場微小流動帶來的影響而更加精準(zhǔn)。但因仿真條件的設(shè)置和理論計算的方法不可能完全一致，仍存在一定的誤差。

圖4 CFD-ACE+和AMESim仿真結(jié)果

利用CFD-ACE+仿真得到的3種節(jié)流結(jié)構(gòu)的液控換向閥閥芯位移隨時間的變化情況如圖5所示，假設(shè)當(dāng)閥芯位移同一波動周期內(nèi)波動幅度與該段波峰值之比小于0.5%時，即為穩(wěn)態(tài)。

從圖5可看出在液控換向閥P-B開啟過程中，3種節(jié)流結(jié)構(gòu)閥芯位移均存在一定的波動，超調(diào)量分別為45.9%、48.9%和48.7%，穩(wěn)態(tài)位移分別為0.64 mm、0.59 mm和0.48 mm。A型超調(diào)量最小，穩(wěn)態(tài)位移最大，達到穩(wěn)態(tài)所需的時間也最短。

圖5 閥芯開啟過程的位移

3.2 液動力分析

滑閥閥芯表面所受到的軸向力可以由閥芯表面所受的壓力對表面積積分得到，閥芯表面軸向力即為液動力[10-12]。圖6給出了3種閥芯開啟過程中所受液動力變化情況，包括瞬態(tài)液動力和穩(wěn)態(tài)液動力，其中瞬態(tài)液動力等于不同瞬間總的液動力與穩(wěn)態(tài)液動力之差?？芍€(wěn)態(tài)液動力的大小均為負值，方向即為閥口關(guān)閉的方向。在閥芯開啟過程瞬間，瞬態(tài)液動力變化較大，設(shè)計閥時須加以考慮。A型閥芯穩(wěn)態(tài)液動力值最小，達到穩(wěn)態(tài)所需的時間最短，但瞬態(tài)液動力波動峰值較大。

圖6 閥芯開啟過程液動力變化

3.3 流固耦合分析

圖7給出了3種節(jié)流結(jié)構(gòu)液控換向閥在開啟初期、最大位移和穩(wěn)態(tài)時流域的壓力場和閥芯所受的馮.米斯應(yīng)力場。馮·米斯應(yīng)力是根據(jù)第四強度理論得到的一種當(dāng)量應(yīng)力，可以用來對疲勞、破壞等進行評價[13]。

圖7 不同時刻壓力場和應(yīng)力場

從圖7可以看出，在閥口節(jié)流區(qū)域靠近閥體的流場均會產(chǎn)生負壓，其中B型閥芯負壓區(qū)域較大。流場的負壓導(dǎo)致產(chǎn)生氣蝕。由于氣蝕，閥芯則容易破壞。在入口靠近閥桿的流場均出現(xiàn)了高壓區(qū)，這會對閥芯形成高壓沖擊，從而影響閥芯的徑向平衡，容易造成閥芯卡死。3種節(jié)流結(jié)構(gòu)閥芯所受馮.米斯應(yīng)力值的分布基本一致，在閥芯臺肩和閥桿連接處馮.米斯應(yīng)力值均最大。圖8為閥芯開啟過程應(yīng)力值的峰值隨時間變化情況。從其中可知，閥芯開啟過程最大馮.米塞斯應(yīng)力值均存在波動，其中A型閥芯穩(wěn)態(tài)值最小，B型閥體內(nèi)芯穩(wěn)態(tài)值和瞬態(tài)值最大。B型閥芯開啟過程所受馮.米塞斯應(yīng)力的峰值為35.5 MPa，閥芯材料(20CrMo)的抗拉、抗剪和抗壓強度都大于200 MPa[14]，根據(jù)換向閥的強度理論，換向閥使用過程中的安全系數(shù)取為2，由此文中所研究的3種節(jié)流結(jié)構(gòu)的閥芯強度均滿足要求，設(shè)計時無需考慮閥芯的疲勞破壞。

圖9給出了3種節(jié)流結(jié)構(gòu)閥芯在開啟初期、最大位移和穩(wěn)態(tài)位移時閥芯的形變及其分布情況。

從圖9可看出，閥芯的徑向形變主要發(fā)生在閥桿上，最大形變分布位置和大小隨時間不斷變化。達到穩(wěn)態(tài)時，閥桿在入口處靠近閥口方向和出口處背離閥口方向產(chǎn)生較大變形。穩(wěn)態(tài)時3種節(jié)流結(jié)構(gòu)閥芯的最大形變分布位置基本相同，B型閥芯最大形變較大，A型閥芯最大形變較小。

圖8 閥芯開啟過程最大馮·米塞斯應(yīng)力值變化

圖9 不同時刻閥芯形變

4 結(jié)論

1) 液控換向閥開啟時在閥口節(jié)流區(qū)域靠近閥體的流場產(chǎn)生負壓區(qū)，在入口靠近閥桿的流場出現(xiàn)高壓區(qū)；

2) 在閥芯開啟同時，瞬態(tài)液動力值變動較大，馮·米斯應(yīng)力最大值主要集中在閥芯臺肩和閥桿連接處，設(shè)計閥時須加以考慮；

3) 穩(wěn)態(tài)時3種節(jié)流結(jié)構(gòu)閥芯的最大變形分布基本相同， A型閥芯最大變形較??；

4) A型閥芯在試驗條件下達到穩(wěn)態(tài)時位移最大，液動力、穩(wěn)態(tài)馮.米塞斯應(yīng)力值最小和位移超調(diào)量最小，時間最短，設(shè)計閥芯時可以考慮采用該節(jié)流結(jié)構(gòu)。

[1] 趙蕾,陳青,權(quán)龍.閥芯運動狀態(tài)滑閥內(nèi)部流場的可視化分析[J] 農(nóng)業(yè)機械學(xué)報,2008(11):142-145，155.

[2] 吳小鋒,干為民,劉春節(jié),等.液壓換向滑閥內(nèi)部結(jié)構(gòu)的健壯性設(shè)計[J].中國機械工程,2015,26(15):2030-2035，2040.

[3] 楊科,金曉宏,肖鵬飛,等.非全周矩形開口滑閥小開口度時流量及液動力特性研究[J].流體機械,2016(7):14-17，62.

[4] 姜濤,黃偉,王安麟.多路閥閥芯節(jié)流槽拓撲結(jié)構(gòu)組合的神經(jīng)網(wǎng)絡(luò)模型[J].西安交通大學(xué)學(xué)報,2016(6):36-41.

[5] 張蕊華,吳啟帆,易健,等.壓電驅(qū)動精密流量閥閥芯結(jié)構(gòu)優(yōu)化設(shè)計[J].潤滑與密封,2016(5):97-101，107.

[6] 姜濤,夏明杰,王安麟.滑移網(wǎng)格法在滑閥三維瞬態(tài)流場解析的適用性分析[J].同濟大學(xué)學(xué)報(自然科學(xué)版),2015(10):1575-1581.

[7] 劉罡,張俊俊,王江勇.基于ADINA液壓滑閥的流固耦合研究[J].機床與液壓,2012(23):158-160.

[8] 朱鈺.液控換向閥內(nèi)流場及動態(tài)特性的數(shù)值模擬[J].哈爾濱工業(yè)大學(xué)學(xué)報,2012(5):133-139.

[9] SINGHAL A K,ATHAVALE M M.Mathematical basis and validation of the full cavitation model[J].Journal of Fluids Engineering,2002,124:617-624.

[10] 何曉暉,孫宏才,程健生,等.基于動網(wǎng)格的液壓閥閥芯啟閉中的液動力分析[J].解放軍理工大學(xué)學(xué)報(自然科學(xué)版),2011(5):491-495.

[11] 閆清東,劉博深,魏巍.基于動網(wǎng)格的沖焊型液力變矩器流固耦合分析[J].華中科技大學(xué)學(xué)報(自然科學(xué)版),2015(12):37-41.

[12] 張曉俊,權(quán)龍,趙斌.內(nèi)流式滑閥壁面壓力分布可視化計算及試驗驗證[J].機械工程學(xué)報,2016(14):196-203.

[13] 何曉暉,欒健,王強.基于CFD的液壓錐閥開啟過程流固耦合分析[J].液壓與氣動,2015(9):122-125.

[14] 李建心,劉效東.液動力對換向閥作用的ANSYS有限元分析[J].機床與液壓,2010,38(17):119-120.

(責(zé)任編輯周江川)

TransientPropertiesofHydraulicOperatedDirectionalValveBasedonFluid-structureInteraction

ZHANG Shunfeng, ZHANG Sheng, WANG Qiang, HE Xiaohui, XU Lei

(College of Field Engineering, PLA University of Science & Technology, Nanjing 210007, China)

Three-dimensional CFD modals were established via CFD-ACE+ software based on fluid-structure interaction and dynamic mesh.The opening process of hydraulic operated directional valves under the spring force, driving force and flow force process were simulated and the results are consistent with the theoretical values and values calculated by AMESim.The transient flow fields, spool displacement, flow force, distribution of von Mises stress, deformation of three different throttle structure of hydraulic operated directional valves in the opening process were compared and analyzed.

hydraulic operated directional valve; fluid-structure interaction; dynamic mesh; dynamic

2017-04-23；

：2017-05-11

：國家自然科學(xué)基金項目(51175511)

張順鋒,碩士研究生,主要從事液壓傳動與控制研究。

王強，教授，主要從事液壓傳動與控制研究，E-mail：wangqiangjs@sohu.com。

10.11809/scbgxb2017.09.032

format:ZHANG Shunfeng, ZHANG Sheng, WANG Qiang, et al.Transient Properties of Hydraulic Operated Directional Valve Based on Fluid-structure Interaction[J].Journal of Ordnance Equipment Engineering,2017(9):151-155.

TH137

：A

2096-2304(2017)09-0151-05

本文引用格式：張順鋒，張勝，王強，等.流固耦合作用下的液控換向閥動態(tài)特性仿真[J].兵器裝備工程學(xué)報，2017(9):151-155.