鄭玉卿,朱西產(chǎn)
棱邊強化薄壁方管軸向壓潰吸能特性?
鄭玉卿1,2,朱西產(chǎn)1
(1.同濟大學汽車學院,上海 201804; 2.湖州師范學院工學院,湖州 313000)
采用Hypermesh 9.0和Lsdyna 971數(shù)值仿真分析了棱邊強化對薄壁方管軸向壓潰過程能量吸收機理的影響?;趩蝹€方管折疊單元能量平衡方程,采用分項能量修正法和影響因子歸一化建立了棱邊強化薄壁方管軸向壓潰能量平衡方程。理論推導出棱邊強化方管平均壓潰力預測公式和吸能增幅表達式,其中,平均壓潰力的仿真結果與預測公式計算結果吻合較好。Q235方管壓潰試驗結果表明,棱邊強化方管平均壓潰力公式也可用于預測原始方管,且效果更好。吸能對比表明:棱邊強化后,棱邊與平板的塑性變形在軸向壓潰過程中存在強耦合關系,兩部分的吸能均有提升。對于截面為56.36mm×56.36mm、厚度為1.0mm的薄壁方管,棱邊約強化4倍,4條棱邊僅占薄壁方管截面9.09%,理論上可使方管總體吸能量提升近30%,可見棱邊強化能顯著提升薄壁方管軸向壓潰吸能水平。
薄壁方管;棱邊強化;能量吸收;平均壓潰力;數(shù)值仿真;壓潰試驗
為減少交通事故乘員傷害數(shù)和碳排放,要求汽車車身結構做得更加耐撞和輕量[1]。當前增強結構截面抗力的方法通常有兩種:其一是增加沖壓件的厚度,即增加局部結構質量,這會導致油耗增加;其二采用高強度或超高強度鋼板,價格昂貴,但優(yōu)勢是利于實現(xiàn)車身輕量化[2],故被廣泛應用。傳統(tǒng)的高強度鋼主要由鐵素體基體組成,先進高強度鋼則由鐵素體和其他相的多相基體組成,多相基體對材料的力學性能、可成形性能甚至吸能特性都有明顯影響[3]。而對于980MPa級以上的超高強度鋼板,由于冷沖壓的成形特性和回彈的制約,需加熱至其伸展性較好的高溫狀態(tài),壓力成形后馬上進行冷卻硬化處理,這種熱沖壓成形技術已被成功開發(fā)并普遍使用,但成本較高[4]。超高強度鋼材料具有超高強度特性優(yōu)勢,也存在塑性延伸性小、易斷裂等缺點[5]。近年來,有研究報告指出高強度鋼板的點焊和激光焊接的熔解熱影響區(qū)會導致母材強度下降的現(xiàn)象[6];同時還有根據(jù)零部件在車身結構中承擔載荷和變形分布的必要性,在同一個零部件中分開制造高強度區(qū)和高延伸區(qū),其工藝成本就更加高昂[7]。
為克服傳統(tǒng)截面抗力增強方法和車身連接工藝的缺點,文獻[8]中論述了帽型結構棱邊強化后在壓彎作用下壓潰側平板的彈塑性屈曲失穩(wěn)過程和全截面抵抗力增強的影響因子,試驗和仿真均表明其截面抗壓和抗彎性能有顯著提高。薄壁方管結構塑性變形特性在汽車碰撞能量吸收系統(tǒng)中起重要作用。漸進式單元折疊是一種有效的吸能崩塌機制,低碳鋼薄壁結構靜動態(tài)軸向壓潰和能量吸收方面的試驗和理論研究已有大量報道,如圓管和方管[9-10]、單帽和雙帽型管[11-12],而高強度鋼對于汽車碰撞安全應用仍是一個挑戰(zhàn),其幾何形狀和材料屬性均會影響車身結構的失效變形模式[13]。本文中選取棱邊強化薄壁方管為研究對象,假設僅僅只將其受焊接和機械加工影響較小的棱邊進行超高應力強化熱處理,探討棱邊應力強化對于在準靜態(tài)軸向壓潰載荷作用下方管分部能量吸收機制的影響,量化分析其軸向壓潰吸能增強效果,為后續(xù)車身核心承載件的輕量化和低成本設計及應用提供技術支撐。
1.1 基本模型
圖1(a)為1/4方管的初始幾何形狀,a為方管的寬度或等效寬度,H為理想化方管折疊單元的半高度。圖1(b)為1/4方管受準靜態(tài)軸向壓潰載荷作用時的典型吸能分析模型,該理想化折疊模型是由Wierzbicki和 Abramowicz[14]給出的,2ψ=π/2是相鄰平面UBC和平面UBA之間的夾角,假設方管變形過程中2ψ和邊長都是不變的。因穩(wěn)定屈曲是期望的屈曲變形吸能模式,根據(jù)文獻[15]中的試驗結果,選取適當?shù)膶捄癖萢/h,可避免發(fā)生Euler失穩(wěn)屈曲模式。文獻[16]中利用系列PVC管,獲取了方管連續(xù)變形的過程,并發(fā)現(xiàn)變形過程中有瞬間為不動的塑性鉸,即水平固定鉸(如AB和BC),但整體上看,整個變形過程中4個棱角處初始為垂直的塑性鉸,又稱傾斜移行鉸(UB和BL)逐漸變成傾斜,水平固定鉸(AB和BC)發(fā)生彎曲變形,其中AB鉸外凸變形,BC鉸內(nèi)凹變形。理想化折疊模型的塑性變形都發(fā)生在塑性鉸上,如果塑性變形沒有擴展,它是可以接受的,但通過試驗和仿真過程觀察,注意到連續(xù)變形過程中傾斜移行鉸UB和BL均在移動,要滿足動力學變形許可,必然存在連續(xù)的面內(nèi)塑性擴展過渡區(qū)[17],也即,水平固定鉸的AB與BC和傾斜移行鉸的UB和BL在B點皆出現(xiàn)小幅彎曲塑性擴展,最后兩種塑性鉸在B點交匯處形成一個環(huán)形殼曲面。
圖1 方管軸向壓潰吸能分析模型
1.2 原始方管壓潰吸能預測
沿用Wierzbicki和Abramowicz的原始方管折疊單元壓潰過程的能量耗散分解法,假設每一個折疊的壓潰距離為2H,將上述1/4方形截面塑性變形區(qū)的能量耗散分為3個部分:環(huán)形殼曲面拉伸的能量耗散W1、水平固定塑性鉸能量耗散W2和傾斜移行塑性鉸耗散W3,根據(jù)能量守恒原則可得
式中Wert為折疊單元在軸向外力擠壓作用下吸收的總能量。
考慮到一般吸能特性模型頂部和底部均有夾持狀態(tài)下的水平固定塑性鉸線,W2能量應該加倍。整個折疊單元的外力功等于1/4折疊單元壓潰過程中內(nèi)部能量耗散總和的4倍,同時將文獻[18]中的公式代入式(1),則方管壓潰過程中一個完整折疊單元能量平衡方程式為
式中:Pmo為原始方管平均壓潰力;Mo為單位長度的塑性極限彎矩;b為圓環(huán)子午線方向半徑;h為方管壁厚;I1,I3均為含ψ的函數(shù)表達式,當ψ=π/4,可求得 I1(π/4)=0.58,I3(π/4)=1.11。
通過一系列數(shù)學推算[19]可得
實際上有效壓潰距離小于2H,Abramowicz和Jones基于a/h=35~50系列軟鋼方管靜態(tài)壓潰試驗得到了每個折疊單元實際有效壓潰距離為2H的0.77~0.78倍[9],因轎車車身薄壁箱型結構的寬厚比與試驗方管接近,結合本文中研究方管的尺寸,有效壓潰距離也取為2H的0.78倍,可得原始方管平均壓潰力修正式為
假設方管軸向壓潰距離為δ,因Mo=Yh2/4,代入式(4),可得整個準靜態(tài)壓潰過程原始方管總體吸能 Wo為
式中Y為材料屈服強度。
值得注意的是,將式(4)預測值與文獻中不同寬厚比的低碳鋼方管靜態(tài)壓潰試驗值對比可知,低碳鋼屈服強度和抗拉強度分別為265和325MPa,預測值比試驗值偏低約17%左右[9],而采用抗拉極限應力或流動應力替代公式中的屈服強度Y求得的平均壓潰力與試驗值更為接近。式(4)和式(5)分別為原始方管軸向壓潰過程平均壓潰力預測公式和吸能預測公式,故本文預測值與試驗值相比可能偏低。
1.3 棱邊強化方管基本特征
棱邊強化方管選用母材為SUS440A鋼板,其化學成分見表1。440A不銹鋼淬火硬化性能優(yōu)良、硬度高,較440B鋼和440C鋼有高的韌性,可用作刃具、量具和軸承。棱邊強化方管截面形狀尺寸參數(shù)見圖2,強化棱邊部用粗實線標出,平板部用細實線標出,長度設為l。
表1 SUS440A鋼化學成分%
圖2 方管橫截面尺寸參數(shù)
豐田汽車公司熱處理試驗表明,該鋼板可通過激光淬火和合理的冷卻控制,獲得超高強度特性。熱處理前后的材料單向拉伸得到的真實應力 真實應變數(shù)據(jù)對比如圖3所示[8],可知平板部母材的屈服強度Yp約為328MPa,棱邊部經(jīng)局部超高強度化調質處理,其屈服強度Yr可達1 300MPa左右,在此引入兩者屈服強度比λ:
圖3 真實應力 真實應變曲線
2.1 數(shù)值模型
采用Hypermesh 9.0和Lsdyna 971軟件對上述原始方管和棱邊強化方管準靜態(tài)軸向壓潰過程進行數(shù)值模擬,可知 λ≈4,設 c=50mm,s=5mm,r=3.18mm,h=1.0mm,長度為 l=400mm。 仿真計算前,先做了一個仿真工況進行仿真精度驗證:選取邊長為1m,厚度為10mm的薄板,對邊簡支、對邊固定且面載荷為1MPa,線彈性CAE分析結果表明,選用16號全積分殼體單元,單元邊長與厚度比為2~12.5,厚度方向積分點設為5個,薄板中部最大撓度仿真值與精確解析解的偏差不超過1%。于是方管劃分為4節(jié)點16號殼體單元,為捕獲軸向壓潰變形過程更多的變形行為,單元平均邊長為2mm,厚度方向設置5個積分點,方管平板部和棱邊部材料模型均選用Mat24多線性彈塑性材料模型,材料特性曲線由圖3中的拉伸試驗數(shù)據(jù)轉換獲得,壓板也劃分16號薄殼單元,材料模式選用Mat20,兩者密度、彈性模量和泊松比均設為7.85×10-9t/mm3,205 800MPa和0.3。方管一端固定,另一端受剛性壓板軸向擠壓,對剛性壓板施加位移載荷,壓板運動方向僅限軸向,壓潰速度為1.0m/s,總壓潰距離設為320mm,計算時間為0.32s,剛性壓板與方管之間采用?Contact_Automatic_Surface_To_Surface 接觸算法,擠壓力則通過設置接觸控制卡中的SPR和MPR設置為1來提取,方管本身還采用了?Contact_Automatic_Single_Surface接觸算法,摩擦因數(shù)均設為0.2,模型總單元數(shù)為22 625,原始方管和棱邊強化方管準靜態(tài)壓潰模型分別見圖4和圖5。
2.2 準靜態(tài)模擬方法
圖4 原始方管準靜態(tài)壓潰模型
圖5 棱邊強化方管準靜態(tài)壓潰模型
Lsdyna顯式積分的最小時間步長,主要由最小單元特征長度和材料中應力波的傳播速度確定。通常顯式積分算法因時間步長較小,只適合瞬態(tài)分析。而準靜態(tài)壓潰時間較長,通常可通過提高加載速率的方法來縮短計算時間,也可采用Lsdyna的質量縮放技術來加快計算過程。與質量縮放技術相比,提高加載速率可更快得到比較準確的準靜態(tài)模擬結果,此時只需消除材料模型中動態(tài)應變率效應的影響[20];當然為了保證模擬計算的準確性,文獻[21]中研究發(fā)現(xiàn),提高加載速率后的動能與內(nèi)能之比值應盡可能小,建議控制在2%以內(nèi),上述方管仿真模型滿足這一要求。
2.3 塑性壓潰過程對比
圖6為原始方管和棱邊強化方管準靜態(tài)軸向壓潰過程中壓潰距離分別為 20,100,200,300和320mm時塑性變形的對比。首先,從整體塑性變形情況來看:壓潰距離為20和100mm時,屬于兩種方管準靜態(tài)軸向壓潰初始階段,均有不同程度的不穩(wěn)定塑性變形發(fā)生,之后均產(chǎn)生了穩(wěn)定的塑性變形,但發(fā)生時間和周期性折疊單元數(shù)不同,其中原始方管較早進入穩(wěn)定壓潰階段,且第2~6個單元形成了周期性折疊單元,而棱邊強化方管進入較晚,第3~5個單元形成了周期性折疊單元;被擠壓的周期性折疊單元高度和寬度明顯不同,原始方管折疊單元外凸部分的高度和寬度均大于棱邊強化方管的,這說明棱邊強化方管的棱邊在彎折過程中,彎折距離更長,將應力強化棱邊周邊的平板卷入更多,導致平板部母材的塑性拉伸和彎曲變形量更大,會加劇平板部能量耗散速率;壓潰距離為320mm時,比較兩種方管第6個折疊單元可發(fā)現(xiàn)原始方管的第6個折疊單元仍是完整和周期性的,而棱邊強化方管的第6個折疊單元是不完整且非周期性的,這是因為棱邊強化方管端部固定約束導致強行壓潰所得的,該現(xiàn)象表明應力強化棱邊的支撐效應使方管折疊單元長度變長,從而無法生成與原始方管相同數(shù)量的周期性折疊單元。上述塑性變形對比分析表明,應力強化棱邊的存在對方管的塑性變形機理確實造成了不小的影響,同時也驗證了本文中選取的合理寬厚比能夠保證方管產(chǎn)生穩(wěn)定的軸向壓潰塑性變形。
圖6 不同壓潰距離時原始方管和棱邊強化方管的塑性變形對比
2.4 棱邊強化對分項吸能的影響
通過仿真動畫演示和估算,觀察棱邊強化方管與原始方管之間變形模式不同點,可歸納為3個方面:(1)平板區(qū)的環(huán)形曲面能量耗散大部分仍是通過環(huán)形面內(nèi)塑性拉伸來實現(xiàn)運動許可的塑性變形模式,當?shù)?個折疊形成后會受到第2個折疊棱邊擠壓的影響,這是因為棱邊應力強化后導致棱邊本身塑性變形困難同時平板易發(fā)生屈曲,就出現(xiàn)了同一棱邊由對折變成Z型錯位拉伸卷曲平板的情形,經(jīng)測定單個折疊單元受棱邊強化影響的平板區(qū)約占整個折疊單元的1/3,故單個折疊單元的平板部耗散能量W1約增大1/3;(2)水平固定塑性鉸(AB和BC),通過平板部和棱邊部的彎曲塑性變形來實現(xiàn)能量耗散,彎曲變形過程中強化棱邊直接參與,故W2會有所上升;(3)對于傾斜移行塑性鉸(UB和BL),主要由棱邊部和平板部的塑性移動彎曲來實現(xiàn)能量耗散,棱邊強化導致棱邊部卷曲困難,只有約1/3棱邊部參與傾斜移行彎曲,因屈服強度比為λ,故W3提高幅度僅約為(λ/3-1)。當然,隨著λ的降低,壓潰過程參與傾斜移行彎曲的棱邊比例會提高,且受影響的平板區(qū)也會隨之減少,故W1,W2和W3的能量增幅均會隨λ的改變而改變。
由2.4節(jié)分析可知,正是棱邊應力強化λ倍導致方管軸向壓潰過程分項吸能量發(fā)生變化,故為數(shù)學表達方便,將W1,W2,W3能量耗散增幅進行影響因子歸一化處理,即統(tǒng)一采用棱邊與平板之間屈服強度比λ來近似表示,增幅修正式分別表示為(λ-1)/9,λ.s,(λ-1)/9。
于是分別將3個修正式代入式(2),可得棱邊強化方管整個折疊單元吸能模型的能量平衡方程為
式中Pmr為棱邊強化方管平均壓潰力。將平均載荷取為
將I1,I3和 λ代入式(9)~式(11),計算可得A1=24.747,A2= 12.56,A3= 11.84,式(8)僅含有 2個未知數(shù)b和H,可分別令1階偏導等于0,即
求解式(12)可得
將式(13)代入式(8),有
從而,得到
和
因方管寬度與厚度比為56.36,近似取有效壓潰距離是2H的0.78倍,可得
同理,假設棱邊強化方管軸向壓潰距離也為δ,可得靜態(tài)軸向壓潰過程棱邊強化方管總體吸能Wr為
式(17)和式(18)分別為棱邊強化方管靜態(tài)軸向平均壓潰力預測公式和吸能公式。將 Y=328MPa,λ≈4,c=50mm,s=5mm 和 h=1.0mm 分別代入式(4)和式(17),可得原始方管和棱邊強化方管的平均壓潰力分別為15.424和20.032kN。
因 a=c+2r=c+4s/π,將其代入式(5)可得
將式(18)減式(19)后再除以式(19),再令棱邊強化區(qū)域占比為ω=s/c,可得棱邊強化方管與原始方管之間的吸能增幅ζ為
由式(20)可知棱邊強化方管吸能增幅ζ與厚度h無關,僅與應力強化棱邊所占截面的區(qū)域占比ω和屈服強度比λ有關。代入棱邊強化區(qū)域占比ω=9.09%和屈服強度比 λ≈4,可得 ζ=0.294,這說明對于截面為56.36mm×56.36mm、厚度為1.0mm和圓角半徑為3.18mm的方管,棱邊強化約4倍,強化區(qū)域僅占截面周長的9.09%,方管吸能水平理論上可提高29.4%。
4.1 平均壓潰力
平均壓潰力是表征薄壁管靜態(tài)壓潰過程中穩(wěn)定壓潰階段的平均抗力水平。圖7描繪了h=1.0mm原始方管和棱邊強化方管在準靜態(tài)軸向壓潰過程中不同壓潰距離對應的壓潰力變化曲線,根據(jù)壓潰力曲線波動情況可判斷原始方管在壓潰距離40mm左右進入穩(wěn)定壓潰階段,之后形成了5個連續(xù)的W型壓潰力波動周期,從300mm開始出現(xiàn)壓潰力陡增現(xiàn)象,說明折疊單元組進入了壓實階段;而棱邊強化方管在初期階段,波動較大,直到壓潰距離120mm左右,才開始進入穩(wěn)定壓潰階段,緊接著出現(xiàn)了3個連續(xù)的W型壓潰力波動周期,最后一個波動周期是不完整的。
圖7 h=1.0mm原始和棱邊強化方管壓潰力曲線對比
上述軸向壓潰過程描述恰好補充解釋了圖6中塑性變形過程和折疊單元形成過程。對兩種方管穩(wěn)定壓潰階段壓潰力取平均值,獲得平均壓潰力曲線,其值分別為13.669和20.717kN,可見通過棱邊應力強化處理,棱邊強化方管仿真平均壓潰力有較大幅度的提升,增幅達到了51.5%。
4.2 能量吸收
因原始方管被壓到300mm時開始出現(xiàn)壓潰力曲線陡增現(xiàn)象,即折疊單元組進入了壓實階段,但是棱邊強化方管并未出現(xiàn),因此需排除這一干擾因素進行能量吸收對比,只處理300mm之前的吸能仿真結果。圖8~圖10分別為原始方管和棱邊強化方管總體、棱邊部和平板部能量吸收曲線對比??梢钥闯?,棱邊強化方管總體、棱邊部和平板部吸能水平都有較大幅度提高,棱邊部吸能提升幅度雖高于平板部,但是兩者吸能值凈增量較接近。
圖8 總體吸能對比
圖9 棱邊部吸能對比
圖10 平板部吸能對比
4.3 仿真與理論結果對比
表2列出了原始方管和棱邊強化方管準靜態(tài)壓潰過程中平均壓潰力和各分部能量吸收對比。由表可知:原始方管理論值大于仿真值,兩者偏差稍大,為11.4%;棱邊強化方管理論值稍低于仿真值,偏差很小,僅為3.4%。原始方管偏差較大的原因是仿真采用了倒圓角方管,而理論公式推導并未考慮倒圓角的影響,通過對直角方管進行CAE驗證分析得到平均壓潰力為 14.606kN,與理論值偏差縮小至5.3%,說明方管棱邊處的倒角形式對平均壓潰力存在一定程度的影響,而棱邊強化方管偏差小是因為預測公式的推導過程重點考慮了應力強化棱邊對分部吸能機理的影響。棱邊強化方管與原始方管對比,理論壓潰力增幅為29.9%,而仿真壓潰力增幅達到了51.6%,造成如此大差異的原因主要是原始方管仿真平均壓潰力偏小所致。吸能對比可知3點:(1)原始方管棱邊部吸能均約占總吸能的1/4,而平板部約占3/4,棱邊強化后,棱邊部吸能占比提高到1/3左右,而平板部約為2/3,說明應力強化后棱邊部吸能水平提高;(2)棱邊應力強化了近4倍,但棱邊部吸能只增加了90.9%,并非預期的4倍關系,同時平板部吸能也增加了36.1%,棱邊部增幅雖遠大于平板部,但吸能增量上仍是平板部稍多一些,兩者凈增量在數(shù)值上較為接近,說明棱邊強化方管在軸向壓潰過程中棱邊部和平板部的吸能方面存在強耦合關系,驗證了平板部在應力強化棱邊影響下會加劇塑性拉伸擴展,從而增加塑性變形吸能效果,同時應力強化棱邊的塑性折彎和移行范圍會隨應力增強幅度提升而下降,而平板部的塑性折彎和拉伸則會加??;(3)棱邊強化方管的棱邊強化區(qū)域僅占整個截面的9.09%,理論吸能水平可提高29.4%,仿真吸能水平可提升50.1%,兩者雖然由于棱邊處的倒角影響存在數(shù)值上偏差,但棱邊強化工藝對于薄壁方管的能量吸收增強效果是顯著的。
表2 原始方管和棱邊強化方管平均壓潰力和吸能對比
4.4 壓潰試驗驗證
原始方管理論預測值與文獻中大量的低碳鋼方管試驗值對比,發(fā)現(xiàn)理論值均偏小,這是因為理論公式是基于方管所有部位材質屬性相同的前提下而導出的,實際上原始方管在加工成形過程中,棱邊部彎折成直角過程存在一定程度的塑性強化現(xiàn)象。為了驗證這一塑性強化現(xiàn)象,裁取了若干塊長度為200mm、寬度為 20或 10mm、厚度均為 1.0mm的Q235平板,經(jīng)冷彎加工成直角的棱邊樣件,試驗樣件如圖11所示,然后對平板和棱邊樣件進行重復的單向拉伸試驗。拉伸速度為5mm/min,拉伸應力應變曲線如圖12所示,由圖可知,折彎直角工藝導致棱邊部屈服強度由276MPa提升至443MPa左右且延伸率有所降低。繼而加工了長、寬、厚分別為60.0,60.3和1.0mm的Q235方管,一側平板部采用焊接封閉,長度均為300mm,進行3次靜態(tài)重復壓縮試驗,如圖13所示,試驗速度為5mm/min,獲得的平均壓潰力試驗值約為14.256kN。
圖11 平板和棱邊拉伸試件
圖12 工程應力應變曲線
圖13 Q235方管軸向壓潰試驗
此時令 s=10mm,可求得 r=6.37mm,c=47.56mm,λ≈1.605代入式(17)求得理論預測值為14.099kN??梢姡紤]棱邊強化的理論預測值與試驗值之間偏差很小,僅為1.1%,而未考慮棱邊強化效果的原始方管理論預測值為13.274kN,與試驗值之間偏差較大,達到了7.4%。此外,未考慮棱邊強化的仿真結果為11.301kN,更是不建議采用,而考慮棱邊強化的仿真結果是13.731kN,與試驗值較接近,因此建議采用。
綜上所述,采用考慮棱邊強化效應的平均壓潰力理論預測公式(式(17))來預測原始單一材質方管的平均壓潰力效果也非常好,因為原始方管在成形過程中彎折棱邊確實存有塑性應變強化現(xiàn)象,這一點也得到了考慮棱邊強化的仿真結果驗證,同時也從側面解釋了之前仿真結果一直偏低的原因。而國際上對于低碳鋼薄壁件平均壓潰力預測,仍一般沿用極限抗拉應力或極限抗拉應力的倍數(shù)來替代屈服應力的做法,但對于中高碳鋼薄壁件,效果似乎并不是很好。
(1)分析了應力強化棱邊對薄壁方管軸向壓潰過程中能量耗散機制的影響,通過修正原始方管單個折疊單元的能量平衡方程,理論導出了棱邊強化方管平均壓潰力預測公式和吸能增幅表達式,量化分析了棱邊強化效應對方管平均壓潰力和吸能水平的影響。
(2)數(shù)值仿真再現(xiàn)了原始方管和棱邊強化薄壁方管軸向壓潰過程的塑性變形特征,并對兩種方管平均壓潰力和能量吸收進行了詳細的對比和討論。理論和仿真壓潰力對比可知兩者吻合較好,原始方管偏差較大為11.4%,而棱邊強化方管偏差相對較小,為3.4%。由原始方管平均壓潰力理論值與文獻中低碳鋼方管的試驗值對比可知,理論預測值偏低17%左右,國際上一般采用抗拉強度來替代屈服應力求解預測值。而Q235方管軸向壓潰試驗結果分析表明,本文導出的棱邊強化方管平均壓潰力公式也可用于預測考慮棱邊塑性強化效應的原始方管平均壓潰力,且預測效果更好,同時也解釋了困擾已久的原始方管理論和仿真結果偏小的原因在于,未考慮原始方管加工過程棱邊折彎塑性強化現(xiàn)象。
(3)吸能對比表明,棱邊強化方管的棱邊部與平板部之間存在強耦合關系,應力強化棱邊部通過增強自身的屈服應力水平來提高平板部的塑性變形吸能能量,從而全面提升棱邊強化方管的總體吸能水平。
(4)對于棱邊強化區(qū)域僅占薄壁方管截面9.09%,棱邊應力強化4倍的薄壁方管,其吸能水平理論上可提升近30%,說明棱邊應力強化對于薄壁方管準靜態(tài)軸向壓潰吸能水平的提升效果是顯著的。制備棱邊超高應力強化方管的關鍵工藝在于棱邊局部熱處理和材質選擇,靜態(tài)軸向壓潰驗證試驗尚待今后繼續(xù)進行。
[1] 崔崇楨,曹立波,顏凌波,等.車與車正面碰撞安全性仿真研究[J].汽車工程,2015,37(6):661-666.
[2] 范子杰,桂良進,蘇瑞意.汽車輕量化技術的研究和進展[J].汽車安全與節(jié)能學報,2014,5(1):1-16.
[3] MALLICK P K.汽車輕量化:材料、設計與制造[M].于京諾,譯.北京:機械工業(yè)出版社,2012.
[4] ALTAN T.Hot-stamping boron alloyed steels for automotive parts[J].Part I: Process Methods and Uses.Stamping Journal,2006:40-41.
[5] DYKEMAN J, MALCOLM S, HUANG G, et al.Characterization of edge fracture in various types of advanced high strength steel[C].SAE Paper 2011-01-1058.
[6] KIM C H,CHOI JK,KANGM J,et al.A study on the CO2laser welding characteristics of high strength steel up to 1500 MPa for automotive application[J].Journal of Achievements in Materials and Manufacturing Engineering,2010,39(1):79-86.
[7] PEREZ-SANTIAGO R,BILLUR E,ADEMAJA,et al.Hot stamping a B-pillar with tailored properties:experiments and preliminary simulation results[C].4thInternational Conference,Hot Sheet Metal Forming of High-performance Steel CHS2, Lulea, Sweden,2013:83-90.
[8] 嵇明,安孫子拓也,岡村俊明.異なる材料強度の稜線を有する箱形斷面の圧縮曲げ強度[C].自動車技術會論文集.とうきょう:日本自動車エンジニア學會,2015,46(3):645-652.
[9] ABRAMOWICZ W,JONES N.Dynamic progressive buckling of circular and square tubes[J].International Journal of Impact Engineering,1986,4(4):243-270.
[10] 鄭玉卿,朱西產(chǎn),董學勤,等.Q235B直縫焊管軸向沖擊性能的理論和試驗研究[J].振動與沖擊,2016,35(20):98-103.
[11] WHITE M D,JONESN.Experimental quasi-static axial crushing of top-hat and double-hat thin-walled sections[J].Int.J.Mech.Sci.,1999,41:179-208.
[12] WHITE M D,JONES N.Experimental study into the energy absorbing characteristics of top-hat and double-hat sections subjected to dynamic axial crushing[J].Proc.Instn Mech.Engrs, Part D:J.Automobile Engineering,1999,213(D3):259-278.
[13] SCHNEIDER F,JONES N.Impact of thin-walled high-strength steel structural sections[J].Proc.Instn Mech.Engrs, Part D: J.Automobile Engineering,2004,218(D3):131-158.
[14] WIERZBICKI T,ABRAMOWICZ W.On the crushing mechanics of thin-walled structures[J].Journal of Applied Mechanics,1983,50:727-734.
[15] REID SR,REDDY T Y,GRAY M D.Static and dynamic axial crushing of foam-filled sheet metal tubes[J].International Journal of Mechanical Sciences,1986,28:295-322.
[16] MENG Q,AL-HASSANI ST S,SODEN P D.Axial crushing of square tubes[J].International Journal of Mechanical Sciences,1983,25(9/10):747-773.
[17] CALLADINE CR.Theory of shell structure[M].Cambridge University Press,1983.
[18] 鄭玉卿,朱西產(chǎn),胡強,等.棱線強化薄壁方管軸向壓潰力學特性[J].同濟大學學報(自然科學版),2016,44(6):949-954.
[19] 余同希,盧國興.材料與結構的能量吸收[M].華云龍,譯.北京:化學工業(yè)出版社,2006.
[20] 桂良進,范子杰,王青春.泡沫填充圓管的動態(tài)軸向壓潰吸能特性[J].清華大學學報(自然科學版),2004,43(5):709-712.
[21] 王青春,范子杰.利用LS-DYNA計算結構準靜態(tài)壓潰的改進方法[J].力學與實踐,2003,25(3):20-23.
Energy Absorption Characteristics of Ridgeline Strengthened Thin-walled Square Tube in Axial Crushing
Zheng Yuqing1,2& Zhu Xichan1
1.School of Automotive Studies, Tongji University, Shanghai 201804; 2.School of Engineering, Huzhou University, Huzhou 313000
A numerical simulation is conducted with Hypermesh 9.0 and Lsdyna 971 to analyze the effects of ridgeline strengthening(RS)on the energy absorption mechanism of thin-walled square tube subjected to axial crushing.Based on the energy balance equation of a folding element in common square tube,the energy balance equations of RSthin-walled square tube in axial crushing are set up by using partial energy correction and influencing factor normalization.The predicting formulas for mean crushing force and the expression of energy absorption increasing margin of RSthin-walled square tube are derived theoretically,in which the results of simulation agree well with the results calculated by predicting formulas.Crushing test results of Q235 square tube show that the predicting formula of mean crushing force for RS thin-walled square tube is also suitable for original square tube with better prediction result.The comparison of energy absorbed indicates that there is a strong coupling relationship of plastic deformation between stress-strengthened ridgeline and plate in axial crushing,so the energy absorbed in both parts all increase.For a specific RSthin-walled square tube with a cross section of 56.36mm×56.36mm and a thickness of 1.0mm,the yield stresses of its four ridge lines are strengthened to as high as 4 times of original ones,so its energy absorbed can increase by about 30%theoretically,though the proportion of four strengthened ridges is only 9.09%of the total perimeter of cross section,demonstrating that RScan significantly enhance the energy absorption capability of thin-walled square tube in axial crushing.
thin-walled square tube; ridgeline strengthening; energy absorption; mean crushing force;numerical simulation;crushing test
10.19562/j.chinasae.qcgc.2017.11.005
?國家863計劃項目(2012AA111302)資助。
原稿收到日期為2015年9月15日,修改稿收到日期為2016年12月28日。通信作者:朱西產(chǎn),教授,E-mail:zhuxichan@ #edu.cn。