首页 > 演讲致辞 > 会议致辞 / 正文
桥式吊车小车运动控制系统的建模及MATLAB仿真(附程序)
2020-11-22 10:05:30 ℃桥式吊车运动控制系统的建模及MATLAB仿真(附程序) 1 简介 桥式起重机是横架于车间、仓库及露天堆场的上方,用来吊运各种物体的机械设备,通常称为“天车”或 “吊车”。它是机械工业、冶金工业和化学工业中应用最广泛的一种起重机械。实际生产中的桥式吊车(天车)类似,是一个MIMO复杂控制系统。
桥式吊车系统由三部分组成:桥架驱动系统,车体驱动系统和重物装吊系统。其工作流程为:先将重物起吊至预先设定好的高度,然后吊车运动将重物运到想要放置的位置上方,最后把重物下放到想要放置的位置上。
2 确定要研究的系统为桥式吊车运动控制系统 桥式吊车系统工作示意图见下图1:
p M m x z F θ mg 图1 桥式吊车工作示意图 对于如上桥式吊车控制系统,首先做如下假设:
1) 吊车的行走运动仅限于吊车一个自由度,即假设桥架不运动,只有吊车在桥架上行走。
2) 吊车行走时吊装重物的绳索长度不变。
图中,x坐标为水平方向,z坐标为垂直方向。重物的摆动是由吊车与重物的运动产生的,可以根据动力学有关规律建立吊车及重物的运动方程式。
1) 在水平方向,吊车和重物整体受力为F(t),由牛顿第二定律得 (1) 2) 在垂直于绳索方向,重物受力为,由牛顿第二定律得 (2) 由吊车在行走时吊装重物的绳索长度不变的假设可得出下面两个关系式:
(3) (4) 式中,为绳索长度。
由(3)可得 (5) (5)代入(1)得:
(6) 同样由式(4)可得:
(7) 将(5)(7)代入(2)得 (8) 又尽量小,所以有如下近似式:
,, 将(6),(8)线性化可得:
(9) (10) 由(9)和(10)计算得 (11) 和 (12) 3) 吊车驱动装置的方程式。吊车由电动机驱动,简化的认为电动机是一个时间常数为的一阶惯性环节,即它产生的驱动力F(t)与其控制电压v(t)之间满足方程式:
(13) 其中K为放大系数。
3 选择系统的输入、输出变量和状态变量 选择5个状态变量分别为:
,,,,;
输入变量为:;
两个输出变量为:,。
4 建立状态空间描述 根据(11)(12)(13)式可得出描述吊车运动系统的状态空间表达式为:
选取适当参数:对一个实际的桥式吊车吊车运动系统,假定具有如下各具体参数:M=1000kg,m=4000kg,l=10m,K=100N/V。将它们代入上面的状态空间表达式得:
5 分析系统的稳定性 用特征值法。在MATLAB中输入以下程序:
eig(A) ans = 0 0 0 + 2.2136i 0 - 2.2136i -1 系统的5个开环特征值不全位于S左平面上,有4个位于虚轴上,所以系统为临界不稳定。
系统输出仿真波形如下图所示:
6 判断系统的能控性 使用MATLAB判断系统的能控性,输入以下程序:
A=[0 1 0 0 0;0 0 -39.2 0 0.001;0 0 0 1 0;0 0 -4.9 0 0.0001;0 0 0 0 -1]; B=[0;0;0;0;100]; C=[1 0 0 0 0;0 0 1 0 0]; rct=rank(ctrb(A,B)) rct = 5 根据判别系统能控性的定理,该系统的能控性矩阵满秩,所以该系统是能控的。
7 采用状态反馈进行系统综合 因为系统是能控的,所以,可以通过状态反馈来任意配置极点。例如将极点配置在:s1=-0.16-j0.16 s2=-0.16+j0.16 s3,s4,s5=-1。
在MATLAB中输入:
P=[-0.16+0.16i,-0.16-0.16i,-1,-1,-1]; K=acker(A,B,P) 求出状态反馈矩阵K:
K = 0.52245 4.8327 -1420.7 -137.21 0.0232 在MATLAB中输入 A-B*K ans = 0 1 0 0 0 0 0 -39.2 0 0.001 0 0 0 1 0 0 0 -4.9 0 0.0001 -52.245 -483.27 1.4207e+005 13721 -3.32 因此综合后系统的状态空间描述为:
采用MATLAB编写m文件进行仿真。
运行仿真程序,得到仿真曲线如下图:
将极点配置在:s1=-0.2 s2=-0.2 s3,s4,s5=-1。计算出K=[0.40816 5.3061 -1438.1 -119.06 0.024]此时输出y的仿真曲线如下:
将极点配置在:s1=-0.16+0.16i s2=-0.16-0.16i s3=-1 s4=-2 s5=-3。计算出K=[3.1347 25.339 -2145.4 553.73 0.0532]此时输出y的仿真曲线如下:
8 实验结论 通过比较3组不同的极点配置下状态反馈系统的输出响应曲线和原系统的输出响应曲线可以看出,不同的极点配置对系统性能有一定的影响,但只要极点都配置在S左平面,就可以保证系统具有一定的动态和稳态性能。
附录:程序 %******************************************* %桥是吊车运动控制系统极点配置设计及仿真 2015.12.30 lwd %******************************************* %建立状态空间表达式 A = [0,1,0,0,0;0,0,-39.2,0,10^(-3);0,0,0,1,0;0,0,-4.9,0,10^(-4);0,0,0,0,-1] B = [0;0;0;0;100] C = [1,0,0,0,0;0,0,1,0,0] D = [0;0] %分析系统稳定性 eig(A) %求A的特征值,通过特征值在S平面分布判断系统稳定性 %系统仿真输出波形 [num,den] = ss2tf(A,B,C,D); sys1 = tf(num(1,:),den); sys2 = tf(num(2,:),den); figure step(sys1) figure step(sys2) %判断系统的能控性 rct = rank(ctrb(A,B)) obs = rank(obsv(A,C)) %如果能空则进行极点配置设计 if 5 == rct %第一次极点配置设计 P = [-0.16+0.16*i,-0.16-0.16*i,-1,-1,-1]; K = acker(A,B,P) A1 = A-B*K [num,den] = ss2tf(A1,B,C,D); sys1 = tf(num(1,:),den); sys2 = tf(num(2,:),den); figure step(sys1) figure step(sys2) %第二次极点配置设计 P = [-0.2,-0.2,-1,-1,-1]; K = acker(A,B,P) A1 = A-B*K [num,den] = ss2tf(A1,B,C,D); sys1 = tf(num(1,:),den); sys2 = tf(num(2,:),den); figure step(sys1) figure step(sys2) %第三次极点配置设计 P = [-0.16+0.16*i,-0.16-0.16*i,-1,-2,-3]; K = acker(A,B,P) A1 = A-B*K [num,den] = ss2tf(A1,B,C,D); sys1 = tf(num(1,:),den); sys2 = tf(num(2,:),den); figure step(sys1) figure step(sys2) end
猜你喜欢
- 2021-10-06 2021反法西斯胜利抗战心得
- 2021-10-06 名著红岩的读书心得读书笔记5篇
- 2021-10-05 弘扬伟大长征精神中心组发言
- 2021-10-02 建设世界重要人才中心和创新高地青年心得体会【五篇】精选
- 2021-05-06 区县义务教育均衡发展实验区县建设方案
- 2021-05-06 疫情防控个人现实表现
- 2021-05-06 学习科学发展观心得体会:坚持“六个必须”
- 2021-05-06 入党申请书思想方面怎么写
- 2021-05-05 校园安全方面演讲稿
- 2021-05-05 2篇基层党员干部政法队伍教育整顿民主生活会个人对照检查材料范文
- 搜索
-
- 重点信访人员处理及稳控措施汇报 04-22
- 「部编版」五年级上语文期末,课内阅读, 10-10
- 成立新党支部请示 10-19
- 十四五规划编制心得体会思路发言5 09-13
- 硕士研究生中期考核个人总结 09-27
- 部队基层管理案例分析及心得体会 10-28
- 学习全国公安工作会议精神心得体会 04-23
- 党员学习四史心得体会3篇 08-11
- 对当前信访维稳工作对策与建议 05-10
- 幼儿园新冠肺炎疫情防控各岗位工作责任 03-25
- 11-25国庆70周年庆典晚会 庆典晚会串词
- 11-25办公室礼仪的十大原则 浅谈办公室的电话礼仪
- 01-17用心灵轻轻地歌唱_心灵的歌唱
- 01-17也许你不是我一生的唯一|也许不是我
- 01-17爱了,请珍惜;不爱,趁早放手|爱就珍惜不爱就放手
- 01-17岁月带走的是记忆,但回忆会越来越清晰|有趣又有深意的句子
- 01-17曾经的美好只是曾经,我只想珍惜身边的人|我只想珍惜你
- 01-18从容不惊 [学会笑眼去看世界,不惊不乍,淡定从容]
- 02-03当代大学生学习态度调查报告
- 02-03常用护患英语会话
- 标签列表