- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
哈工大现代控制理论实验报告一
Harbin Institute of Technology现代控制理论基础上机实验一亚微米超精密车床振动控制系统的状态空间法设计院系: 航天学院控制科学与工程系 专业: 探测制导与控制技术 姓 名: 班 号: 学 号: 指导教师: 史小平 哈尔滨工业大学 2015年5月26日目录一、工程背景介绍及物理描述31.1 工程背景介绍31.2实验目的31.3工程背景的物理描述3二.闭环系统的性能指标5三.实际给定参数6四.车床振动系统的开环状态空间模型6五.状态反馈控制律的设计过程7六. 闭环系统数字仿真的MATLAB编程86.1源程序86.2 运行截图9七. 实验结论及心得107.1实验结论107.2 心得体会11一、工程背景介绍及物理描述1.1 工程背景介绍超精密机床是实现超精密加工的关键设备,而环境振动又是影响超精密加工精度的重要因素。为了充分隔离基础振动对超精密机床的影响,目前国内外均采用空气弹簧作为隔振元件,并取得了一定的效果,但是这属于被动隔振,这类隔振系统的固有频率一般在2Hz左右。这种被动隔振方法难以满足超精密加工对隔振系统的要求。为了解决这个问题,有必要研究被动隔振和主动隔振控制相结合的混合控制技术。其中,主动隔振控制系统采用状态空间法设计,这就是本次上机实验的工程背景。1.2实验目的通过本次上机实验,熟练掌握:1. 控制系统机理建模;2. 时域性能指标与极点配置的关系;3. 状态反馈控制律设计;4. MATLAB语言的应用。四个知识点。1.3工程背景的物理描述上图表示了亚微米超精密车床隔振控制系统的结构原理,其中被动隔振元件为空气弹簧,主动隔振元件为采用状态反馈控制策略的电磁作动器。 上图表示一个单自由度振动系统,空气弹簧具有一般弹性支承的低通滤波特性,其主要作用是隔离较高频率的基础振动,并支承机床系统;主动隔振系统具有高通滤波特性,其主要作用是有效地隔离较低频率的基础振动。主、被动隔振系统相结合可有效地隔离整个频率范围内的振动。 床身质量的运动方程为: ——空气弹簧所产生的被动控制力;——作动器所产生的主动控制力。假设空气弹簧内为绝热过程,则被动控制力可以表示为: ——标准压力下的空气弹簧体积;——相对位移(被控制量);——空气弹簧的参考压力;——参考压力下单一弹簧的面积;——参考压力下空气弹簧的总面积;——绝热系数。电磁作动器的主动控制力与电枢电流、磁场的磁通量密度及永久磁铁和电磁铁之间的间隙面积有关,这一关系具有强非线性。 由于系统工作在微振动状况,且在低于作动器截止频率的低频范围内,因此主动控制力可近似线性化地表示为: ——力-电流转换系数;——电枢电流。其中,电枢电流满足微分方程: ——控制回路电枢电感系数;——控制回路电枢电阻;——控制回路反电动势;——控制电压。二.闭环系统的性能指标 闭环系统单位阶跃响应的超调量不大于5%,过渡过程时间不大于0.5秒()。三.实际给定参数=1200N/m m=120kg=980N/A c=0.2 R=300 L=0.95H四.车床振动系统的开环状态空间模型首先假定为常数,将式两边求关于时间的二阶导数可得: (5)记为 (6)其中。式(6)两边求导数得: (7)由式(6)可得: (8)由式(7)可得: (9)将式(8)和(9)代入式(4)可得:将非线性项视为干扰信号,略去不计,可得线性化模型: (10)五.状态反馈控制律的设计过程(1)根据带入实际参数后的线性化模型列写系统状态空间表达式:令则uy =x(2)设计一状态反馈控制律,使对y的控制满足性能指标由性能指标要求的闭环系统单位阶跃响应的超调量不大于5%,过渡过程时间不大于0.5秒(),取阻尼比为0.8,,进一步的,将闭环极点配置如下:设状态反馈控制律,则解得故加入控制律后的状态方程为uy =x六. 闭环系统数字仿真的MATLAB编程给定初始条件:,,下面进行仿真分析。6.1源程序第一个文件simu046.mfunction dx=simu046(t,x)A=[0,1,0; 0,0,1; -5850,-940,-63];dx=A*x;第二个文件do_simu046.m[t,x]=ode45(simu046,[0,1],[6*10^-5,2*10^-5,-0.8*10^-5]);subplot(3,1,1);plot(t,x(:,1),r-);legend(x_1);grid;subplot(3,1,2);plot(t,x(:,2),b-);legend(x_2);grid;subplot(3,1,3);plot(t,x(:,3),k-);legend(x_3);grid;阶跃响应曲线绘制:A=[0,1,0; 0,0,1; -5850,-940,-63];B=[0;0;-
文档评论(0)