弹簧质量阻尼系统的建模与控制系统设计样本.docx
- 文档编号:17990678
- 上传时间:2023-04-24
- 格式:DOCX
- 页数:34
- 大小:476.66KB
弹簧质量阻尼系统的建模与控制系统设计样本.docx
《弹簧质量阻尼系统的建模与控制系统设计样本.docx》由会员分享,可在线阅读,更多相关《弹簧质量阻尼系统的建模与控制系统设计样本.docx(34页珍藏版)》请在冰豆网上搜索。
弹簧质量阻尼系统的建模与控制系统设计样本
任课教师签字
华北电力大学研究生结课作业
学年学期:
第—年第—学期
课程名称:
线性系统理论
学生姓名:
号:
提交时间:
.11.27
目录2
1研究背景及意义4
2弹簧-质量-阻尼模型4
2.1系统的建立5
2.1.1系统传递函数的计算7
2.2系统的能控能观性分析9
2.2.1系统能控性分析10
2.2.2系统能观性分析11
2.3系统的稳定性分析12
2.3.1反馈控制理论中的稳定性分析方法错误!
未定义书签。
2.3.2利用Matlab分析系统稳定性14
2.3.3Simulink仿真结果16
2.4系统的极点配置19
2.4.1状态反馈法19
242输出反馈法21
242系统极点配置21
2.5系统的状态观测器23
2.6利用离散的方法研究系统的特性26
2.6.1离散化定义和方法26
2.6.2零阶保持器29
2.6.3—阶保持器31
2.6.4双线性变换法34
3.总结36
4•参考文献37
弹簧•质量•阻尼系统的建模与控制系统设计
1硏究背景及意义
弹簧、阻尼器、质量块是组成机械系统的理想元件。
由它们组成的弹簧-质量-阻尼系统是最常见的机械振动系统,在生活中具有相当广泛的用途,缓冲器就是其中的一种。
缓冲装置是吸收和耗散过程产生能量的主要部件,其吸收耗散能量的能力大小直接关系到系统的安全与稳定。
缓冲器在生活中处处可见,例如我们的汽车减震装置和用来消耗碰撞能量的缓冲器,其缓冲系统的性能直接影响着汽车的稳定与驾驶员安全;另外,天宫一号在太空实现交会■对接时缓冲系统的稳定与否直接影响着父会对接的成功。
因此,对弹簧-质量-阻尼系统的研究有着非常深的现实意义。
2弹簧-质量-阻尼模型
数学模型是定量地描述系统的动态特性,揭示系统的结构、参数与动态特性之间关系的数学表示式。
其中,微分方程是基本的数学模型,不论是机械的、液压的、电气的或热力学的系统等都能够用微分方程来描述。
微分方程的解就是系统在输入作用下的输出响应。
因此,建立数学模型是研究系统、预测其动态响应的前提。
—般情况下,列写机械振动系统的微分方程都是应用力学中的牛顿定律、质量守恒定律等。
弹簧-质量-阻尼系统是最常见的机械振动系统。
机械系统如图
2.1所示,
图2-1弹簧-质量-阻尼系统机械结构简图
其中呵、吒表示小车的质量,q表示缓冲器的粘滞摩擦系数,広表示弹簧的弹性系数,表示小车所受的外力,是系统的输入即
比p=Fip,兀⑴表示小车的位移,是系统的输出,BDyiW=,
i=l,2。
设缓冲器的摩擦力与活塞的速度成正比,其中叫"kg,m2=2kg,k1=k3=100N/cm.,k2=300N/cm,cx=c3=3N・s/cm,c2=6N・s/emo
2.1系统的建立
由图2.1,根据牛顿第二定律,分别分析两个小车的受力情况,建立系统的动力学模型如下:
对叫有:
叫肴+c禹+c2(Xx-X2)+kpq+kjtXi-X2)=Fi(t)
对叫有:
m2X2+c3X2+c2(X2-心)+k3X2+k2(X2—XJ=F2(t)
联立得到:
对…:
叫肴+(c】+-c2X2+(k】+k2)Xt-k2X2=F/t)
对叫•m2^2+Cc2+C3)^2~CZ^1+(%+*3)X2、—k2X1=F2W
令=Xpx2=X2,xs=x4=X2,ut=FvU2=F2;
Yi=X[,y2=X2
得出状态空间表示式:
因此,状态空间表示式为:
0
0
10
-
0
0I
0
0
01
0
0
ki+k2
&
ci+ca
C2
、+
1
0
mL
m二
nv
X
m2
k2+kS
C2
C2+C3
0
1
m2
m2
rxij
m2-
由此能够得出
1000
■0100.
已知:
=1,m2=2,kt=k3=100,k2=300,cx=c3=3,c2=6
代入数据得:
00101
000
-4003009
00.5」
1000
.0100.
2丄1系统传递函数的计算
在Matlab中,函数ss2tf给出了状态空间模型所描述系统的传递函数,其一般形式是[num,den]=ss2tf(A,B,C,D,iu),其中iu是输入值。
用Matlab将状态空间表示式表示为传递函数:
在输入1单独作用的情况下
A=[0010;0001;-400300-96;150-2003-4.5];
B=[00;00;10;00.5];
C=[l000;0100];
D=[00;00];
[num,den]=ss2tf(A,B,C,D,1)
运行程序,得到:
num=
0
-0.0000
1.0000
4.5000
200.0000
0
-0.0000
-0.0000
3.0000
150.0000
den=
1.0e+004
*
0.0001
0.0014
0.0623
0.1800
3.5000
在输入2单独作用的情况下:
A=[0010;0001;-400300-96;150-2003-4.5];
B=[00;00;10;00.5];
C=[l000;0100];
D=[00;00];
[num,den]=ss2tf(A,B,C,D,2)
运行程序,得到:
num=
0-0.00000.50004.5000200.0000
den=
1.0e+004*
0.00010.00140.06230」8003.5000
由此可知:
位移心⑴对外力%)的传递函数是:
十4.5s+200
-s4+14s3+623s2+1800s4-35000
位移%)对外力%的传递函数是:
X23s+150
_s4+14s3+623s2+1800s+35000
位移X1(©对外力F2(◎的传递函数是:
Xi(Q3s+150
-s4+14s3+623s2+1800s+35000
位移X2©对外力F2@)的传递函数是:
X2(t)_0.5s2+4.5s+200
_s4+14s3+623s2+1800s+35000
2.2系统的能控能观性分析
在反馈控制理论中只讨论输入量对输出量的控制。
而这两个量的关系唯一地由系统的传递函数所确定。
一个稳定的系统,一定能控。
同时,系统的输出量本身就是我们想要控制的量,对于一个实际的系统来说,输出量当然杲能够被观测到的,因此在反馈控制理论中没有必要设立能控和能观这两个概念。
然而在现代控制理论中,能控和能观是两个重要的基本概念。
我们把反映系统内部运动状态的状态向量作为被控量,而且它们不一定杲实阪上可观测到的物理量,至于输出量则是状态向量的线性组合,这就产生了从输入量到状态量的能控性问题和从输出量到状态量的能观测性冋题。
在现代控制中,分析和设计一个控制系统,必须研究这个系统的能控性和能观性。
状态方程描述了输入u(t)引起状态X(t)的变化过程;输出方程则描述了由状态变化引起的输出Y(t)的变化。
能控性和能观性正是分别分析U⑴对状态X(t)的控制能力以及Y(t)对X(t)的反应能力。
2.2.1系统能控性分析
设线性定常系统的状态方程为
x yp=Cxp+dup 式中anxn矩阵 Bnxr^P车 C——mxn矩阵 Dmxr矩阵 系统能控的充分必要条件为: 能控判别阵U广[BAB...AHB]的秩R(Uc)=n, 用Matlab计算能控矩阵的秩,从而对该系统的能控性进行判 别,程序为: A=[0010;0001;-400300-96;150-2003-4.5]; B=[00;00;10;00.5]; C=[l000;0100]; D=[00;00]; Qc=ctrb(A,B) Rl=rank(Qc) 运行程序,得到: R1= 4 等于矩阵行数,由此能够判断,系统是完全能控的。 2.2.2系统能观性分析 设线性定常系统的状态方程为: Xp=AxCt)+Bu* 卩=Cxp+Du(廿 式中Anxn矩阵 Bnxr矩P车 C——mxn矩阵 Dmxi•矩阵 能观的充分必要条件为: 能观判别阵Uo=[CCA...CAn-1]的秩 R(Uo)=n, 下面,用Matlab计算能控矩阵的秩,从而对该系统的能控性进行判断: A=[0010;0001;-400300-96;150-2003-4.5]; B=[00;00;10;00.5]; c=[l000;0100]; D=[00;00]; Qo=obsv(A,C) R2=rank(Qo) 运行程序,得到: R2= 4 满秩,因此能够判断,该系统是完全能观的。 综上所述,这杲一个既能控又能观的系统。 2.3系统的稳定性分析 2.3.1反馈控制理论中的稳定性分析方法 稳定性杲一个系统能够被采用的最基本的条件,是系统的固有属性。 稳定系统的定义如下: 设控制系统处于某一起始的平衡状态,在外力的作用下,它离开了平衡状态,当外作用消失后,如果经过足够长的时间它能够恢复到起始的平衡状态,则称这样的系统为稳定的系统,否则称为不稳定的系统。 由稳定性的定义可见,稳定性是系统去掉外力作用后自身的一种恢复能力,因此是系统的一种固有特性。 对于线性定常系统,它取决于系统本身的结构和参数,而与初始条件和外界作用无关。 线性定常系统稳定的充分必要条件杲: 闭环系统特征方程的所有特征根为负实数或具有负实部的共馳复数,即所有特征根位于复平面的左半平廂。 只要有一个闭环特征根分布在右半平廁上,系统就是不稳定的;如果没有右半平⑥的根,但有纯虚根,则系统是临界稳定的;在工程上,处于不稳定和临界稳定的线性定常系统是不能采用的叫 在古典控制系统中,我们判断系统的稳定性经常见劳斯-赫尔维茨代数判据.时域分析法、根轨迹法、频域分析法等方法,但那只针对低阶系统。 实际的工业生产中,经常合遇见一些特别复杂的系统。 这时古典控制理论中的方法就有点捉襟见肘了。 1892年俄国学者李雅普诺夫提出的稳定性理论是确定系统稳定性的更一般性理论,它采用了状态向量描述,不但适用于单变量、线性、定常的系统,而且适用于多变量,非线性、时变的系统。 李雅普诺夫理论在建立一系列关于稳定性概念的基础上,提出了判断系统稳定性的两种方法: 一种方法是利用线性系统微分方程的解来判断系统稳定性,称为李雅普诺夫第一法或间接法;另一种方法杲首先利用经验和技巧来构造李雅普诺夫函数,进而利用李雅普诺夫函数来判断系统稳定性,称为李雅普诺夫第二法或直接法。 2.3.2利用Matlab分析系统稳定性 随着计算机技术的发展,在现代控制理论中,我们经常采用Matlab判断系统的稳定性。 对于线性定常系统,典型的系统输入信号类型有脉冲.阶跃、斜坡.加速度、正弦信号。 系统的稳定性是对任何输入信号而言,即若一个系统是稳定的,则其在任何输入信号情况下对应的输出曲线是收敛的。 然而,阶跃信号包含了另外几种常见输入信号的特性,因此我们常经过观察系统的单位阶跃响应曲线判断判断系统的稳定性。 若系统的单位阶跃响应是收敛的,则系统一般是收敛的;否则,是发散的。 在Matlab中输入相应系统的状态空间表示式矩阵来求取系统的特征值: A=[0010;0001;-400300-96;150-2003-4.5]; B=[00;00;10;00.5]; C=[l000;0100]; D=[00;00]; eig(A) 运行程序,得到: ans= -5.7735+22.3859i -5.7735-22.3859i -0.9765+8.0332i -0.9765・&0332i 由此能够知道,经计算得出A阵的所有特征根均在复平面的左半平庖,因此得出该系统是稳定的。 给系统加起阶跃信号: A=[0010;0001;-400300-96;150-2003-4.5]; B=[00;00;10;00.5]; C=[l000;0100]; D=[00;00]; step(A,B,C,D) 结果如下 图2-2阶跃响应曲线 由图能够看出,在阶跃响应下,系统在一定时间内收敛于某一固定值,因此能够判断系统是稳定的,但同时我们也能够看出,系统的调节时间比较长,如果想要减少调节时间,那么需要重新配置极点,对系统进行改进。 下面的章节将对系统进行极点的配置。 2.3.3Simulink仿真结果 根据上述原理,用Matlab中的Simulink组件进行仿真。 根据状态空间表示式,搭建系统模型如下图所示: 我们分别对只有输入1作用下和只有输入2作用下的系统使用 Simulink进行仿真,让其与Matlab图像进行对比 图2-4ul作用时响应曲线 图中,绿色为输出1的曲线,蓝色为输出2的曲线。 经分析: 此曲线与对应Matlab曲线一致,系统稳定,可是超调量较大,调节时间较长。 (2)仅有出作用,系统的输入如下所示: iroofriit0 图2-5u2作用时响应曲线 图中,绿色为输出1的曲线,蓝色为输出2的曲线。 经分析: 同样,此曲线与对应的Matlab曲线一致,系统稳定,可是超调量较大,调节时间较长。 在U]、112共同作用下,系统的输出如下图所示: 图2-6ul.u2同时作用时响应曲线 图中绿色为输出1的曲线,蓝色为输出2的曲线。 经分析: 此曲线与Matlab曲线一致,系统稳定,可杲超调量较大,调节时间较长。 需要逬行极点配置,使系统得到更好的性能。 2.4系统的极点配置 控制系统的性能主要取决于系统极点在根平面上的分布。 因此,在系统设计中,一般杲根据对系统的品质要求,规定闭环极点应有的分布情况。 所谓的极点配置就是,就是经过选择反馈矩阵K,将闭环系统的极点恰好配置在根平廂上所期望的位置,以获得所希望的动态性能。 2.4.1状态反馈法 极点冋题首先解决是否能经过状态反馈来实现给定的极点配置,即在什么条件下才有可能按照规定的要求来配置极点。 其次是,这样的反馈阵K如何确定的问题。 状态反馈示意图 图2-7状态反馈示意图 (1)采用状态反馈配置系统极点条件: 系统工(A,B,C)采用状态反馈,任意配置其闭环系统极点的充要条件为: 系统工SBC)完全能控。 若系统不是完全能控的,就必须按能控性分解,只能任意配置可控的极点。 (2)极点配置的方法: 若原系统S3)可控,则采用状态反馈阵K,有([A-BKM)可控。 设原系统的特征方程为*+%严+…+加+讣0。 其中K=[心虫,…则有: A-BK= 0 0 1… 0•… 0 0 B= 0 ■-仏)+心) -(©+«)•・・ -(陽+&)_ 丄 配置后的闭环特征方程为: S”+(%+«一1)/"+…+(5+人)$+(4)+&)=()• 假设闭环系统希望的极点为几=(人仏,…人),得到: /(/1)=($_入)心_/12)・($_人)=”'++・+爪+几 O 为使系统达到希望性能,对比式 (1)和式 (2)中系数,使之相等,即可求得状态反馈阵K=k,«,rS]。 采用状态反馈配置系统极 点不改变系统可控性,它不能影响系统中不可控部分模块。 2.4.2输出反馈法 图2-8输出反馈示意图 对于完全能控的单变量系统,不能采用输出线性反馈来实现闭环系统极点任意配置。 不能任意配置极点,正是输出线性反馈的基本弱点。 为了克服这个弱点,在经典控制理论中,往往采取引入附加校正网络,经过增加开环零极点的方法改变根轨迹走向,从而使其落在指定的期望位置上。 对于完全能控的单变量系统工S3C),经过带动态补偿器的输出反馈时限极点任意配置的充要条件是: 1.系统完全能观测;2. 动态补偿器的阶数为n-lo 2.4.2系统极点配置 在现代控制理论中是用系统内部的状态来描述系统的,因此经常从系统的状态引出信号作为反馈量。 利用状态反馈只能改变系统能控部分的极点,而不能改变不能控部分的极点,因此利用状态反馈进行极点配置的充分必要条件是系统必须是完全能控的。 对一个可控系统,在采用状态反馈后,能够实现闭环极点的任意配置,即经过状态反馈的方法,使闭环系统的极点位于任意期望 的位置上。 对于 xCt)=十Bug 其中X是状态变量(n维),u是控制信号,这里选取控制信号为u=—kx 因此,X&)=(A-BK)x(t) 系统的稳态响应和瞬态响应特性由矩阵A-EK的特征决定 虽然理论上系统的闭环极点离S左半平面越远越好,可是在工业生产实践中,系统极点离左半平⑥越远,系统的运动状态就变化的越快,这就要求执行机构快速运作,即使再好的执行元件也令短时间内被损坏掉。 因此新的极点的绝对值大约是原系统极点绝对值的3至4倍左右。 取Pl=-15+40i;P2=-15-40i;P3=-3+10i;P4=-3-10i; 利用Matlab进行极点配置,希望能够减小超调量,缩短稳定时间以优化系统。 Matlab程序如下: A=[0010;0001;-400300-96;150-2003-4.5]; B=[00;00;10;00.5]; C=[l000;0100]; D=[00;00]; p=[-l5+40i,-l5-40i,-3+10i,・3・10i]; k=place(A,B,p) step(A-B*k,B,C,D) 运行程序,得到: -234.6522 131.8512 14.4561 6.3957 643.3762 -89.9765 6.7658 36.0878 StepFtesponse x103 From: h (1) From: ln (2) Time(sec) 图2-9稳态响应曲线 由响应曲线能够看出该系统重新配置极点后,具有较快的调节时间,而且也减少了超调量,改进了系统的动态性能与稳态性能。 2.5系统的状态观测器 -tZ}— 图2-10状态观测器示意图 经过状态观测器能够任意配置系统的极点,从而使闭环系统具有期望的稳态和动态性能。 但在工业生产中,系统的状态变量并非都是物理量,或者杲难以测得的量。 这样一来,系统的所有状态变量未必都能够直接测量得到,因此,状态反馈这种控制方式在许多实际控制问题中往往难以直接应用和实现。 状态观测器就是利用系统的外部输入输出信息来确定系统内部的状态,进而,在系统的极点配置状态反馈中,用观测器得到的状态估计值代替系统的真实状态。 下图为状态观测器的结构图: 图2J1状态观测器示意图 使用MATLAB为本系统设置状态观测器,选用极点配置时的 极点,程序如下图所示: A=[0010;0001;-400300-96;150-2003-4.5]; B=[00;00;10;00.5]; C=[l000;0100]; D=[00;00]; p=[-l5+40i,-l5-40i,-3+10i,・3・10i]; Kl=place(A,B,p) A1=A-B*K1 Ll=(pbce(A',C,p))‘ A2=A丄1*C L2=(place(Ar,C,,p))t A3=A1丄2*C sys2=ss(A2,B,C,D) sys2=ss(A3,B,C,D) 运行上面程序,得到: LI= 7.0833 30.0895 -30.5796 15.4167 -41.6552 -96.5401 16&1877 200.0790 A2= -7.0833 -30.08951.0000 0 30.5796-15.4167 01.0000 -358.3448 -1&1877 L2= 396.5401 -400.0790 -9.0000 3.0000 6.0000 -4.5000 3.7432 71200 -21.4563 -3.7432 190.9894 93.5822 115.5037 -24.2083 A3二 -3.7432 7.1200 1.0000 0 21.4563 3.7432 0 1.0000 -655.5795 -119.9176 -18.2856 30.9515 -81.9216 -402.0612 -29.0527 -17.7144 其中L1代表没逬行状态反馈时的状态观测反馈矩阵,L2代表进行了状态反馈的状态观测矩阵。 2.6利用离散的方法研究系统的特性 2.6.1离散化定义和方法 利用数字计算机对线性定常连续系统求数值解是现代科学技术研究中常见的一种方法,它不但方便,而且精确。 由于实际工业生产中线性定常连续系统被控对象需要在线控制等,必须将连续时间系统的状态方程转化为离散系统的状态方程,即将矩阵微分方程化成矩阵差分方程,这就是连续系统的离散化。 根据离散系统的构成设备不同能够将离散系统分为采样控制系统和数字控制系统: "•采样控制系统: 控制系统的构成中选择了采样开关(或者含有幵关特性的设备)。 b.数字控制系统: 控制系统的控制器选择了专用数字计算机。 —般,把系统中的离散信号是脉冲序列形式的离散系统,称为采样控制系统或脉冲控制系统;而把数字序列形式的离散系统,称为数字控制系统或计算机控制系统。 采样控制系统: 采样控制系统是对来自传感器的连续信息在某些规定上的时间瞬时值上取值。 例如,控制器系统中的误差信号能够是断续连续的脉冲信号,而相邻两个脉冲之间的误差信息,系统并没有收到。 如果在有规律的间隔上,系统取得了离散信息,则这种采样称为周期采样;反之,如果信息之间的间隔是时变的,或随机的,则称为非周期采样,或随机采样。 在采样控制系统中,把连续信号转变为脉冲序列的过程称为采样过程,简称采样。 实现采样的装置称为采样器,或采样开关。 用丁表示采样周期,单位为s;f5=1/T,表示采样频率,单位为1/S;气=2irfs=2tt/T
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 弹簧 质量 阻尼 系统 建模 控制系统 设计 样本