华中科技大学电液控制工程四通电液伺服阀控液压缸控制系统仿真摘要：本文通过对电液伺服四通滑阀控液压缸系统进行数学建模和Matlab仿真研究系统的传递函数、响应特性以及波特图。
关键词：四通伺服建模仿真响应特性波特图引言：电液伺服控制系统是电液控制技术最早出现的一种应用形式。
通常所说的电液伺服控制系统，从其构成来说，就是指以电液伺服阀作为电液转换和放大元件实现某种控制规律的系统，它的输出信号能跟随输入信号快速变化，所以有时也成为随动系统。
电液伺服控制系统将液压技术和电气、电子技术有机地结合起来，既有快速易调和高精度的响应能力，又有控制大惯量实现大功率输出的优势，因而在国防和国民经济建设的各个技术领域得到了广泛的应用。
作为电液伺服系统中不可缺少的组成部分，液压动力机构由液压控制元件、执行元件和负载组成，又称为液压动力元件，它的动态特性对大多数液压伺服系统的性能有着决定性影响，因此，其传递函数是分析整个液压伺服系统的基础。
液压动力元件可以分为四种基本形式：阀控液压缸、阀控液压马达、泵控液压缸和泵空液压马达。
四种液压动力元件虽然结构不同，但其特性是类似的，本文通过建立数学模型，分析零开口四通滑阀和对称液压缸组成的液压系统的流量特性、力平衡方程和控制传递函数，获得系统的响应特性。
系统组成和原理：电液伺服控制系统根据输出信号的不同分为电液位置伺服系统、电液速度伺服系统和电液力伺服系统。
本文四通阀控液压缸属于电液位置伺服系统，其原理如右图1所示，四通滑阀控制液压缸拖动带有弹性和粘性阻尼的负载作往复运动。
该液压伺服控制系统的结构框图则如下图2所示。
图1 四通阀控液压缸原理图u i + u g i q L F L建立系统数学模型：流量方程由图1可知，从阀进入液压缸做强的流量除了推动活塞运动外，还要补偿液体的压缩量和管道等的膨胀量，补偿液压缸内、外泄漏，即q 1=A p dx p /dt+V 1/βe (dp 1/dt)+Ci(p 1-p 2)+C e p 1 (1)q 2=A p dx p /dt-V 2/βe (dp 2/dt)+Ci(p 1-p 2)-C e p 2 (2)式中，A p 为活塞面积，x p 为活塞位移，分别为左右进油腔容积，为液压弹性模量，分别为液压缸左右腔压力。
又由于q L =0.5 (q 1+ q 2)可得q L = A p dx p /dt+V 1/4βe (dp 1/dt)+C t p L (3)式子表明：四通阀控液压缸的负载流量包括推动液压缸活塞运动所需的流量、总泄漏量、总压缩量。
力平衡方程忽略活塞与缸体的摩擦力和油液质量的影响，液压缸的输出力与负载力（运动件的惯性力、运动件的粘性摩擦力、弹性负载力和其他负载力）平衡，由图1可得A p p L =M t d 2x p /dt 2+B p dx p /dt+kx p +F L (4)式中，为活塞及由负载折算到活塞上的总质量，为活塞及负载等运动件的粘性阻尼系数，k 为负载运动时的弹簧刚度，为作用在活塞上的其他负载力。
传递函数及系统方框图将式（3）（4）进行拉普拉斯变换之后参照文献【1】可以得到传递函数式X p =X vK q A p −K ce A p 2(V t4βe K ces+1)F L V t M t4βe A p 2s3+(K ce M t A p 2+B p V t 4βe A p 2)s 2+(kV t 4βe A p 2+B p Kce A p 2+1)s+K ce k A p 2(5)式中，K ce =(K C +C t )是包括泄漏在内的总的压力流量系数，阀芯位移X v 是指令信号，负载力F L 是干扰信号。
此系统中，液压动力元件的负载主要是惯性负载，弹性刚度k 很小，即弹性负载很小，可以忽略；而负载粘性摩擦系数B p 一般较小，B p V t s 2/(4βe A p 2)、B p K ce s/A p 2这两项可以忽略；令液压弹簧刚度K h =4βe A p 2/V t ，它是液压缸两腔完全封闭时液体压缩形成的液压弹簧的刚度，K h >>k ，认为k/K h <<1，因此式(5)可以简化为放大器 伺服阀 液压缸负载 反馈装置 图2 系统原理方框图X p=X v K qA p−K ceA p2(V t4βe K ces+1)F Ls(s2ωh2+2ξhωhs+1)(6)本系统采用型号为QDY2-D10的电液伺服阀，阀口面积梯度W=2.5X10-2m，间隙δ=5X 10-6m，传递函数G sv(s)=3060X10-6/(s2/6002+2 X0.5s/600-1)；最大行程L=0.5m，折算总质量M t=1000kg；负载力F L=2500N；供油压力p s=63 X105Pa；反馈传感器增益K f=100V/m；放大器增益K a=0.048A/V；液压缸有效面积A p=8 X10-4m2；油液的弹性体积模量βe=7000 X105Pa；液压油的绝对粘度μ=1.8 X10-2m。
若只考虑中位泄漏，有K ce=πWδ2/32μ=3.41X10-12(m3/s)/PaV t=A p L=4 X10-4m3固有频率和阻尼系数ωh=√4βe A p 2V t M t =67s-1ξh=KceA p√βe M tV t=0.2负载力系数Kce A p2=0.53 X10-5m/sN Vt4βe Kce=0.0042由式(6)可得，只有指令信号X v时Xp X v =K qA ps(s2ωh2+2ξhωhs+1)(7)只有干扰信号F L时Xp F L =−K ceA p2(V t4βe K ces+1)s(s2ωh2+2ξhωhs+1)(8)综上所述可以得到系统的方框图如下F L(s)-u i(s)+ I(s) Q L(s) + X p(s)-u f(s) 0.0483060 X10−6s26002+2X0.5600s+11.25X1051s(s2672+2X0.2s67+1)0.53X10-5X(1+0.042s)100图3 系统方框图Matlab仿真：代码>> w=-8*pi:0.01:8*pi;I=tf([0.048],[0 1 ]); %各环节传递函数Q1=tf([3060*10^(-6)],[1/600^2 2*0.5/600 1]);Q2=tf([1.25*10^5]);Q=I*Q1*Q2;F=tf(0.53*10^(-5)*[0.042 1],[0 1]); %干扰力环节G=tf([0 1],[1/67^2 2*0.2/67 1 0]);X1=F*G;H1=tf([100],[0 1]);X2=feedback(Q*G,H1,-1);X=X2-X1; %闭环传递函数>> XsvL=Q*G-X1;>> Xsv=Q*G %无干扰开环传递函数Transfer function:18.36------------------------------------------------------------------6.188e-010 s^5 + 3.879e-007 s^4 + 0.0002355 s^3 + 0.007637 s^2 + s>> bode(Xsv)>> grid on>>bode(XsvL) %绘制无干扰开环波特图>> figure;>> XsvL %有干扰开环传递函数Transfer function:-1.377e-016 s^6 - 8.962e-014 s^5 - 5.448e-011 s^4 + 0.00409 s^3 + 0.1096 s^2 + 18.36 s-----------------------------------------------------------------------------------------------1.378e-013 s^8 + 9.01e-011 s^7 + 5.539e-008 s^6 + 3.495e-006 s^5 + 0.0005039 s^4 + 0.01361 s^3+ s^2>> bode(XsvL)>> grid on**比较图4和图5可以发现，干扰力主要影响系统的相频特性，而且影响的频段频率较高。
>> X %闭环传递函数Transfer function:-1.377e-016 s^6 - 8.962e-014 s^5 - 5.448e-011 s^4 + 0.00409 s^3 + 0.1096 s^2 + 18.36 s- 0.009731----------------------------------------------------------------------------------------------1.378e-013 s^8 + 9.01e-011 s^7 + 5.539e-008 s^6 + 3.495e-006 s^5 + 0.0005039 s^4 + 0.4226 s^3+ 11.96 s^2 + 1836 s>> figure;step(X)参考文献：【1】电液控制技术，易孟林曹树平刘银水，华中科技大学出版社【2】电液伺服比例阀控缸位置控制系统仿真研究，孙衍石靳宝全熊晓燕，太原理工大学机械电子研究所【3】matlab绘制系统函数波特图，，源于百度。