管奔，臧勇，逄曉男，呂智勇，曲為壯

(1. 北京科技大學(xué) 機(jī)械工程學(xué)院，北京，100083；2. 山東省萊蕪鋼鐵股份有限公司型鋼廠，山東 萊蕪，271126)

彈塑性彎曲是一種重要的金屬加工技術(shù)。在現(xiàn)代工業(yè)生產(chǎn)中，板條、型材和板的彈塑性彎曲被廣泛地應(yīng)用于大型金屬結(jié)構(gòu)及各種日常用品的制造過程[1]，而在金屬條材的生產(chǎn)過程中，作為最后一道塑性變形性質(zhì)的輥式矯直工序更是典型的需要經(jīng)歷多次反彎的復(fù)雜彈塑性彎曲過程[2]。為設(shè)定輥式矯直過程合理的工藝參數(shù)及保證矯后條材的幾何精度，需要對輥式矯直過程所特有的此類連續(xù)彈塑性反彎過程進(jìn)行深入研究和準(zhǔn)確認(rèn)識。針對彈塑性彎曲，國內(nèi)外學(xué)者進(jìn)行了許多基于理論探索和工程應(yīng)用的研究工作[3-5]，大部分工作是以一次彎曲成型的金屬塑性加工為應(yīng)用背景，從而不存在考慮彎曲歷史對截面特性的影響問題。實(shí)際上，金屬條材的矯直是典型的多次彈塑性反彎的力學(xué)過程，矯直理論研究的力學(xué)基礎(chǔ)卻依然是基于不考慮截面彎曲歷史的簡單彈塑性彎曲理論，從而將截面各次反向彎曲作為沒有關(guān)聯(lián)的、各自獨(dú)立的過程進(jìn)行分析，因而造成理論分析與現(xiàn)場應(yīng)用的較大差距[2,6-7]。使得對矯直過程現(xiàn)實(shí)問題的有效分析手段多為各種數(shù)值方法[8-14]。彈塑性彎曲是一種典型的截面不均勻變形狀態(tài)，會(huì)使截面產(chǎn)生宏觀殘余應(yīng)力[15]，所以輥式矯直過程中截面多次彈塑性反彎是一個(gè)不斷產(chǎn)生新的殘余應(yīng)力而同時(shí)截面應(yīng)力又在不斷進(jìn)行重分布的演變過程，在此過程中前次彎曲產(chǎn)生的殘余應(yīng)力遺傳至下一彎曲過程又必將會(huì)對下一彎曲過程的反彎特性產(chǎn)生影響。所以，在金屬條材彈塑性連續(xù)反彎過程中存在一種截面彎曲特性與應(yīng)力分布相耦合的演變過程，這時(shí)經(jīng)典彈塑性彎曲理論已經(jīng)不再適用，彎曲參數(shù)與力學(xué)參數(shù)間簡單的函數(shù)關(guān)系也將不復(fù)存在。因此，要建立合理的輥式矯直過程彎曲力學(xué)模型、精確分析多次反彎的截面彎曲特性，必須建立考慮遺傳變形歷史的截面彈塑性彎曲理論。本文作者以無量綱化的矩形截面為例，基于彈塑性彎曲的基本理論，建立考慮截面反彎特性與應(yīng)力遺傳關(guān)系的彈塑性二次反彎過程力學(xué)模型體系，并對彎曲歷史對截面應(yīng)力分布及彎曲特性的影響進(jìn)行研究，在此基礎(chǔ)上，揭示輥式矯直過程截面彎曲特性的復(fù)雜函數(shù)關(guān)系，并對進(jìn)一步的輥式矯直過程截面更多次反彎的參數(shù)關(guān)系解析函數(shù)形式進(jìn)行論述。

1 矩形截面的彎曲特性

設(shè)矩形截面寬度為B，高度為H，彈性極限彎矩為Mt，彈性極限曲率角為At。彎曲時(shí)，截面彈性區(qū)高度為Ht，截面彎矩為M。

對矩形截面定義如下4個(gè)參數(shù)：

由此，則實(shí)現(xiàn)了整個(gè)彎曲過程曲率、彎矩等各關(guān)系的無量綱化，無量綱參數(shù)ξ，C和是截面不同參數(shù)相對其彈性極限狀態(tài)的比值。由截面高度為±1的無量綱化的矩形截面推導(dǎo)產(chǎn)生的各種關(guān)系適用于任何理想金屬的矩形截面。

2 彎曲應(yīng)力與反彎特性耦合解析

2.1 簡單彈塑性彎曲

若完全不考慮(或不存在)截面的彈塑性彎曲歷程及由其產(chǎn)生的殘余應(yīng)力，則在此狀態(tài)下進(jìn)行的彈塑性彎曲可定義為簡單彈塑性彎曲過程。而金屬條材的輥式矯直截面第一次彎曲過程即屬于此簡單彈塑性彎曲狀態(tài)。

經(jīng)典的彈塑性彎曲理論中彈塑性彎曲和彈復(fù)過程如圖1所示，其截面應(yīng)力分布可分為加載應(yīng)力和卸載后殘余應(yīng)力2種狀態(tài)[16]。設(shè)截面屈服極限為σs，彎曲曲率比為C，則沿矩形截面高度z方向彎曲加載過程應(yīng)力分布方程為

對應(yīng)力分布方程進(jìn)行積分，得彎矩與曲率比的關(guān)系(即M-C關(guān)系)為

卸載后截面彎曲殘余應(yīng)力分布方程為

圖1 簡單彈塑性彎曲截面應(yīng)力演變過程Fig.1 Evolution of section stress during simple elastic-plastic bending process

2.2 二次彈塑性反彎

復(fù)雜彈塑性彎曲即是彎曲截面在簡單彈塑性彎曲基礎(chǔ)上，再經(jīng)過了一次或多次彈塑性彎曲的過程。在復(fù)雜彈塑性彎曲過程中，由于殘余應(yīng)力的存在和作用，截面的變形歷史開始影響彎曲過程中的曲率比與彎矩的關(guān)系。以輥式矯直截面第2次反彎過程為例進(jìn)行分析，其在力學(xué)模型上是最簡單的2次反彎過程，設(shè)截面第1次正向彎曲曲率比為CI，第2次反向彎曲曲率比為CII，現(xiàn)分析其第二次彎曲過程。

2.2.1 應(yīng)力分布

設(shè)第1次彎曲截面殘余應(yīng)力為σcI，第2次彎曲過程產(chǎn)生的彎曲應(yīng)力為σwⅡ，根據(jù)彎曲過程變形的幾何協(xié)調(diào)性和殘余應(yīng)力的疊加性，若不考慮應(yīng)力的屈服條件，其理論彎曲應(yīng)力的基本方程應(yīng)為

第2次彎曲單獨(dú)形成有應(yīng)力在彈性變形區(qū)應(yīng)滿足下式：

殘余應(yīng)力則滿足式(3)。

設(shè) ξII為截面二次彎曲后的彈塑性分界點(diǎn)，應(yīng)力滿足：

在ξII以外為塑性區(qū)應(yīng)力應(yīng)直接取為sσ。

由式(5)可以確定 ξII的位置。不過，由式(3)和圖1(b)的復(fù)雜性可知：ξII和ξI的相對位置(前后二次的彎曲曲率比)不同，截面的應(yīng)力分布不同，具體如圖2所示，所以式(5)所確定的ξII的具體公式也不一樣。

圖2 復(fù)雜彈塑性彎曲加載過程應(yīng)力分布情況Fig.2 Distribution of section stress during loading step of complicated elastic-plastic bending process

(1) ξII＞ξI。當(dāng)后一次彎曲程度相對較低，加載后不能統(tǒng)一前一次彎曲造成的殘余應(yīng)力疊加效應(yīng)時(shí)，截面應(yīng)力分布如圖2(a)所示，稱為1類應(yīng)力分布方式。將式(3)及式(4)帶入式(5)可得：

整理得彈區(qū)比ξⅡ：

在考慮截面應(yīng)力屈服條件的情況下，此時(shí)的截面加載后應(yīng)力分布方程如下：

將式(7)帶入上式并整理得復(fù)雜彈塑性彎曲過程的1類彎矩與曲率比關(guān)系方程(即M-C關(guān)系)如下：

(2) ξII＜ξI。當(dāng)后一次彎曲程度相對較高，加載后完全消除前一次彎曲造成的殘余應(yīng)力效應(yīng)時(shí)，截面應(yīng)力分布如圖2(b)所示，稱為第2類應(yīng)力分布方式。與第1類應(yīng)力分布推導(dǎo)相同，可以獲得：

2.2.2 應(yīng)力分布形式判別

由以上推導(dǎo)可知：復(fù)雜彎曲加載過程的應(yīng)力分布形式與CI和CII的組合關(guān)系有關(guān)。為求得應(yīng)力分布方式的判定方法，可取兩分布方程都適用的邊界情況，則：

整理得兩類應(yīng)力分布邊界函數(shù)為：

實(shí)際上，式(12)就是 ξII=ξI的關(guān)系式。

另外，除上述的彈塑性變形外，第二次彎曲也有可能存在彈性變形的情況。其極限情況是僅在邊部出現(xiàn)塑性應(yīng)力。此時(shí)由式(6)確定的ξII應(yīng)落在邊部，即為1。相應(yīng)的彎曲曲率比即為復(fù)雜彈塑性彎曲狀態(tài)下截面的彈性極限彎曲曲率比，即

由式(13)可以看出：復(fù)雜彈塑性彎曲狀態(tài)下，截面彈性極限彎曲曲率比CtII已經(jīng)不再等于 1，而是 CI的函數(shù)，即與截面的變形歷史有關(guān)。正常情況下， CI＞1，因此CtII＜1，即需要較小的力矩即可使斷面產(chǎn)生塑性變形。

將式(12)和式(13)表示的曲線畫在(CI，CII)平面上，可以得到復(fù)雜彈塑性彎曲過程的應(yīng)力狀態(tài)分布情況(圖3)。上述2條曲線將平面分為彈性變形、1類應(yīng)力和2類應(yīng)力3個(gè)區(qū)域，具體情況則與CI和CII2個(gè)參數(shù)的搭配相關(guān)，即復(fù)雜彈塑性彎曲的應(yīng)力分布由2個(gè)彎曲參數(shù)共同決定。

由此可見：在輥式矯直截面二次反彎過程中，代表截面彎曲歷史的參數(shù) CI對后次彎曲存在遺傳影響作用，不考慮截面變形歷史的彈塑性彎曲理論已經(jīng)無法精確描述截面的彎曲性質(zhì)。

圖3 復(fù)雜彎曲過程在(CI，CII )平面上的應(yīng)力分布Fig.3 Variation of stress distribution with bending parameter combination (CI, CII)

2.3 高次復(fù)雜彈塑性彎曲

在輥式矯直的實(shí)際過程中，金屬條材各截面最終都將經(jīng)歷3次以上的彈塑性連續(xù)反彎過程。由以上二次反彎的彈塑性解析可知：隨著截面反彎次數(shù)的不斷增多，其應(yīng)力遺傳耦合關(guān)系將會(huì)越來越復(fù)雜。

通過對二次反彎的推導(dǎo)，可以對更高次的復(fù)雜彈塑性彎曲的截面彎曲特性函數(shù)進(jìn)行如下描述：若設(shè)輥式矯直過程截面經(jīng)歷n次彎曲，則其截面變形歷史將由CI，CII，…，Cn共 n個(gè)彎曲參數(shù)進(jìn)行定義，其應(yīng)力分布方程可能是一個(gè)n元2n+1段的分段函數(shù)，其應(yīng)力分布特征需采用n維空間方可進(jìn)行描述。

可見：采用解析方法對輥式矯直高次的復(fù)雜彈塑性彎曲過程進(jìn)行徹底的解析較為困難，有賴于采用合理的數(shù)值方法進(jìn)行分析。

3 二次反彎過程的彎曲特性分析

3.1 加載過程M-C曲線

加載過程中的 M-C曲線集中反映著彈塑性彎曲加載過程截面的彎曲特性，對于實(shí)際工程應(yīng)用中合理地設(shè)定彎曲參數(shù)具有重要指導(dǎo)意義。

由以上理論公式(8)和(11)可以得到二次反彎加載過程中彎曲特性(即M-C曲線)，如圖4所示。

圖4 二次反彎過程M-C曲線Fig.4 M-C curves during the second reverse bending process

由圖4的M-C曲線可以看出：CI=1時(shí)，前次彎曲并沒有進(jìn)入彈塑性彎曲狀態(tài)，二次反彎過程的M-C曲線與簡單彈塑性彎曲過程相同。對比CI=1的曲線，前次簡單彎曲若進(jìn)入彈塑性彎曲狀態(tài)，則二次反彎過程截面彈性極限彎矩將按照式(13)的規(guī)律減小，相同曲率比條件下截面所需彎矩更小，即截面更易于發(fā)生彈塑性彎曲。同時(shí)，隨著前次彎曲曲率比的增大，截面二次反彎彎矩減小的趨勢將放緩。而 M-C關(guān)系在整個(gè)(CI，CII)平面上的曲面及等值線投影見圖5。

由以上對復(fù)雜反彎加載過程進(jìn)行的分析可以發(fā)現(xiàn)，經(jīng)歷二次反彎的截面更易于進(jìn)行彎曲。這也就是說，在實(shí)際輥式矯直的彈塑性彎曲過程中，矯直機(jī)第3輥(即截面第 2次彎曲位置)相較于經(jīng)典理論應(yīng)減小彎矩水平，即適當(dāng)減小壓下量，才能達(dá)到截面理想的彎曲曲率要求。

圖5 二次反彎過程M-C關(guān)系在(CI，CII)平面上的曲面Fig.5 M-C relation surface on (CI, CII) coordinate plane during the second reverse bending process

3.2 卸載過程應(yīng)力分布及回彈特性

由于二次反彎加載彈復(fù)后應(yīng)力分布函數(shù)σ及M-C關(guān)系函數(shù)極為復(fù)雜，因此不再列出完整的二次反彎殘余應(yīng)力分布函數(shù)。

以sσ=235 MPa為例，當(dāng)CII=3時(shí)，卸載后截面的殘余應(yīng)力分布狀態(tài)隨CI的變化如圖6所示。由圖6可以看出：前次彈塑性彎曲歷史對截面中心部位的殘余應(yīng)力分布形式具有很大的影響，而在連續(xù)反彎過程中截面邊部的殘余應(yīng)力水平則與本次反彎參數(shù)相關(guān)較大，邊部殘余應(yīng)力處在不斷變化的狀態(tài)，對經(jīng)歷輥式矯直過程多次彎曲的截面殘余應(yīng)力狀態(tài)的分析有賴于其整個(gè)變形歷史的彎曲參數(shù)。

圖6 CII=3的二次反彎卸載后截面殘余應(yīng)力分布Fig.6 Distribution of section residual stress after unloading step of the second reverse bending process for CII=3

同樣，由于彈復(fù)過程的純彈性效應(yīng)，彈復(fù)后截面曲率比Cf符合下式：

回彈比是考察截面回彈特性的重要參照量。該量在板條彎曲加工中用于判斷截面是否能更好地保持受彎時(shí)的形狀[1]。在輥式矯直過程中也是對條材殘留曲率比計(jì)算的重要指標(biāo)。二次反彎過程中回彈比隨前次彎曲參數(shù)的變化如圖7所示。

圖7 二次反彎回彈比曲線Fig.7 Springback ratio curve during the second reverse bending process

由圖7可以看出：前次彎曲歷史的存在會(huì)使得回彈比增大，而且隨著前次彎曲曲率比的增大回彈比也呈增大趨勢，即彎曲曲率比CII與彈復(fù)后曲率比Cf更加接近。這就說明經(jīng)歷二次反彎的截面卸載時(shí)能更好地保持受彎時(shí)的形狀。因此對于輥式矯直過程而言，經(jīng)歷前次彎曲的金屬條材在矯直機(jī)第3輥的2次反彎過程中需要相對于經(jīng)典公式減小壓下量以利于其殘余曲率比的控制。

4 結(jié)論

(1) 理論解析證明輥式矯直多次彈塑性反彎過程中截面存在著與反彎歷史相關(guān)的、復(fù)雜的應(yīng)力分布與反彎特性關(guān)系。

(2) 復(fù)雜彈塑性彎曲的彎曲特性參數(shù)均是包含截面整個(gè)彎曲歷史參數(shù)的復(fù)雜函數(shù)。

(3) 輥式矯直過程經(jīng)歷二次反彎的條材截面彈性極限彎矩減小，截面更易于彎曲。

(4) 輥式矯直過程經(jīng)歷復(fù)雜反彎的截面殘余應(yīng)力狀態(tài)有賴于其整個(gè)變形歷史的彎曲參數(shù)。

(5) 經(jīng)歷二次反彎的截面回彈比增大，截面能更好地保持受彎時(shí)的形狀，應(yīng)有針對性地相對于經(jīng)典公式減小壓下量，控制殘余曲率比。

[1]余同希, 章亮熾. 塑性彎曲理論及其應(yīng)用[M]. 北京: 科學(xué)出版社, 1992: 2-126.YU Tong-xi, ZHANG Liang-chi. Plastic bending theory and its application[M]. Beijing: Science Press, 1992: 2-126.

[2]崔甫. 矯直原理與矯直機(jī)械[M]. 2版. 北京: 冶金工業(yè)出版社,2002: 157.CUI Fu. Straightening and straightening machine[M]. 2nd ed.Beijing: Metallurgical Industry Press, 2002: 157.

[3]Yu T X, Johnson W. Influence of axial force on the elastic-plastic bending and springback of a beam[J]. Journal of Mechanical Working Technology, 1981(6): 5-12.

[4]Schleinzer G, Fischer F D. Residual stress formation during the roller straightening of railway rails[J]. International Journal of Mechanical Sciences, 2001, 43: 2281-2295.

[5]Johnson W, Yu T X. On the range of applicability of results for the springback of an elastic/perfectly plastic rectangular plate after subjecting it to biaxial pure bending[J]. Int J Mech Sci,1981, 23: 631-637.

[6]臧勇, 王會(huì)剛, 崔福龍. 型鋼輥式矯直壓彎撓度的彈塑性解析[J]. 機(jī)械工程學(xué)報(bào), 2005, 41(11): 47-51.ZANG Yong, WANG Hui-gang, CUI Fu-long. elastic-plastic analyses of bending deflection on section roller straightening[J].Chinese Journal of Mechanical Engineering, 2005, 41(11):47-51.

[7]周存龍, 王國棟, 謝東鋼, 等. 輥式矯直過程中板帶彎曲撓度的確定[J]. 太原科技大學(xué)學(xué)報(bào), 2009, 30(1): 48-50.ZHOU Cun-long, WANG Guo-dong, XIE Dong-gang, et al.Effect of entrance/exit leveler roller intermesh to plate flatness[J].Journal of Taiyuan University of Science and Technology, 2009,30(1): 48-50.

[8]薛軍安, 崔麗, 胡賢磊, 等. 輥式矯直過程的接觸傾角與曲率耦合分析[J]. 中國冶金, 2009, 19(2): 23-26.XUE Jun-an, CUI Li, HU Xian-lei, et al. Analysis of contact angle and curvature during roller leveling process[J]. China Metallurgy, 2009, 19(2): 23-26.

[9]周劍華, 吳迪, 趙憲明, 等. 輥式水平矯直對重軌斷面尺寸的影響分析[J]. 鋼鐵, 2009, 44(2): 40-43.ZHOU Jian-hua, WU Di, ZHAO Xian-ming, et al. Effect of horizontal roll straightening on cross section dimensions of heavy rail[J]. Iron & Steel, 2009, 44(2): 40-43.

[10]崔麗紅, 臧勇, 章博, 等. 多輥矯直過程中H型鋼斷面的應(yīng)力演變規(guī)律[J]. 北京科技大學(xué)學(xué)報(bào), 2008, 30(8): 942-946.CUI Li-hong, ZANG Yong, ZHANG Bo, et al. Stress evolvement rule of an H-beam section during multi-roller straightening process[J]. Journal of University of Science and Technology Beijing, 2008, 30(8): 942-946.

[11]周文, 劉學(xué)毅. 高速道岔尖軌矯直的有限元分析[J]. 西南交通大學(xué)學(xué)報(bào), 2008, 43(1): 82-95.ZHOU Wen, LIU Xue-yi. FEM simulation of straightening tongue rail of high-speed turnout[J]. Journal of Southwest Jiaotong University, 2008, 43(1): 82-95.

[12]王福華, 劉安中, 李友榮, 等. 輥式矯直機(jī)上輥壓下規(guī)程的仿真分析[J]. 武漢科技大學(xué)學(xué)報(bào), 2009, 32(4): 351-354.WANG Fu-hua, LIU An-zhong, LI You-rong, et al. Simulation analysis of upper rolling schedule of roll-type straightener[J].Journal of Wuhan University of Science and Technology, 2009,32(4): 351-354.

[13]高燕, 臧勇. 輥式矯直中 H 型鋼斷面畸變的仿真分析[J]. 北京科技大學(xué)學(xué)報(bào), 2006, 28(12): 1157-1161.GAO Yan, ZANG Yong. Finite element simulation of section deflection during H-beam roller straightening[J]. Journal of University of Science and Technology Beijing, 2006, 28(12):1157-1161.

[14]王會(huì)剛, 劉學(xué)江, 劉炳新, 等. 有限差分法在 H 型鋼輥式矯直壓下?lián)隙扔?jì)算中的應(yīng)用[J]. 鍛壓技術(shù), 2005(6): 41-43.WANG Hui-gang, LIU Xue-jiang, LIU Bing-xin, et al. Application of finite difference calculus in bending deflection calculation of H-beam roller straightenning[J]. Forging & Stamping Technology,2005(6): 41-43.

[15]米谷茂. 殘余應(yīng)力的產(chǎn)生和對策[M]. 朱荊璞, 邵會(huì)孟, 譯.北京: 機(jī)械工業(yè)出版社, 1983: 10.Yonetani Shigeru. The engender theory and countermeasure of residual stress[M]. ZHU Jing-pu, SHAO Hui-meng, transl.Beijing: China Machine Press, 1983: 10.

[16]徐秉業(yè), 劉信聲. 應(yīng)用彈塑性力學(xué)[M]. 北京: 清華大學(xué)出版社, 1995: 50-51.XU Bing-ye, LIU Xin-sheng. Applied elastic-plastic mechanics[M]. Beijing: Tsinghua University Press, 1995: 50-51.