控制阀粘滞特性补偿控制方法研究

控制阀粘滞特性补偿控制方法研究

ID:35081977

大小:6.21 MB

页数:80页

时间:2019-03-17

上传者:U-24835
控制阀粘滞特性补偿控制方法研究_第1页
控制阀粘滞特性补偿控制方法研究_第2页
控制阀粘滞特性补偿控制方法研究_第3页
控制阀粘滞特性补偿控制方法研究_第4页
控制阀粘滞特性补偿控制方法研究_第5页
资源描述:

《控制阀粘滞特性补偿控制方法研究》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

'一'-‘-‘乂.^'^'、.、''.-、>^://''三、;:./^>:-1分'言巧-巧r/.’分心:"A--'--一.^''—二:/若^,/Z。,—V‘/'^#尤:々'--'-《''‘'iy.';兴.焉;^知,萨\的:、;、t分::.察:新;;—:-分类号:TP206单位代码;1的35'識.六;无申请号332047密级;21事斗‘:.萝抑?才、:議4嚇:碁養/^參,.:,;-:—心挺券.V、奉托>三户':义户.六,V.:硕去学位论文f'';;,.;.‘'如^:过琴礙.八.''..."-.:^;’7'人、.荀,V,龄冷..肾.余.;.龄.>矿T声細;;到ir%.’-厶‘:.治短齡r\乂盛於,乂'於'.'':‘知.,键邱护V:V'/■,c.中文论文題目:控制阀抱滞特佳补偿控制方法研究y英文论文题目:CompensationControlMethodsforStictionintheControlValvei*1.申请人姓名;.:;,王挺任.、’"';:一-.指导教师:谢磊副教授-—':?专业名称;控制科学与工程、、^研究方向:制阀粘滞特性补偿控制—擦......V、r蟲I糕.,所在院系.:控制科学与工程学院:...馬.、,V堂结J看"巧皆..■巧/扣克.装譜H师‘ ̄.论文提交日期^_〇二_方_全T_i十巧.0齐食刮)’。;朱恕劈J.当’'.:‘.'y、.y巧表-r;Ky一、梦父t^h:专.^.^'遍.獻心載德謂1’'.?.—、'';:,,,V:/皆LV斯於—//,,乂戶\、片尽 验制阀摇滞特化补偿验制方法巧宛戀论文作者签名:班也指导教师签名:今论文评阅人1:^2评阅人:^评阅人3;^评闽人4;慕亦军(浙江工业大学)答辩委员会主席;慕晋辉教摄(中国计量学院计量测试工巷学院)1委员:陈剑教授(浙江大学控制科学与工程学院)委员2;谢磊副教擾(浙江大学控制科学与王程学院)^^3:杨秦敏副教擾(浙江大学控制科学与工程学院)答辩日期;2016.03.09 CompensationControlMethodsforStictionintheControlValve戀\>’/\L/All化orssignature;^’Suervisorssina化re:pg—ExternalReviewers:AnonymousAnonymousAnonymousYiiunCai.ZhejiangUniversityofTechnologyExamininCommitteeChairerson:gpJinhuiCai.ChinaJiliangUniversityExamininCommitteeMembers:gJianChen.ZheiiangUniversityLeiXie,ZhejiangUniversityOinminYan泛Zheian泛University,iDateoforaldefence:March9.2016 浙狂大学研究生学位论文独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。除了文中特别加W标注和致谢的地方外,论文中不包含其他人已经发表或撰写过的研究成果,也不包含为获得浙江大学或其他教育化构的学位或证书而使用过的材料一。与我同王作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。学位论文作者签名:签字日期:>八年參月口日学位论文版权使用授权书本学位论文作者完全了解浙江大学有权保留并向国家有关部口或机构送交本论文的复印件和谣盘,允许论文被查阅和借阅。本人按权浙江大学可W将学位论文的全部或部分内容编入有关数据库进行检索和传播,可采用影印、缩印或扫描等复制手段保存、汇编学位论文。(保密的学位论文在解密后适用本提权书):学位论文作者签名i■側1导师签名:签字曰期:7^/《年夸月0曰签字曰期:年月/曰/| 浙江大学硕±学位论文摘要在工业过程控制系统中,气动控制阀是最常见的终端元件,其性能的好坏将直接与控制系统的调节品质相关。由于气动控制阀长期与被控介质直接接触,在工作一段时间后经常出现非线性特性,其中W姑滞特性最为常见。控制阀的枯滞特性容易诱发控制回路产生振荡现象,从而导致生产过程效率降低甚至生产停止。因此,在检测出控制阀存在接滞特性后,采取适当的补偿控制手段克服控制阀的裕滞作用十分必要。本文在阀读了大量关于控制阀摇滞特性补偿控制的文献基础上,从理论结合实践的角度出发,对控制阀粗滞特性补偿控制进行了较为深入的研究:。本文的主要研究内容如下1、分析了Knocker补偿法参数选取的依据,发现其脉冲幅值的选取相对困一难许多,,往往是通过人王试凑的方式然后提出了种基于黄金分割法的脉冲幅值在线自整定寻优算法,并对所提出的算法进行了仿真研究,验证了该算法的有效性。2、分析了两步补偿法的常规实现方式,发现其实现过程耗时较长,尤其是面对具有大迟延特性的过程系统一,然后提出了种两步补偿法的新型实现方式。所提出的新型两步补偿法无需等待过程系统稳定下来,提高了两步补偿法的应用效率,,并对该实现方式进行了仿真研究验证了该方法的有效性。3、在控制阀粘滞特性单参数模型MPC补偿控制的基础上,将更具代表性的双参数He模型转换为对应的输入约宋条件,利用MIQP求解器实现了控制阀姑滞特性的MPC补偿控制,并利用仿真平台将该方法的补偿效果与常规PID控制效果进行了对比,验证了该方法的有效性。关键词:气动控制阀;粘滞补偿;Knocker补偿法;两歩补偿法;MPC补偿法I 浙江大学硕女学位论文ABSTRACTIntheindustrialprocesscontrolsystem,pneumaticcontrolvalveisthemo巧commonterminalelementwhoseerformancewillbedirectlassociatedwkh化e,pyqualityofcontrolsystem.Thepneumaticcontrolvalvealwaysworkswiththecontroleeium.Soafteraeriodofmesomenonneaharactetcsmaldmdti,lircrisipyaearamongthemstictionisthemo巧commo打?ThestictionintheV过Ivewillcausepp,,oscUlationofcontrolloopseasUyandresiUtineficienclossandevenastandstillof,ytheproductiveprocess.Therefore,afterthedetectionofthestictioni打thevalve,itisnecessarytointroducea打appropriatecompensation化overcomethestiction.Basedttesthtsssothetnmethodfoonhereladreferenceeheifocusencompensaiorthe,stictioninthevalvefromtheoreticalandpracticalaspect.Themaincontentsasfollows:’1、The化lectionbasisoftheKnockercompensatorsparametersisanalyzed.化isfoundthattheselectionofulseamlitudeisrelativelmore出ficultandusuallppyyreliesonmanualerrorexperiment.Thenanoptimizationalgorithmofpulseamplitude-seton-rosesetawttlfuni打呂lineispodbasedonoldencio打l.Andasimulaionsudispgy'tintroducedtodemonstratetheroosedalorihmsusefulness.ppg2-、Theconventionalimlementationof化etwomovecomensationisanazedlypp,ichmaCO巧100cleeecatheaceestetwh:muhimsillinfofthroce化smwihlarey,pypyg巧-delaycharacteiiic.Thena口ewimlementationofthetwomovecomensationppwithoutwatinfortherocesssstemtoestableisroosedwhichcanmroveigybppip,p-theeficiencyofthealicationofthetwomovecome打说tion.Andasimulationpppstudyisintroduced化demonstratetheproposedimplementation、useftilness.3BasedonMPCcomensationcontrolforthesinlearametermodelofstiction^pgp?inthecontrolvalveHemodelofdoublearametersisconverted化thecorresondin,ppginputconstraints,andMIQPsolverisappliedtorealize化eMPCcompensationcontrolofstictioninthevalve*AndthecomensationefectiscomparedtothepconventionalFIDcontrolbysimulati打g,theeffectivene巧oftheMPCcompensationIII 浙江大学硕±学位论文controlisdemonstrated.Keywords:PneumaticControlValve;StictionCompensation;Knocker;Two-moveMPC;IV 浙江大学硕±学位论文目录摘要IABSTRACT虹@录V1绪论1111.引言1.2控制阀概述21.2.1控制阀简介2124.2.控制阀常见特性.361控制阀猎滞特性补偿控制方法研究现状11控制阀粗滯特性经典补偿方法6.3.127.3.控制阀粗滞特性改进补偿方法1.4本文主要肉容122基于黄金分巧法巧Knocker参数在线自塞定寻优算法1521.1引言521.2基于Kano横型的控制阀粗滞特性建模62.3Knocker补偿法特点分析192.3.1Knocker补偿法实现原理19232Knocke:20..r补偿法参数的选取2.3.3仿真研究及分析222.4基于黄金分割法的在线白整定寻优算法232.4.1黄金分割法的基本思想232.4.2Knocker关键参数a的寻优242.4.3仿真研究及分析262.5本章小结273控制闲藉滞持性开环两步补偿法的新型实现方式29312.引言9V 浙江大学硕±学位论文3.2基于He模型的控制阀摇滞特性建模303.3开环两步补偿法特点分析31313.3.1开环两步补偿法实现原理3.32.32常规开环两步补偿法的互次实现334.3.3仿真研究及分析3.4基于控制器输出期望值估计的新型开环两步补偿法%3.4.1新型开环两步补倦法的实现方法%3.4.2仿真研究及分析403.5本章小结414基于MIQP求解问題实现控制闲抱滞特性的MPC补偿控制43434.1引言4.2控制阀粘滞特性模型的约東条件转换444.2.1单参数d模型的约宋条件转换444.2.2双参数He摸型的约宋条件转换464.3MPC补偿法特点分析494.3.1MPC控制算法4942MPC50.3.补偿法的实观方法4.3.3求解MIQP问题524.3.4仿真研究及分析54445.本章小结75结東语巧5.1研究工作总结巧52研60.究工作展望参考文献63致巧67作者简介69作者在攻读巧壬学位巧间的科研成果71VI 浙江大学硕±学位论文1绪论1.1引言随着企业规模的不断扩大,现代工业控制系统中包含成千上万的控制回路,并且回路之间存在着精合作用,这些控制回路的性能表现将直接影响企业的产品品质、生产安全、生产成本、节能环保等指标,因此王业过程控制的重要性不言而喻。然而,由于工业控制系统日兹复杂,许多控制回路的日常维护工作并不及时都是出现了故障才去修理所L义控制回路性能表现不佳是十分普遍的现,往往,1一]象[26王业过回路进行的牡。项持续两年的调查结果显示,在对000个程控制,一的控制回路性能表现良好能分析中,剩余的控制回路的性能,大约仅有王分之一均有待提高。研究表明,如果现代工业控制系统中能有百分之的控制回路性能P]得到改善或者提高其能源利用率,那么将会挤出上亿美元的利润空间。34[]有研究指出,在工业过程控制系统中,那些性能表现较差的控制回路大部分都处于振荡状态。控制回路的振荡现象将直按导致企业的产品品质变差、生产一系列问题安全性降低、设备损耗加快、生产成本提高等。然而,在实际工业生产中,即使控制回路出现了振荡现象,但是因为相关设备并未到维修周期,是不可能停止生产过程而去对设备进行检修的。因此,如何在设备正常运行期间检测6[]控制回路振荡现象产生的原因是十分必要的。控制回路产生振荡的原因主要分为王种:控制器参数整定不当、过程存在噪声/外部干扰、过程设备存在问题。7-9"[]20% ̄/有数据显示,大约r30〇的控制回路振荡问题是由于控制阀出现了非线性特性所引起的,控制阀的非线性特性主要包括回滞、死区、跳变、据滞等,这些一特性均会对回路性能产生不同程度的影响。阀口的单位维修成本非常高,般在400美元至800美元,现代工业生产过程系统中大约有300万个阀口,如果毎个阀口都需要维修一,那无疑将是笔巨额开支。因此,利用过程控制手段来延长阀口的维修周期,保证在王业生产过程中阀口仍然能够正常工作显得栓外重要。在控制阀的各种非线性特性中,粗滞特性最为常见,并且控制阀抱滞特化非9][工业生产常容易诱发控制回路产生振荡。控制阀粘滞特性补偿控制不需要在周期内停止生产过卷,单纯依赖补偿信号或调节PID控制器即可消除或减弱控制1 浙江大学硕±学位论文回路的振荡现象一,是种十分有效的方式。控制阀枯滞特性的补偿控制作为工业过程系统性能改善的一个重要手段王工业生产过程的正常,对石油、化、发电等运转具有重要意又。1.2控制阀概述1.2.1控制闽简介控制阀作为工业过程控制系统中最常见的执行器,其应用质量的好坏将直接与控制系统的调节品质相关。控制阀作为过程控制中的终端元件,相比过去它得到了人们更多的重视。控制阀应用的好坏和许多因素相关,主要包括产品本身质量是否过关、操作人员安装是否得当、操作是否合理、维修是否及时等,同时也和选型息息相关。如果控制阀选型失误,可能会造成系统开开停停,严重的甚至无法投入使用,所W对于工业巧程控制系统的操作设计人员来说,应该充分认识到控制阀在现场的重要性,对其影响过程性能的相关因素要引起足够重视。控制阀的调节开度关系到被控介质的流量特性,这也就意味着控制阀要长期与被控介质直接接触,面对高温、高压、高粧度、强腐蚀等复杂环境,控制阀的选型要十分讲究。依据工业过程的特性,控制阀选型有几个原则:(1)按照工艺条件,选择合适的控制阀结构和椅料;(2)按照工艺对象的特点,选择控制阀的(3)流量特性;按照工艺操作参数,选择合适的控制阀口径尺寸;(4)按照王艺()选过程的要求,选择所需要的辅助装置;5择合涅的执行机构,执行机构的响应速度需要满足工芭要求。本文研究的控制阀类型主要是气动控制阀,这是因为气动控制阀在实际工业过程生产系统中的使用率高这90%W上。气动控制阀主要包括互个部分:阀体、执行器、阀口定位器。阀体是控制阀下部直接与管路相连接的部分,它通过阀抒与执行器栖连接,在执行器的推动下阀体可W改变其与管路间的流通面积,使得经过阀体的流量发生变化,进而使过程输出满足控制目标;执行器是用来驱动阀20 ̄抒运动的装置100KPa,然后,它可W接收的气压信号产生相应的推力,驱?动阀杆运动L乂将420mA的标准电信号转换成2(MOOKPa的标;阀口定位器可准气压信号,同时阀口定位器还能根据阀杆位移反馈信号使输入电信号与位移信2 浙江大学硕±学位论文-号保持良好的线性关系,在数字控制系统中经常使用。图11展示了常见的气动控制阀的工作原理,不同类型的阀口细节设计不同。当电信号U作为控制信号输入到阀口定位器时,定位器可W通过鼓膜的变化将电信号转换成弹黃为,使得进入到执行器上下两室的空气量发生变化,这样阀杆活塞两侧的受为随即发生变化,进而驱动阀杆发生移动。多数情况下,阀口定位器的尺寸保持不变,通常会改变执行器的规格来适应不同需求。执行器越大,其内部的气窒越大,这意味着在阀口定位器规格不变的情况下,要花费更长的时间才能改变执行器上下两室的空气压为差,这样阀杆移动就比较缓慢。反之,执行器越小,则可W更快化驱动阀巧移动。UActuatorL—-翊——_rIJV麵^ValveETl國-圈11气动控制阀工作原理图控制阀在理想工作环境中是呈线性特性的,但在实际作业中,控制阀大部分3 浙江大学硕±学位论文一时间处于较为恶劣的环境,当其工作段时间后,难免出现各种各样的故障,文11-17[]献对控制阀常见的几种故障进行了详细的分析和分类レ,归纳出乂下几种常见故障及其原因:(1)、。阀杆闽患能运动的部分所受阻力增大可能原因:填充原料过于紧、密,,使得阀杆与接触壁的摩擦力增大;阀杆、阀忠的同A度不好导致其受力不均衡;阀杆在长期使用后发生形变;气源含水在温度较低时容易在管内结冰,导致堵塞。(2)控制阀介质泄漏量大:,。可能原因阀杆长短不合适阀口调巧后无法关严;阀体内密封性变差。3一()阀杆仅在某方向运动时表现迟缓。可能原因:气动薄膜执行机构中膜片处破损泄漏;执行机构中密封泄漏。阀杆在往复运动时均表示迟缓。可能原因:阀体内有粘物堵塞;填充原料变质硬化。(4)控制阀在任何开度下都振动。可能原因:支摧不稳;附近有振动源;阀恕与衬套磨损严重。控制阀在接近全闭位置时振动。可能原因;控制阀选大了,常在小开度下使用;单座阀介质流向与关闭方向相反。(5)控制闽的线性特化不佳。可能原因:弹黄反馈性能变差;阀杆与接触壁的摩擦力过大等。122..控制阔常见特隹-2-控制阀常见的非线性特性主要包括死区特性(图1)、回滞特性(图13)、-4--死区加回滞特性(图1)、不录敏区特性(图15)、饱和特性(图16)W及藉滞特性-7)L(围1?乂上控制阀非线性特性图中0P表示控制器输出信号(ControllerOutput),MV表示控制阀输出信号(ManipulatedVariable),巧阀口位置。困-2所示的控制阀死区特性产生的1原因十分复杂,主要是机械元件反应时间和电信号反馈时间存在时间差问题,比如由于摩擦方面的影响,电信号输入了,但是机械运动要先克服静摩擦力再运动,这样就产生了死区。影响死区的因素有很多,但主要还是机械与电的反馈时间快暧的问题。1-3描述了控制阀回滞特牲单独出现时阀口的输入输出特性图,但回滞特性一一-与图12描述的死区特性经常同出现当两者,同出现时阀口输入输出特性如4 浙江大学硕±学位论文一图-41所示。回滞特性是控制阀输出信号总是滞后于输入信号的种特性。很多领域都存在回滞现象,比如介质中的电场或磁场需要逆向偏置的电压才可L乂完全退跑;机械系统中的齿隙也会造成回滞现象,比如变速箱中的齿轮咬合不紧时,正传和反转的机械特性称为反弹(backlash),与回滞特性相似。至1ndeadbanUdX7/W^I^(并OP--图12控制阀死区特性图13控制阀回滞特化iki\dea?>cflMnd>/7deadzonc!y//I)厂鮮I>>OPOP--图14控制阀死区加回滞特性图15控制阀不灵敏区特性一-图15所示的控制阀不灵敏区特性比较容易理解,就是阀口输入信号在某段范围区间时,阀口没有出现相应的线性相应动作。不灵敏区往往出现在零位附近,即输入信号很小时控制阀没有输出。-图16巧示的控制阀饱和特性是指阀口开度有上下限限制,阀口全关时开度00%为0%,阀口全开时开度为1,控制信号超出上下限时阀口无法响应。-7阀枯滞特性在所有非线性特性中最为常见也是本文研究图1所示的控制,一的重点,即阀段死区后即将滑动,。姑滞特性最大的特点是死区加跳变抒在经历5 浙江大学硕±学位论文但此时阀杆发生突然的跳变。A争''SI'广7A资deadbandyA了i---一心—庐屯1-?-i1^巧stckbandiI?I^OPm--图16控制阀饱和特性国17控制阀粘滞特性1.3控制阀摇滞特性补偿控制方法研究现状1.3.1控制阀粘滞特化经典补偿方法控制阀的巧滞特性会引起控制回路发生振燕,所W检测控制阀粗滞特性是否tW存在并对其进行量化分析是非常有意义的。关于控制阀描滞特性的检测方法19-23己经比较成熟,文献[都对粘滞检测进行了相关研究,但是关于控制阀姑滞]特性补偿控制的研究却相对较少。目前巴提出的控制阀粘滞特性补偿控制方法大口4-27Ps-口2-33]1W概可W分为吉类ocker补偿法、控制器调节法。:Kn、两步补偿法K一证nocker补偿法通过叠加系列的短脉冲方波信号到控制器输出信号,保一[24心]控制阀不停运动,避免控制阀停滞在个位置,进而保证被控量来不及发生较大的变化。Knocker补偿信号有王个参数:脉冲幅值、化冲宽度、脉冲周期,口8]这王个参数的选取与过程系统、控制阀枯滞程度相关。利用Knocker补偿法可W快速带离控制阀离开当前描滞位置,当控制系统达到稳态时,如果系统输出与设定值的误差满足控制要求,即可关闭PID控制器,避免阀口发生额外的移动。Knocker补偿法可1^乂有效地消除控制回路的振荡现象,不过是W增加阀口的振动频率,减少阀口的使用寿命为代价。两步补偿法首先在与控制器输出信号相同方向上叠加一个能量较大的补偿一信号,将阀口带离当前指滞位置,然后再反向叠加个补偿信号,将阀口直接移6 浙江大学硕±学位论文29】[动到目标位置。控制阀经过两步移动后到巧目标位置,所W称为两步补偿法。两步补偿法在实现时需要一个重要参数,即阀口到达目标位置所对应的控制器输一出期望值。两步补偿法的般实现方式是驱动阀口进行多次移动(除两步移动外),过程中需要等待系统输出稳定PW,进而获取控制器输出期望值。大多数情况下,两步补偿法需要将回路切换到开环控制,因此在补偿期间不允许设定值变化或过程存在外部干扰。控制器调节法可W保持控制回路处于自动模式,仅根据回路振荡产生的原因PH)来调节控制器参数,,从而消除振荡或减小振荡的幅度和频率。然而为了克服因控制阀枯滞特性所引起的振荡而去调节PID参数,容気影响常规操作下过程输出跟踪设定值的效果或回路抑制干扰的能力,比如去除控制器的巧分作用,通常会造成比较严重的稳态误差。132..控制巧粗滞特性改进补偿方法理想的控制阀粘滞特性补偿器应巧具备L乂下条件;(1)控制阀阀杆动作不应太频繁;(2)过程输出信号波动不应太大;(3)控制信号(控制阀输入信号)不一Kno应太大。许多文献在上节介绍的经典补偿方法的基础上,对cker补偿法和两步补偿法进行了相应的改进。(1)Knocker朴偿法的改进方法a.引入阔值条件虽然Knocker补偿法消除回路振荡的能为很强,但是会导致阀口频繁动作,增加阀口磨损程度工程中一L。在实际,被按量是允许有定的稳态误差的,所乂可Lcke一定阐值的补偿条乂考虑在允许的稳态误差范围内,适当给Knor补偿器设置34[1件当控制阀输入信号变化量Ut超出设定的阔值0时才开后Knocker,例如k()U,一Knoレ补偿器,cker,这样在定程度上可乂适当减小阀口动作的频率改进后的补偿信号表达式为:'asi-W—nt^t+片+T苗[c的|c的pfc=.1M的-1fc()0,otherwise1 浙江大学硕±学位论文b.控制冻结算法一当Knocker补偿控制达到目标效果]^乂后回路中存在的些微弱突变信号有,可能影响控制阀的移动情况,巧使控制阀阀杆移动的距离很微小,但这也是不必27一要的。因此,文献[提出了种控制冻结算法,即当控制回路的稳态误差很小]时,停止PI控制器和Knocker补偿器;如果控制回路的稳态误差较大,则重启PI控制器和Knocker补偿器。基于Knocker补偿法的控制冻结算法的具体实现方式;当Knocker补偿控制介入作用后,如果控制回路的稳态误差小于5%并且保持10秒W上,那么PI,控制器则停止输出控制信号,进而Knocker补偿器的输出也保持为零,这样便冻结了控制信号;如果稳态误差超过5%,并且保持2秒W上,那么则重启PI控制器和Knocker补偿器。此方法化同样适用于PI控制器单独作用。Knocker补偿控制冻结算法在消除回骆振荡W后,面对微弱的突变信号仍然有良好的鲁棒性,即使设定值频繁变化或者过程存在较大干扰,此改进算法也能有效处理。Knocker补偿控制冻结算法的控制回路框图如图-8-91所示,算法流程图如图1所示。Compensator>errori■EnablePIDDisturbance斯u ̄‘的TPIPG产hL-—Enab?le的,的APSP|V||Ut11■■…一———c()?Process?fr=?Va■?lvej^^^邮f戶。ITNoise图-Kno18cker补偿控制冻结算法的控制回路框图8 浙江大学硕dr学位论文Start〔)r,Enable化ePulses一IyesDisablethePIDcontrollerMaximumerroris^2戸<^C<^2during4tJ?^>^^noEnab*lethePIDcontroller圈-9Knocker补偿控制冻结算法流程图1(2)两步补偿法的改进方法a.P瓜控制器与两步法补偿器切换文献口巧提出了一种两步补偿法的改进方法,在该方法中,补偿信号不再叠加到PID控制器输出信号中,控制阀输入信号将在P瓜控制器输出信号和两步1-1法补偿器输出信号之间巧换,如图0所示。首先假设PID控制器巧使在回路因控制阀粗滞特性产生振荡后,仍然能够处理干扰信号和设定值变化,并且当补偿器介入控制巧,PID控制器能够跟踪控制L阀输入信号乂及偏差信号。如果补偿器作用后回路振荡消失,那么控制偏差趋于9 浙江大学硕±学位论文口5P]ID控制零,并且器输出信号化将趋于补偿器输出信号。此时如果过程中,存在干扰信号或者设定值发生变化,补偿器的控制效果会受到影响,误差信号不ID再趋于零,P控制器输出信号也不再趋于补偿器输出信号。那么,将控制阀输入信号切换到PID控制器输出信号,利用PID控制器抑制干扰和跟踪的能力,使控制回路再次处于振荡,然后再利用两步补偿器进行补偿控制。Comensatoron/ofifp——rPIDmodulewithtrackinsinal■gg置^—IwComensator(0iopi|ipypv*1Va一Process-lveet;^()1y ̄■t王i?TrackingsalIgnNoi化L—-图110改进后的两步补偿控制回路框图此改进方法在两步补偿开环实现时,只需要驱动阀口移动四步即可到这目标一12位置,从而估算出控制器输出期望值。式(.给出了四步移动的全过程:第步,)驱动阀口离开当前粗滞位置;第二步,驱动控制阀阀杆反向运动;第H步,在第二步阀杆移动的方向上,驱动阀杆继续运动;第四步,将控制阀带到目标位置并稳定下来。其中,第H步是为了获取过程增益,因为阀抒移动的方向与第二步相同,并且第3步的移动不会受到死区的影响。'<t<t+T的;imUt-<2571(0:0:+Tt含t+T(i)苗1)2;mlm(U=<12/、/、(.f的)t+-別2T<t<T?10:+t+3:T叫(im)巧〇1)3;mimSP-Vt-Ut2T—-—(+)sinat>t+37mr(;^iifii)i/C、pb.闲环两步补偿法两步补偿法一般在实现时需要将控制回路切换到开环控制,而在实际过程控制系,由于,,、统中操作环境较为复杂控制变量要求较为严格很多时候无法将回路切换到开环控制,因此文献[3U提出了闭环两步补偿法来克服控制阀粘滞作用。10 浙江大学硕上学位论文一组短时间的矩形波补偿信号闭环两步补偿法通过在设定值信号上叠加,从而实现与两步补偿法中相同效果的两步移动(控制阀)。其中,矩形波补偿信号参数的设计要非常系统化:基于控制阀粗滞He模型、分析不同阶段的信号、推一导关系式、设计最优化参数。闭环两步补偿法相比般的开环两步补偿法在处理建模误差-、可测量干扰等鲁棒性更好。闭环两步补偿法的实现过程如图1H所不。了^CDt化)tg>【tActpt-图111闭环两步补偿法实现过程-图11-1(a)中虚线表示设定值信号,实线表示过程输出信号1H(b);图-中虚线表示控制阀输出信号示控制器输出信号。图111(b)中A、B、,实线表C、D、E、F所对应的时刻点有W下关系:二+^5t二+0c=<:+13r.()【E二t+o、t二+白p11 浙江大学硕±学位论文0T短其中,T是控制系统的采样周期,s,是过程对象的时间迟延时间矩形波的持续时间。1.4本文主要内容近些年,控制阀姑滞特性朴偿控制的研究受到越来越多国内外学者的重视。在检测出控制回路因控制阀枯滞特性而发生振荡现象后,如何在不中断生产过程的前提下,根据实际情况恰当选用控制阀枯滞特牲模型,并利用描滞补偿器消除回路的振荡现象显得格外重要。目前,国内外研究学者提出的控制阀投滞特性补偿控制方法还很少,所W本文在已有研究王作的基础上,阅读了与控制阀粗滞补偿相关的文献及书籍,,并分析现有控制阀粗滞特化补偿控制方法所存在的不足对控制阀粗滞特性补偿进行了较为全面和《入的研究。本文的主要研究工作内容如下:一章介绍了控制阀摇滞特性补偿控制研究的背景意义第。首先介绍了控制阀的概述信息,对气动控制阀的工作原理进行了简单说明,归纳出控制阀常见的故障及其原因。然后列举了控制阀常见的几种非线性特性,对枯滞特性进行了重点介绍。最后,主要,针对国内外提出的控制阀粗滞特性补偿控制方法进行了综述包括经典补偿控制和改进补偿控制。第二章围绕基于黄金分割法的Knocker补偿信号参数在线自整定寻化算法展开了讨论与研究。首先介绍了控制阀粗滞特性Knocker补偿法的实现原理,并基于控制阀楷滞特性Kano模型进行了仿真研究。然后针对Knocker补偿信号脉冲レ一幅值a难乂选取的问题,提出了种基于黄金分割法的在线自整定寻优算法,分析了黄金分割法在此问题上的优势。最后,对所提出的在线自整定寻优算法进行了仿真,验证了该算法的有效性。第王章控制阀粗滞特性开环两步补偿法研究。首先介绍了两步补偿法的核成一原理,对开环两步补偿法的般实现方式进行了研究,并基于控制阀粗滞特性He模型进行了仿真实验。然后针对开环两步补偿法常规实现方式所存在的不足,提出了一种基于控制器输出期望值的新型实现方式,并详细介绍了新方法的实现步寐。最后,对所提出的开环两步补偿法新型实现方式进行了仿真,通过对比常规开环两步补偿法,验证了新方法的有效柱。12 浙江大学硕±学位论文第四章围绕基于MIQP求解问题实现控制阀投滞特性的MPC补偿控制展开了研究。首先分析了控制阀藉滞特性数据驱动摸型转换为含混合整数的约宋条件问题,分别对控制阀摇滞特性单参数d棋型和双参数He模型进行了输入约束条件的转换。然后讨论了MPC算法在处理有约末问题时的实现过程,分析了MIQP求解器在求解混合整数二次规划问题的注意事项。最后,通过多种情况的仿真实例验证了控制阀粗滞特性模型转换为混合整数约束条件问题的有效性,并将MPC补偿控制效果与常规PID控制效果进行了对比。第五章结宋语。首先总结了本文对控制阀粘滞特性补偿控制方法的研究。然后分析了如何将研究成果更好地应用于实际工业生产控制回路中。最后,提出了一若干有待进步研究的控制阀拱滞特性补偿控制方法。31 浙江大学硕±学位论文41 浙江大学硕±学位论文2基于黄金分割法的Knocker参数在续自整定寻优算法2.1引言通过在控制器输出信号后方叠加补偿信号来消除或减弱因控制阀粗滞特性所引起的回路振荡是目前最常用的手段,而Knocker补偿法可能是目前此类补偿[^Kno一方法中克服控制阀楷滞作用效果最好的方法。cker补偿信号是系列的短脉冲方波信号,其中包括云个参数:脉冲幅值a、脉冲宽度T、脉冲周期hfc,这互一个参数的逸取般没有明确的依据,但是它们的大小决定了控制阀粗滞特性补偿效果的好坏其中又Lレ乂脉冲幅值a最难确定。,义往的经验都是根据理论分析确定脉冲幅值a的初始范围,再通过大量实验反馈,得到脉冲幅值a的最优值,但是这样的参数选取方法效率很低,化很难系统化。一种经验整定方法黄金分割法是,其寻优过程收敛速度较快,能够迅速准确地找到最佳整定参数,目、。在控制器参数自整定寻优中标函数初值和步长的选取都会对结果产生影响。传统的优化算法有很多种,例如随即搜索、函数逼近法、区间消去法和爬山法等,新型的优化算法如PSO算法、混纯优化算法等。但考一虑到计算简单,,无需全局最优解,同时搜索区间为维黄金分割优化算法则更为适用。一维黄金分割法的基础上针对Knocke本章在,r补偿信号脉冲幅值a的选取LIt提出了相应的优化方窠,X误差平方和积分货正(negralSquareError)作为分割搜索时的目标函数,进而实现Knocker补偿信号脉冲幅值a的自整定寻优。这样减少了反复实验的繁琐,使得补偿控制更加自动化,大幅度地提高了对控制阀投滞特住的补偿效率。本章首先分析了Knocker补偿法的特点及其实现原理,针对脉冲幅值a的选取原则进行了分析,并利用仿真实验对Knocker补偿法的效果进行了验证。然后提出了基于黄金分割法的Knocker参数在线自整定寻优算法,在确定脉冲幅值a的初始取值区间后,可W快速搜索到最化值,并通过仿真实验验证了该方法的有效性。15 浙江大学硕±学位论文2.2基于Kano摸型的控制巧摇滞特任建摸控制阀姑滞持性建模也就是建立控制器输出信号(OP)与控制阀输出信号-133639,]MV[()之间的关系,主要分为机理建棋和数据驱动建横两种。机理建模要求大量的物理参数,,对于控制阀来说需要知道阀抒质量、弹黄系数、静摩擦为一、滑动摩擦力等,而且不同阀口的参数不同,参数也存在定的不确定因素,因此机理建模的工作量将非常大,而且不容焉实现。数据驱动建模则只需要相关数据,所需参数也较少,计算量不大,容湯在计算机上进行模拟和仿真,同时能够很好地逼近机理巧型,因此被广泛采用。、双参数控制阀粗滞特性的数据驱动模型有许多种,可W分为单参数模型P,37,39,",42lKP一模型。ano模型叫乍为种双参数数据驱动模型能够很好地描述控制阀的猎滞特性,包括死区、姑滞区、跳变区、滑动区。对于不确定信号或随机信号,Kano模型的处理能为也比较强,其双参数S、/分别表示控制阀枯滞特性曲线中的死区加摇滞区。、跳变区ValvePos扣onElasticForce()iiS'y?^//'??(h)(g)阳^,?b()!1?/〇?/??乙 ̄^?hControllerOutputAirPressure()-图21双参数控制阀祐滞特性输入输出特性16 浙江大学硕±学位论文2-图1中的虚线表示气动控制阀弹黄为与空气压为处于平衡状态时的特性曲线,也就是控制阀没有出现粗滞特性(理想状态)时的状态,此时控制信号输入到控制阀后阀杆位置会沿着虚线移动。当控制阀出现了釉滞特性,控制阀阀杆将)沿着实线移动,不,(a点表示弹黃为与空气压力再次处于平衡状态过由于粗滞作用,即使控制信号增加,阀口仍将停留在原来位置。(b)点阀巧开始移动,此时外部作用力克服了静摩擦力怎,摩擦力从静摩擦为后变成动摩擦力么,则跳变区/的大小为=-J2.1fsfv()阀杆从(b)点直接跳到(C)点。然后阀杆会沿着线k移动,k偏离中央虚线的横向距离为&),这是因为弹糞为与空气压为相差大为瓜。当阀杆到巧(d点时,此,直到空气压力增加到,弹黃为与空气压力差值小于&时阀杆再次停滞,e)/时,才会使弹黄力与空气压力差值等于/〇阀巧重新跳动到(点并继续沿,着线!2移动。当控制信号驱动阀杆移动到(f)点时,空气压为开始减小,此时阀杆准备改变它的移动方向,但由于枯滞作用阀杆会停滞,直到弹黃力与空气压力差值超过)h)&,阀口才会从化点反向跳变務动,然后沿着线移动。(点与(b)点对称,因此空气压为从(f)点到(h)点相差的大小为=S+f.2fs〇口)上述控制阀粗滞特性可^2-2U义用流程图描述,表示控制器输出信号,即阀。口输入信号,y表示阀口位置信号17 浙江大学硕±学位论文ControllerOutputu(t)u阳=10〇心?<>!00>;>;||^]了no—=ut0<'()Tno—?t=—1Auutut()()()^——<0AND1)W=-s;l)stp=1noIIfnoI**"y***"*—^-—<t=■>dd■r5^il()1^5]noI,f ̄nodyt=ut--()()口/)L5s=0=-巧y(ty的UIValvePositiony(t)图2-2控制阀牆滞特性Kano模型流程图18 浙江大学硕±学位论文2.3Knocker朴偿法特点分析2.3.1Knocker朴偿法实现原蕃?^一[1?1班在2002年提出了种有效的控制阀枯滞特性补偿方法该方法麵封,将一系列的短脉冲方波信号壹加到控制器输出信号,利用补偿信号的能量来克服""一Knocke控制阀的粗滞作用。由于这些短脉冲方波信号像敲击样,因此称为r补偿信号。Knocker补偿信号有王个参数:脉冲幅值a、脉冲宽度T和脉冲周期hfc,如图2-3所示。a??hk图2-3Knocker补偿信号Knocker补偿法应用于控制回路的结构如图2-4所示。1 ̄?Knockerr(t\)图2*4加入Knocker补偿器的控制回路框图SP表示设定值,PV表示过程输出;控制器输出信号为Uc(t),Knocker补偿=U器输出信号为Ufc(t),阀口输入信号为It,则w(t)U的+fct。其中Knocker的f()Utレ补偿器输出信号fc(的大小及方向可乂用下述分段函数表示:)au-<si巧nc(〇Uttt+h+Tc(,fc(p))p\M=fc.3的口)〇t>t+h+T、,pk一UUPID其中,t表示上个脉冲开始时间,fc〇的方向由的变化率决定。控p(c(〇19 浙江大学硕±学位论文制器输出信号在叠加Knocker补偿信号之后的变化曲线如图2-5所示。Controlsignalii"Ac户/Ia/-"I—I*t图2-5PID控制信号叠加Knocker补偿信号后的变化曲线Knocke一r补偿器产生的毎个短脉冲方波,其能量大小及方向应该恰好补偿一掉控制阀枯滞作用。如果个短脉冲方波的能量过小无法克服枯滞作用,那么,一阀口会继续停滞在当前位置,补偿无效,远;如果个短脉冲方波的能量过大远超过了艳滞作用的大小,那么阀口移动的距离会超过预期位置,发生过补偿。因此,在诊断出阀口出现粗滞特性并且化滞程度可知的情况下,需要适当地设计Knocke一r补偿信号的参数,保证每个短脉冲方波的能量恰好可W补偿掉控制阀粗滞作用,并且不出现过补偿现象。关于控制阀姑滞特性的诊断W及粗滞程度的检测在许《文献中邑给出相应方法,因此并不作为本文研究的重点,本文涉及的控制阀粗滞特性补偿控制方法均在粗滞已知的情况下进行的。Knocker补偿法的优势很明显,它并不影响PID控制器参数的设计,对控制回路过程对象的特性也没有任何要求,需要的参数信息量很小,算法的思想化十分简单,易于实现。23.2Knocker.补偿法参数巧选取Knocker补偿信号的能量由脉冲幅值CU脉冲宽度T和脉冲周期hfc王个参数决20 浙江大学硕±学位论文定,如果它们选取的不恰当,则会影响到控制阀摧滞特社的补偿效果,也就无法消除或减弱因控制阀摧滞特性所引起的控制回路振荡现象。通常Knocker补偿器要求毎个短脉冲方波的能量足够小,因为在阀杆滑动后,控制回路测量信号不能立即返回给PID控制器,这主要是由于测量变送环节有一ce定的响应时间,因此Knokr补偿器接收到的输入信号(PID输出信号)存在""微小迟延特性,从而导致有错误的短化冲方波信号继续在补偿。如果单位脉冲信号的能量较大,控制阀可能会接着克服新位置的粗滞作用继续移动。所W,""当Knocker补偿信号的单位脉冲方波能量很小时,错误的补偿信号叠加到控制信号后则不足(^驱动阀口离开新的描滞位置,但是单位脉冲方波的能量化不能太小,至少保证阀口能离开当前猎滞位置。Knodcer补偿信号Wfc的和过程输出;Kt)之间的关《可心乂用式口.4)表示y=U24k.命()其中,G表示过程对象的传递函数,Gc表示PID控制器的传递函数,则过程输p出可L乂用W下式子表示y=2.5)(式(2.5)表明过程输出与0T成正比,而且aT决定了Knocker补偿信号中单位短脉冲方波的能量。一考虑在阀杆运动的瞬间,不能叠加过多能量的补偿信号,因此脉冲宽度TL般要相对较短。考虑到实际情况,T的值应该大于采样时间t通常可乂巧置为采样时间/I的整数倍,同时也需要根据阀口大小来进行相应调整。Knocker补偿信号的脉冲周期h最好和控制回路闲环时间常数相关,这样可レ乂保证两次有效脉冲fc补偿信号之间的间隔不会发生太多变化。最简单的方法就是令脉冲周期hfc成为脉一<n<心冲宽度T的函数=riTn为正整数25。,例如hfc,般可W议置由此可乂,看出,脉冲周期hfc和脉冲宽度T的选择还是化较容焉的,但是脉冲幅值a的选巧就要相对困难许多。在hk和T确定W后,a的值会直接影响单位脉冲补偿信号的能量,21 浙江大学硕止学位论文a选得过一如果大,阀口的开度就会很大,这样执行器侧就很难控制。所Wa的一 ̄取值要相对较小,般范围在1%4%,但是也没有具体的理论依据来确定最合适的a值。另外,在a的初始范围确定レ乂后,需要通过反复的实验数据来手动找工作量会很一到最合适的a值,这个大,因为人工试凑的方式存在定的不确定因素。2.3.3仿真研究及分析KnockerATLAB仿真2-4所示控制阀粘滞特性补偿控制M模块框图如图。设定值SP选用幅度为1的阶跃信号,控制阀工作点为0.5,rCO是白噪声干扰,-能量为0.010.005(21)KanoKa,噪声增益为;控制阀粘滞特性图使用模型,no5=化2=化15,/被控对象的传递函数为模型中参数设置为;G.6切口)0.5s+1PI控制器P=0.8/=0.3控制器采用常规,其中参数设置为,Knocker补偿器;T选h的整=选为脉冲宽度T的的脉冲宽度为系统采样时间数倍,T2周期h;脉冲fc3倍=6过反复实验调整冲幅值a的最佳值为化2。仿真时长设置为,;经,脉200秒,在120秒时加入Knocker补偿控制,控制回路姑滞特性补僖前后控制效果如图2-6所示。12I11iIIIII.[1-0-.8\r\Kno=ckera1201ttsf严_r^0-.4^0--.21--SP-PV-020406080100120140160180200奇Timess〔)图2-6控制回路加入Knocker补偿前后效果图2-6LK从图中可乂看出,控制回路在引入nocker补偿控制前,回路因阀口姑22 浙江大学硕止学位论文滞作用而处于振荡状态;在引入Knocker补偿控制后,回路振荡明显消除,过程2-7输出能稳定在设定值附近,控制效果明显好转。阀口的振动情况如图所示。玉 ̄IiIiIIII1[0--.8-0.'.4广S厂0-一.2fI--(V0--02II._■--0.〇-〇-8〇20406080100120140160180200Timess()图2-7控制回路加入Knocker补偿前后阀口振动情况2-7中从图可W看出,当控制回路引入Knocker补偿控制后,阀口的振动颠率明显加快。阀口的这种频繁动作能够使得被控量来不及发生较大幅度的变化,从而消除了控制回路的振荡现象,使得控制回路的稳态误差减小。然而,振动频率的加快显然会导致阀口磨损如剧,从而加速缩短阀口的使用寿命。:.r2.4基于黄金分割法的在线自整定寻优算法L4.1黄金分劑法的基本思巧黄金分割法又称为0.618法,是寻优方法中的经典方法,其算法简单、效果显著一,是许多寻优算法的基础。黄金分割法的本质是种序贯消去法,按照固定比例逐步缩小搜索区间:。黄金分割法的基本思想如下6和区间缩减的相对精度5需要先给定=①初始区间[a,],再确定循环次数n加+1’并依据取点么式获得试探点A、X2的值,即(;。譜。)3-V5— ̄-x=a-bax+()()Vs-i=+— ̄-X口2()作句223 浙江大学硕i学位论文==②计算函数值/l/0:1)和/2〇2)。/<=,则转④⑤,⑨比较/l和&。如果A>/2;如果/i/2,则转;如果/l/2则转⑥。_=乂=X=乂=a_===a66,+6X,’,nn1,④令1,i心2i/i/2/2/沁)=0n*0如果n,则转⑦,;如果则转⑤。====_==打=n_ab乂乂XXa+占乂,乂1,⑤令a,2,2i’i2/2A,/i/(1),n=0n*如果,则转⑦;如果0则转⑤。,===—2—乂=—'@令〇;1£:,6;5:,义0+3>/^6£1,〇+/5121(()/)()2(〇-===n-=12)(&a),//如),/批2),^1,如果n0,则转⑦;如果n单0,)/i/2则转⑨。=2⑦获得极小值点X(a+的/。传统的优化算法有很多,例如随即搜索法、函数逼近法、区间消去法和爬山法等。新型的优化算法化日渐成熟,例如PSO算法、混純优化算法等。但是考一虑到计算简单,不需要全局最优解,同时搜索区间为维的,那么黄金分割法则更为适用。2.4.2Knocker关巧参教a的寻优本节采用黄金分割法来实现Knocker补偿信号关键参数a的在线自整定寻优。为了反映搜索到的脉冲幅值a的控制品质,F,这里选用两个优化目标西数O)和e的。V(x)表示误差目标函数ZSE(IntegralSquareError),即误差平方和的巧分:k+p2=-2.7制化旬)()三^=ifc+l式(2.7)中,fc为开始计算时间,P表示巧£所对应的时间区间,Yi为系统输出,托P为设定值。e(t)表示直接反映控制系统调节品质的目标西数,即被控量偏离设定值的大小:e=-2(t)y.8iysp()24 浙江大学硕±学位论文在Knocker补偿信号的另外两个参数确定レ乂后(脉冲宽度T和脉冲周期/i),fc脉冲幅值a的大小便决定了Knocker补偿信号单位脉冲方波能量的大小。如果能量过大一,会导致阀口开度过大执行器,,侧很难再控制控制回路继续振荡,并且超调量会较大,发生过补偿现象。如果能量太小,则可能无法克服阀口粗滞作用,阀杆位置则不会发生变化,即补偿失效。因此,性能评化函数F〇c)在a的搜索区间大致呈现西曲线状一,这样便可W采用黄金分割法记录毎步的a值,同时将Knocker补偿信号叠加到PID控制器输出信号中^e,通过观察1(对和的变化(〇情况,根据控制要求设置相应的闽值,最终得到可行解。由于在实际工业过程中,控制回路很可能存在噪声干扰,脉冲幅值a可能存在多个极值点区域,6,所W在a趋于稳定后,搜索停止观察目标函数(0是否在阀值范围内。如果不满足相应控制要求,则改变初始区间,重新搜索。本节针对Knocker补偿信号,从单参数变量寻优角度出发,利用黄金分割法在线获取脉冲幅值a的最佳值,具体算法步骤归纳如下:①确定脉冲幅值a的初始区间[的,更新周期r,分割次数n,变量区间,〇2]。的性能评价函数KO),eW及阔值条件0。(〇=口382〇-=②根据取点么式计算试验点和X+0.〇X+2,则1(21),2%0618—.(^2口。1)-<③分别计算相邻两个周期T内的误差平方和和F0:aa,如果2)|5打=nln=(,则+,若p,转⑦,否则转④。④如果7的>1^XV<V,则为较佳点:)沁)2,转獄如果知)〇2),则为较佳点,转⑥。=X==-==⑤令ai,XXXaX77aXii2,2+a2i,〇1)〇2),取,将得i2到的Knocker补偿信号叠加到PID控制器输出信号中,计算该周期下的e(〇.若e(t全部小于0,转⑨否则转⑨。),=X==-==⑧令£l2,XXXaClXVyaX22l,ii+22,〇2)〇l),取i,将得到的Knocker补偿信号叠加到PID控制器输出信号中,计算毎时刻下的e(t).若et全部小于()Q,,转⑨否则转③。⑦停止分割,若e(〇>Q,缩小初始搜索区间[叫,口2],转②,否则转⑨。=⑨a(ai+a2)/2,停止搜索。25 浙江大学硕±学位论文2乂3仿真研究及分析在MATLAB平台上对本章所提出的基于黄金分割法的Knocker参数在线自2-4整定寻优算法进行仿真实验,控制回路框图如图所示。设定值选用幅值为0.2-ss=^eP的阶跃信号,被控对象的传遵函数为G(),控制器采用常规I控制^,P=1495/=化259K其中参数设置为.,,控制阀枯滞特性采用ano模型,粗滞=化2=5参数为5,/化1,过程输出存在白噪声干扰,嗓声能量为0.01,噪声增益为0.005,系统采样时间设置为0.01巧,仿真时间为2100秒。在没有加入Knocker补偿控制时,控制回路由于阀口枯滞作用会使得过程输一出在工作点附近产生明显的振荡-8700,如图2前部分所示。在仿真进行到秒PIK=2=时,在控制器后方加入nocker补偿控制,其脉冲宽度T,脉冲周期hfc6。一节中提出的寻优算法来在线整定脉冲幅值利用上a,黄金分割法的初始区间设置为阳.7入约宋为0<?<1Knocker补偿信号单位脉冲方波,0,控制系统输,]=10=0二=能量的更新周期r秒.01(30.015a〇,闽值条件Q,P,6。化冲幅值的初始值设为零,在t=700秒时开始搜索,900秒レ乂后a稳定在0.215附近,a的搜索过程曲线如图2-9所示。'''os.f-PVa6021022..,[]K=nockeratt700syyy0--.15001000巧002000Timess()图2-8Knocker补偿器黄金分割法寻优过程中的回路补偿效果26 浙江大学硕±学位论文''''=0.6F:v_215f巧。.3尸I^一■""三…■■0.2T0-_.1I〇0500100015002000Timess[]图2-9Knocker补偿器化冲幅值a的搜索过程曲线2-8和图2-9加入Knocker补偿控制后可W看出也就是a从图,在补偿初期,在搜索初期,过程输出的超调量相对较大,对应的a值也较大。当a的拨索值趋于稳定后,其逐渐靠近人工试凑的最优值化22,并最终稳定在区间[0.210,0.220]内,此时,捏制回路的补偿效果也比较理想。仿真结果表巧,利用改进后的黄金分割法对Knocker补偿法参数a的在线寻优算法具有良好的逼近特性,并且收敛速度较快,可W极大程度地缩短手动调节时间,使得Knocker补偿控制迅速生效。25本.幸小结本章详细研究了控制阀枯滞特性Knocker补偿法,介绍了Knocker补偿控制的实现原理,并且基于Kano粗滞模型对Knocker补偿法进行了仿真实验,验证了该方法的有效性,同时分析了谅方法的局限性,即W增加阀口振动频率,缩短。KnockerT阀口使用寿命为代价补偿法有兰个参数,其中脉冲宽度和脉冲周期hfc的选取较为容易,但是脉冲幅值a的选取相对困难许多,也没有任何针对性的优化算法,因此本文提出了基于黄金分割法的脉冲幅值在线自整定寻优算法。黄金分割法相比其他经典/新型的优化算法计算量小,无需全局最优解,同村搜索区一间为维,对于Knocker补偿信号参数a的寻优更为适合。该方法只要先确定脉a的初始取值区间,冲幅值,并将此区间输入到本文所提出的寻优算法中通过设置相应的阔值条件,在线自动搜索出最适合的a值,无需反复实验人工试凑。无27 浙江大学硕±学位论文论是实验研究,还是实践应用,该方法都能大幅度地提高Knocker补偿器的接入效率,进而加快消除或减弱控制回路因阀口粗滞特化而引起的振动现象。最后通过仿真实验,展示了脉冲幅值a的搜索过程,验证了基于黄金分割法在线寻优的有效性。28 浙江大学硕±学位论文3控制阀撫滞特化开环两步补偿法的新型实现方式31.引言29[Knocker补偿法外]控制阀枯滞特性补偿方法中,除了,两步补偿法也比较经典。Knocker补偿法参数的设置与控制阀姑滞程度大小有关,但基本不依赖阀口粗滞特性建模的准确性,这也导致了其无法完全消除回路振荡,只能保证阀口不停留在接滞位置,而且增加了阀口振动的频率,导致阀口寿命缩短。两歩补偿法在消除或减弱控制回路中因阀口粗滞特性所导致的振荡问题有很好的效果。相比Knocker补偿法,两步补偿法避免了增加阀口的振动频率,化不需要调节PID控制器参数,所W在王业应用中两步补偿法得到了很大程度的重视。两步补偿法的核瓜要点是通过驱动阀口进行两步相应移动,从而到达目标位置一主要包括开环两步补偿法。然而,当前各种形式的两步朴偿法都有定弊端,28-%i[]P]和闭环两步补偿法。开环两步补偿法在实现上所需时间较长,准确地说,开环两歩补偿法需要兰步开环阶跃响应,其中两步是用来估算控制器输出期望值的,也就是让阀口到达目标位置。由于很多实际工业过程具有大迟延特性,因此开环两步补偿法在实现时就需要很长的响应时间。闭环两步补偿法需要准确获取过程的动态特化W及阀口的姑滞参数等信息,才能保证阀口能移动到目标位置,l8[]尽管很多辨识方法可レ乂使用,但是过程动态特性W及阀口枯滞特性信息的准确性也无法保证,甚至有些无法获取。一本章主要对开环两步补偿法进行研究,并提出了种新型实现方式,保证开环两歩补偿法的实现过程较快.同时在实现过程中对回路的负面影响较小,并且没有任何中间假设条件,使得开环两步补偿法的补偿性能真正得到提高。此方法主要获取开环两步补偿法中一个重要参数,即使得阀口到达目标位置所需控制器输出信号期望值。而估算该参数所需时间并不像开环两步补偿法常规实现方式那样长,这样使得阀口能迅速到送目标位置,从而让过程输出更快、更准确地跟踪设定值。本章首先分析了传统开环两步补偿法的特点レ乂及常见的兰次实现形式,并对该方法进行了仿真实验。然后提出了基于控制器输出期望值估计的新型开环两步巧 浙江大学硕±学位论文补偿法,介绍了此方法实现原理及步聚,通过仿真实验验证了此方法的有效性。3.2基于He棋型巧捏制巧拖滞特植建摸42f一H]e樓型也是种双参数粘滞模型,它能很好地描述阀口的粗滞行为。相化Kano模型,He摸型结构更加简单,没有那么复杂,应用起来十分方便,在实际工业过程中或实验研究分析中经常被采用。Inutpr的sure"的pmj—wt—t—cuiU.+()u(1j[)]yr,rr、-—?=-U=usi打curnMU£:1v的公(_s)dp()似//的=sn—.=cumjiUig(jcumjiffs)fa3-图1控制阀枯滞特性双参数He模型流程图-图31是控制阀粗滞He模型的流程图,其中尼和/d分别表示阀口处于粗滞状5==-态时静态摩擦力大小和动态摩擦力大小,+/,。中间变量IV是/sd//s&,,cumu表示作用在阀口上全作用在阀口上的残余力但是它无法驱动阀口移动_。cxmiu足够大部为的总和如果_,并且能够克服静态摩擦力尤,阀口将会移动,移动的距离与控制器输出信号U和动态摩擦力/d相关。反之,阀口位置不会变化,cumu将继续保留到下一个控制周期_。H一e模型非常录活,而且也能适用于些特殊情况。例如,当控制阀只有滞=后行为,那么可W设置/;当控制阀只存在静摩擦力时,可W设置为零s/d/d;阀口位置的上下限可L乂分别设置为0%和100%。虽然He模型的逻揖较为简单,但这为实验研究提供很多便捷性,尤其在控制阀接滞特性不明显时。30 浙江大学硕±学位论文基于He模型流程图中的描述,本文对He模型的输入输出行为进行仿真分析,1=化=化2。输入信号为正弦波,幅值为角频率为1,粗滞参数设置为5,,/s/dHe摸型的输入输出仿真图如图3-2所示。1111|0--.5§?<c巧^£〇-^(I>^--0.5^0o!s:.Co打trollerOutput-困3-2He模型输入输出仿真困PMV)(O3.3开环两步抑偿法特点分析3.3.1开环两步补偿法实现原理一第章中介绍了现有补偿控制阀枯滞特性的方法中WKnocker补偿法、,其PID控制器调节法和两步补偿法最具代表性。Knocker补偿法的原理是通过叠加一系列的短脉冲方波信号到控制器输出信号中,补偿信号的能量恰好可W克服控制阀指滞作用一er,保证阀口不断的移动,从而避免阀口停滞在个位置。Knock补偿信号参数的选取在第二章中进行了详细的介绍。虽然Knocker补偿法可W有效并迅速地解决控制回路振荡问题,但却W增加阀口振动频率,减少阀口使用寿Knock一命为代价,这是er补偿法非常大的个弊端。PID控制器调节法根据振荡产生的原因来调节控制器参数(比例、积分),,但只能减弱振荡的幅度或频率很多情况下无法完全消除振荡,而且PID控制器调节法可能在补偿控制阀粗滞特性时对回路跟踪设定值或克服干扰产生负面影响,甚至导致更大的稳态误差。レ两步朴偿法主要有开环两步补偿法和闭环两步补偿法,其中又乂开环两步补31 浙江大学硕±学位论文偿法最为常见且应用最广。开环两步补偿法在实现时需要将闭环控制回路断开,回一使路处于开环控制。开环两步补偿法的实现原理主要分为两步:(1)第步,通过向控制器输出信号中叠加足够大的补偿信号,补偿信号方向与控制器输出信一二号方向致,驱动控制阀移动并离开当前枯滞位置,;(2)第步阀口离开当前二二摇滞位置后,再通过第步补偿信号将控制器反向驱动到目标位置,第步补偿信号方向与控制器输出信号方向相反。这样,驱动控制阀进行了两步移动,阀口最终到达一目标位置并稳定下来。开环两步补偿法在实现过程中需要个重要参数(控制阀到达目标位置时所需控制器输出),此参数需要通过驱动阀口产生额外的移动才能获取。=3*3.2常規开环两步抑偿法的次实现上一节中提到开环两步补偿法实现过程中需要驱动控制阀产生额外的移动来获取m,下面将介绍开环两步补偿法的常规实现方式。,,3-3rt6mtUwt闭环控制回路如图所示,其中()、的、()、的、Ym的和()分别代表设定值、过程输出和噪声/干、反馈误差、控制器输出、控制阀阀杆位置C一扰。是过程输出的静态调节量,是个常数。/〇表示存在枯滞特性的控制阀,采用双参数He摸型。W的r的_e(0m的v(〇的1iI—??C??-^〇(s)?fO)G(s)Qnrc3-3闭环控制回路框图)图(控制阀存在枯滞特性LT-I(neartimeinvat)Gs过程对象采用线性时不变lirian,表巧式如下:() ̄- ̄31TT.()nuT〇S+:(,i〇32 浙江大学硕±学位论文1一4]レ[式(3.リ可乂简化为阶纯迟延模型:=岭)口.2如)控制器采用PI控制:C=1+33而().切()^当控制回路发生了振荡行为,过程输出无法跟踪设定值,通过检测确定是控制阀粘滞引起的则可W将回路切到开环控制。假设设定值rt不变,回路没有,()=0较大的外部干扰w(〇。为了佑算出控制器输出信号期望值的大小,需要口1]组移动中阀口要移动两次先驱动控制阀产生额外两组移动,两组移动间,每一隔时间内要保证回路这到稳态。驱动控制阀移动的方法与两步补偿法实现原理致(驱动控,即在控制器输出信号方向上叠加能量较大的补偿信号脉冲方波身,一制阀离开当前粘滞位置,到达个新位置,再在与控制器输出信号相反的方向上叠加能量较大的补偿信号一,反向驱动控制阀移动到另个新位置。控制阀毎经过一过程输出分别记录为和组移动后到运的稳定位置分别记录为m和m,Y。i22结合He枯滞模型流程围描述的特点,如果驱动阀口产生移动时阀口输入信号先1<增大后减小,那么可乂得到W下关系:K7=(jigg++cT3.4p/(j)jj()其中K表示过程增益,表示回路到这稳态后的过程输出。而m和m,和Ypi22都代表稳态值:,那么可W得到W下两个式子Krri+f+c=35(yi.pid)()K+f+=Onzd)cy23.6p()结合式3:.5和式(3.5就能够得到K和()),/dp.心=,吃37.)(—yc\=--^/d13.巧^(此时,过程增益反被求出,阀口粗滞参数店可^选用上值,通常情况可|^^通过捡933 浙江大学硕±学位论文测手段获取阀口粗滞参数,同时静态调节量C又己知,那么就可W计算出控制器输出期望值TTI:mT-Cssfm=—39ssf正.^()经过上述方法,两步补偿法实现过程中的重要参数m。,便可W获取,于是可W驱动控制阀执行最关键的两步移动,也就是第王次开环阶跃响应,直接将阀口移动到目标位置v。实现方式为在控制器输出信号上叠加第云組补倦信号,先?一个新位置将阀口移动到,避免阀口处于摇滞位置无法移动,再将阀口输入信号レ大小调节为并保持不变,过程输出便可乂准确跟踪没定值。然而,考虑到实脉工业过趕中常规两步补媛法实现的时间会很长,尤其是过程对象有较大迟延特回路在较长一性时,控制段时间内的性能表现会很差。33.3仿其研究及分析为了更加直观地了解常规两步补偿法的实现过程,本文通过MATLAB平台过S-imulink平台搭建如图33对该控制阀补偿法进行仿真研究,通所示的巧环控制回路。过程对象传递函数为38163-25^=0156.46S+1控制器采用常规PI控制器C=02513.11切.(+)()^=—1921设定值r(t)为常数信号=38c.341/干,r,,,静态调节量,噪声扰信号wte==暂不考虑,8.4,,控制阀粗滞模型采用H模型枯滞参数设置为尼/d()3.5243,仿真采样周期为0.5秒,仿真时间为7000秒。回路在上述控制条件下发生了振荡现象,在2500抄时将回路切到开环控制,一引入开环两步补偿法。此时,控制器输出处于上升趋势,因此先叠加组正向的0一,方波补偿信号,幅值为2持续100秒,将阀口从当前位置移动到个新位置却;5一(m)再将补偿信号幅值降为阀口输出信号,反向驱动阀口移动到另个新位i00一置V,补偿信号能量保持不变秒时叠加,保证过程输出达到稳态的。在%I一组正向的方波补偿信号,将阀口从稳态位置移动到,幅值为22,持续100秒34 浙江大学硕±学位论文个新位置切;再将补偿信号幅值降为4(阀口输出信号m2),反向驱动阀口移动一个新位置V祈到男,补偿信号能量保持不变2,直到过程输出这到稳态72。将1、饼3.7,.、的、於代入到式中,得到过程增益馬再将馬代入式(39中2()),得到控制器输出期望值m。当m已知则可1^3ssss,乂量加常规两步补偿法中的第次补偿一信号,经过。,将阀口直接驱动到目标位置1;?段时间后过捏输出跟踪设定值。¥-阀口输入信号与阀口位置信号如图34所示出信号与设定值信号如图3-5,过程输所示。8011111|75-_口■70-■-165-A-,……zk…驻二W為—tmi-巧\7口=50--IIjII1000200030004000500060007000Timess〔)困3-4阀口输入信号)(实线与阀口位置信号(虚线)^:301JJ111000200030004000500060007000Timess〔)图3-5过程输出信号(实线)与设定值信号(虛线)35 浙江大学硕±学位论文-34和3-5王次图图清晰地展示了常视开环两步补偿法的实现过程,其中第、レ实现是两步补偿法的核A,,可乂将阀口直接拉到目标位置振荡现象消失。有几,个关键点需要注意:首先毎次驱动阀口离开当前位置需要足够大的输入信号,保证阀口肯定能克服枯滞作用;其次,义须要估算出使得阀口到达目标位置所需的控制信号大小最后,此方法需要将控制回路切换到开环控制。另外,由于补偿期问控制回路处于开环控制,如果设定值发生变化或者过程存在较大嗓声/外部干扰,控制回路将无法响应这些变化,过程也将无法跟踪设定值,这也是一30通过监常规开环两步补偿法的个弊端。不过这个问题臣经被解决了,文献[]=-值来判断回路是否需要巧换回闭环控制测控制偏差e的:r识Ym的的幅,再利用PI控制器的作用使过程输出重新跟踪设定值或抑制嗓声/外部干扰。3.4基于控制器输出巧望值佑计的新型开环巧步补偿法3乂1新型开环两步补偿法的实现方法两步补偿法的关键是耍获取控制器输出期望值ruM,但是常规两步补偿法获取该值需要花费很长时间,而且在获取过程中容易对控制回路产生更坏的负面影响。尽管在文献[18]中提出了利用哈默斯坦系统辨识方法来解决这个问题,但该方法的复杂度非常高,如何降低复杂度来解决控制阀粗滞问题是非常重要的。另外,两步补偿法只需要参数,而系统辨识方法能提供很多信息,但这些信息一在两步补偿法中却利用不到,因此冗余度化过高。基于上述问题,本文提出了种新型的开环两步补偿法,旨在更简单迅速地获取参数m,s,同时在获取过程中尽可能小地影响控制回路性能。通过分析常规开环两步补偿法的实现过程,可W发现驱动阀口移动所需的控制信号要足够大,但过大的控制信号又会导致阀口移动距离较大,从而使得控制回坡的巧能变差,所W找到能合理驱动阀口移动的控制信号显得非常重要。从控制回路历史振荡数据分析来看,当控制信号处于峰值或峰谷附近时,阀口位置都会发生变化,因此本文提出的新型两步补偿法在驱动阀口产生移动时选用历史报一荡数据中控制器输出信号的极值。该方法在回路切换到开环控制后,输出系列一的控制信号使得回路产生个完整周期的振荡信号,通过分析此振荡周期内各个36 浙江大学硕±学位论文变量间的关系,计算出阀口目标位置Vss的值,进而得到控制器输出期望值m,,。整个过程不需要等待过程输出达到稳态值。预先定义几个变量:历史振荡数据中控制器输出信号的最大值为m,最小值i一一一为m,2,个振荡周期内(假设控制信号有两个峰谷个峰值)第个峰谷到峰值的时间为了1,峰值再到第二个峰谷的时间为。。新型开环两步补偿法具体实现步骤如下:①当控制器输出信号处于上升趋势并在峰谷附近时,叠加补偿信号使得阀口输入信号大小为mi,驱动阀口离开当前猎滞位置。②阀口在新位置保持较短时间后,立即叠加反向补偿信号使得阀口输入信号降为mj,反向驱动阀口移动,控制信号保持m2不变,持续时间。一③再次叠加与控制器输出信号相同反向的补偿信号,使得阀口输入信号大,驱动阀口离开步骤②所在位置,小为m,控制信号保持mi不变持续rz时间。i④重复步骤②中的动作,将阀口输入信号再次降为m,保持r时间不变。2i⑤经过上述操作一过程输,手动驱使控制回路产生个周期的振荡信号,出在此期间出现两个峰值一,个峰谷,分别记为Yd、Y2和Yi。通过分析此振藻周期内各个变量间的关系,结合控制阀枯滞He模型的特点,可计算出阀口目标位置的值,进而估算出控制器输出期望值ms,。一组补偿信号驱动阀口离开步骤④所在位置⑥叠加最后,然后立刻将阀口,输入信号补偿到步骤⑤中估算出的值,最终将阀口移动到目标位置。新型开环两步补偿法的实现过程重点在于手动控制过程回路产生至少一个周期的振荡信号,保证在此期间阀口只在两个位置移动,进而利用各个变量之间的关系计算出控制器输出斯望值m?。针对步骤⑥中的推导过程,下面将进行详细分析。((过程输出y〇、、(。由于过程对象有、Yi、Y2对应的时间点分别记录为012一一—迟延特化白t—0(0,因此在防,n位置保持不变,记为V,在的,)i巧阀i2)一阀口处于另个位置,记为V2。公式(3.巧中的过程对象Gs对y(t求导后可写成W下形式:())(i)=--y似0)3.12()hh37 浙江大学硕±学位论文选用tc)时刻作为y(0的觀始时刻,对上式进行拉普粒斯变换-t=-Ys31o)().3W()^^可レ乂得到KVtpiW〇)r:-()IrSs+s+p(是)是314(.)口t-yK1(〇)p化馬=^+Ss+王Tptett因此[〇,i),—t1〇.15=-ytKV-+tKVe而口)()pLCy(〇)p^)同理,tGt,t时[i2)—tt\!了.-口巧yt=KV+tKi;e。()p2(y(i)2)p-7=£-£=1531又因为,7££,则式(3.和式(.\1〇221)句分别可^写成T.17-p口)t=KV+tKi?e;y(i)pi(y(〇)pi)_扛T3.18t=Kt—()V+KT?ep><2)pi(y(i)pi)3151=/ic将式(.和式3.TV+)(巧,sspss代入W下等式rt〇+了1-T-(XcdT+ssy())I-T-化=0c3.19Iy())()JtjL式(3.19)可W利用控制器C(s)(式(3.3))经过拉普拉斯变换的微分方程求得,最终写成W下形式:38 浙江大学硕±学位论文KU ̄=O+了KOhT牛。了pss\2)pl22)3.20/_立_丘(),、、f、(、r'T_Kl_pt__prtVe+7(/CVlep〇y(〇)pi)pbi)p2)|^j^jt=t1731又因为y(y,再结合式(3.和式.则式3.20等号右侧为零〇)(2))(巧,()-+=KUT+了7VT032s2)zz).psCi馬机\(UK不为零,那么p-=U(了的T+"了03.22ssTl+2)l22)()考虑He模型的变量关系v=mt+3.23i(〇)fa()=t-3Vm.242(i)fa()根据式(3.23)和式(3.24),那么式(3心)可(得到^--7了 ̄mt+mt(〇巧(i巧2i。u=f"知T1+T2了1+了23.25()饥了—+呵八了212。2_=_—t江T+TT+T1212由于饥、:^和那么阀口目标位置1^12/^都已知,,,便能够计算出来。如果控制信号m(X)先增大后减小,根据He枯滞模型控制信号期望值771?为—3771"+.26"ss尼()t反之,如果控制信号m(先滅小后增大,那么ni为)is=-3mVd.27ssssf()此时一,新型开环两步补偿法只需最后次两步实现。一レ3步骤⑥中驱动阀口移动到新位置的控制信号般可乂设置为式.2(巧的值t=-m(比h)心±a(mim2)(3.2巧jwc39 浙江大学硕±学位论文a一t其中为调节系数,保证阀口可W克服姑滞作用w为最后;sw次两步实现开fi=-0始时刻,(s(+了。witch22注意,本文所提出的新型开环两步补偿法如果在实现时控制信号处于下降趋势,上述方法仍然适用,只需修改部分么式中的正负号即可,本文不再做额外介一绍,只K个方向为例。3.2.4仿其研究及分析MATLAB仿真3-338新型开环两步补偿法模块框图如图所示。设定值为,为了与常规开环两步补偿法进行对比,过程对象、控制阀粗滞程度、控制器参数3一工业过程中等均与.3.3节中致,另外考虑实际必定存在嗓声干扰,仿真中在过程输出引入能量为化01的白噪声千扰声增益为0.013,10000,嚷仿真时间为秒。通过检测控制回路历史振荡数据发现=,控制信号最大值m66.7i,最小值一一一m=53.3,,在个振荡周期内(控制信号有两个峰谷个峰值),第个峰2=二==谷到峰值的时间:T307抄:T788。ti,峰值再到第个峰谷的时间2秒在4000秒时,将控制回路巧换到开环控制,并利用本文所提出的开环两步补偿法-对控制阀糖滞特性进行补偿63-7所示,控制效果如图3和围。I![IJIII70--50-45--'1111111L-4.0200030004000目000600070008000900010000^Timet切图3-6阀口输入信号)(实线与阀口位置信号(虛线)40 浙江大学硕±学位论文尸afIIIIII!I1-50V-im4Hy〇yz25--[I1II]1j2000300040005000600070008000900010000Timets](图3-7过程输出信号(实线)与设定值信号(虛线)-5-7对比图3和图3,可W看出本文所提出的新型开环两步补偿法的实现过25004000程大约需要秒,而常规开环两步补偿法的实现过程大约需要秒。考虑新型开环两步补偿法在实现过程中无需等待过程输出稳定下来,因此可W适当缩短模拟振荡的时间,补偿过程可L乂更快。另外,在补偿过程中,新型开环两步补偿法没有对控制回路造成更坏的影响,而常规开环两步补偿法使得过程输出远远偏离设定值。3.5本章小绪本章详细研究了控制阀粗滞特性开环两步补偿法,介绍了常规开环两步补偿法的实现原理,并基于控制阀粗滞特性双参数He模型进行了仿真实验。通过分析常规开环两步补偿法的优缺点一,提出了种新型的开环两歩补偿法该方法可,W缩短对控制阀粗滞特性补偿的时间,并且在补偿过程中几乎不对回路产生额外的负面影响。本文所提出的新型开环两步补偿法在变量关系推导过程中需要保证控制阀只在两个位置变化,因此在控制回路切到开环控制后,需要驱动控制回路产生至少一个周期的振荡信号,这样可W保证阀口位置可控,从而满足上述条件。然后基于控制阀枯滞特性双参数He模型的特点,同时利用拉普拉斯变换获得此振荡周期内各个变量间的关系,进而快速计算出阀口目标位置,这样便可L乂得到两步41 浙江大学硕±学位论文补偿法中最重要的参数m,,。该方法栖比常规开环两歩补偿法最大的优势是补偿过程中无需等待过程输出稳定下来。42 浙江大学硕±学位论文4基于MIQP求解问题实现控制阀撫滞特性的MPC补偿控制4.1引言一模型预测控制(MPC)在王业生产中是种被广运应用的先进控制算法。模fA45型预测控制在很多工业领域中都十分常见]工,例如石油化领域、造纸领域、天然气管道领域、食物加王领域、汽车制造领域、航空航天领域W及冶金领域。[^模型预测控制之所W这样流行,是因为它可W很好地处理有约束问题包括,输入约宋条件问题、状态约宋条件问题W及输出约東条件问题。此外,模型预测,预测未来时域被控过程输出信号的行为控制可W利用己知的内在棋型,在线计[^算出当前时刻最优的控制器输出信号。然而,模型预测控制设计方法的潜在限制是它在处理不确定模型问题时显得为不从知,这主要是因为过程对象本身建模时含有不确定因素,或者传感器动态特性不确定等,又或者执行机构无法准确建模,例如控制阀。基于混合整数二次规划问题(MIQP)来设计模型预测控制器可W用于非线性47t3特性的补偿控制,例如本文研究的控制阀枯滞特性补偿。通过分析控制阀摇L滞特性不同模型的特点,并结合模型预测控制算法框架的特点,可乂针对性地进行控制阀摇滞特性补偿说计L。因此利用模型预測控制可乂有效地解决因控制阀死区特性和跳变特性所导致的输出信号难L乂跟踪问题。在多种控制阀艳滞特楼模型中,本文将选取有代表性的控制阀摇滞特性单参数d模型和双参数He模型来分别进行研究。另外,由于过程变量的多样性将使得操作情况很难满足系统的约宋条件,因此会导致过程输出有过多或不必耍的瀑移。通过理解并研究控制阀梢滞特性的行为特点,将有助于消除或减小其对控制回路的负面影响。口D]本章主要利用MIQP求解问题来实现MPC控制器的最优设计,从而对控制阀粗滞特性的非线性进行补偿。首先分析了控制阀粗滞特性单参数d模型的特点,通过引入逻辑状态变量将该摸型转换为输入约宋条件问题。然后将此方法扩展到更加常见的控制阀枯滞特性双参数He模型中,同样将其转换为可利用的输43 浙江大学硕±学位论文入约東条件问题。之后针对双参数He摸型进行了MPC编程,由于过程对象摸型中包含前一时刻的条件变量一,因此具有定的特殊性。最后通过棋拟仿真验证IP求解问题的MPC控制器设计补偿方法的有效性P了这种基于MQ,并与常规ID控制效果进行了对比,同时讨论了此方法的适用范围。4.2控制阀抱滞特性棋型的约宋条件推换4.2.1单參《d巧型巧巧宋条件转换--d棋型的图41描述了控制阀中执行机构的摧滞特性单参数输入输出行为。输入信号Ut在任何时刻都可能是正值或者负值,即输入信号大小可能增大或者i()减小レ乂如下定义:果ut是正值,阀,这里将其作为横坐标。这两种情况可0)如i()抒移动的正向距离用d表示;Ui)如果Ui(t)是负值,阀巧移动的反向距离用d。表p示。因此d和d的表达式如下:。p—t=-dp(+dU)j的41.()^-d=+-uKtn似(巧)771占控制阀执行机构中的藉滞特性因此可W表这成下述形式;'-么叫>t-m冲d1e[叫的的p()--―U=‘mAu<dt142b私的(的]:似n().的,()pt—U1:otherwise^i)p47此模型可W通过引入[]3个逻辑状态变量来分别表示3个动态区域而输入,的变化量AUt与粗滞死区d的大小关系可W决定阀杆输出的位置。另外需要注|()U=Ut-l时意的是当片)(),说明输入信号处于死区里面。基于MIQP求解问pp题的MPC控制器设计策略需要考虑输入约束。假设控制阀执行机构的非线化单元为整体元件,因此棋型可W写成:44 渐江大学硕-上学位论文"—ut:!()+du(t)>u(t_)——=:<uW<wtrft()u:L4的的().3p()——t1:i=ut心p()i(t)(0tU‘py?々4-图1控制阀粗滞特性d模型的输入输出图一組逻辑状态变量=采用上述简化模型,并引入&,;1,2,3,分别代表输入j变化的3种状态(输入变量为正值、负值或零值)。另外,针对输入变化量Au(fc)的这3个逻辑变量要满足W下关系:3=14.4公0()=1;其中=12mm表示多输入系统的1,,输入维数。MPC可从处理棋型有约宋问题,因此考虑引入逻辑变量来将控制阀粗滞特一组线性不等式約束条件性模型的分段特性餐换为,针对第i个输入信号的约宋条件则可W用如下命题表示(假设输出没有约宋限制);45 浙江大学硕±学位论文W<*^WfcUf()f,7nmf,爪似;<AUfc<AUi()m〇c打i;,S-*>AWfcd4.5ixi()()i--S?Aufe<di2()ii5。一Au佩=0A>-=即"爪0)代表fc010Aufc<001Aufc0。叫()曲,(,代表如(,,代表,)i())i()レ一沮混合整数线性不等式因此,自变量包括连续输入信号变量,上式可乂转换为mUeR和逻辑状态变量5£化化利用文献48中描述的方法可W进行如下定义:[] ̄— ̄^^ufe1imaximax(),,4.6()_="——1加ni加打)将4.4、4.5和4.,最终MPC优化问题的约宋条件可W写成下述巧()()(巧结合起来式:么<AUi(fe)<^Ui扣_max——l-d含Mdf(i加ni)(Ai)Sui— ̄(d-<i-iJ)i{max}i,1-(5ii(fc)+d^巧+6口)4.7馬乂ii咖()折If(fc)+过£>(百£min+df5。,)6>1-叫脚百6诚化(巧)Sw<-扫1S佩献()i柿巧4.2.2双参数化棋型的曲宋条件转换4-2e图是控制阀粗滞特化双参数H模型流程图。它有两个参数,分别是充和,示阀口粗滞静摩擦力的大小。它托怎表,记表示阀口枯滞滑动摩擦力的大小们一==一和5、的关系分别为5+。.充,般来说;/5记//;4尤耍大于记与控制阀描滞特性单参数d模型相似,双参数He模型也和输入信号的变化量相关。不过He模型又引入了中间变量u^,Uf是作用在阀口上的残余外为,它并没有使得阀口产生任何移动,即Uf小于阀口摇滞的静摩擦力。但是Wr的大小化46 浙江大学硕±学位论文He一将作为模型的判断条件之,与输入信号的变化量同时影响阀口的移动情况。〇一_071"是变量_个临时变量,表示阀口当前所受外为的总和它受平衡于阀,cwn口的静摩擦力/。如果_M足够大s,并且可W克服静摩擦为A,那么阀杆位置ut将发生变化tt(^(),u^()由控制器输出信号w()阀口输入信号)与阀口釉滞的动态摩擦力/d的偏移量决定。另外,如果cum_u不能克服静摩擦力尤,那么阀杆一将不会移动,将继续保留到下个控制周期中。Inutpressureutp(),[ ̄—t—tcumjxUr+wu1[()()]"。"IT1',U=-打-.=U-Utscumut1v(0()苗(_/s)()/dv=m—.=cwnjiUsincuurg(_fs)f凸图4-2控制阀粗滞特裡双参数He模型流程图控制阀指滞特性双参数He模型的数学表达式如下:cuml=w——jlit)+Cut1r([()()]'—=ut—si打cwnw(幻()0(_尤),cum_u>\\fs-u=nu■tsicum(r()g(_fs)fa4.(巧=t—U1)v的咕((cumu<,l_lfsj=、liv(t)cumji通过上式可W得到W下关系———=u(tu(t1)cumlUt)jr(.)49.()=cuml>一cuml--Wr(0f'jfsjUr(t)〉fsdfd47 浙江大学硕±学位论文-- ̄=->cumjlcumi<Ut<Wjr0fd./sr()fdfs(=cumicumcumi—t=0u<U.j,_j{)\\fsj而输入信号变化量为—Aut=t—tu(u(1)(4.10()))将(4.10)代入(4.9)中可W得到W下关系:cumu>一Au>-fs的/s/d<-4ownu一Aii<.11_/s的/d/s()-=cumu<*Aut0_\\fs()一=-假设dd利用上文中描述的方法,欲将此聚参数模型转换为a线/s/,’一?=123性不等式约束条件,同样引入组逻辑状态变量,,j,此处只考虑单输&,入过程对象系统。由于此模型的判断条件cumu分为3个区域,所W_y的上限化为3,即取3个逻辑状态变量来分别表示3个区域。满足下列关系:3=14.12公()=1;根据He模型的分段特化,并结合逻辑状态变量,有W下命题(假设输出没有约束限制);M一"k含Mm化()mojc—Au(fe)<AUaxffiS-*Au>4fcd.1i()(^一<—Aud么化)-=S?Aufc03()He一組线性不等式之前进行W下定义将模型转换为:=--1iU婦ufcW()414.()_=--M1馬nmmm)将(4.巧、(4.13)和(4.14)结合起来,自变量包括连续输入信号变量u6和逻辑状态变量5G{0,1},最终MPC优化问题的约宋条件可W写成下述巧式:48 浙江大学硕±学位论文Au-——d>征巧(16加71i)Au-< ̄d(UmaxAul—+d<巧+d饥似)(么)4.15()_Am+d>(巧mm+的&Aul->u(5^i?(3)Am<1-而而1化^()4.3MPC补偿法特点分析4.3.1MPC控渐算法1.预测模型49一3[MPC控制算法,也就是模型预测控制算法,,是种基于模型的控制算法此模型又称为预测模型。MPC可レ乂根据过程对象的历史信息和预测模型,预测过程未来的输出,从而找到最优的未来输入。MPC控制算法只关记过程对象模心、状态方程等棋型掏可型的准确性,对于模型的形式并不要求,例如传递函数乂作为预测模型。此外,无论是线性对象、非线性系统或分布参数系统,只要具有预测功能的信息集合L,化可乂利用预测模型来进行预测控制。预测模型可W描述过程系统未来的动态行为,能够找出未来时刻被控对象输出变化与没定值的偏差,然后W此为依据,调节未来时刻控制信号变化即可使得系统的动态特性得到改善,从而获得比常规PID控制更好的效果。2.滚动优化MPC控制算法的最大特征也就是滚动优化。模型预测控制算法通过未来时刻输出与设定值的偏差这一指标的最优来确定未来时刻控制信号的变化。指标的形式与系统未来行为相关,通常取未来输出与设定值偏差的方差最小。不过也可心乂采用更多形式,,不同形式对应不同的控制性能保证控制性能指标在给定范围内即可。与通常的离散最优控制算法不同的是,模型预测控制的优化并不是全局不变的优化目标,而是滚动式的有限时段的优化目标。毎个采样点的优化目#只一时刻采样点的化化目标又与下一时刻到与从该时刻到未来有限时长有关,而下49 浙江大学硕±学位论文未来有限时长有关,预测时长不变,则优化时刻同时向前推移。因此,模型预测一控制的优化并不是次性离线完成,而是不断在线滚动式优化。这种有限时段优化目标是有局限性的,即理想情况下只能得到全局的次优解,而无法得到最优解,但是这样的滚动形式却能顾及模型失配、时变、外在干扰等不确定性因素,迅速进行补偿,不断利用新的优化目标来与实际情况匹配,保证了实际控制效果的最优。这在复杂的工业过程系统中是十分必要的,因为这样的系统中模型失配、时变、外在干扰等因素无法避免,过程的不确定性会对全局最优解产生很大影响,因此通过在有限时段内建立滚动优化目标反而更加有效。3.反馈枚正一个环节反馈校正化是模型预测控制中非常重要的。预测模型只能作为过程对象动态特性的基础性描述,而实际系统中存在非线性、横型失配、时变、迟延、外在干扰等不确定性因素,因此基于不变模型的预测控制就不可能与实际系统完全匹配,这就需要对基础模型进行在线修正。反馈校正的引入可W很好的提高系统性能,增加滚动优化的优越性。这样,模型预测控制算法通过滚动优化确定了一一系列未来时域的控制信号L乂后,并不把这些未来时刻控制信号次性全部作用到控制器中而是只利用当前时刻的控制信号来预测下一时刻的过程输出,其他,未来时刻的控制信号只作为当前时刻优化出来的控制信号的参考量,从而避免过程系统的输出与设定值偏离过大一。而下时刻则先将实际过程系统的输出反馈给MPC控制器,利用实际输出信息再基于预测摸型进行有限时段的优化。总之,反馈校正力图在预测优化时对过程系统的动态行为做出较准确的判断,利用反馈信息构成巧环优化。4.32MPC.补偿法的实现方法一T对于般LI(线性时不变,L)对象可乂直接利用MATLABMPCToo化0X中的模型预测控制设计工具提供的图形用户界面进行分析一,按照菜单提示要求一步歩进行操作即可,能够迅速地设计出期望的模型预测控制器。但是对于模型具有非线化特性的对象,此方法则不再适用。而本文研究的控制阀粗滞特性摸型是非线性的,因此只能通过编写MPC控制器代码的方式,其中包括死区和跳变实现粗滞补偿控制。50 浙江大学硕±学位论文4-3MPC控制过程如图所示,设计器需要建立被控对象的预测模型,被控包含控制阀摧滞特性的非线性部分。因此,首先需要建立控制阀粗滞特性的非线性模型,本章中采用控制阀粗滞特性单参数d模型和双参数He棋型来进行研究。被控过程中需要引入中间变量U来作为MPC控制器输出信号U经过控制阀(描p滞)后的信号,化就是实际过程对象的输入信号。然后建立控制阀(描滞)和过PC一般需要进行离散化控制程对象的整体预测模型。由于M控制算法,也就是一设置采样时间,这就需要对建立的预测模型进行离敢化。模型离散化的方式般1SZZ有两种(被控过程为S域):()将被控过程从域蒋换到域,再从域转(2),换到时域;将被控过程从S域转换成状态空间方程再从状态空间方程转换到时域。不过MPC也能控制连续系统,但本文不对此进行展开研究。被控过程1IywufcU(fc)irl()1pMPC控制一控制阀 ̄?过程对象器(粗滞) ̄^nI?II图4-3被控过程包含粘滞特性的MPC补偿控制回路一LLMPC优化问题的目标函数般可乂写成X下形式:PM22=■-+R?-Aufc+14.16;式0[y从+0yre/从+0][(;)]()式==il71—式中,M,P表示预测步长表示控制步长,P般大于M;Q表示输出部分的控制一/?加权系数,表示输入部分的控制加权系数;输出部分偏差的平方和从下时刻一时刻开始计算,当前时刻输出已无法改变,可,截止到经过预测时长的未来某一W不计算在内,;为了使输出更好地跟踪设定值目标西数后部分的控制量与输一一入变化量相关,,并从输出前时刻开始计算截止到经过控制时长的未来某时刻;其中最关键的问题是求出y(fc+〇的表达形式,可W利用上文中建立好的被控对象预测模型来找到未来过程输出与控制输入的关系,从而确定目标亟数的最终形式。51 浙江大学硕±学位论文MPC控制算法能够很好的处理有约束问题,在目标函数确定后,根据约束条件并利用相应的求解器即可找到有限时段下的优化目标的次优解,从而获得当前时刻最佳的控制器输出信号。目标函数采用了常见的偏差平方的形式,因此是一二上文将控制阀祐滞模型转化成約束条件后含有二进制变个次规划问题,同时量{0,1},则整个问题变成混合整数二次规划问题(MIQP)。针对上述分析,MPC优化问题最终可レッ写成W下形式:Tmiri+占zzs.t.Cz+d<0417阿(.)z=UdJneZ巨ficdZe(d{Uf其中,Cz+d<0表示混合整数线性不等式约束(4.7)或(4.15),Zc表示输入信号u,?二?Z表示逻辑变量。上式将自变量U中引入进制变量,两者同时作为约束因dAj素。此时测模型、目标函数和约束条件都巴描述清楚,控制阀粗滞特性的补,预偕问题转换为求解MIQP问题,得到的最优解即可作为MPC控制器当前时刻的最佳输出信号,通过不断的滚动优化从而克服或减少因阀口牠滞特性所引起的输出偏移。4.3.3求解MIQP问巧MATLAB工具箱中有提供求解二次规划问题(QP)的方法,但是由于约束条件中含有{0,1}变量,因此需要利用OPTI工具包来求解此类问题。OPTI工具包中可W解决LP、MILP、QP、MIQP、MINLP等问题,本文需要利用其中MIQP求解功能来实现MPC的补偿控制。MIQP的标准形式如下:52 浙江大学硕±学位论文1了了min:xNx+Xxf5s.t.Ax<b斗北AX=b()eqeq<X<UlbuXEZie0勺{4}L上下限其中,自变量取值为整数,约宋条件包含线性不等式约束、等式约束义及或二进制{〇,1}。然而,目标函数(4.16)的形式与(4.1巧的标准形式不同,不能直接得到(4.1巧中所需的矩阵W和/,若通过人工转换的方式计算量较大,尤其当预测步长P和控I工IP求解器打制步长M取较大值时。但是OPT具包中的MQ必须获取矩阵和/,因此无法直接使用该方法求解本文所采用的目标函数优化问题。通过比较分析,OPTI工具包中所提供的混合整数非线性规划(MINLP)求41解器可W直接利用目标函数(.巧,同时约束条件中化可W加入非线性约束,其IPINLP他巧式完全相同。因此,束解MQ问题可^间接转换为求解M问题。MENLP的标准形式如下:minxxfi)s.t.Ax<b<x<U419l.bb()cx<d()CX=deg()巧X£Zi呵巨{(U}可从直接将(4.1巧代入上式中,非线性约東条件为空即可。53 浙江大学硕±学位论文4.3.4仿真研究义分析利用MATLAB平台对本章所描述的基于MIQP求解问题实现控制阀姑滞特性的MPC补偿控制进行仿真研究及分析。本文采用控制阀粗滞特性双参数He模型,因其更具有代表性,相比单参数d模型更能反映阀口的粗滞行为且容易实现。-3HeMATLAB仿真模块框图如图4所示。控制阀楷滞特性使用模型,被控-sG=e被控过程转换到时域并离散化后的形式如下对象为O):^,'。2-〇.2=厂x--U-xCfcr)(fc1)巧+21efc1)了(()p(()■4.2(。刺=xfc-(、(則上式中U为控制阀输出信号,r,为系统采样时间,此处7;取1秒。pHe模型使用MATLAB中S函数实现,MPC控制器使用MATLAB中M函数实现。MPC的M函数中需要将MIQP求解器在当前时刻采样点求得的最优解woace一rks设为全局变量保存到p中,下时刻再作为MIQP求解器的初值,W此一i类推不断更新初值,从而保存毎个采样点上时刻的值。代价函数cost化ficton中也需要将He模型的中间变量Ur和输出信号M设为全局变量保存到workspacep中一,便保存毎个采样点上时刻的值。仿真中MPC预测步长P=6控制步长M=2目标函数输出部分调节权重,,Q=100入部分调节权重R=O.OlnMPC控制(控制阀输入信号),输器输出信号=100=〇U的最大值1^?^,U,由于U中包含3个逻辑变量,逻辑变量只能mai取0或1,因此约宋条件包括线性不等式约宋4.15、等式约宋4.1、上下限约()(巧宋和二进制约束。MPC控制器输出信号U、过程输出信号y、He模型中间变量和输出信号U均在MPC的M函数外作为全局变量賦1次初值,且全部为零值。pMPC采样时间设置为0.5秒,仿真时间设为50秒。==11.当粗滞程度为怎3,时,假设外部干扰不存在记,环境较为理想情70的常数信MPC控制器补偿控制效果如图4-4况下,设定值为号,所示。54 浙江大学硕±学位论文=--1王〇〇1i!1111rr1[?SPAA90--PV80-/\ ̄-万6。>50-f-t/40-30--20--10-1tIn_[IIjIIiI051015202目3035404550Times〔s]4-4MFC控制器补偿控制效果图(无噪声)-4可W看出MPC控制器补偿控制效果较为埋想4,理想条件下从图,,即使控制阀存在一定程度的粘滞,过程输出仍然能够很好地跟踪设定值,并且很快稳定下来。2一.由于般过程系统都存在外巧千扰,为了使得仿真效果更巧近实际系统,,.01益大小在上述条件基础上,引入白噪声外部干扰噪声强度为0,增为化05,MPC控制器补偿控制效果如图4-54-6所示所示,局部细节图如图。=100IIIIiIIIII]…SP90-80A鸟心60-50--/t40--/30-20--/10-/」0---10tIIII[III05101520253035404550Timess()图4-5MPC控制器补偿控制效果(含白噪声)55 浙江大学硕±学位论文■---—II][[[[Iy—SPA己\八'/1I[1III1520253035404550Timess()图4-6MPC控制器补偿控制欢果(含白噪声局部图)一4-5レMPC补偿从图可乂看出,在过程输出加入定程度的白嗓声后,控制效-6可レ果仍然较为理想过程输出在设定值附近变化。从局部图4乂看出,过程输,出在30秒后进入稳定状态,相比无噪声20秒后进入稳定状态较晚。3=6=3.当粗滞程度增大时,,70尤,/d仿真时间设为秒,其他参数不变,MPC控制器补偿控制效呆如围4-7所示。-.=■XOOiiIII!IIr[11l—SP90_80-A心60-^:—4030-_I20--100-^!1I1〇IIjjI[[[III0510152025303540455055606570Timess()图4-7MPC控制器补偿控制效果(粘滞程度较大)一4-7可レ乂看出从图,当控制阀粗滞程度增大倍时,MPC补偿控制效果在过一-5-程输出稳定后相比图4差些,但是误差范围化只在1.51.5,[,]仍然可W满足实际系统的控制要求。56 浙江大学硕王学位论文4.对比常规ND控制器控制效果。仿真时间设置为7000秒,控制阀处于理想线性状态==,其他参数不变,经过多次调节P0.5,10.05,时控制效果较巧,当控制阀出现枯滞特性=6=-,/3,P瓜。,枯滞程度为尤d控制效果如图48所示1〇〇r1111;|90--80--70nmnnnmnnmnnmmmmnmmmnmmmmf60-t50-40-30--20-I—SP10-PV,,,0〇1000200030004000500060007000Times〔s)4-8P图ID控制器补偿控制效果(枯滞程度较大)4-8中可W看出从图,当控制阀出现较大程度的粘滞特性时,过程输出发生4-7振荡,并且不能稳定地跟踪设定值。而采用MPC控制器补偿控制时(图),同样的枯滞程度,但是控制效果明显好于常规PID控制。4.4本章小结本章详细研究了基于MIQP求解问题设计控制阀粗滞特性MPC补偿控制器的方法,通过引入逻辑变量分别将控制阀枯滞特性单参数d模型和双参数He模型转换为线性不等式约宋条件,从而利用MPC能够很好地解决有约宋问题的特点,将过程回路因控制阀粗滞特性所导致的输出无法堆确跟踪设定值的问题克服了。MPC对于模型失配等因素的处理能为较好,解决了无法完全准确建立被控IP问对象模型的难点,需要注。在求解MQ题上意目标函数与标准型在形式上的区别,可W等价转换为求解MINLP问题。通过《个仿真实验,验证了MPC补偿控制方法的可行性,其补偿效果也与控制阀粘滞程度相关。当枯滞程度较小时,MPC可W很好地控制过程输出稳定地跟踪设定值,而当粘滞程度较大时,虽然控制效果有所减弱,但相比常规門D57 浙江大学硕±学位论文控制,其控制性能有很大程度的提高,具体指标也与目标画数的调节参数、MIQP求解器优化程度等栖关。58 浙江大学硕±学位论文5结束语5.1研究工作总结现代工业控制系统规模日渐庞大,复杂度也越来越高,在成千上万条控制回一20% ̄30%路中种普遍现象,其中有约的振荡是由于回路中控制阀出,振荡是现了粗滞特性所导致的。在检测出控制阀存在粗滞特性后,采用合理的补偿控制方法可W克服枯滞作用,从而消除或减弱控制回路的振荡现象,这样保证了工业生产过程的正常运作,维护了企业的利益。粗滞特惟作为控制阀最常见的非线性一特性受到了越来越多的重视,而描滞特性的补偿控制已成为控制王程界的个重要研究课题。控制阀粗滞特性补偿控制的研究可W从W下几个方面入手:(1)如果不巧引入任何外部的枯滞补偿器,可W尝试先分析控制回路振荡产生的原因,利用描述函数判断回路是否存在极限环,再西绕控制器积分作用的调整来消除或减弱振荡现象;(2)在控制器输出信号中叠加不同类型的补偿信号,利用补偿信号的能量来克服控制阀的粗滞作用,从而驱动控制阀适到相应的位置。这种方法不需要调节控制器参数,单纯依赖补偿信号的作用来消除或减弱振荡现象,因此补偿信号的形式将决定最终的补偿效果;(3)可W考虑引入控制指标监测器,并且设置相应的阀值条件,根据控制指标的变化来决定是否引入摇滞补偿器或切换回常规控制一。这是种较为综合性的补偿控制方法,可W根据不同的控制系统吴活调整方窠,补偿过程中回路可能会被切换到开环模式:(4)如果引入控制阀粗滞补偿器,补偿器参数的设置可W结合控制回路历史振荡数据,通过历史振荡数据的分析来快速找到最佳补偿方窠。本文的主要研究工作是在已有的控制阀粗滞特性补偿方法的基础上,提出了相应的改进或扩展方窠レ,弥补了原有补偿方法的不足,使得粘滞补偿方法可义更好地运用到实际工业过程系统中。59 浙江大学硕±学位论文(1)Knocker脉冲幅值在线自整定寻优算法Knocker补偿法可能是目前效果最好的控制阀枯滞特性补偿方法,但是其单位脉冲信号参数的选取往往依据工程经验,缺少有效的选取方法,尤其是脉冲幅一值。本文提出了种基于黄金分割法的脉冲幅值在线自整定寻优算法,用于快速有效地找到最合适的Knocker参数,避免了人工试凑,提高了补偿效率和补偿精度。通过仿真实验,验证了该方法的有效性。(2)开环两步补偿法的新型实现方式两步补偿法的常规实现方式是在补偿器介入控制时,将控制回路巧换到开环模式,然后再利用补偿信号驱动控制阀进行王次除跃响应,阶跃响应间隔巧要求过程系统达到稳态,这样就使得两步补偿法的实现过程耗时很长。本文在分析了因控制阀枯滞特性而产生的振荡信号的历史数据,提出了基于He模型的开环两一种新型实现方式步补偿法的,该方法在回路切换到开环模式后手动驱动回路产一一生至少个周期的振荡信号,通过找到个周期内各变量间的关系,快速获取两歩补偿法的关键参数,,无需过程系统稳定下来。通过仿真实验将结果与常规两步补偿法的补偿效果进行了比较,验证了本文所提方法的有效性。(3)基于MIQP求解问题的MPC补偿控制MPC能有效处理有巧束问题,将控制阀摇滞特性的数据驱动模型转换为相应的输入约宋条件,即可利用MPC消除回路的振荡现象。基于上述研究,本文分别将控制阀枯滞特性单参数d模型和双参数He模型转换为相应的输入约来条I件,P求解器获取当前输入,并采用误差的平方作为MPC的目标函数利用MQe最优值。通过仿真实验,验证了本文所扩展的H模型转换的输入约宋条件的有效性。5.2研究工作展望目前控制阀枯滞特性补偿控制的研究己经取得了一定的成果,但是在研究过一一程中发现还有些枯滞补偿方法可W深入探索,也有许多问题有待进步巧究。未来的研究工作可L乂从W下几个方面着手:(1)控制阀楷滞模型转约束条件本文将控制阀粗滞特性单参数d模型转输入约宋条件扩展到双参数He模型60 浙江大学硕±学位论文转输入条件,还有更多复杂的控制阀粗滞特性数据驱动模型可W研究,例如Kano模型、Choudhury模型,但是这样的模型复杂度更高,需要引入更多的還揖变量才能完全描述模型特性。如果能将更多的控制阀粗滞模型转成相应的约宋条件,那么MPC补偿控制的实用化将更强。(2)系统化的控制阀姑滞补偿方窠针对不同的枯滞模型需要采用不同的补偿方窠,并且随着枯滞程度的变化也要调整相应的方窠因此将现有的控制阀艳滞补偿方窠整合成一,个智能的补偿系统,通过控制阀枯滞特性检测选择更适用的粗滞模型,再依据粗滞模型及枯滞程度智能切换到最佳的补偿器是很有实际意义的。(3)#回路报滞补偿控制本文所研究的控制阀粘滞特性补偿控制都是建立在单回路控制系统的基础上的,而实除工业过程系统中多回路十分常见。因此,如何在多回路稿合性较强的控制系统中设计猎滞补偿器需要进一步巧究。61 浙江大学硕±学位论文62 浙江大学硕±学位论文参考文献[1]M.A.A.S.Choudhury,S.L.Shah,N.F.Thornhill.DiagnosisofProcessNonlinearitiesandValveStictionDataDrivenApproaches.Springer,2008.口]M.A.A.S.Choudhury,N.F.Thornhill,S.L.Shah.Modellingvalvestiction.Cont200513-658trolEngineeringPracice5:641.,,()-mance3T.HSnd.Acontroooerfomonitor.ControEnneerinractce[]gglullpprlgigPi,131-9951:15431551.,()4NFTholTHa.Dttttr..rnhil.lundeecionanddianosisofoscillaioninconol[],gggEn-11loos.Controlineerinractice997510:3431354.pggP,,()[5]0.Taha,G.A.Dumont,M.S.Davies.Detectionand出agnosisofoscillationsin-controlloos.IEEEConferenceDecisionandControl1996310:24322437.p,,()6T.MiaoD.E,Sebor.Automaticdetectionofexcessiveloscillatorfeedback[],gyycontrolloops.IEEEInternationalConferenceonControlApphcations,1999,1:%9-%4.7W.L.Bialkowski.Dreamsversusrealit:aviewfrombothsidesofthea:[]ygpmanufacturingexcellencewithcomeonlythroughengineeringexcellence.PulpandPaCanada194119-27er993:.p^,。)D?目.Ender.ProcessCO打trolerformanceNotasoodasou化ink.Control[巧pgyEngineering,1993,40(10).[9]L.Desborough,R.Miller.Increasingcustomervalueofindustrialcontrol’-i化rinHon知ienceperformancemongeywellsexp?乂icheSymosium2002.p,10]T.Hagglund.Africtioncompensatorforneumaticcontrolvalves.Journalof[pProcessContro2002127-l:89904.,,巧)U刘巨良过程控制仪表.北京:化学工业出版社.1998.[],12何衍庆邱宣振,杨杰,王为国控制阀工程设计与应用.北京:化学工业[],,出版社.2005.[13]M.Bartys,民.Paton,MSyfert,S.D.L.Heras,J.Quevedo.Introduction1:0theDAMADICSactuatorFDIbenchmarkstudy.ControlEngineeringPractice,-596200614:577.,6()63 浙江大学硕±学位论文[M]J.M.F.Calado,J.M.G.S.D.Costa,M.Bartys,J.Korbicz.FDIapproach化theDAMADICSbenchmarkroblembasedonualitativereasonincouledwi也pqgptrP-fuzzyneuralnetworks.ConolEnineerinractice2006146:685698.gg,,()[1气D.Diigteg6r,E.Frisk,V.CocquempotjM.Krysander,M.Staroswiecki.Structuralanalysisoffaultisolabilityin化eDAMADICSbenchmark.ControlEn2006147-ineeri打ractice6:巧608.ggP,,()1J.M.Ko%ielny,M.Bartys,P.Rz)iejewski,J.SadaCosta.Actuatorfault[糾q出sti打gui浊abilitystudy化r化eDAMADICSbenchmarkproblem.ControlEnPacte2006-146:45652.ineerinric6gg,,()17F.Previdi,T.Parisini.Mode^freeactuatorfaultdetectionusinasectral[]gpestimationaproach:thecaseoftheDAMADICSbenchmarkroblem.Controlpp2006146-644EnineerinPractice6:%.gg,,()18MeaB.Huan,DetectionandDianosisofStictioninControlLoos:State.Jlli[],ggpoftheArtandAdvancedMethods.SringerVerlag2010.p,19A.Horch.Asimplemethodfordetectionofstictionincontrolvalves.Control[]En17101221-31ineerinPractice999:12.gg,,()20A.Sina.sr.Asmemethodfordetectnvave巧ctonn[]ghl,TLSalbuyipligliiittrouofPocessontro2005154-oscillaincono:371382lloos.JrnalrCl.gp,,()[21]R.Srinivasan,民?Rengaswam乂S.Narasimha打,R.Miller.ControlLoopPerformanceAssessment.2.HammersteinModelApproachforStictionDiagnosis.Industrial狂EngineeringChemistryResearch,2005,44。7):67-196728.M-22.A.A.S.ChoudhuryM.JainS.L.Shah.Stictiondefinitionmodellin[],,,g,tectonanduant-deiqification.JournalofProcessControl,2008,1834:()232-243.[23]R.Sri打ivasan,化.Rengaswamy.Techniquesforstictiondiagnosisandcomensationinrocesscontrolloos.AmericanControlConference2006.ppp,24THSActoncomensatororneumatccontrovavesournaof.lund.friipfill.Jl[]ggptr2002-Pro说ssConol12:897904.,,巧)Att-tttr25Tl.HSlund.uomaiconlineesimaio。ofbacklashinconolloos.Journa[]ggp64 浙江大学硕±学位论文tr2007-Proc的SConol1717:489499.,,()口6]L.Ivan,S.Lakshminarayanan.Anewun近edapproachtovalvestictiont-quanificationandcompensation.Ind.Eng,Chem.Res2009487:34743483.,,()[27]M.Cuadros,C.Munaro,S.Munareto.Novelmode^freeapproachforstictiont-comensaionincontrolvalves.Ind.En.Chem,Res20125125:84658476.pg,,()口8]民.Siinivasan,艮?Rengaswamy.Stictioncompensationinproce巧controlloops:Aframeworkforintegratingmeasuredandcompensation.Ind.Eng.Chem.Res.200544249164-17:94.,()29R.SrinivasanR.Renaswam.Aroachesforefficientstictioncomensation[],gyppp2008-inrocesscontrolvalves.ComutersandChemicalEnineering,3212:ppg,()218-229.[30]M.Cuadros,C.Munaro,S.Munareto.Improvedstictioncompensationinneumaticcontrolvalves.ComutersandChemicalEnineerin201238:ppgg,,-106114.1WC-3J.ang.losedloocomensationmethodforoscillationscausedbycontrol[]pp-valvestiction.Ind.Eng.Chem.Res2013,5236:1300613019.,()口2]Gerry,J.;Ruel,M.How化measureandcombatvalvestictiononline.2011ISATechnicalConf.HoustonTexasSetember2011.,,,p,[33]M.Mohammad,B.Huang.Compensationofcontrolvalvestictionthroughtr-conollertuninJtr201222111g.ournalofProcessConol:80089.,,巧)34付川.103110:,了维明,控制回發中的阀口迟滞补偿方法自动化仪表,20[],()8-11.口引K.J.AstrSm,T.HSgglund.AdvancedFIDcontrol.ResearchTrianglePark,NorthCarolina:InstrumentSocietyofAmerica,2006.36M.A.A.S.Choudhury,S.L.Shah,N.F.Thornhill,D.S.Shook.Automatic[]detectionanduantificationofstictionincontrolvalves.ControlEngineeringqPratce2006-ci142:13%1412.,,y)[37]M.Kano,H.Maruta,H.Kugemoto,K.Shimizu.PracticalmodelanddetectiontttINIFADtralgorihmforvalvesicion,CSymposiumonynamicsandConolofProce巧SstemsCambrideMAJu-yl57,2004.,g,,y65 浙江大学硕±学位论文38民.Srinivasan,民.民enaswamU.NallasivamV.Raavelu.Issuesi打modelin[]gy,,jg巧ictioninrocesscontrolvalves.AmericanControlConference2008:p,3374-3379.[39]何熊熊,邹涛,杨悦,李信,气动控制阀静摩擦的Kano模型与检测方法的一-进步巧究.化工学报2008597:16911697.,,()40A-.StenmanRGustafssonK.Forsman.Asetatonsetmenibadmehodfor[],,gdetectionofstictionincontrolvalves.InternationalJournalofAdaptiveControl--andSinalProcessin20031779:化5634.gg,,()4--1S.L.Chen,K.K.Tan,S.Huang.TwoLayerBinaryTreeDataDrivenModelfor[]Ch-ValveStiction.Ind.Eng.em.Res2008478:28422848.,,()[42]Q.P.He,J.Wang,M.Potmann,SJ.Qin.Acurvefitingmethodfordetectingvalvestictio打inoscillatingcontrolloops.Ind.Eng.Ghem.艮es,2007,46(13):4549-4560.[4^D.Sebo巧,T.Edgar,D.Melichamp.ProcessDynamicsandControl.Wiley,2011,33:45.()44SJ.inTIA.Badwell.Asurveofindustrialmodelredictivecontrol[]Q,gyp-1:echnology.ControlEngineeringPractice,2003117:132.,()[45A.Ghazzawi,E.AliA.Noud,E.Zafiriou.Onlinetuninstrateibrmodel],ggyred-ictivecontrollers.JournalofProcessControl20013:2652.p,1184,()[46]J.M.Mac均owski.PredictiveControlwithConstraints.PrenticeHall,Inc.2002.47H.ZabiriY.Samudia.Ahbridformulationanddesinofmodelredictive[],yygpcontrolforsystemsunderactuatorsaturationandbacklash.JournalofProcessCon-709tro化2006167:693.,()[48]A.Bemporad,M.Morari.ControlofSystemsIntegratingLogic,Dynamics,andtrt-ConsaAt19935407427ins.uomatica88:.,,巧)[49]史纽梅,冯丽辉,钟伟红,彭云峰.基于模型的预测控制算法研究.最明理王大学学报-:理工2004276:7275版,.,()-MPC0HZabYIPBinthPCont.iri.Samudia.MasederesenceofrolValve口],yQ如ction.ChemicalProductandProcessModeling,2009,4口?)66 浙江大学硕±学位论文致谢值此论文即将完成之际,我要向所有帮助过我的老师、同学、朋友和家人们表示由衷的感谢!本文在谢磊副教授的悉瓜指导和亲切关怀下完成。谢老师知识渊博、敬职敬业、治学严谨,,是我学习和工作中的榜样。谢老师不仅仅在学业上对我耐伯指点而且在生活上对我也关怀备至,经常通过分享自己的人生经历来鼓励我,使我获益良多此我表示对谢磊副教搜最真擎的敬意和由衷的感谢!,在感谢苏宏业教授レ、陈剑教授、徐巍华副教授乂及系里其他老师在学习和生活中给予我的帮助。特别感谢实验窒成员谈斐棋、郭子旭、宋超给予的帮助和鼓励,通、王培宇过与他们的交流讨论,使我得到了很多启发。感谢欧阳权、王凡、李爽、张鹏飞、陈新如等多年来的关怀和帮助。一特别感谢我的家人,在我多年求学生涯中直默默在背后支持我,永远祝福他们!王教巧公M年先j子走4(867 浙江大学硕±学位论文68 浙江大学硕±学位论文作者简介王挺任,男,1991年6月出生于黑龙江省鹤岗市萝北县,中国共产党党员,2009年考入华北电为大学自动化学院,2013年毕业,获王学学击学位。同年由于成绩优异,表现突出,免试推荐到浙注大学控制科学与工程系,在智能系统与控制研究所攻读硕壬学位,师从谢磊副教授。研究方向包括:工业过程控制和控制阀非线性特性研究。69 浙江大学硕±学位论文70 浙江大学硕±学位论文作者在攻读硕壬学位期间的科巧成果发表文章[1]TingrenWang,LeiXie,FeiqiTan,HongyeSu.Anewimplementationof-oen-tttrplooptwomovecomensationmehodforoscillaionscausedbconolpy-148-43valvestiction.IfacPapersOnline,205,8:4338.()发明专利一1谢磊.种因阀口粘滞特性而引起控制回路振荡的补偿方法.申请[],王挺任号:20151008098LX公布号:CN104635497A71

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。
大家都在看
近期热门
关闭