现代控制理论大作业

1、Harbin Institute of Technology现代控制理论基础上机实验报告之一亚微米超精密车床振动控制系统的状态空间法设计课程名称： 现代控制理论 院 系： 航天学院 自动化 班 号： 作 者： 皮永江 学 号： 指导教师： 刘杨、井后华 哈尔滨工业大学2014年6月5日1. 工程背景介绍超精密机床是实现超精密加工的关键设备，而环境振动又是影响超精密加工精度的重要因素。为了充分隔离基础振动对超精密机床的影响，目前国内外均采用空气弹簧作为隔振元件，并取得了一定的效果，但是这属于被动隔振，这类隔振系统的固有频率一般在2Hz左右。2. 实验目的通过本次上机实验，使同学们熟练掌握：a)

2、控制系统机理建模b) 时域性能指标与极点配置的关系c) 状态反馈控制律设计d) MATLAB语言的应用3. 给定的实际参数与数学建模3.0参数与物理模型机床的已知参数k0=1200Nmm=120kgke=980NAc=0.2R=300L=0.95H上图表示了亚微米超精密车床隔振控制系统的结构原理，其中被动隔振元件为空气弹簧，主动隔振元件为采用状态反馈控制策略的电磁作动器。床身质量的运动方程为：ms+Fp+Fa=0 (1)Fp-空气弹簧所产生的被动控制力Fp-作动器所产生的主动控制力假设空气弹簧内为绝热过程，则被动控制力可以表示为：Fp=cy+k0y+pr1-Vr(Vr+Aey)nAe (2)V

3、r-标准压力下的空气弹簧体积y=s-s0-相对位移（被控制量）pr-空气弹簧的参考压力Ar-参考压力下单一弹簧的面积Ae=4Ar-参考压力下空气弹簧的总面积n-绝热系数电磁作动器的主动控制力与电枢电流、磁场的磁通量密度及永久磁铁和电磁铁之间的间隙面积有关，这一关系具有强非线性。由于系统工作在微振动状况，且在低于作动器截止频率的低频范围内，因此主动控制力可近似线性化地表示为：Fa=keIa (3)ke-力-电流转换系数Ia-电枢电流其中，电枢电流Ia满足微分方程：LIa+RIa+EIa,y=ut (4)L-控制回路电枢电感系数R-控制回路电枢电阻E-控制回路反电动势u-控制电压综上得到如下方程组

4、：ms+Fp+Fa=0 (1)Fp=cy+k0y+pr1-VrVr+AeynAe (2)Fa=keIa (3)LIa+RIa+EIa,y=ut (4)3.1如果忽略非线性部分数学建模ms+Fp+Fa=0 (1)Fp=cy+k0y (2)Fa=keIa (3)LIa+RIa=ut (4)y=smy=-(cy+k0y+keIa)my=-(cy+k0y+keIa) 整理 Ia=-1ke(my+cy+k0y)Ia=-1ke(my+cy+k0y)带入4式-Lkemy+cy+k0y-Rkemy+cy+k0y= ut整理得-Lmkey-Rm+Lckey-Rc+Lk0key-Rk0key= ut设状态变量为：

5、x1=y,x2=y,x3=y得到状态方程：x1=x2x2=x3x3=-Rk0Lmx1-Rc+Lk0Lmx2-Rm+LcLmx3-keLmu状态空间表达式：x1x2x3=-Rk0Lm-Rc+Lk0Lm-Rm+LcLm+00-keLmuy=100x1x2x3代入数据：-Rk0Lm=-300*12000.95*120=-3157-Rc+Lk0Lm=-300*0.2+0.95*12000.95*120=-10.52-Rm+LcLm=-300*120+0.95*0.20.95*120=-316-keLm=-9800.95*120=-8.6那么状态空间表达式为：x1x2x3=-3157-10.52-316

6、+00-8.2uy=100x1x2x3显然系统能控，可以采用状态反馈进行任意配置极点。3.2考虑非线性部分数学建模因为系统工作在低速，微位移情况下，那么对于(2)式中pr1-Vr(Vr+Aey)nAe0点泰勒一阶展开nprAe2Vrnyn=1.41Pr=0.4*105Pad=0.3mAr=d34= 0.0707Ae=4Ar=0.2827m2h=0.28mVr=Arh=0.0198m3q=nprAe2Vrn=1.41*4*10000*0.01981.41=1.14107(4)式中EIa,y=5.4y ，ms+Fp+Fa=0 (1)Fp=cy+(k0+q)y (2)Fa=keIa (3)LIa+R

7、Ia+5.4y=ut (4)y=smy=-(cy+k0y+keIa)my=-(cy+k0y+keIa) 整理 Ia=-1kemy+cy+k0+qyIa=-1kemy+cy+(k0+q)y带入4式-Lkemy+cy+k0y+qy-Rkemy+cy+k0y+qy+5.4y= ut整理得-Lmkey-Rm+Lckey-Rc+Lk0+Lq-5.4kekey-Rk0+Lqkey= ut设状态变量为：x1=y,x2=y,x3=y得到状态方程：x1=x2x2=x3x3=-R(k0+q)Lmx1-Rc+L(k0+q)-5.4keLmx2-Rm+LcLmx3-keLmu状态空间表达式：x1x2x3=-R(k0+

8、q)Lm-Rc+Lk0+Lq-5.4keLm-Rm+LcLm+00-keLmuy=100x1x2x3代入数据：-Rk0+RqLm=-300*1200+300*1.141070.95*120=-3107-Rc+Lk0+Lq-5.4keLm=-300*0.2+0.95*1200+1.14107-5.4*9800.95*120=-1.2108-Rm+LcLm=-300*120+0.95*0.20.95*120=-316-keLm=-9800.95*120=-8.6x1x2x3=-3107-1.2108-316+00-8.6uy=100x1x2x3显然系统能控，可以采用状态反馈进行任意配置极点。4.

9、性能指标与理论设计4.0性能指标：闭环系统单位阶跃响应的a) 超调量不大于5%b) 过渡过程时间不大于0.5秒（=2%）。由经验公式得到：pe-1-2100%0.7ts4n0.5(=2%) n8为留有一定余量，设n=10，=22(=45。)得到闭环主导极点为：-10j10,取第三个极点为-100；得到闭环特征多项式为：f*=+100+10+j10+10-j10 =3+1202+2200+200004.1忽略非线性状态控制采用状态反馈控制，设状态反馈矩阵为：K=k1k2k3u=v+KX=v+k1k2k3 x1x2x3f=detI-A+bK =000000-3157+8.2k1-10.52+8.2

10、k2-316+8.2k3=-100-13157+8.2k110.52+8.2k2+316+8.2k3=3+316+8.2k32+10.52+8.2k2+3157+8.2k1f=3+316+8.2k32+10.52+8.2k2+3157+8.2k1f*=3+1202+2200+20000316+8.2k3=12010.52+8.2k2=+8.2k1=20000 k1=2054k2=267k3=-24得出状态控制规律为：v=u+KX=u+-24 x1x2x3系统校正框图如下图：2054267-24yx1x2x3-8.21s1s1su-316-10.52-3157Figure 1系统模拟结构图4.2

11、考虑非线性状态控制K=k1k2k3u=v+KX=v+k1k2k3 x1x2x3f=detI-A+bK =000000-3157+8.2k1-10.52+8.2k2-316+8.2k3x1x2x3=-3107-1.2108-316+00-8.6u=-100-13107+8.2k11.2108+8.2k2+316+8.2k3=3+316+8.2k32+1.2108+8.2k2+3107+8.2k1f=3+316+8.2k32+1.2108+8.2k2+3107+8.2k1f*=3+1202+2200+20000316+8.2k3=1201.2108+8.2k2=22003107+8.2k1=200

12、00 k1=-3.65106k2=-1.46107k3=-24得出状态控制律为：u=v+KX=v+-3.65106-1.46107-24 x1x2x3控制流程简图在此处键入公式。uyVBCAKFigure 2系统模拟结构图5. 数字仿真5.1闭环传递函数校正后系统状态空间表达式：x1x2x3=-20000-2200-120x1x2x3+00-8.6uy=100x1x2x3Simulink仿真框图：Figure 3 Simulink仿真框图系统0初始状态响应曲线：Figure 4闭环系统0初始状态响应曲线系统非0初始状态响应曲线：Figure 5系统非0初始状态响应曲线：5.2闭环传递函数与数字

13、仿真s=C(sI-A)-1Bs=100s-100s-s+120-100-8.6=-8.6s3+120s2+2200s+20000为了更好的观测系统的超调量p，以及调整时间ts，可以将系统闭环传递函数增益调整为1，或者将阶跃响应的幅值调到20000-8.6，那么系统的稳态输出为1。这里采用前者：s=-8.6s3+120s2+2200s+20000s=20000s3+120s2+2200s+20000Matlab仿真程序为()：y=tf(20000,1,120,2200,20000);step(y)；grid;Figure 6系统阶跃响应曲线由图像可知，系统超调p4%，ts0.42s满足性能指标要求。6. 心得与体会本实验借助Matlab实现基于状态空间法设计的超精密车床振动控制。所设计系统超调p4%，ts0.42s（2%误差带）。设计满足指标要求。在设计过程中，先将系统的非线性部分利用泰勒一阶展开转换成线性方程组，建立状态空间表达式，原系统是完全能控的，那么能采用状态反馈进行