戴杰濤　李烈軍　戴文笠　張祖江

(1.華南理工大學(xué) 機(jī)械與汽車工程學(xué)院, 廣東 廣州510640; 2.寶鋼集團(tuán) 廣東韶關(guān)鋼鐵有限公司, 廣東 韶關(guān) 512123)

張力下的板形屈曲臨界載荷和寬度分析*

戴杰濤1李烈軍1戴文笠2張祖江2

(1.華南理工大學(xué) 機(jī)械與汽車工程學(xué)院, 廣東 廣州510640; 2.寶鋼集團(tuán) 廣東韶關(guān)鋼鐵有限公司, 廣東 韶關(guān) 512123)

摘要：針對(duì)現(xiàn)場(chǎng)板帶材的實(shí)際板形瓢曲情況,在分析以往板帶材屈曲按全板寬進(jìn)行計(jì)算的解析模型的基礎(chǔ)上,文中采用壓應(yīng)力區(qū)域作為板形瓢曲變形區(qū)域,以反映現(xiàn)場(chǎng)的實(shí)際板形情況,并據(jù)此建立了張力影響板形瓢曲變形的解析計(jì)算模型,獲得了不同張力大小情況下板帶材的臨界屈曲載荷與寬度的變化情況:臨界屈曲載荷隨張力的增大而增大,而臨界屈曲寬度隨張力的增大而減小,與現(xiàn)場(chǎng)的實(shí)際情況一致.為進(jìn)一步驗(yàn)證解析計(jì)算結(jié)果,建立了張力作用下板帶材瓢曲的有限元計(jì)算模型,仿真計(jì)算結(jié)果與解析計(jì)算結(jié)果基本一致.

關(guān)鍵詞：板帶材;張力;屈曲;有限元

張力對(duì)整體浪形的影響一直是板形瓢曲計(jì)算領(lǐng)域的一個(gè)薄弱環(huán)節(jié),因?yàn)樵趯?shí)際的生產(chǎn)過(guò)程中,帶鋼是一個(gè)整體,隨著張力的增大,板帶材的壓應(yīng)力會(huì)隨之變化,導(dǎo)致板形瓢曲區(qū)域也發(fā)生變化.同時(shí),有限元仿真計(jì)算結(jié)果也表明,在板形瓢曲計(jì)算中,張力越大,屈曲寬度越小.現(xiàn)有張力對(duì)板形瓢曲變形影響的研究主要集中在連續(xù)退火爐內(nèi)張力對(duì)熱瓢曲的影響[1-5],有關(guān)張力對(duì)熱、冷連軋過(guò)程中板形瓢曲變形影響的研究相對(duì)較少[6-9],最常見(jiàn)的方法是在解析研究過(guò)程中忽略張力對(duì)板形瓢曲變形寬度的影響,認(rèn)為在不同張力情況下板形瓢曲變形的寬度保持不變[6,9],計(jì)算結(jié)果不能有效地反映板帶材張力對(duì)板形瓢曲變形寬度的影響.此外,以往對(duì)現(xiàn)場(chǎng)板形瓢曲進(jìn)行解析計(jì)算時(shí),均認(rèn)為鋼板的屈曲區(qū)域是整體寬度的[10-14],但在實(shí)際的生產(chǎn)過(guò)程中,不管有沒(méi)有張力,板帶材的瓢曲區(qū)域均沒(méi)有達(dá)到整體寬度,只是在部分寬度內(nèi)發(fā)生,因此,如果按整體寬度對(duì)板帶材的臨界屈曲載荷進(jìn)行計(jì)算,獲得的臨界屈曲載荷往往要小于實(shí)際的板形臨界屈曲載荷值.以往按整體寬度對(duì)冷軋后薄帶鋼板形瓢曲行為進(jìn)行計(jì)算,獲得的臨界屈曲載荷均很小[8-9,11,14],以寬為1 600 mm、厚為1 mm的鋼板為例,按文獻(xiàn)[11]方法計(jì)算獲得的臨界屈曲載荷為1.512 MPa,折算到板形浪形為0.2 IU(1 IU表示100 m長(zhǎng)的板帶材縱向延伸差為1 mm).按后屈曲計(jì)算,如果其浪形達(dá)到5 IU,則其浪形高度將達(dá)到10 mm以上,而實(shí)際浪形的控制標(biāo)準(zhǔn)為高度在3 mm以內(nèi),相應(yīng)的浪形在5 IU以內(nèi),顯然上述計(jì)算結(jié)果不符合現(xiàn)場(chǎng)的實(shí)際情況.

針對(duì)上述問(wèn)題,文中采用壓應(yīng)力區(qū)作為板形瓢曲變形區(qū)域?qū)ΜF(xiàn)場(chǎng)實(shí)際板形情況進(jìn)行跟蹤分析,建立了張力影響板形瓢曲變形的解析計(jì)算模型,以解決采用全板寬計(jì)算導(dǎo)致臨界屈曲載荷偏小的問(wèn)題;同時(shí)將張力引入解析模型,以解決以往解析模型在考慮張力時(shí)不能有效反映張力對(duì)板形瓢曲寬度影響的問(wèn)題.為進(jìn)一步驗(yàn)證解析計(jì)算結(jié)果,文中還建立了張力作用下板帶材瓢曲變形的有限元仿真分析模型,以驗(yàn)證解析計(jì)算模型的正確性.

1張力對(duì)板形瓢曲影響的解析計(jì)算

1.1　張力作用下板材的載荷函數(shù)

根據(jù)現(xiàn)場(chǎng)鋼板的板形應(yīng)力情況,以板中心為坐標(biāo)原點(diǎn),y表示縱向位置,x表示橫向位置,可以將鋼板的板形載荷定義為(文中在計(jì)算時(shí)均以中浪為例進(jìn)行分析),

(1)

式中:Ny為鋼板沿寬度方向的載荷分布值;N0為邊界不均勻分布載荷尖峰值;b為鋼板寬度;l、m、n為常數(shù),其數(shù)值根據(jù)板形儀上顯示的板形曲線進(jìn)行擬合得到;dN為待定常數(shù).

考慮到板形應(yīng)力應(yīng)滿足自平衡條件[8-12],即

(2)

可以求出式(1)中待定常數(shù)dN:

(3)

設(shè)板材受到的張力為Nz,則板材的載荷為

(4)

1.2　板形瓢曲區(qū)域的分析

在以往的板形瓢曲計(jì)算過(guò)程中,往往認(rèn)為板形的瓢曲區(qū)域?yàn)檎w寬度,而從圖1所示的中浪及板形應(yīng)力圖可以看出,鋼板的起浪區(qū)域不是整體寬度,而是受壓應(yīng)力作用的區(qū)域.同時(shí)從受力角度分析也可以知道,鋼板在受拉應(yīng)力區(qū)域不能發(fā)生板形瓢曲變形.因此,文中在對(duì)板形瓢曲進(jìn)行解析計(jì)算時(shí),將鋼板的瓢曲區(qū)域限定在受壓應(yīng)力區(qū)域.

圖1　在拉、壓應(yīng)力作用下板材的中浪變形實(shí)物圖Fig.1　Photo of plate central wave deformation under tensile and compressive stress

因此,對(duì)于式(4)所示的載荷情況,板材的瓢曲區(qū)域可以通過(guò)下式獲得:

(5)

bw=f (b,Nz/N0)

(6)

1.3　浪形函數(shù)

從上述浪形瓢曲區(qū)域的分析可以看出,發(fā)生浪形瓢曲變形的區(qū)域?yàn)棣?[-bw,bw;-a,a],而非整個(gè)方板區(qū)域Ω[-b,b;-a,a],屈曲區(qū)域的長(zhǎng)度a為浪形的半波長(zhǎng),其具體取值由能量最低原理確定.

根據(jù)現(xiàn)場(chǎng)的實(shí)際浪形模態(tài)特征,取板形瓢曲區(qū)域兩側(cè)邊為固支約束,則撓度函數(shù)W(x,y)應(yīng)滿足如下邊界條件:

(7)

在滿足上述邊界條件的情況下,根據(jù)現(xiàn)場(chǎng)浪形函數(shù)案例的實(shí)測(cè),浪形的撓度函數(shù)可描述成如下數(shù)學(xué)表達(dá)式:

(8)

式中,A為撓度幅值.

撓度函數(shù)W的仿真結(jié)果如圖2所示,顯然此撓度函數(shù)更符合圖1所示的實(shí)際板形情況.

圖2　撓度函數(shù)仿真結(jié)果Fig.2　Simulation result of deflection function

1.4　臨界屈曲載荷的求解

根據(jù)小撓度理論和能量原理[14-15],屈曲時(shí)刻對(duì)應(yīng)的臨界狀態(tài)由式(9)獲得:

T=U

(9)

其中,微彎狀態(tài)下矩形板帶材的彎曲變形能U為

(10)

作用在中面內(nèi)的中面內(nèi)力做功T為

(11)

式中，Nx為鋼板沿長(zhǎng)度方向的載荷分布值，Nxy為鋼板剪切載荷分布值.文中討論的現(xiàn)場(chǎng)鋼板情況下的Nx和Nxy值均為0.

由此可知,式(7)-(9)構(gòu)成了求解屈曲載荷N0的方程組.獲得N0的表達(dá)式

(12)

引入

θ=a/bw

(13)

根據(jù)臨界屈曲載荷最小,對(duì)臨界載荷N0求θ的偏導(dǎo),則有

(14)

從而獲得θ=2/3.

將θ=2/3及式(6)、(13)代入式(12)中,可以獲得

N0=f(Nz,b,h)

(15)

在給定張力Nz的情況下,通過(guò)式(15)獲得臨界屈曲載荷后代入式(6)中,可以進(jìn)一步獲得張力作用下板材的臨界屈曲寬度.

1.5　與整體寬計(jì)算結(jié)果對(duì)比

為了驗(yàn)證采用壓應(yīng)力區(qū)域?qū)挾茸鳛榍鷮挾扔?jì)算結(jié)果的正確性,針對(duì)現(xiàn)場(chǎng)某塊規(guī)格為板寬2b=1 600 mm、厚度h=1.0 mm的鋼板進(jìn)行了跟蹤,此鋼板在軋制過(guò)程中板形儀上顯示鋼板的板形應(yīng)力為4.89 IU,板形應(yīng)力函數(shù)如式(16)所示,在檢測(cè)平臺(tái)上測(cè)量得到的浪高大約為3.3 mm.

Ny=N0x2

(16)

針對(duì)上述現(xiàn)場(chǎng)的實(shí)際鋼板,在張力為0的情況下,按壓應(yīng)力區(qū)域屈曲計(jì)算獲得的臨界屈曲載荷為4.284 MPa,而以整體寬計(jì)算獲得的臨界屈曲載荷為1.512 MPa.為了更好對(duì)比兩者的計(jì)算結(jié)果,采用攝動(dòng)法[14-15]對(duì)后屈曲路徑進(jìn)行計(jì)算(具體計(jì)算過(guò)程在此不作贅述),按壓應(yīng)力區(qū)域屈曲在延伸差為4.89 IU(即內(nèi)應(yīng)力載荷為10.27 MPa)的情況下,獲得的浪高為4.16 mm,而按整體寬計(jì)算獲得的浪高為10.35 mm,顯然按壓應(yīng)力區(qū)域進(jìn)行計(jì)算更符合現(xiàn)場(chǎng)的實(shí)際情況.

1.6　算例分析

針對(duì)1.5節(jié)中的鋼板情況,文中計(jì)算了不同張力情況下板形的臨界屈曲載荷和寬度.

在彈性模量E=210 GPa、泊松比μ=0.3的情況下,板形的臨界屈曲載荷和寬度隨張力的變化如圖3所示.從圖中可以看出:板帶材的臨界屈曲載荷隨著張力的增大而增大,并且近似呈直線增長(zhǎng);隨著張力的增大,板帶材的臨界屈曲寬度迅速降低,達(dá)到一定的值后下降趨勢(shì)趨于平緩.

圖3　臨界屈曲載荷和寬度隨張力的變化曲線Fig.3　Changing curves of critical buckling load and width with outer tension

為了進(jìn)一步驗(yàn)證上述計(jì)算結(jié)果,對(duì)活套上同樣規(guī)格的鋼板進(jìn)行了不同張力調(diào)整試驗(yàn),從試驗(yàn)中可以看出,當(dāng)張力增大時(shí),活套上鋼板的浪形寬度會(huì)變小,這與圖3(a)的計(jì)算結(jié)果一致,同時(shí)隨著張力的增大,本來(lái)有浪形的鋼板會(huì)被拉直(即鋼板的臨界屈曲載荷增大),這與圖3(b)的計(jì)算結(jié)果也一致.這從側(cè)面驗(yàn)證了文中建立的張力對(duì)板形屈曲影響解析模型的正確性.

2張力對(duì)板形屈曲影響的有限元仿真

為了進(jìn)一步分析張力對(duì)板形屈曲的影響,也進(jìn)一步驗(yàn)證解析計(jì)算過(guò)程的正確性,文中應(yīng)用ANSYS有限元軟件對(duì)承受張力的板帶材進(jìn)行了仿真分析[16].

2.1　模型參數(shù)與網(wǎng)格劃分

由于Shell63單元具有彎曲和薄膜兩種功能,可以承受面內(nèi)載荷和法向載荷,同時(shí)其已經(jīng)考慮了應(yīng)力強(qiáng)化和大變形,因此文中采用Shell63單元來(lái)劃分網(wǎng)格.

選取材料性能參數(shù)如下:彈性模量為210 GPa,泊松比為0.3,帶鋼沿各方向的熱膨脹系數(shù)為1.2×10-6/℃,厚度為0.001 m.

這里可以直接建立一個(gè)矩形面積作為帶鋼模型,并令原點(diǎn)為模型的中心點(diǎn),x、y方向分別為帶鋼的寬度方向(橫向)和長(zhǎng)度方向(縱向),面積尺寸初步定為8 000 mm×1 600 mm(長(zhǎng)度×寬度,分析過(guò)程中寬度保持不變),網(wǎng)格數(shù)為30×150.

2.2　約束及邊界載荷

根據(jù)軋制理論可以知道，帶鋼的縱向瓢曲浪形是由于帶鋼縱向延伸在寬度方向上的分布不一致.為了描述這種不一致的縱向延伸分布,文中采用不均勻的溫度場(chǎng)來(lái)描述帶鋼的不均勻延伸,同時(shí)對(duì)載荷離散化后將離散化的數(shù)據(jù)添加到每個(gè)節(jié)點(diǎn)上,這樣就構(gòu)成了中浪求解的載荷條件.在帶鋼的端邊界上以節(jié)點(diǎn)載荷的形式給定張力,其合力與帶鋼的張力大小一致.對(duì)于模型的邊界約束,以中浪為例,給定圖4所示的約束條件.

圖4　張力下板帶材屈曲的有限元模型Fig.4　Finite element model of strip plate buckling under tension

2.3　張力條件下板形縱向瓢曲模態(tài)

圖5給出了張力條件下的鋼板瓢曲浪形模態(tài),從圖中可以看出,在張力條件下鋼板的浪形模態(tài)在寬度方向上不是鋼板的整體寬度,而是部分寬度,這驗(yàn)證了文中對(duì)張力條件下帶鋼采用壓應(yīng)力區(qū)域而非整體寬度進(jìn)行計(jì)算的正確性.

圖5　張力下板帶材的屈曲模態(tài)Fig.5　Buckling modal of strip plate under tension

2.4　張力對(duì)臨界屈曲載荷和寬度的影響

不同張力情況下板帶材的臨界屈曲載荷和寬度的解析計(jì)算結(jié)果與有限元仿真結(jié)果比較如圖6所示,從圖中可以看出:在不同張力情況下,板帶材的臨界屈曲載荷隨著張力的增加而增加,有限元仿真計(jì)算結(jié)果小于解析計(jì)算結(jié)果,這是因?yàn)樵诮馕鲇?jì)算過(guò)程中,固支約束邊界給定在板帶材的實(shí)際邊部,而不是在發(fā)生瓢曲區(qū)域的邊部,這使得實(shí)際發(fā)生瓢曲區(qū)域的邊界約束要比解析計(jì)算的邊界約束弱;隨著張力的增大,板帶材的臨界屈曲寬度減少,有限元仿真結(jié)果大于解析計(jì)算結(jié)果,這同樣是由于有限元計(jì)算過(guò)程中添加的約束條件與解析計(jì)算中添加的約束條件不一致造成的.

圖6　臨界屈曲栽荷和寬度的解析計(jì)算結(jié)果與有限元仿真結(jié)果比較Fig.6　Comparison of critical buckling load and width betweenresolution results and finite element simulation results

3結(jié)論

文中采用壓應(yīng)力區(qū)域而不是全板寬作為板形瓢曲變形區(qū)域?qū)ΜF(xiàn)場(chǎng)實(shí)際板形情況進(jìn)行跟蹤分析,建立了張力影響板形瓢曲變形的解析計(jì)算模型,獲得了不同張力情況下板帶材的臨界屈曲載荷和寬度的變化情況.為進(jìn)一步驗(yàn)證解析計(jì)算結(jié)果,文中還建立了張力作用下板帶材瓢曲變形的有限元仿真分析模型,獲得的仿真結(jié)果驗(yàn)證了解析計(jì)算模型的正確性.

參考文獻(xiàn):

[1]劉順明,尹顯東,鄭曉飛,等.退火爐內(nèi)帶鋼拉伸失穩(wěn)與屈曲失穩(wěn)分析 [J].軋鋼,2013,30(6):22-25.

Liu Shun-ming,Yin Xian-dong,Zheng Xiao-fei,et al.Ana-lysis of tensile instability and buckling instability of strip in the furnace of continuous annealing line [J].Steel Rolling,2013,30(6):22-25.

[2]李會(huì)免,劉建雄,李俊洪.連續(xù)退火爐內(nèi)帶鋼橫向張力分布研究 [J].材料熱處理技術(shù),2010,39(4):148-150.

Li Hui-mian,Liu Jian-xiong,Li Jun-hong.Research of transverse tension distribution on steel strip in continuous-annealing furnace [J].Material & Heat Treatment,2010,39(4):148-150.

[3]戴杰濤,李烈軍,張祖江,等.中厚板寬板軋制中邊浪的產(chǎn)生機(jī)理與抑制措施 [J].華南理工大學(xué)學(xué)報(bào):自然科學(xué)版,2014,42(12):64- 69.

Dai Jie-tao,Li Lie-jun,Zhang Zu-jiang,et al.Generation mechanism and restrain measure of edge wave on the medium plate rolling [J].Journal of South China University of Technology:Natural Science Edition,2014,42(12):64- 69.

[4]胡廣魁.連續(xù)退火爐內(nèi)鋼帶穩(wěn)定運(yùn)行張力分析 [J].寶鋼技術(shù),2010(5):47-50.

Hu Guang-kui.Tension analysis for stable strip running in a continuous annealing furnace [J].Baosteel Technology, 2010(5):47-50.

[5]戴江波.冷軋寬帶鋼連續(xù)退火線生產(chǎn)線上瓢曲變形的研究 [D].北京:北京科技大學(xué)機(jī)械工程學(xué)院,2005.

[6]林振波.冷帶軋機(jī)板形判別模型的有限條分析 [D].河北:燕山大學(xué)機(jī)械工程學(xué)院,1993.

[7]邊宇虹,劉宏民.求解帶材軋后大撓度屈曲變形的一個(gè)通用方法 [J].機(jī)械工程學(xué)報(bào),1994,30(增刊):21-27.

Bian Yu-hong,Liu Hong-min.Universal method analysing the large deflection buckling deformation of rolled strip [J].Chinese Journal of Mechanical Engineering,1994, 30(Suppl):21-27.

[8]楊荃,陳先霖.軋制帶材的瓢曲生成路徑 [J].北京科技大學(xué)學(xué)報(bào),1994,16(1):53-57.

Yang Quan,Chen Xian-lin.The deforming route of buckled waves of rolled strip [J].Journal of University of Science and Technology Beijing,1994,16(1):53-57.

[9]楊荃.冷軋帶鋼屈曲理論與板形控制目標(biāo)的研究 [D].北京:北京科技大學(xué)機(jī)械工程學(xué)院,1992.

[10]Poplawski J V,Jr Seccombe D A.Bethlehem’s contribution to the mathematical modeling of cold rolling in tandem mills [J].Iron Steel Engineering,1980,57(9):47-58.

[11]常鐵柱.薄寬帶鋼板形斜向和橫向屈曲變形行為研究 [D].北京:北京科技大學(xué)機(jī)械工程學(xué)院,2009.

[12]Roberts W L.Cold rolling of steel [M].New York:Marcel Dekker Inc,1978.

[13]Tarnopolskaya T,Gates D J.Analysis of the effect of strip buckling on stability of strip lateral motion with application to cold rolling of steel [J].Journal of Dynamic Systems Measurement and Control,2008,130(1):1-7.

[14]戴杰濤,張清東.冷軋薄板中浪板形缺陷的屈曲及后屈曲理論與軋制試驗(yàn)研究 [J].機(jī)械工程學(xué)報(bào),2011, 47(2):44-50.

Dai Jie-tao,Zhang Qing-dong.Analysis and experimental on central buckling and post buckling of thin cold-rolled sheet [J].Journal of Mechanical Engineering,2011,47(2):44-50.

[15]戴杰濤,張清東,秦劍.薄寬冷軋帶鋼局部板形屈曲行為解析研究 [J].工程力學(xué),2011,28(10):236-242.

Dai Jie-tao,Zhang Qing-dong,Qin Jian.Analysis of local buckling for thin cold-rolled strip [J].Engineering Mechanics,2011,28(10):236-242.

[16]卿偉杰,楊荃.冷軋帶鋼整體和局部屈曲及后屈曲的有限元分析 [J].北京科技大學(xué)學(xué)報(bào),2000,22(4):377-380.

Qing Wei-jie,Yang Quan.Study on cold rolled strip glo-bal and local buckling,post-buckling using the finite element method [J].Journal of University of Science and Technology Beijing,2000,22(4):377-380.

Analysis of Critical Load and Width of Plate-Shaped Buckling under Tension

DaiJie-tao1LiLie-jun1DaiWen-li2ZhangZu-jiang2

(1. School of Mechanical and Automotive Engineering, South China University of Technology, Guangzhou 510640, Guangdong,

China; 2. Shaoguan Iron and Steel Co.Ltd., Baosteel Group, Shaoguan 512123, Guangdong, China)

Abstract:Based on the actual plate and strip buckling situation in the site and by analyzing the analytic model of plate and strip buckling, which used to be calculated as a whole, this paper takes the compressive stress region as the plate-shaped buckling deformation region to reflect the actual plate shape situation in the site, and then establishes an analytic calculation model of the outer tension affecting the plate-shaped buckling deformation. Thus, the changes of the critical buckling load and critical width of plate and strip are obtained under different outer tensions, that is, with the increase of tension, the critical buckling load increases while the critical width decreases, which accords well with the actual situation in the site. In order to further verify the analytic calculation result, a finite element calculation model of the plate and strip buckling is established under outer tensions, and the simulation result is consistent with the analytic calculation one.

Key words:plate and strip; tension; buckling; finite element

中圖分類號(hào)：TG335.5

doi：10.3969/j.issn.1000-565X.2015.09.017

作者簡(jiǎn)介：戴杰濤(1982-)，男，博士后，工程師，主要從事冶金機(jī)械力學(xué)行為和板帶鋼板形的機(jī)理分析與控制研究.E-mail： djt66008@163.com

*基金項(xiàng)目：中國(guó)博士后科學(xué)基金資助項(xiàng)目(2014M562168)

收稿日期：2015-03-10

文章編號(hào)：1000-565X(2015)09-0107-06

Foundation item： Supported by the China Postdoctoral Science Foundation(2014M562168)