茍遼, 楊寧都
(四川凌峰航空液壓機械有限公司,四川 廣漢618300)
錐式單向閥在液壓系統(tǒng)中的應(yīng)用很廣泛,其功能是使液流只能沿一個方向流通而不得向反方向流動。在某些場合錐式單向閥也可起到直動式溢流閥的作用[1]。如圖1所示,錐式單向閥由閥座、閥芯、后蓋和彈簧組成。通常,閥座和后蓋為固定件,閥芯安裝在閥座腔內(nèi),彈簧安裝在閥芯腔內(nèi)。閥座用于安裝閥芯,閥芯在彈簧力的作用下與閥芯組成硬密封,保證壓力油只能向一個方向流動。后蓋和閥芯給彈簧提供一個初始預(yù)壓力,保證彈簧的初始彈力,閥芯靠彈簧的彈力與閥座緊密貼合。
無壓力油作用時,閥芯在彈簧彈力的作用下處于關(guān)閉狀態(tài)。有壓力油作用時,錐式單向閥閥芯在壓力油作用下開啟,壓力油通過閥芯與閥座之間的縫隙從閥芯前端流動至彈簧腔。此時,閥芯與閥座在油液的完全潤滑下工作,二者之間的摩擦力很小。因此,閥芯、閥座和彈簧組成低阻尼振蕩系統(tǒng),根據(jù)其特性,單向閥工作時其部位的壓力容易產(chǎn)生較大的超調(diào)量,使液壓系統(tǒng)的局部壓力瞬時增大、瞬時減小,這樣的壓力突變對液壓系統(tǒng)的穩(wěn)定性造成極為不利的影響并產(chǎn)生噪聲,同時可能造成液壓系統(tǒng)強度較弱的部位破裂。所以,有必要研究單向閥的動態(tài)特性,本文從錐式單向閥低阻尼系統(tǒng)的階躍響應(yīng)方面研究其動態(tài)特性。
圖1 錐式單向閥結(jié)構(gòu)示意圖
目前,國內(nèi)針對單向閥類產(chǎn)品開展的相關(guān)研究有:劉上等[2]進行了單向閥流路系統(tǒng)自激振蕩特性研究,建立了描述該系統(tǒng)動態(tài)過程的非線性動力學(xué)模型,通過數(shù)值計算證明了系統(tǒng)在一定條件下出現(xiàn)極限環(huán),系統(tǒng)是局部穩(wěn)定,大范圍不穩(wěn)定;余武江等[3]對單向閥自激振動機理進行了研究,采用線性分析法提出了單向閥的臨界穩(wěn)定曲線,通過單向閥三維動態(tài)流場仿真分析,利用數(shù)值研究了不同工作壓力、不同質(zhì)量流量條件下的單向閥動態(tài)穩(wěn)定特性;張生昌等[4]利用VOF法建立了混輸泵出口單向閥內(nèi)流場的氣液兩相流CFD模型,采用UDF及動網(wǎng)格技術(shù),研究了介質(zhì)不同含氣率下彈簧剛度對出口單向閥開啟特性的影響;王細波等[5]針對單向閥打開過程中持續(xù)震顫現(xiàn)象,采用AMESim分析軟件建立單向閥數(shù)值計算模型,對不同流量條件和氣體介質(zhì)條件下的工作特性進行仿真計算,得出入口壓力在較低范圍內(nèi)時閥門發(fā)生振蕩,在入口壓力達到一定值后,閥門處于穩(wěn)定狀態(tài)。
本文在此背景下,根據(jù)錐式單向閥的動力學(xué)運動方程、流體流動的連續(xù)性方程和流體流過節(jié)流窗口的流量方程,建立其低阻尼振蕩系統(tǒng)的運動微分方程,通過拉普拉斯變換建立其傳遞函數(shù),分別在原模型的基礎(chǔ)上增大系統(tǒng)阻尼系數(shù)值和彈簧剛度值,在MATLAB/Simulink環(huán)境下對其數(shù)學(xué)模型進行階躍響應(yīng)計算仿真,研究在阻尼系數(shù)和彈簧剛度變化的情況下系統(tǒng)階躍響應(yīng)品質(zhì)的各評判參數(shù),為工程實際應(yīng)用提供一定的參考意見。
式中:ΔP=ps為進油壓力;A為閥芯承受液壓力的有效面積;ks為液動力系數(shù),ks=CdCvπdmsin 2β;x0為閥口開度為零時的彈簧預(yù)壓縮量;x為閥芯位移;m為運動件的質(zhì)量,其大小為閥芯質(zhì)量及1/3彈簧質(zhì)量之和;c為閥芯運動黏性阻尼系數(shù);k為彈簧剛度。
當(dāng)壓力油進入錐式單向閥時,閥芯在壓力油作用下相對閥座開啟形成通油窗口,流體通過通油窗口流動時,閥芯產(chǎn)生運動。運動過程中的閥芯承受6種力,分別為由壓力油產(chǎn)生的液壓力、自身慣性力、與閥芯運動速度成正比的阻尼力、與彈簧壓縮量成正比的彈簧力、閥芯與閥座間的摩擦力、閥芯運動時產(chǎn)生的液動力[6]。
此外,閥芯在開啟過程中會產(chǎn)生兩種液動力,分別是瞬態(tài)液動力和穩(wěn)態(tài)液動力。當(dāng)閥芯在開啟過程中,閥芯運動導(dǎo)致主閥口開度x不斷變化,閥腔內(nèi)液體的流速也不斷變化,液體的慣性會在閥芯上產(chǎn)生作用力,即瞬態(tài)液動力。
閥芯開度的變化速率是影響瞬態(tài)液動力的因素,而閥芯開度不影響瞬態(tài)液動力大小。一般而言,瞬態(tài)液動力對閥芯的作用力非常小,因此在其運動方程中對這部分力忽略不計[7]。
為計算方便,在ps=ps0,x=x1領(lǐng)域內(nèi)對式(2)進行小偏差線性化處理,得[3]
當(dāng)壓力油在單向閥中流動時,其液體流量連續(xù)性方程以壓力流量連續(xù)性方程為基礎(chǔ)。對于閥芯而言,其液體從進油口進入,通過閥口開度流出。由閥芯開度處流出的流量遵循液體壓力流量連續(xù)性方程。
閥芯運動產(chǎn)生位移后會導(dǎo)致閥芯下腔容積變大,變大的容積需要液壓液填充,由此導(dǎo)致的影響體現(xiàn)在流量連續(xù)方程中為QV=A1dy/dt。
對式(6)進行拉氏變換得
理論上,當(dāng)錐式單向閥的閥芯在壓力油的作用下相對閥座打開一定的開度時,閥芯頭部與閥座之間形成一個同心環(huán)型間隙。但是在實際的工作過程中,這個同心環(huán)型間隙在外界因素(如振動等)的影響下,不一定是同心環(huán)型間隙。因此,錐式單向閥開啟后其流量方程選擇薄壁孔流量公式:
合并式(4)、式(7)和式(10),得
根據(jù)式(11),建立錐式單向閥低阻尼振蕩系統(tǒng)的方塊圖,如圖2所示。
將圖2所示的錐式單向閥低阻尼振蕩系統(tǒng)方塊圖演化為由典型環(huán)節(jié)組成的方塊圖,其方塊圖如圖3所示。
根據(jù)圖3所示的方框圖,可得錐式單向閥低阻尼振蕩系統(tǒng)的傳遞函數(shù)為
圖2 方塊圖
圖3 轉(zhuǎn)化為典型環(huán)節(jié)的方塊圖
對該傳遞函數(shù)表達式的分母進行簡化,得
分析思路:為分析錐式單向閥低阻尼振蕩系統(tǒng)中阻尼系數(shù)和彈簧剛度對系統(tǒng)階躍響應(yīng)的影響,在其原模型的基礎(chǔ)上,分別增大系統(tǒng)阻尼系數(shù)和彈簧剛度,在MATLAB/Simulink環(huán)境下仿真計算系統(tǒng)的階躍響應(yīng)。
在MATLAB/Simulink環(huán)境下,根據(jù)傳遞函數(shù)建立錐式單向閥低阻尼振蕩系統(tǒng)的框圖,給錐式單向閥低阻尼振蕩系統(tǒng)的原始模型輸入階躍信號,進行仿真計算[8-9],得到圖4所示的仿真曲線。
如圖4所示,錐式單向閥低阻尼振蕩系統(tǒng)在階躍信號的激勵下,其上升時間約為0.35 s,最大峰值約為1.28 MPa,調(diào)整時間約為1.7 s,穩(wěn)定時間為3 s。
增大錐式單向閥低阻尼振蕩系統(tǒng)的阻尼系數(shù),其他參數(shù)不變,進行仿真計算,得到如圖5所示的仿真計算曲線,其上升時間約為0.42 s,最大峰值約為1.15 MPa,調(diào)整時間約為1.6 s,穩(wěn)定時間為2.5 s。
增大錐式單向閥低阻尼振蕩系統(tǒng)的彈簧剛度,其他參數(shù)不變,進行仿真計算,得到如圖6所示的仿真計算曲線,其上升時間約為0.51 s,最大峰值約為1.32 MPa,調(diào)整時間約為3.2 s,穩(wěn)定時間為4.5 s。
圖4 仿真曲線1
由表1可知,在保證其他參數(shù)不變的情況下,增大系統(tǒng)的阻尼系數(shù),上升時間會增大,但是系統(tǒng)的調(diào)整時間和穩(wěn)定時間都相對降低,且系統(tǒng)的壓力峰值也相對降低。
在保證其他參數(shù)不變的情況下,增大系統(tǒng)的彈簧剛度,上升時間、調(diào)整時間和穩(wěn)定時間都相對增加,壓力峰值也相對升高。
圖5 仿真曲線2
圖6 仿真曲線3
表1 仿真參數(shù)
綜上分析,在工程應(yīng)用中可根據(jù)實際需求,增大或減小錐式單向閥低阻尼振蕩系統(tǒng)的阻尼系數(shù)或彈簧剛度來達到應(yīng)用目的。
錐式單向閥工作時,閥芯、閥座和彈簧組成一個低阻尼振蕩系統(tǒng),本文在忽略溫度對壓力油的影響、主閥的卡緊力、流量的脈動等因素的條件下,建立了錐式單向閥低阻尼振蕩系統(tǒng)的數(shù)學(xué)模型,并在MATLAB/Simulink環(huán)境下研究了阻尼系數(shù)和彈簧剛度對其階躍響應(yīng)的影響,結(jié)論如下:
1)在原錐式單向閥結(jié)構(gòu)的基礎(chǔ)上,增大低阻尼振蕩系統(tǒng)的阻尼系數(shù),系統(tǒng)階躍響應(yīng)的峰值由原系統(tǒng)的1.28 MPa降低至1.15 MPa,說明增大阻尼系數(shù)有助于降低系統(tǒng)峰值。2)在原錐式單向閥結(jié)構(gòu)的基礎(chǔ)上,增大低阻尼振蕩系統(tǒng)的彈簧剛度,系統(tǒng)階躍響應(yīng)的超調(diào)量由原系統(tǒng)的1.28 MPa增大至1.32 MPa,說明增大彈簧剛度,有助于增大系統(tǒng)峰值。3)系統(tǒng)阻尼系數(shù)和彈簧剛度的變化,會影響系統(tǒng)階躍響應(yīng)的各時間參數(shù),在實際工程應(yīng)用中可根據(jù)實際需求調(diào)整阻尼系數(shù)和彈簧剛度。