保捱科技网
您的当前位置:首页冷带轧机轧辊偏心补偿控制的应用研究

冷带轧机轧辊偏心补偿控制的应用研究

来源:保捱科技网
摘要    对轧辊偏心问题的研究,在我国尚处于起步阶段,成功实例很少.随着各行各业对板带材质量要求的不断提高,轧辊偏心对轧出厚度的影响己变得不容忽视,成了现代板带材轧制领域研究的热点问题。本文以唐山恒通精密薄板有限公司的1270mm可逆冷带轧机液压微调    厚控系统为研究背景,深入剖析了轧辊偏心对板带材厚度的影响机理,提出了1270mm可逆冷带轧机利用快速傅立叶变换法进行轧辊偏心补偿的实施方案。针对补偿算法在具体实现中的一些技术问题做了深入研究,并对该控制方案进行了仿真试验,证明了控制方案的可行性。依据控制方案选取了相应的外围辅助元件,将控制方案应用到于实际生产中,使板带精度得到了较大的提高。/针对液压板厚控制系统内位置闭环采用常规PI  D控制方法参数调整范围小,鲁棒性差的缺点,本文设计了增强型PID自动位置控制算法,提高了伺服阀阀芯位置反馈环节的调节作用,使系统的动态响应和稳定性有了明显提高。    为了进一步改善液压板厚控制系统的响应性能,本文提出了缸位移微分负反馈对位置内环的功能进行补偿。仿真结果表明,加入缸位移微分负反馈使系统的响应特性得到了明显提高。此外,本文还用机理分析的方法建立了液压板厚控制系统的动态数学    模型,并对数学模型进行了必要的简化。仿真结果表明,用机理方法建立的数学模型可以精确评估系统的各项性能。)夕办一关键词冷轧机;板厚控制轧辊偏心;虚厄面蔽行派履快速傅立叶变换卜      /                      燕山大学工学硕士学位论文Abstract    There are fewer successful examples in roller eccentricity researchdomain, which is in the preliminary stage in china. With the increaserequest of the every walk of life to strip quality, influence which theeccentricity affect in strip quality has no longer been neglect and becomethe hotspot problem in study of realm of modern strip rolling.Bas    ing on hydraulic micro-adjustment AGC of reversing cold strip millin Tangshan HengTong high-quality strip limited company, this paper,analyze the effect principle which roller eccentricity affect strip gauge,bring up eccentricity compensation control implementary print which applyFFT method. The paper carefully researches some technical question ofcompensation arithmetic, simulate and testify its feasibility. According asthe practical of control blue print, the author select corresponding peripheryassistant element and apply the control blue print in practical product, raisethe quality of cold strip.Ai    m at the less adjustment range and robustness of general PID controlin hydraulic AGC system, the paper design enhanced PID automaticposition control arithmetic, raise the adjustment effect of the core ofservovalve position, evidently improve the dynamic response and stability.I    n order to improvement the response characteristic, the paper bring upnegative differential feedback of hydraulic cylinder displacement tocompensate the function of position inner loop, emulational result indicatethat the addition of negative differential feedback of hydraulic cylinderdisplacement raise evidently the response characteristic of system.I    n addition to, the paper constructs the dynamic model of hydraulicAGC system through the mechanism analysis method and simply it.Emulational result makes clear that the mathematic model is practical toanalyze accurately the every characteristic of systemKeywords: cold mill; AGC; roller eccentricity; enhanced PID automaticposi          tion control arithmetic; negative differential feedback ofhydr            aulic cylinder displacement; FFT11                                        第I章绪论第1章绪论冷轧薄板(    带)广泛用于汽车制造、食品包装、家用电气、机械、轻工、仪表、通讯和军事工业等各领域,是非常重要的通用钢材之一,从产量上看约占冷轧板产量的20%左右。为了生产出更高精度的板带材产品,占领国际市场,各国相继投入了大量的人力、物力,开发研制了多种现代化大型轧机。具有代表性的有日本日立公司研制的工作辊可轴向移动的HC轧机,德国施罗曼一西马克公司研制的连续可变凸度CVC轧机,日本三菱公司研制的成对轧辊交叉式的PC轧机,德国曼内斯曼一德马克公司的UPC轧机等。这些轧机配以现代化控制手段,可实现多功能精密轧制,生产出微米级厚度精度的板带材产品。在我国,上千台套装备较为落后的国产轧机如何在投资较少的情况下    应用一些较先进的控制技术,以提高带材质量川,还有许多理论和技术问题需要研究,轧辊偏心补偿控制便是其中要重点研究的问题之一。1.1选题的背景、目的及意义    本文是以横向课题“唐山恒通1270mm冷带轧机液压微调厚控系统的研制”为基础撰写的。唐山恒通公司始建于1996年,原有一台1100mm轧机和一台1000mm平整机,现在又从国外引进了2台1270mm冷带轧机。这2台1270mm冷轧机原系国外的二手设备,仅引进了机械设备、电机等电气设备,无传动及厚控系统部分。轧机厚控微调部分原为电动压下,并预留了液压压下的位置。引进后传动部分以及液压AGC部分由国内研制。车L机在控制上采用了一些独特的方法,技术含量较高,与我校先前研制的几套液压AGC厚控系统存在一些差别。作者有幸参加了整个控制系统的设计、安装和调试工作,在学习硬件设计的基础上,协助导师进行了软件的编程和调试;并根据现场出现的新问题及操作人员提出的具体建议对系统不断地进行改进。由于轧辊偏心不仅对带钢厚度均匀性产生直接的影响,而且会使常规调节器的调节质量恶化,对带钢厚度均匀性产生间接的不良影响,如果不消除轧辊偏心对板厚的影响,很难轧出更高精度的板带材。针对这种情况,本文对轧辊偏心问题进行了深入的探讨和研究,制定出了车L辊偏心总体控制方案,改进了软件的算法,解决了其在具体实现中的一燕山大学工学硕士学位论文些关键问题,将快速傅立叶变换成功地应用到液压板厚控制系统中,极大地提高了板厚控制精度。为了改善系统的鲁棒性,本文设计了增强型PID自动位置控制算法,很容易使系统的动态特性、稳定性达到要求,有效地降低了现场调试和操作的难度。应用缸位移微分负反馈使系统动静、态特性得到了明显改善,满足了生产要求。经过改造后的轧机,取得了巨大的经济效益,无论其轧制速度和板厚控制精度都有了很大的提高。每月最高可轧制带钢9000吨,年轧制量可达10万吨以上,板厚控制精度达到了微米级,完全达到了设计能力。为今后改造和设计同类轧机提供了可靠的实践检验和设计依据.1.2板带轧机厚度控制发展综述    板带材轧制生产,20世纪70年代以来,取得了飞跃的发展。生产效率己经达到相当高的水平,产品质量也有很大的提高,如在尺寸精度上,冷轧板材己经达到纵向和横向都具有较高的精度P1;在产品性能上,已经初步作到按预期目标控制生产。这样快速的发展是在轧制理论的指导下,在机、电、仪表和控制设备的基础上,借助电子计算机而取得的。板带材轧制理论近年来己经发展成为一门比较完整的技术科学,因此对轧制理论的发展进行必要的讨论和回顾,进一步估计它的发展趋势和指出必须解决的问题,将具有一定的现实意义,以下将从三个方面进行论述。1.2.1轧制的早期历史在1    8世纪早期,轧机已普遍使用于英格兰和欧洲[2-41。大约此时瑞典的伟大力学天才波尔希姆写过关于轧机的文章,他自己设计了一架十分相似于现代劳特轧机的轧机,只是他采用了四个轧辊,并用支持辊驱动。    一台真正预示以后轧机形状的“轧制铅板的机器”于1728年由英格兰引入法兰西。同年在英格兰佩恩发表了一项专利,用轧机将锻制棒材轧成所需要的形状。可是,佩恩的想法似乎并末变为现实。1759年,一项关于“抛光及轧制金属”的专利被授予英格兰的布洛克莱,确切地说其成果就是使用者可按照需要对轧辊进行刻槽0 17“年另一个英国人拍内尔接受了一项有关刻槽的轧辊的专利,利用成对的齿轮箱使轧辊协调一致地转动。这以前,轧辊都是单独驱动的,不相等的转动速率引起轧辊过度磨损,也使得在每个轧辊的两侧有必要安置导板。z                                       第1章绪论在此同一时期内,这些热轧机的一般外貌开始改变成现代的形状[[21例如,铸造的牌坊和轧机每边的单独螺丝,都在1783年普莱菲尔德所提出的英国专利说明中得到了描述。此时的热轧机主要有两个辊(一个置于另一个之上)的单机架组成,此种常用的薄板热轧机被称为“两辊回递式轧机”。一个轧制车间可使用一台以上的这种轧机,每台轧机完成的工作阶段,但每台轧机的操作原则仍然一样。18世纪也出现了金属在连续机架中轧制的连轧机。见之于记录的第一    台真正连轧机,于1766年在英国为福特取得了专利,它用以热轧线材。科克肖特和克劳谢大约在1790年于英国设菲尔德附近创立了一台四辊式连4L机,它大约为1524mm长,不到609mm高,其生产能力可能是每天一吨或两吨[5-61。他们于1798年发表的专利,是关于轧制厚板和薄板的连轧机的。在同一年黑兹尔丹对一台线材轧机增加了机械导板。在1    9世纪初,英国的工业正方兴末艾,对钢铁提出空前的需要,因此,轧机的发展很快,在钢铁生产中也很重要。伯肯肖在1820年开始使用第一台铁轨轧机,生产长4572 ̄ 5486mm的鱼腹式熟铁轨。1831年在美国首次轧出T形轨(基本设计与今日所用的钢轨相同)。1849年佐雷在巴黎轧出第一根工字钢梁。在19世纪中叶,第一台可逆式轧板机于英国帕克盖特工厂投入生产,    并于1854年用以轧制制造“远东”号蒸汽船的板材[7-8]。在1848年,德国伦德斯多夫的德伦发明了万能轧机,他在大约七年以后建立了第一台这种形式的轧机。所有轧辊在同一垂直平面的四辊轧机,于1872年在英国沃林顿为布莱克利所采用,以精轧熟铁块,作为轧制铁轨的原料。热连续轧钢机到1890年左右有一个大的发展。1892年在波希米亚的    特普利兹建立了一台半连续热轧带钢轧机,它具有机械齿轮传动的连续式二辊精轧机列[9-101.据报道,己轧制出宽1270mm、厚2-3mm、长18288mm的薄板扳。该轧机采用包括两台三辊轮机的粗轧机列和包括五个机架的精轧机列.轧辊尺寸为625 X 1498mm,机架间距为2743mm。每一机列由746KW的原动机传动。这一时期人们开始把辊轧的实践经验总结为艺术或者技巧加以掌握,    但是没有能进行较深入的科学概括[[21。虽然也提出了若千假设,但是因为缺乏必要的物理模拟试验,一些假设的根据是不足的,有些甚至是不正确燕山大学工学硕士学位论文的。很多方面的阐述都没能正确反映实际过程。例如,有关轧制压力的计算,在1859年格日绍夫提出单位压力在变形区入口处有一个极大值,而且沿着人口到出口的方向逐渐降低,因而它与轧件在变形区的厚度成比例。轧制理论也才刚刚开始建立,远没有达到应用的程度[10-u[;单回路调节的自动控制理论,尚未应用于控制轧机这类较复杂的机器。因此,轧制理论从成熟程度来看,还处于萌芽阶段,从理论水平来看,也只是处于艺术阶段。过去有一个习惯说法,轧钢是一门艺术,大概就源于此。1.2.2冷轧的发展历史    虽然自1660年以后,己经在二辊轧机上进行了板带钢的热轧,但是,直至十八世纪末,板带钢的冷轧却似乎并末成功实现。作为一个生产过程,板带钢冷轧的发展,只是在具有较小的中辊(工作辊)的劳特三辊冷轧机发展以后,才真正得到推动。但小工作辊的优点已更早地为人们所认识,因为波尔希姆描述了一台用于板带热轧的四辊轧机[I-21,采用细长的熟铁工作辊和大的铸铁支持辊,这是由于“小轧辊比大轧辊具有大得多的延伸能力”。冷轧起初只是应用于二辊和双二辊轧机上,后来美国也采用了威尔莫    特和曼式的四辊和六辊的多辊轧机。用于冷轧带钢的第一台四辊轧机在近至1923年才首先被勒德隆钢公司采纳作试生产之用.轧辊辊颈轴承的改进,也促使冷轧机的应用日益增加。滚动轴承早在    1890年首先用于二辊冷轧机,在1909年用于多辊冷轧机的支持辊和在1926年用于四辊轧机。    二辊式可逆冷轧机于20世纪20年代首先用于德国(早在1917年在专利文献中已有所披露),而四辊可逆式则用于1932年.美国第一台这样的冷轧机于1933年安装于加里。带钢连续冷轧的首次记录要追溯到大约1    904年,当时,韦斯特利奇伯格刚公司安装并开动了一台四机架二辊连轧机,每个机架单独用速度可调的直流电机驱动。具有机架间张力和张力卷取机的真正的连轧机操作,大约于1915年在匹兹堡的英里斯一贝利钢公司和苏必利尔钢公司安装的轧机上才得到发展。第一台四机架四辊冷轧机由美国轧机公司于1926年在巴特勒工厂投入生产[12-13]在可逆式轧机和连轧机的操作中,随着所轧带刚长度增加和轧制速度    第1章绪论愈来愈高,如何处理轧件的问题必须引起注意.冷带卷取机似乎领先于热带卷取机大约十年,第一台冷卷取机大约在1893年由施米茨公司建于德国。该卷取机设计较好,具有楔形可胀缩的部分,并且由轧机来驱动,通过打滑的皮带来控制张力。大约1920年,苏必利尔钢公司和威斯汀豪斯电气公司联合发展的单独电力驱动的卷取机,是通过控制电流来维持恒张力的。    卷取机和电力驱动的共同发展,首先带来卷取机与轧机之间的,然后引起连轧机邻接机架之间的带钢张力控制的发展。第一台连轧机操作时,机架之间的带钢是松弛的,以后技术进步到利用活套下垂或活套张紧,最后大约于1920年,进展到在机架之间采用大的张力。在20世纪30年代,热轧钢带的冷轧,由一种比较专门化的小规模的    生产过程,发展到在冷轧板带和扁钢生产中占头等重要的地位。冷轧带钢的最大有效宽度,自1925年以后得到迅速增长,而且在一定宽度条件下,其最小厚度也得到减小p41    30年代建立的典型薄板轧机是三机架的2133.6mm连轧机,此轧机采用直径520mm的工作辊,直径1422mm的支持辊,驱动电机总容量5110KW,轧制速度高达2.75m/s。为了生产镀锡板(所轧带钢的最终厚度级别为0.254mm),在本世纪三十年代,五机架连轧机很快投入使用.这样一种轧机于1941年安装于欧文工厂,它采用533/1346X 1219mm的轧辊和8280KW的驱动电机。它能在高达19m/s的速度下轧出带钢。1960年开始引用六机架镀锡板轧机,装备了更大的功率及稍大的工作轧辊(直径一般为533mm,此轧机为584mm)。在这种轧机上,采用双电动机驱动,并且部分地在计算机控制下操作。大量物理模拟试验为轧制理论发展成为一门技术科学打下了必要的基    础。但是在这时,轧制理论还没有能形成一门相对的、系统的科学,因为这些物理模拟试验,无论在试验方法上、结果分析上还是理论概括上,都处于孤立、静止和割裂的状况.这个时期还很难用数学来帮助解决试验中的问题,试验结果也难以直接指导生产。轧制生产及工艺还需要依靠机械学、电力学、金属学、力学以及工作者的经验和技艺.轧制理论尚处于发展阶段;从理论水平来看,它处于物理模拟阶段17一1,]5。70年代起,液压厚控技术(    AGC)开始应用,使板厚控制技术产生了重燕山大学工学硕士学位论文重大变革[51。液压AGC的响应速度比电动AGC快2个数量级以上,因此,使AGC的内环执行机构几乎可以近似为一个比例环节(相对于AGC的响应速度),以实现可变等效刚度的控制效果。随着液压技术与计算机技术的结合,板厚控制技术开始向着大型化、高速化、连续化的方向发展(19-211    随着控制理论新成果在冷轧宽带刚轧机上的应用,冷轧带刚产品的精度也进一步提高0 1990年8月投产的日本新日铁八蟠制铁所的冷连轧机生产线,其板厚控制系统由31个回路的自动板厚控制AGC和8回路的自动张力控制ATC构成[9-10]高性能交流传动装置的应用以及液压压下装置的能力最大限度的发挥,使板厚精度控制在士0.5%以内,成为这一阶段冷连轧机计算机网络控制的典型代表。    冷轧薄板发展如此迅速的主要原因是[1"31,(    1)钢材在热轧过程中的温降和温度分布不均给生产带来了难题,特别是在轧制厚度小而长度大的薄板带产品时,冷却上的差异引起的轧件首尾温差往往使产品尺寸超出公差范围,性能出现显著差异。现代热连轧机,目前设计可能轧出最小厚度为1.2mm,但实际生产中很少生产1.8mm或1.5m,以下的热轧板卷.而冷轧则不存在热轧温降与温度不均匀的弊病,可以得到厚度更薄、精度更高的冷轧带钢和冷轧薄板。现代冷连轧宽带轧机和双机架二次冷轧可生产厚度为0.10-0.17mm的冷轧薄板,作为镀锡原板,即使不经二次冷轧也可生产0.2-3.5mm厚的冷轧薄板。现代可逆式冷轧机可生产0.15-3.5mm厚的冷轧板。这些产品用热轧方法是不可能生产的.从厚度精度上看,现代热轧厚度精度通常为土50u,m,而现代冷连轧板厚精度高达土5pm,比热轧厚度精度高10倍。(    2)目前热轧工艺技术水平尚不能使刚板表面在热轧过程中不被氧化,也不能完全避免由氧化铁皮造成的表面质量不良。因此热轧不适于生产表面光洁程度要求较高的板带钢产品。热轧板表面粗糙度热轧状态为20,um.酸洗后为25u,m。而冷轧板表面清洁光亮,并可根据不同用途制造不同表面粗糙度的钢板。冷轧板按表面粗糙度分为3种:一种是无光泽的钢板,其表面租糙度为3"-1印m,一般适用于作冲压部件,并且当需涂喷刷漆时这种钢板附着性较强;第二种是光亮板,其表面粗糙度大于0.知m,这种钢板主要作为装饰镀铬用厚板等;第三种是压印花纹钢板,采用表面具有70-12伽m凸凹的平整辊平整钢板,这种钢板用于仪表壳及家具装饰等。第t章绪论这样的表面质量热轧是无法满足的。(    3)冷轧钢板的另一突出优点是性能好、品种多、用途广。通过一定的冷轧变形程度与冷轧后的热处理恰当的配合,可以在比较宽的范围内满足用户的要求。如汽车用薄板几乎全部须经冲压成型,这样深冲性能就成为薄板生产和使用的核心问题。热轧薄板的塑性应变比R仅可达到0.8-0.95。而冷轧第一代沸腾钢汽车板万为1.0.1.2,第二代08A1钢为1.4-1.8,而第三代冷轧汽车板万为1.8^-2.8,这是热轧无法达到的。由于上述原因1    21,冷轧薄板的生产得到迅速发展。早期的冷轧板轧制速度不到1 m/s,而今已达41.6m/so钢板的宽度1905年是406mm, 1925年是914mm,而今最宽可达2337mm o钢卷重量也从几吨发展到60吨。一座现代化的冷轧厂年产量可以达到250万吨。从世界范围来看:日本1960年冷轧板产量100万吨,而30年后的今    天生产能力己达2000多万吨,增加了20倍。原西德1950年年产冷轧板16.3万吨,而1974年达821.3万吨,增加了50多倍。1.2.3冷轧带钢在我国的发展情况我国冷轧宽带钢的生产开始于1    960年,首先建立了1700mm单机可逆式冷轧机,以后陆续投产了1200mm单机可逆式冷轧机,MKW1400mm偏八辊轧机、1150mm二十辊冷轧机和1250mmHC单机可逆式冷轧机等,70年代投产了我国第一套1700mm连续式五机架冷轧机,1988年建成了2030mm五机架全连续冷轧机.现在我国投入生产的宽带钢轧机有18套,窄带钢轧机有418套。在这30多年中,我国冷轧薄板生产能力增加了20多倍,生产装备技术水平已由只能生产低碳薄板而发展到能生产高碳钢、合金钢、高合金钢、不锈耐热冷轧薄板、镀锌板、涂层钢板、塑料复合薄板和硅钢片等。但随着四化建设的发展,无论在数量和品种质量上都远远满足不了四化建设发展的需要,为此我们必须增建新轧机,改造现有冷轧机,大力发展冷轧生产[22-2311.3偏心补偿控制    在带材轧制过程中,有很多因素会引起轧件的厚度偏差【14-711。对于变化比较缓慢的扰动[24.261,如来料厚度、支撑辊油膜、轧辊热膨胀等一般采取相应的补偿方法来克服其影响,并取得了令人满意的效果。但轧辊偏心7                                         燕山大学工学硕士学位论文是一种高频干扰,不仅用一般的方法难以凑效,而且使压力AGC产生相反作用的调节[[181随着各行各业对板带材质量要求的不断提高,轧辊偏心轧出厚度的影响己变得不容忽视,成了现代板带材轧制领域研究的热点问题。    车L辊偏心可归纳为两种基本类型,即由辊身和辊颈的不同轴度产生的偏心及轧辊本身具有的椭圆度造成的偏心。第一种偏心主要是由机械作用形成的,它包括以下几个方面:轧辊轴心和轴套之间的不一致轴向(如滑动)的残余应力导致的轧辊弯曲;轧辊打磨过程中造成的磨口偏心;轧辊的不均匀磨损:轴承和辊颈之间的键销作用形成的辊颈弯曲。第二种偏心主要是由辊身的不对称应力,轧辊圆周柔量的不一致以及轧辊的加工形成的。当然,实际情况下轧辊的偏心可能是两种情况的综合[z7-z81通常,偏心信号混杂在各种扰动和随机信号之中,偏心补偿的关键就    是设法获取偏心信号或其主要成份,进行控制。由于液压压下系统响应快,所以大部分现代轧机均采用液压压下来进行调节。对轧辊偏心的研究国外始于70年代,我国在这方面的研究才刚刚起步,但大都处于理论仿真阶段,成功实例很少。偏心补偿领域较成功的技术是从轧制力传感器输出信号中提取偏心分量进行补偿。如各种数字滤波器偏心补偿装置以及傅立叶变换偏心补偿装置等[30-331在冷带轧钢轧制系统研究过程中,偏心控制问题的研究尚应深入发    展,还一些问题急待解决。首先是如何把偏心控制方法应用于实际生产之中。其次是如何解决偏心信号中的高次谐波成分的补偿。第三是随着对带材质量要求的不断提高,在补偿支撑辊偏心的同时必须对工作辊的偏心进行补偿。最后,如果轧机运行速度较高,则轧辊的偏心频率也相应的增高,此时液压压下控制装置的频率特性就不一定能够满足要求,尤其是相位延迟问题,这种情况下必须解决。目前在轧辊偏心补偿研究领域又有人提出了一些新的方法,如神经网    络自适应滤波在偏心控制中的应用,利用小波消噪原理来设计消除轧辊偏心的小波分析器等,但大都处于理论研究阶段,如何将其用于实际工程,有待大家进一步研究。第1章绪论1.4本文主要研究内容(    I)在参阅国内外大量文献的基础上,结合偏心补偿控制问题,进行了板带轧机厚度控制发展的文献综述。(    2)本文以唐山恒通精密薄板有限公司的1270mm可逆冷带轧机液压微调厚控系统为研究背景,深入剖析了轧辊偏心对板带材厚度的影响机理,提出了1270mm可逆冷带轧机利用快速傅立叶变换法进行轧辊偏心补偿的实施方案。针对轧机运行过程偏心信号的抽样、截断,在控制中过程中偏心补偿信号的给定方法及轧机外围辅助元件的选取这些在实际应用中的技术问题进行了深入的研究。经过理论分析,实验室仿真和大量的现场实验最后将以快速傅立叶变换为基础的偏心补偿控制顺利地运用到了实际生产之中。(    3)由于液压板厚控制系统中常规型PID控制方法完全依靠PID控制器参数的调节作用,参数调整范围较小,系统鲁棒性较差。本文采用增强型PID自动位置控制算法来对其进行了改进,加强了伺服阀阀芯位置反馈环节的调节作用,相对减弱了PID控制器的调节作用,很容易使系统的动态特性、稳定性达到要求,有效地降低了现场调试和操作的难度。根据生产的要求,应用缸位移微分负反馈和来对位置内环的功能进行补偿,改善了系统的动态性能,减小了系统的相位滞后,降低了轧制成品的厚度波动。(    4)在对液压厚控系统元件性能认真分析的基础上,用机理分析的方法建立了系统的动态数学模型,并对数学模型进行了必要的简化。仿真结果表明,用机理方法建立的数学模型可以精确评估系统的各项性能。燕山大学工学硕士学位论文第2章板厚控制基础2.1机座的弹跳曲线和弹跳方程    轧制时,在轧制力的作用下,轧机工作机座(轧辊及轴承、压下装置和机架等)会产生一定量的弹性形变。工作机组的总变形量可达2-6mm,以致对轧制成品的精度造成影响,在轧制薄而宽的带材来说川,影响更为明显。    设轧件进入辊缝前,原始辊缝等于so:轧件进入辊缝,开始轧制后,在轧制力尸的作用下,工作机座会产生弹性形变,使辊缝增大,使与轧件接触的轧辊形状呈凹形。弹性变形的结果使实际压下量减小,轧件的出口厚度大于空载辊缝值。轧后轧件厚度为[19-201式中h一轧后轧件厚度(简称轧件厚度)(m );So一轧辊空载辊缝(m       );.      f一机座弹性变形(m ),可见,要得到厚度为h的轧件,应把空载辊缝S。调整到比轧件厚度h小    h=SO                          + f                      (2-1)一个机座弹性变形.f的量值上。2.1.1弹跳曲线机座的弹性变形与轧制力有关。图2-    1表示了机座弹性变形量f与轧制力尸之间的关系曲线,称为机座弹性变形曲线或称弹跳曲线.机座弹性变形曲线直线段的斜率,称为机座的纵向刚度模数[171。轧机的纵向刚度模数可表示为[171△尸              KS一tga=万式中K:一轧机纵向刚度模数((N/二):      AP一轧制力的变化量(N),O      f一弹性变形的变化量(M).(2-2)在轧制力较小的时候,机座弹性变形与轧制力成曲线关系,这是由于    机座各零件之间的接触面凹凸不平和轧辊的非线性接触变形造成的。当轧第2章板厚控制基础机座弹性变形曲线      图2-1curve of mill base distortionFig.2-1 Spring制力达到一定数值后,机座弹性变形与轧制力就成直线关系。由于通常都不在低负荷区轧制,故一般都把弹性曲线看成线性。机座的弹性变形除与轧制力有关外,还与机座的纵向刚度模数有关。    机座的纵向刚度模数越大,由此造成的轧件纵向厚差也就越小,即轧机控制轧件纵向厚差的能力越强。2.1.2弹跳方程从图2-1可以看出,轧后轧件厚度h可以近似表示为[211h=So +4S=S。十二KS    (2-3)式中h一板带材的实际出口厚度(,);so一空载辊缝值(M)      ;K、一轧机纵向刚度模数(    N1m)c      尸一轧制力(N)o式(    (2-3)称为机座的弹跳方程,它表示了轧件厚度与空载辊缝、轧制力和轧机纵向刚度模数之间的关系。由于机座各零件的非线性接触变形不稳定,每次换辊之后都有变化,故弹跳曲线的非直线部分是经常变化的,式(2-3)很难在实际中应用。因此在实际生产中,必须首先进行压靠,即将轧辊预压靠到一定的压力PO,将此时轧辊辊缝值记为零,称为人工零位。如燕山大学工学硕士学位论文图2-2人工零位时弹跳曲线Fig.2-2 Spring curve of manpower zero图2一所示,曲线A'就是预压靠曲线。将压靠力为PO时的轧辊辊缝记为人工零位,实际值为一负值,设为弓,当实际的空载辊缝为So时,辊缝指示器读数则为瓦二S+So。曲线A与压靠曲线A'完全对称,为此辊缝下的机座弹跳曲线。考虑到轧辊液体轴承油膜厚度变化的影响,精确的弹跳方程可表示为P一Po一(s一h=So+Ks式中S。一人工零位后的空载辊缝指示值(P      O一人工零位时的预压靠力(N);(      5—液体摩擦轴承油膜厚度(m );      8a一对应于人工零位时的油膜厚度(m).轧辊液体摩擦轴承油膜厚度变化量s-so与轧制力和工作辊转速有    关,可按下式计算[21式中n—工作辊转速(m;)s/一‘(yP一rn.)勾岭(2-4)(2-5)第2章板厚控制基础n。一人工零位时的工作辊转速(;)s/mk—常数,决定于轴承的结构参数,由试验确定。图2-3示出了油膜厚度随工作辊转速和轧制力变化的关系曲线,可见随着压力的增大,油膜厚度减小,随着转速的增大,油膜厚度也相应的增大。£氛、叱轧制力减小方向n/rpm图2-3油膜厚度与工作辊轧速及轧制力示意图Fig.2-3 Sketch map of oil film thickness about work roll speed and press2.2轧件的塑性变形与塑性方程在轧制过程中,轧件在轧辊间承受轧制力的作用而产生塑性变形。与    此同时,轧件也给轧辊以大小相等的作用力(也称为轧制力))122-231。在轧辊直径及圆周速度相等、轧件机械性能均匀的情况下,轧制过程相对轧件的水平中心线是上、下对称的。如果除轧辊对轧件的作用力外无其它的作用力,则根据轧件的平衡条件,两轧辊作用在轧件上的作用力尸必大小相等、方向相反、作用在同一直线上,并且按对称条件,力的作用线必与轧辊中心线相平行。以上说的只是一个理想情况,实际轧制时轧件的受力情况非常复杂。    影响轧制力的因素主要包括以下四个方面:(    1)相对压下量在其它条件一定的情况下,随着相对压下量的增加,轧制力增大.    (2)轧件的塑性变形抗力金属的塑性变形抗力是度量其抵抗变形能力的力学指标。对轧制而言,随着轧件的塑性变形抗力增加,使其产生塑性变形所需的轧制力增大。(    3)张力采用张力轧制是冷轧板带的工艺特点之一。采用张力轧制使燕山大学工学硕士学位论文轧制力显著降低,并且张力愈大,轧制力愈小。不论前张力或张力均使轧制力减小,但后张力的作用效果更大。(    4)接触表面摩擦系数轧辊与轧件间接触表面的摩擦系数对轧制力的影响也很大。摩擦系数愈大,轧制力亦愈大。在一定的轧件宽度和轧辊半径条件下,轧制力P是所轧带材宽度B,    来料厚度H,出口厚度h,摩擦系数P,前、后张力今几以及变形抗力QS等的函数P                   = P(B, H, h, ju,T f,Tb, v,二〕(2-6)当把除h以外的所有变量视为常量时,轧制力P与轧件厚度h的关系曲线成为轧件的塑性变形曲线,如图2-4所示,其关系式称为轧件的塑性方程图2-4轧件塑性变形曲线Fig.2-4 Spring curve of workpiece表示为(2-7)尸二O(                         h)轧件塑性变形曲线的斜率,称为轧件的塑性刚度系数,可表示为入M=一二了仁71                门8P        (2-8)式中K、一轧件的塑性刚度系数(Nlm)e轧件塑性在一般轧机使用范围内基本是线性的,方程可近似表示为P=坑(劫+a󰀀)(2-9)第2章板厚控制基础轧件塑性刚度系数K、反映了轧件的软硬程度,它与一系列因素有关,在塑性曲线的近似直线的范围内,可表示为KM=馆刀                        (2-10)式中16一轧制力直线段延长线与横坐标轴的夹角。一般说来,坯料厚度H和张力T愈小,摩擦系数k愈大和变形抗力v,    愈大,塑性系数KM愈大。2.3轧机弹塑曲线将弹跳曲线与塑性曲线绘制在同一图上,即可得轧机的弹塑曲线图    印-h图),两曲线交点O的横坐标即为轧后的轧件厚度h。其中曲线A由式(2-4)确定,称为轧机的弹跳曲线,曲线B由式(2-7)确定,称为轧件的塑性曲线。曲线A的斜率即是轧机刚度KS,曲线B的斜率KM称为轧件的塑性刚度,它表征了使轧件产生单位变形所需的轧制力(1-51    曲线A, B的交点O决定了实际轧制力的大小,从而决定了该道次的出口厚度h,因此,一切影响曲线A, B交点位置的因素,都将影响轧件的出口厚度。板带材的厚度控制就是力图在交点位置变化了的情况下,维持交点的横坐标h不变。图2-                             5轧机P-h图Fi                      g.2-5 Figure of Rolling Mill P -h2.4  ¥L件厚度波动的原因影响轧件厚度的因素很多,主要有空载辊缝、轧机刚度模数,轴承油    膜厚度等(241。因此,无论是分析造成轧件厚度波动的原因,还是阐明板厚控制的基本原理,都应从这几个方面入手[171(    1)空载辊缝的变化轧辊的偏心、磨损和热膨胀等都会使实际的空载燕山大学工学硕士学位论文图2-6空载辊缝变化对轧件厚度的影响Fig.2-6 Influence in gauge due to gap vary辊缝so发生变化,从而使轧件的厚度产生波动.如图2-6所示,当空载辊缝S。变化Sol或Sot时,弹跳曲线的位置将由A移到At或A2,轧件厚度将由h变化到hi或h2-(    2)轴承油膜的变化有弹跳方程(2-4)可见,轧辊轴承油膜厚度的影响机理,与空载辊缝变化对轧件厚度的影响机理是相同的。随着油膜厚度的增加,轧件厚度变薄.    (3)轧件纵向刚度模数的变化在轧制过程中,由于轧辊的磨损和热膨胀沿辊身长度方向分布的不均匀,将使辊间的接触状况发生变化,造成辊系的弹性变形量波动,即轧机的纵向刚度模数发生变化。另外,轧件变形抗力的波动,也会通过影响变形区工作辊的弹性压扁,而使轧机的纵向刚度模数发生变化.如图2-6所示,当轧机的纵向刚度模数增加,弹跳曲线相应地由A变为A3时,实际轧出厚度h减小。可见,提高轧机的纵向刚度模数,有利于轧出更薄的板带材。在实际轧制过程中,以上诸因素对轧件实际轧出厚度的影响不是孤立    的,而往往是同时起作用。所以,在进行厚度控制时,必须综合考虑各因素的影响。第2章板厚控制基础2.5纵向厚差方程轧件轧后的纵向厚差方程口7    一,9],是建立各种厚度控制方式数学模型的基础。对弹跳方程式(2-4)两边取增量得Ah一AS.+竺一△:方程KS          (2-11)将轧制压力方程式(2-6)按泰勒级数展开,忽略去高次项可得轧制力的增量(2-12)”一aph    △A“h ・+                 P“・aB            AB+           8aHAp研・H+     apf<Aft +...将式(2-12)写成如下形式△尸二一K󰀀, A万十△F(2-13)式中△if,二长AB +于-AH +于- AU +...aB                aH aft‘,  、一_ap_ap__apKA,为塑性曲线在工作点处的斜率,AF是外千扰量BH, ,u, v:等引起的轧制力变化量.将式(2-I1)和式(2-13)联立可得昭。一△s十(2-14),.KM偏1甲.KK,空载辊缝变化量AS。可表示为战=                   -Ax,一AX,一Ax,十气(2-15)式中Ax,一压下位移增量(m),使辊缝变小时为正,与ASo符号相反;帆一支撑辊的偏心量(      m),偏心在工作辊一侧时为正,反之为负;      AX,一轧辊热膨胀量(m),热胀(使辊缝减小)为正,反之为负;Ax      。一轧辊的磨损量(m),恒为正;将式(2-15)代入式(2-14)可得一公,一公。一Ax‘十Ax。十,.KM1,尸一(2-16)燕山大学工学硕士学位论文式(2-16卿为纵向厚差方程。2.6各种因素对纵向厚差的影响系数由(    (2-16)式可以清楚地看出各种因素变化对纵向厚差的影响,并可定量地进行分析和计算【”]。在式(2-16)中,令Ax, =Ax, =Axm=os=Ax;二o,可得压下位移增量4x,对纵向厚差4h的影响系数  KSKS + Km(2-17)同理,可得4x, ,  4x, ,  Ox.和△氏对4h的影响系数  K,sKS + K,,,(2-18)(      2-19)    K,sKS   + KM将式(2-12)与式(2-16)联立,忽略掉其它因素对板厚的影响只考虑△万对Oh的作用,可得aP          竺=axOH  KS + KM(2-20)同理可得△Ta . OT,、如和△。,对Oh的影响系数翌鳖+K丝些+K竺玺+K护「巧一+K ?一城 从一崛从-如从一奴(2-21)        一-----一-  -M       一 衬      -K一K-K一K(2-22)(2-23)-A,-材(2-24)第2章板厚控制基础上述公式中轧制力对各因素的偏导数拳、头、任冬、景 ̄、于、.子二可___一,.-・-心,,--一,一r一ah    ax aTo  OT, ap  aas根据相应的轧制力公式求得。通过式(2-17)一(2-24)可得各种因素对纵向厚差的影响,.        、…_.,,_.。__。,二。-. ap  ap  ap  ap  ap  ap,2.7厚度控制的基本方式调整压下是厚度控制的最主要和最有效的方式,它通过改变空载辊缝    的大小来消除各种因素的变化对轧件厚度的影响[21。令△(5 =0,则由式(2-11)与式(2-13)联立可得[16-201Ah=一三s(SoA    KS+K",(2-25)式(2-25)相应的结构图如图1-2所示。图2-                        7调整压下控制结构图Fi                g.2-7 Control block diagram of adjusting press当△F,中的某些因素引起出口厚度变化Ah时,可通过调整压下来改变    辊缝,即让AS。二((1 +-" )Ah,即可维持厚度不变。由此机理构成的厚度控Ks’                    K,.。,_,                          ,,二 ̄一_一,、L,_工二二“,一一制系统成为自动厚度控制系统,简称AGC。通过式(2-17)^-(2-20)可得各种因素对纵向厚差的影响,采用调整压下方法就能对板厚进行精确的控制。2.8轧机等效刚度    对于制造好的一台轧机来说,轧机的结构确定后,轧机固有的刚度Ks是基本不变的,因此,把轧机的固有的刚度称为轧机的自然刚度[251    在液压压下的轧机中,由于液压压下具有调节速度快,反应灵敏等优点,使得因轧制力等因素波动而产生的轧机弹性变形可以通过调整辊缝及时得到补偿。这种情况下,虽然轧制力的波动可能很大,但纵向厚差可以燕山大学工学硕士学位论文控制得很小,这相当于提高了轧机的纵向刚度模数。生产中,对于因坯料厚度波动或其机械性能不均等因素所引起的板厚偏差,要求轧机刚度越大越好,刚度越大,则出口板厚偏差越小。实际上辊缝是在不断变化的,通过自动控制系统不断调节液压压下,消除轧机的全部弹性形变,使带钢出口厚度不变。由此引出了轧机的可变刚度K。的概念。不考虑油膜厚度的影响(AS=0),由式(2-11)得(2-26)Ah=ASO+州KS                  轧机固有刚度可定义为人。=- L△尸      Ah‘        一-(2-27)式中K:一轧机固有刚度(N1m);      AS—为空载辊缝增量(M);AP一为轧制力波动(      N);Ah一为轧机负荷辊缝变化(      M).对于具有快速响应特性的液压压下系统,变刚度控制系统通过对轧制    力的测量,控制液压压下系统动作,在一定程度上补偿了轧制力变化引起的轧机弹跳,相当于改变了轧机刚度,称为轧机等效刚度Ks. 4.式中AP—轧制力波动;Ks                       ,, =△习Ah                    (2-28)Ah一为负荷辊缝变化      在空载辊缝增量方程式(    2-15)中,令Ax, =AX, =AX. = 0,即只考虑压下位移增量的影响,则有ASO                         =-Axp               (2-29)设a为补偿系数,压下位移增量Ax,按下式给定△x,=a-△尸          KS        (2-30)则将式(2-26)、式(2-29)和式(2-30)联立,可得调整压下后的纵向厚差为(2-31)而定义等效纵向刚度模数[3]第2章板厚控制基础。=.K,,KS1一a(2-32)    轧机等效纵向刚度模数K,表示轧机厚度控制能力,从式((1-27)可知,改变补偿系数a的大小,就可改变等效纵向刚度模数Ks,,实现不同的控制方式【17-191. a与Kss的关系见下表表2                  -1 a与KSA之间的关系Ta                  ble 2-1 Relation between a and KSAa取值口=1  KS  A值KSA=.控制方式恒辊缝控制硬特性  0<a<1口二0  口<0a二一刀KSA>KSKSA =KSKSA<KSKSA=0自然刚度控制软特性控制恒压力控制    实践证明,要使补偿系数a=1是很困难的。为了使控制系统稳定,目前可达到的是a=0.9左右,即轧机等效纵向刚度模数最大等于其自然刚度模数的10倍。采用硬特性控制时,使纵向厚差减小,但使板形变坏。恒辊缝控制(    最硬)使纵向厚差最小(为零),但板形最坏。采用软特性控制时,使纵向厚差加大,但板形变好。恒压力控制(最软)使纵向厚差最大,但板形最好。一般在轧制的前几个道次采用硬特性控制,在轧制的后几个道次采用恒压力控制,以达到平整轧件,改善板形的目的。2.9本章小结    本章首先从机座的弹跳方程和轧件的塑性方程入手,运用轧机的弹塑曲线详细分析了轧件厚度波动的原因,通过纵向厚差方程的推导得出了各种因素对纵向厚差的影响系数,最后对变刚度控制方法的原理进行了剖析。燕山大学工学硕士学位论文第3章1270mm冷带轧机的组成及控制模型3.1液压厚控系统在20世纪60年代,液压伺服系统开始得到广泛应用,当时它只用于    在一些低负荷的系统之中,这些系统的自然振荡频率相对来说都比较小,对液压系统的要求也不高[261。直到70年初液压伺服系统才被应用在轧机板厚控制系统中。随后20年的实际应用证明,在板厚控制系统中引用液压伺服系统时非常成功的。一个完整的液压板厚控制系统一般有六部分构成:    (    1)电液伺服阀电液控制系统中将电信号转换为机械信号的元件以力矩马达应用最广,力矩马达与液压放大元件一起,组成电液伺服阀。是系统中的控制元件。(    2)传送管路传送管路如果太长,系统的滞后作用将加大。可通过伺服阀与液压缸的直接安装来消除。    (3)液压缸液压缸是伺服系统的执行元件.在系统中的作用是迅速、准确地复现伺服阀的动作〔系统输入)。(    4)轧机包括轧机的动态特性,自然振荡频率,刚度和阻尼效应。在系统中属于被控对象,即负载.    (5)回路即泄路,产生备压,在压力闭环中必须考虑。(    6)传感器它包括位移、压力和速度传感器等测量设备。    在现代轧机中,板厚板形控制系统的主要环节均采用电液伺服控制系统,这主要是由于电液伺服控制系统较之机械、电气元件做执行机构的系统具有明显的优点.主要表现在[271.(    1)液压系统结构简单,对于轧机设计,尤其是对旧轧机的改造,带来方便。(    2)力矩惯量比大,使其在同样输出功率的条件下较机械和电气机构具有更小巧的结构,容易实现。(    3)液压系统惯量小,相对来说响应速度快,这就为高速轧制的快速跟踪调节提供了可能性。鉴于此,唐山恒通公司决定将引进的二手轧机原始的电动压下,改为    第3章1270mm冷带轧机的组成及控制模型液压压下。我有幸参加了电控部分设计、安装和调试工作.以下将对1270mm冷带轧机的组成及工作原理做详尽的介绍。3.2  1270mm冷带轧机的组成与工作原理1    270mm冷轧机的总体组成如图3-1所示。控制系统主要包括:上位机、下位机、模拟量控制部分、操作台及外部设备等。电气控制系统的主要作用是通过1270mm冷带轧机进行带材厚度的控制128-291控制}              叫夏刃勤电动压下#Ir"LI I}}  Di人}纂匕-+}fik一凡        f1AR恻一}9L}}模1}_尸侧厚仪一之斗        卿嘶一赘3era       Te&jx7 I T 下罕h」K乡辊CTf叫?I飞\厂一、卷取机¥r}园4传动图3-          1  1270mm冷带轧机控制系统组成图Fig.3-1 Making up of the control system for 1270mm cold strip3.2.1厚控系统硬件的总体结构和主要功能    上位机是操作者与厚控系统进行人机对话的窗口,用于系统运行状态的监视、原始数据的输入以及轧制规程的计算等。系统监视的状态主要包括:缸位移值、轧制压力、轧件出、入口厚度、轧制速度以及工作状态等;输入数据包括:钢卷号、材质、带材宽度、入口厚度、成品规格等。下位机实施各种控制,包括BISRA补偿、偏心补偿、油膜厚度补偿和加减速补偿以及厚度监控等,并上传显示数据.上下位机通过RS232标准串行口进行通讯。模拟量控制部分主要完成各种现场信号的预处理及实现位置和轧制力闭环控制:执行完成摆辊和压靠的功能;接收和向计算机传送数据;燕山大学工学硕士学位论文等等。模拟控制部分的设计质量将直接影响控制系统的稳定性,可靠性以及控制的精度。作为现场操作人员对系统进行全部功能控制的操作台,集中了系统的    全部监视功能和控制功能。台面上设有手动控制和自动控制功能转换开关。手动状态下,各参数设定值由操作人员手动控制来完成;自动控制状态下,操作人员只要在上位机中输入带材的各原始参数,即可由计算机自动完成设定。此外操作台还安有厚度环接入与切除按钮,方便维修和调试。操作台安装有摆辊、压靠、四位手柄开关以及厚度拨码设定开关等。摆辊功能用于预摆初始电动压下位置,四位手柄开关的四位分别为同上、同下、左上右下、左下右上,正常轧制时操作人员可以根据监视的轧制状况手动微调辊缝位置。厚控系统的总体结构框图如图3-2所示图3-                      2厚度自动控制系统总体框图Fi                  g.3-2 General block diagram ofAGC system厚控系统的原理框图如图3-    3所示。该系统在控制方式上使用了二套的控制回路:上位机沁叫下位机伺服阀卜叫轧机卜叫测厚仪缸位移轧制力显示图3-                          3厚控系统原理框图Fi                        g.3-3 Principle diagram for AGC(1)缸位移闭环一S环(2)恒压力控制闭环一尸环第3章1270mm冷带轧机的组成及控制模型    操作者可以根据轧制工艺选择不同的闭环形式进行轧制。S环精度(1u m),调整范围((0-12.5mm)a P环工作范围为(0-5000)KNa缸位移闭环闭环一般在轧制的前几个道次使用,压力闭环主要在压靠时和轧制的后一、二个道次使用,以用于消除轧辊偏心等因素的影响。图3-    2总体结构中的基本通讯关系如图3-4所示,图中虚框内即为模拟量闭环控制部分,从组成单元上划分,它有3个部分,即快速模拟量反馈控制系统、操作台以及现场操作箱。在正常运行状态下,上位机、操作台以及现场操作箱是操作者实施控制的界面,在操作台上可完成系统全部的基本功能,快速模拟量反馈部分是控制系统的核心。图3-4厚度控制系统各单元通讯关系图Fig.3-4 Communication for unit ofAGC    将图3-3中的模拟量闭环控制系统分离出来并具体化如图3-5所示。正常运行状态下模拟量闭环控制系统的给定量(对应二种闭环形式的缸位图3-          5模拟量闭环控制系统原理图Fig.3-5 Principle diagram of analog closed-loop control system燕山大学工学硕士学位论文移量和压力值)是由下位机通过D/A给出的。当然根据需要也可以由模拟量通过拨码盘给定。对应上述二种闭环形式,实际的输出分别由缸位移传感器和压力传感器检测,将实际反馈量与给定量比较,由相应闭环状态的调节器进行控制,经功率放大器后驱动伺服阀,使缸位移值或轧制力保持与给定值相同。可见,模拟量闭环控制系统完成的是一个自动位置控制或自动压力控制。除上述主要功能外,还自动完成来自下位机的弹跳量补偿以及手动控制方式等。自动控制状态下,厚度监控补偿量是由下位机直接在内部运算后加到响应给定值上的。3.2.2厚控系统主要控制方法随着对高精度板带材的需求日益增加,原始的电动压下系统已越来越    不适应市场需求。基于电液伺服系统的液压板厚控制系统以其控制的高精度和快速性在轧钢领域受到了广泛的应用。因此1270mm冷带轧机厚控系统采用液压压下来代替以前的电动压下,即采用液压AGC系统.液压AGC系统主要用来清除带材纵向厚度偏差,提高成材率,改善自动化程度。它具有以下特点:惯性小,响应快,精度高,安全可靠,操作方便等.厚控制系统的主要通过调节液压压下改变辊缝来达到控制的目的。液    压压下是在电动压下基础上进行微调,利用电动压下的大行程对液压压下系统进行预摆辊缝初始位置,使液压微调尽可能在最小的液柱高度上进行控制,以提高液压系统刚度。电控系统部分主要由下位机指挥,将来自现场的各种数据进行综合、处理之后将控制信号经D/A输出,经功率放大后驱动液压伺服机构,对液压缸进行位置控制,进而实现对带材厚度进行控制[28-291    1270mm液压冷带轧机AGC厚度自动控制系统组成主要有:(    1)手动控制闭环功能状态置于手动控制时的原理如图3-6所示当“压下”按钮按下时++15v电源通过W, Pi和凡构成回路,“。点具有一正电位,经功率放大的后驱动伺服阀使轧辊压下:反之“抬辊”按钮按下时,-15v通过凡.码和R3构成回路,u。点为负电位使轧辊抬起,压下与抬起的速度分别由W,. W2调节.(    2)压力闭环控制当闭环功能状态量置于压力闭环(简称尸环)控制时的自动控制原理图如图3-7所示,给定压力P是两侧压力之和,将尸值除第3章1270mm冷带轧机的组成及控制模型图3-                            6手动控制原理操左侧传动侧Fi                g.3-6 Control principle for manual operation以2后分别送入两侧的压力闭环控制系统。轧制力和几(备压)由压力传感器检测,并将经二次仪表处理后的输出信号送入轧制力仪表板(简称尸仪表板).P仪表板对两侧的信号进行处理之后,分别反馈到各自的PID调节器输入端,形成两侧各自的压力闭环控制。当P轨迹为零时,自动调节的结果应为PI=P2=1/2P.在调整板形需要的情况下,调节P轨迹,可以人为使两侧的压力不同,但始终保持PI+P2=Po P轧制调整可在操作台和现场操作箱上进行。仪表板将尸、,P:处理后的信号,一方面送入计算机的A/D,另一方面对二者求和P=PI+P2、求差△P=PI -P2,进行反馈并显示。操作侧计给模给传动侧图3-                        7压力闭环控制原理图Fi              g.3-7 Principle diagram of pressure close-loop control    压力闭环的主要作用是在摆零时将压力给定由手动或自动设定到某一值,对轧辊进行压靠,然后对两侧的辊缝值和缸位移值进行拨零,作为辊缝和缸位移设定的初始值。有时也可以在压力闭环状态下,进行带钢的平27                                          燕山大学工学硕士学位论文整,或在轧制的最后一两个道次,改为压力闭环轧制,以消除偏心等影响,提高带钢的板形指标。(    3)缸位移闭环控制当闭环功能置于缸位移闭环(简称s环)控制时的自动控制原理如图3-8所示。S环控制原理与压力闭环基本上相同,这里的缸位移给定值S由下位机通过D/A给出,即最优轧制表的缸位移给定值。两侧的自动位置控制系统均以S为基准进行闭环控制。S,    , SZ是指系统摆零之后的实际缸位移值,但为了充分利用缸位移传感器的可测范围,一般摆零时缸位移二次仪表输出总是偏向正或偏向负,不一定正好在零位,经S仪表板检测后送入下位机,下位机给出一平衡信号,使摆零时两侧的偏差为零,此后,再调整S轨迹时可使S, :A凡,产生操作侧计境模给传动侧图3-                  8缸位移闭环控制原理图Fig.3-8 Principle diagram of hydraulic cylinder position close-loop Control人为的缸位移偏差,以手动控制板形。缸位移闭环的位置控制范围较大,正常可达1    2.5mm左右。当来料带钢较厚时,一般在轧制的前几个道次均采用缸位移闭环进行轧制,后几个道次采用压力闭环。缸位移闭环是整个厚度控制的基础。(    4)测厚仪监控系统在轧机出口侧直接用测厚仪测出轧件厚度,然后反馈调整压下装置,改变空载辊缝,消除厚度偏差。因为压下位置闭环控制和轧制压力变化补偿,都不能消除轧辊磨损、热膨胀对空载辊缝的影响以及位移传感器和压力传感器本身的误差对轧出厚度的影响,因此系统中必须加有测厚仪监控系统。测厚仪监控系统位于系统的外环,属于前馈补偿控制[(2)92                                          第3章1270mm冷带轧机的组成及控制模型由2.5节式(2-17)知可得厚差4h与压下调节量dx,,之间的传递函数(3-1)式(    3-1)与(2-17)相差一个负号,原因是式(2-17)中的Oh是AXp所造成的厚度变动量,而式(3-1)中的4h是要消除的厚度差。要使h减小,x;应该增大。从式(  (3-1)可知,要消除厚差劫,必须将压下装置移动((1十且)劫的距______二_、.K.__一KS’            离,即要移动比厚差Ah大共sAh的距离。KM                        由于空载辊缝改变量ASo =-dxp,所以式(2-17)可改写为叮=硫从    KSKS + Km(3-2)式中17一压下效率系数。从(3-2)可见,当空载辊缝改变AS。时,它所引起的实际厚度的变化量Ah要小于ASo。这说明,压下位置的改变量没有完全反映到轧出厚度的变化上,而其中的一部分表现为轧机弹性变量的形变。当轧件变形抗力口:较大而轧件厚度h很薄时,由于轧件塑性刚度系数KM很大,则压下效率很低。此外,考虑到轧机结构的,测厚仪一般安装在离辊缝较远的位置,    通常L=150-1750mm。因此,厚度测量点与轧制点不在一处,滞后时间为(3-3)式中公—测量滞后时间((S);1一轧辊中心线到测厚仪的距离(    m);v一轧制速度(      m/S)e由于有时间滞后,所以采用比例控制很难保证控制系统的稳定性。而    要保证稳定控制,则需要降低放大倍数,这将影响稳态精度。通常为提高系统的稳定性,测厚仪监控系统是断续采样的,有的还设    置死区,以避免系统频繁动作。采样周期根据滞后时间和系统时间常数确定。由于存在时间滞后,难以纠正因入口厚度骤变而引起的轧后厚差。它只能纠正低频、缓变的偏差,并且必须经过几次传递时间之后,才能得到燕山大学工学硕士学位论文                          .目二二二二....二......二.....二纠正[21。因此它位于系统的外环,用作厚度监控。(    5)倾辊控制为了保证板形良好,在来料为楔形情况下,要求轧辊两侧具有一定的倾斜,即所谓的倾辊控制。常规的倾辊控制是将两侧的缸位移偏差值与倾辊设定值比较后进行PID控制,将PID输出结果以互为相反的值分别作用于两侧伺服放大器1151。这种情况下,两侧的位置控制采用的是同一个PID控制器。本文作者采用的是具有作用的倾辊控制方法,其控制结构图如图3-9所示。图中‘Sx.为倾辊控制设定值,Sx’为缸位移给定值,C. I Cd分别为操作侧和传动侧控制器。该控制模式是采用两个各自的PID控制器,分别控制两侧的缸位移值,而将倾辊设定值互为相反地加入到两侧PID控制器的输入端。这种控制模式的特点是调整两侧图3-9倾辊控制结构图Fig.3-9 Block diagram of tilt roll control的PID控制器参数时互不影响,当两侧轧机特性不完全相同时,仍可通过调整两侧PID控制器参数使两侧闭环响应特性大致相同。而常规倾辊控制模式,两组PID控制器参数调整相互影响,给实际调试带来困难。当采用作用的倾辊控制模式时,系统还应设置大行程同步控制器,    系统两侧的压下状态总是受到大行程同步控制器的监测,当两侧压下偏差超过设定值时,则启动大行程同步控制器,两侧偏差最大值。而在正常轧制两侧压下偏差较小时,同步控制器是不工作的。(    6)油膜厚度补偿本轧机支承辊为油膜轴承,油膜厚度变化对辊缝影响较大,其值与轧制压力和轧制速度有关,为补偿这部分影响设置了油膜厚度补偿环节[30-321。在轧制过程中根据轧制压力和轧辊转速的实测值,计算出油膜厚度的变化量,然后调节压下对油膜厚度的变化单独进行补偿,消除其对轧出厚度的影响。30                                          第3章1270m。冷带轧机的组成及控制模型在纵向厚差方程(    2-16)中,令Ax== Ax,=Ax.=4F=0和Ah=0,可得补偿油膜厚度变化所需要的压下位移调节量为Ax                      ,=一△S                   (3-4)负号说明,如果油膜厚度增加,则压下位移量应该减少。上式中的△s是实际轧制条件下相对设定轧制条件的油膜厚度变化量    可以有式(2-5)时得到,实际轧制时通常用折线代替图2-3中的曲线,并以表格的形式存放在内存储器中。(    7)轧辊偏心补偿支承辊的偏心对空载辊缝和轧后厚度造成周期性的影响。而且当轧后厚度增加时,轧制力减小,当轧后厚度减小时,轧制力增加,这将会引起轧制力闭环系统误动作,进一步加大厚差,所以需要单独给予补偿[21。这一部分在后面还要着重介绍,这里不再详细论述。    该厚度控制系统的最大特点是在普通缸位移闭环基础上,设有倾辊控制、油膜补偿控制、增强型PID位置自动控制算法、偏心补偿控制等功能。有效地改善了轧机的动态响应性能,提高了板厚控制的精度。该控制系统采用了双闭环的控制方式,外环由秒流量,前馈,测厚仪    监控,速度、偏心,弯棍力等补偿算法构成。外环各功能主要有由下位机来实现的,即有下位机接受内环模拟部分传送来的数据,经过下位机的运算,然后将输出值再传给内环模拟部分,有这部分来实现具体的控制功能。内环采用模拟量部分包括压力闭环和位置闭环,在轧制时根据具体情况进行选择。轧机的执行机构是伺服阀,由它去控制轧机的压下装置,其中直接控制部分的控制量和反馈量均为模拟量.由于控制的精度很高(/    en级),要求下位机的工作量不可能太大,并且控制系统的响应速度必须很高。由下位机和模拟量控制系统共同构成的控制系统则可以很好的把这种优点显示出来,很好的完成控制的要求。在本系统中,模拟量控制系统是它的最为基本的环节,是本设计的特点,它在控制中的作用不可替代。3.3软件的设计思想,来完成板厚控制系统的各种    软件根据硬件电路的特点和设计的要求功能。下位机采用了c语言进行编程控制。它完成的主要功能有摆辊,压靠,压力闭环控制,位置闭环控制,厚度闭环监控,油膜厚度监控,偏心燕山大学工学硕士学位论文图3-1          0厚控程序整体流程设计Fig.3-10 Design of thickness program about the whole lfow补偿,状态采集,与上位机串行通讯,A/D采样等。如图3-10所示,它描绘出了厚控程序的大致框架。下面将对厚控程序所完成的主要设计思想进行简单的介绍[33-341(    1)初始化包括给定值的初始化,通讯端口的初始化(传输速率、奇偶校验等),轧制压力值的初始化,AM卡的初始化,D/A卡的初始化(轧制力和缸位移在初始化时的设定值是不同的)等.(    2)摆辊程序由于液压缸的最大缸位移只有25mm,在压靠时,或者在穿带后进行轧制时,可能存在这种情况:当缸位移达到最大值时,仍不第3章1270mm冷带轧机的组成及控制模型能满足压靠所要求的压力值,或者在轧制时,不能很好的完成所要求的压力闭环或位置闭环控制。原因在于:在轧制前,没有合理的设定初始辊缝值和缸位移的初始值。例如:在轧制的过程中,每一道次之后,压力值或者缸位移值都要增大,此时,极有可能出现这种情况,缸位移己达到了最大值,仍不能满足压靠或者轧制时的要求,这必将对生产造成严重影响。因此有必要提出摆辊的概念。    摆辊的思想是:在压靠前,或轧制前,通过计算机设定一个合理的值,采用电动压下使整个缸体向下移动。使刚体稳定在某一位置。在这个位置,通过液压控制,可以使缸位移在给定的工作行程内伸缩自如,即液压缸完全缩回去时可方便的完成穿带:当使工作辊完全接触时,仍留有一定裕量,以便实现所要求的最大压力值.在现场时,首先使缸体伸出12mm(即为所要留下的摆辊辊缝值),然后采用电动压下,使其下移到两个工作辊刚刚接触,然后令缸体收回,完成摆辊。可见,当液压缸完全收缩回去,12mm的间隙足以穿带;当工作辊完全接触时,仍留有很大的裕度,足以完成压力闭环或者缸位移闭环的要求。    摆辊时,操作台按钮必须置在位置闭环(若为压力闭环,则报错),非压靠状态(若为压靠,则报警),系统设定状态。这就要求程序能够不断的检测按钮的状态,并能够及时预报,提示错误信息。满足条件后,才可以进行摆辊.若检测到停止摆辊的信号,应停止并返回。在状态检错部分,当操作台按钮不是所要求的状态,向上位机发送错    误代码,再由上位机完成错误提示,以及解决问题的办法。需要说明的是,关于错误代码的解释,在编程时,上位机和下位机之间己经有了定义。为了减小下位机的工作量,下位机只传送错误代码。摆辊时,在S仪表板上,必须不能有压靠值,需清零;然后要使缸体伸出,达到所设定的值。在摆辊之前,可能存在这样的情况。如果电动压下己将液压缸压下很    多,即所留下的辊缝值小于12mm(缸体完全收缩回去时),比如说:只有7mm。这是如果进行摆辊,由于7mm不会给缸体留下足够的伸缩空间,当工作辊接触时,刚体仍有5m,没有伸出。这样的话,液压缸势必将一直向下伸出,工作辊上必然承受很大的压力。当这个压力超过了一定的范围,会对工作辊和轧机造成重大的损坏,出现事故。解决问题的方案就是:在液压缸向下伸出的过程中,不断的检测液压    燕山大学工学硕士学位论文缸的压力信号(于工作辊之间的压力相接近),当两侧的压力之和超过了一定值时,就进行事故处理程序。这时,要完成的工作有:停止摆辊,使液压缸完全收缩回去;向上位机发送错误信号,并由上位机提示错误信息和解决方案。(    3)压靠程序进行压靠是为了消除弹跳曲线中非直线段的影响,更好的利用用纵向厚差方程,可参考2.1.2节中的人工零位法。在压靠时首先要保证轧制状态为手动和压力闭环;否则的话将进行错误处理。然后从拨码盘上读入压靠值,让程序检验压靠值的设定是否合理,若不合理,必须进行错误处理,然后重新设定,待条件满足后就可进入压靠程序。在装上工作辊后,即第一次用新辊进行轧制时,必须进行压靠,储存    压靠值;换辊后,也必须进行压靠,这种情况下必须注意刷新压靠值。在轧制的过程中,虽然时常要进行抬辊,穿带,换卷等操作,只要仍然是原辊进行轧制(没有更换工作辊),就可不必进行压靠操作,只要从计算机中调入己知的压靠值即可。    进行压靠时,压力值在压力闭环中给定完成,此时计算机将不断的检测位移传感器的值(压靠值),并进行数字量反向反馈(清压靠值)。(这个功能在S仪表板上实现)。当轧辊的压力达到所设定的压力时,通过计算机的回馈作用可以使缸位移的闭环输出为零(此时的压靠值就是人工零位点)。但是,压靠的测量值不可能十分稳定,有一定程度的波动,将此值用A/D卡采入计算机后,转化为整形量,又将从D/A卡输出,送到压靠值清零接口。这样的话,就避免不了出现这种情况:计算机的输出始终在跟踪所测量的压靠值,即使是在十分稳定的情况下,这两个值也几乎是不可能完全相等的。何况传感器还存在精度问题,A/D卡和D/A卡还存在量化的误差等等。因此,必须容许两者存在有一定范围的误差。因此要编制相应的误差处理程序。在压靠的过程中,控制程序是在不断的进行循环,误差允许处理程序在每一次循环中都将被运行。在进行动态的比较之后,可以使两侧的压靠值在一定的范围之内波动。(    4)位置闭环控制厚制系统采用了计算机和模拟部分联合控制的方法,闭环功能主要有模拟环节来实现。这样设计,一是因为模拟量的控制响应速度比较快,可以保证轧制的精度,二是为了减小下位机的工作量,保证下位机具有足够快的处理速度。因此,在位置闭环控制中,程序主要第3章1270mm冷带轧机的组成及控制模型完成的功能有:向上位机传送厚度设定值,轧件出口厚度值,缸位移值(用以绘制曲线);检测手柄和快速手柄的状态,并做出响应[3513.4板厚控制模型的建立3.4.1液压伺服系统动态模型的建立    自动厚度控制(AGC)系统的内环是由液压压下系统组成的自动位置控制(APC)系统。o APC系统由位移传感器测量液压缸位移量构成反馈闭环控制.APC系统作为AGC系统的内环,其性能指标直接决定着AGC的指标。1270mm可逆冷轧机单侧液压压下APC系统的组成如图3-11所示。位里给定一控制器曰1功放卜洲伺服阀+              位里反馈图3-11                          APC系统示意图Fi                    g.3-11 Sketch map for the APC system下面分几步来求液压伺服系统各部份的传递函数,以最终建立其动态    数学模型[261(    1)电液伺服阀由机理分析和频率响应试验可知电液伺服阀的传递函数具有如下形式a  ),s2 +2凸J 2s式中C(s)一为伺服阀滑阀输入电流(MA)1      ,(s)一为伺服阀滑阀位移(M);。,一为伺服阀等效固有频率(      Hz)古—阻尼系数。      其中固有频率和阻尼系数由制造商提供,‘:)一、{旦〔丛二)-l 'D)        了t一}了  军-5(3-5)(3-6)燕山大学工学硕士学位论文式中1。一为伺服阀额定电流(mA);K。一为电流放大器的增益系数;            UX一位移变送器输出(my );U,一输入控制信号(      mv)e伺服阀的流量q与伺服阀的位移成正比,同时与其输入输出端口的压力差有关。。一。。0c_95O0P0式中P一伺服阀内液体压力(Pa);AP一液体压力差(    Pa);(3-7)q    。一伺服阀额定流量(m'1s)a(    2)传送管路传送管路的动态行为不能通过一个简单的时间差分方程来描述,它是关于时间和空间的函数。虽然用差分方程描述是最好的近似,但如果和轧机其它元件的动态特性相匹配仍有一些困难。根据传输管路的特点,我们用差分方程而不是偏微分方程来描述其特性。根据质量和动量守恒定律,可得以下二个等式。一q.一普Ps式中(3-8)q一管路内液体流量(m'/s);A,一管路截面积(m2);L,一管路长度(M);q。一进入液压内的液体流量(m'/s);八一液体等效容积弹性模量(Pa);s—微分算子。P一:一令qos+Rp・。式中几一液压缸内液体压力(Pa);(3-9)      R,一压力降系数;P一液体的密度(      kg/m');第一个等式是根据传到送管路中油的惰性得到的,第二个等式代表了油的摩擦损失。等效弹性模量,包括液体弹性模量、管路的塑性模量和阻塞气体的作用,可以通过下式计算第3章1270mm冷带轧机的组成及控制模型上_鱼+兰式中dpfl,    dp几管路的内径(m );E—        管路材料塑性系数;(3-10)管路厚度(M);tp—        被阻塞气体体积与总的液体体积比;气体的弹性系数s—微分算子            液体的等效弹性模量是随着压力而变化的(压力越大弹性模量越大)。气体的阻塞作用不可避免降减慢系统的响应速度。    (3)液压缸液压缸的压力由四个主要的因素决定:液压缸的流量90,泄漏流量9;,液体的可压缩性而造成的体积变化9。和活塞的运动速度9.=综合这四方面的因素,可得下式。一击〔。一,一K,(“-Pb)]式中A。一液压缸的截面积(脚,:)      K,一流量泄漏系数;几一液压缸的备压(Pa)      ;(3-11)Y。一开始液压的容积(      m3 );      x一液压缸的位移(M).因为液压缸的位移x在等式(3-10)中位于分母的位置上,可知式(3-10)是非线性。在现在的液压厚控系统的设计中,管路应做的尽可能的短,或者将伺服阀直接安装在液压缸上。这样就可以忽略由于管路摩擦造成的压力损失。(    4)回路由于回路的存在,轧制时必须考虑液压缸的备压。备压计算公式如下。一。。・PL.奇xs2・R,xs式中A,一回路的横截面积(m2);L。一回路的总长度(      M);Pe      。一初始备压(pa);(3-12)燕山大学工学硕士学位论文R,一压力降系数。(  5)位移检测环节包括传感器与二次仪表等,这部分的传递函数可表示为全up=K,4.X0兀S+1S);式中Tx一位移检测中引入的时间常数((3一13)Kj一位移传感器及相应二次仪表的放大系数(、而子)。设K,,=90儒,Kb—1宁爪pPLP.。AP综合式(3-5卜(3-13)Kd= R, + A0 "L,S ,得到液压伺服系统的动态模型如图3-12所示[201    图3-12对液压系统各部份特性都做了详细的研究,从动态模型的可以看出液压伺服系统是高度非线性的,但用它来研究液压缸的动态特性仍显不足,因为它没有把轧机部分包括在内,不包含控制功能,很难估计液压图3-1        2液压伺服系统动态模型Fig.3-12 Dynamic model of the hydraulic servo system伺服控制系统中控制器的作用和其动态特性。为了对液压板厚控制系统的性能进行精确的分析,建立一个切实可行的控制模型是非常重要的[[371第3章1270mm冷带轧机的组成及控制模型3.4.2板厚控制系统动态模型为了便于分析,如图3-    13所示[[261,对轧机液压伺服系统的动态物理模型作了适当的简化处理。伺服阀传递函数可等效为一个积分环节和一个r一一一一一一一一一一一一一一WT心今      K  IUT,S+1 I一I T,S,TtS匕-------------------_  --- j准11:图3-1  3常规型PID控制模型Fig.3-13 General PID Control model惯性环节,其时间常数分别为乃、T,;设位置环控制器为K,、伺服阀控制器为K3、伺服阀增益设为K4,伺服阀阀芯位置反馈增益设为Kf,惯性主阀流量增益设为系数K2,液压缸积分时间常数设为T2,液压缸位置反馈系数设为KF。在以上设定情况下,可得常规型PID控制模型。图3-1    3所示模型仍只是液压伺服系统的的模型,没有把轧机系统,各种干扰,各种控制算法和补偿算法考虑在内。将这些因素考虑在内可得液压板厚控制的一般物理模型如图3-14所示,图中把PID调节器和轧机控制器Ki放在一起,合为调节器K},伺服阀传递函数中的常量也做了相应的线性简化。得到的动态数学模型一般包括以下四部分[38-391(    I)轧机机架,待轧带材(厚度及塑性系数均可能有变化),来自带材自身因素以外的各种干扰。    (2)具有完备动力能源的调节装置和执行机构组成的液压系统,作为控制系统的内环.一般包括:液压缸、伺服阀、负载限定阀、和位置、压力变送装置。(    3)控制工作站,作为整个控制系统的给定和监控机构,是人机对话的窗口。通常指上位机,操作台,及一些外部设备。    (4)控制系统的外环,实施各种控制算法和干扰补偿,包括秒流量控制,测厚仪监控,BISRA补偿、偏心补偿、油膜厚度补偿和加减速补偿等方法。39                                      燕山大学工学硕士学位论文图3-1                            4液压板厚系统控制模型Fi                    g.3-14 Modelof hydraulic gage control system一般由下位机来实现。3.5本章小结    本章从液压板厚控制系统入手,对1270mm冷带轧机的组成,控制系统各部分的功能,工作原理,软件的设计思想进行了详细的分析,在参阅大量文献的基础上,推导出了轧机控制系统的动态数学模型,对精确分析液压板厚控制系统的各项性能有重要意义。第4章系统动态性能分析及改进第4章系统动态性能分析及改进4.1动态特性分析    合适的设计、控制方式的选择和优良元器件的选用是使液压板厚系统达到最优动、静态性能的关键。对液压板厚控制系统的设计标准和动、静态性能进行分析,设计和改进各自的系统以适应不同板厚控制项目的需要,是非常重要的[19]4.1.1车L机的负载设计特性    一个设计优良的控制系统首先应满足一定的负载特性。自动位置控制系统作为板厚控制系统的内环,必须对微弱变化信号具有快速响应的能力,而且对与大行程信号,动力部分应能提供足够的能源。即在轧制力和辊缝允许的变化幅度内,实时做出调整,以保证轧制过程的正常进行。这可以从其动、静态负载特性来进行分析[[40-42[静态负载特性一般是指辊缝的变化速度位于3mml    s和5mm/s之间时,系统的跟随性能。低速变化值,可以确保一定的板厚控制精度;上限值,可以防止超调过大,但相应地减弱了系统的响应速度.动态负载特性一般是指辊缝以一定的正弦规律变化时,系统的动态跟    随性能。辊缝变化的幅值和频率均应在一定幅度内;幅值变化可根据实际需要来定,一般大于0.5mm,小于0.75mm o频率变化设计时应大于0.2Hz,小于支撑辊最大旋转频率的二倍。低频变化值,也是为了保证足够的板厚控制精度,最大频率,是为了减弱轧辊偏心对轧制的干扰(轧辊偏心是一种周期性高频干扰信号)。4.1.2元件的选择对厚控系统性能的影响    自动位置控制系统作为板厚控制系统的内环,一般包括三部分[20[:控制元件、执行元件,能源供给部分和调节装置。控制元件指伺服阀,执行元件一般指液压缸、负载限定阀等,能源供给能力主要与伺服阀的额定流量有关,调节机构指调节器,以及位置、压力变送等装置。其中伺服阀、液压缸,压力变送器,位置变送器的正确选择,对于板厚控制系统的设计至关重要,在某种程度上它决定者设计的成败。轧机控制系统的响应性能41                                       燕山大学工学硕士学位论文直接和伺服阀、变送器和执行机构根随轧制力变化的能力有关。(    1)伺服阀的选择对系统性能的影响伺服阀当尺寸太大时,因为时间常数较大,将减弱系统的响应速度。尺寸太小,额定流量也较小,当流量较大时,又会出现饱和现象,不能满足所需。图4-1是当供油压力为Ps=21 Mpa,阻尼等于0.01时,伺服阀的阶跃输入分为0.025mm, 0.05mm,0.1 Omm(曲线1, 2, 3分别对应0.025mm, 0.05mm, 0.10mm),额定流量为631/min时的响应对比曲线。C六L、.0卜・-,・--一;一‘、〔.舍二二二t1--二...二二二,甲,一‘-------件〔I0                             0.05                                    0.1        0.15       0.2 t1s图4-1伺服阀响应曲线图                        Fi            g.4-1 Comparison diagram of step response of sevrovalve    由图4-1可见当阶跃输入信号太大时伺服阀将会出现饱和现象,从而降低系统的响应速度.为了克服上诉缺点,可以同时安装两个小尺寸的伺服阀,一个用于稳定状态,另一个用于瞬时响应,即小流量时一个伺服阀工作,当流量增大到一定值时,两个伺服阀同时工作。这样既可以满足控制系统的瞬时响应特性,又可以不受伺服阀饱和度的影响。为了减小管路对系统性能的影响,伺服阀安装时可采用直接安装的形式.1    270m二冷带轧机,根据实际的控制要求,采用单伺服阀,直接安装方式即可。伺服阀额定流量为6311min,这样即经济,又可满足控制的需要。(    2)液压缸的选择对系统动态特性的影响这里主要讨论液压缸的行程选择问题。液压缸的工作行程主要取决于液压缸的类型,对于液压缸工第4章系统动态性能分析及改进作行程小于25mm的,通常认为是短行程液压缸;工作行程大于250mm的,认为是大行程液压缸。轧制时的所需实际工作行程有轧机的大小决定。图4-2是液压缸I作行程分别为2.5mm, 25mm, 250mm的阶跃响应曲线(1, 2,3分别对应工作行程2.5mm, 25mm, 250mm)a!/今/,2/二=一灭几,,一分产,'3 -,了了犷,以,.一        -—」十0                      0.05      0.1      0.15     0.2     t/s图4-                    2液压缸阶跃响应曲线对比图Fi            g.4-2 Comparison diagram of step response of sevrovalve其中伺服阀采用的是阻尼为0.    01, Ps=21 Mpa,单行程,直接安装方式。从图中可看出,250mm的液压缸响应曲线出现频率约98Hz的波动。由于轧机的惯性,摩擦力,传感器、伺服阀响应时间和液压缸油的压缩过程会在系统中形成死区,对大行程液压缸的影响要远远大于短行程液压缸。此外行程太大的液压缸固有频率与轧机的阻尼振荡频率相接近,快速轧制时容易产生谐振,不过当轧机在低速运行时大行程液压缸的频率响应性能却很好。1270mm冷带轧机工作行程小于25mm,因此选择一个短行程的液压缸即可。(    3)压力变送器、位置变送器的选择对系统响应时间的影响当伺服阀和执行机构选定以后,只有通过正确选择反馈机构,即通过变送器的选择来减小系统的响应时间。而且液压板厚控制系统的分辨率、精度直接与变送装置的响应特性有关。所选变送器的分辨率、精度是由系统所能满足的最小板厚变化量决定的,即整个轧制过程执行机构的跟踪能力决定。响应时间应满足系统动态性能,即能对各种干扰和控制信号做出及时响应.当燕山大学工学硕士学位论文然随着变送器分辨率与精确性的提高,响应时间也会增大,使系统的响应速度减慢。因此选择变送器应综合考虑这两方面的因素,使其综合性能达到最佳。对于短行程的液压缸,位置变送器一般选择线性电压差动形式的(LV DT),分辨率理论上可达无限大,全行程误差不超过0.1%,时间常数可达0.5ms。本系统所选位置变送器采用了美国ANALOG DEVICES的LVDT专用集成电路,电路的长时间稳定性、抑制温漂的能力、在工业现场的抗千扰能力都很强,非常适合于冶金工厂环境条件下使用。对于长行程液压缸位置反馈元件的选择要复杂一些,还应增加一些辅助设备。目前,在分辨率、精度、和响应时间之间的最佳折衷方案是选择线性旋转变压器,将它与数字变送模块相结合。分辨率可达0.0005mm,精度可达0.002mm,响应时间0.25ms。压力变送器选择的自由度相对大大一些,一般要与测压元件性能相一致。此外还要综合考虑它分辨率、精确性及安装位置。4.1.3板厚控制系统的振动性    液压伺服系统虽然是高度非线性的,但在轧机一般使用的压下范围内基本是线性的,因此从研究线性系统的稳定性入手研究轧机系统的稳定性仍有很大意义。因为线性多自由度动态系统的运动完全可用系统各个变量对各种信号响应的综合来进行分析,因此我们从研究单自由((SDOF)度机械系统的运动入手,对多自由度机械系统的稳定性进行分析[21。单自由度机械系统的动态性能可用下述差分方程表示[43-451mx(                        t) + cx(t) + kx(t)二F(t)     (4-1)式中从一物体质量(KS);c—阻尼;            k一刚度系数(Nlmm);F(    t)一时间域内系统的输入;x(    t)一系统在时间域内响应输入产生的位移(mm).将式(4-1)进行拉普拉斯变换,可得(                      ms' +cs+k)X(s)=F(s)                 (4-2)设SDOF系统的传递函数为H(s),则有下式H(s)=X一厂1      5一」ms' +cs+k(4-3)第4章系统动态性能分析及改进将((4-3)式改写为典型振荡环节的形式,设阻尼系数为者,自然振荡频率为田。。可得下式H(s)=X-F5一夕(4-4)式中。。一、压,古=丁.,设阻尼系数为。,阻尼自然频率为。d,式丫脚                乙VKI脚(4-4)可进一步化为控制系统中零、极点的形式I          H(s)=式中P=。十J )od(系统的极点):’=Q一jWd (P的复共扼极点);X一F5一5=.一—卫x(  s一p)(s一P*)(4-5)多自由度线性振荡系统可看作有单一频率的单自由度线性系统的综    合。每一个单自由度系统都有一个独一无二的自然频率、质量、阻尼系数存在。因此多自由度系统的传递函数可表示为式中N一系统中总的变量个数;P.,  P:一第;个变量的共辘极点:,一变量序号(1,2,3... );P',,一第k个变量对第.1个输入的响应。(4-6)式完全表征了多自由度振荡系统的特性。利用留数定理,由(4-5)式可得          1A                          A'= ̄一一一2n一一一=—                          十一=H(s)(4-7)Pa=-和.(阻尼项);。‘一。n抓-v(阻尼振荡频率)。E'-     m,Hik (S)=r-I is一Pl As一P).二Pi -Pk(4-6)X-月5、}5(  :一p)(s一p')  s一p  s一P.燕山大学工学硕士学位论文A,  A'一系统的留数;通常认为是复向量。系统的对输入函数F(t)瞬时响应,可经过拉氏反变换来求得,即下式X (t)=L-'夏X(s))一F(t)[L-' {-二+:}]=F(t)(Ae-+A'e0*')s一p            s一p(4-8)口+J田d,P .=v一jWd,因此可得又因为P=                  x(t) = F(t)e' [AeI"'' + A'e-''I' )                 (4-9)依据线性系统的可列可加性,可得到多自由度系统中第k个分量对第.l个输入的响应函数:xk(t)=F, (t)艺e" [A,e'"0`+A,e-'"0t)(4-10)    从式(4-9)和式(4-10)可以看出,当输入降为零时,系统的瞬时响应总会衰减到零。原因可参看等式(4-9),在方括号里面,自然指数的幂均为虚数,在括号外自然指数的幂为实数。虚指数幕导致系统在稳定状态以阻尼振荡频率Wd振荡,可以通过欧拉公式e"' =cosy+siny来进行分析・设A=a十声,有((4-9)式可得:Ae'吻‘  +A'e-'=2a cos Codt一2,8 sin o)dt(4-11)它是一个正弦函数,角频率为。‘,振幅为了4a'+4尸,初相角为ar tan(-a /脚.再考虑实指数部分,上文己经说明v=一和,,由于任何实际系统都不可能有负的角频率,口只能是负值。因此在没有外部干扰的情况下,系统的输出总会衰减到零,衰减的速率由Q的值决定。因此使机械系统产生持续的自激振荡只可能有一种途径,即外加控制信号及控制系统设计不当所造成的,下面通过反馈系统来进行说明。图4-3是反馈系统的典型框图。图4-                      3反馈系统框图Fi                    g.4-3 Block of feedback system图4-3所示系统的传递函数为章--)『)式中G, (s)一为是执行机构本身的传递函数,同等式(4-7)G2       (s)一反馈机构的传递函数。等式(4-13)仍可化为式(4-7)的形式—X (s)二F(s)[式中部分分式的数目是由式(4-12)所具有的极点数目所决定的,它与反馈函数G2 (S)有关。因为机械系统反馈函数的形式由反馈机构各自的特殊性所决定,U2的特征也不,丁能相同。在系统中,a2有可能是一个正值,系统有位于s平面右半部的极点,整个系统变的不稳定。而在等式(4-4)中所有的极点都位于:平面的左半部,系统是无条件稳定的。液压板厚控制系统的振动问题可概括为以下三种[    [45-50].    (I)有外部输入F(ja))引起的,它于系统的传递函数这种情况,系统的振荡幅度主要由F(r)的幅值和系统的静态阻尼决定。在轧机上表现为各种外部干扰信号,如轧辊偏心信号,轧制时速度的波动,带材塑性系数和厚度的变化。(    2)轧制频率与轧机的固有频率相接近如轧制力的变化频率与轧机的的阻尼振荡频率co"相接近,可引起轧机的谐振,振动幅度由执行机构本身的传递函数决定。轧机的自然频率一般在60到100Hz之间,操作频率正常情况下应小于20Hz.    (3)用于信号反馈的机械执行机构和电气结构所造成的如用于液压缸位移检测的位置变送器,压力传感器及模拟变送电路的选择或设计的不合理。    在冷带轧机轧制过程中,这三种因素都可能存在,因此在考虑消除轧制时的振动性问题时,必须综合考虑三者的影响。4.24.2.1板厚控制系统性能的改进方法增强型PI  D自动位置控制算法在实时控制系统中,最常用的控制器是PI    D控制器。1270mm冷带轧机厚控系统的内环亦采用了模拟PID控制器。常规PID控制系统原理框图第一X一尸系统动态性能分析及改进G,(s)H(s)=1+G2(s)q(S)(4-12)B                     B"+—s一(Q2+jm2)  s一(Q2一jm2)十…](4-13)47                                         燕山大学工学硕士学位论文如图3-13所示。系统由模拟PID控制器和被控对象组成[["]. PID控制器是一种线性控制器,控制原理图如图4-4所示。比例分仁俪不乌u t衬   }}}}一 c(t)微分图4-                    4模拟PID控制系统原理框图Fi              g.4-4 Block diagram of simulate PID control principle从图中可看出,它根据给定值r(t)与实际输出值c(t)构成控制偏差e(                       t) = r(t)一c(t)                  (4-14)将偏差的比例(P)、积分(I)和微分((D)通过线性组合构成控制量,对被控对象进行控制,故称PID控制器。其控制规律为二,八_。r.....l e(t)dt.TDde(t),.t,,=n n1匕v)宁一-r-I(4-15)界                  dt或写成传递函数形式G(s)二U(s)一、尸(1 +于十TD s)二乙Ls)                 I'ST, TD S'十KpT,s+1Ts      (4-16)式中Kp一比例系数T,      一积分时间常数      几一微分时f ul常数简单说来,PID控制器各校正环节的作用如下:(    1)比例环节即时成比例地反映控制系统的偏差信号e(t),偏差一旦产生,控制器立即产生控制作用,以减少偏差。(    2)积分环节主要用于消除静差,提高系统的无差度。积分作用的强弱取决于积分时间常数T, , T,.越大,积分作用越弱,反之则越强。(    3)微分环节能反映偏差信号的变化趋势(变化速率),并能在偏差信第4章系统动态性能分析及改进号值变得太大之前,在系统中引入一个有效的早期修正信号,从而加快系统的动作速度,减小调节时间。为了使轧制的带材,更加符合实际的要求,在板厚控制系统的外环,    通常应加有偏心补偿,油膜补偿等控制算法。这就要求位置闭环具有较好的动态特性。1270mm冷带轧机位置闭环采用的是电液伺服控制系统,它的动、静态特性对板厚控制起着决定性的作用。AGC系统功能的进一步增强,主要取决于电液伺服系统的响应特性II1,研制具有快速响应和高精度的电液伺服控制系统是板厚控制的关键技术之一。作者在恒通1270mm可逆冷轧机上对电液伺服系统采用了增强型PID位置控制算法,使电液伺服系统具有更好的动、静态特性I2o1在3.    4.2中已得常规性PID控制模型,在此基础上本文提出增强型PID位置控制算法。。为了便于比较,现将二者共同画在图4-5中。从图4-5可rU-兰孙祠、:I  K4T,S+I}1 日      Kz    _T12Sl!】IJKf(a)常规型PID控制模型给定{二万-吕二召 一_ K4于石iT,S+I曰Ts -门一州K2门T;S(                          b)增强型PID控制模型图4-                5常规PID控制与增强型PID位置控制模型(        a) Model of general PID control (b) Model of enhanced PID controlFi          g.4-5 Model of general PID control and enhanced PID control以看出,常规型PID控制和增强型PID位置控制模型的主要区别在于PID49                                                 燕山大学工学硕士学位论文控制器作用的位置不同。常规PID控制的作用是进行液压缸位置偏差的校正,其输出是伺服阀控制回路的参考信号,而增强型PID位置控制的作用是进行伺服阀阀芯的位置偏差校正。图4-5(a), (b)结构图中外环部分是相同的,比较二者的动、静态特性,仅需讨论虚线框内部分的传递函数即可。设PID控制器的积分时间常数为乙,微分时间常数为几,比例系数为朴。,伺服阀控制部分增益合并为双K=KK4lT3),可得增强型PID位置控制算法内环传递函数Y        (s)一(T,TDSZ+K,T,s+1)K_U(s)不s2 (Ts + 1) + KK f(T,Tos2 +T,.凡s+1)奇(TTDS2・KpT,S・”KTK, T, f・S'‘命+TTD)S2 +TKpS+1.    、____‘_,_‘,T__、__T._,一一,_一二(4-17)                ’凡K‘一’‘KKf中凡一般较大,因此该条件很容易满足.当PID控制器只含比例作用(T; =TD=0)时,可重新推得传递函数KfY(S)二KfU(S)一T,一   S2+-二一一S+1共:,+2}. S+1KKP      Kf    KKPKf叭叭由式(4-17)可知,当(七共丁十不几)凡>二分-时,系统即可稳定。实际系统二一,,二1                                                      1(4-18)“中142)}一W T P}K :f}, oj}“内环”效成二阶“荡环节”的固“““”率;    乙洲.KKPf可知,当阀芯位置反馈系数Kf增大时,口。增大,氛变小,即可通过改变伺服阀阀芯位置反馈系数来改变W}, }.的大小,从而改变系统的动态性能。当T.几不为零时,亦可得到同样的结果.常规型PI    D控制情况下的内环传递函数为U(                s)不s        s(不s+l)+KK,氛一喜下三I一           I }n为内环等效成二阶振荡环节时的阻尼系数。Y          (s)=  TTD兰土T,TDs+I、一一.一K—二第4章系统动态性能分析及改进街(T,TDSZ +T,TDS“,“‘奇“’+ 1 IS+1)(4-19)可知这种情况下内环闭环传递函数含有一积分环节,而外环也含有一个积分环节,使整个位置闭环系统不容易稳定,只有通过细致调节PID控制器参数才能使系统稳定。因此,增强型PID控制加强了伺服阀阀芯位置反馈为的调节作用,相对减弱了PID控制器的调节作用,通过调节为和PID控制器参数,可使系统的动态特性、稳定性达到要求;常规型控制则完全一阵于一0           0.02     0.04     0.06     0.08     il.1  WS图4                        -6两种PID控制效果比较Fi                g.4-6 Comparison diagram two kind of PID control依靠PID控制器参数的调节,才能满足稳定性要求,这种情况下PID控制器参数调整范围较小,鲁棒性较差。图4-6为两种PI17控制情况下,阶跃0.025mm输入时,位置闭环(外环)输出Sx的响应曲线。由图可知,位置输出达稳态值60%所需时间,常规型约为0.033s,增强型控制约为0.022s.因此,增强型PID控制具有良好的动态响应特性Isi-saI4.2.2缸位移微分负反馈    为了进一步研究液压控制系统的动态性能,对系续的各项性能进行分析是非常重要的。由3.4.2节图3-14和图4-5可得如留4-7单机架冷轧机控制系统的一个检测模型。    图4-7所示的模型,可用于板厚控制系统各项性能的分析和检测。下面就通过这个模型来比较系统加有缸位移微分负反馈与不加微分负反馈时燕山大学工学硕士学位论文外部干扰调节器给定115带材K,󰀀I待轧带材礴L出带材cKu伺服阀管路液压缸I T's IK,位置/轧制力轧制力变送I   1+Ts#LL &#d1L位置变送1+T微分反馈图4-                    7液压AGC系统动态检测模型Fi            g.4-7 Dynamic examine model of the hydraulic AGC system系统对阶跃和正弦信号的响应能力。(    1)系统对阶跃信号的响应能力图4-8是缸位移变送环节加有微分负反馈进行补偿与不加微分负反馈时,系统的阶跃响应对比曲线,图中,1‘尼150.05户一、从-一:l'2尹    如      『呀,,’-,0.040.03卜--一t-----------.----------咨-0.02卜-尸盯』-.520.01_:_00.      05           0.1t/s图4-8阶跃响应曲线对比图Fig.4-8 Comparison diagram of step response curve为增加缸位移微分负反馈时的阶跃响应曲线,2为原系统的阶跃响应曲线。虽然不加补偿时的响应曲线己有较合理的上升时间和超调,仍然可以通过增加缸位移微分负反馈提高其性能。通过曲线I和2的对比可知,增加反馈补偿后,系统响应速度加快,超调量减少,系统很快趋于稳定。当通过第4章系统动态性能分析及改进改变系统的增益来改善没有补偿系统的上升时间与超调时,将受到伺服阀饱和特性的,而使系统变得不稳定。当增加缸位移微分负反馈后,可以将系统增益值大幅度提高,从而改变系统的动态响应性能。通过将输出的微分信号反馈到系统输入端,可以增大系统阻尼,改善系统动态性能。本系统采用的补偿方法是缸位移微分负反馈。因轧机APC系统可以简化为二阶环节,下面以二阶环节为例来进行定量的分析[(471。如图4-9所示R(s) E(s)必门川        -其-}一      C(s)下        万I s(s+[SW,) I    I I图4-          9微分负反馈控制的二阶系统Fig.4-9 Negative differential feedback control of second-order由图4-9,系统的开环传递函数为G(s)=山,    12} +K, ro󰀀 s[s /(2枷。+K,心)+I]K= ̄一一竺1一(4-20)式中,开环增益为2咨+K,    m}(4-21)相应的闭环传递函数为二,、jo}玉Wks)=种下,,万丁一一一,一..s-+‘百,        田󰀀s +田。式中(4-22)参=者十合Kp.(4-23)由式〔    4-20)^-(4-23)可见,微分负反馈会降低系统的开环增益,从而加大系统在斜坡输入时的稳态误差;但不影响系统的自然频率,并可增大系统的阻尼比。在设计测速反馈控制系统时,可以适当增大原系统的开环增益,以弥补稳态误差的损失,同时适当选择测速反馈系数K,,使阻尼比参在0.4到0.8之间,从而满足给定的各项性能指标。微分负反馈控制对系统输入端噪声有滤波作用,对系统组成元件没有过高的质量要求,使用场53                                      燕山大学工学硕士学位论文合广泛。微分负反馈可以改善系统动态性能,但会增大稳态误差。为了减小稳态误差,必须加大原系统的开环增益,而使K,单纯用来增大系统阻尼。(    2)系统对正弦信号的适应能力实际轧制时待轧带材厚度的波动、材料塑性系数的变化以及输入的偏心补偿信号,可通过系统对正弦信号的适应能力来进行模拟。此外由于轧机的自身因素也会对轧制造成干扰,这些干扰频率一般都大于IOHz。因此,轧制时的频率响应分析,即系统跟随正弦信号的能力,比起阶跃响应分析来说更重要。下面仍通过前面提到的测速反馈来改进系统的响应能力。图4-10是缸位移变送环节加有微分负反馈进行补偿与不加微分负反馈时的波特比较图。    图4-10中曲线a为不加微分负反馈时的系统的幅频、相频特性,b为加有微分负反馈时的系统的幅频、相频特性。由图可见,缸位移微分负反馈改善了系统的动态性能,减小了系统的相位滞后,降低了轧制成品的厚度波动。    因此说合理的设计液压板厚控制系统,能够充分发挥液压系统的潜能,增强系统的鲁棒性,改善其响应性能.4.3本章小结本章从论述轧机的动、静态负载特性入手,首先分析了液压板厚控制(电)划竖(灿弓)划毕频率(                             rad/sec)图4-1                            0伯德对比图Fi                      g.4-10 Comparison of bode diagram第4章系统动态性能分析及改进系统中典型元件的选择对轧机性能的影响和轧机的振动问题。然后采用增强型PID自动位置控制算法,加强了伺服阀阀芯位置反馈环节的调节作用,相对减弱了PID控制器的调节作用,很容易使系统的动态特性、稳定性达到要求;应用缸位移微分负反馈进一步有效地改善了轧机的动、静态特性,提高了系统抗高频干扰的能力,减小了系统的相位滞后,满足了生产的实际要求。燕山大学工学硕士学位论文第5章偏心补偿控制5.1轧辊偏心的起因从广义上讲,由于轧辊和轧辊轴承形状的不规则而造成的轧辊偏心会    导致轧件厚度周期性变化。这些不规则形状包括以下几个方面,如图5-1所示[551270=娜一1800    (a)偏心(b)双周期必2701800                                                             1800270=多90o(                c)三周期(d)实际形状1                -辊身转动中心2一轴颈转动中心3-辊身和轴颈转动中心图5一1轧辊横断面类型                          (    a) Eccentricity  (b) Double periods  (c) Three periods  (d) Practicality shapeFi                          g.5-1 Transecttype of roll    (1)相对于辊身或轴承与轴承之间,支撑辊轴承偏心。。(    2)工作辊和支撑辊身最大直径D和最小直径d之间的差值而引起的轧辊椭圆度。    (3)相对于轧辊轴承和辊身,支撑辊轴承旋转套筒的偏心。(    4)辊身圆周的不均匀性。(    5)轴承圆周的不均匀性。(    6)支撑辊轴承旋转套筒的不均匀性。第5章偏心补偿控制    (7)轧辊轴承中轴颈的不均匀性。由于这些不规则变化而引起的周期性辊缝变化通常称为轧辊偏心。    其起因可归纳为以下五类:    (t)设计不当这主要是指轴承设计性能,例如油膜轴承中的键会造成支撑辊轴承旋转套筒圆周上的局部不均匀。(    2)修磨不当修磨设备的缺陷及操作和维护方法不当而造成修磨误差。(    3)装配不当当拆卸和装配轴承以及在轧辊和轴承座之间安装轴承时,常会产生划伤和异物介入等造成的误差.    (4)轧辊和轴承变形当载荷超过轧辊和轴承的设计能力后,轧辊和轴承就会发生塑性变形和金属破坏。镶套式支撑辊过载后,支撑辊发生几何变形就会引起辊套相对心轴产生移动,此外,轧辊和轴承圆周上发生不均匀磨损也会造成轧辊偏心。(    5)辊缝控制不当根据辊缝控制类型的不同,轧辊偏心的影响可能被抑制、放大或不发生变化。5.2轧辊偏心对轧制力和板厚的影响    理论分析和实验结果表明,轧辊的偏心反映在辊缝、轧制力和带钢厚度上,是一种高频周期波。偏心信号的频率取决于轧制速度,幅值取决于轧辊偏心量及其安装的相对位置。轧辊偏心的存在导致直接辊缝的周期性变化,从而造成带钢厚度的波动.在配有常规压力AGC系统的轧机上,轧辊偏心的存在还会导致压力AGC系统误动作,使调节质量恶化[56-58]5.2.1轧辊偏心对轧制力的影响    当控制系统保持缸位移位置恒定时,辊缝将随着轧辊偏心发生周期性的变化。辊缝的变化可表示为轧辊总偏心量的函数。旋转的轧辊其偏心将引起轧制力周期性的变化[59]图5-    2说明零辊缝时轧制力的变化,即辊缝中没有轧件时机架受力情况。显然轧制力周期非正弦变化,它的主频与支撑辊旋转频率一致[59]。轧制过程中,轧制力波动图的形状多少要受到轧件厚度、硬度以及带刚张力变化的影响。但主频周期取决于支撑辊旋转频率。由于上下支撑辊存在直径差,因此轧辊偏心会出现偏摆现象。如图5-3所示燕山大学工学硕士学位论文‘、久90魂J八1]75 ,两支撑辊中辊径较小的轧辊旋转一圈比另一支撑辊旋状一圈所要的时    间要少,这个时间可表示为偏摆周期T b (60-651T6=了xVmm)式中D,, DZ -分别为较大和较小辊直径(轧机速度(MM/S)频率几为式中厂一上支撑辊旋转频率(Hz)儿一下支撑辊旋转频率(Hz)          零辊缝位置闭环条件下,通过记录轧辊旋转过程中轧制力波动很容易发现偏摆现象。在一个偏摆周期中,轧制力从最小值升到最大值,然后又降到最小值。如图(5-3)所示,当上下支撑辊相互到达偏心相加的角位置时,轧制力波动最大;而当到达偏心相减的位置时,轧制力波动最小。的60丹0          I        2        3        4        5       Nsa              =主频曲线b一实际轧制力曲线图5-        2零辊缝时轧辊偏心引起的轧制力波动Fig.5-2 Rolling Force lfuctuate due to roll eccentricity in zero gap‘飞、气V几一支撑辊旋状一周所需时间Tb                    一偏摆周期图5-                          3轧机偏摆现象示意图Fi              g.5-3 Sketch map of partial sway phenomenon of mill万D,DZ2;    r}。2一。,}(5-1).                   f6 =I!1 - f 2l                 (5-2)第5章偏心补偿控制5.2.2轧辊偏心对轧件厚度的影响    在测厚仪式的反馈厚度自动控制系统中,由于测厚仪一般安装在距离引起厚度波动的辊缝较远的地方,使整个系统存在一定时间的延迟,所以实际轧出厚度波动不能在反馈信号中得到及时的反映「66-68]。为了解决这一问题,通常采用厚度自动控制系统((GM-AGC)。即在轧制过程中,检测出轧制力P和空载辊缝So,根据弹跳方程h二SO +习K:计算出成品厚度ha在厚度自动控制系统(    (GM-AGC)中,轧辊偏心是通过对轧制力和辊缝的改变来影响产品厚度的。在一个机架中,每个轧辊的偏心量对合成轧辊偏心量。都有影响。因此,辊缝会随着轧辊旋转而发生变化。在有轧制力的情况下,由于轧辊偏心的存在而使辊缝减小时,将产生附加轧制力AP ,这个力将使轧机发生弹性形变AS,AS,=AP 1K.轧件厚度将减少Ah}(5-3)Ah,=APIKM因为轧辊偏心量e等于:(5-4)e=留.+Ah.(5-5)有式(5-3)至(5-5)可以推导出轧辊偏心。和相应的轧制力AP之间的关系其表达式如下:AP=eKMKsAKM + KS,(5-6)如2.8节所述,轧机有效刚度KS,a等于:(5-7)Ks,=Ks I(1一a)从式((5-3). (5-6). (5-7),可以算出轧件厚度变化量从与合成轧辊偏心量。的比率或轧辊偏心在轧件上的表现系数Ahe=e   1+(K,(1一a)lKs)1_(5-8)燕山大学工学硕士学位论文式(    5-3)至(5-8)中Ks. KM所代表物理量的意义与前面章节相同。设液压机构刚度设为KH,可得偏心作用的物理模型,见图5-4e图5一              41w心作用物理模型Fi    g.5-4 Physical model of eccentricity effect车L辊偏心、轧制力、和轧件厚度在P-h图中的关系见图5-5。从图5-5中可以看出,当有偏心e存在时,实际板厚减少了4h.,但由于同时轧制图5-              5轧辊偏心、轧制力和轧件厚度在P-h图中关系Fi    g.5-5 Relation among of roll eccentricity, rolling force and gauge in P-h diagram力增大,压力AGC控制系统此时会认为板厚增大了Ah2,因而在控制器的作用下系统朝着使板厚减少的方向动作。结果使常规控制器的质量进一步恶化。尤其是当实现轧机刚度可调时,随着轧机刚度的提高轧辊偏心量将完全反映在轧件上。而当轧机刚度减小时,尽管偏心影响被消除,但入口第5章偏心补偿控制厚度波动的影响仍很严重。依据2.    8节所述。下面就变刚度控制中三种典型情况下偏心对板厚的影响进行论述。(    1)恒压力控制在这种模式下Ks,, = 0,轧制力保持不变。这时,轧机的等效刚度模数为零,纵向厚差最大,在轧制过程中调整压下使轧制力恒定,故又称为等压控制。一旦轧辊偏心引起辊缝趋于闭合和张开,控制系统立即发出辊缝补偿信号以维持恒定的轧制力。这样辊缝中工作辊圆周表面之间的物理距离就可保持不变。只要有充足的响应时间,这种控制模式就能有效地消除偏心对轧件厚度的影响,即式(5-8)中的O, h= 0.(    2)自然刚度控制此种模式a = 0 )  KS,, = Ks, 4h =0PIKS。在这种情况下,辊缝执行机构的位置应保持恒定,即在轧制过程中不进行压下调整,空载辊缝So恒定。轧机合成偏心的一部分会表现在轧机上。轧辊偏心表现在轧机上的部分大小取决于轧件塑性系数KM和轧机自然刚度K:的比值。随着这个比值的增加,表现在轧件上的偏心会相应减小。    (3)恒辊缝控制此种模式a=1, Ks,=00・这时,轧机的等效纵向刚度模数为无穷大,空载辊缝调整量AS,=一△P/KS。由于偏心会造成辊缝的变化,变化幅度与偏心幅度大小相同,方向相反。系统误认为是轧件厚度的变化而进行调节,最终将偏心的变化最终全部表现在轧件上,即有, hO= e,这样,在轧机有效刚度为无穷大的控制模式中轧辊偏心量会全部表现在轧机上.1    270mm轧机在轧制的前几个道次采用位置闭环控制,主要用来消除厚度偏差,在后一两个道次采用压力闭环控制,这时的入口厚度偏差较小,以压力闭环控制来消除轧辊偏心等的影响,保证具有良好的板形。在位置闭环情况下,轧辊偏心对出口带材的厚度影响较大,对轧制高精度带材来讲,必需设法消除轧辊偏心对带材厚度的影响,以此提高板厚控制精度[4)一般情况下,轧辊偏心引起的带材厚差超出了带材厚度公差范围.例    如当两轧辊偏心最大值为50 fan,机架刚度系数为200 t l mm,材料的塑性刚度为500 t / mm,偏心使被轧带材产生14 fan以上的厚度偏差,而厚度公差一般为1 ^-5 fan,因此,若不消除偏心,便无法满足带材精度的要求。以上分析表明,轧辊的偏心对带刚厚度的影响可归纳为两个方面:    (    1)辊缝的周期性变化造成带刚厚度的周期性变化。燕山大学工学硕士学位论文                      (2)使常规压力厚度控制系统的调节质量恶化.5.3快速傅立叶变换法在偏心补偿中的应用轧辊偏心信号为一复杂的高频波,其变化周期与主传动转速成正比,    其中起主导作用的基波份量近似为正弦波[68-701。通常,偏心信号混杂在各种扰动和随机噪声之中。偏心补偿的关键就是设法获取偏心信号或其主要成份(通常是其基薄份量),来进行控制[[71-721。在常规压力AGC系统中,轧辊的偏心会引起带材出口厚度的波动,而补偿厚度波动一般采用动态改变轧制力或改变压下的方法来实现。因为电动压下机构响应慢,很难即时对轧辊进行偏心补偿,因此大部分现代轧机均采用液压压下来进行调节。对轧辊偏心问题的研究,国外始于三十年代,但大都处于仿真阶段,    成功实例很少,我国近几年才开始这方面的研究。偏心补偿较成功的技术是从轧制力传感器输出信号中提取偏心份量进行补偿,如各种数字滤波器偏心补偿装置以及傅立叶变换补偿装置等。本课题应用轧辊偏心量e和轧制力变化量峨之间的关系,首先,将轧    机预压靠运行,通过测量轧制力的变化而算得轧辊偏心的变化量,在轧机正常运行时将偏心变化量离散取反后加到位置内环之中,以达到偏心补偿的目的。鉴于快速傅立叶变换在数字信号处理方面的成功应用,本文首先系统地归纳了偏心信号的特点,指出了它在轧制过程中对辊缝、轧制力和带钢厚度造成的不良影响,在此基础上提出了1270mm可逆冷带轧机利用快速傅立叶变换法进行轧辊偏心补偿的实施方案.针对轧机运行过程偏心信号的抽样、截断,在控制中过程中偏心补偿信号的给定方法及轧机外围辅助元件的选取这些在实际应用中的技术问题进行了深入的研究,并通过仿真验证了软件算法的正确性及控制方案的可行性。依据控制方案选取了相应的外围辅助元件,并进行了大量的现场实验,根据现场中出现的实际问题对控制方案进行了必要的改进.经过理论分析,实验室仿真和大量的现场实验,最后将偏心补偿控制顺利地运用到了实际生产之中,使板带材轧制精度得到了很大提高.5.3.1傅立叶级数正弦余弦函数在区间(t    o , t。十T)内构成了一个完备的正交函数集第5章偏心补偿控制。依据正交函数性质,在这个正交区间内,用这个完备正交函数集的线性组合可以来表示任意函数f (t),如下式[73-74]f (t)=艺(a󰀀 cosnmot+b󰀀 sin n coot)=ao+艺(a󰀀 cosncoot +b; sin mot)(5-9)这就是任意函数在有限区间(to , to +T)内三角函数形式傅立叶级数展开式。当f(t)是以T=2习WO为周期的周期函数,则式(5-9)在整个区间(-00,。)内都成立。各分量的系数可分别由式(5-12)求出二二一丈+Tf (t) cos nwotdtfod‘+r          -(5-10)nrootdt6n=式(5-10)和(5-11)中令to=-T12,则有一。;a。一争lr}z f (t) dt・,一景TfZf(t) cosnmot、一争r/zf(t)s‘二%tdt逼近的均方误差最小。在傅立叶级数中,当n取有限值时,则正交函数集是不完备的,式(5-9)仅为At)的近似逼近,可以证明,由式(5-12)-(5-14)所确定的份量系数,其让N     = oo,式((5-9)还可合并成另一种形式[[74_75,            f (t)=Ao +A, cos(wot + }P,) + Az. cos(2root + 4P2) +…                                + A_ cos(ncoot +t)󰀀) +=Ao+艺A󰀀m cos(nmot)不难得出上下式(5-15)和式(5-9)有下列关系内.么一丙f、f (t) sin no),t+Tsin z nmotdt(5-11)(5-12)(5-13)(5-14)(5-15)燕山大学工学硕士学位论文Ao =a,(5-16)(5-17)(5-18)(5-19)A󰀀,󰀀=了a'.+时a󰀀=A󰀀󰀀, Cos汽b,=-A󰀀󰀀, sin ip󰀀Jg9'k=—一b"      口行          (5-20)式(5-25)所示的无穷三角亦称为傅立叶级数(傅氏级数),第一项A。称为周期函数f(t)的恒定份量(或直流份量);而第二项A,󰀀, cos(root十Pi)称为一次谐波(或基波份量),其周期或频率与原周期函数的相同;其它各项统称为高次谐波,即2次,3次,4次,:谐波。这种把一个周期函数展开或分解为具有一系列谐波的傅立叶级数称为谐波分析。指数函数集{    e'""0`) (n=0, 11, 12,…)在区间(to , to +T)内也是完备正交函数集.T=2。。为基频份量周期・任意函数f(    t)在区间(to,  t。十T)内,可以用{(e'""lO}这个完备的正交集来表示:            f (t)=Fo+F,e'-o +Fze'z-v+…+F e''woj+…+F_                    ,e-'"0‘十几z e-}2"0`十…十及,e一”肠‘+一Y F e'""0`(5-21)式中e'""0`是实变量t的复函数,上式中各份量系数可按式(5-13)计算,仍令,则有t。二一习2聪一[::le凡=当f(t)为周期信号时,式(5-22)在整个区间(一。,。)内成立。式(5-21)是指数形式的傅立叶级数。傅立叶级数的三角形式与指数形式在本质上是同一的,下面将从傅立叶级数的三角形式推导出其指数形式,从而说明二者的关系,为应用快速傅立叶变换做好准备。由式(5-9)可得厂一鱿f (t)("'"0'0` )' dte'--'(e'""0` )' dtf (t)e-'""I'dt=TrI"0'0`e-J"tdt几f(t)e-'""^`dt‘,一,2)第5章偏心补偿控制f (t)=as+艺(a" cos ncoot+b" sin ncopt)=・。・鑫【二(e ,"d},+2     e-'" O,r卜)+"   b・(     e'"a` 2-j=e-'"va )]・。・E}󰀀s, 0" -2一+'b"E!.+2.jb"一 ’”‘,定义厅= -n,则n=一厅,代入上式,得声a_一i          b__‘、r ,_a_十jb_・、.病、,了(t)一“。十丢卜飞升)e'一’+充(一寸书)e".._".        (-')-14)蕊二丁              ‘有式((5-13)和式((5-14)不难得出a-" = an, b-" = -b,故r,   ,a .一jb.)e,""0`=At)一。。十全(琴典)e'"0V十乙气一一代不一一.(a" - jb"・。・󰀀_,卢一   2}' -,’    2(5-25)若令F=则有a,一jb󰀀2    对所有n¥。:且Fo=asf(t)=艺F e'"0'0`式的关系。按c。的定义,可以得出(5-26)上式即傅立叶变换的指数形式。从中可知傅立叶变换的三角形式与指数形;一争f f (t)(cos一‘一,sinnwo t)=其中Fo一粤rf(t)dt一asT              、    将非周期信号看成周期为无穷大的周期信号,应用指数形式傅立叶级数,可以得到傅立叶正变换和傅立叶反变换,它们构成一个傅立叶变换对,其典型表示式如下1                 r f(tle-'""0`                 (5-27)T-                                0一、燕山大学工学硕士学位论文正变换F(co)一皿f (t)e-'0' dt反变换f(t)=(5-28)F(co)e'"dco(5-29)傅立叶变换对记为f(                          t) E--), F(w)    (这里信号f t)一般为连续时间信号,而F(co)称为f(t)的频谱密度函数,或简称At)的频谱,t为时间变量,w为角频率。    为了表示一个周期函数分解为傅氏级数后包含那些频率份量和各份量所占的“比重”,用长度与各次谐波振幅大小相对应的线段,按频率的高低把它们依次排列起来,就得到图5-6所示的图形。这种图形称为f(t)的频0}                     wa   3w,  5wu 7wo。图5一幅度频谱                          Fi                        g.5-6 Amplutide pedigree谱(图).这种频谱,只表示出各谐波份量的振幅,所以称为幅度频谱。由于各次谐波的角频率是。。的整数倍,所以这种频谱是离散的,有时又称为线性频谱.通过图形我们可以清楚地看到各次谐波在信号中所起作用的大小,在偏心信号中基波分量起主导作用,如图5-2所示。5.3.2用离散傅立叶变换变换对连续信号进行谱分析离散傅立叶变换(    Discrete Fourie Transform)简称DFT,目前己经成为在计算机上研究傅立叶变换的专有名词。我们从傅立叶积分变换的数值积分,可以推导出DFT. DFT分析的对象是离散信号,并且是有限长序列。变换‘肠                                     第5章偏心补偿控制得到的结果也是一个有限长离散序列,是频谱的离散值。DFT在数字技术中是一个重要的概念,研究它的性质及快速算法在数字信号处理中起着很重要的作用(741给定有限长离散信号x(    n)0<n<_N一1月>_N可定义x(n)的离散傅立叶变换X(                          k)(5-30)洞艺间x(n)e-J( N ).n.k式中k一自然数(0, 1, 2, -""N-1)    --令W,v=e'J(N ),则X (k)=DFT[x(n)]=Ix:x(n        )jV, ,k(5-31)式中k一自然数((0, 1, 2, """N-1)n一自然数(    0, 1, 2, """N-1)式中X(    k)是x(n)的离散傅立叶系数,也就是我们需要的“频谱”。如果对离散傅立叶系数x(k)作离散傅立叶反变换,可获得原序列x(n) ox(n)二,N 1 K}X.0(k)WN0k(5-32)式((5-31)和式((5-32)形成离散变化对,记作:x(n) F-), X(k) o工程实际中,经常遇到连续信号x(    t),其频谱函数X(动亦是连续函数。为了利用DFT对x(t)进行频谱分析,应先对x(t)进行时域采样,得到x(n) = x(nT),又因为计算机的存储量有限,必须用门函数进行截断。离散傅立叶变换,可以看成是作了截断与抽样处理的连续傅立叶变换的样本,下面通过傅立叶变换的数值积分计算来详细说明如何通过DFT来对连续信号进行谱分析。当用数值积分公式来计算傅立叶变换,实际上是用定积分来计算积分    区间为无限的广义积分,在这无疑是引起误差的,甚至会引起很大误差。这种积分区间由无限变为有限相当于将信号x(t)作了截断,这样作是不可燕山大学工学硕士学位论文避免的。不仅在数值积分中X(t)要做截断,就是在离散傅立叶变换中,x(t)也要作截断,这是计算傅立叶变换的共同问题。用定积分来计算广义积分引起的误差是可以计算的,可设函数x(t)是因果信号,即x(t)=o, t<o,则X (m)=f x(t)e--dt=rx    (t)e-'"dt截断所引起的误差为(5-33)“〔。)二〔x(t)。一‘dt在允许误差内。(5-34)对于实际信号,这种误差是可以在允许范围之内,这是因为〔96),(1)当To足够大时,x(t)可以是十分小,因此󰀀n (m)也十分小,有可能(2)有许多非周期信号是时限信号,即x(t)=0,  t>To这时傅立叶正变换实际上就是定积分运算,误差Ax()o) = 0 a傅立叶正变换的数值积分公式可用式(5-33)来推导:AT.艺x(nAT》一,ATX(ro)一rX(t)。一*-(5-35)其中To为截断长度;N为样本数;采样周期AT                       =Tp/N                    (5-36)显然此时X(叫仍是G)的连续函数,如图5-7(c)所示.对X(oc)在区间[[0, coo I上等间隔采样,采样点数仍为N点,采样间隔为F(频率分辨率)。如图5-7(e)所示。参数wo, To, N和F的关系如下式F=4?o=-止一N         N-AT(5-37)由于N-AT =T6。所以(5-38)                       F=2兀将to=kF和式(5-37)代入(5-35)得一,2N'!mX (kF)=AT.艺x(nT)e加璐第5章偏心补偿控制令X (k)二X(kF),x(n)=x(nT),则X(k)一△T.艺x(n)ex(n)=x(nT)一1-2}=4T     . Dh!'[x(n)]____(5-39)juznk,FN[ 1梦X(k沁]=x(X(、k*,)e' LN '0k_一N     =ko(5-40)4TIDFT[X(k)]通过式((5-39)可知,连续信号的频谱特性可以通过对连续信号采样并进行DFT再乘以△T的近似方法得到。因此说,离散傅立叶变换可以看成是作了截断与抽样处理的连续傅立    叶变换的样本,分析截断和抽样对频谱的影响,可以了解对连续信号采样式应注意的问题,即如何决定采样间隔时间和怎样对其进行截断。现用广义函数6(t)形成的冲击序列,对连续信号x(t)作理想抽样。这种理想抽样信号仍属于连续信号,可对它作连续傅立叶变换,可用图解分析其频谱变化,截断和抽样的影响是有以下三个方面引起的。    (1)时域抽样引起频谱周期变化变换对X(t) H X(oc)如图5-7(a)所示・对x(t)作理想抽样,即将x(t)乘以单位冲激序列-5,11(t) I SIT (t)是冲激序列,其频谱8.0 (w)亦是冲击序列。则有下式x, (t)=x(t) - 6"T (t)=x(t)艺6(t -nAT)=艺f (nt1T)S(t一nAT)9- (W)=。。Y .Y(。一noco )(5-41)(5-42)式中。0=2习AT变换对soT (t) H凡。(。)如图5-7(b)所示根据卷积定理,时域两函数相乘的傅立叶变换是这两个函数频谱的卷积即m)x(t) -凡T (t) H上X(*凡o(})=X. (})2厂(5-43)这是x(    动和周期冲激频谱凡。(。)作卷积,其结果是将频谱x(动以周期。二二2到OT延伸为周期变化的频谱。式(5-43)的傅立叶变换对表示于图燕山大学工学硕士学位论文傅氏变换对<片=令(  a)c二二;>(  b)普<二=合(  c)心二=令(  d)心==)(  e)图5-      7连续信号进行的谱分析Fig.5-7 Spectrum analysis of continuous signal5-7(c),它说明时域抽样对于频谱的影响。设x(t)的截止频率为口.,则理想时域作离散抽样必须满足。,=2川AT                 ? 2o).=m0(5-44)否则,由于抽样频率不满足抽样定理会产生“混迭效应”,抽样后的频谱第5章偏心补偿控制上X(动* 8_.Jro)就和原来频谱X(动不完全一样。产生了对X佃)的第一次27r、”’修改。    (2)时域截断引起频谱的变化时域抽样所得到的无限多个样本值是不适于计算机作计算的,可以乘以门函数STo (t)进行截断,得到一个有限长时间函数X2 (t)=X(t) "56T (t) - STo (t)U(t + To)一U(,一T.),其频谱为G:。(口)门函数g,-O(t)二2’        2-sin些STO (t) -- GTO (W)TOS.( 2O)(5-45)伽)为抽样函数,如图5-7(d)所示,时域乘以门函数,对于频域为和(口)作卷积X(m) *凡o (ro)l * GTO (W)=X2(rw)  (5-46)X(t) -SAT(t)'8TO(t)。。牛1牛这一傅立叶变换对如图5-7(e)所示,经过和GTO (W)卷积后,频谱出现了皱波,这是对原有频谱X(co)的第二次修改。显然,TO越大,皱波越小。(    3)频域抽样引起时域信号周期变化为了得到离散频谱,应对图5-7(e)所示频谱作理想抽样,乘以冲激序列8,* (w),它对应在时域中亦为冲击序列STO (t ),设STo (t)=艺S(t - nTo)则Se. (oc)=△。Ys(。一kAr))式中△。二2TrITO根据时域卷积定理,频谱乘以戈. (动对应于时域作卷积z (t)=NO .15.T (t) .9 TO (01 *乓o (t)频率抽样亦必须满足抽样定理,即(5-48)(5-47)          (5-49)t)的样本值不完否则,由时域卷积会产生混迭效应,z (t)的样本值就会和X(24w ->几燕山大学工学硕士学位论文全一样,会再一次在时域产生修改。    从图5-7可以看出截断和抽样对傅立叶变换的影响。以上分析表明,截断和抽样对于原来的傅立叶变换对都会产生某种程度的修改。截断和抽样如果处理不当,则会引起畸变,离散的样本值则不能代表原来的信号和频谱。5.3.3傅立叶积分变换的离散变换对    下面,我们从采样定理出发来讨论如何从傅立叶变换得到它样本值的离散变换对1941:    (1)截断时域截断为有限长度To的有限信号。频谱也截断为有限长,当作带限信号。其截止频率为aim,带宽(    2)样本数频域抽样间隔为△co,时域抽样间隔为△T,则时域样本数频域样本数wo=2叽                         (5-50)N,                       =To /AT                   (5-51)从=wo                       /OW                    (5-52)频域及时域离散化,其变量关系为t=n4T                                                         oc=kAro                    (5-53)(    3)采样个数的确定设时域及频域抽样都各自满足抽样定理,时域作N:点抽样,则有To = N,AT;时域作NZ点抽样,则有co。二Ni 4co。有时域抽样定理m:二2习OT _>2 m󰀀=。。,可取。。=2AT,即NZOW = 2可OT,它要求                      TO-)c4=2可Nz                   (5-54)同理从频域抽样定理可得,To _<2 习Ow,可取To= 2△。,即N, AT = 2习△。,要求Oc                       o-AT=2习N,                   (5-55)若取N, =N2 =N,则有AT                         - Aco=2到N(5-56)由式(5-35),令co =kAm可得X (kAco)=OT工x(nAT)e-jATev n k第5章偏心补偿控制.........曰..曰...且口X(k)二L -T 0 NS'-'f (n)eN    _0(5-57)式中n—自然数((0, 1, 2, """N--1)k—自然数(    0, 1, 2, """N--1)k, N取值同前,有式(5-40),可得下式山0        xkn)=二甲了I                            V I‘夕「              k- N2tx(X(kk,j)ej _2N ,rnk(5-58)式(5-57)和式(5-58)表示了傅立叶变换的离散变换对,记为x(n) <--+ X(k) e图5-7所示傅立叶变换的离散变换对如图5-8所示。在轧机偏心补偿控制中N值取为1024,可有效保证采样的精度。图5-8离散傅立叶的离散变换对Fig.5-8 Discrete transform couple of Discrete-Time Fouriertransform5.3.4离散傅立叶变换和傅立叶级数的关系有5.    3.3节可知,一个理想的抽样信号[751,它的傅立叶变换同样是一个周期冲激序列,它们各自的样本值在一个周期内组成了DFT变换对。分析对任何意信号f (t)都是适用的,因此任何一个周期序列都可以通过傅立叶分析建立它的离散谱,也就是通常所说定义的离散傅立叶级数(DiscreteFourier Series,简称DFS) o周期为N的周期序列,可以表示为    s(                         n)=z(n+rN)                    (5-59)这种周期序列,可以分解为一系列正弦一余弦序列或复指数序列,但它与连续傅立叶级数不同,周期序列的分解,仅包含N个分量。从图5-6和图5-7我们可以清楚的看出这一点。这是因为指数份量。j(等)-n(k+rN)是周期为N燕山大学工学硕士学位论文的周期函数,而且仅包含N个不同的数值j (} )-n(k,rN)e‘。/(N)n一。(N)-k(n+,N)J(N )}'kj因此(5-60)只需N个份量就可以把序列表示出来Y        (n) = I翼X(k)e(          5-61)在展开式中,n,k意义同前,各份量的周期为叼k,角频率为k.2习Na在k=1时为基波,其周期为N,基波角频率为22习N,。,,IN > e・,(N)为基频序列,它是一个随,作周期变化的序列.X(k)是周期为叼k份量的系数。式中还乘以1(N,它可以使DFS和DFT的变换形式统一。经过推导可得式(5-61)中的系数洲艺间.z(k)x(n)e一,(NM-k(5-62)由上式可知,X(k)是一个周期为N的周期序列,和DFT一样可令风=。WN -e-J ( N)则DFS变换对可以写为洲艺词    --双n)W Nmkk”的谐波由于x(的是周期性离散信号,离散的特点决定了不存在周期N<1k                      而周期性的特点决定了k=0, 1, """N-1,变化时只存在有限的N次谐波而不能象连续信号那样有无限多个谐波份量,因为当k>N时,并不出现更高次的谐波而是又重复出现了k=0, 1,,二N-1次谐波,这是由时域及频域双重的周期性决定的,它使x(n)及X(k)都是N项级数取和而构成n个采样点的序列。实际上可知离散傅立叶变换正是在一个周期内的离散傅立叶X-x(5-63)_上nS -1X(k)W -k"N完:a      --(5-64)第5章偏心补偿控制级数。离散傅立叶变换是有周期性特点的离散变换对。这个结论可通过以下推导得到。设有一有限长序列x(    n)・、・,一}“昙of     x- t、  _!x(n)0<_n5N一1n   ?N9所示。它们之间的关系如下可将x(n)延拓成为一周期序列x(n),如图5-r-N O! N-1 N      2N- I图5-    9有限长序列的周期延拓Fi  g.5-9 Period ductibility of finite serialx(n)一艺x(n+rN)OSn<-N一1(5-65)      (5-66)定义周期序列x(的从。二0到n=N-1第一个周期的范围叫做“主值区间”其它n值N一2(也可以规定为而周期序列在主值区间的值x(n)称为“主值序列”。即x(n)是x(n)的主值序列,而Y(n)是x(n)的周期延拓。因此可有这样的结论:DFT是DFS的主值序列,在主值区间的DFS就是DFT,因此可以说当DFT的频率分辨率与DFS的基波频率相等时,它们之间的系数关系是一一对应的。既有DFT各次谐波的的幅值、幅角可得对应DFS各次谐波的的幅值、幅角。因为偏心信号周期性的特点,用DFT及其逆变 一<     一燕山大学工学硕士学位论文换即可得到我们所要偏心信号。5.3.5快速傅立叶变换法的优点    直接计算DFT,计算工作量相当大,特别当数据点数N很大时,使用计算机也颇费时间。DFT的计算是大量的复数乘法及加法。由DFT的定义式及其矩阵矩形可知,DFT的运算需要将输入的N点数据x(n),与复系数卿‘依次两两相乘再求和。每计算出一个X(k)值,需进行N次复数乘法和N-1次复数加法。一共有N个X(k)需要计算,因此计算X(k)需要作NxN=N2次复数乘法和Nx(N-1)次复数加法。而实际数据处理中,一般N值可以相当大。例如N =2"  =1024。全部DFT的工作量需要N2二220 =1048576次复数乘法,Nx (N一1) =1047552次复数加法。都是上百万次复数运算。而且还需要用计算机计算出全部复数变换系数WNki又是N2个复数加法。此外,还需要将输入数据和中间运算数据存储起来,一般所需存储容量与运算次数成正比,当N很大时,用计算机直接计算DFT就要耗费极长的计算时间和占据整个计算机的存储容量,这不仅无法实时处理数据,也很不经济[741由于这个原因,1    965年美国的Cooley和Tukey为了计算DFT而提出的FFT发表以后,引起了世界各国广泛的注意,20年来发展极快。目前不仅有各种FFT的计算程序,而且有各种FFT专用硬件及FFT处理机作为商品出售。FFT之所以受到如此重视,首先式FFT将DFT的计算速度提高了    叼1092 N倍。使许多信号的处理工作能与整个系统的运行速度协调,它的应用就从某些数据的事后处理及系统的模拟研究,而进入到数据的实时处理,另外在提出FFT之前,在通讯、雷达及其它领域,数字信号处理在速度、成本方面都赶不上模拟系统,FFT的出现使得用数字系统分析频谱,其优越性超过了模拟系统,为广泛应用数字方法处理信号打开了崭新的局面。    FFT算法的基本原理在于利用DFT原始变换矩阵的多余性,化简了原始变换矩阵。把计算长度为N序列的DFT应用蝶形计算结构逐次地分解成计算较短序列的DFT。现在常用的有两种基本的FFT,第一种称作时域分组(或叫时间抽选),它将x(n)(标号n常表示和时间有关的量)逐次分解成较第5章偏心补偿控制短序列作蝶形流程计算。第二种称作频域分组(或称频率抽选)是将X(k) = DFT[x(n))分解成较短序列(标号k常表示和频率有关的量)作蝶形流程计算。只有对DFT有充分的理解,才能理解FFT的运算结果,因为FFT只不过是DFT的一种特殊算法,是一种在计算机上DFT的快速计算法。要讨论FFT算法,必须首先了解DFT变换系数Wn    N的周期性和对称k性。下面将分别进行讨论。5.3.6离散傅立叶变换系数的周期性和对称性j u n.k是一个复指数周期序列,它对于序数n和k都具有周期性WNmk、-特点,其周期为N。即存在着衅‘=衅+              r-N)-k=衅(k+rN)=Wink+rN)            (5-67)虽然共有n-k个系数,但只有N个系数是的,其余系数只不过是这N个系数的重复值。x(    n)的DFT可以看成x(n)的Z变换在单位圆上的均匀抽样,而其变换系数WnNk也可以看成是单位圆上均分这一个单位圆的复单位矢量。对于N=8,其分布如图5-10所示i771。从图可以看出,衅‘不仅具有周期性,图5    -10矢量WNk的分布Fig.5-10 Distributing of vector of衅也具有对称性。在单位圆上:W"0,,WN二一1。于是有(        5-68)Wmk.W万=-WNmkWN (              nk+z)也即,凡是相位差180’的单位矢量它们一正一负互相对称这样对于N=8”  燕山大学工学硕士学位论文的八个系数,存在着                  W8a = -Wea        WSb = -W¥zWs                   e      = -W, B          W=一        g = -WBs这样一来,64个系数由于周期性只有8个是不同的,其余的值都和这8个相同,再根据对称性分析,这8个系数实际上只有四个系数不同,所以由这64个系数所组成的系数矩阵中,存在着大量系数雷同。这就造成DFT计算中,存在着不必要的重复计算。5.3.7快速傅立叶变换及相关程序信号流图可以形象地表示出FFT的整个运算过程,因此,它就成为分    析和表示FFT的重要方法。在运用信号流图作FFT分析时,是将长序列的DFT运算分解成较短序列的DFT运算,并在运算中利用变换系数衅k的对称性和周期性,产生蝶形运算结构,然后用信号流图将它表示成蝶形计算流程图,这种蝶形计算流程成为一种FFT的基本运算形式。可以利用这种蝶形流程图将较长序列的DFT运算分级进行蝶形运算,从而建立整个蝶形流程的FFT信号流程。有了这种FFT信号流图便可以编写计算程序,从而使用计算机作FFT信号流程[73]    应用蝶形流程表示FFT运算时有两种基本形式,在5.3.5节己详细论述,由于频率抽选与时间抽选计算原理相同,本课题主要针对时间抽选的方法进行研究。以下针对第一种方法进行详细分析。首先选取N为2的整数幂,即    N=2M, M为整数。将X (k) =艺x(n)W"‘的计算过程中,首先将x(n)分成两个州2点序列:一个序列由x(n)的偶数点组成(n=0, 2, 4"二),另一个序列由x(n)的奇数点组成(n=l, 3, 5"二)。这样可得到(n)WN k+艺x(n)衅‘X(k)=艺x(n) -衅‘=艺x月为侧盈月为奇盆(5-69)为了便于分析,应作变量置换:n为偶数,令n=2r; n为奇数,令。=2r+lo可得下式n,、    ,r,2,、公2r + 1FN2r+1)-k=A tK)=乙x(Gr厂、+乙x(r_0                                r=0第5章偏心补偿控制式中-I(N)2W,2=。=e  2=叽,于是上式可写成71)WNKt二(2;十lwN‘一。(*卜W󰀀 H (k)(5-X(k)=乞x(2r)W}k、,二。          万丝_        1式中G(k,一鳌二(2r)W,",为原始序。二(。。。。、偶数点、N12 A DFT'.0万‘称为原始序列x(n)中。为奇数点时叼2点DFTH(k)=艺x(2r +1)W,'-0            2这两个变换式是符合DFT定义的。我们可以先求得叼2点的DFT, G(k)及H(k),然后利用分组关系式X(k) = G(k)十叮H(k)来求得所要计算的X (k)。以下以N=8为例,画出它的时域分组信号流图,如图5-11所示。其    中G(k)及H(k)的DFT以框图表示,利用分组关系式(5-91)表示出X(k)的计算流程.    由图5-11可以看到一个N点的DFT己被分解成两个叼2点的DFT,但是,我们看到G(k)和H(k)只有叼2个点,k=0, 1.…(叼2一1),而X(k)却是N个点。从关系式X (k) = G(k)十叫H(k)中只能计算出k=0,  1.…(叼2一1)的X(k)值,如果要从G(k)即H(k)计算出全部X(k)的值,需要利用G(k), H(k)具有周期性的特点。它们是叼2点的DFT,是以叼2为周期的周期序列,有如下关系式:_对于分组关系中的力口权系数(WNk+Z当k>舒艺间N-1                  x(2r)(W,') ''t+W N乞二(2r + 1)(W,2 )'-k一,(N)(5-70).=0                                    N_,G(k +叼2)=G(k)H(k +叼2)二H(k)  时N工Z(5-72)燕山大学工学硕士学位论文二(k"2)_二*二夸_二*kW(NWNkWN-W,分组关系式73)X(k+ 2)=G(k)-W󰀀H(k)形流程图,如图5-12所示。G(o)(5-74)依照这个过程在进行第二次,第三次。经过化简就可得到一个8点FFT蝶XM(p)十叫Xm (4)(5-75)Xm(p)一WNX. (9)如果选其中p" 9相邻两级蝶形运算元素的序号,m表示蝶形运算得级数。取N=2"',则经过M次时间抽选N点DFT都变成了各级的对偶运算。在这个完整的蝶形图上可以看到每一级蝶形运算都是一组运算,    每一蝶形对偶结的运算不用于计算其它结点,而它本身的运算也和其它运算无关,从而加快了运算速度。傅立叶逆变换原理上与快速傅立叶变换相同,这里不过多讨论。    由给定N=8的蝶形计算流图5-12可以看到,输入数据x(n)的标号编排不是自然顺序,而是以x(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)的次序进入计算机存储单元。这是我们为了缩短序列长度N而将n按奇偶分组所造成。也就是FFT运算中存在着标号编排问题。这在实际安排程序时必须注意,但却是不难解决的,只要把十进制表示为二进制数码以后倒置,码位置倒置的x(0)x(2)x(4)x(6)x(l邓x(5x(7X(o)X(1)X(2)X(3)X(4)X(5)X(6)X(7)图5    -11  DFT的时域分组流图Fig.5-11 Flow chartof  time domain of DFT第5章偏心补偿控制x(n)序列顺序就会按奇偶分组。这样做的好处是可按即位运算来安排流程,即运算时保证同一对对偶蝶形运算的结果仍存放于原来这对数据的地址中,不需要另外安排存储器来存储各级中间运算数据。在FFT运算中,一般以安排作即位运算的方式居多,因为这样做可以节省内存。X(o)XMX(2)X(3)X(4)X(5)X(6)X(7)H(                         3)图5-                    12八点DFT快速蝶形算法流图Fi        g.5-10 Flow diagram of fast butterlfy arithmetic of eight point DFT    下面给出快速傅立叶变换及其逆变换的C语言子程序[771,程序共分码位置倒置子程序,加权系数计算子程序,蝶形运算子程序等三部分。其中计算加权系数时采用的是实部与虚部分别计算的方法。这样编程条理清晰,利于实现,可轻松与主程序中其它部分进行联接。void kkfft(pr,pi,n,k,fr,fi,l,il)i  nt n,k,l,il;doubl  e pr[],pi[],fr[],fi[];  {int it,m,is,ij,nv,nu;doubl  e p,poddr,poddi;/*采样值由主函数给出中/f    or (it=0; it<=n-1; it++)        {m=it; is=0;f        or (i=0; i<=k-1; i十+)/*码位置倒置混序子程序*/燕山大学工学硕士学位论文      {j=m/2; is=2*is+(m-2*j); m--j;)f      r[it]=pr[is]; fi[it]=pi[is];}      pr[0]=1.0; pi[0]=0.0;/*计算加权系数*/p=6.283185306/(1.0*n);pr[1]=cos(p); pi[1]=-sin(p);               /*I=0,进行傅立叶变换*/if (1!=0) pi[l]=-pi[l];                    /*1=1,进行逆傅立叶变换*/for (i=2; i<=n/2-1; i++)  { pr[i]=pr[i-1]*pr[l]-pi[i-1]*pi[l];pi    [i]=pi[i-1]*pr[l]+pr[i-1]*pi[l];}  /*计算加权系数的实部*//*计算加权系数的虚部中//*蝶形运算*/for (nu=1; nu<=k; nu++){m=pow(    2,(k-nu)); nv=pow(2,nu);    for (it=0; it<=(m-1)*nv; it=it+nv)f        or (j=0; j<=(nv/2)一1; j++)          {poddr=pr[m*j]*fr[it+j+nv/2]-pi[m*j]*fi[it+j+nv/2];poddi            =pr[m*j]*fi[it+j+nv/2]+pi[m*j]*fr[it+j+nv/2];f          r[it+j+nv/2]=fr[it+j]-poddr;/*傅立叶变换的实部*/f          i[it+j+nv/2]=fi[it+j]-poddi;/*傅立叶变换的虚部*/f          r[it+j]=fr[it+j]+poddr;/*傅立叶变换的实部巾/f          i[it+j]=fi[it+j]+poddi;/*傅立叶变换的虚部巾/}          }    if (1!=0)f  or (i=0; i<=n-1; i++)    毛fr[i]=fr[i]/(1.0*n);/*逆傅立叶变换的实部*/f      i[i]=fi[i]/(1.0*n);                 /*逆傅立叶变换的虚部*/}      if (il!=0)f  or (i=0; i<=n-1; i++)*计算傅立叶变换的模及幅角*//{     pr[i]=sgrt(fr[i]*fr[i]+fi[i]*fi[i]);/*计算傅立叶变换的模+/第5章偏心补偿控制if (fabs价[i])<0.000001 *fabs(fi[i]))/*计算傅立叶变换的幅角巾/  {if ((fi[i]*fr[i])>0) pi[i]=90.0;el    se pi[i]=-90.0;}  elsepi  [i]=atan(fi[i]/fr[i])*360.0/6.283185306;}    return;5.4实际生产中应注意的问题    在轧辊偏心补偿中正确使用快速傅立叶变换,在轧机上还必须安装一些辅助设备。1270mm冷带轧机在轧辊上安装有与支持辊同轴的光电脉冲发生器,使用光电脉冲发生器的两组脉冲信号,一组是零点定位脉冲,即相应于轧辊圆周上某一预先固定的点,每转发出一个脉冲,此脉冲作为检测、补偿和控制的初始定位信号;另一组是圆周定点脉冲,即沿轧辊圆周每转发出100个脉冲,以保证采样和控制的准确定点。在轧辊偏心补偿控制中利用快速傅立叶变换时,数据的采集工作很重    要。传统的采样方式一般固定抽样间隔时间,但轧辊的转速在实际中很难固定为恒定值,这种抽样方式便失去了意义。针对这个问题,1270mm板厚控制系统采用等间距抽样来实现,如图5-13所示,将某一信号f (t)等间隔图5-                            13采样分析图Fi                        g.5-13 Analysis diagram of sampling分成N等份,然后提取各点处的值。不管函数的周期增大还是缩小,只要是N不变,各抽样点得到的值就不会改变。在轧辊偏心信号的提取中,即可利用光电脉冲发生器的第二组脉冲信号,在一定时间内对轧辊偏心值等吕3                                       燕山大学工学硕士学位论文间距采样。无论轧辊是否匀速转动,对采样值都不会产生影响(751因为要利用快速傅立叶变换,因此抽样点数应当为2    N个,其中N取整数值.本课题采用的N值为10,即抽样值共为1024个,这样大的抽样点数很容易满足其谱分析范围。抽样点数不足时补零,大于1024时可删去多余数据。截断时间应等于或稍大于轧辊的偏摆周期。此时经傅立叶变换得到的谱分辨率正好等于其基频频率值,便于对信号进行分析。通过下位机,将计算得到的参数模型转换成与第二组脉冲序列相对应    的离散点模型。偏心控制环节将前述环节处理出来的偏心信息转换成为调节液压推上机构的电压信号,以补偿轧辊偏心的影响,动作频率约为5mso在偏心控制器中运用了两级中断,中断A与零点定位脉冲相对应,中断B与圆周定点脉冲相对应,中断A的优先级高于中断B的优先级,两级中断的运用不仅确保了控制的精确性,而且有效地避免了偏心控制器的误动作.中断B的作用使得本系统不仅可以在正常轧制状态下补偿偏心影响,而且能够在加减速阶段对偏心干扰加以抑制.1    270mm冷带轧机采用上述控制方案进行偏心补偿,取得了很好的效果。实际运行时轧制力,出口厚度变化如图5-14.之9087858280(a)轧制力变化情况鹰0.3100.3000.305偏心补偿投入0.2950.2900  24681012141618                  (b)板厚变化情况图5-                14偏心补偿效果图(a) Things of rolling press vary ( b) Things of strip gauge varyFi    g.5-14 Effect diagram of eccentricity compensation94                                       第5章偏心补偿控制5.5本章小结    本章首先详细分析了轧辊偏心的的特点,掌握了其性质,指出了它在轧制过程中对辊缝、轧制力和带钢厚度造成的不良影响。然后进一步分析了轧辊偏心对轧制力和对轧件厚度的影响。在参阅大量文献的基础上,提出了利用快速傅立叶变换来对轧辊偏心信号进行处理,从而提取其主要成份(基频份量)的方法,重点对其在实际中应用中的一些关键问题提出了解决方法。经过偏心补偿后的控制系统在实际应用中取得了很好的效果。燕山大学工学硕士学位论文结论    自2001年9月份起,协助导师完成了唐山恒通精密薄板有限公司1270mm四辊可逆冷带轧机液压微调厚控制系统的研制、安装和调试工作。系统于2002年三月份正式投入运行,在轧制速度和精度等方面完全达到了控制要求,取得了巨大的经济效益。每月最高可轧制带钢9000吨,年轧制量可达10万吨以上,板厚控制精度达到了微米级,完全达到了设计能力。为今后研究更高精度的同类轧机提供了可靠的实践检验和设计依据,对推动轧机生产有积极的意义。本文以实际课题为背景,取得了以下研究成果。(    1)9入剖析了轧辊偏心对板带材厚度的影响机理,提出了1270mm可逆冷带轧机用快速傅立叶变换法进行轧辊偏心补偿的实施方案。很好地解决了偏心信号的抽样、截断,偏心补偿信号在正常轧制时的给定方法,轧机相应外围辅助软件的选取这些在实际应用中的技术问题.最后将以快速傅立叶为基础的偏心补偿控制成功地用到了实际生产中,使板带材轧制精度得到了很大提高。(    2)利用增强型PID自动位置控制算法,加强了伺服阀阀芯位置反馈环节的调节作用,相对减弱了PID控制器的调节作用,很好地克服了常规型PID控制方法系统鲁棒性差的缺点,有效地降低了现场调试和操作的难度。为了进一步增强液压板厚控制系统的响应性能,应用缸位移微分负反馈,改善了系统的动态性能,减小了系统的相位滞后,满足了实际生产的要求。    (3)用机理分析的方法建立了液压板厚控制系统的动态数学模型,并对数学模型进行了必要的简化。仿真结果表明,用机理方法建立的数学模型可以精确评估系统的各项性能。需要说明的是,在冷带轧机液压板厚控制系统的研究过程中,偏心控    制问题的研究尚需深入发展,还有一些问题急待解决。主要包括以下二方面:(    1)本课题只对轧辊偏心分量中的基波分量进行了补偿,如何解决偏心信号中高次谐波成份的补偿问题。(    2)如何将一些较先进的理论,如神经网络滤波,小波变换理论应用到轧辊偏心补偿控制领域上来,其中重点要解决其实时控制问题。参考文献参考文献1李虎兴,葛愁琦,陈贻宏.冷轧窄带钢生产.北京;冶金工业出版社,1995:66-702  Tezuka. T, Yamashita. T, Application of a New Automatic Gauge Control System for theTandem     Cold Mill. Conference Record一IAS Annual Meeting,2001,2:955-9603傅作宝.冷轧薄钢生产.北京:冶金工业出版社,1996:1-104  V. Ya. Osadchij, A. L. Vorontsov. Influence of Eccentricity on Tube Bending duringDr    awing. Intermet Inzhiniring,2001,(4):48-495  H.E. Merritt, Hydraulic Control Systems. New York:Wiley,1967:238-2416  Kondo. Katsuya, Kubo. Toshimasa, Estimation of Rolling Parameters by SignalPr    ocessing of Measured Data. Tetsu-To-Hagane/Journal of the Iron and Steel Institute ofJ    apan,2002,88(6):314-3197王一译冷轧机新型控制系统国外电气自动化1991,(2):11-198  Donald J.Faplano. Control of Strip Thickness in Hot Rolling. Iron and SteelEngi    neer, 1985,62(11):34一39  Hajime Watanaba,Chief. Control Equations for Dynamic Characteristics of Cold RollingTa    ndem Mills. Iron and Steel Engineer,] 974,51(2):59-6410 Nakatsuji Takeshi, Mori Atsunobu. The Tribological Effect of Mechanically ProducedMi    cro-dents by aM ico rDiamond Pyramid on Medium Carbon Steel Surfaces inRol    ling-sliding Contact. Kobe City College of Technology, Department of MechanicalEngi    neering,2002,36(6):663-674I1 Kimura, Hiroshi. Recent Development of Aluminium Rolling Technology. SumitomoKei    kinzoku Giho/Sumitomo Light Metal Technical Report, 1989,30(3):42-4912 Shibata Masamichi, Goto Masao, Ohta Atsuhiko. Development of Long Life rollingBea    rings  For  Use  in  the  Extreme  Conditions.  ASTM  Special  TechnicalPubl    ication,2002,(5): 8-1013 Chen Weiguo, Chen Zhangwei, Ding Fan. Analysis and Modification of Hydraulic RollerBendi    ng Force Control System.. Kang 1'ieh/Iron and Steel (Peking),2002,36(11): 37-4114魏衡江.轧钢过程自动化技术的新进展(下).冶金自动化.北京,1996,20(3):1-415魏衡江,轧钢过程自动化技术的新进展(上).冶金自动化北京,1996,20(2):5-1016 Wang Wenjing, Liu Zhiming, Li Qiang. Anti jamming Technique Research in Dynamic已7                                          燕山大学工学硕士学位论文St    ress Test Under Random Loading. Proceedings of SPIE一The International Society forOpt    ical Engineering,2001,(10): 15-1717彭天乾.微米级冷轧带刚厚度控制系统.冶金自动化北京,1996,20(6):1-518连家创,刘宏民板厚板形控制J匕京:兵器工业出版社,1996:20-2519金学俊.液压AGC在板带轧机上的应用.液压气动与密封,2000,24(8):45-4720朱义国,毛志忠,刘赤兵.带材轧制中轧辊偏心控制问题的综述.基础自动化,    1998,(3):1-421 Remn-Min Guo. Analysis of Dynamic Behaviors of Tandem Cold Mills UsingGe    neralized Dynamic  and  Control  Equation.  IEEE  Transactions  on  IndustryAppl    ications,2000,36(3):842-85322 Kazuma  Gumi.  New    Gage  Control    System   for  Tandem.  Iron  and   SteelEngi    neer, 1994,71(12):42一623 Yuan-Liang Hsu, Chang-Pin Liang, Song-Jau Tsai. An Improvement of HAGC Responsef    or CSC No. I HSM. IEEE Transaction on IndustyAr pplications,2000,36(3):854-86024张殿华,王军,褚恩辉.东北大学国家重点实验室冷连轧机厚控系统.冶金自动化.北京,    1999,23(1):24-2725陈建华,张殿华.压力AGC有关问题的综述.钢铁研究,1999,(9):36-3926王军,王国栋.厚度计型变刚度控制研究.轧钢,2001,18(6):3-527 Remn-Min Guo. Evaluation of Dynamic Characteristics of HAGC System. Iron and steelEngi    neer, 1991,68(7):52-62.28 Birol Yucel, Okumus Ertugrul, Kara Gokhan. The Response of Twin Roll Cast AluminumSt    rips to Thermo-mechanical Processing. Proceedings of Sessions, TMS AnnualMe    eting,2002,15(l 1):l 107-111329 Katori. Hideo, Yoshitani. Naoharu, Ueyama. Takatsugu, Nyu. Shuichi. Application ofTwo-    degree-of-freedom Control System to Automatic Gauge Control. Proceedings of theAmer    ican Control Conference, 1992,1(6):24-2630叶光平.液压式自动厚度控制.系统与装置.2000,24(4):30-3331 Bilkhu. T. S. Dynamic Control of Tension, Thickness and Flatness for a Tandem Cold Mill.AI    SE Steel Technology,2001,78(10):49-5432周易,徐开兴.四辊单机架冷轧厚度控制系统.冶金自动化.北京,自动化系统与装置,    1995,(2):20-24参考文献33 Hayashi Chihiro, Utakoji Masaru, Yamada Tateo. Renpei. Mandrel Mill HydraulicSc    rewdown Control. Sumitomo Metlas,1990,42(4):139-14734 Li Mouwei, Zhang Yng, aZhang Dazhi. Ways of Increasing Efficiency nd aProductPr    ecision of Cold Tandem Strip Mill at Anshan Steel. Kng aT'iehIlron ndaSt    eel,2001,36(11):32-3635 J.J. Fitzpatrick, J.G. Lenard, The Wear of Tool Steel Rolls during Hot Rolling of lowCa    rbon Steels. Iron nd aSteelmaker,2001,28(4):45-4936 Wng aZhenxinag, liu jinajun. Error Analysis for AGC Model nd aResearch for the NewF  unction of Mill Characteristic Cuvre, 7`h IFAC Symposium on Automation in MMM,1992:    356-36137 Nikiforuk,RN. The Large Signal Response of a Loaded High-Pressure HydraulicSe    vromechanism. The Institution of MechanicalEngi neers, Proceedings, 1965,180:32-3638 Zhang Disheng, Ling aQihong. The Application of Dynamic System Optimization ndaAdapt    ive Algorithm in Hot Rolling. ACTA Automatica SINICA,1992,18(5):573-57939 K. Aistleitner, A. Kugi, nd aB. Manhartsgruber. On-line Identification of DeformationRe    sistncae  and  Mill  Stretch  Coefficient  for  Cold  Rolling  Mills.  inPr    oc.METAL'98,1998,2(5):12-1440胡寿松.自动控制原理.北京:国防工业出版社,2002:90-9541 M. Bodson nd aS.C. Douglas, Adaptive Algorithm for the Rejection of PeriodicDi    sturbances with Unknown Frequency. Automatica,1997,33,(12):2213-222142 C. Canudas de Wit, P. Lischinsky. Adaptive Friction Compensation with Partially KnownDyna    mic Friction Model. Int. J. Adapt. Contr. SignalPr ocess, 1997,(11):65-8543 M. Espana, L. Prlay. On the Global Dynamics of Adaptive Systems: A study of na    ElementayEr xample. SIAM J. Contr. Optimi,1993,31(5):1143-116644 Vladmir B. Ginzburg. Selection of Optimum mill Configurations for cold rolling. ALSESt    eel Technology, 1999,(12):34-3945赵静一,王益群,王永昌四辊轧机弯辊电液伺服控制系统动态特性的研究.中国机械S程,    1998,9(4):6-1046 Monaco G. Dynamics of Rolling Mills-Mathematical Model nd aExperimental ResultAI    SE Year Book,1977:399-41047 NesslerGL, nd aCoryJ r. Cause nd aSolution of Fitfh Octave Backup Roll Chatter on 4-HCol    d Mills nd aTemper Mills. AISEYear Book, 1989:482-486g9                                          燕山大学工学硕士学位论文48 Brown D. L. Grinding Dynamics,University of Cinci-nati Doctoral Thesis, 197649 Remn-Min Gun. Evaluation of dynamic characteristics of HAGC system. Iron and steel    Engineer, 1991,68(7):52-6250 E.Teoh, G. Goodwin, W. Edward, and W. R. Davies. An Improved Thickness Contorllerf  or a Rolling Mill. in Proc. 9a' IFAC Congress, Budapest, Hungary, 1984:101-10851 Liang Qihong, Zhang Disheng, Geng Qingbo, Chen Jindi. Techniques for ReducingVa    riation Within one Plate. Journal of Beijing Institute of Technology, 1996,5(2):148-15452 Feldmann, F, Kerkmann, M. Process Optimization for a 6-high UniverslRa eversing ColdMi    ll. MPT Metallurgical Plant and Technology Intemational,2001,24(4):118-12753 Lawrence K.Saunders. Identifying Control Modernization Requirements for a Tandem    Cold Mill. Iron and steel Engineer,1982,69(9):37-4254陶永华,尹怡欣.电气自动化新技术丛书一新型PID控制及其应用.机械工业出版社,    1998:1-1055金兹伯格(Ginzburg V B.)著,姜明东,王国栋译.高精度板带材轧制理论与实践J匕京,    冶金工业出版社,2000:205-21056 Patula E J. Experimental Verification of Roll Eccentricity in Sleeve Motion. Iron and Steel    Engineer, 1975,(1):55-5857赵嘉蓉,杨节,徐木庆.消除轧辊偏心的引起的轧制力滤波器.武汉冶金科技大学学报,    1996,19(4):407-41158 K.Aistleitner, L.G. Mattersdorfer, W. Haas and A.Kugi. Neural Network for Identification    of Roll Eccentricity in Rolling Mills. J. Mater. Process. Technol,1996,60:387-39259 S.S. Garimella, K.Srinivasan. Application of Repetitive Contorl to EccentricityCompe    nsation in Rolling," ASME J. Dynam, Syst., Meas., Contr,1996,118:657-66460 Kugi. Andreas, Haas. Werner, Schlacher,. Kurt. Active Compensation of Roll Eccentricityi    n Rolling Mills. IEEE Transactions on IndustyrA pplications,2000,36(2):625-63261 Garcia. Daniel F, Lopez. Jose M, Obeso. Faustino. Multiprocessor-based Equipment forRe    al-time Diagnosis of Roll-eccentricity Influence in Finishing Hot Strip Mills. IEEEI    ntenrational Symposium on IndusrtilaEl ectronics, 1997,3(7):7-1162 Garcia. Daniel F, Lopez. Jose M, Suarez. Francisco J. Novel Rela-time Fuzzy-based    Diagnostic System of Roll Eccentricity Influence in Finishing Hot Strip Mills. IEEETr    ansactions on IndustyrA pplications, 1998,34(6):1342-134963 Choi. S.Q Johnson. M.A, Grimble. M.J. Polynomial LQG Control of Back-up-rollEc    centricity Gauge Variaitons in Cold Rolling Mills. Automatica,1994,30(6):975-992匆参考文献64 Malnov. aA.L. The Decreased Effect of Eccentricity of Back-up Rolls on the Thickness of    a Strip During Cold Rolling at Reverse Mills.乃azheloe Mashinostroenie,1994,(8):25-2865 Iwashita. Tetsuyuki,  Suzuki. Yoshihiro, Sakai. Toshio. New Roll EccentricityCompe    nsation System. Furukawa Review, 1993,(8):117-11966 Carlos Canudas de, Laurent Praly. Ad叩tive Eccentricity Compensation. IEEE Transactionon     Control System Technology,2000,8(5):757-75667 Daniel F. GarcIa, Josd M. Lbpez, Frnciasco J. Sufirez. A Novel Real-time Fuzzy-BasedDi    agnostic System of Roll Eccentricity Influence in Finishing Hot Strip Mills. IEEE    IndustyApplr ication Society Annual Meeting New Meeting, Louisiana, 1997:2117-212468祝诗平,王学勤,苏玉刚.小波消噪原理及在轧辊偏心补偿控制中的应用.农业机械学报,    2000,31(4):99-10269李育苗.改进的快速富里叶变换及其在轧钢机偏心补偿中的应用.自动化学报,1    990,16(2):151-15470刘丁,马宁,郑岗.神经网络自适应滤波在轧辊偏心控制中的应用.重型机械,    2000,(5):3-671楼雅斌.消除轧辊偏心影响的自适应数字系统的分析与研究.鞍钢技术,2000,(3):27-2972马宁,刘丁,郑岗随机逼近发法对轧辊偏心信号的检测用与控制.重型机    械,1999(6):31-3373张彦仲,沈乃汉.快速傅里叶变换及沃尔什变换J七京:航空工业出版社,1989:101-13074易克初,田斌,付强.语音信号处理.北京:国防工业出版社,,2000:150-16575陈逢时子波变换理论及其在信号处理中的应用.北京:国防工业出版社,1998:200-21076张贤达.通信信号处理.国防工业出版社北京:国防工业出版社,2000:180-19277王国栋,刘相华.金属轧制过程人工智能优化出版社.北京:冶金工业出版社,    2000:20-25燕山大学工学硕士学位论文攻读硕士学位期间承担的科研任务与主要成果在攻读硕士学位期间协助导师完成了1270mm冷带轧机液压板厚控制系统的的研制,安装与调试工作,并发表了以下学术论文:1.杨景明,曹智文,冯雅丽,车海军.1270mm冷轧机液压厚控系统中的一些改进措施.冶金自动化.    北京,己录用2.杨景明,曹智文,液压板厚自动控制系统的设计和性能分析.自动化仪表.上海,已录用    3.冯雅丽,杨景明,车海军,曹智文.1270mm冷轧机液压厚度自动控制系统.中国自动化学会第1    7届青年学术会议,北戴河,2002:434-438.致谢    本课题是在导师杨景明教授的亲切关怀和悉心指导下完成的,在论文完成之际,对导师所给予我的关怀教诲表示最为深切的谢意。导师以渊博的学识和前瞻的眼光,为学生开阔了研究视野,丰富了专    业知识。导师追求真理、敬业爱业、忘我的投入精神,使我在课题的进行过程中受益匪浅,他虚怀若谷,永远进取和谦逊无私的崇高品质是我终生学习的榜样。师恩厚重,学生将永远铭记在心。感谢冯雅丽、赵亚军及其他同学所给予的帮助和支持。    最后,衷心感谢电气工程学院自动化系全体老师多年来的辛勤培养和    教诲。

因篇幅问题不能全部显示,请点此查看更多更全内容