现代控制理论课程设计-大型桥式吊车行车控制系统的状态空间设计

《现代控制理论课程设计-大型桥式吊车行车控制系统的状态空间设计》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、.现代控制理论课程设计题目：大型桥式吊车行车控制系统的状态空间设计班级：学号:姓名：指导教师：15一、项目背景随着我国制造业的发展，桥式起重机越来越多的应用到工业生产当中。在工厂中搬运重物，机床上下件，装运工作吊装零部件，流水线上的定点工作等都要用到起重机。起重机中种数量最多，在大小工厂之中均有应用的就是小吨位的起重机，小吨位的桥式起重机广泛的用于轻量工件的吊运，在我国机械工业中占有十分重要的地位。但是，我国现在应用的各大起重机还是仿造国外落后技术制造出来的，而且已经在工厂内应用了多年，有些甚至还是七八十年代的产品，无论在质量上还是在

2、功能上都满足不了日益增长的工业需求。桥式起重机是横架于车间、仓库及露天堆场的上方，用来吊运各种物体的机械设备，通常称为“天车”或“行车”“吊车”。它是机械工业、冶金工业和化学工业中应用最广泛的一种起重机械。实际生产中的桥式吊车（天车）与倒立摆类似，是自动控制最为经典的实验模型，是一个MIMO复杂控制系统，可以作为控制理论算法研究的理想实验平台。桥式吊车系统由三部分组成：桥架驱动系统，小车驱动系统和重物撞吊系统。其工作流程为：先将重物起吊至预先设定好的高度，然后小车运动将重物运到想要放置的位置上方，最后把重物下放到想要放置的位置上。二、

3、问题提出桥式吊车系统工作示意图见下图1：pMmxzFθmg图1桥式吊车工作示意图2.1被控对象的动态分析对于如上桥式吊车控制系统，首先做如下假设：15①吊车的行走运动仅限于小车一个自由度，即假设桥架不运动，只有小车在桥架上行走。②小车行走时吊装重物的绳索长度不变。图中，x坐标为水平方向，即小车运动自由度，z坐标为垂直方向，即重物运动自由度。重物的摆动是由小车与重物的运动产生的，可以根据动力学有关规律建立小车及重物的运动方程式。(1)小车的运动方程式为：（1）式中M为小车质量（包括了重物吊装机等）；为小车的位移；F为驱动装置所产生的对小

4、车的驱动力；p为重物通过绳索对小车产生的拉力；θ为重物的摆动角度（即绳索与垂直方向的夹角）；(2)重物的水平运动方程为：（2）式中，m为重物质量（包括了吊钩等）；为重物的水平方向位移。(3)重物垂直运动方程为：（3）式中，为重物垂直方向位移。由小车在行走时吊装重物的绳索长度不变的假设可得出下面两个关系式：（4）（5）式中，l为绳索长度。由（1-1），（1-2），（1-3）可得出：（6）（7）由（1-4）可得：代入（1-6）得：15同样由式（1-7）可得：又θ(t)尽量小，所以有如下近似式：，，经线性化和计算可得：(一)(二)（4）小车

5、驱动装置的方程式。小车由电动机驱动，简化的认为电动机是一个时间常数为Td的一阶惯性环节，即它产生的驱动力F(t)与其控制电压v(t)之间满足方程式：(三)其中K为放大系数。2.2选择系统的输入、输出变量和状态变量选择5个状态变量分别为：，，，，控制量为：　　两个输出变量为：　　，我们只对重物的水平位移即输出变量作分析。2.3　建立状态空间描述根据㈠㈡㈢式可得出描述小车运动系统的状态空间表达式为：15选取适当参数:对一个实际的桥式吊车小车运动系统，假定具有如下各具体参数：M=1000kg,m=4000kg,l=10m,K=100N/V。

6、将它们代入上面的状态空间表达式得：基于以上讨论，可以把桥式吊车小车运动过程自动控制系统表述为，一个对小车驱动电机的合适控制，使小车系统在桥架上从初始位置精确地运动到一个事先规定的位置，并保证具有良好的动态特性，特别是重物的摆动角度尽量小，小车运动到规定的位置时重物的摆角为零。2.4被控对象的特征值(8)作为被控对象的的起重机系统，其对应的（开环）特征值，可将式（8）带人对应的（开环）特征方程求出，即(9)求解式（9），得（开环）特征值(10)代人实际参数后，求得特征值为（11）2.5调节对象（起重机系统）自身动态分析由式（10）知，此

7、调节对象5个（开环）特征值中，有两个位于坐标原点，两个位于虚轴，一个位于负实轴，将这5个（开环）特征值的分布与系统结构图结合起来分析可知：15描述的装置驱动的特征，由于该装置系一串联接入的一阶惯性环节，因此其对应的特征值将为负实数并可单独给以分析。描述的是小车的动力学特征，这是因为在与之间，也就是在与之间相当于存在两个相互串联的积分环节，且无反馈支路存在，显然，两个位于坐标原点的特征值将是与此相对应的。这样一对共轭虚数特征值描述的将是吊钩的无阻尼振荡的动力学特征。这是因为闭环负反馈子系统对应的传递函数为，显然，其对应的一对极点即为。3

8、.建模及仿真①用特征值法。在MATLAB中输入以下程序：A=[01000;00-39.200.001;00010;00-4.900.0001;0000-1];eig(A)ans=000+2.2136i0-2.2136i-