張子騫 楊會林 田永利

1.東北大學(xué)，沈陽，110004 2.天津大學(xué)，天津，300072

0 引言

近年來，高精度薄壁金屬管材在工業(yè)、軍事、醫(yī)療、航空航天、石油化工等領(lǐng)域中的需求急劇上升，為使管材具備理想的直線度和圓度，必須在出廠或使用前對其進(jìn)行精整矯直。目前，薄壁管材的矯直生產(chǎn)一般通過人工經(jīng)驗或反復(fù)試矯完成，生產(chǎn)效率較低、產(chǎn)品質(zhì)量穩(wěn)定性較差，且針對的理論研究尚少。基于經(jīng)典管材矯直理論確定的工藝參數(shù)在生產(chǎn)實踐中，表現(xiàn)出矯直精度低、截面畸變，甚至無法對工件有效矯直，并不適合薄壁管材的矯直生產(chǎn)，因此亟待對薄壁管材的矯直機(jī)理和變形過程進(jìn)行深入的理論研究。

薄壁管材的矯直是針對已彎曲的管材施加更大程度的反向彎曲來實現(xiàn)的，與管材彎曲理論相似，彎曲外側(cè)和內(nèi)側(cè)之間存在一個應(yīng)變中性層［1］。隨著彎曲變形程度的不斷增大，中性層將偏移管材中心層，對偏移量的解析及其對矯直過程影響的分析是后續(xù)深入研究薄壁管材矯直理論的基礎(chǔ)。目前，國內(nèi)外針對中性層偏移問題的處理有兩種方法:一種是基于梁彎曲的工程理論，忽略中性層的偏移［2］;一種是由管材精彎理論建立的偏移模型，其理論基礎(chǔ)為平面應(yīng)變條件下平板塑性彎曲的Hill理論［3-5］。這兩種方法在處理普通板材和管棒材的彎曲與矯直問題時具有較好的適用性，但不適合薄壁管材的矯直，同時這些模型對材料性能的變化不敏感。為此，本文針對薄壁管材的結(jié)構(gòu)特點，確定了平面變形的應(yīng)力應(yīng)變關(guān)系，基于此對應(yīng)變中性層的偏移量進(jìn)行了解析，分析了其變化趨勢和對矯直過程的影響，為深入研究薄壁管材的矯直機(jī)理和變形過程提供理論依據(jù)。

1 基本假設(shè)

在建立中性層偏移模型時需引入以下假設(shè):

(1)材料是均質(zhì)、連續(xù)及各向同性的，忽略矯直變形區(qū)的Baushinger效應(yīng)，即材料具有拉壓一致性，并認(rèn)為應(yīng)力中性層與應(yīng)變中性層重合。

(2)任意垂直于中心軸的管材截面在矯直過程中始終為平面，且與中心軸保持垂直;相鄰兩截面間無傾斜扭曲，并忽略管材截面尺寸的變化。

(3)薄壁管材矯直變形過程中的塑性區(qū)域遵循體積不變原則。

(4)由于管壁很薄，故忽略管材矯直過程中的壁厚變化，同時徑向正應(yīng)力遠(yuǎn)小于其他應(yīng)力，可忽略不計。

(6)矯直工作過程平穩(wěn)，管材在工作輥面上做純滾動，在矯直輥縫內(nèi)的變形區(qū)可視為等彎矩純彎曲過程［6］。

2 矯直變形分析

薄壁管材通常指管徑與壁厚之比遠(yuǎn)大于20的管材，其在矯直過程中易出現(xiàn)起皺、扭曲和破裂等缺陷，為避免缺陷產(chǎn)生，通常采用減小單次反彎變形量，增加反彎次數(shù)的方法，一般采用8輥或10輥斜輥矯直機(jī)。工作輥分為上輥與下輥，在空間上交錯布置于管材兩側(cè)，與被矯直管材成一定夾角，如圖1所示［6］。具有初始曲率的管材依次通過中間各組具有等曲率輥形的矯直輥，進(jìn)行連續(xù)彎曲變形實現(xiàn)矯直。

圖1 上下矯直工作輥空間布置

這里采用切塊法對薄壁管材在輥縫內(nèi)的矯直變形區(qū)進(jìn)行分析，與經(jīng)典管材彎曲理論在圓柱坐標(biāo)系下選取微元體的方法不同，管材在矯直時具有初始曲率，管殼表面為空間曲面，故這里選用空間子午坐標(biāo)系，在輥縫內(nèi)矯直變形區(qū)的任意位置取微元體單元ABCD，如圖2所示。該單元體由2個夾角為dθ的管材子午面與2個夾角為dφ的平行圓包圍而成，子午面與某一基準(zhǔn)子午面夾角為θ，平行圓處的管殼表面法線與彎曲軸線夾角為φ，因管壁較薄，取微元體的厚度為管壁厚度t。同時規(guī)定任意一點的坐標(biāo)方向如下:X軸方向為平行圓的切線方向，Y軸方向為管材子午線的切線方向，Z軸方向為管材外表面的法線方向。

薄壁管材的矯直變形與普通管材彎曲變形的相似之處:變形區(qū)以中性層為界，X方向外側(cè)材料主要受切向拉伸而伸長，內(nèi)側(cè)材料主要受切向壓縮而縮短。

圖2 薄壁管材矯直區(qū)應(yīng)力應(yīng)變狀態(tài)模型

薄壁管材的矯直變形與普通管材的彎曲變形也有很大不同:在Z方向上，普通管材彎曲理論認(rèn)為外側(cè)壁厚減小，內(nèi)側(cè)壁厚增大;薄壁管材因壁厚較小，矯直工藝本身要求矯直前后壁厚基本無變化，故矯直變形程度相對較小，因此矯直過程造成的厚度變化是十分微小的，由假設(shè)(4)知，Z方向上的變形可忽略不計，即εr=0;同時Z方向上的正應(yīng)力遠(yuǎn)小于其他應(yīng)力，有σr=0，故矯直區(qū)的三維應(yīng)變狀態(tài)簡化為平面應(yīng)變問題，則Z平面上的切應(yīng)力分量 τθr=τφr=0，切應(yīng)變分量 γθr=γφr=0。

薄壁管材的矯直與普通管材的彎曲的另一個重要區(qū)別是加工工藝不同。普通管材在彎曲時，沿軸線運動通過彎曲模，管材本身并不繞軸旋轉(zhuǎn)，因此分析時可將管材子午線切線方向即Y方向的應(yīng)力與應(yīng)變忽略;薄壁管材在圖1所示的空間交錯的工作輥中螺旋前進(jìn)，管材除了沿軸線直線運動外，還繞Y軸旋轉(zhuǎn)，再者因管材厚度小，Y方向上抵御失穩(wěn)的能力差，因此Y方向的應(yīng)力與應(yīng)變需要計入。同時，由于薄壁管材矯直變形程度相對較小，截面上橫剪力τθφ雖然存在，但由假設(shè)(2)、(6)可認(rèn)為τθφ是足夠小量，而由橫向剪應(yīng)變引起的橫截面任何變形可以忽略［7］，因此 τθφ=γθφ=0。

經(jīng)上述分析，矯直變形區(qū)的應(yīng)力應(yīng)變狀態(tài)與相關(guān)幾何關(guān)系如圖2所示。其中，r為薄壁管材半徑，Rw為薄壁管材幾何中心軸反彎半徑，管材彎曲矯直時，沿彎曲切線方向(X向)作用的應(yīng)力、應(yīng)變?yōu)棣姚群挺纽龋毓懿淖游缇€切線方向(Y向)作用的應(yīng)力、應(yīng)變?yōu)棣姚蘸挺纽眨瑥澢鷳?yīng)變中性層偏移角為α，t為管壁厚度。

3 矯直區(qū)彈塑性應(yīng)力應(yīng)變狀態(tài)

3.1 應(yīng)力平衡微分方程

對微元體單元ABCD列寫應(yīng)力平衡微分方程，在X方向上由于假設(shè)τθφ=0，因此

在Y方向上，考慮σφ的同時還要考慮σθ在Y方向上的投影，如圖3a所示，σθ在平行圓徑向的投影分量為σθrdφdθ，則該力在Y方向上的分量為 σθrdφdθcosφ，如圖3b所示。因此，對Y方向列寫應(yīng)力平衡微分方程，有

又因為 sin(φ+dφ)=sinφcosdφ+cosφsindφ ≈sinφ+dφcosφ，將其代入式(2)并化簡整理得

圖3 σθ的投影關(guān)系圖

Z方向上的正應(yīng)力σr=0，因此不考慮該方向的應(yīng)力平衡方程。

3.2 應(yīng)變分析

薄壁管材具有初始曲率，故可假設(shè)微元體ABCD初始處于壓縮區(qū)，則反彎后微元體處于拉伸區(qū)。由于管壁較薄，若想實現(xiàn)矯直其初始彎曲程度不能太大，可認(rèn)為矯直前薄壁管材的幾何中心軸與中性層重合，則管材的初始長度l0為

式中，R0為薄壁管材矯直前的初始曲率;θ0為薄壁管材矯直前的初始中心角。

微元體ABCD處初始纖維長度lφ0為

薄壁管材在矯直時需反向彎曲，中性層向彎曲軸向方向偏移α角，此時中性層的長度lw為

其中，θw為薄壁管材矯直時的反彎中心角，此時微元體ABCD處的纖維長度lφw為

因此，Y方向的真應(yīng)變?yōu)?/p>

變形前后中性層的長度保持不變，矯直區(qū)的中性層長度等于管材初始長度，即

因此

3.3 彈性區(qū)應(yīng)力應(yīng)變關(guān)系

為使薄壁管材得到矯直而截面形狀基本不變，必須保證中性層附近的某個區(qū)域為彈性變形區(qū)域。由于彈性變形區(qū)域較小，且彈性變形量通常較小，因此可認(rèn)為彈性區(qū)內(nèi)各變形層的初始長度均為管材的初始長度，則彈性區(qū)Y方向的真應(yīng)變可近似表示為

在彈性區(qū)域中，金屬遵循線彈性變形的基本原理，其本構(gòu)關(guān)系遵循廣義胡克定律:

其中，μ為泊松比，將式(10)代入式(9)得

將式(11)代入式(3)并整理得

該一階線性非齊次微分方程的通解為

將式(14)～式(16)代入式(13)并化簡整理得

將式(17)代入式(11)并化簡整理得

又由假設(shè)(1)可知，當(dāng)φ=α?xí)r，σθ=0，將其代入式(18)可得

將式(19)代入式(17)、式(18)，得彈性區(qū)X、Y方向的主應(yīng)力:

3.4 塑性區(qū)應(yīng)力應(yīng)變關(guān)系

薄壁管材在矯直過程中，除中性層附近的彈性區(qū)域外，其余均為塑性變形區(qū)域，遵循金屬塑性變形相關(guān)規(guī)律。由假設(shè)(3)知，在塑性變形區(qū)有εθ+εφ=0，則 εr=(εθ+εφ)/2。平面應(yīng)變狀態(tài)下，由增量理論知σr=(σθ+σφ)/2。因此，塑性區(qū)的等效應(yīng)力、等效應(yīng)變?yōu)?/p>

由假設(shè)(5)知

塑性變形相對于彈性變形而言，變形程度較大，故Y方向塑性應(yīng)變由式(8)表示。將式(8)、式(22)、式(23)代入式(24)并整理得

4 中性層偏移數(shù)學(xué)模型

薄壁管材在矯直過程中，彈性區(qū)域與塑性區(qū)域必然存在一個屈服分界面，當(dāng)管材的變形量增大到一定程度，材料發(fā)生初始屈服，由假設(shè)(5)可知，彈性區(qū)的等效應(yīng)力與等效應(yīng)變的關(guān)系為

當(dāng)彈性區(qū)外層發(fā)生初始屈服時，由Mises屈

服準(zhǔn)則有

式中，σs為屈服應(yīng)力。

將式(27)、式(26)代入式(28)得

若此時屈服面偏移中心軸的角度為β，將式(20)、式(21)代入式(29)得

因此

而對于塑性區(qū)的應(yīng)力應(yīng)變狀態(tài)，在初始屈服時同樣應(yīng)滿足式(28)，即

5 數(shù)學(xué)模型分析與仿真實驗

由式(33)可見，中心層偏移角α與薄壁管材的工藝結(jié)構(gòu)參數(shù)(初始相對彎曲半徑R0/r、反彎相對彎曲半徑Rw/r)有關(guān)。通常，矯直工藝將來料管材的R0/r最大值規(guī)定為定值，這樣α便隨Rw/r的變化而變化;同時，α還與材料的力學(xué)性能有關(guān)，式(33)中的量綱一參數(shù)S、T就是由材料的力學(xué)性能決定的。不同金屬的力學(xué)性能與結(jié)構(gòu)參數(shù)(薄壁管材半徑r=20mm，管材壁厚t=1mm)如表1所示［8-9］，計算了不同Rw/r下的偏移角值，如圖4的空心點所示;基于此擬合了α隨Rw/r的變化曲線，如圖4的實線所示。為驗證偏移模型的正確性采用ANSYS/LS－DYNA對典型材料的薄壁管材矯直變形過程進(jìn)行了動態(tài)仿真，測量了相同條件下的中心層偏移角，在圖4中以虛線表示。

表1 不同金屬力學(xué)性能與結(jié)構(gòu)參數(shù)表

圖4 不同典型金屬材料薄壁管材α隨R w/r的變化曲線

由圖4可見，中心層偏移角的仿真實驗值與理論計算值存在一定差距，這主要是由仿真過程未計入管材的初始變形、建立有限元模型時網(wǎng)格劃分較大、理論計算模型并未涉及約束形式等原因造成的，但兩者差距并不是很大，基本控制在5°以內(nèi)，同時仿真實驗值與計算值隨Rw/r的變化趨勢大體相同，因此可認(rèn)為式(33)的計算結(jié)果近似準(zhǔn)確。

由圖4可見，在R0/r一定的前提下，中心層偏移角α隨Rw/r的減小而增大，這說明當(dāng)初始彎曲程度相同時，對于同一管材半徑，矯直反彎半徑越小，管材的變形程度越大，中性層向內(nèi)偏移越大;對于同一反彎半徑，管材的半徑越大，中性層偏移越大，這些定性分析結(jié)果與文獻(xiàn)［3-4］得到的結(jié)果相似。

對比圖4a和圖4b可知，不同材料的薄壁管材在矯直時的偏移角相差很大，這說明α與材料的力學(xué)性能有關(guān)，材料的塑性越好，偏移程度越大，這是文獻(xiàn)［3-4］中模型所無法體現(xiàn)的。由圖4a可見，對于碳素結(jié)構(gòu)鋼、合金結(jié)構(gòu)鋼等材料的管材，由于材料塑性變形能力相對較弱，管材初始彎曲半徑較大，為防止管材在矯直過程中發(fā)生破壞，其反彎半徑也相對較大，因此中性層偏移較小，一般均在2mm以內(nèi)，故在分析計算時可認(rèn)為中性層不發(fā)生偏移。由圖4b可見，對于不銹鋼、鋁合金與銅合金等材料的管材，由于材料塑性變形能力相對較強(qiáng)，其初始彎曲半徑通常較小，若想實現(xiàn)矯直，其反彎半徑也相對較小，中性層偏移較大，故在分析計算時需考慮中性層的偏移。

上述分析僅針對普通規(guī)格薄壁管材的反彎矯直，并不適用于微小直徑及超大直徑薄壁管材的矯直與彈塑性彎曲。

6 結(jié)論

①本文建立的薄壁管材矯直過程的中性層偏移與相關(guān)理論模型，可對彈塑性變形區(qū)的應(yīng)力應(yīng)變狀態(tài)、初始屈服面與應(yīng)變中性層位置進(jìn)行理論估算與分析。仿真實驗證明，該理論模型是近似正確的，為進(jìn)一步深入研究薄壁管材矯直機(jī)理和變形過程提供了理論依據(jù)。②薄壁管材矯直時，中心層的偏移量不僅與工藝結(jié)構(gòu)參數(shù)有關(guān)，還與材料的力學(xué)性能有關(guān)，它隨相對反彎半徑的增大而減小，隨金屬塑性變形能力的增強(qiáng)而增大。③對于碳素結(jié)構(gòu)鋼、合金結(jié)構(gòu)鋼等材料的薄壁管材，矯直過程的中性層偏移量較小，在工程計算分析時可忽略不計;對于不銹鋼、鋁合金與銅合金等金屬材料，中性層偏移量較大，分析計算時必須考慮。④本文建立的理論模型及其方法可推廣應(yīng)用到薄壁型鋼的矯直及相關(guān)領(lǐng)域的分析研究中。

［1］Stechpwicz F.Bending with Upsetting of Copper Tube Elbows［J］.Journal of Materials Processing Technology，2000，100:236-240.

［2］崔甫.矯直原理與矯直機(jī)械［M］.北京:冶金工業(yè)出版社，2002.

［3］鄂大辛，郭學(xué)東，寧汝新.管材彎曲中應(yīng)變中性層位移的分析［J］.機(jī)械工程學(xué)報，2009，45(3):307-310.

E Daxin，Guo Xuedong，Ning Ruxin.Analysis of Strain Neutral Layer Displacement in Tube － bending［J］.Chinese Journal of Mechanical Engineering，2009，45(3):307-310.

［4］余同希，章亮熾.塑性彎曲理論及其應(yīng)用［M］.北京:科學(xué)出版社，1992.

［5］Li Longyuan，Kettle R.Nonlinear Bending Response and Buckling of Ring－stiffened Cylindrical Shells under Pure Bending［J］.International Journal of Solids and Structures，2002，39:765-781.

［6］張子騫，張柏森，楊會林，等.管棒材等曲率矯直力模型可視化設(shè)計［J］.東北大學(xué)學(xué)報(自然科學(xué)版)，2012，33(3):409-413.

Zhang Ziqian，Zhang Baisen，Yang Huilin，et al.The Visual Design on the Force Model of Tubes and Bars Equal Curvature Straightening［J］.Journal of Northeastern University(Natural Science)，2012，33(3):409-413.

［7］王俊奎，張志民.鈑殼的彎曲與穩(wěn)定［M］.北京:國防工業(yè)出版社，1980.

［8］成大先.機(jī)械設(shè)計手冊［M］.5版.北京:化學(xué)工業(yè)出版社，2009.

［9］劉鳴放，劉勝新.金屬材料力學(xué)性能手冊［M］.北京:機(jī)械工業(yè)出版社，2010.