摘 要：建立板材激光彎曲的三維非線性準(zhǔn)靜態(tài)彈塑性熱力耦合有限元模型。使用有限元軟件MSC.Marc對(duì)中厚船舶鋼板的激光彎曲成形過(guò)程進(jìn)行數(shù)值模擬。計(jì)算了船舶鋼板激光彎曲成形過(guò)程的溫度場(chǎng)和變形場(chǎng)，并進(jìn)行相應(yīng)的實(shí)驗(yàn)驗(yàn)證。模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合較好。用建立的模型對(duì)中厚船舶鋼板的激光彎曲成形過(guò)程中鋼板的幾何效應(yīng)進(jìn)行數(shù)值模擬，討論了一定工藝條件下鋼板幾何參數(shù)與彎曲角度之間的關(guān)系，為在將來(lái)實(shí)際生產(chǎn)中制定和優(yōu)化鋼板激光彎曲成形的工藝參數(shù)提供理論依據(jù)。 關(guān)鍵詞：船舶鋼板；激光彎曲；幾何效應(yīng)；數(shù)值模擬；有限元 中圖分類號(hào)：TG365.25 1 引 言 板材激光彎曲成形是近些年提出的一種新的板材塑性成形的加工工藝方法。激光彎曲成形利用激光束掃描金屬板材誘發(fā)的內(nèi)部非均勻分布的熱應(yīng)力，使板材發(fā)生局部塑性屈服，從而產(chǎn)生彎曲變形[1]。 早期的研究主要采用傳統(tǒng)的實(shí)驗(yàn)研究方法，這種方法證據(jù)充分、說(shuō)服力強(qiáng)，針對(duì)影響因素能快速直觀地得出結(jié)論。許多學(xué)者采用實(shí)驗(yàn)方法研究了板材激光彎曲成形過(guò)程，并得出一些結(jié)論，總結(jié)了一些規(guī)律[2-8]。但是，由于板材激光彎曲成形過(guò)程是一個(gè)復(fù)雜的熱彈塑性變形過(guò)程，其影響因素眾多，如：激光功率、光斑大? ⑸杪肪丁⑺俁群痛問(wèn)⒒蝗忍跫⒓諧痔跫宀牡募負(fù)緯嘰繅約叭任鐨圓問(wèn)土ρ閱懿問(wèn)齲俠砩杓坪陀嘔す饌淝尚喂ひ詹問(wèn)檬笛櫸椒ń斐紗罅康氖笛橄摹６捎檬的Ｄ餳際酰岷仙倭康氖笛檠櫓ぃ詡撲慊賢瓿杉す饌淝尚喂痰姆掄媯槐鼗ǚ汛罅康娜肆ξ锪Χ允滴锝釁飾觶行Ъ跎倭聳笛橄摹Ｒ虼耍的Ｄ夥椒ㄒ殉晌昀窗宀募す饌淝尚喂痰鬧饕芯渴侄蝃9],國(guó)內(nèi)外學(xué)者對(duì)板材激光彎曲成形已進(jìn)行了大量的數(shù)值模擬研究［10-18］。 板材激光直線掃描彎曲成形過(guò)程中，板料的幾何尺寸對(duì)成形有比較大的影響。板料的幾何尺寸主要包括板料的長(zhǎng)度L（垂直于激光掃描路徑方向的尺寸）或L＇（自由端距激光掃描路徑的距離）、寬度b（沿激光掃描路徑方向的尺寸）和厚度t。盡管有許多學(xué)者對(duì)一些材料的板材激光彎曲成形過(guò)程中板材的幾何尺寸效應(yīng)開(kāi)展了一些實(shí)驗(yàn)[4-8]和數(shù)值模擬研究[14-18]，但基本上都是針對(duì)薄板，對(duì)厚度超過(guò)3 mm的板材激光成形的幾何效應(yīng)的數(shù)值模擬的報(bào)道還很少。本文在驗(yàn)證了激光彎曲成形數(shù)值模擬結(jié)果準(zhǔn)確性的基礎(chǔ)上，由模擬結(jié)果結(jié)合一定的實(shí)驗(yàn)數(shù)據(jù)，研究中厚船舶鋼板激光彎曲成形過(guò)程中板材的幾何尺寸與彎曲角度之間的關(guān)系，討論這些參數(shù)對(duì)激光彎曲成形規(guī)律的影響，為制定激光彎曲成形的合理的工藝參數(shù)組合提供依據(jù)。 2 板材激光彎曲成形模型的建立 2.1 彈塑性有限元熱力耦合基本方程 有限變形彈塑性有限元方程的矩陣形式可描述為 式中：[K]為切線剛度矩陣，{△u}為節(jié)點(diǎn)位移增量，為等效外載荷，{Fnr}為Newton-Raphson失衡力。 假設(shè)環(huán)境溫度恒定，考慮板材與周圍空氣的對(duì)流換熱，溫度計(jì)算的有限元基本方程為式中：[C]為熱容矩陣，[Kt]為熱傳導(dǎo)矩陣，{T}為節(jié)點(diǎn)溫度，為節(jié)點(diǎn)溫度對(duì)時(shí)間的導(dǎo)數(shù)，{Q}為熱流矢量。 聯(lián)立有限變形彈塑性有限元方程（1）和溫度計(jì)算的有限元方程（2），得到熱力耦合計(jì)算的基本方程為式中{F}是力矢量，包括節(jié)點(diǎn)力矢量和熱應(yīng)變引起的力矢量。 2.2 有限元模型的建立 基于有限元軟件MSC.Marc建立板材激光彎曲成形過(guò)程三維彈塑性熱力耦合計(jì)算模型，考慮了材料的熱物理性能和力學(xué)性能參數(shù)與溫度的非線性關(guān)系。模型中采用不同尺寸的試板，材料為船舶用低碳鋼，建立的三維熱力耦合模型如圖1所示。激光沿板材寬度方向直線勻速掃描，由于溫度和變形集中在掃描線附近，故在掃描線附近區(qū)域加密網(wǎng)格，遠(yuǎn)離掃描線的區(qū)域保持稀疏的網(wǎng)格。激光參數(shù)為功率P=1000 W，光斑大小為8×8 mm2，掃描速度為7.5 mm/s。材料對(duì)激光的吸收系數(shù)為0.65。 (1) 熱源 激光束作為外加表面熱流載荷輸入。通過(guò)用戶子程序描述激光束的連續(xù)掃描，使激光熱載荷以熱流密度的形式施加到相應(yīng)單元上。采用熱流均勻分布的方形光斑，其能量密度分布I為式中A為材料對(duì)激光的吸收系數(shù)，P為激光束的輸出功率，d為激光光斑寬度。 (2) 邊界條件 固定激光源，板沿Z的負(fù)方向以勻速v平動(dòng)。約束板材一端單元的所有節(jié)點(diǎn)X和Y方向的自由度，以避免板材發(fā)生剛體位移，這與板材激光彎曲成形過(guò)程的實(shí)驗(yàn)一致。板材經(jīng)激光束掃描后在空氣中自然冷卻，只考慮板材表面與空氣的對(duì)流換熱，而不考慮熱輻射。對(duì)流換熱邊界條件為式中h是表面對(duì)流換熱系數(shù)。 (3) 初始條件(4) 材料熱物性參數(shù)材料的熱物性參數(shù)和力學(xué)性能參數(shù)均為溫度的函數(shù)。對(duì)未知溫度范圍內(nèi)的材料性能按線性插值法確定。 3 模擬與實(shí)驗(yàn)結(jié)果的比較 通過(guò)實(shí)驗(yàn)來(lái)驗(yàn)證模擬結(jié)果的有效性，實(shí)驗(yàn)用船舶鋼板的尺寸為：長(zhǎng)L=300 mm，寬b=150 mm，厚t分別為4、5、6 mm，采用激光位移傳感器測(cè)量了下表面距自由端5 mm處一點(diǎn)Y向（即厚度方向）位移隨時(shí)間的變化。圖2(a)和(b)分別表示厚度為4 mm和6 mm鋼板激光彎曲Y向位移的模擬和實(shí)驗(yàn)之間的比較。 圖2表明，變形的模擬和實(shí)驗(yàn)在彎曲的整個(gè)過(guò)程都吻合得很好，最終Y向位移的模擬值分別為2.220 mm和0.979 mm，都經(jīng)過(guò)三次實(shí)驗(yàn)驗(yàn)證，測(cè)量的平均值分別為2.222 mm和1.027 mm，計(jì)算誤分別為-0.09%和-4.67%。通過(guò)實(shí)驗(yàn)證明，所建立的模型比較精確。 4 模擬結(jié)果與分析 4.1板材厚度對(duì)彎曲變 板材在厚度方向上所形成的溫度梯度是鋼板產(chǎn)生彎曲變形的根本原因，因此板材的厚度對(duì)激光彎曲變尤為突出。圖3表示在圖示工藝參數(shù)條件下，鋼板的厚度與角度的關(guān)系。從圖中可以看到，在該工藝條件下，隨著板厚的增加，彎曲角度先增大，隨后迅速減小，當(dāng)板厚大于激光光斑寬度時(shí)，彎曲角度減小速度減緩，這時(shí)的彎曲角度已經(jīng)很? Ｓ紗絲芍卑搴癯臣拗凳保薹ń屑す獬尚巍Ｊ導(dǎo)噬希宀腦膠瘢莧惹奈露確逯翟降停擾蛘土吭叫。捎諼律鸕那αΦ南陸抵翟叫。蓖淝璧耐淝匾蒼醬螅宀暮穸仍黽釉斐砂宀牡哪謨?chǔ)力减小，两烦C嫻淖酆獻(xiàn)饔檬垢職逅芐越嵌妊桿偌跣 ５卑宀暮穸刃∮詮獍嚦磯鵲囊話耄絳跣“宀暮穸齲牧轄⑸芻蟶賬穡淝嵌確炊跣 Ｋ裕雜謖庖蛔楣ひ詹問(wèn)宀暮穸仍詡す夤獍嚦磯鵲?.5倍～1倍之間，可以產(chǎn)生比較大的彎曲角度。4.2 板材寬度對(duì)彎曲變 圖4表示板材在長(zhǎng)度和厚度尺寸一定的條件下，不同的寬度對(duì)激光彎曲角度的影響。從圖可以看出，板材的寬度對(duì)彎曲變較大。隨鋼板寬度增加，彎曲變形最開(kāi)始增加很快，而后逐漸變緩，到達(dá)一定寬度后，變形基本不再增加。這是因?yàn)榧す饧訜崾茄刂宀牡膶挾确较蜻B續(xù)進(jìn)行的，通常激光束的直徑相對(duì)板寬很小，使得同一時(shí)刻被加熱的材料也很少。在特定的工藝下被加熱處產(chǎn)生塑性變形時(shí)，其它在寬度方向上沒(méi)加熱到的材料阻礙加熱區(qū)域材料的變形，具有一定的冷端抑制作用。鋼板越寬，冷端抑制作用越明顯。又因?yàn)槔涠藢?duì)加熱過(guò)程中彎曲的阻礙作用大于對(duì)冷卻過(guò)程中的正向彎曲的阻礙作用，所以鋼板越寬，背向激光源的反向彎曲角度越小，而朝向激光源的最終彎曲角度越大。但是當(dāng)板材的寬度達(dá)到某一值后，冷端抑制作用達(dá)到最大，使曲線逐漸趨于平緩。4.3 板材長(zhǎng)度對(duì)彎曲變 在激光彎曲成形中，板材一端夾持，夾持端的長(zhǎng)度對(duì)彎曲成形影響不大，而激光掃描線距板材自由端的距離L＇對(duì)激光彎曲角度有一定的影響。圖5為在一定的工藝條件下，當(dāng)鋼板厚度和寬度尺寸不變時(shí)，激光掃描線距板材自由端的距離對(duì)激光彎曲成形角度的影響。圖中的曲線表明，激光掃描線到自由端的距離越大，所獲得的彎曲角度越大。和寬度對(duì)成差不多一致，激光掃描線距自由端的距離越大，自由端對(duì)它的約束作用就越大，彎曲角度也越大。但是，當(dāng)L＇達(dá)到某一長(zhǎng)度之后，自由端對(duì)它的約束作用不再增加，曲線逐漸趨于平緩。 5 結(jié) 論 建立了船舶鋼板激光彎曲成形的三維非線性準(zhǔn)靜態(tài)彈塑性熱力耦合有限元模型，計(jì)算了船舶鋼板激光彎曲成形過(guò)程的溫度場(chǎng)和變形場(chǎng)，通過(guò)實(shí)驗(yàn)驗(yàn)證了模型的有效性，并用建立的模型對(duì)中厚船舶鋼板的激光彎曲成形過(guò)程中鋼板的幾何效應(yīng)進(jìn)行了數(shù)值模擬研究，得出以下結(jié)論： 1）模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合較好。 2）一定的工藝條件下，彎曲角度隨鋼板厚度的增加，先增加，隨后迅速減小，最后逐漸變緩，超過(guò)一定極值后，不能發(fā)生彎曲變形；對(duì)于本文的激光工藝參數(shù)組合，板材厚度在激光光斑寬度的0.5倍～1倍之間，可以產(chǎn)生合適的彎曲角度。 3）鋼板彎曲角度隨鋼板的寬度和激光掃描線到鋼板自由端的距離的增加先迅速增加，隨后逐漸變緩，超過(guò)一定值后，彎曲角度不再變化。 參考文獻(xiàn) [1] M. Geiger, F. Vollertsen. The Mechanisms of Laser Forming. Annals of the CIRP, 1993,42(1):301-304. [2] Gareth Thomson, Mark Pridham. Material property changes associated with laser forming of mild steel components. Journal of Materials Processing Technology, 2001,118:40-44 [3] M. Merklein, T. Hennige, M. Geiger. Laser forming of aluminium and aluminium alloys—microstructural investigation. Journal of Materials Processing Technology, 2001, 115:159- 165 [4] K.C.Chan, J.Liang. Laser bending of An Al6013/SiCp aluminium matrix composite sheet. Journal of Materials Processing Technology, 2000,100:214-218. [5] Wu shichun, Zheng Jinsong. An experimental study of laser bending for sheet metals. Journal of Materials Processing Technology. 2001,110:160-163. [6] 管延錦，孫勝，季忠. 板料激光成型技術(shù)的實(shí)驗(yàn)研究. 光學(xué)技術(shù), 2000,26(3):260-262. [7] 陳敦軍，吳詩(shī)惇，向毅斌，李淼泉. 鈦合金板材激光彎曲成形的研究. 航空學(xué)報(bào)，2001,22(2):187-189. [8] 王秀鳳，王秀彥，林道盛. 板料激光彎曲的實(shí)驗(yàn)研究. 鍛壓機(jī)械，1999, 3:8-10. [9] 張立文，鐘琦，呂波，裴繼斌，張國(guó)梁. 板材激光彎曲成形數(shù)值模擬的研究現(xiàn)狀及展望. 金屬成形工藝，2003,21(1): 25-29 [10] An.K.Kyrsanidi, Th.B.Kermanidis, Sp.G.Pantelakis. Numerical andexperimental investigation of the laser forming process . Journal of Materials Processing Technology.1999,87:281- 290. [11] Z. Hu, R. Kovacevic, M. Labudovic. Experimental and numerical modeling of bucking instability of laser sheet forming. International journal of machine tools & manufacture. 2002,42:1427-1439. [12] Zhang Liwen, Zhong Qi, Pei Jibin, Zhang Guoliang, Xia Yuanliang. FEM Simulation of Laser Forming Process of Shipbuilding Steel Plate. 2nd International Conference Thermal Process Modeling and Computer Simulation. Nancy,France, 2003. [13] 季忠,吳詩(shī)惇. 板料激光成形時(shí)的形變場(chǎng)研究. 塑性工程學(xué)報(bào),1998,5(2):33-38 [14] An.K.Kyrsanidi, Th.B.Kermanidis,Sp.G.Pantelakis. An analytical model for the prediction of distortions caused by the laser forming process. Journal Materials Processing Technology, 2000,104:94-102. [15] Zhong Ji, Shichun Wu. FEM simulation of the temperature field during the laser forming of sheet metal. Journal of Materials Processing Technology, 1998,74: 89-95. [16] 李緯民, M.Geiger, F.Vollertsen. 金屬板材激光彎曲成形規(guī)律的研究. 中國(guó)激光, 1998, 25 (9): 859-864. [17] 李緯民，盧秀春，劉助柏. 激光彎曲工藝中板材厚度的影響規(guī)律. 中國(guó)有色金屬學(xué)報(bào), 1999,9(1):39-43. [18] 季忠，吳詩(shī)惇. 板料激光彎曲變形的幾何和能量效應(yīng)研究. 航空精密制造技術(shù), 1996, 32(4):21-23.