胡均平，胡 騫

（中南大學機電工程學院，長沙 410083）

液壓沖擊器零位負開口配油法油壓脈沖特性研究

胡均平，胡 騫

（中南大學機電工程學院，長沙 410083）

針對國內液壓沖擊器普遍采用零位正開口配油閥而存在的高、低壓油互混的問題，探究采用負開口配油時沖擊器的動態(tài)特性，以尋求配油閥最佳的零位負開口型式。依據(jù)液壓沖擊機構的基本原理，結合流體縫隙流動時流態(tài)轉變機理，建立動態(tài)非線性數(shù)學模型，并采用Stateflow模擬頻繁換向過程，搭建了液壓沖擊器Simulink仿真模型，對采用零位負開口配油方式的液壓沖擊器工作壓力和流量特性等進行了研究，并設計實驗進行了實測分析，實驗所得結論與仿真結果相符。在此過程中又著重分析了不同的零位負開口量對沖擊器活塞腔室壓力的影響，為液壓沖擊器零位負開口配油方式的設計提供了參考依據(jù)。

液壓沖擊器；配油閥；零位負開口量；壓力脈沖

液壓沖擊器是破碎鑿巖機械的重要執(zhí)行機構，現(xiàn)今國內對其核心部件配油閥的研究，大多是針對采用零位正開口滑閥作為配油閥的結構型式進行的［1－2］，而對于采用零位負開口閥的情況則研究甚少，且這些研究大都認為負開口閥配油會導致活塞工作腔產生瞬間封閉，形成有害的油壓沖擊［3］。然而在工程上，采用負開口閥作為配油閥的案例卻早已出現(xiàn)，例如瑞典Atlas－Copco公司生產的Cop1238液壓沖擊器［4－5］，該沖擊器能避免正開口配油帶來的類似高、低壓油互混的問題［6－7］，使液壓油的壓力能利用率得到了進一步提高。

基于這一現(xiàn)狀，本文擬探究采用負開口閥配油時液壓沖擊器的動態(tài)特性，通過建立Matlab Simulink仿真平臺，對液壓沖擊器采用負開口配油閥時的性能參數(shù)進行仿真研究，重點闡述不同的零位負開口量對各腔室壓力、輸入壓力油流量和進油口壓力等的影響。

1 工作原理及負開口閥配油特點

1.1 配油原理與結構參數(shù)

本文采用雙控式的配油方案，如圖1所示，并以Atlas Cop1238液壓沖擊器作為仿真結構參照，表1列出了其主要結構參數(shù)。

圖1 雙控式配油原理圖Fig.1 Schematic of dual control oil distributing system for hydraulic impactor

表1 Atlas Cop1238沖擊器主要結構參數(shù)Tab.1 Parameters of Atlas Cop1238 impactor

圖1所示狀態(tài)為活塞已完成上一次沖擊，閥結束沖程換向，整個系統(tǒng)處于回程開始階段。此時壓力油與前腔相通，后腔通回油，在前后腔壓力差的作用下，活塞將向左作回程運動?；钊爻探Y束后，閥也已回程換向，使壓力油通往后腔，前腔則通回油，這樣活塞必向右沖程加速，直到打擊釬桿完成一次沖擊［8］。

1.2 負開口閥配油特點

負開口閥也叫正重疊閥（Overlap Valve），就是當閥芯處于零位的時候閥芯凸肩與閥座凸肩的重疊量為正。

圖2展現(xiàn)出負開口閥閥芯在負開口區(qū)間附近運動時，閥芯凸肩與閥座凸肩的重疊關系及各腔室連通情況，閥芯箭頭指向為閥芯運動方向，油腔內箭頭為油液流動方向。由（b）到（c）的過程，閥芯正好處于負開口區(qū)間。

當閥芯正好處于負開口區(qū)間的時候，原本被閥道溝通的腔室此刻完全被隔開，成為各自獨立的腔室，壓力也被隔斷，高、低壓油不會混到一起。

圖2 閥芯運動于負開口區(qū)間結構示意圖Fig.2 Structure schematic of overlapped valve

2 負開口配油法數(shù)學模型

2.1 閥芯運動于負開口區(qū)間油液特征

沖擊器配油閥換向頻率高達50～60 Hz，且一般要求在1～2ms的時間內完成大開口量的油路切換動作，因而流經配油閥的油液常處于節(jié)流狀態(tài)［9］。閥芯運動于負開口區(qū)間的時間雖然很短，但油液節(jié)流情況卻最為復雜，有必要對閥芯運動于負開口區(qū)間及其附近的過程作進一步細分，如表2所示。

正常節(jié)流狀態(tài)時，閥芯還未進入負開口區(qū)間。表中參數(shù)作為區(qū)別正常節(jié)流和強烈節(jié)流狀態(tài)的指標，與雷諾系數(shù)Re直接相關［10］：

式中：Q表示流經閥道的油液流量；dv為閥凸肩的直徑；v為油液流速；μ為油液運動粘度（取為18×10－6m2/s）；取Re＝100。

表2 閥芯運動于負開口區(qū)間及其附近的細分方法Tab.2 Subdivision for valve spool movement in the negative open range

強烈節(jié)流狀態(tài)時，閥芯開口量處于Amin和Ak之間，Amin的計算方式與Ak類似，只是雷諾系數(shù)取為Re＝1 200。隨著閥開口量At的減小，At對壓力的影響會越來越大。此時需要對油液情況作進一步處理，基于液壓流體力學的薄壁小孔和環(huán)形縫隙分析理論，結合動態(tài)節(jié)流分析方法，作如下合理簡化［11］：

式中：Δpd為節(jié)流狀態(tài)壓降；ρ為油液密度（取900 kg/m3）；cd為流量系數(shù)（這里取為0.83）；A＾t為前一時刻閥開口量；uv為閥運動速度；Δpr為高壓腔壓力上升量；c為沖擊波在液體中的傳播速度（取為1 000 m/s）；Δv為流速下降值；tc為沖擊波傳播時間；t為閥口關閉時間。由于tc＝2l/c，t＝s/uv，l為油道長度（等效為0.37 m），s為閥芯剩余開口長度，不難得到s＝（Ak－Amin）/π/d，Δv＝Q/Ak，所以式（4）可化為：

腔室封閉狀態(tài)時，原本被閥道溝通的腔室成為各自獨立的油腔?；钊那昂笄粚⒊霈F(xiàn)瞬間的封閉，若此時活塞回程減速速度不為零，將導致后腔的油液受到進一步壓縮。這種狀態(tài)持續(xù)時間非常短，不能簡單認為會產生有害的壓力沖擊。

首先必須將油液本身的可壓縮性考慮進來［12］：

式中：KL為液壓油的體積彈性模量（取1.6×109N/m2）；Δpc為壓力上升量；ΔV為油液被壓縮的體積；V0為壓縮前的體積；Δyp為閥芯進入負開口區(qū)間時活塞走過的位移；A為油液作用面積。

另一方面，還應考慮油路的流量補償作用，它也能減緩油壓突變［13］。設油道的體積彈性模量為Ke（各段綜合等效為3.5×108N/m2），則：

式中：Δpe為壓力改變量；ΔVe為補償油液的體積；VE與油道內油液體積相關（取為500 mL）。

不難發(fā)現(xiàn)，式（6）和式（9）都遵循關系：

λ為比例系數(shù)，這就好比在油路中安裝了一個剛度很大的蓄能器。即，上述兩方面的原因可等效為在活塞油腔和配油閥油腔間存在一個可吸收壓力沖擊的蓄能器。

2.2 運動體微分方程

目前對于液壓沖擊器活塞和閥芯運動規(guī)律的研究方法主要有線性模型和非線性模型兩種。前者簡明易懂、計算簡便，但不能用于十分精確的設計性研究；后者考慮了壓力、流量等的非線性變化，更接近實際情況［14］。何清華［8］提出的三段分析法綜合了上述兩種方法的優(yōu)勢，能實現(xiàn)更為精確的計算機仿真。下面本文基于三段分析法來指導仿真。

根據(jù)運動體（包括活塞和閥芯）的動力學平衡、流量平衡及氣體狀態(tài)方程，可得液壓沖擊機構運動體聯(lián)合運動的基本微分方程組：

（1）活塞動力學平衡方程

實際流量Qt＝Qi＋Qe，其中Ph為蓄能器工作壓力，Pah為蓄能器充氣壓力，Vh為蓄能器氣室容積，Vah為蓄能器充氣容積，k為氣體絕熱系數(shù)（氮氣取1.4），kv＝50（KL＋Ke）/（KL·Ke），為補償流量系數(shù)（算得1.8 ×10－7），B1～B9為與結構參數(shù)有關的量，p1，p2，p3對應前述壓力。

3 負開口配油法仿真與結果分析

3.1 Stateflow狀態(tài)轉換

液壓沖擊機構一個周期的耗時約0.02 s，而工作狀態(tài)細分后卻達十幾個，且每個狀態(tài)都有各自不同性能參數(shù)的變化規(guī)律。為此，本文采用Matlab中的Stateflow模塊來生成狀態(tài)轉換監(jiān)控邏輯，解決活塞和配油閥復雜的狀態(tài)循環(huán)轉換問題［15］。圖3就是在Stateflow模塊中建立的液壓沖擊機構狀態(tài)轉換流程圖。

內嵌的initial函數(shù)將閥芯運動于各區(qū)間的特征參數(shù)，以及諸如狀態(tài)指標Ak、油壓改變量Δpd、Δpr、Δpe等中間變量的關系式通過m函數(shù)文件初始化到模型中，以供另一內嵌的函數(shù)Fcn＿pn調用，計算不同階段各腔室的壓力和流量，所得結果作用于活塞和閥芯子模塊。

圖3 基于Stateflow的狀態(tài)流程圖Fig.3 States transfer for control valvemovement

3.2 Simu link仿真模型

為構成一個穩(wěn)定的閉環(huán)系統(tǒng)，還需將Stateflow模塊嵌入到Simulink平臺中與其他模塊連接，共同完成參數(shù)的輸入、輸出，誤差的修正和反饋。

整個模型以沖擊能量Ei、沖擊頻率f、系統(tǒng)工作油壓pi和基本結構參數(shù)（包括閥芯零位開口量）為輸入量，通過對Stateflow模塊中的兩內嵌函數(shù)循環(huán)調用，計算出腔室壓力和流量，輸入到piston和valve子模塊中，又將輸出的活塞速度和閥芯速度反饋回模型輸入端，不斷修正活塞位移和閥芯位移，使系統(tǒng)達到穩(wěn)定。如圖4所示。

圖4 液壓沖擊機構simulink仿真模型Fig.4 Simulationmodel for hydraulic impactor

3.3 仿真結果及油壓脈沖分析

模型采用ODE45算法求解，依常用工況令Ei＝300 J，f＝50 Hz，pi＝15MPa。初選仿真負開口量z0＝0.3 mm，得圖5（a）所示穩(wěn)定后的活塞和閥芯速度曲線。

圖5 負開口量z0＝0.000 3 m時仿真結果曲線圖Fig.5 Simulation curves at z0＝0.000 3 m

圖5（a）中，以回程速度方向為正方向，up為活塞速度曲線，活塞每個運動周期存在明顯的回程加速、回程減速、沖程加速和打擊停頓四個階段；uv為閥芯速度曲線，曲線圖像精確反映出閥芯在回程加速最末時和沖程最末時迅速進行回程和沖程換向，達到理想配油狀態(tài)。

圖5（b）是進油口壓力p和實際流量Qt的變化曲線。p的變化曲線表明，盡管泵輸出穩(wěn)定的壓力pi＝15 MPa，但沖擊器入口處油壓卻并未穩(wěn)定在這個值，而是在13.4 MPa～14.1 MPa這個區(qū)間內呈周期性變化。同樣地，進油口流量Qt也并不是穩(wěn)定于壓力油源輸入的流量Qt＝90 L/min，而是在70～90 L/min的范圍內呈周期性變化。由此可知，在高頻換向的工作狀態(tài)下，沖擊器高壓膠管的補償流量作用對進口油壓的影響不可忽略。

圖5（c）和圖5（d）是系統(tǒng)穩(wěn)定后壓力p1，p2變化情況的仿真曲線，圈出部分為出現(xiàn)油壓脈沖區(qū)間：前腔壓力p2在半周期時出現(xiàn)壓力驟降，之后立即回歸平穩(wěn)；而后腔壓力p1在對應時刻出現(xiàn)壓力脈沖。由圖5（c）可知，系統(tǒng)油壓為15 MPa時后腔壓力脈沖峰值油壓為16.8 MPa，超出系統(tǒng)壓力12%左右，并未出現(xiàn)劇烈沖擊。

因此，即便閥芯在負開口區(qū)間時活塞前后腔室出現(xiàn)了瞬時封閉，但由于高壓膠管的補償流量作用和壓力交替變時油液可壓縮性，使得油壓脈沖大幅減緩。

4 實驗研究與油壓脈沖特性分析

4.1 實驗方案與結果分析

為驗證文中采用零位負開口配油閥所建立模型的合理性，筆者在原有的正開口配油閥沖擊器的基礎上設計制作了相匹配的零位負開口量為0.3 mm的配油閥閥芯。為進一步探究油壓變化的特性，引入壓力傳感器測量沖擊器活塞腔室壓力變化，實測數(shù)據(jù)經DAQ系統(tǒng)采集，最后輸入到計算機中進行波形顯示，如圖6所示，實驗臺的液壓泵輸出穩(wěn)定流量90 L/min。

圖6 腔室壓力測量實驗Fig.6 Experiment for pressuremeasure of chamber

圖7 負開口量z0＝0.000 3 m時實測曲線與仿真曲線圖Fig.7 Comparison of simulation curves and test curves at z0＝0.000 3 m

圖8（a） 負開口量z0＝0.000 1～0.000 5 m時p1曲線圖Fig.8（a）Curve of p1at z0＝0.000 1～0.000 5 m

圖8（b） 負開口量z0＝0.000 1～0.000 5 m時p2曲線圖Fig.8（b）Curve of p2at z0＝0.000 1～0.000 5 m

圖7是實測的活塞腔室壓力變化曲線與仿真曲線的對比，為清晰地反映各自變化情況，將曲線在圖中的時間錯開了0.002 s。忽略由于高頻測試帶來的噪聲信號和測量誤差，可知所建模型仿真的曲線與實測曲線匹配良好：它們都在經歷負開口區(qū)間時波動劇烈，而在穩(wěn)定區(qū)間，實測值在仿真曲線上下穩(wěn)定波動，它們的有效值基本保持一致。

4.2 不同負開口量對油壓脈沖的影響

負開口量的大小直接決定閥芯處與負開口區(qū)間的時間，也就直接了活塞腔室處于封閉狀態(tài)的時間。圖8（a）～（b）為不同的負開口量對壓力p1、p2油壓脈沖的影響對比。

圖8（a）為定義負開口量z0從0.0001 m到0.0005 m內變化時所得p1曲線，圖8（b）為對應的p2曲線。該圖反映出，當z0取值逐漸增大時p1的壓力峰值也會隨之增大，而且增長越來越劇烈，呈現(xiàn)二次增長趨勢；于此同時，p2的谷值則會隨之越來越小，并在沖程近半時出現(xiàn)壓力驟降的瞬間。

由此可知，z0取值越大時油液處于封閉狀態(tài)的時間就越長，油液壓縮性和膠管蓄油效應越不明顯，又由于閥芯在活塞回程和沖程過程中各要經歷一次負開口區(qū)間，故會出現(xiàn)壓力脈沖峰值更大和谷值壓力更小的現(xiàn)象。

4.3 油壓脈沖特性與系統(tǒng)壓力關系

進一步探索取不同值時油壓脈沖特性與系統(tǒng)穩(wěn)定壓力的關系，如圖9所示。

圖9 不同z0值對p與后腔穩(wěn)定壓力p1關系的影響Fig.9 Effects of relation between p and p1of different value of z0

圖9中，負開口量增大時，系統(tǒng)穩(wěn)定壓力雖有輕微增大，但大體不受影響，且波動趨勢不變。當活塞處于沖程階段時，由圖中可看出：后腔壓力與系統(tǒng)壓力基本保持一致，增減性也一致。由此可見，油壓的脈沖峰值與系統(tǒng)的最大穩(wěn)定壓力關系不大。

綜合上述各圖像的關系曲線可得表3所示后腔壓力p1和系統(tǒng)穩(wěn)定壓力隨負開口量變化的改變情況。表3更簡明地反映出，系統(tǒng)穩(wěn)定壓力和后腔穩(wěn)定壓力并不隨z0的變化而變化，z0主要影響p1壓力脈沖的最大最小值。表3中，峰壓率表示壓力峰值超過穩(wěn)定壓力的百分比。由表還可以看出：當z0≤0.000 3 m時，壓力峰值不會太大，峰壓率未超過18%，這是由于油液的可壓縮性和管道的油液補償特性在這一瞬間減緩了壓力沖擊，同樣也不會有壓力驟降產生，對機構本身不會產生有害作用，加之這個負開口區(qū)間避免了油液互混，能提高油液利用效率，于施工有利；而當z0＞0.000 3 m時，壓力峰值會迅速增大，超出油液可承受范圍，峰壓率將超過35%，而且會出現(xiàn)瞬間的壓力驟降，可能導致空穴產生，這是要避免的，不能用于工程設計。

表3 p1和系統(tǒng)壓力隨負開口量變化的改變情況Tab.3 Effects of p1of different value of z0

5 結 論

（1）負開口配油法縫隙節(jié)流過程復雜，在閥芯經歷負開口區(qū)間時需建立準負開口區(qū)間這一劃分，進行動態(tài)節(jié)流分析，進一步完善非線性動態(tài)模型。

（2）由于將油液本身存在可壓縮性和油路中的蓄能作用，在液壓油處于強烈節(jié)流狀態(tài)且壓力很大時，需將這些因素考慮進來。實測結果表明，這樣做是完全合理的。這不僅適用于液壓沖擊器，對于其他類似動態(tài)流體分析同樣適用。

（3）閥芯負開口量的大小是影響油液脈沖特性的關鍵因素，后腔的壓力峰值將隨著負開口量的增加而迅速增加，但若保證負開口量在一定范圍就不會造成有害的油液沖擊。通過對后腔壓力隨負開口量變化的研究得出，當負開口量取在z0≤0.0003 m時，活塞腔瞬間封閉不會導致有害的壓力沖擊和空穴，這樣的負開口閥配油既能避免高、低壓油液互混，又能確保沖擊器工作可靠。

［1］楊國平.液壓沖擊器換向閥設計方法［J］.建筑機械，2007（7）：75－77.

YANG Guo-ping.Design method of directional valve for hydraulic percussion device［J］.Architecture and Mechanism，2007（7）：75－77.

［2］趙宏強.液壓沖擊器配油閥設計［J］.鑿巖機械氣動工具，2005（2）：29－32.

ZHAO Hong-qiang.The design of distribution valve on hydraulic breaker［J］.Zaoyan Jixie Qidong Gongju，2005（2）：29－32.

［3］HE Qing-hua.Theoretical analysis and design/calculation formulate for hydraulic impact mechanism［J］.Transactions of NFsoc，1995，5（1）：161－121.

［4］Melamed Y，Kiselev A，Gelfgat M，etal.Hydraulic hammer drilling technology：developments and capabilities［J］.Journal of Energy Resources Technology，2000，122（3）：1－7.

［5］胡銘，董鑫業(yè).阿特拉斯公司鑿巖鉆車與液壓鑿巖機介紹［J］.鑿巖機械氣動工具，2011（4）：37－49.

HU Ming，DONG Xing-ye.Introduction of Atlas rock drilling rigs and hydraulic rock drills［J］.Zaoyan Jixie Qidong Gongju，2011（4）：37－49.

［6］朱鈺.液控換向閥內流場及動態(tài)特性的數(shù)值模擬［J］.哈爾濱工業(yè)大學學報，2012（5）：133－139.

ZHU Yu.Numerical simulation of flow field and dynamic characteristics of pilot operated directional valve［J］.Journal of Harbin Institute University，2012（5）：133－139.

［7］馬富銀，楊國平，吳偉蔚.液壓沖擊器流場分析［J］.振動與沖擊，2012，31（12）：73－78.

MA Fu-yin，YANG Guo-ping，WU Wei-wei.Fluid field analysis for a hydraulic impactor［J］.Journal of Vibration and Shock，2012，31（12）：73－78.

［8］何清華.液壓沖擊機構研究·設計［M］.長沙：中南大學出版社，2009：5－10.

［9］Lisowski E，Rajda J.CFD analysis of pressure loss during flow by hydraulic directional control valve constructed from logic valves［J］.Energy Conversion and Management，2013，65：285－291.

［10］Chiang L E，Elìas D A.A 3D FEMmethodology for simulating the impact in rock-drilli ng hammers［J］.International Journal of Rock Mechanics＆Mining Sciences，2008，45：701－711.

［11］章宏甲，黃誼.液壓傳動［M］.北京：機械工業(yè)出版社，2000：35－41.

［12］Jia Zhen-yuan，Ma Jian-wei，Wang Fu-ji，et al.Hybrid of simulated annealing and SVMfor hydraulic valve characteristics prediction［J］.Expert Systems with Applications，2011，38：8030－8036.

［13］Kilic E，Dolen M，Koku A B，et al.Accurate pressure prediction of a servo-valve controlled hydraulic system［J］.Mechatronics，2012，22：997－1014.

［14］Liu Wei，Yang H Y，Xu B，et al.Simulation study on the control valve in hydraulic operating mechanism of the high voltage circuit breaker［C］.ASME International Mechanical Engineering Congress and Exposition，Boston，MA，United states，Proceedings，2009，5：151－158.

［15］胡均平，郭勇，張政華，等.液壓打樁錘主控閥換向性能的動態(tài)特性分析［J］.中南大學學報（自然科學版），2012（7）：2622－2628.

HU Jun-ping，GUO Yong，ZHANG Zhen-hua，et al.Dynamic analysis ofmain control valve reversing performance for hydraulic pile hammer［J］.Journal of Central South University（Science and Technology），2012（7）：2622－2628.

Oil pressure pulse features of a hydraulic im pactor with an overlapped oil distributing valve

HU Jun-ping，HU Qian
（College of Mechanical and Electrical Engineering，Central South University，Changsha 410083，China）

Due to the fact that high pressure oil can mix together with low pressure oil using underlapped oil distributing valves widely used in domestic hydraulic impactors，the dynamic features of a hydraulic impactor with an overlapped oil distributing valve were studied and the optimal neutral position type valve was explored.According to the operating principle of hydraulic impactors and the transition principle of seam flow states，a nonlinear dynamic mathematicalmodel was established，and combining with Stateflow to simulate a frequent reversing process，the dynamic simulation model of a hydraulic impactor was obtained.Theworking pressure and oil flow features of a hydraulic impactor with an overlapped oil distributing value were investigated.A test was designed and its results were analyzed.The test results agreed well with those of simulation.The effects of different negative neutral open values on the dynamic chaber pressure of hydraulic impactorwere analyzed emphatically.The study resultsprovided a reference for designing the optimal neutral position type value of hydraulic impactors.

hydraulic impactor；oil distributing valve；negative neutral open value；oil pressure pulse

TH137

A

10.13465/j.cnki.jvs.2014.24.026

國家自然科學基金（51175518）

2013－08－09 修改稿收到日期：2013－12－12

胡均平男，博士，教授，博士生導師，1964年7月生

胡騫男，碩士生，1989年9月生