《龙门吊车重物防摆的鲁棒-PID-控制方案.docx》由会员分享，可在线阅读，更多相关《龙门吊车重物防摆的鲁棒-PID-控制方案.docx（32页珍藏版）》请在课桌文档上搜索。
1、龙门吊车重物防摆的鲁棒PID控制方案周耀龙144173248摘要龙门吊车作为一种运输工具，广泛应用在现代工厂、安装工地和集装箱货仓等的装卸与运输作业。他在离地面很高的轨道上运行，具有占地面积小、省时省工的优点。关键词龙门吊车；双闭环；PlD龙门吊利用绳索一类的柔性体代替刚体工作。由于惯性，运动过程中会使吊重产生摇摆，不利于起重机的快速对位。文中采用拉格朗日方程的方法建立了龙门起重机的动力学模型，并用MATLAB仿真功能验证了数学模型的有效性。然后设计了防摇摆的双闭环PID控制方案，并合理选择参数，使控制具有较强的鲁棒性。并用Simulink进行仿真实验，验证控制方案的合理性。一、前言桥式起重机
2、或门式起重机广泛用于车站、码头、仓库、工厂等场所搬运物料，是工厂、铁路、港口及其他部门实现物料搬运机械化的重要设备。尤其是轨行式集装箱龙门式起重机是集装箱堆场的主要装卸机型，作为现代物流装备之一其应用得到逐步推广。当起重机小车或大车运行时，控制起重机的起吊重物相对于小车中心竖直线的偏摆幅度,可以减小吊重的晃动程度，从而实现起重机的快速对位，如吊具与集装箱对位，起吊集装箱与底盘车对位等。以及集装箱在堆场的准确码放，以提高装卸作业效率。控制吊重摇摆程度的方法一种是采取防摇措施，主动控制小车或大车使吊重从静止运动到目标位置过程中始终保持吊重在较小范围内摆动，另一种是采取减摇措施，当吊重偏角较大时，被
3、动控制小车或大车使吊重摇摆的幅度在最短时间内衰减到规定范围内。起重机的这种主动防摇和被动减摇问题可归结为起重机的快速对位问题。起重机吊重防摇控制系统就是使吊重的摆动能得到迅速衰减，在较短内使吊重相对于小车的中心竖直线处于微动状态（即在规定的微小角度内摆动），以利于吊具在工作空间准确对位和集装箱等吊重准确、快速码放，达到起重机快速对位的目的，从而极大地提高起重机的装卸作业效率，明显改善装卸作业安全状态，减轻操作人员的工作强度，消除操作人员之间的经验差别，减少甚至消除快速对位对操作人员经验的过分依赖性。二、系统建模(一)机理建模龙门吊车利用绳索一类的柔性体代替刚体工作，以使得吊车的结构轻便，工作效
4、率高。但是，采用柔性体吊运也带来了一些负面影响，例如吊车负载一一重物的摆动问题一直是困扰吊车装运效率的一个难题。为研究吊车的防摆动控制问题，需要对实际问题进行简化、抽象。吊车的“搬运一一行走一一定位”过程可抽象为如图的模型：图1吊车系统的物理抽象模型图中，小车的质量设为加O,受到水平方向的外力口的作用，重物的质量为加，绳索的长度为L,对重物的快速吊运与定位问题可以抽象为：小车在受到外力尸作用时，使得小车在最短的时间4由A点运动到B点，且摆(),为系统允许的最小摆角。该问题为多刚体、多自由度、多约束的质点系动力学问题，若应用牛顿力学来分析过于复杂，因此采用拉格朗日方程，将力学体系中运动方程从以力
5、为基本单位的牛顿形式，改变为以能量为基本概念的分析力学形式。拉格朗日的普遍形式为：d/肛Tdt两*qk式中，T为质点系的动能，/为质点系的广义坐标，k为质点系的自由度数，匕为广义力。由此可见，拉格朗日方程把力学体系的运动方程从以力为基本概念的牛顿形式，改变为以能量为基本概念的分析力学形式。(二)系统建模实际中的吊车系统受到多种干扰，如小车与导轨之间的干摩擦、风力的影响等。为了便于分析，需对实际系统进行进一步的简化。简化为图所示的物理模型：图2龙门吊车的物理模型重物通过绳索与小车相连，小车在行走电机的水平拉力Fl(N)的作用下载水平轨道上运动，小车的质量为加。(kg),重物的质量为机(kg),绳
6、索的长度为L(m),重物可在提升电机的提升力F2(N)的作用之下进行升降运动；绳索的弹性、质量、运动的阻尼系数可忽略；小车与水平轨道的摩擦阻尼系数为D(kg/s)；重物摆动时的阻尼系数为D(KgmVs),其他扰动可以忽略。取小车的位置为王，绳子长度为当，摆角为再作为系统的广义坐标系，在此基础上对系统进行动力学分析。由上图的坐标系可知，小车的位置和重物的位置坐标为:4)=MJ1=0xm=x1-x2sinx3l=s所以小车和重物的速度分量为：F3J4=。Ixfn=xx-x2sinx3-x2x3cosx3l.1=x2s-x2smx3系统的动能为:T=mv2+mv22wjo2,n=-m(3+/)+m(
7、2+y)2。/nowo2Jm-(m0+m)+m(x+考后-2xix2sinx3-2x1x2x3cosx3)22此系统的拉格朗日方程组为:dfT.Tj7n奈五)一成二F以(史1)-21=F2+mgcosx3dtx2x2d,T、dT_()mgx?SinVX3、出x3x3综合以上公式的系统的方程组为:(m.m0)x1-mx2sinx3-wr3cosx3-2mr2x3cosx3nr2x3sinx3I+Dx1=F1jmx2-mxxsinx3-nxy-mgcosX3=F2Ilmxx3+2nc2x1x3-mxix2cosx3mgx2sinx3+v=0上式是考虑绳子长度变化的情况下的二自由度龙门吊车的运动系统
8、的动力学模型。对于绳子长度不变的情况下，可将上述模型进一步简化，将上式中的：X2=X2=0消去F2,令F=Fl,x2=l=常数，得到绳长不变的情况下运动系统的数学模型为：fI(m+m)X-mlx3cosx3+mlx3sinx3Dx=FIm-r,vclcosx3+mgsinx3+vx3=0(三)模型简化由上式可见，龙门吊车的运动系统的动力学模型为非线性微分方程组。为了便于应用经典控制理论对该控制系统进行设计，必须将其简化为线性定常的系统模型。考虑到实际吊车运行过程中摆动角较小,一般不超过10,且平衡位置为=x3=0,将模型在0=0处进行线性化，此时有如下近似结果：sinjcos1;62sin0考
9、虑到摆动的阻尼系数较小，可以认为R,所以上式可以简化为：Jl(Jnm0)x-加/9+Dx=FIml9-ITix+mg9=0进一步简化为：F=m0X+mg9+Dxx=9+g9对上式进行拉式变化可得：HF(S)=(加+Dy)X(S)+mg(s)IIyX(S)=(/f+g)G)F(s)由上面系统的传递函数形式模型，可得图所示的定摆长吊车运动系统动态结构图，下图就是其中的一种表达形式：图3定摆长吊车运动系统动态结构图（一）去掉反馈环节图4定摆长吊车运动系统动态结构图（二）同理，也可将上述模型转化为状态空间形式，对式进行变换，每个式子只保留一个二次倒数项，可得：L=-2X.9+LFIm。m。m09=_D
10、（m+9+叫mJ取乐.O,O为系统的状态，工,O为系统的输出，则系统的状态空间描述方程为:x=Ax+Buy=cx式中X=x,%,9,91，=F,y=x,91ID一/O D 一知 O O O O l 6 A 6e- A1000uc-o01OJ当不考虑系统的阻尼系数D,即令D=O,则有：iYF(5)=mo(s)+mgQ(s)=(ls2g)Q(s)写成传递函数的形式：iX(S)E+NIQ(S)rQ()IhF(5)m0Is2+(m+mQ)g所描述的系统如图所示:图5系统结构图（四）模型验证数学模型建立完毕，为了检查所建立的模型和实际模型是否具有相同的一些必要性质，要对所建立的模型进行验证。下面给一些物
11、理量赋初值。设绳长为1m,小车质量m0=50kg,重物质量m=5kg,在时间t=0时，给小车一个F=IN的恒力，小车开始位于x=0处，且摆角0二0。根据经验，我们知道，小车将在恒力的作用下做匀加速直线运动，位置不断增加，为一个抛物线曲线。初始状态=0为重物相对小车摆动的一个极限取值，而在恒力作用下，也将使重物相对小车的摆动角存在另外一个极限取值所以，整个摆动运动就是重物在小车的一侧0-，角度之间做往复摆动。所以，在恒力F作用下，小车向前移动，X不断增加，负载的重物在00，区间内摆动，b的取值与初始力F的大小有关。下面利用MATLAB中的Simulink模块进行仿真实验。为了使程序更具有可读性，
12、使用Simulink模块中的封装子程序功能。程序中有两部分，即使用传递函数建立的简化模型和使用微分方程建立的实际模型。F值为输入，X和值为输出，为了便于比较两种模型的输出值的差异，将两种模型的X值和0值放在同一个输出框口中。程序如下：图6利用子系统封装后的模型框图150s2TransferFcnGainTransferFcn 1s2s2+9.8nea图7简化模型的框图图8实际模型的框图Fcn为：(u7-9.8*u8*u3*u4-u8*u6*u5*u5*u3)/(u9+u8*u3*u3)Fcnl为：(u7-9.8*u8*u3*u4-U8*u6*u5*u5*u3)/(u9+u8*u3*u3)*u(
13、4)/u6-9.8*u3/u6图9小车位置坐标X的值Scope1图10重物摆角0从中可以看出，在IN恒力作用下，负载不断的在eOg区间内摆动，小车的位置不断增加。这一结果符合前述的实验设计，故可以在一定程度上确认：该吊车系统的数学模型是有效的。同时，我们也可以看出实际模型和简化模型的曲线基本上是重合的。因此，我们认为近似模型在一定条件下可以表述原系统模型的性质。三、龙门吊防摆的PID控制与仿真实验在设计控制系统过程中，我们假设知道了受控对象和控制器的模型以及它们的各种定常参数。但是，由于存在种种不确定因素，例如参数变化、未建模动力学特性、未建模时延、平衡点的变化、传感器噪声、不可干预的干扰输入
14、等，所以建立起来的对象模型并不能精确的表示实际的物理系统。如何在模型不精确或者存在参数变化的前提下，有效地控制被控对象，尽可能地减小实际系统中这些因素对控制系统品质带来的影响，使系统仍能保持期望的性能，是我们一直面临的问题。在自动控制领域内，控制系统的设计是建立在比较抽象的而且繁琐的数学基础上的,这使得实际工程掌握和运用这些方法较为复杂。PID控制器以其结构简单、稳定性好、工作可靠、调整方便、各个控制器参数具有明显的物理意义而成为工业控制主要和可靠的技术工具。当对一个系统和被控对象不完全了解或者不能通过有效的测量手段来获得系统的参数时，最适合采用PID控制。对于吊车系统的重物防摆控制要求，双闭
15、环PID防摆控制虽具有很好的消摆和定位效果，但针对其绳长和载荷常常不确定，要求所设计的控制系统应具有较强的鲁棒性。下面利用鲁棒PID控制理论对龙门起重吊系统进行设计防摆控制设计，并用MATLAB7.O/Simulink软件进行仿真试验，得到了相应的仿真结果并对结果进行分析。(一)鲁棒PID控制理论下图为PID控制结构框图，典型PID为滞后一超前校正装置。图IlPlD控制结构框图由图可见，PID控制器是通加对误差信号e(t)进行比例、积分和微分运算，其结果的加权，得到控制器的输出u(t),该值就是控制对象的控制值。PID控制器的数学描述为：=KPe+:JOeS+Td式中(以为控制输入，e为误差信
16、号，ZY以为输入量，。为输出量。下面对PID中常用的比例P、比例一积分PI、比例一微分PD和比例一积分一微分PID四种调节器作一简要分析，从而对比例、微分和积分作用有一个初步的认识。(1)比例调节器一比例的作用比例调节器的传递函数GC(S)=Kp,U=Kp(t),即在PID控制器中使Ti!1,Td!O0根据前面所学，为了提高系统的静态性能指标，减少系统的静态误差，一个可行的办法是提高系统的稳态误差系数，即增加系统的开环增益。显然，若使KP增大，可满足上述要求。然而，只有当KP!1,系统的输出才能跟踪输入，而这必将破坏系统的动态性能和稳定性。(2)比例积分调节器一积分的作用在PID调节器中，当T
17、d!0时，控制输出u(t)与e(t)具有如下关系：u(t)=Kpe(t)1e(t)dt首先，通过比较比例调节器和比例积分调节器可以发现，为使e(t)!0,在比例调节器中，!1,这样若e(t)存在较大的扰动，则输出u(t)也很大，这不仅会影响系统的动态性能，也使执行器频繁处于大幅振动中；而若采用PI调节器，如果要求e(t)!O,则控制器输出u(t)由n得到一个常值，从而使输出e(t)稳定于期望的值。其次，从参数调节个数来看，比例调节器仅可调节一个参数Kp,而PI调节器则允许调节参数KP和Ti,这样调节灵活，也较容易得到理想的动、静态性能指标。但是，因PTiSlPl调节器归根到底是一个迟后环节。根
18、据前面介绍的迟后校正原理，在根轨迹法设计中，为避免相位迟后对系统造成的负面影响，零点元靠近原点，即Ti足够大；在频域法设计中，也要求转折频率元且远离!c。这表明在考虑系统稳定性时，Ti应足够大。然而，若Ti太大，则PI调节器中的积分作用变小，会影响系统的静态性能，同时，也会导致系统响应速度的变慢。此时可通过合理调节KP和Ti的参数使系统的动态性能和静态性能均满足要求。采用PI控制，系统的稳态误差为零；且当方减少时，系统的稳定性变差；当Ti增加时，系统的响应速度变慢。(3)PD和PID调节器一微分的作用当PID调节器的Ti!1时，校正装置成为一个PD调节器，这相当于一个超前校正装置，对系统的响应
19、速度的改善是有帮助的。但在实际的控制系统中，单纯采用PD控制的系统较少，其原因有两方面，一是纯微分环节在实际中无法实现，同时，若采用PD控制器，则系统各环节中的任何扰动均将对系统的输出产生较大的波动，尤其对阶跃信号。因此也不利于系统动态性能的真正改善。实际的PID控制器的传递函数如下式：GC（S）=KP（I+式中般大于10。显然，当N!1时，上式即为理想的PID控制器。（二）PID控制器参数整定所设计的PID控制器的控制质量如何，很大程度上取决于其三个参数Kp、7；和7；。因此需要对K。、7；.和7；三个参数进行整定。PID控制器的参数整定的方法有很多，概括起来有两大类：一是理论计算整定，它主
20、要是依据系统的数学模型，经过理论计算确定控制器的参数，这种方法所得到的计算数据还必须通过工程实际进行调整和修改才能在使用。二是工程整定方法，它主要依赖工程经验，直接在控制系统的试验中进行，而且方法简单、易于掌握，在工程实际中被广泛采用。PID控制器参数的工程整定方法通常采用ZiegIeLNiChOIS整定公式。该整定公式是一种针对带有时延环节的一阶系统而提出的实用经验公式。此时，可将系统设定为如下形式：在实际的控制系统中，大量的系统可用此模型近似，尤其对于一些无法用机理方法进行建模的系统，可用时域法和频域法对模型参数进行整定。（三）鲁棒PID控制与灵敏度我们常常希望所设计的系统在不确定参数在一
21、定范围内变化时仍能正常工作。若控制系统是稳健的并有很强的适应能力则是鲁棒控制系统。鲁棒控制系统具有以下特点：灵敏度低；在参数的允许变动范围内能保持稳定；、当参数发生较激烈的变化时，能够恢复和保持预期性能。鲁棒可以视为是系统对那些未加考虑的影响因素的灵敏度，这些影响因素主要包括干扰、测量噪声、和未建模动态特性等。当系统按照设计去完成任务时，它应该能够克服这些因素的影响。灵敏度是控制系统分析与设计的基本问题之一，是用来表征控制系统能受参数变化影响程度的量。当参数只在小范围内摄动时，常用系统灵敏度来度量系统的鲁棒性。系统灵敏度的定义为y=T!Tda/a其中a是参数，r是系统的传递函数。经典的PID控
22、制器的传递函数为Ge(S)=Kp+Klsi+KDS由于它具有较强的鲁棒性，能够在大范围内适应不同的工作条件，同时具有简单易用的优点，因此得到了广泛的应用。在下文中，我们将基于系统根轨迹，采用系统灵敏度来度量控制系统的鲁棒性。(四)鲁棒双闭环PID控制器的设计鲁棒PID控制器设计的基本任务是：确定控制器的结构和参数，以获得最佳系统性能。而在设计中我们关心的是行走小车的定位和重物的摆角控制，因此采用双闭环PID控制，取外环为位置环，内环为摆角环。内环的的设计有较强的跟随性能，可以使吊车在准确定位的同时，摆动也衰减至零，从而达到防摆的目的。为了提高系统的性能，考虑到对象为非线性不稳定系统，以及反馈校
23、正具有以下特点：(1)削弱系统中非线性特性等不希望有的特性的影响；(2)降低系统对参数变化的敏感性；(3)抑制扰动；(4)减小系统的时间常数。所以，对于系统内、外环拟采用反馈校正控制。综上所述，设计出控制系统的结构如图所示。吊乍系统图12吊车防摆控制系统结构图Xr(S)给定1.内环(摆角)控制器设计假设所采用的伺服电机的机电时间常数较小，可以将其等效为比例环节。设70=50kg,Ks=2SNV(电机环节)，重物质量机与绳长/在不同的情况下可以变化，它们的标称值分别取?=5依，1=1*所以，内环系统未校正时的传递函数为(5)Ks280.56U(S)tnls1+(am0ni)g50/+539S2+
24、10.78确定控制器的形式对于内环反馈控制器HI(S)可以有PD,PI和PID三种可能的结构形式，绘制各种控制器结构下的“系统根轨迹”，加以比较并从中选出一种比较适合的控制器结构。图3-3为各种控制器结构下的根轨迹。各种控制器的开环传函的传递函数分别为：0.56K.52+10.29RootLocus-1-0.500.5RealAxis(seconds,1)图13比例(P)控制器下的根轨迹0.56*(KP+5PI:S-S2+10.78RootLocus-1.5-1-0.500.5RealAxis(seconds1)14比例积分(PI)控制器下的根轨迹0.56*(K+KS)PD:7IrS2+107
25、8RootLocusRealAxis (seconds-1)图15比例微分(PD)控制器下的根轨迹0.56YKP+KQ+%PID:LS2+10.78二SPUoOeS) sKSKF+(7+砸；为阻尼系数，=LKSKD2-Jf/IKSKP+(叫+廊)g控制器参数的鲁棒性设计为了保障系统控制具有较好的鲁棒性，即对于绳长/和重物质量小的变化不敏感，需要对内环控制器的参数进行鲁棒性设计。由灵敏度公式知，当某个参数变化时，系统轨迹(如伯德图、根轨迹、奈奎斯特轨迹等等)变化较小，就是说系统对该参数灵敏度较低，即鲁棒性较强。令T=e*,=/可以得到系统对绳长/的灵敏度为MTT虹Ulm0fc2-KsKds-Ks
26、Kp+(m04-m)g2S?+2+显同理可以得到系统对负载质量机的灵敏度为ST=arnmm0fe2+(m)gmg/(m)s2+2zl5+0：为了使系统对参数变化有较低的灵敏度，一般要求在系统参数变化时系统轨迹变化不超过5%o那么在此条件下，确定系统固有参数(绳长和载荷质量)允许的变化范围。用公式表示两变量的鲁棒性设计要求，即为Ky05STX40.051z,m由系统对绳长/的灵敏度公式可以得到H=1M2占为了保持内环系统的快速响应并无超调，取=l,则按绳长的鲁棒性设计要求可以得到即绳长变化范围为10%,即0.9m到LInk由系统对载荷质量m的灵敏度公式可以得到T10x9.8/50S=；为保证内环
27、的跟随性能，使响应时间尽量短，转折角频率3应该选得较大；然而当“选得过大时，系统稳定性变差。为此，可以取3=IOmdk,则按载荷的鲁棒性设计要求可以得到颔5.1m即载荷的变化范围为510%,即0.82kg到30.5kg将&=1和Q=TOmd/S带入内环传递函数公式并取整得到KQ=36,Kp=159综上可知，当内环控制器取KD=36,KP=I59时，内环将具有抑制“绳长变化1佻（即0.9m到1.1m）,载荷变化510%（即0.82kg到30.5kg）”的能力。下面将用仿真实验检验这个结果2.外环(位置)控制器设计鉴于内环调节时间对于外环来说较小，为简化外环系统的设计，可以将内环等效成一个增益为K
28、O的比例环节，则由前述内容可以知道，KK_SKKKP+2828159+(50+5)9.8=0.0056由内环传递函数公式可知这种近似应满足条件：S1+2g7S优由内环设计知=1,fl=0rad/s,S为外环响应频率范围，可以取为外环的剪切角频率,。为满足上式，不妨取“5倍系数”，则有计算得至UCQ095rad/S,即外环的剪切角频率不超过0.95%/S外环的简化设计注意到由摆角到位移的传递函数：迎！二吃冬。该传递函数分子没有一次项，这样的系统容易不稳定。为了便于设计，需要对该环节进行简化。由线性系统的性质，可以将该环节分解为两个并联的环节：比例环节/和积分环节g/。分别对这两个环节进行控制，所
29、得的结果与直接控制环节时是等价的。即可以将对图设计而得到的控制器参数，直接用于图所示的系统中。图17外环控制器设计系统框图图18简化后外环控制器设计系统框图对于二次积分环节gM,本身有两个不稳定的零极点，采用PD控制器能够将此环节校正到稳定状态。为了消除闭环零点对系统动态性能的影响，将控制器放在反馈通道。外环控制器的鲁棒性设计设所采用的反馈调节器的传递函数为2G)=l+为了调节前向通道的增益，起到快速、准确定位的作用，在前向通道内设置一个比例环节调节器G(s)=K图中，二次积分环节闭环传递函数为T=见母-一2JCr(s)/+gJ+2+武其中：庆=KlKOg比例环节的闭环传递函数为:TG)=由于
30、I(S)对系统的影响较小，为了简化设计，可以将石作为系统传递函数进行设计，而把THS)当作系统叠加的扰动进行处理。为保持系统始终稳定，达到较好的鲁棒性，应满足：T1()l+()对频域内所有角频率都成立。为了满足剪切角频率和吊车定位无超调的要求，选定=1,Q=a.6radls,则根据式显=KlKOg、=KiWK函以及验证公式3+ ln(4-9)c可以解得Kl=10.3,K2=2.7,经过优化选择，最终选定参数M=12,K2=2.5o所以有0O672Z,0,16+(l+0,0672Z)2二产+.5s+0.5625为了验证外环设计过程中的简化条件，可以做出不同绳长时HO)和l+7()的幅值图如图3-
31、6所示。上面为l+()伯德图的幅值特性，下面的三条曲线分别为I=1.1,1.0,0.9机时怩()伯德图的幅值特性，这四条曲线在任何频率是都在l+70)的下面，即幅值在任何频率是都小于上者。所以设想成立，整个外环系统可以用心($)代替。在MATLAB中编写如下画图程序来绘制伯德图的幅值特性:胸mPMIOT100101102Frequency(rad/sec)图19不同绳长时的幅频特性曲线num=11.51.125;den=l1.50.5625;numl=O.06048;denl=O.15121.06048;num2=0.0672;den2=0.1681.0672;nun3=0.07392;den
32、3=0.18481.07392;W=Iogspace(-1,2);mag,pha=bode(num,den,w);mag1,pha1=bode(num1,den1,w);mag2,pha2=bode(num2,den2,w);mag3,pha3=bode(num3,den3,w);magdB=20*logl0(mag);magdB1=20*IoglO(mag1);magdB2=20*logl0(mag2);magdB3=20*logl0(mag3);semiIogx(w,magdB);holdonsemiIogx(w,magdBl);holdonsemiIogx(w,magdB2);holdo
33、nSemilogx(w,magdB3);gridonxlabel(Frequency(radsec)ylabel(,MagnitudedB,)title(,外环简化条件验证BOdeDiagram,)接着检验内环的近似条件，画出外环单位负反馈开环传递函数的伯德图，如图，可以得到c=QASradst既满足剪切角频率又满足内环环节的简化条件。在MATLAB中编写如下画图程序来绘制伯德图：num=336O3292.8;Clen=50148849916585.6O;hd=tf(num,den);W=Iogspace(-1,2);Gm,Pm,Wcg,Wcp-margin(hd);mag,pha=bode(
34、num,den,w);magdB=20*logl0(mag);semilogx(w,magdB,r,);xlabel(Frequency(radsec)ylabel(,MagnitudedB,)title(,BodeDiagram,)gridonBodeDiagramOO - H Y Pm6ez-50-60-70IOT100101Frequency(rad sec)图20外环单位负反馈开环传递函数的伯德图ans二11.691269.00872.10380.4809此外，为使电机输出的最大控制力限定在一定的范围内，可在摆角给定和电机电压给定前分别加上饱和限幅环节，以限定电机的最大输出力矩。由于最
35、大输出量得到了限制，可使系统在暂态过程中的超调得以减小，因而稳定性也有所提高。四、 仿真实验（一）Simulink仿真实验根据上面的设计，可以建立如图3-9所示的龙门吊车控制系统的Simulink仿真程序。其中，系统初始状态为零，小车质量为50kg,小车的期望位置为4m载荷质量和绳长变化时仿真的结果分别如图3-10.11所示。图21系统仿真结构图上图中cranemo加1展开后如下:fdjIKuOTOrce内w5IIM提MIransrerrcni在MATLAB中编写如下画图程序来绘制m不同时的响应曲线：%m-0.5kgplot(time,position)holdonplot(time,angl
36、e)gridonxlabel(Time(s)ylabel(position(m),Angle(degree)axis(O12-68)title(质量(m)不同时的响应曲线(L=Im)holdon%m-l.Okgplot(time,position)plot(time,angle)holdon%m-5.Okgplot(time,position)plot(time,angle)holdon%m=20kgplot(time,position)plot(time,angle)holdon%m=30kgplot(time,position)plot(time,angle)holdon%m=50kgpl
37、ot(time,position)plot(time,angle)holdon%m=IOOkgplot(time,position)plot(time,angle)图22重物质量m不同时的响应曲线（L二Im）在MATLAB中编写如下画图程序来绘制L不同时的响应曲线:%L=0.5mplot(time,position)holdonplot(time,angle)gridonxlabel(Time(s)ylabel(position(m),Angle(degree),)axis(012-68)title(,L(m)不同时的响应曲线(m=5kg)holdon%L=0.9plot(time,posit
38、ion)plot(time,angle)holdon%L=1.0mplot(time,position)plot(time,angle)holdon%L=1.Implot(time,position)plot(time,angle)holdon%L=2mplot(time,position)plot(time,angle)图23不同绳长L时的响应曲线(m=5kg)(二)总结由仿真结果可见，摆角和位置在IOS左右时就可以达到稳定，同时，可以看出负载质量与绳长两个可变参数在“鲁棒设计”的参数变化范围之内时，系统控制的动态性能满足要求，即系统灵敏度变化不超过5%o在仿真实验过程中，也给出了参数变化超
39、出我们设计的参数变化范围的情况，从仿真结果可以看出，控制系统一样可以克服这些参数的大范围扰动，保证系统的控制性能。综上所述，采用“鲁棒双闭环PID控制”的吊车系统防摆控制方案可以在定位完成的同时，能基本消除行走过程中的摆动，实现防摆控制的目的。该控制方案的优点是系统具有较强的鲁棒性，随着系统参数，如绳长、载荷的变化，系统的性能变化较小(灵敏度较低)。五、 参考文献1高丙团，陈宏钧，张晓华.龙门吊车系统的动力学建模J.计算机仿真，2006,23(2):50-52.2刘为，瞿军，袁湛,等.龙门吊车防摇模糊控制研窕J.机械制造，2009,47(2):16-18.3郭晨海.龙门吊中的力学问题J.物理教学，2010(1):53-53.
