一种极薄带钢高速精整机组振动计算方法
【专利摘要】本发明公开了一种极薄带钢高速精整机组振动计算方法，包括下述步骤：(1)初步设定第i#辊辊径Di；(2)建立辊子振动模型，并根据辊子振动模型计算第i#辊的临界线速度Vircr；(3)判断第i#辊的临界线速度Vircr与机组最大运行速度Vmax的大小关系，若Vircr>Vmax成立，执行步骤(4)；若Vircr>Vmax不成立，则令Di＝fεDi，其中系数fε>1，转至步骤(2)；(4)建立带钢振动模型，并根据带钢振动模型计算第i#辊与第(i+1)#辊之间带钢第m阶临界运行速度Vimcr。本发明建立了辊子振动模型和带钢振动模型，形成了比较完善的极薄带钢高速精整机组振动计算方法，以避免机组共振。该方法应用于带钢精整处理，为新机组的开发和控制提供了便利。
【专利说明】一种极薄带钢高速精整机组振动计算方法
【技术领域】
[0001]本发明属于带钢精整处理【技术领域】，涉及一种极薄带钢高速精整机组振动计算方法。
【背景技术】
[0002]经冷轧机组轧制后的带钢，必须经过精整处理加工，才能得到高质量的合格产品。精整机组主要进行重卷、修边等工序。精整处理是成品带钢的最后一道工序，这就要求在精整处理过程中不得产生新的缺陷，带钢表面不准产生擦划伤和塑形变形。
[0003]现代化冷连轧机组不断向高速、自动化方向发展，且运行速度越来越高。为了匹配产量，精整机组的设计速度也越来越高。目前，国内运行速度最高的精整机组是中国重型机械研究院股份公司为武汉钢铁集团公司设计的镀锡准备机组，速度达到1000m/min，在国际上，这也是运行速度最高的精整机组之一，达到国际领先水平。
[0004]机组共振是制约高速精整机组稳定运行的关键因素，高速精整机组一旦发生共振，就会严重影响产品质量，还有可能损坏设备。从已发表的文献得知，目前国内外对高速精整机组振动的研究很少，只停留在一般的概念和生产经验描述上，缺乏深入的理论分析，没有提出具体的计算分析机组振动的原理、方法和步骤，无法应用于机组设计及生产控制。
【发明内容】
[0005]本发明的目的在于提供一种极薄带钢高速精整机组振动计算方法，使得在机组设计和生产控制时可以根据该方法设计和控制机组，以避免机组共振。该方法应用于带钢精整处理，效果明显，为新机组的开发和控制提供了便利。
[0006]一种极薄带钢高速精整机组振动计算方法包括下述步骤:
[0007]步骤一、初步设定第i#辊辊径Di ；
[0008]步骤二、建立辊子振动模型，并根据辊子振动模型计算第i#辊的临界线速度V..? ircr ，
[0009]步骤三、判断第i#辊的临界线速度Vira与机组最大运行速度Vmax的大小关系，若Vira>Vmax成立，执行步骤⑷^Vira>Vmax不成立，则令Di = fEDi，其中系数fE>l，转至步骤-* ;
[0010]步骤四、建立带钢振动模型，并根据带钢振动模型计算第i#辊与第(i+1)#辊之间带钢第m阶临界运行速度ViM。
[0011]进一步地，所述辊子振动模型的建立方法为:
[0012]将辊子简化为两边简支的力学模型进行分析，辊子横向振动时的临界线速度Vircr,单位为m/min,按下式计算
[0013]Vim=SOM2Dl^
[0014]其中，E表示弹性模量，单位为Pa ;G表示辊子重量，单位为kg ;Li表示支座间距，单位为m ;1表示辊子惯性矩，单位为m4，辊子惯性矩I的计算公式
[0015]I =-^{ρι — D0 )+10
[0016]其中，Di为辊外径，单位为m ；D0为辊内径，单位为m ；10为辊轴及辐板惯性矩，单位为m4 ；
[0017]进一步地，所述带钢振动模型的建立方法为:
[0018]截取两个辊子之间长度为ai;宽度为b的一段带钢，a,和b的单位均为m，将带钢和辊子作为一个整体；由于辊子为弹性体，带钢振动力学模型可简化为位于文克勒弹性基础上的简支矩形板，以简支矩形板的长度边为X轴，简支矩形板的宽度边为y轴，带钢的两对边(X = 0，X = at)简支,其余两对边(y = 0，y = b)自由,在两简支边受有大小相等、方向相反的拉力Ni, Ni单位为N。
[0019]进一步地，所述第i#辊与第(i+1)#辊之间带钢第m阶临界运行速度Vimra，单位为m/min,按下式计算
[0020]Vimcr = SODiQffl
[0021]其中，Di为辊外径，单位为m ； ωπ表示带钢第m阶固有振动频率。
[0022]本发明的优点是:
[0023]本发明建立了辊子振动模型和带钢振动模型，可以计算辊子发生振动时的临界线速度和带钢发生振动时的临界运行速度。在设计时，可以通过调整辊径Di使临界线速度大于机组最大速度；在机组传动控制时，可以通过控制机组运行速度和张力避免带钢发生振动。从而，形成比较完善的极薄带钢高速精整机组振动计算方法，计算快速稳定。
[0024]以下将结合附图和实施例对本发明做进一步详细说明。
【专利附图】
【附图说明】
[0025]图1是辊子振动力学模型图。
[0026]图2是带钢振动力学模型图。
[0027]图3是极薄带钢高速精整机组振动计算方法流程图。
[0028]图4是带钢临界运行速度Viira随带钢长度a变化曲线图。
[0029]图5是带钢临界运行速度Viira随带钢宽度b变化曲线图。
[0030]图6是带钢临界运行速度Viira随带钢厚度h变化曲线图。
[0031]图7是带钢临界运行速度Viira随带钢张力T变化曲线图。
【具体实施方式】
[0032]实施例1:
[0033]本实施例所述的极薄带钢高速精整机组振动计算方法，包括以下步骤:
[0034]步骤一、初定第i#辊辊径Di。
[0035]步骤二、建立辊子振动模型，并根据辊子振动模型计算第i#辊的临界线速度Vira。辊子振动模型的建立方法如下:
[0036]如图1所示，辊子简化为两边简支的力学模型进行分析，辊子横向振动时的临界线速度Vira?，单位:m/min,按下式计算:
[0037]
ν.= 30;τ2Α 擇")
[0038]其中，E表示弹性模量，单位为Pa ;G表示辊子重量，单位为kg A表示支座间距，单位为m ;1表示辊子惯性矩，单位为m4，辊子惯性矩I的计算公式:
[0039]/=吾(《)+/0
[0040]其中，Di为辊外径，单位为m ；D0为辊内径，单位为m ；10为辊轴及辐板惯性矩，单位为m4。
[0041]步骤三、判断第i#辊的临界线速度Vira与机组最大运行速度Vmax的大小关系，若Vira>Vmax成立，执行步骤⑷^Vira>Vmax不成立，则令Di = fEDi，其中系数fE>l，转至步骤
--O
[0042]步骤四、建立带钢振动模型，并根据带钢振动模型计算第i#辊与第(i+1)#辊之间带钢第m阶临界运行速度ViM。
[0043](I)带钢振动模型的建立
[0044]如图2所示，截取两个辊子之间长为％ (单位:m)的一段带钢，将带钢和辊子作为一个整体进行研究。由于辊子为弹性体，其力学模型可简化为位于文克勒弹性基础上的简支矩形板，其两对边(X = 0，X = aj简支,其余两对边(y = 0，y = b)自由，带钢宽度b(单位:m)。在两简支边受有大小相等、方向相反的拉力Ni (单位:N)。
[0045](2)位移振型函数模式
[0046]对于带钢固有振动，位移振型函数设为:
[0047]
W {x,,V) = Y,M (y)sin(m=l,2,3，---，n)C2)
ai
[0048](3)薄板横向振动方程
[0049]
O4IF c4W M M C1WL.η,^vtr f 一 λm
-十 ￡---? -1--；--(X Vv ----— — U\3 J
OX cx-oy- cyD ex'
__1、1.4k."
[0050]式中，—D= J^,
[0051]其中，k表示文克勒弹性基础模量，单位:Kg/m3，其取值取决于辊子刚性；w表示带钢固有振动频率，单位:rad/s ;r表示带钢密度,单位:Kg/m3 ;h表示带钢厚度,单位:m ^表示泊松比，V = 0.3。
[0052]对式(3)求导，得到
[0053]
? 厂 /\4{\
ο Iy I mm tr (.ηιπf Λ,
^——> |;# I vjsin ——'X(4)
fir I , J ^J
[0054] f \2 f \
iWπιπ \ %r w, X , mm,产、
_;τ = - - Ym I f)sie -- X(ο)
& 负.-、,J ^ ^ J
[0055]
O4IF ,, ftff f 、.(mm\,,、
——~-ι, V sm1- X(6)
^ 4WlfI
?ν、a.J
[0056]式(4)、(5)、(6)代入式(3)简化得到
[0057]
C W — 2(4 C (V) + [M4^4+ fi^]" I r?, (>>) = O(7)
I Λ, J K ^ J dK a, J —
[0058]式(7)为常系数常微分方程，其解为
[0059]Ym (>!) = eji'x(8)
[0060]式中，Idi11应满足特征方程
[0061]
/Γ/ / \2*.,4 J mKI nm 4 Ni ιηπ Α
Pm — 2 — pm 十 ? — ~~ OC 十—- 二 Oi ))
IJ [I J 0Kai) _
[0062]求解式(9)，得到:
[0063]
A-=JR7^lRf(10)
II?.J I DKa, J
[0064]
ιηπ\Ni ( Μπ/I 1、
P%m— 4十 4 α ~7Γ\1 1Z
?, J $ D{a, }
[0065]
I=滅-f^￥f(12)
[0066]
β — M L^f-T(|3)
w IUj i ο UJ
[0067]故有一般解
[0068]
K (ν) = AmBfh^ + Bmeih'}' + C?/Vr + Dn/^(14)
[0069]从而，得到振型函数
[0070]if (,τ,_ν) =[jweA-y +BmS0tmr +Cmefl,-y + Dme0*~r(15)
[0071]式中，系数Am、Bm、Cm、Dm由简支矩形板两自由边的边界条件确定。
[0072](4)简支矩形板边界条件
[0073]简支矩形板简支边的混合边界条件为
[0074]
Ηλ-?.? =0
< (O2w')(16)
Id / ?
\λ ~(ha
[0075]振型函数W (x，y)恒满足上式。
[0076]在自由边界上各点的弯矩、扭矩、剪力均为零，得到矩形板自由边的静力边界条件
r ( '2 \
OVn
_ I % ? _ Immf I I
I V Γ~— KJ
^ 2 ^ 2 I
[0077]Ir^J、(17)
i^i+(2 —=0
I^ ￡a.cy,U,<.
[0078](5)求解固有振动频率的矩阵方程
[0079]对式(15)求导得到:
[0080]l~^T ^m fiirnelhmy +^,,JLe0imf +CM/4cAeJ +
α> (18)
[0081]
^r = Fjbj--Λ.*.ν +BmBfilmy +Cmefiim,y + Dme04m}' l1-?-^-^-1 jsin-^-.v^^(19)
[0082]^^Ia,ηβι,η^fk^ +B?Jl?eA-'r +CmfilHeP1-"y + Dmfilnep^v ]sin^.t
呀(20)
[0083]
^?J\me,kmr +BnJ2me^v 七C,▲?/、’γ +Dmfi,meA-y ] -?—) sin^-x
α\- 'CV LJ I a, I a,
-Lv-yJ.(21)
[0084]式(18)、(19)、(20)、(21)代入式(17)得到
[0085]y = 0 时，
AJls + K0L + CJlt + DtJl -(Am + Bm + Cm + Dm): O(22)
[0086]y = b 时，
++ CJil,/一 *DaflleA-b
r \2
-v I (Ame/h"b + B Sfitmh + CmZ1-* + Dme"',) = O(23)
KaJx>
[0087]y = 0 时，
AmfiL + + CmfiL.+ ~(2~ v)^—j ( 'H + Smfi2m + Cmfi3m + H ) = 0
(24)
[0088]y = b 时，++ Cm^ + D,Xne^
[0089]
—(2 — )" (A,An/-h + Β,Λ,,,^ + CJl,me^h + D,J4,,-h) = 0(25)
[0090]将式(22)、(23)、(24)、(25)写成矩阵表达式:
[0091][T] {A} = {0} (26)
rd\\ ■.■ ?4、
[0092]式中,[Γ]=:..: , {A} = (Am, Bm, Cm, Dm)To
V 1 …^44 /
[0093]采用迭代法求解，若某一 ωπ值使得行列式|τ| = 0，则为所求第m阶固有振动频率。
[0094](6)第m阶带钢临界运行速度Vimcr
[0095]带钢振动源为辊子转动，因此按照带钢固有振动频率ωπ可求得发生振动时，带钢第m阶临界运行速度Vimra,单位:m/min
[0096]Vimcr = SODiWm (27)
[0097]其中，Di为辊外径，单位为m ； ωπ表示带钢第m阶固有振动频率。
[0098]综上，本实施例建立了辊子振动模型和带钢振动模型，可以计算辊子发生振动时的临界线速度和带钢发生振动时的临界运行速度。在设计时，可以通过调整辊径Di使临界线速度大于机组最大速度；在机组传动控制时，可以通过控制机组运行速度和张力避免带钢发生振动，从而，形成比较完善的极薄带钢高速精整机组振动计算方法。
[0099]实施例2:
[0100]通过采用本发明所提出的极薄带钢高速精整机组振动计算方法对某镀锡准备机组工艺参数进行计算分析。本机组带材规格为:宽度b = 700?1300mm,厚度h = 0.15?
0.55mm,机组最大运行速度Vmax = 1000m/min,最小卷取张力Tmin = 6000N,最大卷取张力Tmax=18000N。主要计算结果如表I。
[0101]表I机组典型规格辊子临界线速度
[0102]辊径，mm[?20~[?60~[230~[300
临界线速度，m/min 2142 2556 5276 9239~
[0103]可见，辊子临界线速度都远远大于机组最大运行速度，机组运行时，辊子不会发生振动。
[0104]带钢临界运行速度Vimra分别随带钢长度a、带钢宽度b、带钢厚度h、带钢张力T变化曲线见图4、图5、图6、图7。从图中可见，带钢临界运行速度Viira随带钢长度a的增加而减小、随带钢宽度b的增加而减小、随带钢厚度h的增加而减小、随带钢张力T的增加而增加。
[0105]在机组设计时，可以根据计算结果设计辊径Ditl在机组控制时，可以控制机组速度快速通过带钢临界运行速度Vimra，或者改变带钢张力T以使机组能够以设定的运行速度生产，避免机组发生共振。否则，机组发生共振，影响产品质量，严重时，损坏机械设备。
[0106]上面结合附图对本发明的实施方式作了说明，但本发明并不限于上述实施方式，在本领域的普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。
【权利要求】
1.一种极薄带钢高速精整机组振动计算方法，其特征在于包括下述步骤: (1)初步设定第i#辊辊径Di； (2)建立辊子振动模型，并根据辊子振动模型计算第i#辊的临界线速度Vira； (3)判断第i#辊的临界线速度Vira与机组最大运行速度Vmax的大小关系，若V1T>Vmax成立，执行步骤⑷^V1T>Vmax不成立，则令Di = LDi，其中系数fE>l，转至步骤⑵； (4)建立带钢振动模型，并根据带钢振动模型计算第i#辊与第(i+1)#辊之间带钢第m阶临界运行速度Vimra。
2.根据权利要求1所述的极薄带钢高速精整机组振动计算方法，其特征在于，所述辊子振动模型的建立方法为: 将辊子简化为两边简支的力学模型进行分析，辊子横向振动时的临界线速度Vira，单位为m/min,按下式计算
i Pf
\τ ?α w- η I
^ircr 30/Γ Di I (I γ 其中,E表示弹性模量,单位为Pa ；G表示棍子重量,单位为kg ；Li表示支座间距,单位为m ；1表示辊子惯性矩，单位为m4，辊子惯性矩I的计算公式
I=^(A^D/)+/o 其中，Di为棍外径,单位为m ;D(!为棍内径,单位为m ；10为棍轴及辐板惯性矩,单位为m4。
3.根据权利要求1所述的极薄带钢高速精整机组振动计算方法，其特征在于，所述带钢振动模型的建立方法为: 截取两个辊子之间长度为ai;宽度为b的一段带钢，％和b的单位均为m，将带钢和辊子作为一个整体；由于辊子为弹性体，带钢振动力学模型可简化为位于文克勒弹性基础上的简支矩形板，以简支矩形板的长度边为X轴，简支矩形板的宽度边为y轴，带钢的两对边(X = O, X = at)简支,其余两对边(y = O, y = b)自由,在两简支边受有大小相等、方向相反的拉力NyNi单位为N。
4.根据权利要求1所述的极薄带钢高速精整机组振动计算方法，其特征在于，所述第i#辊与第(i+1) #辊之间带钢第m阶临界运行速度Vimra,单位为m/min，按下式计算
Vimcr = SODi ωηι 其中，Di为辊外径，单位为m ; CO111表示带钢第m阶固有振动频率。
【文档编号】B21B38/00GK104128376SQ201410283031
【公开日】2014年11月5日 申请日期:2014年6月23日 优先权日:2014年6月23日
【发明者】孙亚波, 刘渭苗, 张康武, 任玉成 申请人:中国重型机械研究院股份公司