哈工大现代控制理论实验报告

1、现代控制理论基础上机实验报告之二基于降维观测器的亚微米超精密车床振动控制院 系航天学院专业自动化姓名李蒙班号1004102指导老师王述一哈尔滨工业大学2013年6月5日、工程背景介绍超精密车床隔振系统的作用超精密机床是实现超精密加工的关键设备，而环境振动又是 影响超精密加工精度的重要因素。为了充分隔离基础振动对超精 密机床的影响，目前国内外均采用空气弹簧作为隔振元件，并取 得了一定的效果，但是这属于被动隔振，这类隔振系统的固有频 率一般在2Hz左右。这种被动隔振方法难以满足超精密加工对隔 振系统的要求。为了解决这个问题，有必要研究被动隔振和主动 隔振控制相结合的混合控制技术。隔振系统的物理描述

2、上图表示了亚微米超精密车床隔振控制系统的结构原理，其 中被动隔振元件为空气弹簧，主动隔振元件为采用状态反馈控制 策略的电磁作动器。机床质量 空气弹簧粘性阻尼系数机床位移s一空气弹簧刚度系数地基位移0 主动隔振系统作动器（不表示参数）上图表示一个单自由度振动系统，空气弹簧具有一般弹性支 承的低通滤波特性，其主要作用是隔离较高频率的基础振动，并 支承机床系统；主动隔振系统具有高通滤波特性，其主要作用是 有效地隔离较低频率的基础振动。主、被动隔振系统相结合可有 效地隔离整个频率范围内的振动。床身质量的运动方程为：ms + F + F = 0（1）f空气弹簧所产生的被动控制力；f作动器所产生的主动控制

3、力。假设空气弹簧内为绝热过程，则被动控制力可以表示为：F = cy + k y + p 1 -V /(V + A y)nA(2)p 0 rr r eeV r标准压力下的空气弹簧体积；y = s-s0相对位移(被控制量)；p r空气弹簧的参考压力；A r参考压力下单一弹簧的面积；A： = 4 Ar 参考压力下空气弹簧的总面积；n -绝热系数。电磁作动器的主动控制力与电枢电流、磁场的磁通量密度及 永久磁铁和电磁铁之间的间隙面积有关，这一关系具有强非线 性。由于系统工作在微振动状况，且在低于作动器截止频率的低频范围内，因此主动控制力可近似线性化地表示为：F = kI(3)k e 力-电流转换系数；I

4、 a电枢电流。其中，电枢电流匕满足微分方程：LI + RI + E(I ,y) = u(t)(4)L 控制回路电枢电感系数；R控制回路电枢电阻；E 控制回路反电动势；u 控制电压。基于降维观测器的振动控制系统设计针对亚微米超精密车床的振动控制系统，我们可以采用全状态反 馈法设计控制规律。但是在工程实践中，传感器一般只能测量基座和 床身的位移信号，不能测量它们的速度及加速度信号，所以后两个状 态变量不能获得，换句话说全状态反馈很难真正实现。为了解决这个问题，本实验设计一个降维(2维)状态观测器，用来 解决状态变量、 的估计问题，从而真正实现全状态反馈控制。二、实验目的通过本次上机实验，使同学们熟

5、练掌握：降维状态观测器的概念及设计原理；线性系统分离原理的内涵；进一步熟悉极点配置及状态反馈控制律的设计过程; MATLAB语言的应用。三、闭环系统性能指标要求所设计系统：闭环系统单位阶跃响应的超调量不大于5%；过渡过程时间不大于0.5秒（ = 0.02）。闭环系统渐近稳定降维观测器渐近稳定已知一个车床的实际参数如下：k0 = 1200N/m，k = 980N/A，m = 120kg，c = 0.2，R = 300 Q，L = 0.95H。五、建立控制系统的开环状态空间模型假定s为常数，将式y = s-s0两边求关于时间的二阶导数可得: y = s = -L（F + F ）=一（cy +-k

6、y p 1-V （V + A y）nA + k I ）（5）m0 rr r ee e a记为：(6)y = - J-(cy + k y + +k I )m0e a其中=pi-W /(V+Ay)nA。rr r ee对式(6)两边继续求导得：(7)结合车床质量运动方程及电磁作动器原理，由式(6)可得:(8)则:my +cy + ky + =0a(9)综合以上各式:-CO - . eco_my+cy +妇+ _人河+以+ *+冲,力预,)aCOCOLmy + (Lc + Rtti)y + (Lk + Rc)y + Rk y + L + R k E(I ,y)(= -k U t) 00e ae将非线性

7、项况+R kEQ y)视为干扰信号，略去不计，可得 e a即:线性化模型为：Lmy +(LC + Rm)y + Lk + Rc)y + Rk y = -k 00e考虑到以上微分方程，.可取状态弯量为：(10)可得系统开环状态方程为：T % 了 Rkx =-0 x、3 Lm 1Lk + Rc Lc + Rm kr-r nx1x2x3=0R-010Lk + Rc5L LmLm由此得开环系统的状态空间表达式为:y = h0r i x101Lc + RmLmx2 +x30kLmu(11)x0 02x3Lm 2 Lm 3 Lm代入车床如下实际参数k = 1200N/m，k = 980N/A，m = 12

8、0kg，c = 0.2，R = 300 Q， 0eL = 0.95H。一 010 一x 0 一1001x+023157.910.5315.8x3一8.6得开环系统的状态空间表达式为:x1.x2.x3x0.1x2x一 010 一-0 -A=001，B=0-3157.9- 10.5-315.8-8.6即:C= 10六、降维观测器方程的推导过程由开环系统状态空间表达式可以看出，状态变量气是可观测 的。由于系统能观，rank（C）=1，故其余的2个状态分量只需用2维的降维观测器进行重构即可。降维观测器的设计过程如下：设受控系统记为：0 = （A，B，C），由于对降维观测器的要求 是保证闭环系统渐近稳定

9、和降维观测器渐近稳定，故只需要保证观测器极点在根平面左侧即可。考虑到需要突出主导极点的作用，则所选极点应该距离虚轴足够远。此处，选两个极点为-300和-400取变换矩阵为T，其中1 0 0CT-1 =010CL 0J0 0 1（由于只需保证第一行为C且1|=0，为方便计算，此处取变换阵为单位阵），1 0 0T =0100 0 1一 010 一A = T-1AT =001_至全一-3157.9- 10.53-315.8-0 -BB = T-1B = B =0=-1-8.6BL 2C = CT = C =1 0 0，经过此变换，系统被按照状态能检测分为了两部分（可检测部分和不可检测部分)。降维观测

10、器期望特征方程为：/*(!) = (X + 300)(1 + 400)二处 + 700人 +1200000设& = ，得观测器特征多项式：此+ 10.53)2f (X) = det(V-(A-GA)二处 + (315.79 + 炉)人 + (315.79 + g 221211故315.79 + 7 = 7001315.797 + F + 10.53 = 12000012解得：G =384.21 一g2-1340.2则w = (A-GA)w+ -GIT+ (JT-OT)y +(B -GB )ux23-A-_w+ Gy_L2384.21-1340.2221222122111A-X 二23计拿.-A

11、-A-二 w+ Gy得：-384.21 -A-F-14896010 w =1329.7-315.8w+930950y +u 8.6则:A-X =yw + 384.21yw -1340.2y2再变换到, 得:TA A 一1 X Xy元 + 384.21_L_七-1340.2y七、基于降维观测器的状态反馈控制律设计-&兀根据性能指标。=e i-C2 x 100% 5%，解得 q0.69，所以 p0 46.360。根据性能指标七子 0.5，解得知n Z 8。n留出裕量，取& = 0.8，o =12.5，则：列=10， J1-&2 = 7.5。为 此得两共轭极点为七=-10 + 7.5j，52 = -

12、10-7.5j，取第三个极点为 七=-100。于是得出系统期望特征多项式为： f（人）=。+10 + 7.5j）（人 +10-7.5j）。+100）=膈 + 120X2 + 2156.25人 +15625 （12）设状态反馈控制律为：u =伙1k2 k3 x2 + rx3则闭环系统的状态空间表达式为:.X1.X2.X300-3157.9 - 8.6k110-10.53 - 8.6k201-315.8 - 8.6kX0-1X2XL 3则此时闭环系统的特征多项式为f *(X) = X3 + X2(315.8 + 8.6k3) + X (10.5 + 8.6 k2) + 319 + 8.6k 71（13）将式（12）与式（13）比较可得:W9.7 k = 249.5 k =-22.8l 3315.8 + 8.6k = 120-22.8- w + 384.21y w： -1340.2y+r-Ar- 384.21 一*-148960-0 一w=1329.7-315.8w+930950y+-8.6u10001x+0-10.5-315.8-8.60 x = Ax + Bu =0-3157.9u通过数字仿真，可以知道系统是稳定的，并且满足动态性能要求。心得体会：通过设计进一步理解掌握了状态空间法的设计步骤，充分将所学 的知识应用到