大型桥式吊车行车控制系统的状态空间设计.docx
《大型桥式吊车行车控制系统的状态空间设计.docx》由会员分享,可在线阅读,更多相关《大型桥式吊车行车控制系统的状态空间设计.docx(16页珍藏版)》请在冰豆网上搜索。
大型桥式吊车行车控制系统的状态空间设计
^Vbrtfvwest"UniversityforO^ationaRties
现代控制理论课程设计
题目:
大型桥式吊车行车控制系统的状态空间设计
班级:
学号:
姓名:
指导教师:
、项目背景
随着我国制造业的发展,桥式起重机越来越多的应用到工业生产当中。
在工厂中搬运
重物,机床上下件,装运工作吊装零部件,流水线上的定点工作等都要用到起重机。
起重机中种数量最多,在大小工厂之中均有应用的就是小吨位的起重机,小吨位的桥式起重机广泛的用于轻量工件的吊运,在我国机械工业中占有十分重要的地位。
但是,我国现在应用的各大起重机还是仿造国外落后技术制造出来的,而且已经在工厂内应用了多年,有些甚至还是七八十年代的产品,无论在质量上还是在功能上都满足不了日益增长的工业需求。
桥式起重机是横架于车间、仓库及露天堆场的上方,用来吊运各种物体的机械设备,通常称为“天车”或“行车”“吊车”。
它是机械工业、冶金工业和化学工业中应用最广泛的一种起重机械。
实际生产中的桥式吊车(天车)与倒立摆类似,是自动控制最为经典的实验模型,是一个MIMC复杂控制系统,可以作为控制理论算法研究的理想实验平台。
桥式吊车系统由三部分组成:
桥架驱动系统,小车驱动系统和重物撞吊系统。
其工作流程为:
先将重物起吊至预先设定好的高度,然后小车运动将重物运到想要放置的位置上方,最后把重物下放到想要放置的位置上。
二、问题提出
桥式吊车系统工作示意图见下图1:
2.1被控对象的动态分析
对于如上桥式吊车控制系统,首先做如下假设:
1吊车的行走运动仅限于小车一个自由度,即假设桥架不运动,只有小车在桥架上行走。
2小车行走时吊装重物的绳索长度不变。
图中,x坐标为水平方向,即小车运动自由度,z坐标为垂直方向,即重物运动自由度。
重物的摆动是由小车与重物的运动产生的,可以根据动力学有关规律建立小车及重物的运动方程式。
(1)小车的运动方程式为:
MXM(t)=F(t)-psinr(t)
(1)
式中M为小车质量(包括了重物吊装机等);Xm为小车的位移;F为驱动装置所产生的对小车的驱动力;p为重物通过绳索对小车产生的拉力;9为重物的摆动角度(即绳索与垂直方向的夹角);
(2)重物的水平运动方程为:
mXm(t)二psin=(t)
(2)
式中,m为重物质量(包括了吊钩等);Xm为重物的水平方向位移。
(3)
重物垂直运动方程为:
式中,Zm为重物垂直方向位移
由小车在行走时吊装重物的绳索长度不变的假设可得出下面两个关系式:
(5)
(6)
(7)
Zm(t)=ICOS,(t)
式中,I为绳索长度。
由(1-1),(1-2),(1-3)可得出:
Mxm(Umxm(t)讦(t)
v»v>
mxm(t)cos^(t)mzm(t)si(t)二mgsin^(t)
由(1-4)可得:
»«**2
Xm(t)二XM(t)Tcosv(tp(t)lsinl(t)
(Mm)xM(t)-mlcosw(t)mlsin二⑴二(t)=F(t)
同样由式(1-7)可得:
2
Zm(t)=—lcos日日(t)—lsin日(t)日(t)
又B(t)尽量小,所以有如下近似式:
2
sin日(t)乏日(t)cos日(t)^1sine(t)日(t)^0
经线性化和计算可得:
(4)小车驱动装置的方程式。
小车由电动机驱动,简化的认为电动机是一个时间常数为Td的一阶惯性环节,即它产生的驱动力F(t)与其控制电压v(t)之间满足方程式:
TdF(t)F(t)二Kv(t)
其中K为放大系数。
2.2选择系统的输入、输出变量和状态变量
选择5个状态变量分别为:
X^Xm,X^Xm,X3",X4",X5二F
控制量为:
U=V
两个输出变量为:
%=XM,y2
我们只对重物的水平位移即输出变量%二Xm作分析
2.3建立状态空间描述
根据㈠㈡㈢式可得出描述小车运动系统的状态空间表达式为:
10000y=Ix
[00100一
选取适当参数:
对一个实际的桥式吊车小车运动系统,假定具有如下各具体参数:
M=1000kg,m=4000kg,l=10m,K=100N/V将它们代入上面的状态空间表达式得:
■0
1
0
0
01
[
-01
0
0
-39.2
0
-3
10
0
0
0
0
1
0
x+
0
0
0
-4.9
0
-4
10
0
卫
0
0
0
一1一
100一
x
u
基于以上讨论,可以把桥式吊车小车运动过程自动控制系统表述为,一个对小车驱动
电机的合适控制V(t),使小车系统在桥架上从初始位置精确地运动到一个事先规定的位置
Xm=r,并保证具有良好的动态特性,特别是重物的摆动角度讯t)尽量小,小车运动到
规定的位置时重物的摆角为零。
2.4被控对象的特征值
(开环)特征方程det('l-A)=0求出,即
1
T(10)
1,2=0,3,4二-j
代人实际参数后,求得特征值为
%2=0,*4=±j2.21sT,=_sT
2.5调节对象(起重机系统)自身动态分析
由式(10)知,此调节对象5个(开环)特征值中,有两个位于坐标原点,两个位于虚轴,一个位于负实轴,将这5个(开环)特征值的分布与系统结构图结合起来分析可知:
1
5丄描述的装置驱动的特征,由于该装置系一串联接入的一阶惯性环节,因此其对应
5Ta
的特征值将为负实数并可单独给以分析。
,描述的是小车的动力学特征,这是因为在Xi与X2之间,也就是在Sa与Sa之间相当于
存在两个相互串联的积分环节,且无反馈支路存在,显然,两个位于坐标原点的特征值将是与此相对应的。
(mAmB)g
'3,4二-J\a43二-J—-这样一对共轭虚数特征值描述的将是吊钩
V叫1
的无阻尼振荡的动力学特征。
这是因为闭环负反馈子系统对应的传递函数为,
(1/S2)1
2=二显然,其对应的一对极点即为打4。
(1(1/s2)a43)S2*433,4
3.建模及仿真
1用特征值法。
在MATLAB^输入以下程序:
A=[01000;00-39.200.001;00010;00-4.900.0001;0000-1];
eig(A)
ans=
0
0
0+2.2136i
0-2.2136i
-1.0000
系统的5个开环特征值不全位于S左平面上,有4个位于虚轴上,所以系统为临界不稳定。
2也可以采用MATLAB/Simulink构造系统开环控制系统的仿真模型,如下图所示。
其输出仿真波形如下图所示:
由上图可判断系统是不稳定的
③在MATLAB^输入以下程序,可得原系统的阶跃响应:
A=[01000;00-39.200.001;00
title('x2');
010;00-4.900.0001;00
00-1];
subplot(5,1,3);
B=[0;0;0;0;100];;
plot(t,x(:
3));grid
C=[10000];
D=0;
xlabel('t(s)');ylabel('x3(t)');
sys0=ss(A,B,C,D);
title('x3')
t=0:
0.01:
50;
subplot(5,1,4);
[y,t,x]=step(sys0,t);
plot(t,x(:
3));grid
subplot(5,1,1);
plot(t,x(:
1));grid
xlabel('t(s)');ylabel('x4(t)');
xlabel('t(s)');ylabel('x1(t)');
title('x4')
title('x1');
subplot(5,1,5);
subplot(5,1,2);
plot(t,x(:
3));grid
plot(t,x(:
2));
grid;
xlabel('t(s)');ylabel('x5(t)');
xlabel('t(s)');ylabel('x2(t)');
title('x5')
x2
2
XJyL2x
八Cd
、
r\、
\\f\f\
AJ
/VVV
/\/
x4
o
di
X
o
di
X
八/Xf\H
\f\C\
Pta.
从5个状态变量的波形可看出系统是不稳定的。
4.判断系统的能控性
使用MATLAB^」断系统的能控性,输入以下程序:
A=[01000;00-39.200.001;00010;00-4.900.0001;0000-1];B=[0;0;0;0;100];
C=[10000;00100];
rct=rank(ctrb(A,B))
ret=
5
根据判别系统能控性的定理,该系统的能控性矩阵满秩,所以该系统是能控的。
5.采用状态反馈进行系统综合
因为系统是能控的,所以,可以通过状态反馈来任意配置极点。
例如将极点配置在:
s仁-0.16-j0.16s2=-0.16+j0.16s3,s4,s5=-1。
在MATLAB^输入:
A=[01000;00-39.2010A-3;00010;00-4.9010A-4;0000-1];B=[0;0;0;0;100];
P=[-0.16+0.16i,-0.16-0.16i,-1,-1,-1];
K=acker(A,B,P)
K=
1.0e+003*
0.00050.0048-1.4207-0.13720.0000
因此,求出状态反馈矩阵为:
1
根据公式求出输入变换系数:
C(BK-A)B
在MATLAB^输入:
B=[0;0;0;0;100];
C=[10000];
K=[0.54.8-1420.7-137.20.0232];
t=C*inv(B*K-A)*Bt=
2
得:
1=0.5
采用MATLAB/Simulink构造系统状态反馈控制系统的仿真模型,如下图所示
当输入为1时仿真波形如下
当输入为8时:
其仿真波形为:
Timeoffset:
0
在MATLAB^输入以下程序可得加入状态反馈后,系统的各状态变量阶跃响应:
A=[01000;00-39.200.001;00010;00-4.900.0001;0000-1];
B=0.5*[0;0;0;0;100];
C=[10000];
K=2*[0.54.8-1420.7
0.0232];
D=0;
A仁A-B*K;
sys0=ss(A1,B,C,D);
t=0:
0.01:
50;
[y,t,x]=step(sys0,t);subplot(5,1,1);plot(t,x(:
1));gridxlabel('t(s)');ylabel('x1(t)');title('x1')