基于動態軋制力的四輥冷軋機七自由度振動特性研究

基于動態軋制力的四輥冷軋機七自由度振動特性研究

劉浩然1，張瑩2，時培明3，趙紅旭3

(1. 河北省 特種光纖與光纖傳感實驗室，秦皇島066004;2.天津 一重電氣自動化公司，天津300457； 3.燕山大學電氣工程學院，秦皇島066004)

摘要：針對四輥冷軋機座輥系系統，考慮軋輥振動位移對軋制力影響，建立基于動態軋制力的非線性七自由度振動模型；用平均法求解該模型獲得系統的幅頻特性方程，利用奇異性理論獲得系統的分岔行為。通過現場數據分析各參數對幅頻特性及對系統分岔、混沌特性影響?？蔀樗妮伬滠垯C的動力學研究、抑制振動提供理論參考。

關鍵詞：四輥冷軋機；動態軋制力；平均法；分岔；混沌

中圖分類號:TH113.1;TB123文獻標志碼：A

基金項目：國家973基金資助項目(2011CB0132005)

收稿日期：2014-11-12修改稿收到日期：2015-05-20

收稿日期：2014-07-22修改稿收到日期：2014-12-18

Vibration characteristics of a cold rolling mill with a model of seven degrees of freedom and their influence on dynamic rolling force

LIUHao-ran1,ZHANGYing2,SHIPei-ming3,ZHAOHong-xu3(1. College of Information Science and Engineering, Yanshan University, Qinhuangdao 066004, China；2. Tianjin first heavy Electric Automation Co.,Ltd, Tianjin 300457, China；3. Institute of Electrical Engineering, Yanshan University, Qinhuangdao 066004, China)

Abstract:Considering the influence of roll vibration on the rolling force, a nonlinear seven degrees of freedom vibration model was established for the stand rolls system of a four-roll cold mill. The amplitude-frequency characteristic equation of the nonlinear multiple degrees of freedom model was derived by using the average method, further, the bifurcation behavior of the system was revealed in light of the singularity theory. The influence of various parameters on the amplitude-frequency characteristic, and the bifurcation and chaos characteristics of the system were analysed with practical data. The results provide a theoretical reference for studying the kinematics and inhibiting the occurrence of vibration of four-roll cold tandem mills.

Key words:four-roll cold mill; dynamic rolling force; average method; bifurcation; chaos

軋機座系統(尤其高速冷連軋機)在軋制過程中，輥系普遍存在的垂直振動會造成板帶厚度波動及表面振痕，嚴重影響板帶質量，降低生產效率，并威脅設備的生產安全，為亟需解決之問題[1-2]。

Younes等[3]據實驗所得工藝參數建立軋機座線性垂直振動模型，且將軋件等效成具有線性剛度的彈性元件。馬維金等[4]提出六自由度集中質量模型，利用MATLAB軟件計算模型的固有頻率，并與現場測試結果對比驗證模型的正確性。楊旭等[5]通過建立兩自由度垂向系統非線性自激振動模型，分析非線性參數變化對系統穩定性及振動特性影響，并利用數值方法驗證理論推導結果的正確性。侯東曉等[6]通過分析非線性剛度、阻尼等參數對系統振動影響知，非線性剛度變化會引起激勵幅值跳躍。因冷連軋機軋制過程的復雜性，多自由度模型更能正確反映軋機的工作狀態，因此本文考慮軋輥振動位移對軋制力影響，建立七自由度非線性垂直振動模型；并用平均法對該系統解析求解；據某廠實際參數，通過仿真驗證分析系統的幅頻特性、分岔及混沌特性，為抑制軋機輥系間發生振動提供理論依據。

1軋機座輥系七自由度非線性垂直振動模型

1.1垂直振動力學模型

對軋機垂直系統一般采用集中參數模型，即將垂直模型視為若干慣性元件與彈性元件組成的質量彈簧系統。本文將四輥軋機垂直系統簡化為七自由度彈簧質量系統計算模型，見圖1。其中，m1～m7分別為上機架、液壓缸、上支撐輥系、上工作輥系、下工作輥系、下支承輥系及下機架等效質量；k1，c1分別為上、下機架間等效剛度及阻尼；k2，c2分別為上機架與液壓缸間等效剛度及阻尼；k3，c3分別為液壓缸與上支承輥間等效剛度及阻尼；k4，c4分別為上支承輥與上工作輥間等效剛度及阻尼；k5，c5分別為上、下工作輥間等效剛度及阻尼；k6，c6分別為下工作輥與下支承輥間等效剛度及阻尼；k7，c7分別為下支承輥與下機架間等效剛度及阻尼；k8，c8分別為下機架與地面間等效剛度及阻尼；ΔP為軋制力變化量。

1.2垂直振動數學模型

據垂直振動動力學模型得軋機機械結構振動方程為

(1)

考慮軋機工作過程中上下輥振動位移、速度時刻變化，因此本文采用動態軋制力模型[7]，即

P(x)=P(x0)+ΔP(Δx)

(2)

式中：P(x0)為未發生振動時軋制力；ΔP(Δx)為動態軋制力變化量，Δx=Δx4-Δx5為上下工作輥位移變化量之差。

ΔP(Δx)=b1Δx+b2Δx2+b3Δx3

(3)

以上工作輥為研究對象，考慮動態軋制力下用軋制力變化量代替上下工作輥間等效剛度及阻尼，據圖1模型得上工作輥的振動方程為

P(x0)+ΔP(Δx)

(4)

令x3=x30+Δx3，x4=x40+Δx4，x5=x50+Δx5，且x30、x40、x50為穩態位移，Δx3、Δx4、Δx5為位移變化量，則式(4)轉化為

k4(x30+Δx3-x40-Δx4)=

P(x40-x50)+ΔP(Δx4-Δx5)

(5)

P(x40-x50)=-k4(x30-x40)

(6)

將式(3)、(6)代入式(5)，得關于位移變化量的非線性振動方程為

k4(Δx3-Δx4)-b1(Δx4-Δx5)-

b2(Δx4-Δx5)2-b3(Δx4-Δx5)3=0

(7)

為便于表達，令x3=Δx3，x4=Δx4，x5=Δx5，x=Δx4-Δx5，則式(7)可轉化為

b1x-b2x2-b3x3=0

(8)

2七自由度非線性振動方程振動響應分析

2.1平均法求解系統非線性振動響應方程

冷連軋機軋制過程中機座的垂直輥系系統會受前后機座及周圍環境影響，從而影響板帶表面質量。設軋機受周期性外部擾動Q，則式(1)可寫為

(9)

式中：Q為外部擾動矩陣。

當動態軋制力中非線性項為高階小量時，可認為式(9)的解具有形式為

(10)

由平均法可解出

(11)

式中：

(12)

設上工作輥受外部擾動為q4cos(ωt)，由式(12)得

f4(aj,θj,φj)=

K45a5cosφ5)-(C43a3sinφ3+C44a4sinφ4+

b2(a4cosφ4-a5cosφ5)2+b3(a4cosφ4-

(13)

則(aj,θj)的微分方程為

(14)

由平均法定義，可導出系統的平均化方程為

(15)

式中：

b2(a4cosφ4-a5cosφ5)2sinφ4+b3(a4cosφ4-

a5cosφ5)cosφ4+b2(a4cosφ4-a5cosφ5)2cosφ4+

(17)

系統(9)與激勵同頻的解為

xsj=asjcos(ωt-θsj),(j=1,2,…,7)

(18)

對應平均化方程(15)的非零平衡點(asj,θsj)，滿足超越方程為

U4(asj,θsj)=0，S4(asj,θsj)=0

(19)

由式(19)導出系統的幅頻特性曲線為

(20)

2.2系統非線性幅頻特性分析

以某四輥軋機實際參數為例，機架輥系各部等效質量見表1，等效剛度見表2，等效阻尼見表3。其中，q4=1.5×105N；b1=-2.291 8×1010N/m；b2=-9.769 5×1012N/m；b3=-9.345 2×1016N/m。

表1　機架輥系各部分等效質量(kg)

表2　機架輥系各部分等效剛度(GN/m)

表3　機架輥系各部分等效阻尼(MN·S/m)

因動態軋制力代替上下工作輥間等效阻尼及剛度，故K44=k4，C44=c4。c4，k4取不同數值時，由仿真獲得上工作輥振動方程的幅頻特性曲線，見圖2、圖3。由圖2看出，阻尼因子不僅影響振幅，亦影響共振區域，且阻尼因子增大振幅減小，共振區域亦減小。由圖3看出，剛度主要影響共振振幅，即剛度增加系統共振振幅逐漸增大。

圖2　不同阻尼下幅頻曲線 Fig.2 The amplitude-frequency curves under different damping

圖3　不同剛度下幅頻曲線 Fig.3 The amplitude-frequency curves under different stiffness

3系統非線性分岔特性分析

式(20)兩邊同除以J1，得系統分岔響應方程為

(21)

式中：l=-(J2-J3ω2)/J1；h=(J4ω4+J5ω2+J6)/J1；

g=(J7ω2+J8)/J1。

將式(21)兩邊同乘以a4，得

(22)

據奇異性理論，此時系統的轉遷集為

(1)分岔點集：B0(Z2)=B1(Z2)=φ(φ為空集)

(2)滯后點集：H0(Z2)={h=0}，H1(Z2)={h=l2/3,l≤0}

(3)雙極限點集：D(Z2)={h=l2/4,l≤0}

(4)轉遷集：∑=B0∪B1∪H0∪H1∪D

系統轉遷集將平面分成四個區域(1)、(2)、(3)、(4) 見圖4，轉遷集臨界點與各區域對應的分岔拓撲見圖5。圖5中(1)反映系統在分岔點處于穩定狀態，而(2)、(3)、(4)則表示在分岔點出現滯后、跳躍現象，造成軋機產生振動。因此應將開折參數選在系統的穩定區域。

圖4　轉遷集 Fig.4 Transition

圖5　分岔拓撲結構 Fig.5 Bifurcation topology

4七自由度非線性系統混沌特性分析

冷連軋機軋制工藝的復雜性決定整個軋制過程會呈現不同的運動狀態。軋機輥系系統全局分岔見圖6、圖7。由兩圖看出，系統出現周期、陣發性混沌及倍周期分岔等多種不同運動狀態。外擾力取不同值時相軌跡及Poincare截面見圖8～圖12。

圖6　外擾力q 4變化時分岔圖 Fig.6 Bifurcation with changing of external incentive force q 4

圖7中，當外擾力q4<3.3×105時，Lyapunov指數位于零位以下，系統處于穩定的周期性運動狀態。圖8中，q4=1×105時相軌跡由三條閉合曲線組成，Poincare截面由三個孤立點組成，反映出此時系統處于周期運動狀態。

圖7　Lyapunov指數 Fig.7 Lyapunovexponent

圖8　q 4=1×10 5時相軌跡與Poincare截面 Fig.8 The phase trajectories and Poincare with q 4=1×10 5

圖9　q 4=3.3×10 5時相軌跡與Poincare截面 Fig.9 The phase trajectories and Poincarewith q 4=3.3×10 5

圖10　q 4=5×10 5時相軌跡與Poincare截面 Fig.10 The phase trajectories and Poincarewith q 4=5×10 5

圖11　q 4=6.5×10 5時相軌跡與Poincare截面 Fig.11 The phase trajectories and Poincarewith q 4=6.5×10 5

圖12　q 4=8×10 5時相軌跡與Poincare截面 Fig.12 The phase trajectories and Poincare with q 4=8×10 5

外擾力q4處于3.3×105～4×105區間時系統發生陣發性混沌，由圖7 Lyapunov指數曲線在零位以上、以下交替出現，即系統運動狀態在混沌運動與周期運動狀態之間切換。由圖9看出，q4=3.3×105時相軌跡與Poincare截面的圖像表述均滿足系統陣發性混沌特征。外擾力q4處于4×105～6×105區間時系統處于穩定的周期運動狀態。由圖10相軌跡與Poincare看出，系統相軌跡由八條閉合曲線組成，Poincare截面由八個孤立點組成，此時系統處于周期8運動狀態。當q4=6.5×105時系統由陣發性混沌退化為周期3運動。由圖11看出，相軌跡反復纏繞，Poincare截面出現混沌吸引子，且圖7Lyapunov指數位于零位以上，證明系統此時處于混沌運動狀態。當q4=8×105時系統再次發生陣發性混沌，由圖12看出，Poincare截面出現混沌吸引子，Lyapunov指數大于零，滿足系統混沌狀態特征。

5結論

(1)考慮工作輥在軋制過程中振動因素對軋制力影響，建立基于動態軋制力的七自由度非線性振動模型。

(2)用平均法求解七自由度非線性振動方程獲得上工作輥的幅頻特性方程。分析幅頻曲線知，阻尼因子增大共振振幅、區域逐漸減??；剛度增大共振振幅增加。

(3)利用奇異性理論討論系統分岔特性，確定使系統穩定運行的開折區域。分析分岔圖、Lyapunov指數、相軌跡及Poincare截面知，外擾力變化會使系統呈現不同運動形式。

參考文獻

[1]瞿志豪,王康健,魏立群,等．冷連軋機顫振的一種振因研究[J]．機械設計與研究,2006,22(1):89-92.

QU Zhai-hao,WANG Kang-jian, WEI Li-qun,et al. A study of the chatter cause on the cold tandem rolling mill[J]. Machine Design and Research, 2006, 22(1): 89-92.

[2]閆曉強.連軋機振動控制重要進展綜述[J]. 冶金設備, 2013(2): 55-57.

YAN Xiao-qiang. Important progress on tandem rolling mill vibration control[J]. Metallurgical Equipment, 2013(2): 55-57.

[3]Younes M A, Shahtout M, Damir M N. A parameters design approach to improve product quality and equipment performance in hot rolling[J]. Journal of Materials Processing Technology, 2006, 171(1): 83-92.

[4]馬維金,李鳳蘭,王俊元,等．四輥軋機的六自由度垂直振動模型研究[J].中國機械工程,2011,22(24):2962-2965.

MA Wei-jin, LI Feng-lan, WANG Jun-yuan,et al. Modelling and validation of vertical vibration with 6 degrees of freedom for 4-high hot strip mill stand[J]. China Mechanical Engineering, 2011, 22(24): 2962-2965.

[5]楊旭,李江昀,童朝南. 冷軋機垂向輥系非線性振動建模與穩定性分析[J]．振動、測試與診斷,2013,33(2): 303-307.

YANG Xu, LI Jiang-yun, TONG Chao-nan. Nonlinear vibration modeling and stability analysis of vertical roller system in cold rolling mill[J]. Journal of Vibration, Measurement & Diagnosis, 2013, 33(2): 303-307.

[6]侯東曉,陳浩,劉斌,等軋機輥系垂直非線性參激振動特性分析[J]．振動與沖擊,2009,28(11):1-6.

HOU Dong-xiao, CHEN Hao, LIU Bin,et al. Analysis on param etrically excited nonlinear vertical vibration of roller system in rolling mills[J]. Journal of Vibration and Shock, 2009, 28(11): 1-6.

[7]侯東曉,朱月,劉浩然,等. 基于動態軋制力的冷軋機非線性振動特性研究[J]．機械工程學報,2013,49(14):45-50.

HOU Dong-xiao, ZHU Yue, LIU Hao-ran,et al. Research on nonlinear vibration characteristics of cold rolling mill based on dynamic rolling force[J]. Journal of Mechanical Engineering, 2013, 49(14): 45-50.

第一作者董得義男，博士，副研究員，1979年生

通信作者李志來男，研究員，1965年生

第一作者張大海男，博士生，1988年生

通信作者費慶國男，教授，1977年生