李欣濤，龍連春

(北京工業(yè)大學(xué)機(jī)械工程與應(yīng)用電子技術(shù)學(xué)院，北京 100124)

1 引 言

板殼結(jié)構(gòu)在面內(nèi)壓力作用下的屈曲破壞是常見的破壞形式。實(shí)際結(jié)構(gòu)有時(shí)會出現(xiàn)整體或局部溫度升高的情況。由于熱膨脹產(chǎn)生熱應(yīng)力，使材料的力學(xué)性能發(fā)生變化，因此溫度的變化會影響結(jié)構(gòu)的屈曲承載能力。在實(shí)驗(yàn)中，強(qiáng)激光因具有方向性好、能量密度高等特點(diǎn)，常作為熱源用于材料的局部快速加熱，因而在數(shù)值模擬中研究者們常采用激光作為熱源。

針對具有局部缺陷的薄壁結(jié)構(gòu)的屈曲問題，Raju等人[1]采用有限元方法研究了局部對稱損傷圓板的熱屈曲，用材料力學(xué)性能的局部下降程度描述其影響，得到簡支和固支條件下圓板的屈曲溫度。Jin等人[2]運(yùn)用數(shù)字圖像相關(guān)方法測試了圓形鋁板在熱載荷作用下的屈曲，實(shí)驗(yàn)測得的屈曲溫度與理論結(jié)果接近。Shaterzadeh等人[3]針對含有圓形開孔的對稱與非對稱復(fù)合材料層合板，總結(jié)了開孔尺寸、邊界條件、鋪層方向等對熱屈曲的影響。Yaghoobi等人[4]對功能梯度材料薄板在恒定、線性以及非線性溫度載荷作用下的屈曲進(jìn)行了分析，得到了其在6種不同邊界條件下的屈曲溫度。尹益輝等人[5-6]利用Berger和Galerkin的近似法，導(dǎo)出了固支圓薄板在連續(xù)激光與橫向力聯(lián)合作用下的軸對稱大撓度的表達(dá)式，討論了溫度場與橫向力的聯(lián)合作用對圓板大撓度的作用機(jī)制。鄧可順等人[7-9]分析了結(jié)構(gòu)熱屈曲問題的非線性性質(zhì)，提出了考慮材料性質(zhì)參數(shù)隨溫度變化的熱屈曲試探解法，實(shí)現(xiàn)了在力載荷和溫度載荷的共同作用下結(jié)構(gòu)的熱屈曲有限元分析。Morimoto等人[10]通過解析和數(shù)值模擬方法研究了功能梯度材料板在局部受熱作用下的熱屈曲，得到了材料的不均勻性、長寬比以及熱載區(qū)域大小對臨界屈曲溫度的影響。Hilburger等人[11]采用數(shù)值模擬與實(shí)驗(yàn)測試相結(jié)合的方法，分析了初始缺陷對光滑薄壁圓柱殼屈曲和后屈曲響應(yīng)的影響。Guan等人[12]采用熱彈塑性有限元模型，分析了金屬板在激光輻照和預(yù)加載荷作用下發(fā)生彎曲時(shí)的力學(xué)特性和變形過程，認(rèn)為熱應(yīng)力和預(yù)加應(yīng)力集中是導(dǎo)致金屬板產(chǎn)生變形的主要原因。

局部激光輻照下薄板的穩(wěn)定性受邊界條件、薄板尺寸、激光功率密度、激光輻照時(shí)間等因素的影響。迄今為止，對激光輻照下薄板屈曲問題的研究多集中于理想線彈性屈曲分析，計(jì)算得到的當(dāng)結(jié)構(gòu)發(fā)生屈曲時(shí)的臨界載荷與實(shí)驗(yàn)結(jié)果相差較大，而采用考慮結(jié)構(gòu)的初始缺陷、材料非線性和變形非線性等因素的屈曲分析方法時(shí)，所得結(jié)果更接近實(shí)際情況。

非線性有限元計(jì)算可以采用多種方法，如牛頓-拉普森法、修正牛頓-拉普森法和準(zhǔn)牛頓-拉普森法，但是對于非穩(wěn)定結(jié)構(gòu)，如屈曲問題、接觸問題等，就必須采用弧長法。這是由于在結(jié)構(gòu)非線性有限元分析的求解過程中，當(dāng)出現(xiàn)結(jié)構(gòu)失穩(wěn)或材料軟化時(shí)，傳統(tǒng)的牛頓-拉普森法無法通過極值點(diǎn)，而改進(jìn)的弧長法能夠解決此類問題。本研究采用弧長法對面內(nèi)受單向均勻壓縮的薄板進(jìn)行非線性屈曲分析，綜合考慮結(jié)構(gòu)初始缺陷、幾何非線性、材料非線性和溫度變化對材料參數(shù)的影響，分析激光輻照時(shí)間、激光光斑半徑和薄板厚度對薄板結(jié)構(gòu)屈曲承載力的影響。

2 數(shù)值計(jì)算模型

模型采用尺寸為200 mm×200 mm×2 mm的正方形薄板，激光輻照區(qū)域位于薄板中心，圓形光斑。薄板頂邊和底邊采用固支約束，左、右兩邊不約束任何自由度。薄板材料為鋁合金，材料參數(shù)隨溫度的變化如表1[13]所示，材料的塑性特性如表2所示，其中ρ為密度，μ為泊松比，T為溫度，E為彈性模量，α為熱膨脹系數(shù)，cp為定壓比熱容，λ為熱導(dǎo)率。數(shù)值計(jì)算中，對于表1和表2未列出的其他溫度點(diǎn)下的材料參數(shù)，采用線性插值方法獲取。

表1 材料參數(shù)Table 1 Material parameters

表2 塑性參數(shù)Table 2 Plastic parameters

3 局部激光輻照下薄板的熱屈曲分析

連續(xù)激光輻照金屬表面時(shí)，金屬表層吸收部分激光能量并轉(zhuǎn)化為熱，通過熱傳導(dǎo)在物體內(nèi)擴(kuò)散，使局部溫度上升；局部溫度升高會使材料的性能發(fā)生變化，同時(shí)結(jié)構(gòu)內(nèi)的任意微元因溫度變化而引起的熱膨脹會受到周圍相鄰微元的限制，產(chǎn)生熱應(yīng)力，當(dāng)熱應(yīng)力達(dá)到一定值后，繼續(xù)輻照會導(dǎo)致薄壁結(jié)構(gòu)發(fā)生屈曲。本研究采用溫度-位移耦合求解方法，分析激光功率密度(F)、光斑半徑(R)和薄板厚度(h)對薄板結(jié)構(gòu)熱屈曲時(shí)間的影響，計(jì)算中考慮了結(jié)構(gòu)的幾何非線性和材料非線性。

3.1 激光功率密度對結(jié)構(gòu)熱屈曲的影響

以薄板厚度2 mm、激光功率密度30 W/cm2、光斑半徑20 mm為例，假設(shè)激光被全部吸收，通過有限元分析，得到薄板中心的溫度和應(yīng)力隨時(shí)間變化的規(guī)律，如圖1、圖2所示?？梢?，中心點(diǎn)溫度隨輻照時(shí)間的增加近似呈線性升高；中心點(diǎn)應(yīng)力隨輻照時(shí)間的增加先增大后減小，輻照約3.7 s時(shí)中心點(diǎn)應(yīng)力達(dá)到最大。通過提取底邊的支反力，繪制出支反力-輻照時(shí)間曲線，如圖3所示。從圖3中可以看出：當(dāng)激光功率密度為20 W/cm2時(shí)，輻照6 s左右，支反力開始下降，結(jié)構(gòu)發(fā)生屈曲；當(dāng)激光功率密度為30 W/cm2時(shí)，激光輻照3.7 s左右，結(jié)構(gòu)發(fā)生屈曲；當(dāng)激光功率密度達(dá)到60 W/cm2時(shí)，輻照2 s左右，結(jié)構(gòu)即已發(fā)生屈曲。圖4為發(fā)生熱屈曲的時(shí)間(tB)隨激光功率密度(F)變化的關(guān)系，可見，隨著激光功率密度的增加，薄板發(fā)生熱屈曲的時(shí)間變短。

圖1 中心點(diǎn)溫度變化曲線Fig.1 Temperature-irradiation time curve at the central point

圖2 中心點(diǎn)應(yīng)力變化曲線Fig.2 Stress-irradiation time curve at the central point

圖3 不同功率密度的激光輻照下支反力與輻照時(shí)間的關(guān)系Fig.3 Support reaction force vs.irradiation time for different laser power densities

圖4 發(fā)生熱屈曲時(shí)的輻照時(shí)間與激光功率密度的關(guān)系Fig.4 Laser power density vs.irradiation time of thermal buckling

3.2 激光光斑半徑對結(jié)構(gòu)熱屈曲的影響

保持激光功率密度30 W/cm2不變，改變輻照光斑半徑R，分析2 mm厚薄板的熱屈曲行為。圖5給出了R分別為15、20、25、30和35 mm時(shí)支反力隨輻照時(shí)間變化的曲線，熱屈曲時(shí)間(tB)與輻照光斑半徑(R)的關(guān)系如圖6所示。當(dāng)R=15 mm時(shí)，薄板在8 s左右發(fā)生屈曲；當(dāng)R增大到35 mm時(shí)，輻照2 s時(shí)薄板便發(fā)生熱屈曲；薄板發(fā)生熱屈曲的時(shí)間(tB)隨輻照光斑半徑(R)的增大而縮短。

圖5 不同光斑半徑條件下支反力與輻照時(shí)間的關(guān)系Fig.5 Support reaction force vs.irradiation time for different spot radii

圖6 發(fā)生熱屈曲時(shí)的輻照時(shí)間與光斑半徑的關(guān)系Fig.6 Spot radius vs.irradiation time of thermal buckling

3.3 薄板厚度對結(jié)構(gòu)熱屈曲的影響

將激光功率密度30 W/cm2和光斑半徑20 mm保持不變，改變薄板厚度，分析薄板厚度對結(jié)構(gòu)熱屈曲的影響。對于厚度(h)分別為0.5、1.0、1.5、2.0和2.5 mm的薄板，其支反力隨輻照時(shí)間的變化曲線如圖7所示?？梢钥闯?， 0.5 mm厚的薄板輻照0.3 s左右便發(fā)生熱屈曲，而2.5 mm厚的薄板輻照8 s左右才發(fā)生熱屈曲。發(fā)生熱屈曲時(shí)的輻照時(shí)間(tB)與薄板厚度(h)的關(guān)系如圖8所示。

圖7 不同厚度薄板的支反力與輻照時(shí)間的關(guān)系Fig.7 Support reaction force vs.irradiation time for plates with different thicknesses

圖8 發(fā)生熱屈曲時(shí)的輻照時(shí)間與薄板厚度的關(guān)系Fig.8 Plate thickness vs.irradiation time before thermal buckling

4 局部激光輻照下薄板的屈曲承載力分析

采用第2節(jié)中的計(jì)算模型，將邊界條件和加載情況變?yōu)椋罕“宓走叢捎霉讨Ъs束，頂邊施加垂直于頂邊并指向薄板中心的初始載荷，約束其他5個(gè)方向的自由度。

分析局部激光輻照下薄板的屈曲承載力時(shí)，需要綜合考慮薄板的整體溫度分布、材料性能參數(shù)變化、局部熱應(yīng)力等問題，本研究將此過程分為3步：(1) 采用溫度-位移耦合法分析局部激光輻照下薄板的溫度場和應(yīng)力場；(2) 線性屈曲分析；(3) 將第1步所得的不同輻照時(shí)間下薄板溫度場和應(yīng)力場的分析結(jié)果作為初始狀態(tài)，將第1步中的一階變形作為初始缺陷，采用弧長法分析薄板的屈曲承載力。

4.1 激光輻照時(shí)間對結(jié)構(gòu)屈曲承載力的影響

首先進(jìn)行線性屈曲分析和溫度-位移耦合求解分析。舉例說明：當(dāng)激光功率密度為30 W/cm2、光斑半徑為20 mm、薄板厚度為2 mm時(shí)，線性屈曲分析的一階模態(tài)云圖如圖9所示，激光輻照5 s時(shí)整個(gè)薄板的溫度分布如圖10所示，其中U為相對位移。

圖9 線性屈曲的一階模態(tài)Fig.9 First-order modal of the linear buckling

圖10 激光輻照5 s時(shí)薄板的溫度云圖Fig.10 Temperature nephogram of plate after laser irradiation of 5 s

采用弧長法分析局部激光輻照下薄板的屈曲承載力。提取底邊支反力和頂邊上所有節(jié)點(diǎn)的位移之和繪制曲線，圖11為激光輻照0、1、2、3、4和5 s時(shí)薄板底邊支反力-位移曲線。激光未輻照時(shí)，結(jié)構(gòu)的最大屈曲承載力約為9.30 kN；當(dāng)輻照時(shí)間(t)增加到5 s時(shí)，最大屈曲承載力約為8.70 kN。從圖12所示的屈曲承載力-激光輻照時(shí)間曲線可以看出，結(jié)構(gòu)整體屈曲承載力隨著激光輻照時(shí)間的增加而減小。

圖11 不同輻照時(shí)間下支反力與位移的關(guān)系Fig.11 Support reaction force vs.displacement after different irradiation times

圖12 屈曲承載力與激光輻照時(shí)間的關(guān)系Fig.12 Buckling bearing capacity vs. irradiation time

4.2 激光光斑半徑對結(jié)構(gòu)屈曲承載力的影響

當(dāng)激光功率密度為30 W/cm2，激光輻照光斑半徑(R)分別為15、35、55、75和95 mm，輻照2 s時(shí)，2 mm厚薄板的結(jié)構(gòu)屈曲承載力如圖13所示，屈曲承載力與輻照光斑半徑的關(guān)系如圖14所示。當(dāng)輻照光斑半徑為15 mm時(shí)，屈曲承載力約為9.25 kN；當(dāng)光斑半徑增大到95 mm時(shí)，屈曲承載力下降至7.90 kN左右。從圖14中可以看出，結(jié)構(gòu)整體屈曲承載力隨輻照光斑半徑的增大而降低。

4.3 薄板厚度對結(jié)構(gòu)屈曲承載力的影響

當(dāng)激光輻照功率為30 W/cm2、光斑半徑為20 mm、輻照時(shí)間為2 s時(shí)，厚度(h)分別為0.5、1.0、1.5、2.0和2.5 mm的薄板的支反力-位移曲線如圖15所示。0.5 mm厚薄板的屈曲承載力在0.50 kN左右，而2.5 mm厚薄板的屈曲承載力則達(dá)到18.00 kN，結(jié)構(gòu)整體屈曲承載力隨薄板厚度的增加而增大，如圖16所示。

圖13 不同光斑半徑的激光輻照2 s后支反力與位移的關(guān)系Fig.13 Support reaction force vs.displacement after laser irradiation of 2 s with different spot radii

圖14 激光輻照2 s后屈曲承載力與光斑半徑的關(guān)系Fig.14 Buckling bearing capacity vs.spot radius after laser irradiation of 2 s

圖15 激光輻照2 s后不同厚度薄板的支反力與位移的關(guān)系Fig.15 Support reaction force vs.displacement for plates with different thicknesses after laser irradiation of 2 s

圖16 激光輻照2 s后屈曲承載力與薄板厚度的關(guān)系Fig.16 Buckling bearing capacity vs.plate thickness after laser irradiation of 2 s

5 結(jié) 論

(1) 隨著激光輻照時(shí)間的增加，薄板中心點(diǎn)的溫度逐漸升高，中心點(diǎn)應(yīng)力先升高后下降。

(2) 隨著激光功率密度的增加，薄板發(fā)生熱屈曲的時(shí)間呈線性下降；當(dāng)激光輻照光斑半徑增加時(shí)，薄板發(fā)生熱屈曲的時(shí)間呈二次函數(shù)關(guān)系下降；增加薄板厚度可以大幅度延長薄板發(fā)生熱屈曲的時(shí)間。

(3) 增加激光輻照時(shí)間和輻照光斑半徑都會降低整個(gè)薄板的屈曲承載力，且均近似呈線性關(guān)系；薄板厚度對結(jié)構(gòu)整體屈曲承載力的影響更為顯著，隨著薄板厚度的增加，屈曲承載力近似呈指數(shù)增大。因此，增加薄板厚度或采用加筋方法對提高薄壁結(jié)構(gòu)抗熱屈曲能力具有顯著的效果。

[1] RAJU K K,NAIDU N R,RAO G V.Thermal buckling of circular plates with localized axisymmetric damages [J].Comput Struct,1996,60(6):1105-l109.

[2] JIN T L,HA N S,GOO N S.A study of the thermal buckling behavior of a circular aluminum plate using the digital image correlation technique and finite element analysis [J].Thin-Walled Struct,2014,77:187-197.

[3] SHATERZADEH A R,ABOLGHASEMI S,REZAEI R.Finite element analysis of thermal buckling of rectangular laminated composite plates with circular cut-out [J].J Therm Stresses,2014,37(5):604-623.

[4] YAGHOOBI H,TORABI M.Exact solution for thermal buckling of functionally graded plates resting on elastic foundations with various boundary conditions [J].J Therm Stresses,2013,36(9):869-894.

[5] 尹益輝,陳裕澤,鄧禮輝.連續(xù)激光作用下圓板大撓度的理論分析 [J].強(qiáng)激光與粒子束,1999,11(6):682-685.

YIN Y H,CHEN Y Z,DENG L H.Large deflections of a circular plate irradiated by CW laser [J].High Power Laser and Particle Beams,1999,11(6):682-685.

[6] 尹益輝,陳 剛,陳裕澤,等.熱力聯(lián)合作用下圓板屈曲臨界狀態(tài)分析 [J].強(qiáng)度與環(huán)境,2000(2):45-50.

YIN Y H,CHEN G,CHEN Y Z,et al.Analysis of buckling critical loading of circular plate loaded by combined heat and lateral force [J].Structure & Environment Engineering,2000(2):45-50.

[7] 鄧可順,張亞輝.考慮材料性質(zhì)參數(shù)隨溫度變化的熱屈曲試探解法 [J].大連理工大學(xué)學(xué)報(bào),1999,39(3):358-362.

DENG K S,ZHANG Y H.Trial and error method of thermal buckling for complex structures [J].Journal of Dalian University of Technology,1999,39(3):358-362.

[8] 鄧可順.復(fù)雜結(jié)構(gòu)熱屈曲有限元分析 [J].計(jì)算結(jié)構(gòu)力學(xué)及其應(yīng)用,1994,11(2):130-139.

DENG K S.Finite element analysis of thermal buckling for complex structures [J].Computational Structural Mechanics and Applications,1994,11(2):130-139.

[9] DENG K S,JI Z,DAVIES A W,et al.Thermal buckling of axially precompressed cylindrical shells irradiated by laser beam [J].AIAA J,2000,38(10):1789-1794.

[10] MORIMOTO T,TANIGAWA Y,KAWAMURA R.Thermal buckling of functionally graded rectangular plates subjected to partial heating [J].Int J Mech Sci,2006,48:926-937.

[11] HILBURGER M W,STARNES J H Jr.Effects of imperfections on the buckling response of compression loaded composite shells [J].Int J Non-Linear Mech,2002,37(4/5):623-643.

[12] GUAN Y J,SUN S,ZHAO G Q.Finite element modeling of laser bending of pre-loaded sheet metals [J].J Mater Process Technol,2003,142(2):400-407.

[13] 孫承緯.激光輻照效應(yīng) [M].北京:國防工業(yè)出版社,2002:28-32.

Sun C W.Effect of laser irradiation [M].Beijing:National Defense Industry Press,2002:28-32.