首页 > 领导讲话 > 党务讲话 / 正文
桥式吊车小车运动控制系统的建模及MATLAB仿真(附程序)
2020-11-22 05:40:45 ℃桥式吊车运动控制系统的建模及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
猜你喜欢
- 2024-01-13 2024乡镇优秀组织委员优秀党务工作者先进事迹材料
- 2024-01-08 2024年度优秀党务工作者事迹材料(范文推荐)
- 2024-01-07 (3合集)关于党务论谈材料汇编
- 2023-12-27 2024年度优秀党务工作者事迹材料(范文推荐)
- 2023-09-26 2023关于党务干部培训方案6篇
- 2023-07-15 党建党务知识测试题(范文推荐)
- 2023-07-15 市文化和旅游局优秀党务工作者先进事迹材料
- 2023-07-14 2023年某局优秀党务工作者简要事迹材料
- 2023-07-05 优秀党务工作者先进事迹材料()
- 2023-07-02 2023年全市公安机关党务干部培训会议动员讲话
- 搜索
-
- 关于进一步规范干部借调工作的实施办法 04-08
- 2021国家开放大学电大本科《国际私法》 11-03
- 组织生活会专题学习心得体会 09-07
- “党建+安全”主题党日活动方案 03-31
- 镇退役军人服务站情况汇报 11-04
- 纪委监委意识形态工作经验交流:发挥专 05-12
- 看《即将到来的对华战争》观后感 07-24
- 劳务派遣人员薪酬管理办法 04-07
- 民族团结进步创建工作情况总结 04-10
- 2020年全面从严治党履行“一岗双责”情 07-09
- 11-25国庆70周年庆典晚会 庆典晚会串词
- 11-25办公室礼仪的十大原则 浅谈办公室的电话礼仪
- 01-17用心灵轻轻地歌唱_心灵的歌唱
- 01-17也许你不是我一生的唯一|也许不是我
- 01-17爱了,请珍惜;不爱,趁早放手|爱就珍惜不爱就放手
- 01-17岁月带走的是记忆,但回忆会越来越清晰|有趣又有深意的句子
- 01-17曾经的美好只是曾经,我只想珍惜身边的人|我只想珍惜你
- 01-18从容不惊 [学会笑眼去看世界,不惊不乍,淡定从容]
- 02-03当代大学生学习态度调查报告
- 02-03常用护患英语会话
- 标签列表