浙江大学学报(工学版), 2019, 53(7): 1265-1273 doi: 10.3785/j.issn.1008-973X.2019.07.004

机械与能源工程

小球式旋转直驱压力伺服阀卡滞机理研究

陆亮,, 夏飞燕, 訚耀保,, 原佳阳, 郭生荣

Spool stuck mechanism of ball-type rotary direct drive pressure servo valve

LU Liang,, XIA Fei-yan, YIN Yao-bao,, YUAN Jia-yang, GUO Sheng-rong

通讯作者: 訚耀保,男,教授. orcid.org/0000-0002-2671-6811. E-mail: y-yin@tongji.edu.cn

收稿日期: 2018-08-23  

Received: 2018-08-23  

作者简介 About authors

陆亮(1983—),男,助理教授,从事激振流体力学基础理论与机电液控制技术研究.orcid.org/0000-0002-9403-330X.E-mail:luliang829@tongji.edu.cn , E-mail:luliang829@tongji.edu.cn

摘要

针对小球式旋转直驱压力伺服阀(BRDDPSV)静态测试卡滞问题,建立阀芯运动全局函数,包括基于缝隙流理论建立倾斜阀芯径向力模型,基于Coulomb摩擦理论建立阀肩触壁静摩擦-滑动摩擦模型. 理论解析曲线合理复现了静态测试卡滞问题:偏心驱动下阀芯逆时针旋转倾斜,右侧阀肩触壁,初始静摩擦导致阀芯卡滞,逐渐提升的电流水平克服摩擦形成阀芯运动超调. 为了保证电流指令与控制压力的近似比例特性,阀芯回拉复位,形成重复的正向驱动阀芯卡滞. 基于阀肩不触壁原则,获得阀芯是否卡滞阈值条件. 研究结果表明:增大阀芯与阀套初始半径间隙或减小小球偏离阀芯轴线的初始偏心量,均可以提高阀芯不卡滞的输出压力阈值;对于21 MPa系统压力及0~8 MPa输出压力的实际需求,在不改变其他参数的情况下,将初始半径间隙和初始偏心距分别调整为5.1 μm和0.2 mm,可以在维持原有性能的基础上获得阀芯运动全局不卡滞的最优解.

关键词: 小球式旋转直驱 ; 压力伺服阀 ; 阀芯卡滞 ; 防卡滞参数优化

Abstract

Aiming at the sticking problem of ball-type rotary direct drive pressure servo valve (BRDDPSV), the global motion function of the spool was established, including the radial force model of inclined spool based on gap flow theory and the static-sliding friction between the spool shoulder and sleeve wall based on Coulomb friction theory. Sticking problem was reasonably reappeared by theoretical analysis curve: spool rotates and tilts counterclockwise by eccentric force, then the right shoulder of spool contact sleeve wall, initial static friction leads to spool sticking, gradually rising current overcomes friction to form spool movement overshoot. The valve core is pulled back to form a repetitive positive drive spool block in order to ensure the approximate proportion of current command and control pressure. The threshold condition of the sticking problem was obtained based on the principle of non-contact between spool shoulder and sleeve wall. Results show that increasing the initial radius gap between spool and bush or decreasing the initial eccentricity of ball from spool axis can increase the output pressure threshold of non-stuck. For the actual demand of 21 MPa system pressure and 0-8 MPa output pressure, the initial radius clearance and the initial eccentricity were adjusted to 5.1 μm and 0.2 mm respectively. The optimal solution of global non-sticking can be obtained by maintaining the original performance.

Keywords: rotary ball direct drive ; pressure servo valve ; spool stuck ; anti-stuck parameter optimization

PDF (1792KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

陆亮, 夏飞燕, 訚耀保, 原佳阳, 郭生荣. 小球式旋转直驱压力伺服阀卡滞机理研究. 浙江大学学报(工学版)[J], 2019, 53(7): 1265-1273 doi:10.3785/j.issn.1008-973X.2019.07.004

LU Liang, XIA Fei-yan, YIN Yao-bao, YUAN Jia-yang, GUO Sheng-rong. Spool stuck mechanism of ball-type rotary direct drive pressure servo valve. Journal of Zhejiang University(Engineering Science)[J], 2019, 53(7): 1265-1273 doi:10.3785/j.issn.1008-973X.2019.07.004

Moog[1]发明了喷嘴挡板型两级电液压力伺服阀,Howard[2]对两级阀设置一个与阀芯相连的拉簧直接作用于力矩马达完成机械反馈,改善了阀的死区,提高了压力阀的输出特性. Atchley[3]发明射流管压力伺服阀,相对于喷嘴挡板阀的双油路管道,降低了加工难度. 压力伺服阀由于动态响应快、控制精度高,在民用和军用伺服系统中广泛使用[4-6],但两级阀前置级组件复杂,制造成本高,且易被油液污染颗粒堵塞[7-9].

1970年以来,因电磁技术、功能材料的快速发展,出现了利用马达、磁致伸缩材料等直接驱动功率滑阀的单级电液伺服阀[10-14](direct drive valve,DDV),如MOOG公司的D633系列直接驱动伺服阀采用直线力矩马达驱动[15],但由于大体积的直线力矩马达与阀体同轴布置,导致整阀空间占用较大. 20世纪80年代末期,小体积旋转直接驱阀(rotary direct drive valve,RDDV)出现[16-20],该类阀采用与阀芯轴线垂直布置的高功率密度旋转电机,取代体积较大的直线力矩马达,使得整阀的占用空间大大减小.

研究的小球式旋转直驱压力伺服阀通过2处偏心设置,将电机的旋转运动转化为阀芯的直线运动进而控制压力,在减小空间占用的同时,结构简单、不易堵塞、泄漏量小. 静态特性测试表明,当输入斜坡电流信号时,样机控制压力重复出现无响应又突然超调的问题;因阀芯运动开度与控制压力近似比例,可以推测阀芯运动出现阶段性卡滞,并在突然释放后超调运动,本文称为阀芯运动的卡滞与猛冲现象或简单叙述为阀芯卡滞.

液压阀芯的卡滞问题是液压阀的主要故障之一. 冀宏等[21]针对固体颗粒造成的液压卡滞问题进行一系列研究,分析液压滑阀中μm级颗粒的沉积情况以及颗粒侵入阀芯阀套配合间隙后的力学行为. 刘新强等[22-23]基于热力耦合变形理论,计算滑阀工作过程中的节流发热以及阀芯阀套的变形,分析因结构热变形而发生的阀芯卡死问题.

本文所研究的旋转直驱压力伺服阀,特殊的小球式偏心驱动结构使得滑阀阀芯受到倾覆力矩和不均衡径向液力作用,进而造成阀芯运动过程中出现倾斜甚至阀肩触壁的问题. 基于缝隙流理论建立倾斜阀芯径向力模型,基于Coulomb摩擦理论建立阀肩触壁静摩擦-滑动摩擦模型,进而建立阀芯运动全局函数,以求完整复现阀芯运动卡滞机理,提出防止卡滞的参数优化设计方法.

1. 工作原理与阀芯卡滞

1.1. 工作原理

小球式旋转直驱压力伺服阀由电机、阀芯、阀套、壳体、压力传感器及反馈电子控制器构成,如图1(a)所示. 电机输出轴偏离电机轴线e0,末端固结驱动小球,阀芯末端加工柱形孔,柱形孔内径与小球外径的名义尺寸相等. 将小球置于柱形孔内,且球心位于阀芯轴线上方h0h0>0),则该小球-柱形孔运动副可以实现将电机的旋转运动转化为阀芯的水平直线运动和阀芯绕轴线的转动,小球作用在阀芯上的作用力可以分解为水平驱动力Fx和转动力矩Tβ. 小球-柱形孔运动副的结构如图1(b)(c)所示.

图 1

图 1   小球式旋转直驱压力伺服阀结构原理

Fig.1   Structure and working principle of BRDDPSV


该小球式旋转直驱压力伺服阀的滑阀副为两凸肩、三通的正开口结构形式,负载口置于供油口和回油口之间,阀芯水平运动形成供油口和回油口间联动的节流作用,从而控制中间负载口的输出压力,最终目的是使得该阀能够输出与输入电流指令成比例的压力. 该阀的结构与工质参数如表1所示.

表 1   小球式旋转直驱压力伺服阀结构工质参数

Tab.1  Parameters of BRDDPSV structure and working medium

参数 参数值
电机电阻Rc 23.3
阀芯质量mv/g 5
阀芯总长度l/mm 24
阀芯直径Dv/mm 6
阀肩长度L/mm 6.4
小球偏离阀芯轴线距离h0/mm 0.5
阀芯与阀套间的初始半径间隙δ/μm 3
滑阀预开口量U/mm 0.1
负载容腔容积V/mL 10
工作介质 15#液压油
介质温度θ/°C 25±0.5

新窗口打开| 下载CSV


1.2. 静态测试

设计要求21 MPa系统压力经过该压力伺服阀后,输出0~8 MPa的控制压力[24]. 为了测试该小球式旋转直驱压力伺服阀的静态特性,搭建图2所示的测试台. 该测试台可以提供21 MPa油源压力,性能测试设备包括测试PC机和信号发生器,内置伺服阀动静态特性测试软件;数字控制PC为便携式笔记本;控制PC机通过信号采集、输出及放大设备与伺服阀和信号发生器相连. 压力测量采用量程为0~25 MPa的固态压阻式压力传感器.

图 2

图 2   压力伺服阀静态特性测试实验台

Fig.2   Test bench for static characteristic test of pressure servo valve


设计目的是获得指令电流、阀芯运动开度与控制压力的近似比例关系,但实际上静态测试效果不理想. 如图3所示,在阀芯开启瞬间,出现压力死区,紧接着出现压力上升,并恢复至拟定值(①-②);恢复后进入新的压力死区,并且在电流增至一定程度后,压力出现超调(②-③),因控制压力与阀芯运动开度近似比例,故称为阀芯卡滞-猛冲区域;压力超调后又拉回至拟定值附近,出现新一轮“卡滞-猛冲”(③-④). 在阀芯闭合阶段(⑤-⑥),出现卡滞问题,但严重程度相对开启阶段较弱.

图 3

图 3   静态测试实验曲线

Fig.3   Experimental curve of static test


控制压力静态测试结果的复杂变化是源于阀芯的不规则运动导致的开度复杂变化,阀芯的不规则运动是源于偏心驱动运动过程中耦合液动力、阀芯倾斜不均衡径向力与壁面接触摩擦力的共同作用. 考虑在整阀数学模型的基础上,细化阀芯动力学模型中的阀腔液动力、缝隙流动径向力、阀肩触壁摩擦力模型,以揭示阀芯运动卡滞机理.

2. 卡滞机理建模复现

2.1. 数学建模

小球式旋转直驱压力伺服阀的系统框图如图4所示. 该系统包括2个闭环:一是电机转角闭环,目的是提高电机抗扰动能力;二是控制压力闭环,目的是使输出压力pL与指令电流信号i成比例. 经PID控制器运算后的信号输入电机,电机末端的小球将力和力矩传递至阀芯,滑阀的水平位移在阀口处形成节流口进而输出一定的压力,该压力由压力传感器测得后反馈至输入端控制电机电压. 该阀的数学模型由电机驱动方程、动力转换方程、阀芯动力学方程、负载压力方程及反馈控制方程组成.

图 4

图 4   小球式旋转直驱压力伺服阀系统框图

Fig.4   System block diagram of BRDDPSV


2.1.1. 电机驱动方程

电机转角限制在±30°内. 当输入电机的电压为u时,电机的输出电磁力矩为Tem,转过相应的角度α,进而推动阀芯运动. 阀芯作用在电机转轴上的阻力矩为Tf. 电机特性测试表明,Temα近似呈抛物线关系,可以表示为

$ {T_{{\rm em}}} = {k_{{\rm t}}}u - {k_{{\rm m}}}{\alpha ^2}. $

式中: k t为电流力矩系数, k m为转角系数. 电机转子的动力学方程为

$ {J_{{\rm r}}}\frac{{{{{\rm d}}^2}\alpha }}{{{\rm d}{t^2}}} = {T_{\rm em}} - {B_{\rm r}}\frac{{{\rm d}\alpha }}{{{\rm d}t}} - {T_{\rm f}}. $

式中:Jr为转动惯量,Br为阻尼系数.

2.1.2. 小球-柱形孔运动副动力转换方程

小球-柱形孔运动副的动力转换关系如图1(b)所示. 建立图1(b)所示的坐标系,由于电机和阀套均固定在阀体上,以阀套轴线为X轴,电机轴线为Y轴,垂直于XY轴且过两轴交点的直线为轴. 小球所在轴线的初始位置位于YOZ平面内.

当电机输出转角为α时,由驱动小球和柱形孔的传动方式可知,阀芯水平位移xv和绕阀芯轴线的转角βv可以表示为

$ {x_{\rm v}} = {e_0}\sin \alpha , $

$ {\beta _{\rm v}} = \arcsin \;({{{e_0}(1 - \cos \alpha )}}/{{{h_0}}}), $

因此电机阻力矩为

${T_{\rm{f}}} = {F_x}{e_0}\cos \alpha + ({T_\beta }\cos \beta _{\rm{v}}/{h_0}){e_0}\sin \alpha. $

2.1.3. 阀芯动力学方程

一般情况下,阀芯运动处于水平状态,除驱动力外主要考虑运动阻尼与液动力;本文的小球式直驱阀受偏心驱动,不可避免地发生倾斜,须细化考虑偏心导致的不对称缝隙流动径向力以及阀肩发生触壁后形成的摩擦力.

图5所示为XOY平面内,阀芯倾斜后的受力示意图. 偏心驱动点位于阀芯右侧轴线上方,故最初于右侧阀肩与阀套壁面接触时出现摩擦力f,阀芯运动除正向开启和反向闭合外,还可能因卡滞后超调运动并拉回,摩擦力方向存在变化;事实上因阀芯的倾斜运动,于左侧阀肩同样可能出现摩擦力. 仅左侧阀肩存在压差流,故径向力Fc出现于左侧阀肩,且与阀芯的倾斜运动趋势相关,Fc的作用方向存在变化.

图 5

图 5   阀芯受力示意图

Fig.5   Schematic diagram of force acting on spool


由牛顿第二定律可知,阀芯动力学方程如下.

阀芯在水平驱动力Fx作用下克服阻尼力、液动力及摩擦力运动,控制阀口开度xv

$m\frac{{{{\text{d}}^{\text{2}}}x_{\rm v}}}{{{\text{d}}{t^{\text{2}}}}}{\text{ = }}{F_x} - {B_{x{\text{v}}}}\frac{{{\text{d}}{x_{\text{v}}}}}{{{\text{d}}t}} - F{\rm _s} - f.$

式中:Bxv为阀芯水平运动黏性阻尼,m为阀芯质量,Fs为阀芯受到的液动力,f为阀芯与阀套接触时的摩擦力.

阀芯在驱动力矩Tβ作用下克服转动阻尼Bβv,使阀芯绕轴线转动βv

${J_{\beta {\text{v}}}}\frac{{{{\text{d}}^{\text{2}}}\beta_{\rm v}}}{{{\text{d}}{t^{\text{2}}}}} = {T_\beta } - {B_{\beta {\text{v}}}}\frac{{{\text{d}}\beta _{\text{v}}}}{{{\text{d}}t}}.$

式中:Jβv为均质阀芯的转动惯量.

无支承的阀芯在FxFc作用下倾斜转动γv

${J_{\text{z}}}\frac{{{{\rm d}^2}{\gamma _{\rm{v}}}}}{{{\rm d}{t^2}}} = {F_x}{h_0} - {F_{\rm c}}\left(a - \frac{L}{2}\right) - {B_{\gamma {\rm{v}}}}\frac{{{\rm d}{\gamma _{\rm{v}}}}}{{{\rm d}t}} - {k_{\rm{n}}}(l - a) - f\frac{D}{2}.$

式中:aLlD为阀芯的结构尺寸,分别表示左阀肩至质心的距离、左阀肩的长度、阀芯长度、阀肩直径; $B{\gamma _{\rm{v}}} $为阀芯倾斜运动的黏性阻尼;Jz为均质阀芯绕质心倾斜的转动惯量;γv为阀芯倾斜角度;kn为阀芯与阀套接触挤压的正压力.

阀芯在Fc和重力的作用下克服径向阻尼力,使阀芯上下运动yv

$m\frac{{{{\text{d}}^{\text{2}}}{y_{\rm v}}}}{{{\text{d}}{{{t}}^{\text{2}}}}}= {F_{\rm c}} - mg - {B_{y{\text{v}}}}\frac{{{\text{d}}{{{y}}_{\text{v}}}}}{{{{\text d{t}}}}} - k_{\rm n}.$

式中:yv为阀芯径向位移.

1)液动力. 阀芯动力学方程中的液动力Fs是指阀腔内流动的油液作用在阀芯内侧端面上的作用力. 液动力主要为稳态液动力,表达式[26]

$\begin{aligned} & F_{\rm s}= \\ & 2C_{\rm v}C_{\rm d}W\cos \theta x_{\rm v}\left( {p_{\rm s}\! -\! p_{\rm L}} \right) \!-\!\! 2C_{\rm v}C_{\rm d} W\cos \theta \left( {U - x_{\rm v}} \right)p_{\rm L} = \\ & 2C_{\rm v}C_{\rm d}W\cos \theta p_{\rm s}\frac{{x_{\rm v}(U - x_{\rm v})(U - 2x_{\rm v})}}{{{{(\sqrt 2 x_{\rm v} - U)}^2} + 2(\!\!\sqrt 2 - 1)x_{\rm v}U}}. \end{aligned} $

式中:θ为射流角;Cv为速度系数;Cd为流量系数;W为面积梯度,WDU为滑阀预开口量;ps为系统压力;pL为负载压力.

由式(10)可以看出,液动力随ps的增大而增大,且当0<xv<U/2时,Fs>0;当U/2<xv<U时,Fs<0;当xv=U/4时Fs取极大值,xv=3U/4时Fs取极小值. 在阀芯启闭过程中,液动力的大小和方向均可能发生改变.

2)径向力. 阀芯动力学方程中的Fc是指在阀肩与阀套的缝隙处,由于泄漏流量在缝隙中存在压差,阀肩周向缝隙宽度不等导致阀肩径向分布大小不等的压力,这些压力可以等效替代为一集中力的形式,以便表达和计算. Fc由缝隙流压力分布在圆周方向积分得出.

对于呈微小夹角的倾斜平板而言,单宽流道中流体的作用力[25]

$ F = \int\limits_0^L {p{{\rm{d}}}x = \left\{ {{p_1} - \frac{{{h_2}}}{{{h_1} + {h_2}}}\left( {{p_1} - {p_2}} \right)} \right\}} L .$

式中:p1p2分别为缝隙入口和出口压力,h1h2分别为入口和出口缝隙高度,L为缝隙长度.

对于无支承的阀芯来说,取其绕质心逆时针旋转为正方向,向上浮动为正方向,则任一时刻阀芯倾斜γv、浮动yv时的阀肩和阀套位置关系如图6所示.

图 6

图 6   阀芯倾斜时阀肩和阀套的位置关系

Fig.6   Position relationship between valve shoulder and sleeve on spool leaning condition


图6可知,

$\begin{array}{*{20}{l}} {{h_1} = \delta + \left( {{e_1} - y_{\rm v}} \right)\cos \theta ,\;{h_2} = \delta + \left( {{e_2} - y_{\rm v}} \right)\cos \theta .} \end{array}$

式中:e1e2分别为阀芯倾斜导致的缝隙入口和出口的偏心值,δ为阀芯阀套间的半径间隙. 对式(12)沿阀肩周向积分,可得缝隙中的流体在阀肩上的作用力,取向上为正方向,则

$ F_{\rm c} = \frac{{2\pi \delta LR \Delta p ({e_{\text{1}}} - {e_{\text{2}}})}}{{{{({e_1} + {e_2} - 2y_{\rm v})}^2}}}\left[ {1 - {{{\left( 1 - {\frac{{{e_1} + {e_2} - 2{y_{\rm{v}}}}}{{2\delta }}} \right)}^{ - {1}/{2}}}} } \right].$

式中:R为阀肩半径;Δp为阀肩两端压差,对于左阀肩,Δp=ps.

对于左端阀肩,有

$\begin{array}{l} {e_1} = (a - L){\gamma _{\rm v}}, \; {e_2} = a{\gamma _{\rm v}} . \end{array} $

由于右阀肩两端均为回油压力,对于右端阀肩,有Δp=0,则Fc=0.

代入阀芯结构参数,可得Fcγvyv的关系,如图7所示.

图 7

图 7   缝隙压差流引起的径向力与阀芯径向位移和倾斜角度的关系

Fig.7   Differential pressure slit-flow produced radial force vs. spool radial displacement and leaning angle


图7的径向力求解结果可见,当阀芯逆时针旋转时,作用在左阀肩的径向力向上,反之向下,即径向力的作用可以减弱阀芯倾斜的趋势. 倾斜角度对径向力起主导作用,径向位移对径向力影响较弱.

3)摩擦力. 当阀芯倾斜后,在逆时针旋转力矩下发生阀肩与阀套触碰挤压,如图8所示.

图 8

图 8   阀肩触壁摩擦受力示意图

Fig.8   Sketch of friction when valve shoulder touching sleeve wall


由于阀芯和阀套间的接触为小变形,可以用等效弹簧阻尼来建立接触模型,kn由Hertz定律[27-28]确定.

$\left. \begin{array}{l} {k_{\rm n}} = {F_{\rm k}} + {F_{\rm d}} = {K_{\rm n}}{\xi ^m} + C\xi \displaystyle\frac{{{\rm d}\xi }}{{{\rm d}t}}; \\ \xi = \left\{ \begin{array}{l} \left| {(l - a){\gamma _{\rm v}} + {y_{\rm v}} - \delta } \right|,\,\,(l - a){\gamma _{\rm v}} + {y_{\rm v}} - \delta > 0; \\ 0,\;(l - a){\gamma _{\rm v}} + {y_{\rm v}} - \delta < 0. \\ \end{array} \right. \end{array} \right\} $

式中:Fk为等效弹簧力;Fd为等效阻尼力;Kn为等效弹簧刚度;ξ为接触挤压深度,可以由γvyv表示;C为阻尼系数,与恢复系数和接触刚度等有关,

$ C = 0.75\left( {1 - {\lambda ^2}} \right){K_{\rm n}}{\xi ^{1.5}}/{v_0}, $

其中λ为恢复系数,通常取0.5;v0为碰撞接触点的初始相对速度.

由Hertz接触理论可知,当阀芯与阀套碰撞接触时,接触表面近似为椭圆,Kn可以表示为

$ \left.\begin{array}{l} {K_{\rm n}} = \displaystyle\frac{4}{{3\pi \left( {{\sigma _i} + {\sigma _j}} \right)}}{\left[ {\displaystyle\frac{{{R_{\rm i}}{R_{\rm j}}}}{{{R_{\rm i}} - {R_{\rm j}}}}} \right]^{\frac{1}{2}}},\\ {\sigma _{\rm i}} = \displaystyle\frac{{1 - v_{\rm i}^2}}{{\pi {E_{\rm i}}}},\\ {\sigma _{\rm j}} = \displaystyle\frac{{1 - v_{\rm j}^2}}{{\pi {E_{\rm j}}}}. \end{array}\right\} $

式中:RiRj分别为阀套和阀芯半径,vivj分别为阀套和阀芯的泊松比,EiEj分别为阀套和阀芯的弹性模量.

根据Coulomb摩擦模型可知,当阀芯与阀套接触时,产生的摩擦力f可以表示为阀芯水平运动速度的分段函数:

$ f = \left\{ \begin{array}{l} {\mu _{\rm n}}{k_{\rm n}}{\mathop{\rm sgn}} \left(\displaystyle\frac{{{\rm d}{x_{\rm v}}}}{{{\rm d}t}}\right),\;\displaystyle\frac{{{\rm d}{x_{\rm v}}}}{{{\rm d}t}} \ne 0;\\ u(t),\;\displaystyle\frac{{{\rm d}{x_{\rm v}}}}{{{\rm d}t}} = 0, \left| {u(t)} \right| < {f_{\rm{s}}};\\ {f_{\rm s}},\;{\text{其他}}. \end{array} \right. $

式中:μn为阀芯阀套间的摩擦系数,ut)为合外力,fs为最大静摩擦力.

在获得理论分析结果之前,通过Coulomb摩擦模型可以初步推测,当阀芯倾斜至右侧阀肩触壁后,若驱动力不足以克服摩擦力,将出现静摩擦,阀芯卡滞;因控制策略要求输出与指令电流匹配的阀口开度及输出压力,则会在卡滞状态下不断提高电流水平,一旦驱动力克服摩擦力,静摩擦变为滑动摩擦,阻力变小,则进一步形成阀芯运动超调;在同样的控制要求下,过度的开度将形成反向电流而拉回阀芯. 可见摩擦是阀芯卡滞-猛冲的主要原因,是否形成摩擦以及是否是静摩擦或滑动摩擦,则依赖于驱动力、液动力、径向力耦合作用下阀芯姿态的变化.

2.1.4. 负载压力方程

正开口滑阀副的预开口量为U,滑肩直径为D. 阀芯水平运动位移xv后,在供油口和回油口处形成联动节流作用. 当供油口输入体积弹性模量为E的恒压ps油源时,pL可以由阀腔内液体流动的连续性及容积V的负载腔内油液压缩性求得:

${C_{\rm d}}\pi D{x_{\rm v}}\sqrt {\frac{{2({p_{\rm s}} - {p_{\rm L}})}}{\rho }} - {C_{\rm d}}\pi D(U - {x_{\rm v}})\sqrt {\frac{{2{p_{\rm L}}}}{\rho }} = \frac{V}{E}\frac{{{\rm d}{p_{\rm L}}}}{{{\rm d}t}}.$

一般来说,对于不考虑油液可压缩性的稳态过程,有

${p_{\rm L}} = \frac{{{x_{\rm v}}^2}}{{2{x_{\rm v}}^2 - 2{x_{\rm v}}U + {U^2}}}{p_{\rm s}}.$

由式(20)可以看出,当ps一定时,pLxv的增大而单调增大,可以通过控制压力变化情况大致推断阀芯的水平运动情况.

2.1.5. 反馈控制方程

该小球式旋转直驱压力伺服阀采用闭环压力反馈的控制方法,如图4所示. 将pL与指令电流i的偏差信号进行比例KP、积分KI运算(PI),输出电机的控制信号u. 该压力闭环反馈的控制方法常规,且易于搭建模拟电路. 压力反馈系数及电机转角反馈系数分别为kf1kf2.

控制方程可以表示为

$u = {k_{\rm pwm}}\left[ {\left( {{k_{\rm b}}i - {k_{{\rm f}1}}{p_{\rm L}}} \right)\left( {{K_{\rm P}} + {K_{\rm I}}\frac{1}{s}} \right) - \alpha {k_{{\rm f}2}}} \right].$

式中:kpwm为驱动信号放大比例系数,kb为输入电流信号调理比例系数.

2.2. 理论分析

借助Simulink工具对式(1)~(21)进行计算分析,如图9所示. 理论结果较合理地反映了实验曲线的特征,理论曲线的压力死区和超调量相对实验值较小,可能源于实际油液阻尼和摩擦系数大于理论选取值.

图 9

图 9   控制压力理论计算结果与实验结果对比

Fig.9   Control pressure comparison between theoretical computation and experimental test


图 10

图 10   理论求解阀芯水平运动速度曲线

Fig.10   Theoretical curve of spool horizontal velocity


为了进一步说明在复杂耦合力及反馈控制作用下的阀芯运动特性,图10给出阀芯水平速度曲线,图11给出阀芯右端面中心点的径向位移. 结合图9可以看出多次卡滞过程中控制压力、阀芯水平运动速度及右端面径向位移的对应关系,即当阀芯发生卡滞时,右侧阀肩触壁,水平速度为零,控制压力恒定;当卡滞解除,过调的指令电流导致阀芯“猛冲”,速度超调,控制压力超调;为了保证指令电流与控制压力的近似比例关系,阀芯拉回,速度为负,控制压力向拟定值恢复.

图 11

图 11   理论求解阀芯径向位移曲线

Fig.11   Theoretical curve of spool radial displacement


综上所述,小球式旋转直驱压力伺服阀阀芯卡滞机理主要在于开启过程中偏心驱动下阀芯不可避免地发生逆时针旋转倾斜,导致右侧阀肩触壁,初始静摩擦导致阀芯卡滞;不断提升的电流水平克服静摩擦,一旦卡滞解除则形成较低水平的滑动摩擦或不再触壁,阀芯“猛冲”,但因指令电流与输出压力近似比例控制的需要,驱动力反向迅速拉回阀芯;指令电流与控制压力的关系恢复正常后,开始新一轮的正向驱动、逆时针倾斜、阀芯卡滞、驱动力超调、卡滞解除、压力超调而又阀芯拉回的情况. 在阀芯闭合过程中,因驱动力作用点位于阀芯右侧的轴线上方,虽然出现了一定程度的卡滞,但严重程度低于开启阶段.

3. 防止卡滞参数优化

3.1. 压力参数匹配

阀芯在阀套内出现首次卡滞,在于右侧阀肩与阀套壁面接触,则以接触卡滞即h2=0作为条件进行模型求解,令下式为零:

${h_2} = \delta - (l - a){\gamma _{\rm v}} - {y_{\rm v}}.$

计算结果如图12所示. 由式(22)为零获得现有结构参数条件下的系统压力与控制压力关系曲线,如图12中的“点-线”所示,左侧为理论获得的不发生卡滞区间,右侧为理论上的卡滞区间. 对该假设开展3组实验,ps分别为21、16和8 MPa,对应的阈值控制压力pL0分别为3.2、4.1和8 MPa,实验阈值点落在理论获得系统压力vs.阈值控制压力曲线附近,但存在一定的误差,3种系统压力工况下的控制压力实验与理论误差分别为14.3%、12.5%和13.8%,误差较大的原因可能在于读表误差和现场实验的较大干扰因素,但试验与理论结果的基本规律是一致的;可见理论假设符合预期. 理论判断不发生卡滞是基于不发生首次卡滞,而实验是从整个阀芯运动过程中获得的控制压力阈值点,但通过图911的分析可知,首次卡滞最严重,后续卡滞在原理上与首次相同且具重复性,由此可知避免首次卡滞,能够在全局上避免卡滞.

图 12

图 12   阀芯卡滞压力阈值理论与实验对比

Fig.12   Pressure thresholds comparison between theory and experiment


图12的分析结果可得如下规律:系统压力越低,不发生卡滞的控制压力裕度越大,但对于实际系统,不能以牺牲系统压力为代价获得更大的不卡滞区间,需要从结构参数上获得更可行的参数优化方法.

3.2. 结构参数优化

阀结构参数优化比降低系统压力更可行,但须考虑由此带来的负面影响. 针对阀芯与阀套间的初始半径间隙δ及驱动力作用点偏离阀芯轴线的初始距离h0两个可调的重要结构参数进行考虑,系统压力设定为所需的21 MPa,控制压力为所需的0~8 MPa.

1)阀芯与阀套间的初始半径间隙. 由式(22)可知,δ是影响阀芯卡滞与不卡滞的临界判断条件何时达成的重要参数. 图13针对δ在现有的3 μm基础上进行调整,获得不同δ情况下,阀芯是否发生卡滞的区域划分,每条曲线的左侧为阀芯不卡滞区间.

图 13

图 13   不同初始半径间隙下的阀芯卡滞与否区间划分

Fig.13   Section partition for stuck or not on different initial radius gap conditions


图13可以看出,增大间隙可以显著增大阀芯不卡滞区域的范围. 例如,将原先的3 μm间隙增大至9 μm后,可以实现系统压力21 MPa,控制压力在0~8 MPa全局范围内阀芯均不发生卡滞,但在一定程度上导致泄漏量增大,泄漏量的表达式为

$q = \frac{{\pi d{\delta ^3}{p_{\rm s}}}}{{12\eta L}}(1 + 1.5{\varepsilon ^2}).$

式中:ε为阀芯在阀套内的偏心系数,取最大偏心系数ε=1考虑,当间隙由3 μm增至9 μm后,泄漏量能够在原来基础上最大增加26倍.

2)驱动力作用点偏离阀芯轴线的初始距离.

h0是重要的设计参数,若h0=0,该小球-柱形孔运动副卡死不工作,即h0存在下限. 令式(22)为零,获得不同h0下的卡滞判断,如图14所示,每条曲线的左侧为阀芯不卡滞区间.

图 14

图 14   不同初始偏心量下的阀芯卡滞与否区间划分

Fig.14   Section partition for stuck or not on different initial eccentricity conditions


图14可以看出,随着h0的减小,阀芯不发生卡滞的范围逐渐增大,当h0减小至0.08 mm时,理论上可以保证当系统压力为21 MPa时,控制压力在0~8 MPa全局范围内阀芯不卡滞. 由式(5)可知,减小h0会增大Tf,增加电机负荷,引起电机驱动能力不足.

3)全局不卡滞的关键参数匹配关系. 单一改变δh0,尽管都能够在理论上获得全局(系统压力为21 MPa,控制压力为0~8 MPa)不卡滞的设计依据,但都不可避免地带来较大的负面影响,如过度增加泄漏或电机驱动功率.

以全局(系统压力为21 MPa,控制压力为0~8 MPa)不卡滞为优化设计目标,将δh0均作为变量,以式(22)为零建立关键结构参数的阈值关系,如图15所示,曲线右侧为全局不卡滞区间. 在考虑泄漏损失和功率损失的基础上,权衡关键参数的选取. 举例说明,当初始偏心距由原先的0.5 mm减小至0.2 mm时,只需将δ增大到5.1 μm即可正常使用,此时泄漏量仅比原来增加3.9倍.

图 15

图 15   全局不卡滞的结构参数匹配(ps=21 MPa,pL0=8 MPa)

Fig.15   Structural parameters matching for no-stuck(ps=21 MPa,pL0=8 MPa)


表1所示,除初始偏心距与阀芯阀套初始半径间隙外,还可以针对其他参数的调整作更多选择,获得更合理的改进办法.

4. 结 论

(1)针对小球式旋转直驱压力伺服阀静态测试卡滞问题,在建立电机驱动方程、动力转换方程、阀芯动力学方程、负载压力方程及反馈控制方程全局数学模型的基础上,完善阀芯运动过程液动力、径向力及摩擦力模型,理论结果复现了静态特性阀芯卡滞问题.

(2)研究表明,阀芯卡滞机理主要在于偏心驱动下阀芯不可避免发生逆时针旋转倾斜,导致右侧阀肩触壁,初始静摩擦导致阀芯卡滞,逐渐提升的电流水平克服摩擦形成阀芯运动开度、运动速度与控制压力超调;为了保证电流指令与控制压力的近似比例特性,阀芯回拉复位,形成重复的正向驱动阀芯卡滞.

(3)以不发生右侧阀肩首次触壁为原则,建立伺服阀压力与结构参数阈值的关系式,获得不发生卡滞的参数区间. 在保证全局正常工作压力的前提下,建立不发生卡滞的关键结构参数匹配关系,为优化设计提供了依据.

参考文献

MOOG W C. Electrohydraulic servo mechanism: US2625136 [P]. 1953-01-13.

[本文引用: 1]

HOWARD C T. Flow control servo valve: US 2790427 [P]. 1957-04-30.

[本文引用: 1]

ATCHLEY R D. Valve: US, US3017864 [P]. 1962-01-23.

[本文引用: 1]

MASKREY R H, THAYER W J

A brief history of electrohydraulic servomechanisms

[J]. Journal of Dynamic Systems Measurement and Control, 1978, 100 (2): 110

DOI:10.1115/1.3426352      [本文引用: 1]

BRADLEY D A. Mechatronics: electronics in products and processes [M]. New York: Routledge, 2018.

EL-ARABY M, EL-KAFRAWY A, FAHMY A

Dynamic performance of a nonlinear non-dimensional two stage electrohydraulic servo valve model

[J]. International Journal of Mechanics and Materials in Design, 2011, 7 (2): 99- 110

DOI:10.1007/s10999-011-9150-x      [本文引用: 1]

SHARMA S K, JR C E S, GSCHWENDER L J, et al

Stuck servo valves in aircraft hydraulic systems

[J]. Lubrication Engineering, 1999, 55 (7): 27- 32

[本文引用: 1]

YIN Y B, YUAN J Y, GUO S R

Numerical study of solid particle erosion in hydraulic spool valves

[J]. Wear, 2017, 392-393: 174- 189

DOI:10.1016/j.wear.2017.09.021     

訚耀保. 极端环境下的电液伺服控制理论及应用技术[M]. 上海: 上海科学技术出版社, 2012.

[本文引用: 1]

AUDREN J T, MERLET E, MELEARD J, et al. Valve control device: US, US7026746 [P]. 2006-04-11.

[本文引用: 1]

SONG M G, HUR Y J, PARK N C, et al

Design of a voice-coil actuator for optical image stabilization based on genetic algorithm

[J]. IEEE Transactions on Magnetics, 2009, 45 (10): 4558- 4561

DOI:10.1109/TMAG.2009.2021399     

俞军涛, 焦宗夏, 吴帅

基于液压微位移放大结构的新型压电陶瓷直接驱动阀设计及仿真

[J]. 机械工程学报, 2013, 49 (2): 151- 158

DOI:10.3969/j.issn.1674-8530.2013.02.013     

YU Jun-tao, JIAO Zong-xia, WU Shuai

Design and simulation study on new servo valve direct driven by piezoelectric actuator using hydraulic amplification

[J]. Chinese Journal of Mechanical Engineering, 2013, 49 (2): 151- 158

DOI:10.3969/j.issn.1674-8530.2013.02.013     

WU S, JIAO Z, YAN L, et al

Development of a direct-drive servo valve with high-frequency voice coil motor and advanced digital controller

[J]. IEEE/ASME Transactions on Mechatronics, 2014, 19 (3): 932- 942

DOI:10.1109/TMECH.2013.2264218     

YANG Z, HE Z, LI D, et al

Direct drive servo valve based on magnetostrictive actuator: multi-coupled modeling and its compound control strategy

[J]. Sensors and Actuators A Physical, 2015, 235: 119- 130

DOI:10.1016/j.sna.2015.09.032      [本文引用: 1]

LEI H, HU X

Research on application of D633 series direct drive servo-proportional control valves

[J]. Machine Tool and Hydraulics, 2007, 35 (9): 171- 172

[本文引用: 1]

JOHNSON D D. Direct drive servo valve with rotary valve: US, US 4800924 A [P]. 1989-01-31.

[本文引用: 1]

LAUX K. Motor-to-spool coupling for rotary-to-linear direct drive valve: US, US5263681 [P]. 1993-11-23.

LEONARD M B. Rotary servo valve: US, US 5954093 A [P]. 1999-09-21.

原佳阳, 訚耀保, 陆亮, 等

旋转直接驱动式电液压力伺服阀机理及特性分析

[J]. 机械工程学报, 2018, 54 (16): 186- 194

YUAN Jia-yang, YIN Yao-bao, LU Liang, et al

Analysis of rotary direct drive electro-hydraulic pressure control servo valve

[J]. Chinese Journal of Mechanical Engineering, 2018, 54 (16): 186- 194

ROBERT D, JOHN W. Direct dive valve-ball drive mechanism: US, US4672922 A [P]. 1987-06-16.

[本文引用: 1]

冀宏, 崔腾霞, 陈晓明, 等

滑阀间隙中方形微米颗粒的旋转现象

[J]. 兰州理工大学学报, 2017, 43 (5): 51- 55

DOI:10.3969/j.issn.1673-5196.2017.05.009      [本文引用: 1]

JI Hong, CUI Teng-xia, CHEN Xiao-ming, et al

Rotation phenomenon of square micron particles in clearance of spool valve

[J]. Journal of Lanzhou University of Technology, 2017, 43 (5): 51- 55

DOI:10.3969/j.issn.1673-5196.2017.05.009      [本文引用: 1]

刘新强. 液压滑阀热力耦合形变与污染卡滞机制的研究[D]. 兰州: 兰州理工大学, 2013.

[本文引用: 1]

LIU Xin-qiang. The research of thermal-mechanical coupling and sticking mechanism for contamination of hydraulic spool valve [D]. Lanzhou: Lanzhou University of Technology, 2013.

[本文引用: 1]

刘晓红, 柯坚, 刘桓龙

液压滑阀径向间隙温度场的CFD研究

[J]. 机械工程学报, 2006, 42 (b05): 231- 234

[本文引用: 1]

LIU Xiao-hong, KE Jian, LIU Heng-long

CFD research on temperature field in radial clearance of hydraulic spool valve

[J]. Chinese Journal of Mechanical Engineering, 2006, 42 (b05): 231- 234

[本文引用: 1]

LI B, JIAO Z, WU S, et al. Research on a new direct drive electro-hydraulic proportional pressure valve for aircraft braking system [C] // IEEE International Conference on Aircraft Utility Systems. Beijing: IEEE, 2016: 1192–1196.

[本文引用: 1]

WANG D

The force analysis of the phenomind of hydraulic block on hydraulic spool valve

[J]. Journal of Shenyang University: Natural Science Edition, 1994, (4): 76- 78

[本文引用: 1]

AMIRANTE R, MOSCATELLI P G, CATALANO L A

Evaluation of the flow forces on a direct (single stage) proportional valve by means of a computational fluid dynamic analysis

[J]. Energy Conversion and Management, 2007, 48 (3): 942- 953

DOI:10.1016/j.enconman.2006.08.024      [本文引用: 1]

FU G

An extension of Hertz's theory in contact mechanics

[J]. Journal of Applied Mechanics, 2007, 74 (2): 373- 374

DOI:10.1115/1.2188017      [本文引用: 1]

HUTCHINGS I, SHIPWAY P. Tribology: friction and wear of engineering materials [M]. London: Butterworth-Heinemann, 2017.

[本文引用: 1]

/