板带轧机液压压下-垂直振动特性研究

侯东晓， 方 成， 陈善平， 阎 爽

(1.东北大学秦皇岛分校 控制工程学院， 河北 秦皇岛 066004； 2.京能秦皇岛热电有限公司， 河北 秦皇岛 066004)

板带轧机中的辊系振动问题一直是影响钢铁企业正常轧制生产的技术难题，辊系异常振动的发生不仅影响着轧制产品的精度与表面质量，严重时还可能引发钢带断裂并损坏设备[1-3].

针对轧机振动问题，学者们从不同的角度进行了研究[4-8].Yun等[9]和Hu等[10-11]考虑轧制角与咬入角对轧制稳定性的影响，建立了水平-垂直-扭转耦合模型.Johnson等[12]考虑了接触非线性对工作辊-支承辊界面动态特性的影响，指出接触非线性会产生高频谐波振动.Swiatoniowski等[13-14]考虑轧件变形的非线性弹性力和参激振动，研究了轧机辊系垂直振动特性.侯东晓等[15]考虑压下缸和平衡缸多个缸体间结构弹性约束影响，建立了分段非线性弹性约束下轧机振动模型，并研究了该模型的分岔特性.学者们在此基础上进一步考虑多个液压缸的弹性约束影响：刘浩然等[16]研究了多非线性弹性约束下的轧机振动；刘飞等[17]进一步研究了液压缸在分段弹性力与非线性摩擦力共同作用下的轧机振动特性；刘彬等[18]针对含液压缸非线性弹簧力的轧机振动模型，通过设计相应的吸振器控制轧机辊系振动行为.

目前轧机液压系统非线性因素主要考虑了液压缸结构中缸体弹性刚度的非线性，液压缸在轧制过程的工作状态也较为复杂，除了缸体机械结构非线性因素影响外，缸体中液压油流量的非线性也可能对轧机振动造成影响，目前对该问题还尚未有较为深入的研究.

因此本文考虑液压压下缸中液压缸伺服电磁阀非线性流量影响，首先建立了非线性流量影响下液压缸的分段非线性刚度，然后建立了板带轧机液压压下-垂直振动的非线性动力学模型.通过求解得到该非线性系统的幅频方程，并研究了轧机振动系统在不同开折参数下的分岔行为.最后以实际轧机参数为例，研究了不同参数对轧机幅频特性的影响，为进一步抑制轧机振动提供理论参考.

1 非线性流量下分段非线性弹性力

在轧制过程中，板带轧机液压压下装置主要负责调节和保持辊缝大小，减小偶然冲击对轧件表面质量的影响，液压压下系统结构示意图如图1所示.

为了避免轧机在工作过程中产生的振动使液压缸无杆腔内压力油回流，一般在四通伺服电磁阀前设置单向阀来阻止此现象的发生.当活塞杆向下移动时单向阀打开，油液进入无杆腔；活塞杆向上移动时单向阀关闭，防止油液流出.此时油液流量将表现为不同的运动状态.液压压下缸的等效刚度可视作由活塞杆刚度与液压油刚度串联构成，在上述非线性流量的影响下液压压下缸的等效刚度可表示为

(1)

式中：k1(x)为活塞杆上移时无杆腔中液压油等效刚度；k2(x)为有杆腔等效刚度；k3(x)为活塞杆下移时无杆腔液压油等效刚度；e1为稳态轧制时液压压下缸无杆腔的初始压缩量.

图1 轧机液压系统示意图

在轧制过程中，液压缸中的液压油始终受力，产生形变.液压油体积变化情况可以通过体积模量表示，液压油受到的压力与体积变化关系可以表示为

(2)

式中：βe为液压油体积模量；ΔP为末态与初始压力差；ΔV为末态与初始体积差；V0为初始体积.

当活塞杆产生向下位移x时，无杆腔内的压力变化有

(3)

式中：L1为无杆腔长度；A1为无杆腔截面积.

采用如下的四通阀非线性流量方程[19]：

(4)

式中：Cd为四通阀的阀口流量系数；W为四通阀开口梯度；xv为四通阀阀芯位移；PL为液压缸无杆腔压力；Ps为伺服阀供油端压力；Pt为伺服阀回油端压力；ρ为液压油密度.

压下缸在工作过程中由缸内压力油提供轧制力，故阀芯位移xv≥0.考虑到管路压力损失对液压缸刚度影响较小[20]，可将其忽略，并考虑电磁阀响应速度，可得在响应时间内无杆腔内流入的液压油体积为

(5)

当活塞杆产生向下位移x时，可得由非线性流量产生的位移变化Δx与压力变化ΔPn：

(6)

(7)

由式(6)和式(7)可得无杆腔液压油等效刚度k3(x)为

(8)

将式(8)所表示的无杆腔等效刚度在稳态轧制处x=e1泰勒展开为

(9)

此时可得非线性弹性力

(10)

式(10)对位移进行积分，可得到非线性弹簧的弹性势能.由于弹性势能具有对称性，忽略奇次幂项，可得

(11)

对式(11)求导得到液压缸无杆腔位移x

(12)

同理有

(13)

此时可得到非线性流量影响下液压压下缸的分段非线性弹性力为

(14)

2 轧机辊系分段非线性动力学模型

四辊板带轧机主要由工作辊、支撑辊、轧机牌坊和传动装置等组成.由四辊轧机垂直方向结构可建立4自由度模型.其中工作辊与支撑辊关于轧件呈上下对称分布，在模型建立中可将4自由度的模型简化成2自由度的模型[21].并考虑上述非线性流量下液压压下缸的分段非线性作用力影响，可建立如图2所示的四辊轧机辊系分段非线性动力学模型.

图2 轧机辊系分段动力学模型

由于支承辊质量远大于工作辊质量，可将其简化为等效质量m[22]，k与c分别是轧件的等效线性刚度与阻尼，k2(x)为轧机液压压下缸有杆腔的等效刚度，k(x)为考虑非线性流量影响下液压压下缸分段非线性刚度，e1和e2分别为稳态轧制时液压压下缸内无杆腔与有杆腔内液体的初始弹性变形(e1>0，e2<0)，F(t)为周期性外扰力.

此时可得轧机辊系分段非线性动力学方程为

(15)

式中，x表示辊系垂直振动位移.将式(15)两侧同时除以m，可将方程化简为

(16)

此时的分段非线性弹簧力可写为

(17)

式中：γ11=k11/m；γ13=k13/m；γ21=k21/m；γ23=k23/m；γ33=k33/m.

3 轧机辊系的幅频特性方程

假设轧机辊系受到周期性外部扰动，该扰动为F(t)=Fsinωt.将系统中非线性影响因素考虑为弱非线性项，引入小参数ε，令ω2=(1+εσ)·ω02.式(16)可写为

(18)

(19)

式中，σ为频率调谐因子.

采用平均法，设式(18)有如下形式的解：

x=a(t)cosψ，ψ=ωt+φ(t).

(20)

式中：a(t)，φ(t)为时间t的慢变函数.

将式(20)代入式(18)中，可得出

(21)

(22)

(23)

其中，

式(23)即为板带轧机液压压下-垂直振动系统的幅频特性方程.

4 轧机辊系分岔特性分析

由于板带轧机模型参数或多或少会存在一定的误差，模型参数的轻微变化有可能使得轧机表现出不同的分岔特性.而分岔现象往往与轧机失稳状态密切相关，如果失稳状态不及时加以控制，将导致轧机辊系出现不可预测的振动，从而导致产品质量降低，严重时甚至会破坏轧机结构的稳定性.因此可采用奇异性理论研究轧机振动系统参数变化对轧机辊系振动状态的影响.

将式(23)在a=a0处进行泰勒展开，并略去高次项后可得

c1a6+c2a5+c3a4+c4a3+c5a2+c6a+c7=0.

(24)

4.1 分岔特性对比

在对轧机振动系统采用奇异性理论进行分析时，往往需要对原振动方程进行化简.采用不同的化简方式会对转迁集的分布产生较大影响.彭荣荣等[23]将余维为3的3-参数分岔方程，令其3个分岔参数分别为零，得到了3个二维转迁集，并对其分岔行为进行讨论.刘浩然等[16]令一个4-参数分岔方程的参数两两为零，得到了6个二维转迁集，并分别讨论其分岔行为.

本文通过将余维为4的4-参数分岔方程的参数分别为零，得到了4个三维转迁集，并在转迁集分隔成的空间内分别讨论其分岔图.避免了多个参数同时为零时结果与原分岔方程差异较大的情况发生.

4.2 分岔方程计算

令

(25)

将式(25)代入式(24)中进行线性变化，可得轧机液压压下与垂直振动系统的分岔方程：

(26)

理论与实践结合能有效促进工程建设发展。公路工程建设施工中，以理论指导实践非常重要。目前我国公路工程软基施工中存在缺乏理论指导的问题，应加强对软基处理技术的研究，促进工程技术的提升。施工中要遵循因地制宜的施工原则，不同的地区对工程施工技术有不同的要求。在施工前必须对周围环境进行细致考察。形成合理完善的施工方案。

α=(α1，α2，α3，α4)∈Rk;

由式(26)可得余维数为4的非奇异矩阵A，其中Gα，Gαz，…和gμ，gμz，…为G(z，μ，α)与g(z，μ)的一阶、二阶以及高阶偏导数.

(27)

可得出非奇异矩阵A行列式的值为

由奇异性理论可知，判定普适开折是否成立的充分必要条件是detA(0，0，0)≠0.因此G(z，μ，α)是g(z，μ)的一个4-参数普适开折，它与式(24)等价.μ是分岔参数，随着μ的变化，轧机辊系的分岔曲线将呈现出不同的拓扑结构，它的变化会导致辊系在某分岔形式下的振动幅度发生变化.由于该分岔方程所含参数较多，故分组讨论随着参数变化导致系统在三维空间上产生的分岔行为.系统转迁集由分岔集、滞后点集和双极限点集三者并集构成，可通过以下公式得出.

分岔集B={α∈Rk|，存在(x，μ)使得在(x，μ，α)处有G=Gα=Gμ=0}.

滞后点集H={α∈Rk|，存在(x，μ)使得在(x，μ，α)处有G=Gx=Gxx=0}.

双极限点集D={α∈Rk|，存在(xi，μ)(i=1，2)，x1≠x2，使得在(xi，μ，α)处有G=Gx=0}.

分岔集为

BR=∅，BI=∅.

其中下角标R和I表示常规转迁集和由于约束条件所产生的新增转迁集.

滞后点集为

HR={30m4+6α2m+2α3=0}，HI=∅.

双极限点集为

DR=∅，DI=∅.

系统的转迁集为

J=HR∪DR∪BI∪HI∪DI.

此时系统的转迁集图像如图3所示.由转迁集分隔出的各个子空间内的分岔图形如图4所示.分岔图中横坐标代表分岔参数，纵坐标代表系统振幅.

图3 α1=0时轧机振动系统的转迁集

图4 α1=0时图3中不同区域内分岔图

由图3与图4可知，转迁集将轧辊系统的振动情况分隔成7个空间，在不同的区域内轧辊振动的分岔图曲线拓扑性质不同.因为轧辊的振幅总是一个大于零的值，即a≥0，在分岔方程中的约束条件为z≥β，β=h2/7.当轧机振动系统的参数在转迁集所分隔的区域内时，分岔图的拓扑性质不发生变化.

2) 当α2=0时.

分岔集为

BR=∅，BI=∅.

滞后点集为

HR={6m5+4α1m3+2α3m+α4=0}，HI=∅.

双极限点集为

DR=∅，DI=∅.

此时系统转迁集如图5所示，分岔图形如图6所示.

图5 α2=0时轧机振动系统的转迁集

图6 α2=0时图5中不同区域内分岔图

3) 当α3=0时.

分岔集为

BR=∅，BI=∅.

滞后点集为

双极限点集为

DR=∅，DI=∅.

此时轧机振动系统转迁集如图7所示，分岔图形如图8所示.

4) 当α4=0时.

分岔集为

BR=∅，BI=∅.

滞后点集为

HR={24m3+8α1m+3α2=0}，HI=∅.

图7 α3=0时轧机振动系统的转迁集

图8 α3=0时图7中不同区域内分岔图

双极限点集为

DR=∅，DI=∅.

此时轧机振动系统转迁集如图9所示，分岔图形如图10所示.

图9 α4=0时轧机振动系统的转迁集

图10 α4=0时图9中不同区域内分岔图

通过上述四种情况可以看出，转迁集将轧机振动系统分成若干个子空间，在同一子空间内任意选择点振动系统的分岔图均有着相同的拓扑性质，在不同的子空间内取点轧机振动系统分岔图的拓扑性质不同.当轧机振动系统不同参数发生变化时会出现不同的分岔行为，通过研究分岔行为可以避免轧机振动系统进入不稳定的参数区域.同时可以得出，改变分岔参数μ可使轧机辊系处于不同的振动状态，通过调整分岔参数可让轧辊的振动处于稳定状态，这能够为降低轧机辊系的垂直振动提供理论参考.

5 仿真研究

根据某厂1780热连轧机图纸数据进行仿真，取轧机的实际结构参数以及液压压下系统参数如表1所示.

表1 1780热连轧机参数

由式(23)中的幅频方程可以得到如图11所示的板带轧机液压压下-垂直振动系统的幅频曲线，从图中可以看出，当外扰力频率增加时，轧机辊系的振幅出现了非线性系统中特有的“跳跃”与“滞后”现象，当外扰频率从曲线左侧开始增大时，振幅沿曲线从点A经过点B，C移动至点D，达到最大振幅，随着外扰力频率继续增加超过点D时，辊系的振幅突然降至点F，随后经过点G，轧辊的振幅在点D出现了跳跃现象；当外扰频率从右侧开始减小时，轧辊振幅沿曲线从点G经过点F到达点E，当频率继续降低小于点E时，轧辊振幅升至点C并沿着CB方向运动经过点A.因此当外扰力频率经过DE段时，系统的振幅是不稳定的，在实际工作中需要避免.

图11 振动系统幅频响应曲线

同时可以发现，由于液压压下系统的分段非线性弹簧力，当外扰频率到点B时，振动系统的非线性刚度发生改变，幅频曲线出现拐弯.这种现象说明，在不同的振幅区段，由于振动系统非线性刚度变化会导致轧辊处于不同的振动状态.

图12为伺服电磁阀不同响应时间下的系统幅频曲线.从图中可以看出，随着电磁阀响应时间变小，幅频曲线的拐弯特性变弱，在拐弯点附近的非线性刚度变化变小，轧机辊系的不稳定振幅频率区域变小；当电磁阀的响应时间变长，轧辊的不稳定振幅频率区域变大，幅频曲线的滞后现象明显.

图12 响应时间对幅频特性的影响

图13为不同外激励下振动系统幅频曲线的变化情况.从图中可以看出，随着外激励增加，在同一扰动频率下的轧辊振幅将变大，共振频率范围也将变宽.当外激励使振幅低于e1时，幅频曲线中的拐弯特性消失，轧辊处于稳定状态.

图13 外激励对系统幅频特性的影响

图14为不同阻尼比下轧机振动系统的幅频曲线.可以看出，随着阻尼比的增加，轧辊的振动幅值降低，当阻尼比增加到一定程度时，轧机振动系统的幅频曲线拐弯现象消失，因此可通过增加轧机系统的阻尼比降低轧辊振动.

图14 阻尼比对系统幅频特性的影响

由图11～图14可以看出，伺服阀响应速度、外激励与阻尼比的变化均会对振动系统的幅频曲线产生影响，通过改变参数(如提高系统的阻尼比和降低外激励幅值，提高伺服电磁阀响应速度等)可以使轧机辊系的振动减弱，避免在工作中出现共振现象.

6 结 论

1) 依据四辊轧机辊系机械的特点和液压系统非线性特性，考虑板带轧机液压系统电磁伺服阀非线性流量对垂直振动的影响，推导出液压系统非线性流量与液压压下缸非线性弹簧力的关系，建立了板带轧机液压压下-垂直振动动力学模型，为研究非线性因素对轧机辊系振动影响提供了一种新的模型依据.

2) 分析了轧机振动系统在不同参数下的分岔拓扑结构，得到4组不同的转迁集与分岔图.可以看出，随着分岔参数m取值的变化，分岔图呈现出不同的拓扑结构，导致轧机辊系处于不同的振动状态.因此可通过调整工艺参数以及分岔参数m避免轧机处于发生振动的频率区.

3) 研究系统参数对板带轧机液压压下-垂直振动系统幅频特性中的拐弯与跳跃现象的影响，提高伺服电磁阀的响应速度可以有效减小非线性系统的跳跃现象，减小振幅不稳定频率区域.并且通过提高系统阻尼和降低外激励幅值可以有效降低系统振幅，避免轧机出现共振现象，从而可以提高轧机辊系工作的稳定性.