气动系统建模仿真设计

气压控制伺服系统的数学建模及仿真模型建立

关于气动伺服系统的数学建模，主要是通过分析系统的运动规律，运用一些己知的定理和定律，如热力学定律、能量守恒定律、牛顿第二定理等，通过一些合理而必要的假设和简化，推导出系统被控对象的基本状态方程，并将其在某一工作点附近线性化，从而获得的一个近似的数学模型。虽然这些模型不是很准确，但还是能够反映出气动伺服控制系统的受力和运动规律，并且借此可以分析出影响系统特性的主要因素，给系统的进一步分析和控制提供依据和指导。

另外，利用Simulink 工具包可以不受线性系统模型的限制，能够建立更加真实的非线性系统，同时其模型分析工具包括线性化和精简工具。因此，本文在数学模型的基础之上，利用Simulink 对所研究的气压力控制系统尝试建立一个非线性数学模型，并对该模型进行计算机仿真。

由于气动系统的非线性，如气体的压缩性较大，且在气缸的运动过程中容腔中气体的各参数和变量是实时变化的，所以对气动系统的精确建模是比较困难的。所以为了建立系统的模型，我们对控制系统作一些合理的假设，来简化系统的数学模型。假设如下：

(1)气动系统中的工作介质—空气为理想气体； (2)忽略气缸与外界和气缸两腔之间的空气泄漏； (3)气动系统中的空气流动状态为等熵绝热过程； (4)气源压力和大气压力恒定；

(5)同一容腔中的气体温度和压力处处相等。 1) 比例阀的流量方程

在实际的伺服控制系统中气体的流动过程十分复杂，气动元件研究中使用理想气体等熵通过喷管的流动过程来近似代替。一般计算阀口的流量时采用Sanville 流量公式：

k

k s d k s

d

s

m P P P P k RT k P q 1212A +⎪⎪⎭

⎫ ⎝⎛-⎪⎪⎭⎫

⎝⎛-= 0.528<

s

d

P P ≤1 )1(212A 1

1

+⎪

⎭

⎫

⎝⎛+=-k RT k k P q k s m 0≤s

d P P

≤0.528

其中：P s—为阀口上游压力；

P

d

—为阀口下游压力。

0.528为临界压力比。当阀口上、下游的压力比小于等于0.528时，气体通过阀口的流量达到最大值，即气体以音速流动，此时下游压力的降低不会使质量流量再增加，出现了所谓的“壅塞”现象，这种现象使气体流经阀时具有很强的非线性，也是以空气作为传动介质系统中的固有特征。当控制阀上、下游压力之比大于0.528小于1时，通过阀口的气体质量流量不仅取决于阀的结构，而且还取决于阀口上、下游压力，此时通过阀口的气体流动状态为亚音速流动[26]。

由于气动元件部的结构比较复杂，不同于渐缩喷管。这使流动的音速和亚音速分界点不是压力比为0.528的点。为解决这个问题，流量计算的新的发展趋势是通过临界压力比b来计算描述气动元件的过流能力，并用这个参数来计算经过比例阀的流量。

因此，比例阀进出气口的流量方程为：

1

1

1

1

1

1

Q

2

1

d v

s

m

k

d v s

s

P

c x P b

P

P

c x P b

k P

ω

ω-

⎧

⎪

≤≤

⎪⎪

=⎨

⎪

⎛⎫

⎪≤≤

⎪

⎪+

⎝⎭

⎩

(1-1)

2

2

1

1

2

2

1

Q

2

1

d v

m

k

d v

P

c x P b

P

P

c x P b

k P

ω

ω-

⎧

⎪

≤≤

⎪⎪

=⎨

⎪

⎛⎫

⎪≤≤

⎪

⎪+

⎝⎭

⎩

(1-2)

其中：

Cd—流量系数

ω—阀口面积梯度

x

v

—阀芯位移

Ps、P

—气源压力、大气压力

P

1

、P

2

—气缸左、右腔压力

利用Simulink工具对进气口式进行建模，如图1-1所示，对其子系统封装

如图1-2所示。

同理可对出气口进行建模并封装子系统。 2) 压力微分方程

根据质量守恒定律，假定工作介质为连续的，储藏到某控制体中去的质量的储藏率应该等于流入的质量流量减去流出的质量流量。即：

()dM d V dV d M M V dt dt dt dt

ρρ

ρ••

-=

==+∑∑入出 将气体状态方程代入上式并化简可得：

1dM dV V dp p dt RT dt k dt ⎛⎫

=+ ⎪⎝⎭

假定T 1=T 2=T ，忽略温度变化的影响，将气缸两腔参数分别代入上面公式，得：

图1-1 阀进口流量方程

图1-2 阀进口流量方程封装模块

111

1

1

()

m dx

k RTQ P A dP dt dt

V -=

(1-3)

222

2

2

()

m dx k RTQ P A dP dt dt

V +=

(1-4)

其中：

A 1、A 2—气压缸左、右腔面积 V 1、V 2—气缸左、右两腔体积

Qm 1、 Qm 2—气缸进出左、右两腔的流量 x —气缸活塞位移

用Simulink 对(3)式建模如图1-3所示，子系统封装如图1-5所示。同理对(4)式进行建模如图1-4所示，子系统封装如图1-6。

图1-3 气缸左腔流量压力方程