连铸结晶器振动位移系统中非线性处理

李 强 方一鸣,2 李建雄 郑会成

1.燕山大学工业计算机控制工程河北省重点实验室，秦皇岛，066004 2.国家冷轧板带装备及工艺工程技术研究中心，秦皇岛，066004

0 引言

连铸结晶器的非正弦振动是发展高效连铸的关键技术之一[1-2]。利用伺服电机的单方向、变角速度转动来驱动结晶器非正弦振动是一种新方式[3]，具有传动系统简单、结构紧凑、使用寿命长、节能降耗、便于维护等优点。

在工业实际中，连铸结晶器振动系统存在初始时刻偏心轴机械零位的偏移、外部时变扰动等问题，因此常采用结晶器振动位移闭环控制来提高结晶器的跟踪精度，保证结晶器振动的平稳性。伺服电机驱动的连铸结晶器振动位移系统中，系统输出量为结晶器振动位移，执行机构为伺服电机。根据文献[3]的分析可知，结晶器的振动位移曲线近似正弦函数曲线且为中心对称，而伺服电机为单方向转动，电机角位移恒大于0，因此在区间[0，∞)内，伺服电机角位移与结晶器振动位移不是一一对应，增加了控制器设计的复杂性。文献[4]采用T-S模型将输出非线性系统转化为线性时变系统，并考虑系统的约束条件来设计控制器；文献[5]采用反馈线性化的方法对模型进行输入输出线性化，实现系统的解耦和非线性模型的线性化；文献[6]通过构建输出与输入的函数关系，采用极值搜索方法设计控制器。但伺服电机近似为线性系统，结晶器振动位移为非线性，采用上述方法求解跟踪控制问题仍然比较复杂。

针对输出方程为非线性周期函数的问题，文献[7]采用分段函数的处理方法，以结晶器振动周期为界，对不同的区间段，采用不同的控制策略建立结晶器振动位移输出到转角的映射关系。分段函数方法已应用于多种控制器的设计[8-10]。但在实际中，由于采样周期的存在易造成分段区间判定条件的误差，使得分段函数变化提前或滞后，从而造成计算误差，影响控制效果。

本文在分段函数的基础上提出一种基于误差的非线性处理算法，通过建立偏心轴转角误差与结晶器振动位移的映射关系，来保证伺服电机转角到结晶器振动位移非线性关系的唯一性。为验证本文非线性处理算法的有效性，设计了伺服电机驱动的连铸结晶器振动系统位移环PID控制器，并通过仿真和实验研究，对比分析正弦和非正弦两种给定信号情况下的跟踪效果，验证本文所提方法的有效性。

1 伺服电机驱动的连铸结晶器模型分析

1.1 系统模型

伺服电机驱动的连铸结晶器振动系统装置见图1，系统的控制过程为：结晶器振动位移检测装置将结晶器振动位移反馈到数控系统，位置控制器根据期望轨迹信号和实际检测信号分析计算出控制律，控制伺服电机转动，通过减速器、偏心轴和连杆机构实现结晶器振动。控制过程见图2。

图1 连铸结晶器振动系统装置图Fig.1 Diagram of continuous casting mold vibration system

图2 连铸结晶器振动位移系统控制框图Fig.2 Structure drawing of the control system for the continuous casting mold vibration displacement system

伺服电机模型为

(1)

式中，θn为电机转角；ω为转子的角速度；Ud、Uq分别为定子电压的d、q轴分量；Id、Iq分别为定子电流的d、q轴分量；L为定子绕组等效电感；ψf为转子永磁体产生的磁链；Rs为定子电阻；p为电机的极对数；J为转子转动惯量；B为黏性摩擦系数；TL为负载转矩。

由于实际中电流环的采样周期很小，针对d-q轴电流分别采用2个PI控制器来稳定跟踪误差，控制器参数值由伺服控制器自整定，因此，本文主要基于式(1)的前两个方程设计位置环控制器，即

(2)

在实验室的理想情况下，模型参数不考虑外部的影响，均看作常值。

机械传动机构主要包括偏心轴、连杆机构及结晶器振动台，偏心轴连杆结构见图3。

图3 偏心轴Fig.3 Eccentric shaft

根据文献[3]分析可得机械传动部件机理模型：

S≜hsinθm

(3)

式中，h为结晶器振动系统振幅；θm为偏心轴角位移，θm=θn/i；i为减速齿轮减速比。

通过分析式(2)、式(3)可知，在进行位移反馈控制时，构建结晶器振动位移与伺服电机转角之间的关系能有效简化控制器设计。但结晶器振动位移为周期函数，而偏心轴角位移与伺服电机角位移为线性关系，偏心轴转角与结晶器位移间的正弦函数关系在[0，∞)内的逆解不唯一。因此，有效解决结晶器振动位移与偏心轴转角间一一对应的关系，能够简化控制器设计。

1.2 现有方法分析

(4)

图4 偏心轴转角转化过程Fig.4 Conversion process of the eccentric shaft angle

本文在文献[7]分段函数的基础上提出了一种非线性处理算法，通过建立偏心轴转角误差与结晶器振动位移的映射关系，保证伺服电机转角与结晶器振动位移非线性关系的唯一性。

2 连铸结晶器振动位移系统非线性处理算法

针对偏心轴转角与结晶器振动位移的非线性周期函数问题，以及k值变化可能提前或滞后的特点，提出一种基于误差的非线性处理算法。

图5 正弦函数在象限区间中的对应关系Fig.5 Corresponding relation between reference circle and sine function

假设1 结晶器振动位移S及期望轨迹S*是连续可微的，一阶导数和二阶导数存在且有界。

假设2 结晶器振动位移S及期望轨迹S*的相位差小于π/2。

(5)

证明：

首先，根据实际工况和图5所示，结合结晶器振动系统模型，可以得出：

(6)

则综合上述分析可得

(7)

根据式(6)、式(7)可以得出

其次，证明结晶器角位移输出误差的连续性[12]。

假设结晶器振动位移实际信号值滞后于期望信号值，相应的反正弦函数值如图6所示，图中，taj表示期望信号为极值的时刻；tcj表示实际信号为极值的时刻；在taj和tcj时刻之间，为便于分析，选取实际信号与期望信号值相等的时刻为tbj。则根据式(5)的表达式：

图6 利用反正弦角度值求解Δθm的过程Fig.6 Solution process of Δθmbased on the arcsine value

(8)

(9)

(10)

(11)

(12)

(13)

同理可以证明，式(5)在时刻taj、tbj、tcj(j=1,2,…)是连续的。

对于结晶器实际振动位移超前期望信号等情况，可采用类似的方法证明式(5)的连续性，本文不再赘述。证毕。

根据推论1可知，结晶器角位移输出误差可以通过结晶器振动位移和期望位移及其导数的关系来确定。伺服电机角位移是单调递增的，可以推出伺服电机角位移误差：

Δθn=iΔθm

(14)

通过推论1可以看出，在非线性处理算法过程中，可以根据结晶器位移及其导数建立偏心轴角位移误差与结晶器振动位移的映射关系。该算法涉及对结晶器位移求导。工业控制中，求导计算很难实现，且现场复杂环境容易对信号造成干扰，因此，需要构造含有滤波功能且不基于对象模型的微分器。为保证微分器的跟踪性能和滤波性能，本文基于双幂次函数设计了非线性微分器。

引理1 下面的系统相对于原点是有限时间收敛的：

(15)

农田是保证农民经济收入的基础，加强小型农田水利工程建设能够保证农作物的增产和人畜饮水，从而促进农民经济收入和农村经济的发展，但我们也不能一味的追求经济利益而忽略了赖以生存的生态环境，下面我们就现在小型农田水利工程中存在的问题进行分析并提出相关的方法以供农业部门的参考。

证明:

选择Lipschitz连续的Lyapunov函数：

(16)

则有

(17)

给定任意实数σ>0和系数(c1，c2)，假设存在常数g满足如下方程：

(18)

则可以推出

(19)

由于c1，c2>0，且0<α1<1，0<α2<1，所以k<0，则系统平衡点(0，0)是全局有限时间稳定平衡点。

根据引理1可知，基于双幂次函数设计的非线性微分器能够保证在有效时间内跟踪给定量，并同时得出数值的导数。

定理1 下面所示的非线性微分器：

(20)

对于二阶连续可导的信号v(t)，存在μ>0(其中，μα2>2)和Γ>0，当t>ξΓ时，满足

x1-v(t)=O(ξμα2)

证明:

即

根据引理1可知，Lyapunov函数是Lipschitz连续的，则一定存在其Lipschitz常数CL[14]，因此可得

D+(V∘Z)(τ)=

v″(τ)]T]-[z2-a1sig(z1)α1-a2sig(z2)α2]T}≤

通过坐标变换，对于t>ξΓ，可得‖[e1ξe2]‖≤ξμα2pημ即|e1|≤ξμα2pημ，|e2|≤ξμα2-1pημ。

3 仿真和实验

3.1 控制系统的仿真

为验证本文所提非线性处理算法的有效性，首先采用仿真技术对系统的动态性能进行分析。利用MATLAB/Simulink搭建系统仿真模型，对比文献[7]中基于分段函数的非线性处理方法，采用工业中常用的PID控制，分别对正弦和非正弦两种给定信号进行了对比分析。

伺服电机驱动的连铸结晶器振动试验台中机械部分参数如下：减速比i=5.114 5，结晶器振幅h=3 mm；永磁同步电机的参数：PN=20.4 kW，IN=45 A，nN=1 500 r/min，L=4.6 mH，ψf=0.96 Wb，B=0.004，Rs=0.14 Ω，p=3,J=0.054 7 kg·m2。微分器参数a1=a2=2，ξ=0.1，α1=0.4，α2=1.6。

系统整定后的直轴电流PID控制器参数Kp1=12.982 V/A，τ1=2 ms；交轴电流PID控制器的参数Kp2=12.982 V/A，τ2=2 ms；位置环PID控制器参数Kp3=200 V/A，τ3=100 ms。

(1)给定信号为正弦信号S*=hsinωt，其中，ω=πf/30；f为连铸结晶器非正弦振动频率；仿真中取f=13/6 Hz。正弦振动时，负载TL=(5.133 5+ 6.498 5sinωt)N·m。

结晶器启动位置的不确定性使得结晶器存在一定的初始偏差，取初始时刻偏心轴机械零位偏移d=0.2rad。仿真结果如图7所示。可以看出，与分段函数法相比，在控制器参数相同的情况下，本文方法具有较好的跟踪效果，跟踪精度较高，动态性能好。

(2)给定的位移信号为德马克非正弦波形:

(b)位移跟踪误差曲线图7 正弦信号时的结晶器轨迹及误差曲线Fig.7 Mold tracking and tracking error curves with sine signal

非正弦振动的负载为

TL=(5.133 5+6.498 5sin(ωt-Asinωt))N·m

式中，σ为波形偏斜率。

仿真中取f=13/6 Hz，σ=0.24。初始时刻的偏心轴机械零位偏移为0.2 rad。仿真结果如图8所示。可以看出，给定结晶器振动位移为非正弦波形时，本文的非线性处理算法在跟踪速度和跟踪精度方面优于分段函数处理算法，具有较好的控制效果。

(a)位移跟踪轨迹曲线

(b)位移跟踪误差曲线图8 非正弦信号时的结晶器轨迹及误差曲线Fig.8 Mold tracking and tracking error curves with non-sinusoidal signal

3.2 控制系统的实验

为验证本文所提非线性处理算法的有效性，在连铸结晶器振动实验台上进行了相关实验。实验平台由数控部分和机械部分组成，如图9所示。数控部分由西门子的S7-400PLC、FM-458等模块和伺服驱动控制器S120组成，机械部分由减速器、连杆机构及连铸结晶器振动台组成。

1.伺服电机 2.减速器 3.联轴器 4.偏心轴 5.连杆 6.结晶器振动台 7.伺服驱动控制器 8.PLC图9 伺服电机控制柜和连铸结晶器振动台Fig.9 Servo motor control cabinet and vibration platform of continues casting mold

连铸结晶器振动位移控制系统中，伺服电机的电流环由西门子的伺服驱动器S120通过PI控制实现，PI控制器参数可通过S120伺服驱动器自整定。在此基础上，利用本文提出的非线性处理算法对结晶器实际位移跟踪值做映射处理，设计位移环PI控制器。在装有西门子Step7软件的PC机上编写数据采集程序，结晶器的位移通过位移传感器采集到PLC，再上传到上位机中进行处理，通过WinCC监控画面显示。

图10所示为结晶器振动期望轨迹为正弦信号时，闭环控制下结晶器振动位移和速度曲线,图中曲线分别为结晶器正弦振动时的期望轨迹与相应的实际曲线。通过图10可以看出，基于本文非线性处理算法的PI闭环控制能够有效减小结晶器初始相位偏差并将其保持在一定的范围内，实现结晶器振动位移的渐近跟踪控制，表明期望轨迹为正弦信号时，所设计非线性处理算法有效。

图11所示为结晶器振动期望轨迹为非正弦信号时，闭环控制下结晶器振动位移和振动速度曲线,图中曲线分别为结晶器非正弦振动时的期望轨迹与相应的实际曲线。通过图11可以看出，采用非线性处理的PI闭环控制可以实现非正弦情况下结晶器振动位移的跟踪控制，表明所设计的非线性处理方法可以有效解决周期函数逆解不唯一的问题。

4 结论

本文通过建立偏心轴转角误差与结晶器振动位移之间的非线性映射函数关系，来保证伺服电机转角与结晶器振动位移非线性关系的唯一性，并通过设计的PID控制器来验证非线性处理方法的有效性。分别采用MATLAB仿真及伺服电机驱动的连铸结晶器振动实验台，对正弦和非正弦位移给定信号进行实验研究。仿真结果表明，与分段函数方法相比，基于本文非线性处理算法设计的PID控制器在跟踪速度和精度上均有提高，实验结果表明所设计的方法有效。

(a)位移跟踪曲线

(b)速度跟踪曲线图10 结晶器正弦振动时的位移和速度曲线Fig.10 Mold displacement and velocity curves under sine vibration

(a)位移跟踪曲线

(b)速度跟踪曲线图11 结晶器非正弦振动时的位移和速度曲线Fig.11 Mold displacement and velocity curves under non-sinusoidal vibration