阀控马达控制系统仿真

1.应用背景

在海上风机吊装工程中，经常需要进行速度控制，如原动机调速、吊臂垂直及回转装置的速度控制等。在该项目的工程设计及应用当中的电液位置伺服系统也经常采用速度局部反馈回路来提高系统的刚度和减小伺服阀等参数变化的影响，提高系统的精度。

电液速度控制系统按控制方式可分为：阀控液压马达速度控制系统和泵空液压马达速度控制系统。阀控液压马达系统一般用于小功率系统，而泵控马达系统一般用于大功率系统。在本次的实验中，主要针对阀控马达速动控制系统的校正前后变化，通过MATLAB的simulink对其进行仿真比较分析。

2.电液速度控制系统原理

首先给出阀控液压马达速度控制系统的实际物理模型：

如图１所示，该系统由伺服放大器、电液伺服阀、液压马达、测速电动机等组成。测速电机轴与负载机轴相联，用于检测负载轴的速度，检测到的速度信号与指令信号差（误差信号）经伺服放大器进行功率放大，产生的电流用来控制电液伺服阀的阀芯位置，电液伺服阀输出压力油驱动液压马达及负载旋转。

根据所建立的物理模型，可以建立相对应的闭环控制系统原理的方框图：

伺服放大器

伺服阀

液压马达

测速机

速度信号电压

+- 图2 阀控马达速度控制系统方框图

3. 系统各环节数学模型

3.1 伺服放大器

伺服放大器输出电流ΔI 与输入电压Ue 近似成正比，其传递函数可用伺服放大器增益Ka 表示：

（1）

但通常的速度控制系统采用积分放大器，对原系统加以校正才能稳定工

作。校正后的积分放大器增益Ka 表示为：

（2）

式中：Ue 为积分放大器额定电压，V ；

Ka 为积分放大器增益，A/V 。 3.2 伺服阀

伺服阀的流量增益为： (3)

式中：为伺服阀流量增益，m 3/(s*A)； 为伺服阀空载流量，m 3/s;

为伺服阀额定电流，A 。

伺服阀传递函数为：

(4)

式中：Q 0为伺服阀流量，m3/s；

ΔI为电流增量，A；

为伺服阀固有频率，rad/s;

为伺服阀的阻尼比。

3.3液压马达

负载总惯量为：

(5)

式中：为负载总惯量，kg*m2;

为液压马达惯量，kg*m2;

为负载的惯量，kg*m2。

因而求得液压的固有频率：

(6)

式中：为液压固有频率，rad/s；

为液压马达排量，m3/rad；

为液压油的有效体积弹性模量,Pa；

为总压缩体积，m3。

及阻尼比：

(7)

式中：为液压马达的阻尼比；

为阀的流量压力系数,m3/ (s·Pa)。

最终可以求得液压马达的传递函数：

(8)

式中：是液压马达负载的传递函数；

为液压马达的的增益，m3/ S·A。

3.4测速机

速度传感器（测速机）的数学模型为：

(9)

式中：为传感器电压，V；

为速度传感器增益，（Vs/rad）。

4.仿真系统动态方块图

根据所给出的系统原理方框图以及所求得的各环节数学模型，可以得出仿真系统经过积分环节校正后的系统动态方块图如下：

K a s

K sv

s 2

ωsv2

+2ξh

ωsv

s+1

K s

s2

ωh2+

2ξh

ωh s+1 K fs s

Ur

Uf Ue

+

-

θm

图3 积分环节校正后的速度控制系统方块图

5.问题描述

现有一阀控液压伺服速度控制系统，其原理图如图1，考虑伺服阀的动态，伺服参数为：=3060e-6,=600,=0.5,液压缸的参数为=1.25e6，=388，=0.94，放大器增益=0.05，试用simulink对该系统进行动态分析，并比较校正前后系统稳定性，分析校正前后的变化。

6.解题步骤

（1）根据给定的参数，利用simulink对该系统进行未加积分环节校正的

系统建模，如下图所示：

图4 未校正系统的simulink模型

给上述模型中的变量赋值，编写matlab程序，之后对上述闭环模型进行仿真，结果出错，得不出闭环仿真结果，说明闭环系统不稳定。

（2）绘制未校正系统的伯德图。Matlab程序如下：

clear all;clc;

Wsv=600;zuni1=0.5;Ksv=3060e-6;Ka=0.05;i=3;Kf=0.175;

Wh=388;zuni2=0.94;Kh=1.25e6;

sys=tf(Ka*Ksv*i*Kf*Kh,conv([1/Wsv^2 2*zuni1/Wsv 1],[1/Wh^2

2*zuni2/Wh 1]));

figure(1);margin(sys);

grid;

figure(2);subplot(121);

pzmap(sys);

grid on;

subplot(122);

nyquist(sys);

figure(3);subplot(121);

impulse(feedback(sys,Kf*i));

grid;

subplot(122);

step(feedback(sys,Kf*i));

运行上述程序后，得到未校正开环系统的奈氏图：