用于调整和控制的方法及装置

申请号 CN200680049429.3 申请日 2006-12-18 公开(公告)号 CN101346676B 公开(公告)日 2011-02-09
申请人 ABB公司; 发明人 佩尔-埃里克·莫登; 蓬图斯·贝里斯滕; 里卡德·林德奎斯特;
摘要 本 发明 涉及用于调整和控制具有变化的材料流率的工业过程的方法及装置。其中装置包括:用于将激励 信号 加到 控制器 输出信号 (U(t))上的加法器(10b);用于响应于所述 激励信号 来测量所述特性的测量系统(6);基于模型的调整单元(10c),该调整单元适于:基于所述特性的所述测量结果(z(t))和来自于控制器的输出信号(u(t)),估计描述变化的材料流率作用的过程模型结构的至少一个具有未知值的参数的值;基于所述参数的估计值,计算用于描述从控制器输出到控制器输入的动态特性的模型;基于用于描述从控制器输出到控制器输入的动态特性的所述模型,执行控制器的基于模型的调整。
权利要求

1.一种用于调整控制器(10)的方法,该控制器控制具有变化的材料流率的工业过程的特性(p(t)),
该方法包括:
-注入被加到控制器输出信号(U(t))上的激励信号
-响应于所述激励信号,接收所述特性的测量结果(z(t));
-选择过程模型结构,该模型结构描述变化的材料流率的作用,且是时间连续的,该过程模型结构包括至少一个具有未知值的、独立于变化的材料流率的参数;
-基于所述特性的所述测量结果(z(t))和来自于控制器的输出信号(u(t)),估计所述至少一个参数的值;
-基于所述参数的估计值,计算用于描述从控制器输出到控制器输入的动态特性的模型,并基于此执行控制器的基于模型的调整。
2.根据权利要求1所述的方法,其中该方法包括:
-选择过程模型结构,该过程模型结构按每个输入/输出对包括至多一个被所述过程模型结构的变化的实际的过程增益所影响的模型参数。
3.根据权利要求1所述的方法,其中该方法包括:
-选择过程模型结构,该过程模型结构按每个传动装置包括至多一个被传动装置的变化的实际响应时间影响的模型参数。
4.根据权利要求1所述的方法,其中该方法包括:
-在过程模型中包括由于数据处理和通信导致的任何时延;
-选择所述过程模型结构,使得所述时延的变化按每个模型输入不影响多于一个模型参数。
5.根据权利要求1所述的方法,其中该方法包括:
-使所述过程模型结构由两部分组成:一部分描述传输行为,其中与传输有关的滞后和与传输有关的时间常数与材料流率成反比;另一部分独立于材料流率。
6.根据权利要求1所述的方法,其中所述控制器用于控制所述材料 在轧机系统中的扎制处理,所述轧机系统包括使所述材料通过的轧辊,用于测量通过轧辊之间的辊隙的材料的平整度的测量辊,以及将经过测量辊的材料卷起的缠卷机,以及其中,该方法包括:
-为所述模型的每个部分导出一阶加滞后模型;
-引入传输行为参数(ξ),其定义与传输行为的纯一阶响应有关的纯滞后度;
-在纯滞后传输行为的情形中,使得滞后由已知距离(Ld)来描述,其中该已知距离(Ld)为从辊隙至测量辊的距离;
-在纯一阶传输行为的情形中,使得相关时间常数由已知距离(Lm)来描述,其中该已知距离(Lm)为从辊隙至缠卷机的距离;
以及
-估计过程增益(K)、独立于流的与传输有关的时间常数(T)、独立于流的与传输有关的滞后(D0)和传输行为参数(ξ)这些参数中的至少一个参数。
7.根据权利要求1所述的方法,其中该方法包括:
-用预测误差方法估计至少一个参数。
8.根据权利要求7所述的方法,其中该方法包括:
-在大的预测误差的情况下,根据从预测误差大小和所估计的先前预测误差的标准差之间的比值所计算的信度因子来减小测量结果和回归量的作用;
-从而实现保护以避免异常值。
9.根据权利要求8所述的方法,其中该方法包括:
-以基本采样速率对过程输出进行采样,应用具有相同采样速率的滤波器,重新采样到较慢的采样速率,此较慢的采样速率被用于来自过程操作的值的所述记录,在预测计算和参数估计更新中将该采样、滤波和重新采样的作用包括进来。
10.根据权利要求1所述的方法,其中该方法包括:
-执行识别实验,意味着利用影响所述工业过程的至少一个被操纵变量的有意图的激励来对过程进行操作;
-在所述识别实验期间在线估计至少一个参数。
11.根据权利要求1所述的方法,其中该方法包括:
-在线计算控制器的所述的基于模型的调整。
12.根据权利要求10所述的方法,其中该方法包括:
-实时地向用户显示结果;
-从而使得用户能够一旦获得满意的结果就决定结束识别实验。
13.根据权利要求10所述的方法,其中该方法包括:
-在过程操作期间在线估计至少一个参数,同时用所述控制器控制所述至少一个参数;
-在线执行用于控制器的所述基于模型的调整计算;
-在线将调整结果应用于所述控制器;
-从而实现真正的自适应控制。
14.根据权利要求13所述的方法,其中该方法包括:
-在通过控制器的在线更新进行操作期间估计比识别实验期间估计的参数更少的参数;
-允许来自于识别实验的参数估计结果限定至少一个参数的值,并将该参数或这些参数视为在通过控制器的在线更新进行操作期间是已知的。
15.一种用于调整控制器(10)的装置,该控制器控制具有变化的材料流率的工业过程的特性(p(t)),其特征在于:
-用于将激励信号加到控制器输出信号(U(t))上的加法器(10b);
-用于响应于所述激励信号来测量所述特性的测量系统(6);以及
-基于模型的调整单元(10c),其适用于:基于所述特性的所述测量结果(z(t))和来自于控制器的输出信号(u(t)),估计描述变化的材料流率的作用的、时间连续的过程模型结构的至少一个具有未知值的参数的值,所述参数独立于变化的材料流率,基于所述参数的估计值,计算用于描述从控制器输出到控制器输入的动态特性的模型,基于用于描述从控制器输出到控制器输入的动态特性的所述模型,执行控制器的基于模型的调整。

说明书全文

用于调整和控制的方法及装置

技术领域

[0001] 本发明涉及用于调整和控制具有至少一个传动装置的工业过程的方法及装置,以影响所生产产品的质量,其中生产速度或材料流率是变化的,从而产生了变化的控制条件。本发明优选地针对辊轧机应用中的自动调整和自适应控制,例如用于采用任何数量机械的或其他的传动装置来改进轧制产品的平整度。本发明还可应用于变化的材料是诸如板、带、织物(web)或流体的其他过程。

背景技术

[0002] 通过辊轧机的工作辊的辊隙分布(profile)和轧制带的厚度分布来确定轧制产品,如带的平整度。带平整度于是可受到对影响轧机及其工作辊隙分布的不同控制装置的操纵的影响。这种传动装置可以是机械装置,如工作辊弯辊装置、中间辊弯辊装置、倾斜装置、中间辊窜动辊装置、顶凸传动装置,或诸如工作辊冷却/加热传动装置之类的热装置。
[0003] 在金属冷轧的平整度控制中使用许多影响平整度分布的传动装置。在标准的解决方案中,平整度偏差借助于轧机矩阵而被映射到传动装置空间,其从各传动装置描述静态平整度响应。此分解产生许多控制环路,每个控制环路用于一个传动装置。这些环路被配备有PI控制器。当前这些控制器的调整是基于对每个环路模型的离线识别。此外,由于例如变化的辊轧速度所引起的模型中已知的变化在参数调度方式中被考虑。使每个环路的过程增益(process gain)不确定的两个因素是:其对轧制材料的依赖性,以及所采用的轧机矩阵和轧机实际行为之间可能的差异。
[0004] 当轧制带时,重要的是总是保持期望的平整度分布。与期望平整度的偏差可能导致损失很大的带破裂以及所生产的带卷报废。因此平整度控制系统的任务是使实际的平整度分布尽可能地接近期望的平整度分布,这对控制系统的计算速度和准确度提出了高的要求。
[0005] 无论变化的条件如何,为减小质量变化并将生产质量保持在规范以内,都需要对工业过程的控制器进行很好地调整。特别是,对高生产率的期望常常会对将过程控制得足够好以避免生产不合格材料以及避免由于破裂而导致生产中断的能提出挑战。造纸机中的纸张破损或轧机中的带断裂是导致高昂生产损失的实例。
[0006] 控制器的调整常常是基于从实验数据(例如从步骤测试)找出简化模型的例程,结合自动查找良好控制器调整的方法,并假定此模型足够好地代表了过程的行为。为使此例程成功,关键是在实验期间模型公式能够捕获过程的实际行为,且所获得的模型在对其中可能出现条件变化的过程的正常操作期间保持有效。通过适当地关注健壮性,调整方法可以允许所采取的标称行为有一些变化。如果过程的动态特性中有本质的变化-例如具有变化的生产速度-则已知的这种变化通过参数调度来处理。如果变化在调整实验期间出现,若应用标准方法则模型识别会被严重干扰。
[0007] 常用手段是黑箱识别,即估计模型的离散时间公式中的参数,在与实际控制中使用的采样周期相同的采样周期中表示。然而,在某些过程的动态特性随生产速度而改变的情形中,由于模型的真实参数值会变化,因此估计会被速度变化干扰。这包括:如果以每单位时间进行采样,此为最普通的情形,则随后动态特性中随速度变化的部分会给出变化的离散时间模型参数,如果有时是由于实用原因按每材料流量进行采样,对于时间特性反比于速度的动态特性,还给出恒定的离散时间模型参数值,则随后动态特性中不随生产速度变化的部分会给出变化的离散时间模型参数。
[0008] 对于PI控制环路的调整,在平整度控制的标准解决方案中,变化的辊轧速度阻止了用黑箱识别方法来确定模型。采样按每单位长度来执行,故传输行为的模型对采样数据来说是不变的,但对于以此方式采样的数据,传动装置动态特性不能得到不变的模型,因此整个离散时间模型会变化。此外,采样周期会因为变化的下采样倍数以及与此相关的变化的前置滤波而变化。
[0009] 为获得准确的控制,应当基于过程如何响应被操纵变量的改变而很好地调整控制器。辊轧过程的增益依赖于未清楚得知的许多参数。对于平整度控制,相关增益由轧制的材料、采取的轧机矩阵与真实矩阵的实际一致性以及其他事项所影响。
[0010] 在将初始控制任务,即使用若干传动装置保持生产材料的整个宽度上的平整度,分到若干传动装置测量环路中后,当前的控制策略是基于试运行期间的标准控制环路调整。这通常以对每个环路进行单独的离线模型参数估计和对该模型的调整来完成。为使该行为更加有效率,应当在线估计相关的模型参数并实时显示给用户,以决定何时结束该行为。而且,估计应当以不被速度变化所干扰的方式来执行。
[0011] 一个问题是,如果过程或材料改变,则控制会变得不准确,即便之前是准确的,这将导致产品质量低下或出现废料。材料改变例如会导致特性如材料的厚度、宽度或硬度的改变。在材料改变后,控制需要适应新材料的特性。控制适应新材料越快,生产质量恢复得越快。
[0012] 为避免该问题,自适应控制得以应用。自适应控制通常应用于黑箱模型。然而,因为时不变传动装置相关的动态特性,由于速度变化而具有变化的采样周期的离散时间模型会有变化的参数值。因而基于这种模型的黑箱识别的自适应控制将无效。
[0013] 与用于材料流率变化的工业过程的控制器调整有关的问题,如辊轧期间的平整度测量,是对受控特性的测量的采样速率依赖于流率,如辊轧速度,这导致采样间隔的变化。估计时间离散模型的参数是简单且已知的。然而,由于采样间隔的变化,不能使用时间离散模型来为具有变化材料流率的过程调整控制器。

发明内容

[0014] 本发明的目的是解决与调整控制器相关的问题,其中控制器用于材料流率变化的工业过程,变化的材料流率导致变化的控制条件。
[0015] 根据本发明的一个方面,此目的通过权利要求1所限定的方法实现。
[0016] 该方法包括:
[0017] -注入加到控制器输出信号上的激励信号
[0018] -响应于所述激励信号,接收过程特性的测量结果;
[0019] -选择过程模型结构,该过程模型结构包括至少一个具有未知值的参数,该模型结构描述变化的材料流率的作用;
[0020] -基于特性的测量结果和来自于控制器的输出信号来估计参数的值,[0021] -根据参数的估计值,计算用于描述从控制器输出到控制器输入的动态特性的模型,并基于此执行控制器的基于模型的调整。
[0022] 本发明使用工业过程的现有知识和已知数据。根据本发明,描述变化的材料流率的作用的模型被使用。已知模型的一些参数。剩余的那些未知的或预计会变化的参数,典型地为过程增益,通过识别来在线地确定。这与调整规则如lambda调整规则结合形成调整方法。当被应用于在线控制时,其形成了自适应控制,该自适应控制考虑过程和材料的不确定性和变化。当被应用于自动调整时,它给出迅速且可靠的调整控制器的途径。本发明允许控制任何类型的传动装置。
[0023] 通过使用连续时间模型公式并考虑有关过程的已知内容来选择模型结构,使得未知参数独立于材料流率以及独立于其他变化的参数。“独立”意味着:不管变化的材料流率或变化的其他参数如何,真实值保持合理的恒定。从而需要通过估计来确定的参数的数量被最小化,结果是该估计变得更可靠。基于模型的调整方法于是被应用于根据所估计的参数、已知参数和关系以及实际应用的任何信号处理所确定的总模型。
[0024] 本发明的方法创建了改进的、稳定的和健壮的控制系统。采用所提出的发明将大大减少与目前解决方案的材料流率变化所引起的变化的条件相关的质量控制问题,特别是与变化的材料参数和缓慢改变的过程参数相结合的质量控制问题。
[0025] 本方法基于离散时间黑箱模型识别实现优于自适应控制器的优点,通过估计而得出的参数将在更长的时间段保持恒定,例如独立于生产速度,这使得仅仅应用自适应更加可行。这进而提供与非自适应控制比较起来实际上能够提高控制性能的基础。提高性能的益处通常在于相同质量平时的更高的生产速度,和/或在相同生产速度时的更高的生产质量,这意味着废品更少或价格更优或二者皆有,和/或生产停顿更少。用作对最初调整的支持,该方法还能减少试运行次数。
[0026] 本发明可应用的实例是金属冷轧中的平整度控制。通过具有四阵列形式的传感器的测量辊来测量张力分布。测量系统每旋转一周取得四次完整的张力分布样本,并将该样本转换成平整度分布。因此存在按每材料流的量进行采样的情形,或更准确地,按轧制材料的长度进行采样的情形。使用所呈现的来自传动装置的静态响应,测量能够被映射到多个控制环路。不同控制环路中的实际动态响应主要包括:传动装置伺服器的响应,或者在热装置的传动装置情形中的加热或冷却响应,以及与带传输相关的响应。与传动装置有关的部分被视为时不变的,与传输有关的时间特性(滞后(dead-time)、时间常数)反比于辊轧速度。此类型过程的离散时间模型的黑箱识别是无用的,因为结果仅对于该操作点有效。另一方面,本发明的方法提供在整个操作范围都有效的估计,且易于被用户解释。
[0027] 此外,本发明使操作员能够完全使用自动模式,在保证同样质量的情况下,以生产废品更少和/或辊轧速度更高的形式来提高轧机的产出。
[0028] 根据本发明的实施例,本方法包括选择过程模型结构,该过程模型结构是时间连续的,且包括独立于变化的材料流率的参数。本实施例使得尽管采样间隔可变但仍能够估计用于模型的参数。
[0029] 根据本发明的另一个方面,其目的通过可直接载入处理器或计算机的内部存储器计算机程序产品来实现,该程序产品包括当程序在计算机上运行时,用于执行根据所附的一组方法权利要求的方法步骤的软件代码部分。计算机程序在计算机可读介质上提供或通过网络如因特网来提供。
[0030] 根据本发明的又一个方面,当程序用以使计算机执行根据所附的一组方法权利要求的方法步骤,且该程序在计算机上运行时,该目的通过其上记录有程序的计算机可读介质来实现。根据本发明的再一个方面,该目的通过权利要求19所限定的装置来实现。
[0031] 该装置包括:加法器,用于将激励信号加到控制器输出信号上;测量系统,用于响应于所述激励信号来测量所述特性;以及基于模型的调整单元,适于基于所述特性的所述测量结果和来自于控制器的输出信号,对描述变化的材料流率的作用的过程模型结构的至少一个具有未知值的参数的值进行估计,基于所述参数的估计值,计算用于描述从控制器输出到控制器输入的动态特性的模型,并基于所述用于描述从控制器输出到控制器输入的动态特性的模型,执行控制器的基于模型的调整。附图说明
[0032] 为更好地理解本发明,将参考以下附图。
[0033] 图1示出了具有一个轧机机座的辊轧机、控制装置、传动装置、平整度测量装置和平整度控制系统的概略图。
[0034] 图2示意性地示出了本发明。
[0035] 图3主要示出了自适应平整度控制的框图
[0036] 图4示出了自适应平整度控制系统中不同方法步骤的流程图
[0037] 图5更详细地示出了图4的方法步骤“更新所需模型参数的在线估计”。
[0038] 图6在图形中示出了作为预测误差和所估计的预测误差的标准差之比的函数的信度因子(belief factor)。

具体实施方式

[0039] 本发明涉及用于调整自动控制器的方法,该方法可用于一次性(one-shot)调整且对于自适应控制也同样良好。在本发明的典型应用中,工业过程在某种程度上具有时不变的动态特性,以及在某种程度上具有反比于生产速度或材料流率的特性,诸如滞后和时间常数的特性。可以根据与过程行为和生产速度的影响有关的实际知识来选择参数表示,使得所估计参数的真实的但未知的值将为常数,其独立于速度和其他变化的参数。从而,估计不会被速度变化所干扰,且要估计的参数更少。结果是,短时间内获得更准确的估计,且在自动调整对话期间或自适应控制的任何时间段中,根据所计划的调度,这种估计不会被运行正常生产的过程所干扰。实际的控制器调整使用结合了包括过程速度的影响的过程行为的先前知识以及结合了通过所估计的参数获得的知识的模型。
[0040] 在控制系统中执行的任何信号处理/传送/映射的作用也并入用于调整的模型中。因此,除对于已知变化的参数调度外,本发明还对需要由估计来获取的未知变化提供自适应。与先前已知的方法相比,结果是控制得更好、获得得更快并且准确性更高。这使得不合格产品更少,且在保证质量同时生产速度能够更高。
[0041] 在金属冷轧的平整度控制的实例中,已知的是,传动装置的动态特性以及任何由于通信和处理所导致的时延均独立于辊轧速度。速度可以变化,其会影响描述传送行为所需的滞后和/或时间常数。而且,测量系统可以将滤波器应用于基本测量,并可以对测量结果进行下采样,因此实际采样周期变为基本采样周期的整数倍。可根据辊轧速度而有意地改变滤波器的传递函数和下采样倍数。速度、滤波器和采样倍数的这些变化是已知的。
[0042] 对每个控制环路所要估计的参数可以选择为:
[0043] 时间常数Ta,描述传动装置响应,
[0044] 滞后D0(独立于辊轧速度的总滞后的部分)
[0045] 增益K,以及
[0046] 因子ξ,描述传输行为在何种程度上为纯滞后或时间常数。
[0047] 根据需要,这些参数真实值独立于速度和变化的信号处理。以长度来表示,纯传输滞后为从辊隙至测量辊的距离Ld;对于传输响应,纯时间常数为从辊隙至缠卷机的距离Lm。为了以时间来表示这些距离,可以简单地将这些距离除以速度。距离为常数且已知。
[0048] 由于对于不同的传动装置来说传输响应不同且不能预先得知,故假定对每个单独的传动装置来说,纯滞后的程度为ξ,纯时间常数的程度为(1-ξ)。这定义了传输行为参数ξ。独立于速度的滞后D0表示从测量到因传动装置运动而导致的辊隙变化的实际开始的总的纯时延,即,处理时间、通信、处理时延以及传动系统中任何纯时延的平均和。
[0049] 根据本发明,在试运行或重新调整期间,在自动调整会话中,当这些参数(K,Ta,D0,ξ)通过估计得以确定时,除了增益K之外,这些参数均被期望保持为常数。由于增益会随变化的材料以及所假定的轧机矩阵与真实的轧机矩阵的可能偏差而改变,因此自适应控制被应用,其中增益是用于估计的唯一参数,同时,已知参数的变化以参数调度的形式予以考虑。
[0050] 如上所述,根据本发明的自适应平整度控制按每个环路采用模型的先有知识,并通过递归识别来在线确定其增益。结合lambda-调整规则(或一些其他的调整规则),这形成用于每个环路的自适应控制,其在每个新卷材初始辊轧期间被应用。模型结构、传动装置的动态特性、依赖于速度的传输行为以及独立于速度的滞后被用作自适应控制的先有知识。
[0051] 其中所包括的参数(Ta,D0,ξ)优选地在递归识别实验中被确定,对于每个环路都被执行一次,激励信号在开环路中被提供给传动装置或在闭环路中被加到控制器的输出上。传输行为所涉及的已知因子为辊轧速度和轧机的物理距离。作为调整的先有知识来考虑的其他因子是所应用的采样速率和滤波。这在单个卷材的辊轧期间也可以改变。采样周期为用于基本测量采样的周期的倍数。此倍数可以改变,且基本采样周期反比于辊轧速度。移动平均滤波器被应用,其长度可以改变。该移动平均滤波器应用于基本采样速率。
[0052] 用于参数值的初始确定的递归识别与在自适应控制期间所用的递归识别相当类似,除其所处理的参数的数量无需估计之外。它们使用模型结构的知识以及变化的采样和滤波,并对模型的连续时间公式中的参数进行估计。从而,所估计的参数独立于实际使用的速度、采样和滤波。对于某些情形,与标准的教科书式的方法相比,递归识别还被提供有一些健壮性的增强,如异常值(outlier)保护、自适应死区(adaptation dead zone)以及对所估计的参数值的限制。
[0053] 图1-4通过示出如何将本发明应用于平整度控制而对本发明进行了阐明。图2示出了一般情况。
[0054] 如图1所公开,平整度控制系统1被集成在包括具有若干传动装置3和轧辊4的轧机机座2的系统中。开卷机(未示出)向轧机机座2馈送带5,以及使该带5通过轧机机座2,借此,该带5经过平整度测量系统6,如“Stressometer”系统,并在缠卷机7上被卷起。该轧机机座可控制轧辊4的倾斜、弯曲和/或移位,从而其改变辊隙的形状,该辊隙就是当带5通过轧辊4之间时受到作用的地方。还可以有影响辊隙的加热或冷却传动装置或任何其他类型的传动装置。该辊轧过程的最终产品是具有期望平整度的轧制带5。在串接的轧机中(未示出),进入的带可来自在前的轧机机座,而不是开卷机。轧机之间的这种变化从平整度控制这点来看并不重要,平整度控制系统可相同。如图3中所示,平整度控制系统1被设计成包括多个先进的组成单元,具有所需的全部功能性。
[0055] 如图2中所公开,主要示出了本发明的流程图,根据本发明的方法包括选择模型结构。本发明进一步的步骤用调整会话8或自适应控制9两种模式都可以执行。在这两种情形中,基于来自过程操作的数据,激励可被应用于过程且参数估计为其他情况下是未知的参数提供值。随着模型结构被给出,先前未知的参数通过估计被确定,以及利用所代入的已知参数和关系,得以对过程模型进行设计。然后,如果执行了任何信号处理,则信号处理的作用也被添加以形成完整的模型。该完整的模型被用于控制器的调整,其中可应用任何恰当的基于模型的方法。随着新的数据在新的采样时刻变得可用,模型确定和基于模型的调整被重复进行。
[0056] 在调整会话期间,参数估计和调整的结果被显示给用户,从而使得能够判断和决定何时结束会话。如果用户决定的话,则调整会话的最终结果会被应用于控制器。调整会话能够在开环情况下执行,也能够利用控制器在闭环操作中执行。在自适应控制期间,控制器工作于闭环操作,且调整结果出现时被应用于该控制器。无论控制是否会被激活,以及无论控制是否会是自适应的,都通过所编程的条件和/或操作员的决定来确定。
[0057] 在图3中,基于本发明的调整方法,公开了自适应平整度控制系统。单元控制器10可被解释为多变量控制器或多个单环控制器。此处描述的是后者,假设控制任务要被分成等于轧机传动装置的数量的多个控制环路,每个控制环路具有自己的单输入单输出控制器。
[0058] 控制器输出u(t)被馈送给传动装置3或传动装置伺服器(未示出)。传动装置的动作影响轧机机座2中的辊隙,从而带的局部特性p(t)正好通过辊隙。测量系统6通过辊隙和缠卷机7之间的测量辊11来感测带5上的张力分布s(t)。利用Stressometer系统,测量辊每旋转一周s(t)被采样四次。平整度测量系统6将原始测量结果转化为平整度分布z(t)。这可以包括滤波和重采样。矢量z(t)的尺度通常在16至64之间。将矢量z(t)与平整度参考值r(t)相比较以产生差值,即平整度误差e(t)。此平整度误差e(t)被映射p p到传动装置空间以产生所谓的“参数化误差”e(t)。u(t)和e(t)的尺度小于或等于轧机传动装置3的数量。
[0059] 调整包括对加到控制器输出的激励信号进行注入(injection),以及利用所收集p的u(t)和e(t)的值来为每个控制环路估计过程模型的参数。这样做时,使用关于过程的先前所知,故仅需要估计未知参数。未知参数的估计值随后与先前已知的参数一同使用,以形成描述从控制器输出到控制器输入的整个动态特性的完整模型,使得该完整模型能够被用于通过基于模型的调整规则来调整控制器。激励信号被生成直到参数估计足够准确。
[0060] 此调整作给出了要使用的为结果的控制器参数。为每个控制环路执行参数估计、完整模型的形成以及调整。参数估计和调整被在线执行,调整结果也被在线应用以提供自适应控制。自适应控制中有少数参数要估计是有益的。可少到每个控制环路一个参数,即过程增益。不同卷材的真实值可以不同,但同一卷材的真实值始终相同。因此,自适应控制应当在卷材的最初辊轧期间被应用,直到得到增益值。此后,激励信号被切断(设置为零),且不对该控制环路进行进一步的估计,但调整规则仍被用于进行参数调度,这是因为模型会随先前知识如辊轧速度、滤波器传递函数和采样周期被更新而得以更新。
[0061] 在试运行或重调整期间,同样的单元方案应用于自动调整会话,除了被馈送给控制器的控制器参数随后不会因为新的估计参数值而被连续地更新之外。取而代之的是,估计的进展将显示给用户,以允许判断何时认为结果是满意的并结束会话,以及停止应用激励信号。在自动调整会话中,可以使控制器处于手动模式,仅将激励信号馈送给相关的轧机传动装置。
[0062] 图4的流程图公开了被重复执行的、用于基于本发明的平整度控制系统的实例的步骤序列。方框12中所描述的第一行为是从测量系统6得到平整度测量结果。在此情形中测量辊11每旋转一周给出四个测量结果,且测量系统6以此速率或通过整数因子下采样的速率传输新的、可能被滤波的测量结果。
[0063] 随后的行为13是形成或确定作为平整度参考值的偏差的平整度误差。测量结果是矢量,并被从相应的参考矢量中减去,以形成平整度误差。
[0064] 第三步骤14是形成或确定参数化误差,即,被映射到控制环路空间的参数化误差。控制环路的数量小于测量结果矢量的尺度,因此误差被投影为对应于控制环路的数量的更小尺度,使用轧机矩阵,按每控制器输出对带变形进行建模。
[0065] 第四步骤15是更新所需模型参数的在线估计。使用连续时间模型,由于如单独的图中所描述的更少数的参数被在线估计,因此模型的一些部分及模型结构可视为已知。
[0066] 进一步的步骤16是将参数估计与已知部分结合以形成完整的模型。当形成此最终模型时,所估计的参数、视为已知的恒定参数以及可以以已知方式进行变化的参数,如生产速度、采样周期和测量结果滤波(如参数调度)所有这些都要被考虑。
[0067] 另一步骤17是应用调整规则来得到被更新的控制器调整。任何合适的基于模型的调整方法都能够被应用。对于PI控制器,可使用结合了被简化到其所需模型形式(一阶加滞后(first order plus deadtime)的模型的lambda调整。对调整的输入为规定在性能和健壮性之间的所需折衷的最终模型和可调整的参数。
[0068] 第七步骤18是应用该控制,即,利用在该环路对参数化误差起作用的被重新调整的控制器来计算控制器输出。该控制器(如PI控制器)保持像任何正常的控制器一样工作,其仅得到在工作中被更新的调整参数(在PI控制器的情形中的增益和重置时间)。
[0069] 对于每个控制环路,方法步骤15到18被重复,如果用于环路的估计参数或用于调度的参数(速度、采样周期、测量结果滤波器)没有变化,则方法步骤16和17可跳过。
[0070] 下一步骤19是:如果控制环路空间与传动装置空间不同,则从控制环路空间到传动装置空间进行重新映射。如果传动装置被映射为更少数量的控制环路,则控制器输出也被重新映射为传动装置的实际数量。
[0071] 下一步骤20是向传动装置或传动装置伺服器设置点馈送控制器输出。一般每个传动装置都有伺服器,在此情形中控制器输出被馈送给这些伺服器而不是直接馈送给传动装置。
[0072] 图5的框图示出了根据本发明的连续时间动态模型中物理参数递归估计的实施实例的关键行为或步骤。建模的过程具有影响该过程且对于控制该过程可用的输入,且该过程具有输出。其可以是更大过程的一部分,其中该部分被选择,以与分配给它的控制器一起形成控制环路,可能是很多环中的一个环。
[0073] 图5的框图示出了在新的测量结果变得可用时所执行的行动。可以不必为每个抽样都执行这些行动,而是可根据一些标准更少地执行这些行动。本方案指出了一个这样的标准,自适应死区21,意味着当数据中的信息太少时,不执行参数更新。
[0074] 在可以以任何采样速率出现的采样时刻获得测量结果22,过程输入u的变化与这些采样时刻相关联。因此在已过去的许多时刻,在23中y和u成对被获知。这些样本y和u中的一些被集中在要被用于离散时间模型的回归矢量 (t)中,以预测当前的测量结果y(t)。目前的模型被表示为连续时间公式,其一些参数可以被看作是已知参数,而其他的参数将被估计。在时间t,这些参数的先前估计是矢量θ(t-1)。在步骤24中形成完整的连续时间过程模型,该模型假设这些估计为真实值,使用合适的采样周期,该模型在步骤25中被转换为所需的离散时间预测模型。这么做时,任何测量结果滤波器26的作用都会被考虑。
[0075] 在步骤27中,以实际测量结果与预测值的差形成预测误差ε(t)。为了对估计进T行更新,还需要关于估计参数的预测值的梯度ψ(t)。该梯度可通过解析而得出,也可通过步骤28中的数值微分获得。后一手段意味着为一系列其中估计的参数被扰动很小的模型计算预测值。预测值的差除以参数中的扰动得到关于该参数的所需偏导数的近似值。
[0076] 现在,随着预测误差ε(t)和转置预测梯度ψ(t)被得出,其余的参数更新可用标准递归识别来决定。然而,框图中指出了两个重要的实用扩展,异常值保护29和参数估计的限制30。
[0077] 如29中所示,异常值保护对当前预测误差的大小与其统计结果进行比较,并计算充分减小大误差作用的因子。参见图6,其中显示了此信度因子fb的实例。该函数包括对预测误差方差的在线估计,其中要对预测误差与当前大小进行比较。在简单的情形中,通过对均方预测误差进行低通滤波来执行此方差估计。
[0078] 当想要的参数更新达到所允许的任何参数范围之外时,需要进行一些抵消来使估计留在范围之内。这可通过多种不同途径来实现。此处所指出的一个途径是假设真实参数值位于规定的限制之内。如步骤30中所示,指向外部的更新通过用异常值模拟、通过利用公共因子来按比例缩小预测误差和预测梯度来处理。因子的值可被选择成使被更新的估计正好位于允许区域的边界,或者在内部的一些部分或距离上。不需要参数的真实值在允许范围之内的可替选的方式是将违反的参数估计设置在该限制上,并对剩余参数重新进行估计。这可以被有效地执行为剩余参数的尔可夫(Markov)估计,将原始估计视为被白噪声干扰的测量结果,其协方差等于与递归估计相关联的协方差矩阵P。
[0079] 在由于异常值检测或延伸至所允许区域之外的估计而导致的可能的减少之后,预T测梯度ψr(t) 和预测误差εr(t)被用于确定参数估计更新行为31的方向和大小。估计器增益K(t)乘以预测误差εr(t),该增益在步骤32中通过先前时刻的预测梯度和协方差矩阵P来确定。在使用遗忘因子λ(forgetting factor)的情形中,计算为
[0080] K(t)=P(t-1)ψr(t)[λ+ψr(t)TP(t-1)ψr(t)]-1。
[0081] 在对协方差矩阵P的更新行为33中,要考虑几个实际的方面。遗忘这一方面普遍都知道。是否被执行为均匀遗忘(uniform forgetting)或定向遗忘(directional forgetting)或依赖于卡尔曼(Kalman)滤波器方法或在P轨迹上寻址目标会成为要体验的事项。与遗忘有关的还有协方差终结(wind-up)的危险。当信息很少时,在应用合适的自适应死区的情况下,可通过制止更新来进行根本的抵消。协方差矩阵通常被以因数分解的T形式进行更新,例如通过将其写为P=LDL,并更新L(单位对的下三角矩阵)和D(对角矩阵)而不是更新P本身。
[0082] 以下为基于预测误差的递归识别的基本算法,比较: 和Stoica的9.5章,“System Identification”(Prentice Hall,1989)。
[0083]
[0084]
[0085]
[0086]
[0087] 有关参数矢量的预测 及其梯度ψT(t)的表达式将依赖于模型公式和有关扰动特征的假定。对于连续时间模型识别的情形,其中当前预测和响应包括等分布抽样数据,以下方式是可能的。根据利用使用实际采样周期的估计 获得的模型来计算包含相应的离散时间预测模型参数的矢量 所涉及的测量样本值和被操纵的变量(那些预测所需的变量)被集中在矢量 中。随后预测被形成为
[0088]
[0089] 此处,基础假设是:可根据所知的离散时间模型的“arx”结构来对扰动进行建模。(arx=具有控制的自动回归。)“arx”模型常常写为(e(t)代表白噪声):
[0090]
[0091] 用于此模型的最佳的提前一步(one-step-ahead)的预测器为上述(2)中的类型,其中
[0092]
[0093]
[0094] 采用此方式,为形成所需的离散时间参数矢量所要采取的步骤包括:基于所估计的参数和其他已知值,形成连续时间模型,然后以众所周知的方式将其转为离散时间,并将结果模型表示为传递函数的形式,以获得参数a1… b1… 和k。这通过以下实例,即具有部分依赖于传输流率v的滞后的一阶过程来说明。在此实例中,连续时间模型为[0095]
[0096] 在真实情形中,L和流率或速度v可能会被视为已知参数。剩余的参数:增益K、时间常数T以及滞后的独立于速度的部分D0,也被视为已知,或组成所估计的参数矢量0的部分,这意味着这些参数是通过递归估计获得其值的参数。总之形成了模型(5)。采样周期h的对应的离散时间模型为:
[0097]
[0098] 其中α=e-h/τ,β=e-(kh-D)/τ,D=D0+L/v,且k为 圆整到最接近的更大的整数。因此,对应的arx模型具有(3)的形式,其中k如上所给出,且
[0099] a1=-α b1=K(1-β) b2=K(β-α) na=1 nb=2 (7)[0100] 通常情形中,连续时间模型会被最方便地用公式表示为状态空间模型,其首先被转换为离散时间状态空间模型,然后被转换为传递函数形式一仍根据以上所述的实例。该方式允许在完整的模型中将已知参数和未知参数进行任何结合。例如,如果存在应用于原始测量结果的滤波器以得到用于识别的可用采样值,则已知滤波器会在参数a1… b1…的表达式的推导中予以说明。
[0101] 尤其对于简单的模型结构来说,所需预测梯度ψT(t)可以通过解析得出。可替选地,预测梯度也可通过以下方式的数值微分获得。小偏差逐个地被引入每个估计的参数(0的元素),并对每种情况对预测进行计算。对于每个估计的参数,预测中的结果偏差随后除以参数偏差,以形成相应的梯度元素。
[0102] 如果要估计离散时间模型参数,如模型类型(3)的a1… b1… 则从一阶实例的(6)和(7)可看出:任何实际参数增益(K)、时间常数(T)或独立于速度的滞后(D0)的变化会影响估计的参数a1、b1和b2中的至少两个。由于真正的基本参数增益(K)、时间常数(T)或独立于速度的滞后(D0)的变化被期望是彼此独立的,且因为任何用户都可以解释它们的值的含义,所以对它们进行直接估计要好于对a1、b1和b2进行估计。用户将能够判断所估计的值是否是可靠的。所涉及的对连续时间模型参数的估计的计算更复杂,但是当所开发的算法所需的软件实施完成时,利用如今容易得到的硬件是可以解决该问题的。
[0103] 控制器的调整基于过程模型。采样以及特别是应用于测量结果的滤波器的作用也要加以考虑。任何适于所用模型和控制器类型的基于模型的调整方法都可被应用。对于PI控制,一种很简单的方法是已知的、针对具有特定时间常数lambda(λ)的闭环系统的lambda调整。用户通常规定lambda因子或减速(slowdown)因子λf,而不是lambda本身。在其众所周知的基本形式中,lambda调整仅被应用于具有或不具有滞后的一阶系统。在将模型简化为一阶之后,lambda调整也可被应用于更高阶模型。作为简单的模型降阶的实例,具有三个时间常数和滞后的三阶模型
[0104]
[0105] 因此可以利用一阶加滞后模型来近似:
[0106]
[0107] 随后,如果控制器采样周期为Ts,则应用于此模型的lambda调整给出PI控制器增益k以及重置时间TI:
[0108] TI=τ1
[0109]
[0110] 这为具有一个、两个或三个时间常数以及滞后的模型概括了一个lambda调整的变量(在公式(10)中允许T3和T2为零)。(10)中的lambdaλ可以根据lambda因子λf确定为
[0111]
[0112] 如果存在未被包括在过程模型中的测量结果滤波器,则通过将其包括在要被降为一阶的动态模型中来对其进行考虑。例如,具有nf个相等系数的FIR滤波器(有限冲激响应)的情形可以因此利用一阶时间常数Tf=(nf-1)Tb/2来近似,其中Tb为滤波器采样周期。
[0113] lambda调整的准确公式还有其他变量,例如涉及在公式中包含采样周期Ts、λ和λf之间的关系、合并过程的应用,以及同样以降低模型阶数的方式应用高阶模型。在此不对这些变量进行详述。此外,lambda调整仅被用作实例,本发明所描述的方法对于基于模型的调整的其他变量也同样应用良好。
[0114] 要注意的是,虽然以上描述了本发明的示范性实施例,但在不背离所附权利要求限定的本发明的范围的情况下,可以对所公开的解决方案进行一些变化和修改
QQ群二维码
意见反馈