气动系统建模仿真设计Word格式文档下载.docx
- 文档编号:6633825
- 上传时间:2023-05-07
- 格式:DOCX
- 页数:27
- 大小:365.47KB
气动系统建模仿真设计Word格式文档下载.docx
《气动系统建模仿真设计Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《气动系统建模仿真设计Word格式文档下载.docx(27页珍藏版)》请在冰点文库上搜索。
这使流动的音速和亚音速分界点不是压力比为0.528的点。
为解决这个问题,流量计算的新的发展趋势是通过临界压力比b来计算描述气动元件的过流能力,并用这个参数来计算经过比例阀的流量。
因此,比例阀进出气口的流量方程为:
(1-1)
(1-2)
Cd—流量系数
ω—阀口面积梯度
xv—阀芯位移
Ps、P0—气源压力、大气压力
P1、P2—气缸左、右腔压力
利用Simulink工具对进气口式进行建模,如图1-1所示,对其子系统封装如图1-2所示。
图1-1阀进口流量方程
图1-2阀进口流量方程封装模块
同理可对出气口进行建模并封装子系统。
2)压力微分方程
根据质量守恒定律,假定工作介质为连续的,储藏到某控制体中去的质量的储藏率应该等于流入的质量流量减去流出的质量流量。
即:
将气体状态方程代入上式并化简可得:
假定T1=T2=T,忽略温度变化的影响,将气缸两腔参数分别代入上面公式,得:
(1-3)
(1-4)
A1、A2—气压缸左、右腔面积
V1、V2—气缸左、右两腔体积
Qm1、Qm2—气缸进出左、右两腔的流量
x—气缸活塞位移
用Simulink对(3)式建模如图1-3所示,子系统封装如图1-5所示。
同理对(4)式进行建模如图1-4所示,子系统封装如图1-6。
图1-3气缸左腔流量压力方程
图1-4气缸右腔流量压力方程
3)气缸活塞的力平衡方程
根据牛顿第二定律可得到气缸的力平衡方程如下:
P1A1-P2A2-Ff=m
+Ky+F(1-5)
Ff—作用在气缸上的摩擦力
F—作用在气缸上的的外力负载
m—气缸上运动部件的质量及负载质量总和
K—负载弹簧刚度
根据力平衡方程(5)式在Simulink中建立模型如图1-7所示,进行子模型封装如图1-8所示。
图1-7气缸力平衡方程
图1-8气缸力平衡方程封装模块
4)气缸摩擦力模型
摩擦力是影响气动伺服控制系统性能的重要因素,摩擦力的大小、方向取决于滑动摩擦副的材料、表面粗糙度、润滑条件、受力大小及温度等因素。
气缸的摩擦力对气动伺服系统的影响最大,特别是气缸低速运动时更为明显,所以研究摩擦力的影响因素对系统的建模至关重要。
气缸摩擦力是非线性的,通常将气缸摩擦力分为动摩擦力和静摩擦力,其中动摩擦力又分为库伦摩擦力和粘性摩擦力。
当气缸在静止时摩擦力较大,而它一旦开始运动时,摩擦力随着速度增加急剧下降,在达到一定速度,即临界速度后又随着速度的上升而增加。
这一摩擦特性产生了气缸在低速运动时的爬行现象,同时影响气动伺服定位系统的性能。
当前工程上位置控制系统中应用较为广泛的气缸摩擦力模型是Stribeck摩擦模型,其摩擦力与速度关系曲线如图1-9所示,摩擦力首先随着速度的增加而降低,到一定速度后又随着速度的上升而下降,其公式为:
Fs—静摩擦力
Fc—库仑摩擦力
u—活塞速度
us—Stribeck分离速度
δ—待定系数,介于0.5到2之间
图1-9气缸Stribeck模型摩擦力与速度关系曲线
Stribeck摩擦模型较好地描述了低速下的摩擦力的行为,用一个衰减指数项体现了负斜率摩擦现象。
但是Stribeck模型没有考虑到摩擦滞后、变化的临界摩擦力等非线性因素带来的影响,在速度穿越零时,摩擦力发生突变,并且突变值较大,在力控制系统中直接反馈到的变量是力,摩擦力的突变会导致反馈力发生突变,进而引发系统高频振荡,不符合实际情况。
实际情况中,摩擦力还具有时间依赖性,即摩擦记忆的特性。
摩擦记忆就是接触表面间相对运动速度发生改变时,摩擦力滞后一段时间才会发生变化的现象,而LuGue模型较好的考虑了这一方面的因素,加入了摩擦力的记忆特性,避免了因为摩擦力突变而产生的高频振荡现象。
因此本仿真模型中采用LuGue模型,LuGue模型不仅考虑了Stribeck速度负斜率影响,并且能反映预滑动位移、摩擦滞后、变化的临界摩擦力和粘性滑动等非线性特性,是目前较为完善的一个模型,具有较高的精度。
LuGre模型将摩擦的接触面看成是在微观下具有随机行为的弹性鬃毛,摩擦力由鬃毛的挠曲产生,其摩擦力模型为:
(1-6)
(1-7)
(1-8)
v—摩擦表面的相对速度
Z—粘滞状态下相对运动表面间的相对变形量
a0—移动前的微观变形量z的刚度
a1—dz/dt的动态阻尼
a2—粘性摩擦系数
根据(1-6)、(1-7)(1-8)三个方程表述的摩擦力模型在Simulink中建模如图1-10所示,然后进行子系统封装。
图1-10气缸LuGre模型摩擦力方程
由LuGre模型作出气缸在低速时的摩擦力与速度的关系如图1-11所示。
此模型中的摩擦力具有记忆特性,在速度过零点时不会发生突变,而是有一定的滞后,在速度增加到反方向的某一个值时才缓慢的回到零,不会产生高频振荡。
并且摩擦力随速度变化关系也满足Stribeck负效应,符合摩擦力变化趋势,比较适合应用于气压力控制系统仿真模型中。
图1-11气缸LuGre模型摩擦力与速度关系曲线
上面已经对气压力控制系统的4个方程进行了建模,将4个子模型联系起来就可以完成对整个系统的建模。
●气压力控制系统的线性化
气压力伺服系统为比较复杂的非线性系统,特性也比较复杂,对其进行控制会比较困难,因此对其进行线性化,虽然线性化不能准确的给出实际系统模型,但它对系统的定性分析提供了一种有效的手段。
下面针对系统的数学模型在某一工作点对系统进行线性化处理。
首先对阀的流量方程(1-1)式(1-2)式进行线性化处理,由Sanville流量公式知,阀的流量方程仅是阀芯位移xv和气缸中气体压力P1和P2的函数,在工作点分别对这些变量取一阶偏微分即可得出微分方程的近似线性化方程:
(2-1)
(2-2)
式中:
然后对压力微分方程进行线性化处理,对(1-3)式(1-4)式进行拉氏变换得出:
从而
(2-3)
同理
(2-4)
其中
,
气缸的力平衡方程:
A1P1-A2P2-F-Ff=m
+Ky
在摩擦力模型中,有一部分与速度成正比的粘性摩擦力,因此线性化过程中可将摩擦力模型简化为
,则力平衡方程变为:
A1P1-A2P2-F-Fj=m
+Bp
+Ky
进行拉氏变换,得
(2-5)
将式(2-1),(2-2),(2-3),(2-4)代入式(2-5),得
可此求得由阀芯位移到气缸活塞位移的传递函数为:
=
在力控制系统中,被控制量是力,将输出力由力传感器转换为反馈电压信号与指令电压信号相比较,得到偏差信号,此偏差信号经过控制器输入伺服阀,使伺服阀到气压缸的流量发生变化,从而使输出力向着减小误差的方向变化。
在力控制系统中,输出力Fg为:
Fg=P1A1-P2A2-Ff=m
+Ky+F
将上式进行拉氏变换,得
又已知电压到阀芯位移的传递函数为二阶振荡环节,即
ωv—伺服阀固有频率
ζ—伺服阀阻尼比
K0—伺服阀增益
综合各部分的传递函数假设,系统的开环传递函数可由下式表示:
Kf为其他部分增益之积
必须指出,在以上分析中,特别是对一些关系式的推演过程,没有考虑气流的泄漏影响;
另外,还忽略了连接管道的分布阻力和管道柔度的影响,即我们采用的是集中参数模型,把管路阻力归并到控制滑阀口处,把弹性变形归并到气缸的活塞位移和气体的容积变化。
这种分析和分析液压伺服控制系统一样,也是在控制阀阀芯位移和气缸活塞位移变化在中间平衡位置附近的小扰动变化围进行的,即以阀的稳态特性的线性化为基础的。
在此引入气压弹簧的概念,假定一个理想的无摩擦无泄漏的气压缸,两个工作腔充满压力气体并被完全封闭。
由于气体具有可压缩性,当活塞受外力的作用时,活塞可以在气压缸移动,活塞的移动使气动缸的一腔压力升高,另一腔压力降低。
根据等熵的假定条件,体积弹性模数
与稳态时的腔工作压力
成正比,即
。
则有
则气压弹簧刚度Kh满足
,得
同液压弹簧一样,气压弹簧只有在动态时才有意义,在稳态时不存在。
假设气缸在初始位置处于平衡位置,即AP10=AP20,则
当活塞处在中间位置时,l1=l2=l/2,此时
上面的式子表明,气压弹簧刚度是活塞位置和工作点压力的函数,最低刚度出现在活塞行程的中间位置,此时气压固有频率最低。
当活塞偏离中间位置时,气压弹簧刚度增大,固有频率将增加。
由传递函数可知,气压系统与液压系统的传递函数具有相同的形式,其动态特征参数也很相似。
明显的差别就是可压缩工作介质,体积弹性模数
完全取决于稳态时的腔工作压力
和气体状态变化指数,即根据等熵的假定条件,
,因此,
的提高受到限制,初始工作压力过高,不仅带来安全问题,且系统元件密封液不易解决。
常规工业中使用的气体压力很低,因而气压伺服系统的固有频率和刚度都很低,和液压系统相比,响应速度慢,延滞时间长。
在系统设计时,应在工艺允许的条件下,尽量采用高的供气压力和尽可能短的连接管道,以提高伺服系统的输出刚度。
●系统仿真分析
在力控制系统开环传递函数中
=
其中Kf为其他部分增益之积
由此可见,系统传递函数由比例环节,二阶微分环节,积分环节和两个振荡环节共同组成的。
二阶微分环节和振荡环节的转折频率分别为负载固有频率ωm,系统固有频率ω0,以及伺服阀固有频率ωsv,并且ωm<
ω0。
下面分析下各个参数对系统传递函数的影响
a)
为负载固有频率,它随着负载弹簧刚度的增大而增大,随负载质量的增大而减小。
b)
为气压弹簧与负载弹簧并联耦合的刚度与负载质量形成的动力机构的固有频率。
它不仅与负载有关,还与气压弹簧刚度有关,气缸两腔面积越大,压强越大,气压弹簧刚度越大,并且气压弹簧刚度还受到活塞位置的影响。
c)
为动力机构的阻尼比。
粘性阻尼越大,负载质量越小,系统阻尼比越大。
负载弹簧刚度越大,气缸两腔面积越大,压强越大,系统阻尼比越小,系统阻尼比也受到活塞位置的影响。
d)
为系统增益。
负载弹簧刚度越大,伺服阀及控制器增益越大,系统增益越高。
气缸两腔压力、面积越大,系统增益越低。
系统增益也因活塞位置的不同而不同。
由上面分析知,系统的传递函数会随着活塞位置的变化而变化,所以我们在分析系统稳定性的时候,要选取系统最不容易稳定的点进行分析,使这一点稳定,系统才能稳定。
以下分具体情况进行讨论。
1)负载固有频率ωm大于伺服阀固有频率ωsv
系统的伯德图如图3-1所示,在伺服阀固有频率ωsv处斜率变为-60dB/10oct,通过负载固有频率时斜率变为-20dB/10oct,过了ω0时斜率又恢复为-60dB/10oct。
由于这种情况下ωm较大,负载弹簧刚度也一般很大,大于气压弹簧固有频率,因此ωm与ω0距离较近,且斜率一直为负值,因此ω0处的谐振峰值不会高于ωsv处幅值,因此谐振峰值不是导致系统不稳定的原因。
由伯德图可以看出,此时相角穿越频率略小于伺服阀固有频率ωsv,但是相角穿越频率处的幅值为正值,幅值裕度为负,系统不稳定,而系统增益是导致不稳定的原因。
此时穿越频率较大,快速性较好,而降低系统的穿越频率有利于系统的稳定性,同时快速性也能满足要求。
因此只需采用比例调节使幅值穿越频率降到小于相角穿越频率,使系统的幅值裕度和相角裕度为正值,系统稳定性较好,系统快速性受到的影响也不大。
随着系统各个参数变化,系统增益也发生变化,因此比例系数也要相应的发生变化。
校正后的系统伯德图如图3-2所示。
图3-1气压力伺服系统开环伯德图
图3-2比例校正后的气压力伺服系统开环伯德图
以一个仿真系统为例,负载分别为惯性负载和弹性负载,m=1kg,K=1000000N/m,控制系统的输出力为F=2000N,选取单出杆气压缸,缸径D=100mm,活塞杆直径为d=20mm,气缸行程为l=300mm,伺服阀的固有频率ωsv=500rad/s,阻尼比ζsv=0.5,此时ωm>
ωsv。
按上面的Simulink模型进行仿真,力响应曲线为图3-3,此时系统不稳定,对此进行比例控制,比例系数为Kp=0.01此时力响应曲线如图3-4所示,系统稳定。
图3-3力响应曲线图
图3-4比例校正后的力响应曲线图
下面分析各个参数在这种情况下对稳定性的影响。
a)质量负载m的影响
根据传递函数的公式知,m的大小影响负载固有频率,系统固有频率及阻尼比,但是对系统增益没有影响。
m的增大使负载固有频率和系统固有频率减小,使ωm向ωsv靠近,并且使系统的阻尼比减小,谐振峰值增加。
因此,在其他条件不变的情况下,增大m不利于系统的稳定。
但是m的增大如果在一定围,即负载固有频率不低于伺服阀固有频率,则系统可以通过比例调节达到稳定。
b)负载弹簧刚度K的影响
根据传递函数的公式知,K的大小影响负载固有频率、系统固有频率、阻尼比及其系统增益。
K的增大使负载固有频率、系统固有频率增加,并且距离靠近,影响可以近似抵消,使斜振频率远离伺服阀固有频率,但是系统阻尼比减小,由于谐振峰值不是影响稳定性的主要原因,对系统影响较小。
K值越大,系统增益越大,但是系统增益与K的关系并不是线性的,K值越大,增益变化越慢。
总体来说,K的增加对系统的影响是多方面的,在负载固有频率不低于伺服阀固有频率的围,总体影响较小。
c)气缸两腔压力及面积的影响
气缸两腔压力及面积影响系统固有频率、阻尼比及系统增益。
气缸两腔面积、压力越大,系统固有频率越大,阻尼比越小,系统增益越小。
在负载固有频率不低于伺服阀固有频率的前提下,负载弹簧一般较大,系统固有频率与负载固有频率距离较近,阻尼比的降低不会对稳定性造成太大的影响,而系统增益的降低幅度也很小,总体来说对系统稳定性影响不大。
d)活塞位置及行程的影响
活塞的位置影响系统固有频率,阻尼比及系统增益。
活塞行程越长,越靠中间,系统固有频率越小,阻尼比越大。
在负载固有频率不低于伺服阀固有频率的前提下,负载弹簧刚度一般较大,而引起的固有频率ω0及其增益部分的
变化较小,因此活塞行程越长,越靠近中间系统的增益
越小,总体来说有利于系统的稳定。
e)摩擦阻尼Bp的影响
阻尼影响系统的固有频率,Bp越大,系统的阻尼比越大,系统越稳定。
负载固有频率越大,他与伺服阀固有频率间距越大,谐振峰值对系统的影响较小,此时阻尼比的影响也较小。
但是负载固有频率与伺服阀频率较接近时,谐振峰值对系统稳定性影响稍大,需进一步降低比例系数,此时增大Bp,系统稳定性变好,能提高比例系数,提高系统的响应速度,得到较好的响应特性。
因此,在这种情况下,增益的变化对系统稳定性的影响较大,因此主要考虑参数变化对系统增益的影响。
2)负载固有频率小于伺服阀频率
此时系统的伯德图如图3-5所示,在负载固有频率ωm处斜率变为+20dB/10oct,通过动力机构固有频率时斜率变为-20dB/10oct,过了伺服阀固有频率ωsv时斜率变为-60dB/10oct。
如图所示,相角穿越频率略大于伺服阀固有频率,并且由于ω0与ωsv距离越远,相角穿越频率越靠近伺服阀固有频率ωsv,相角穿越频率幅值裕度仍为负值,因此考虑像上面一样采用比例控制将系统增益降低,使得幅值裕度和相角裕度为正,比例校正后的系统伯德图如图3-6所示。
校正后系统的幅值裕度和相角裕度都为正值,因此系统是稳定的。
但是如图所示,校正后系统的谐振峰值越过了零分贝线,这会使系统产生超调和震荡,并且谐振峰值越大,调整时间越长。
如果继续降低比例系数将谐振峰值降到零分贝线一下,会使穿越频率大大降低,系统响应过慢,因此只采用比例调节并不能达到较好的控制特性。
图3-5气压力伺服系统开环伯德图
图3-6比例校正后的气压力伺服系统开环伯德图
以上面的仿真系统为例,将负载改为质量m=100kg,刚度K=1000000N/m。
此时ωm<
ωsv,采用比例调节P=0.01,将此模型进行simulink模型仿真,得到的力响应曲线如图3-6所示。
图3-6力响应曲线图
由上面的力响应曲线图可知,系统是稳定的,但是开始的时候有超调和震荡,系统响应特性不好,因此只采用比例调节不能达到较好的控制效果。
下面分析各个参数在这种情况下影响。
m的增大使负载固有频率和系统固有频率减小,阻尼比减小,谐振峰值增加。
由于在这种情况下谐振峰值是导致系统响应特性差的主要因素,因此,m越大,系统的动态特性越不好。
K的增大使负载固有频率、系统固有频率增加,而气压弹簧刚度不变,因此负载固有频率与系统固有频率靠近,有利于系统的稳定性,但是会使系统阻尼比减小,对谐振峰值的影响有二者共同决定。
气缸两腔面积、压力越大,系统固有频率越大,阻尼比越小,系统增益越小,气压弹簧刚度越大。
减小系统增益有利于系统的稳定性,但是固有频率的增大和阻尼比的减小会使得谐振峰值增大。
活塞行程越长,越靠中间,系统固有频率越小,阻尼比越大,系统的增益越小,总体来说有利于系统的响应特性。
e)阻尼Bp的影响
Bp越大,系统的阻尼比越大,谐振峰值越小,系统的超调量越小,调节时间越短,增大Bp会使系统的响应特性变好。
因此,在这种情况下,谐振峰值对系统响应特性的影响较大,因此主要考虑参数变化对谐振峰值的影响。
为了使系统达到较好的特性,考虑像液压系统一样采用双惯性环节,在ωc与ωm之间加入转折频率为ω1的双惯性环节,因此必须先降低增益将ωm降到零分贝线一下,如图3-6所示。
当负载质量或刚度较大时,引起的谐振峰值也较大,采用较大转折频率的双惯性环节并不能完全将谐振峰值降到零分贝线以下,因此必须降低双惯性环节的转折频率。
考虑到阀的相角滞后,ω1处的相角低于-180o,过低的转折频率可能会导致系统本身不稳定,如图3-7所示。
为了使系统有一定的稳定裕度,必须将ω1处的幅值降到零分贝线以下,这就会导致系统的增益进一步降低。
并且从稳定裕度的角度考虑,ω1越小,系统增益越低,响应速度越慢。
所以,采用双惯性环节很难达到较好的控制结果,在仿真系统中尝试加入双惯性环节,将ω1从小到大进行调试,并随之调整增益,不能得到较好的控制特性。
仿真结果证明,在气压力控制系统中,双惯性环节是不可行的。
图3-7采用双惯性环节力控制系统伯德图
考虑采用PID控制,由于微分控制可以较好的改善系统的动态特性,能减小超调和振荡,使系统趋于稳定,缩短调节时间,减小稳态误差,允许增大比例系数,提高控制系统的精度。
但要注意的是,微分时间常数的选取很重要,过大或过小时,系统的响应时间仍然较长,超调量仍然较大。
只有微分时间常数选的比较合适的时候,才能得到较好的控制特性。
为了减小系统的超调和振荡,在比例控制的基础上加入微分控制,构成PD控制,加入微分控制后,系统产生小幅度的高频振荡,引进了高频干扰。
对于惯性较大的实际系统,标准数字PD控制器的积分项需要改进,并且,当瞬时偏差较大时,输出地控制量也很大,会很容易造成溢出的情况。
因此采用不完全微分数字控制器来解决上述问题。
在标准的数字PD控制器的计算式中,加入一个惯性环节,就构成了不完全微分数字控制器。
不完全微分数字控制器不仅可以平滑微分环节所产生的瞬时脉动,而且还能加强微分环节对全部控制过程的影响。
相对于标准数字PD控制,不完全微分数字PD控制器的控制品质较好。
理想微分容易引入高频干扰。
而不完全微分数字PD控制,由于微分作用能够缓慢的持续多个采样周期,使得一般的工业执行机构能够很好的跟踪微分作用输出。
并且由于不完全数字控制器计算式中含有一个一阶惯性环节,具有数字滤波的作用,因此,它的抗干扰能力也较强。
加入一阶惯性环节的不完全微分数字PD控制器的传递函数为:
由传递函数可知,不完全微分数字PD控制器的传递函数与滞后校正的传递函数一样,其中T是超前环节的转折频率,α为滞后超前比,它有滞后校正的优点,利用它的高频衰减特性,在保持系统稳定性的前提条件下,提高系统的低频增益,改善系统的稳态性能,或者在保证系统的稳态精度的条件下,降低系统高频部分的增益,来保证系统的稳定性。
气动力控制系统中,负载质量和刚度越大,阻尼比越小,提高放大系数的限制因素之一是增益裕量,而不是相位裕量,因此采用滞后校正是满足要求的。
而不完全微分数字PD控制器不仅具有滞后校正的上述优点,而且在减小超调和振荡上的效果更好。
在上面的系统中采用不完全微分数字PD控制器,引入T=0.4的一阶惯性环节,取P=0.01,D=0.002,得到较好的动态特性,如图3-8所示。
图3-6采用不完全微分数字PD控制器的力响应曲线图
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 气动 系统 建模 仿真 设计