- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
亚微米超精密车床振动控制系统设计
现代控制理论基础
上机实验报告之一
亚微米超精密车床振动控制系统
的状态空间法设计
一:工程背景介绍
超精密机床是实现超精密加工的关键设备,而环境振动又是影响超精密加工精度的重要因素。为了充分隔离基础振动对超精密机床的影响,目前国内外均采用空气弹簧作为隔振元件,并取得了一定的效果,但是这属于被动隔振,这类隔振系统的固有频率一般在2Hz左右。
这种被动隔振方法难以满足超精密加工对隔振系统的要求。为了解决这个问题,有必要研究被动隔振和主动隔振控制相结合的混合控制技术。其中,主动隔振控制系统采用状态空间法设计。
二:实验目的
通过本次上机实验,使同学们熟练掌握:
1. 控制系统机理建模;
2. 时域性能指标与极点配置的关系;
3. 状态反馈控制律设计;
4. MATLAB语言的应用。
三:工程背景的物理描述
图1
图1表示了亚微米超精密车床隔振控制系统的结构原理,其中被动隔振元件为空气弹簧,主动隔振元件为采用状态反馈控制策略的电磁作动器。此为一个单自由度振动系统,空气弹簧具有一般弹性支承的低通滤波特性,其主要作用是隔离较高频率的基础振动,并支承机床系统。主动隔振系统具有高通滤波特性,其主要作用是有效地隔离较低频率的基础振动。主、被动隔振系统相结合可有效地隔离整个频率范围内的振动。
经物理过程分析得出床身质量的运动方程为:
(1)
——空气弹簧所产生的被动控制力。
——作动器所产生的主动控制力。
假设空气弹簧内为绝热过程,则被动控制力可以表示为:
(2)
——标准压力下的空气弹簧体积;
——相对位移(被控制量);
——空气弹簧的参考压力;
——参考压力下单一弹簧的面积;
——参考压力下空气弹簧的总面积;
——绝热系数。
电磁作动器的主动控制力与电枢电流、磁场的磁通量密度及永久磁铁和电磁铁之间的间隙面积有关,这一关系具有强非线性。
由于系统工作在微振动状况,且在低于作动器截止频率的低频范围内,因此主动控制力可近似线性化地表示为:
(3)
——力-电流转换系数;
——电枢电流。
其中,电枢电流满足微分方程:
(4)
——控制回路电枢电感系数;
——控制回路电枢电阻;
——控制回路反电动势;
——控制电压。
四:闭环系统的性能指标要求
要求闭环系统单位阶跃响应的超调量不大于5%,过渡过程时间不大于0.5秒()。
五:车床振动系统的开环状态空间模型的建立
首先假定为常数,将式两边求关于时间的二阶导数可得:
(5)
记为:
(6)
其中。
对式(6)两边求导得:
(7)
由式(6)可得:
(8)
由式(7)可得:
(9)
将式(8)和(9)代入式(4)可得:
即:
将非线性项视为干扰信号,略去不计,可得线性化模型为:
(10)
令状态变量为:
,,
可得系统开环状态方程为:
由此得开环系统的状态空间表达式为:
(11)
假设某一亚微米超精密车床隔振系统的各个参数为:
,,,,,。
代入式(11)得开环系统的状态空间表达式为:
六:状态反馈控制律的设计
根据性能指标,解得,所以。
根据性能指标,解得。
留出裕量,取,,则:,。为此得两共轭极点为,,取第三个极点为。于是得出系统期望特征多项式为:(12)
设状态反馈控制律为:
则闭环系统的状态空间表达式为:
则此时闭环系统的特征多项式为: (13)
将式(12)与式(13)比较可得:
, 最终解得: ,至此状态反馈控制律设计完毕。
七:闭环系统的数字仿真
1. 闭环系统的单位阶跃响应仿真
由以上设计过程,借助matlab画出系统的simulink仿真图如图2:
图2
得出系统此时的阶跃响应曲线如图3,图4:
图3
图4
由图4分析可见系统输出最大值为-0.0005585,稳态值为-0.0005503.超调量为:1.49%,此时已进入2%误差带,调整时间小于0.5秒。设计满足要求。
2. 闭环系统的全状态响应仿真
假设存在某一初始振动状态,,。
根据闭环齐次状态方程:
编写matlab程序如下:
第一个文件simu046.m
function dx=simu046(t,x)
A=[0,1,0; 0,0,1; -15625,-2156.25,-120];
dx=A*x;
第二个文件do_simu046.m
function do_simu046
[t,x]=ode45(s
您可能关注的文档
最近下载
- FLEX SYSTEM EN2092 1GB ETHERNET调试文档模板.pdf VIP
- 2024年9月中国汽车产销报告.pdf VIP
- 2024年11月广东深圳市福田区选用机关事业单位辅助人员和社区专职工作者143人笔试备考题库及答案解析.docx VIP
- 流体力学 课件 第6章 气体射流.pptx
- 武汉理工大学材料科学基础课件(全套).pdf
- 风貌改造工程施工组织设计.doc
- 药品生物检定技术 全套教案.doc
- 浅谈双碳战略下的自来水厂碳核算碳减排经验分享_供水技术讲座资料课件.pptx
- 2024广州铁路局机考.pdf
- 11《动物的眼睛》教学设计-2023-2024学年科学二年级下册青岛版.docx VIP
文档评论(0)