高崇一， 杜國君， 李建雄， 胡發(fā)科

(1. 燕山大學(xué) 河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點(diǎn)實(shí)驗室，河北　秦皇島　066004;2.燕山大學(xué) 電氣工程學(xué)院自動化系，河北　秦皇島　066004)

?

慣性邊界下帶鋼的非線性振動分析

高崇一1， 杜國君1， 李建雄2， 胡發(fā)科1

(1. 燕山大學(xué) 河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點(diǎn)實(shí)驗室，河北秦皇島066004;2.燕山大學(xué) 電氣工程學(xué)院自動化系，河北秦皇島066004)

摘要：根據(jù)連軋機(jī)軋制過程中帶鋼與軋輥的運(yùn)動機(jī)理，將相鄰兩機(jī)架間的帶鋼簡化為軸向運(yùn)動的Euler梁，軋輥簡化為定軸轉(zhuǎn)動的慣性元件，建立Euler梁在慣性邊界下的非線性振動力學(xué)模型。基于哈密頓原理建立軸向運(yùn)動Euler梁的縱向和橫向非線性振動微分方程，利用Kantorovich時間平均法簡化運(yùn)動方程和邊界條件，并采用修正迭代法求解運(yùn)動方程。最后通過數(shù)值計算獲得了Euler梁非線性振動的幅頻響應(yīng)曲線，并討論慣性邊界條件下的軸向運(yùn)動速度、長度和軋輥的轉(zhuǎn)動慣量對Euler梁振動特性的影響。研究結(jié)果可為控制和分析連軋過程中帶鋼垂直振動提供重要的理論參考。

關(guān)鍵詞：Euler梁；慣性邊界；帶鋼；垂直振動；修正迭代法

軋機(jī)振動普遍存在于軋制生產(chǎn)中，復(fù)雜的振動現(xiàn)象是導(dǎo)致帶鋼表面缺陷的重要因素，并且影響設(shè)備關(guān)鍵部件的使用壽命。關(guān)于軋機(jī)振動問題的研究涉及軋機(jī)系統(tǒng)各個部件的不同振動形式，如軋輥的垂直振動、水平振動、軸向串動、帶鋼的橫向和縱向振動、主傳動系統(tǒng)的扭轉(zhuǎn)振動和軸向振動等[1]。其中帶鋼的橫向(垂直)和縱向振動方面的研究相對較少，而連軋中帶鋼的振動又是軋制過程中不可避免的。因此，對帶鋼振動的研究成了鋼鐵行業(yè)亟待解決的重要問題。

軋制過程中，帶鋼的垂直振動對軋機(jī)和產(chǎn)品質(zhì)量的影響相對較大。若忽略軋機(jī)輥系的垂直振動，可將帶鋼簡化為軸向運(yùn)動梁[2-4]，帶鋼垂直振動可等效為運(yùn)動梁的橫向振動，從而可以利用運(yùn)動梁理論分析帶鋼振動。劉明哲等[3]建立了運(yùn)動板帶的二維動力學(xué)模型和三維動力學(xué)模型，并利用數(shù)值方法分析了軋制過程中運(yùn)動帶鋼的固有特性和穩(wěn)定性。Sun等[4]考慮兩機(jī)架間帶有時變張力的運(yùn)動帶鋼的非線性振動模型，并采用多尺度法求解。Bor等[5]在冷連軋過程中分析了含有非線性參數(shù)的自激振動現(xiàn)象，并給出了近似解析解。Miranker[6]首次推導(dǎo)出了運(yùn)動弦線橫向振動的運(yùn)動微分方程，并對Euler梁進(jìn)行了相應(yīng)分析。Mote[7]分析了軸向變速運(yùn)動系統(tǒng)的橫向振動穩(wěn)定性問題。陳樹輝等[8]對軸向運(yùn)動Euler梁非線性振動內(nèi)共振進(jìn)行了研究。Suweken等[9-10]研究了考慮兩端邊界條件下的軸向運(yùn)動Euler梁的振動特性。

帶鋼振動的研究大多是基于簡支、固支或帶有扭轉(zhuǎn)彈簧的邊界條件[11]，而忽略了軋輥轉(zhuǎn)動慣量的影響，與軋制過程的實(shí)際偏離較大。本文考慮軋輥的轉(zhuǎn)動慣量對帶鋼垂直振動的影響，并且基于軸向運(yùn)動梁理論[12-13]，在不考慮剪切變形的情況下，可將帶鋼振動等效為Euler梁振動，建立軋輥與Euler梁相互作用的動力學(xué)模型，求解慣性邊界條件下軸向運(yùn)動梁的非線性振動。本文采用修正迭代法[14-15]求解軸向運(yùn)動Euler梁振動方程，并利用Matlab對系統(tǒng)進(jìn)行數(shù)值仿真，研究系統(tǒng)各參數(shù)對帶鋼非線性振動特性的影響，為帶鋼連軋生產(chǎn)中的優(yōu)化控制提供理論參考。

1模型的建立

圖1　機(jī)架間運(yùn)動帶鋼及其等效的Euler梁力學(xué)模型Fig.1 Moving strip between stands and mechanical model of Euler beam

根據(jù)哈密頓原理，對軋制過程中的力學(xué)模型進(jìn)行分析，建立數(shù)學(xué)模型。軸向運(yùn)動Euler梁的動能T1為：

(1)

式中:ρ為帶鋼密度，A為Euler梁的橫截面積，u,t與u,x0分別為函數(shù)u(x0,y0,t)對變量t與x0的一階偏導(dǎo)。后文中u,x0x0表示u(x0,y0,t)對x0的二階偏導(dǎo)，依次類推，其他類似符號表示含義類同。

Euler梁的形變勢能U為：

(2)

式中:E為彈性模量，I為慣性矩。

軋輥的動能T2為：

(3)

式中:r為軋輥半徑。

由哈密頓方程：

(4)

可得出慣性邊界下Euler梁的運(yùn)動方程：

EAw,x0w,x0x0=0

(5)

(6)

邊界條件為：

x0=0或x0=l時，

(7)

w(0)=w(l)=w,x0x0(0)=

w,x0x0(l)=0;wmax=φm

(8)

多數(shù)情況下，帶鋼橫向振動引起的軸向運(yùn)動動能相對較小，可令u,t=u,tt=0，則式(5)可簡化為：

(9)

若忽略帶鋼軸向位移對垂直振動的影響，得出的運(yùn)動方程與文獻(xiàn)[4]一致。本文數(shù)學(xué)建模時，將軸向位移與軋輥慣性邊界條件對帶鋼垂直振動的影響加以考慮，所得運(yùn)動方程更加精確。

(10)

(11)

邊界條件整理為：

x0=0或x0=l時，

(12)

(13)

為使以下運(yùn)算簡便，引入無量綱量

將上式代入運(yùn)動方程(10)、(11)和邊界條件(12)、(13)，進(jìn)行無量綱化，則Euler梁的運(yùn)動方程整理為：

(v02-1)φ,xx-φ,xφ,xx=0

(14)

φ,xxxx+Sv02φ,xx-ω2φ-

(15)

邊界條件為：

x=0和x=1，

(16)

φ(0)=φ(1)=φ,x(0)=

(17)

2非線性方程的解

由于求解方程(14)～(17)過程較為復(fù)雜，因此，本節(jié)將采用修正迭代法進(jìn)行求解。

2.1一階近似解

首先，將方程(15)中所有非線性項略去后，可寫為

φ1,xxxx-ω2φ1=0

(18)

求出方程(18)的級數(shù)解：

φ1(x)=a0M0(x)+a1N0(x)+

a2I0(x)+a3K0(x)

(19)

式中:

然后，將φ1(x)代入邊界條件 (17)可解得近似的一階頻率ω1=16.71，可求出系數(shù)

a0=0，a1=0，a2=μ1φm，a3=μ2φm

其中:

從而得出

(20)

接下來將φ1(x)代入到方程(14)中得出：

(21)

由邊界條件(16)得：

2.2二階修正迭代解

下面進(jìn)行二階迭代，將求得的φ1(x)和φ1(x)代入方程(15)并簡化后得：

φ2,xxxx-ω2φ2=αφ1,xx+βφ1,x2φ1,xx

(22)

式中:

其中:

(n=2,3…)

(n=2,3…)

(n=3,4…)

由級數(shù)解的性質(zhì)，方程(22)的解可寫成如下形式：

其中:A0=B0=C0=C1=D0=E0=F0=0

ζ1和ζ2為待定系數(shù)。

將φ2(x)代入邊界條件 (17)，可得：

Dξ=0

(24)

其中:

由detD=0可求出振頻ω2的解析表達(dá)式。由式(24)可求出系數(shù)ζ1和ζ2，從而可求得振幅φ2的解析表達(dá)式，二階修正迭代解得以確定。

3算例與討論

對文獻(xiàn)[16]中的軋機(jī)系統(tǒng)進(jìn)行數(shù)值模擬，帶鋼厚度h為2 mm；軋制速度為5 m/s，則Euler梁無量綱軸向速度v0為1×10-3；軋輥轉(zhuǎn)動慣量J無量綱量為K，當(dāng)K=5，且長度l分別為1.0 m、1.5 m、2.0 m、2.5 m、3.0 m、3.5 m時，考慮振幅φm對振頻ω的影響，其幅頻響應(yīng)曲線如圖2所示，圖中φm為帶鋼振幅無量綱量φ的最大值。在軋輥轉(zhuǎn)動慣量不變的情況下，Euler梁的長度對幅頻響應(yīng)特性影響較大，它將改變振動性態(tài)。從圖2中可以看出，當(dāng)l=1.0～3.5 m時，Euler梁的振動性態(tài)逐漸由硬化型過渡到軟化型，當(dāng)梁長l取值為1.0 m、1.5 m、2.0 m時振動性態(tài)表現(xiàn)為硬化型，即隨著振幅增大頻率逐漸增大，且隨著梁長增大，硬化程度逐漸降低；當(dāng)梁長l取值為2.5 m、3.0 m、3.5 m時，曲線表現(xiàn)為遞減趨勢，非線性表現(xiàn)為軟化型，即隨著振幅增大頻率逐漸減小，并呈現(xiàn)出軟化程度增強(qiáng)的趨勢。隨著長度的遞增，幅頻響應(yīng)曲線振動性態(tài)呈現(xiàn)出較均勻的變化，其中硬化型與軟化型的分界在梁長l為2.0～2.5 m區(qū)間?；谏鲜龇治?，可以得出，在軋制速度和軋輥轉(zhuǎn)動慣量一定的情況下，機(jī)架間距對帶鋼振動的幅頻特性具有較大影響。根據(jù)不同機(jī)架間距下帶鋼振動幅頻響應(yīng)曲線，可對相關(guān)參數(shù)進(jìn)行適當(dāng)調(diào)節(jié)，以減小帶鋼振動對產(chǎn)品質(zhì)量的影響。

圖3為不同振幅下K值與頻率的關(guān)系曲線，當(dāng)Euler梁的長度l=2.0 m時，軸向速度v0=1×10-3，且振幅分別為0.05、0.10、0.15、0.20、0.25時，考慮K值對振頻的影響，其中K值變化體現(xiàn)出軋輥轉(zhuǎn)動慣量的變化。隨著K值的增大，Euler梁的振頻逐漸增大，K值在由0增大到5的過程中，曲線呈現(xiàn)出大幅度遞增趨勢，說明此區(qū)間的K值對振頻的影響較大，并且每條曲線在K=1.5附近近似交于一點(diǎn)，此時對應(yīng)的振頻(ω=16.71)與基頻一致。而隨著K值繼續(xù)遞增，曲線又趨于平緩。隨著K值的無限增大，Euler梁的邊界條件無限趨近于固支狀態(tài)，振頻趨近于固支條件下的頻率，此時近似為軋輥不轉(zhuǎn)動，將帶鋼兩側(cè)軋輥視為固定支架的情況。另外，由圖3可以看出，振幅φm越小，隨K值變化曲線越早趨于平緩，振幅φm越大，K值對頻率的影響范圍越大，影響越明顯。

圖4為不同K值(軋輥轉(zhuǎn)動慣量)下的速度與頻率關(guān)系曲線，當(dāng)梁長l=2.0 m，振幅φm=0.15，且K值分別為5、10、15、20、25時，考慮速度v0對振頻的影響。從圖中可以看出，不同K值下的振動頻率隨軸向速度增大而減小，并且，隨著v0的不斷增大，振頻減小率呈現(xiàn)出增大的趨勢，即軋制速度越大，速度變化對帶鋼振動頻率的影響越大，且頻率減小的越快。當(dāng)K值為15、20、25時，v0-ω2關(guān)系曲線近乎重合，可以看出軋輥轉(zhuǎn)動慣量無限大的情況下，邊界條件趨近于固支狀態(tài)。

圖2 不同梁長下的幅頻響應(yīng)曲線Fig.2Theamplitudefrequencyresponsecurvesunderthedifferentlengthsofthebeam圖3 不同振幅下K值與頻率關(guān)系曲線Fig.3Kandfrequencycurvesunderthedifferentamplitudes圖4 不同K值下速度頻率關(guān)系曲線Fig.4VelocityandfrequencycurvesunderthedifferentK

通過圖3和圖4還可以看出，軋輥的轉(zhuǎn)動慣量對帶鋼振動的影響尤為突出，在轉(zhuǎn)動慣量較小時，對帶鋼振動影響較大，研究帶鋼振動時需要考慮軋輥轉(zhuǎn)動慣量對其影響；在轉(zhuǎn)動慣量較大時，可以近似忽略其對帶鋼振動的影響。

圖5　不同速度下幅頻響應(yīng)曲線Fig.5 The amplitude frequency response curves under the different velocities

圖5為不同速度下的幅頻響應(yīng)曲線，當(dāng)Euler梁的長度l=2.0 m，K=5，軸向速度v0分別為0×10-3、1×10-3、2×10-3、3×10-3、4×10-3時，振幅φm對振頻的影響。從圖中可以看出，隨著振幅φm的增大，Euler梁的振頻逐漸增大。當(dāng)軸向速度v0為0或1×10-3時，兩條曲線幾乎重合，可以看出當(dāng)速度很小時，振幅φm對頻率的影響較大；而隨著速度v0增大，曲線漸漸趨于平緩，說明軋制速度越大，帶鋼振動的振幅φm對頻率的影響越小。

4結(jié)論

模擬軋制過程中軋輥與帶鋼的運(yùn)動，將其近似為慣性邊界條件下的Euler梁，并建立其力學(xué)模型和數(shù)學(xué)模型，采用修正迭代法求解方程，并通過數(shù)值計算方法，討論了相鄰兩機(jī)架間帶鋼的振幅、長度、軋制速度和軋輥的轉(zhuǎn)動慣量對帶鋼振頻的影響。

(1) 在軋制速度和軋輥轉(zhuǎn)動慣量一定的情況下，機(jī)架間距大小對帶鋼振動性態(tài)的影響。當(dāng)機(jī)架間距較小時，即l≤2.0 m，隨著振幅的增大，振頻響應(yīng)逐漸增強(qiáng)，非線性表現(xiàn)為硬化型，隨著機(jī)架間距由小逐漸變大，帶鋼振動性態(tài)的非線性硬化程度逐漸減弱。隨著機(jī)架間距繼續(xù)增大，幅頻曲線呈現(xiàn)出振頻隨著振幅的增大而逐漸減小的趨勢，即振動性態(tài)表現(xiàn)為軟化型，并且，隨著機(jī)架間距的增大，軟化程度逐漸加強(qiáng)。

(2) 在機(jī)架間距一定的情況下，軋輥轉(zhuǎn)動慣量對帶鋼振頻的影響。軋輥轉(zhuǎn)動慣量較小時，其值的大小對帶鋼振動性態(tài)影響較大，隨著轉(zhuǎn)動慣量的無限增大，邊界條件趨于固支狀態(tài)，軋輥轉(zhuǎn)動慣量對帶鋼振動的影響可近似忽略。

(3) 在機(jī)架間距和軋輥轉(zhuǎn)動慣量一定的情況下，軋制速度對帶鋼振頻的影響。軋制速度較小時，帶鋼振頻受振幅變化影響較大，隨著軋制速度的不斷增大，振頻受振幅影響逐漸減小。

以上所得結(jié)論均基于所建立的具有慣性邊界條件下Euler梁非線性振動模型，可為工程實(shí)際提供一定的理論依據(jù)。

參 考 文 獻(xiàn)

[1] 侯福祥, 張杰, 曹建國, 等. 帶鋼冷軋機(jī)振動問題的研究進(jìn)展及評述[J]. 鋼鐵研究學(xué)報, 2007, 19(10): 6-11.

HOU Fu-xiang, ZHANG Jie,CAO Jian-guo, et al. Review of chatter studies in cold rolling[J]. Journal of Iron and Steel Research, 2007, 19(10): 6-11.

[2] Han S M, Benaroya H, Wei T. Dynamics of transversely vibrating beam using four engineering theories[J]. Journal of Sound and Vibration, 1999, 225(5): 935-988.

[3] 劉明哲, 孫建亮, 彭艷, 等. 連軋過程運(yùn)動帶鋼穩(wěn)定性研究[J]. 冶金設(shè)備, 2009,31(5): 24-28.

LIU Ming-zhe, SUN Jian-liang, PENG Yan, et al. Stabilityresearch of moving strip in the rolling process[J]. Metallurgical Equipment, 2009,31(5): 24-28.

[4] Sun J L, Peng Y, Liu H M. Non-linear vibration and stability of moving strip with time-dependent tension in rolling process[J]. Journal of Iron and Steel Research, International, 2010, 17(6): 11-15.

[5] Bar A, Swiatoniowski A, Interdependence between the rolling speed and non-linear vibrations of the mill system[J]. Journal of Materials Processing Technology, 2004,155-156: 2116-2121.

[6] Miranker W L. Thewave equation in a medium in motion[J]. IBM Journal of Research and Development, 1960, 4: 36-42.

[7] Mote C D Jr. Stability of systems transporting accelerating axially moving materials[J]. Journal of Dynamic Systems Measurement and Control Transactions, 1975, 97: 96-98.

[8] 陳樹輝, 黃建亮. 軸向運(yùn)動梁橫向非線性振動的內(nèi)部共振研究[J]. 力學(xué)學(xué)報, 2005, 37(1): 57-63.

CHEN Shu-hui, HUANG Jian-liang. Study of the laterally nonlinear vibration of axially moving beams[J]. Acta Mechanics Sinica, 2005, 37(1): 57-63.

[9] Suweken G, Van Horssen W T. On the transversal vibrations of a conveyor belt with a low and time-varying velocity, Part II: the beam like case[J]. Journal of Sound and Vibration, 2003, 267(12): 1007-1027.

[10] 王波, 陳立群. 微分求積法處理軸向變速黏彈性梁混雜邊界條件[J]. 振動與沖擊, 2012, 31(5): 87-91.

WANG Bo, CHEN Li-qun. Treatinghybrid boundary condition of an Aaxially accelerating viscoelastic beam via a differential quadrature scheme[J]. Journal of Vibration and Shock, 2012, 31(5): 87-91.

[11] 王洋, 許飛, 李勇, 等. 帶鋼-輥?zhàn)?柔性支承混雜系統(tǒng)動力學(xué)建模與耦合振動研究[J]. 振動工程學(xué)報, 2013, 26(4): 599-607.

WANG Yang, XU Fei, LI Yong,et al, Dynamic modeling and coupling vibration analysis of hybrid systems consisting of strip rolls and flexible supports[J]. Journal of Vibration Engineering, 2013, 26(4): 599-607.

[12] 楊曉東,陳立群. 黏彈性軸向運(yùn)動梁的非線性動力學(xué)行為[J]. 力學(xué)季刊, 2005, 26(1):157-162.

YANG Xiao-dong, CHEN Li-qun.Nonlinear dynamical behaviors of an axially moving visco-elastic beam[J]. Chinese Quarterly of Mechanics, 2005, 26(1):157-162.

[13] Chen L Q, Yang X D. Steady-state response of axially moving viscoelastic beams with pulsating speed: comparison of two nonlinear models[J]. International Journal of Solids and Structures, 2005, 42: 37-50.

[14] 李東, 劉人懷. 修正迭代法在波紋圓板非線性振動問題中的應(yīng)用[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 1990, 11(1): 10-15.

LI Dong, LIU Ren-huai. Application ofmodified iteration algorithm in non-linear vibration problem of circular plate with shallow[J]. Applied Mathematics and Mechanics, 1990, 11(1): 10-15.

[15] 杜國君. 夾層圓板的大幅度振動[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 1994, 15(5): 435-445.

DU Guo-jun. Largeamplitude vibration on circular sandwich plate[J]. Applied Mathamatics and Mechanics, 1994, 15(5): 435-445.

[16] 高崇一, 杜國君, 張忠健. 考慮多間隙影響的軋機(jī)主傳動系統(tǒng)扭振分析[J]. 機(jī)械工程學(xué)報, 2014, 50(3): 130-136.

GAO Chong-yi, DU Guo-jun, ZHANG Zhong-jian, Torsionalvibration analysis of the main drive system of a rolling mill with impact of multi-clearance[J]. Journal of Mechanical Engineering, 2014, 50(3): 130-136.

基金項目：河北省教育廳科學(xué)研究計劃重點(diǎn)項目(ZD2015077)；河北省高等學(xué)校創(chuàng)新團(tuán)隊領(lǐng)軍人才培育計劃項目(LJRC013);國家自然科學(xué)基金資助項目(61403332)

收稿日期：2014-09-02修改稿收到日期：2014-11-06

通信作者杜國君 男，博士，教授，博士生導(dǎo)師，1961年12月生

中圖分類號:TH113.1;TG333.11

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2016.01.002

Nonlinear vibration of strip steel under inertial boundary conditions

GAO Chong-yi1, DU Guo-jun1, LI Jian-xiong2, HU Fa-ke1

(1. Hebei Provincial Key Laboratory of Mechanical Reliability for Heavy Equipment and Large Structures, Yanshan University, Qinhuangdao 066004, China;2. Department of Automation, College of Electrical Engineering, Yanshan University, Qinhuangdao 066004, China)

Abstract:According to the movement mechanism of strip steel and rollers in rolling process of a tandem mill, the strip steel between two neighboring stands was simplified to an axially moving Euler beam and rollers were simplified to inertia components rotating about a fixed axis, the nonlinear vibration mechanical model of the beam under inertial boundary conditions was established. The longitudinal and lateral vibration equations of the beam were derived with Hamilton principle. Kantorovich averaging method was used to simplify the motion equations and the inertial boundaries, and the modified iteration method was employed to solve the motion equations. Finally, based on the results of numerical calculation, the amplitude-frequency response curves of the beam were obtained, and the influences of beam axial velocity, beam length and rotational inertia of roller on the vibration performance of the beam under conditions of inertial boundary were discussed. The results provided an important theoretical reference for controlling and analyzing the vertical vibration of strip steel in a continuous rolling process.

Key words:Euler beam; inertial boundary; strip steel; vertical vibration; modified iteration method

第一作者 高崇一 女，博士生，1986年8月生