CN103076743B - 卷取机的张力模糊pid控制方法 - Google Patents

Abstract

本发明公开了一种卷取机的张力模糊PID控制方法，该方法采用电流环和速度环双闭环控制，对电流环和速度环的反馈环节分别设置速度和电流调节器，并采用六个模糊控制器控制电流环和速度环的PID参数，每个模糊控制器采用两输入一输出，输入量均为设定张力的偏差值E以及偏差变化率Ec，输出量分别为电流控制器和速度控制器需要的P、I、D参数；在电流调节器与速度调节器之前分别串入3个相应的模糊控制器，在卷取机运行的过程中，不断调整隶属度函数，对电机转速和电枢电流进行控制，实现对卷取机张力的自动控制，使卷取机的张力实现自动调节始终保持恒定；通过simulink仿真，证实该方法响应快，调整能力强，鲁棒性好，有效地改善了控制效果。

Description

卷取机的张力模糊PID控制方法

技术领域

本发明涉及一种卷取机张力控制方法，通过对模糊控制器对卷取机的电流环及速度环的PID参数控制，实现卷取机张力间接控制。本发明可广泛应用于纤维、造纸、塑料薄膜、电线、印刷品、磁带等轻工业相关领域以及金属制品行业。

背景技术

在冶金行业中，卷取机是将热轧或冷轧钢材卷取成卷筒状的轧钢车间辅助设备，用卷取机将钢材弯曲成卷，从而为增大原材料重量、提高轧制速度、减小轧件头、尾温差提供了有力的条件，由此导致了产品产量与质量的提高；此外，成卷的轧材便于运送，这是各种形式卷取机的共同特点和作用。

卷取机一个重要的控制就是张力控制，在冶金行业中，对薄钢板的卷取就会碰到类似问题。如果张力过大，钢材在卷取的过程中，便会使钢卷的内卷产生相对滑动，由滑动造成的表面划痕，极大地影响了钢材的表面质量，更为重要的是，张力过大会造成钢卷内部应力变大，最终使钢卷的内孔内陷，或者使整个内层卷突出，形成塔形，最终使钢卷的等级降低。如果卷取的张力过小，钢卷会在其自身重力的情况下松散、突出，无法对其包装。而且在松散过程中，由于钢卷的相对滑动，常常会在钢材表面产生大面积的划痕。

解决上述问题，工业现场传统方式通过变频器对卷筒电机进行控制。通常采用PLC对电机进行PI（比例、积分）调节，或者PID（比例、积分、微分）调节。传统的PID控制算法，其设计简单、实施容易，具有很强的鲁棒性，特别对二阶系统具有最优的调节特性，因此在工业控制中一直占有主导地位。但是，张力系统中，常常有很多干扰因素，PI或者PID调节系数，如K_P(比例系数)K_I（积分系数）K_D（微分系数）一旦设定，如果没有人为对其进行调节，这些参数是固定不变的，它们不会随着受控参量的变化而变化，即不会随着环境的改变而改变。但实际上，在张力控制过程中，张力的大小是随着卷径的变化而实时变化的，当卷径较小时，卷辊线速度较小，张力不大。随着卷径的不断扩大，卷辊的线速度变快，张力变大。因此设置静态参数往往达不到最佳控制效果，PID控制常常达不到要求精度。

纵观国内外的张力控制设备的发展，它的研究重点是对张力控制的算法，目前已经有一些控制算法，神经网络控制算法，不依赖对控制对象的数学模型，但是，由于神经网络需要大量的样本且需要对网络样本进行训练，而训练过程比较繁琐，不适宜卷取机这类对控制实时性有较强要求的系统；分析卷取机这一类系统，它是多参量、非线性的，不宜使用精确的数学方法控制。这样的系统，恰恰最适合使用模糊控制技术进行控制。模糊控制是一种基于规则的控制，它直接采用语言控制规则，出发点是模糊语言和相关专家的知识，在设计中不需要建立被控制对象的精确的数学模型，控制机理和策略易于接受与理解，设计简单，便于应用。

发明内容

针对上述现有技术存在的缺陷或不足，本发明的目的在于，提供一种卷取机的张力模糊PID控制方法。本申请首先对卷取机工作特性进行分析，完成卷取机主电机的选型；然后建立卷取机的数学模型，完成卷取机的工作仿真；设计卷取机张力控制的模糊控制器，该模糊控制器不是建立在数学分析基础上，而是根据工程实际运行的样本确定各参数和控制规则，最后将模糊控制器应用于实际卷取机系统，根据控制结果调整输入、输出隶属度函数，直至调整到理想效果。

为了实现上述目的，本发明采用如下的技术解决方案：

一种卷取机的张力模糊PID控制方法，该方法采用电流环和速度环双闭环控制，对电流环和速度环的反馈环节分别设置PID调节器，即速度调节器和电流调节器，并采用六个模糊控制器分别控制电流环和速度环的PID参数的调整，每个模糊控制器采用MISO的形式即两输入一输出，输入量均为设定张力的偏差值E以及偏差变化率Ec，输出量分别为电流控制器和速度控制器需要的P、I、D参数；在电流调节器与速度调节器之前分别串入3个相应的模糊控制器，在卷取机运行的过程中，不断调整隶属度函数，对电机的转速和电机电枢电流进行控制，从而实现对卷取机张力的自动控制，使卷取机的张力T实现自动调节始终保持恒定；具体操作如下：

用户根据公式T=σ₀bh设定所需张力T：其中，b为带钢的宽度，单位：mm；h为带钢的厚度，单位：mm；σ₀为带钢的单位张力，单位：N/mm²；

根据公式M_T=T*(D/2)计算得到张力力矩M_T，单位：Nm；其中，T为张力设定值单位：N，D为卷筒卷径D，单位：m；计算动态力矩M_D，单位：Nm；通过实验得到空载力矩M_o，单位：Nm；将张力力矩M_T、动态力矩M_D、空载力矩M_o相加得到主电机的力矩极限值M_m；

在主电机M上安装脉冲编码器2测得卷筒的转速n，通过S辊的线速度V1计算卷筒的线速度V2，将卷筒的线速度V2与卷筒的转速n经过除法器得到卷筒卷径D；

S辊的线速度V1与S辊的直径D_S经过除法器得到转速给定值，转速给定值与转速n经过比较器得到速度偏差；卷筒卷径D与带钢带宽共同作用的转动惯量经过比例调节器，所述比例调节器用于根据卷筒卷径D与带钢带宽的变化引起的动态惯量变化来调节比例系数；将速度偏差和比例系数输入速度调节器得到速度整定值，速度整定值作为电流给定值；将电流给定值和力矩极限值Mm送入电流调节器，经电流调节器通过如下公式得到主电机M的电枢电流值，使主电机M根据得到的电枢电流I_d运转；

Mm=C_mφI_d

其中，M_m为主电机的力矩极限值；φ为磁通量，单位：Wb；C_m为电动机力矩常数，取2.68N/m·A；I_d为电枢电流，单位：A。

本发明还包括如下其他技术特征：

所述S辊的线速度V1的计算如下：

S辊的驱动电机M1通过减速器1连接S辊，在S辊的驱动电机M1安装脉冲编码器1检测S辊的转速n1，利用下式计算得到S辊的线速度V1：

其中，D_S——S辊的直径，单位m；

n₁——S辊的实际转速，单位r/min；

π————圆周率。

所述比例调节器根据系统卷径带宽的变化引起的动态惯量变化，来调节速度条件器的比例系数的计算公式如下：

其中，

K_n——比例调节器的比例系数；

h——频率宽度，为了跟随特性和抗干扰性，h取5；

β——电流反馈系数，单位：V/A取0.0095V/A；

C_e——主电机的电势转速比，单位：V/(r.min^-1)取0.185V/(r.min^-1)；

T_m——主电机机电时间常数，单位：s，空卷时T_m=0.148s，满卷时T_m＝0.9s；

——惯性时间常数，单位：s，取0.0124s；

α——转速反馈系数，单位：s取0.001s；

R——电枢回路等效电阻。单位：欧姆；

所述的模糊控制器的设置如下；

步骤1：对每个模糊控制器的输入和输出均设置七个语言值，采用Gauss型函数，并确定函数的边界点和中心点，完成所有控制语言的隶属度函数；

步骤2：建立模糊规则表，即确定输入和输出之间的控制规则；

步骤3：对输出的模糊语言利用模糊规则表解模糊。

所述的步骤3采用最大隶属度法解模糊。

在解模糊过程中，如果推理出来的论域有多个元素同时出现最大隶属函数的值，就取这些元素的平均值作为该隶属函数的清晰值。

与现有技术相比，本发明具有如下优点：

1、以人的丰富实践经验和思维过程，构建的模糊规则进行推理和判断。它无需对被控对象建立精确的数学模型，只需做模糊描述即可。易于实行和操作。

2、模糊控制是通过引入模糊逻辑语言变量，通过各变量间的模糊关系进行推理，从而使计算机控制进入那些基于精确模型无法控制的禁区。

3、更精确的实现了张力控制，使设备更加简单，人性化，减少了生产的废品率，降低了成本。

附图说明

图1是本发明的卷取机张力控制系统的结构原理图。

图2是模糊PID控制原理图。

图3是常规的PID张力控制数学模型。

图4是本发明的模糊PID张力控制数学模型。

图5是常规PID控制响应曲线速度波形图。

图6是本发明进行模糊PID控制的响应曲线速度波形图。

以下结合附图和具体实施方式对本发明进一步解释说明。

具体实施方式

1卷取机的建立

1.1卷取机间接张力控制的实现原理

卷取机卷绕带钢时，S辊的线速度设为v1，卷筒上的线速度设为v2，在开始卷取时，v2显然大于v1的，只有这样，开始卷取时才能使带钢紧紧缠绕在卷筒上。否则便会使带钢松散。根据胡克定律F=-kx得到带钢的张力T的计算公式：

其中A——带钢横截面积，单位：mm²；

K——带钢的弹性系数，单位：N/mm；

l——卷筒与辊之间的距离，单位：mm；

t——卷取机建立张力的时间，单位：s；

K₀——

从公式可知，只有当v2大于v1时，被卷取的带钢便会产生弹性形变，从而产生张力T；v2与v1差值为固定值时张力恒定不变。v2与v1中的任一个发生变化都会导致张力不稳定，造成卷取出现质量问题。因此，在卷取机工作过程中应保持v2与v1的速度差为定值。

其中，C——常数，由卷取机系统本身特性决定；

E——主电机的电动势；单位：V；

V——卷筒的线速度；单位：m/s；

I_d——主电机的电枢电流。单位：A；

从上式可知，要求张力大小T不变，那么线速度v2是恒定的，而电动势不变，张力T可以用电枢电流线性表示；如果主电机的电动势E不变，而I_d与V的比也为定值，那么张力T也为一定值；根据公式

其中，C_m——电动机力矩常数，单位：N/m·A；

D——卷筒卷径，单位：m；

φ——磁通量，单位：Wb；

当磁通φ不变时，Id随着卷筒卷径D的变化而变化，并且要与卷筒卷径D的变化保持一致，使张力T达到一个恒值，或者保持I_d不动，而磁通φ随着卷筒卷径D变化，从而可以通过控制主电机的电枢电流I_d实现对卷取机张力的控制。

因为卷取过程中，卷径D不断增大，而线速度不变，但这样便会增大力矩，为了保持卷取机张力T为恒定值，须随着卷筒卷径D的变化调节电机速度，降低转速，达到最终的精确控制。因此，可以通过控制主电机的转速实现对卷取机张力的控制。

综上，通过对电机电枢电流或者对电机的转速的控制，就能实现对卷取机张力的间接控制。因此，卷取机采用电流环和速度环双闭环控制。设计双闭环控制能够获得理想的启动过程。使用两种负反馈分别对电流环和速度环进行调节，对这两个反馈环节设置PID调节器。

1.2卷取机分析

本发明需要设计一种间接张力控制卷取机，应用最大力矩式间接控制，通过对整个卷取张力系统的物理方程进行动态、静态分析，得出一系列参数，达到间接控制卷取机张力的目的。

如图1所示，本发明的方法所应用的卷取机包括可编程控制器S7-300PLC、西门子6SE70变频器、现场总线Profibus-DP、脉冲编码器1、脉冲编码器2、减速器1、减速器2、主电机M、S辊的驱动电机M1和卷筒。

该卷取机具体工作原理如下：

用户根据需要设定卷取机的张力值，张力设定值T（单位：N）与卷筒卷径D的实测值根据公式M_T=T*(D/2)计算得到张力力矩M_T。卷筒卷径D、卷筒转速n的变化产生动态力矩M_D；通过实验得到空载力矩M_o；将张力力矩M_T、动态力矩M_D、空载力矩M_o相加得到主电机的力矩极限值M_m。

主电机M通过减速器2驱动卷筒转动卷取带钢，减速器2选用卧式两级圆柱斜齿轮，在卷筒转动过程中，将来自S辊的带钢经过导向辊缠绕在卷筒上，在主电机M上安装脉冲编码器2测得卷筒的转速n，根据胡克定律F=-kx，只有卷筒的线速度V2与S辊的线速度V1的差值是定值，被卷曲的带材产生固定的形变，才能保证张力恒定。由张力设定值T与带材弹性系数、带材横截面积、卷筒与辊之间距离可建立卷筒的线速度V2与S辊的线速度V1的差值。通过S辊的线速度V1计算卷筒的线速度V2，也可直接在导向辊上安装速度传感器测出卷筒的线速度V2，将卷筒的线速度V2与卷筒的转速n经过除法器得到卷筒卷径D。

S辊的驱动电机M1通过减速器1连接S辊，在S辊的驱动电机M1安装脉冲编码器1检测S辊的转速n1，利用下式计算得到S辊的线速度V1：

其中，D_S——S辊的直径，单位m；

n₁S辊的实际转速，单位m/min；

π——圆周率。

S辊的线速度V1与S辊的直径D_S经过除法器得到转速给定值，转速给定值与转速n经过比较器得到速度偏差；为了保持主电机转速尽可能与给定转速一致，达到系统无静差，转速环节应被校正为II型系统，因此，在卷筒卷径D与带宽共同作用的转动惯量需要经过比例调节器，比例调节器根据系统参数的变化调节比例系数K_n，其计算公式如下：

其中，

K_n——比例系数；

h——频率宽度，为了跟随特性和抗干扰性，h取5；

β——电流反馈系数，取0.0095V/A；

C_e——主电机的电势转速比，取0.185V/(r.min^-1)；

T_m——主电机机电时间常数，单位：s，空卷时T_m=0.148s，满卷时T_m＝0.9s；

——惯性时间常数，单位：s，取0.0124s；

α——转速反馈系数，单位：s，取0.001s；

R——电枢回路等效电阻，单位：欧姆；

将速度偏差和比例系数输入速度调节器得到速度整定值，速度整定值作为电流给定值；将电流给定值和力矩极限值Mm送入电流调节器，经电流调节器得到主电机M的电枢电流值，主电机M根据得到的电枢电流值运转。

所述的速度调节器和电流调节器均为PID控制器，速度调节器和电流调节器设置在变频器中。

所述的现场总线Profibus-DP实现PLC与变频器的通讯。

1.3主电机的选择

在卷取机张力控制中，设计两个PI调节器：电流调节器和速度调节器。其中电流调节器是当卷径D在增加过程中，要建立张力电流i_d1、空载补偿电流i_d2、动态力矩补偿电流i_d3，由此完成系统的动态补偿，才能保证系统卷取电机的电枢电流恒定。

主电机M作为卷取机的动力源，必须产生足够的力矩以驱动系统。从图1可知卷筒有三个力矩：张力力矩M_T、动态力矩M_D、空载补偿力矩M_o，因为卷取机工作时三个力矩相加得到力矩极限值Mm，使其作为速度环的输入，选取的主电机工作应同时满足理论力矩极限和电流极限。卷取机的主电机力矩平衡方程为：

Mm=M_T+M_D+M_o

Mm=C_mφI_d

Mm——主电机的电磁力矩，单位：Nm；

M_T——产生张力所需的张力力矩，单位：Nm；

M_D——动态力矩，主电机加速为正值，主电机减速为负值，单位：Nm；

M_o——空载力矩，单位：Nm；

φ——磁通量；单位：Wb

C_m——电动机力矩常数，由电动机自身性能决定，单位：N/m·A；

I_d——电枢电流。单位：A

通过上式得到主电机的力矩极限值Mm即电磁力矩值Mm，进而得到电枢电流I_d；选取满足电枢电流I_d的合适主电机。本发明中C_m为2.68N/m·A。

1.3.1张力力矩M_T的计算

产生张力所需的张力力矩是在卷筒启动的时候对带钢所产生的力矩，此张力在建立时与所卷取的带材有关。本发明选取带钢的宽度b为1000mm，厚度h为0.5mm，经测定，带钢的单位张力σ₀=140N/mm²，用户依据张力计算公式T=σ₀bh计算出负载的理论张力T为35000N其中，b为带钢的宽度，单位：mm；h为带钢的厚度，单位：mm；σ₀为带钢的单位张力，单位：N/mm²；。卷取机的张力力矩M_T=T*(D/2)，其中，T为张力设定值，单位：N；D为卷筒的卷径，单位：m。

1.3.2动态力矩MD的计算

卷取机在工作过程中，卷径D的变化直接影响带钢的线速度V2、动态力矩补偿电流i_d3；卷径D随着卷筒不断卷取而逐渐增大，转动惯量也随着增大，那么卷筒上的线速度V2变化必然跟不上辊的线速度V1的变化，从而产生动态力矩M_D，必须对卷筒的变化的转动惯量进行补偿。动态力矩M_D的计算过程如下：

式中GD₁ ²——主电机轴转动惯量；单位：N·m²

GD₂ ²——机械传动转动惯量；单位：N·m²

GD₃ ²——卷取机卷筒转动惯量；单位：N·m²

n——卷筒转速；r/min

i——传动比。

或者根据如下公式计算M_D：

式中，

GD²GD₁ ²、GD₂ ²、GD₃ ²三个值之和；单位：N·m²

g——重力加速度；单位：m/s²

v——卷筒线速度；单位：m/min

D——卷筒实时卷径；单位：m

h——带材厚度；单位：m

b——带材宽度，单位m；

S——带材的填充率，与带材致密程度有关，一般为0.85~0.9；

ρ——带材的比重，单位kg/m³，与带材材料有关；

GD₁ ²为主电机轴转动惯量，与电机本身参数有关；GD₂ ²为机械传动（主要为减速器及其连接装置）转动惯量；GD₁ ²、GD₂ ²视为已知定量；GD₃ ²的计算公式如下：

GD₃ ²＝2πρSb(r⁴-r₀ ⁴)

式中，

ρ——带材的比重，单位kg/m³，与带材材料有关；

S——带材的填充率，与带材致密程度有关，一般为0.85~0.9；

b——带材宽度，单位m；

r——带卷半径，单位m；

r₀——卷筒半径，单位m；

由公式可知，GD₃ ²由带材的宽度、密度、卷筒的直径或半径相关，并不恒定。将GD₁ ²、GD₂ ²、GD₃ ²三个值相加，通过上述公式，即可得到卷取机动态力矩M_D。

1.3.3空载力矩M_o的计算

空载力矩M_o主要克服机械摩擦产生的力矩，卷取机自动过程中克服的最大静摩擦力所需的电流可通过点动方式进行测量，先将电流限幅至零，逐渐加大，直到电机启动，这样便可以得到卷取机空载所需的补偿电流，进而得到空载力矩M_o。

在本实施例中，卷筒内径D=508mm；卷筒外径D=1500mm；卷筒宽度700～1250mm；带钢厚度0.2mm~1.2mm；卷取速度传动比i=16.88；工作电压380v；卷取速度max130m/min；上述已经计算出主电机所需最大力矩值Mm和电枢电流I_d，选取的主电机M应满足最大力矩和最大电流要求：主电机参数220V，额定功率90KW，主电机转速750r/min，励磁回路电阻Rx=0.08Ω，励磁回路电感I_∑=0.004H，电枢回路电阻R_a=0.02Ω，过载倍数λ=1.5。

2卷取机模糊PID的设计

2.1自整定设计思想：

PID参数自整定设计思想是先找出PID三个参数Kp、Kl、K_D与偏差|E|和偏差变化率|E_C|之间的模糊关系，在工作中不断检测偏差值|E|和偏差变化率|E_C|，再根据模糊控制原理对参数Kp、Kl、K_D进行在线校正，以满足不同的偏差值|E|和|E_C|时对控制参数的不同要求，而使被控对象有良好的动、静态性能，而且计算量小，便于微机实现。

（1）在偏差值|E|较大，即系统处于响应的初始段，为了使系统响应具有良好的快速跟踪性能，并避免因开始时偏差的瞬间变大，可能引起微分过饱和，而使控制作用超出许可范围，应取较大的Kp和较小的K_C，同时为了避免系统响应出现较大的超调，需对积分作用加以限制，通常取Kl=0。

（2）当偏差值|E|为中等大小时，即系统响应处于中间阶段，为了使系统具有较小的超调，应取较小的Kp，适当的Kl和K_D，以保证系统的响应速度。

（3）当偏差值|E|较小时，即系统处于响应处于结束阶段，为使系统具有较好的稳定性，应取较大的Kp和Kl，同时为避免系统在设定值附近出现震荡，并考虑系统的抗干扰性能，K_D值可取大些，通常为中等大小，当|Ec|较大时，K_C值可取小些。

建立参数Kp、Kl、K_D与偏差值|E|和偏差变化率|E_C|间的二元函数关系：

通过上式计算PID参数Kp、Kl、K_D，即可得到控制器在第n个采样时刻的输出量U（n）：

其中，Kp为比例系数；Kl为积分系数；K_D为微分系数；e（n）是第n个时刻的输出采样值，e（n-1）是第n-1个时刻的输出采样值。

2.2模糊控制器的设计：

把模糊控制理论与PID控制算法结合起来构成模糊PID控制器，如图2所示的自整定模糊PID控制器。本发明的卷取机张力控制的模糊PID自动控制方法，采用六个模糊PID控制器，分别用于控制电流环和速度环的PID参数的调整。模糊控制器设置在可编程控制器PLC中，在电流调节器与速度调节器之前分别串入相应的模糊控制器。每个模糊控制器均为两输入一输出（MISO）结构，输入为卷取机设定的张力偏差值E与张力偏差值变化率E_C。将输入的具体值模糊化，采用mamdani算法对模糊过程进行推理，输出量分别为电流控制器和速度控制器需要的P、I、D参数。

2.2.1确定隶属度函数

每个模糊控制器均为两输入一输出（MISO）结构，输入为卷取机设定的张力偏差值E与张力偏差值变化率E_C，将输入的两个精确量利用隶属度函数模糊化，作为PID动态过程中现象差异的中间过渡，在实际卷取机系统中，不断调整隶属度函数，直到恒张力控制效果理想为止，本发明中最终确定了使用Gauss型隶属度函数，统计表明，自然界众多现象都服从正态分布，因此采用Gauss型做为隶属度函数是合适的。

将期望张力设定在1000N，由工程实际基本可以确定，张力偏差值E输入的基本论域为【-270N,270N】，张力偏差变化率EC基本论域为【-60N/0.1S,60N/0.1S】。上述输入量的越界值均作为边界值处理。论域内的精确量通过隶属度函数转变为语言值{负大、负中、负小、零、正小、正中、正大}，用英文字母缩写为【NB、NM、NS、ZO、PS、PM、PB】。输入输出的隶属度函数参见表1-表5。

表1输入变量E（张力偏差值E）的隶属度函数

语言值	[下边界，中点，上边界]
		NB	[-405，-270，-135]
NM	[-315，-180，-45]
		NS	[-225，-90，45]
ZO	[-135，0，135]
		PS	[-45，90，225]
PM	[45，180，315]
		PB	[135，270，405]

表2输入变量E_C（张力偏差变化率E_C）的隶属度函数

语言值	[下边界，中点，上边界]
		NB	[-90，-60，-30]
NM	[-70，-40，-10]
		NS	[-50，-20，10]
ZO	[-30，0，30]
		PS	[-10，20，50]
PM	[10，40，70]
		PB	[30，60，90]

表3输出变量P的隶属度函数

语言值	[下边界，中点，上边界]
		NB	[15，30，45]
NM	[25，40，55]
		NS	[35，50，65]
ZO	[45，60，75]
		PS	[55，70，85]
PM	[65，80，95]
		PB	[75，90，105]

表4输出变量I的隶属度函数

语言值	[下边界，中点，上边界]
		NB	[1.5，3，4.5]
NM	[2.5，4，5.5]
		NS	[3.5，5，6.5]
ZO	[4.5，6，7.5]
		PS	[5.5，7，8.5]
PM	[6.5，8，9.5]
		PB	[7.5，9，10.5]

表5输出变量D的隶属度函数

语言值	[下边界，中点，上边界]
		NB	[-0.2，0.4，1]
NM	[0.2，0.8，1.4]
		NS	[0.6，1.2，1.8]
ZO	[1.0，1.6，2.2]
		PS	[1.4，2.0，2.6]
PM	[1.8，2.4，3.0]
		PB	[2.4，2.8，3.4]

2.2.2模糊规则表的建立

在每个PID控制器中，比例系数Kp值决定系统的响应速度。增大Kp能提高响应速度，减小稳态误差；但Kp值过大会产生较大的超调，甚至使系统不稳定；减小Kp可减小超调，提高稳定性，但Kp过小会减慢响应速度，延长响应时间。因此调节初期应适当取较大的Kp值以提高响应速度，而在调节中期，Kp则取较小值，以使系统具有较小的超调并保证一定的响应速度；而在调节后期，再将Kp值调到较大值来减小静差，提高控制精度。

表6Kp模糊控制规则表

在PID控制器中，积分控制Kl主要用来消除系统的稳态误差。由于饱和非线性等原因，积分过程有可能在调节初期产生积分饱和，从而引起调节过程的较大超调。因此，在调节过程初期，为防止积分饱和，其积分作用Kl应当弱一些，甚至可以取零；而在调节中期，为了避免影响稳定性，其积分作用Kl应该比较适中；最后在过程后期，则应增强积分作用Kl，以减小调节静差。

表7Kl模糊控制规则表

微分环节的调整主要针对大惯性过程引入的，微分环节的系数K_D的作用在于改变系统的动态特性。系统的微分环节系数能反映信号变换的趋势，并能在偏差信号变换太大之前，在系统中引入一个有戏的早期修正信号，从而加快响应速度，减少调整时间，消除振荡，最终改变系统的动态性能。因此，微分环节的系数K_D值的选取对调节动态特性影响很大。微分环节的系数K_D过大，调节过程制动就会超前，致使调节时间过长；微分环节的系数K_D过小，调节过程的制动就会落后，从而导致超调增加。根据实际过程经验，在调节初期，应加大微分作用，这样可得到较小甚至避免超调；而在中期，由于调节特性对K_D值得变化比较敏感，因此，K_D值应适当小一些并应保持固定不变；在调节后期，K_D值应减小，以减小被控制过程的制动作用，进而在调节过程初期由于K_D值较大所造成的调节过程的时间延长。

表8K_D模糊控制规则表

2.2.3解模糊化

模糊控制器通过模糊规则表，得到一组语言化的控制输出，需要对这些语言的模糊量精确化，也就是解模糊过程，这里采用最大隶属度法进行数据清晰化，即选取隶属度最大的元素，把其作为清晰量的方法。如果推理出来的论域（即输出的模糊子集）中，有多个元素同时出现最大隶属函数的值，那么就取这些元素的平均值作为该隶属函数的清晰值。

完成上述工作后，对卷取机张力控制设计已经完成。

3.卷取机模型建立

以下通过建立卷取机数学模型，来说明本发明的方法的有效性。卷取机常规PID数学模型如图3所示，在该模型的基础上，添加模糊PID控制器得到如图4所示的模糊PID数学模型。

对数学模型中的各参数进行选取：

晶闸管滞后时间常数取T_S=0.0017s;

电流环的滤波时间常数T_oi=0.002s;

晶闸管的比例系数其中取电流调节器的限幅值U_ct=8V，电压U_d0=220V，那么K_S=27.5。

转速反馈系数其中选取U_m=8V，则反馈系数α=0.011;

对于速度调节器分为满卷时和空卷时其中有：

电流环电流调节器的时间常数τ_i=T_n，其中取速度环滤波时间常数K_d=5，则

对于转动惯量，在电机空卷过程中，卷筒和卷取电机的转动惯量折算到电机的轴上的惯量为GD²=1962N·m²，转动惯量

当卷取机工作时，卷入卷筒上得带钢和卷取电机的转动惯量折算到电机的轴上的惯量为GD²=11772N·m²,，转动惯量其中带钢卷取后的直径即卷筒卷径D为1.5m。

电机的电势转速比：

电机力矩电流比：

因此电机在卷筒空卷时：

空卷时T_m＝0.148s

卷筒满卷时：

满卷时T_m＝0.9s

1.6动态仿真

为了验证本发明方法的有效性，通过Simulink建立数学模型，进行卷取机张力控制常规PID动态仿真，得到如图5所示的速度波形图。从图5看出，常规PID的超调量与调节还不是很理想，还有继续优化的空间。

利用模糊PID控制器，采用simulink对卷取机数学模型进行卷取机张力控制仿真，得到如图6所示的速度波形图。在电流调节器与速度调节器之前各串入相应的模糊控制器进行仿真，用以检验设计效果，仿真数学模型如图4所示。

参见如图5和图6，通过常规PID控制与模糊PID控制速度环波形对比，可明显看出，模糊PID超调量与快速量比常规PID控制明显优化了很多，系统运行更加快速及稳定，达到了预期效果。

Claims (1)

1.一种卷取机的张力模糊PID控制方法，其特征在于，该方法采用电流环和速度环双闭环控制，对电流环和速度环的反馈环节分别设置PID调节器，即速度调节器和电流调节器，并采用六个模糊控制器分别控制电流环和速度环的PID参数的调整，每个模糊控制器采用MISO的形式即两输入一输出，输入量均为设定张力的偏差值E以及偏差变化率Ec，输出量分别为电流控制器和速度控制器需要的P、I、D参数；在电流调节器与速度调节器之前分别串入3个相应的模糊控制器，在卷取机运行的过程中，不断调整隶属度函数，对电机的转速和电机电枢电流进行控制，从而实现对卷取机张力的自动控制，使卷取机的张力T实现自动调节始终保持恒定；具体操作如下：

用户根据公式T=σ₀bh设定所需张力T：其中，b为带钢的宽度，单位：mm；h为带钢的厚度，单位：mm；σ₀为带钢的单位张力，单位： N/mm²；

根据公式M_T=T*(D/2)计算得到张力力矩M_T，单位：Nm；其中，T为张力设定值单位：N，D为卷筒卷径D，单位：m；计算动态力矩M_D，单位：Nm；通过实验得到空载力矩M_o，单位：Nm；将张力力矩M_T、动态力矩M_D、空载力矩M_o相加得到主电 机的力矩极限值M_m；

在主电机(M)上安装第二脉冲编码器(2)测得卷筒的转速n，通过S辊的线速度V1计算卷筒的线速度V2，将卷筒的线速度V2与卷筒的转速n经过除法器得到卷筒卷径D；

S辊的线速度V1与S辊的直径DS经过除法器得到转速给定值，转速给定值与转速n经过比较器得到速度偏差；卷筒卷径D与带钢带宽共同作用的转动惯量经过比例调节器，所述比例调节器用于根据卷筒卷径D与带钢带宽的变化引起的动态惯量变化来调节比例系数；将速度偏差和比例系数输入速度调节器得到速度整定值，速度整定值作为电流给定值；将电流给定值和力矩极限值M_m 送入电流调节器，经电流调节器通过如下公式得到主电机(M)的电枢电流值，使主电机(M)根据得到的电枢电流Id运转；

Mm＝CmφId

其中，M_m 为主电机的力矩极限值；φ为磁通量，单位：Wb；Cm为电动机力矩常数，取2.68N/m·A；Id为电枢电流，单位：A；

所述的模糊控制器的设置如下；

步骤1：对每个模糊控制器的输入和输出均设置七个语言值，采用Gauss型函数，并确定函数的边界点和中心点，完成所有控制语言的隶属度函数；

步骤2：建立模糊规则表，即确定输入和输出之间的控制规则；

步骤3：对输出的模糊语言利用模糊规则表解模糊；

所述的步骤3采用最大隶属度法解模糊；在解模糊过程中，如果推理出来的论域有多个元素同时出现最大隶属函数的值，就取这些元素的平均值作为该隶属函数的清晰值；

S辊的驱动电机(M1)通过第一减速器连接S辊，在S辊的驱动电机(M1)安装第一脉冲编码器(1)检测S辊的转速n1，利用下式计算得到S辊的线速度V1：

其中，DS——S辊的直径，单位：m；

n₁——S辊的实际转速，单位：r/min；

π——圆周率；

所述比例调节器根据系统卷径带宽的变化引起的动态惯量变化，来调节速度调节器的比例系数，给速度调节器叠加一个5％的给定值，使速度饱和；饱和的速度使电流环限幅输出，当速度饱和后，力矩极限接管系统控制；比例系数的计算公式如下：

其中，K_n——比例调节器的比例系数；

h——频率宽度，为了跟随特性和抗干扰性，h取5；

β——电流反馈系数，单位：V/A，取0.0095V/A；

C_e——主电机的电势转速比，单位：V/( r.min^-1)，取0.185V/( r.min^-1)；

T_m——主电机机电时间常数，单位：s，空卷时Tm＝0.148s，满卷时Tm＝0.9s；

T_∑n——惯性时间常数，单位：s，取0.0124s；

α——转速反馈系数，单位：s，取0.001s；

R——电枢回路等效电阻，单位：欧姆。