郭鵬程 葉 拓 曹淑芬 徐從昌 李落星, 2
1.湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙,4100822.湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院,長(zhǎng)沙,410082
高應(yīng)變速率下鑄態(tài)AM80鎂合金的變形行為及數(shù)值模擬
郭鵬程1,2葉 拓1,2曹淑芬1徐從昌1,2李落星1, 2
1.湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙,4100822.湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院,長(zhǎng)沙,410082
AM80鎂合金;動(dòng)態(tài)沖擊;應(yīng)變速率敏感性;Johnson-Cook本構(gòu)方程
研究表明:汽車油耗與自重成線性關(guān)系,在其他條件不變的前提下,若汽車自重降低10%,則油耗至少可降低5%,因此,減輕汽車自重、降低能耗成為各大汽車生產(chǎn)商提高競(jìng)爭(zhēng)力的關(guān)鍵[1]。在汽車輕量化的浪潮中,雖然鋼仍然占主導(dǎo)地位,但以輕質(zhì)結(jié)構(gòu)材料為代表的鎂合金,由于密度低、比強(qiáng)度和比剛度高,機(jī)加工、阻尼減振和電磁屏蔽性能優(yōu)異等,是實(shí)現(xiàn)汽車輕量化最理想的結(jié)構(gòu)材料。鎂合金零部件在汽車上的應(yīng)用可以降低汽車自重和油耗,提高加速性能和駕駛舒適度。在最近十年,鎂合金在汽車工業(yè)中的應(yīng)用得到了很大程度的發(fā)展。
鎂合金結(jié)構(gòu)件在實(shí)際服役過程中,除了承受靜態(tài)載荷之外,還要承受高速?zèng)_擊載荷(如沖擊、碰撞環(huán)境)[2]。這種高速?zèng)_擊載荷導(dǎo)致構(gòu)件在極短的時(shí)間內(nèi)失效,具有不可預(yù)見性,且后果極其嚴(yán)重[3]。隨著鎂合金在汽車上應(yīng)用水平的不斷提高,尤其是鎂合金在汽車前端等部件上的使用,使得鎂合金在沖擊載荷作用下的斷裂規(guī)律及失效形式顯得尤為重要。因此,研究鎂合金在高速?zèng)_擊載荷下的動(dòng)態(tài)響應(yīng)和斷裂失效行為對(duì)有效防止鎂合金在高速?zèng)_擊過程中的斷裂,提高鎂合金構(gòu)件在汽車和航空航天中的可靠性和使用壽命,充分發(fā)揮鎂合金材料的潛在塑性具有重要意義[4]。
AM系鎂合金是一種新型鎂合金,它不僅具有優(yōu)良的鑄造性能,而且其伸長(zhǎng)率、韌性和抗沖擊能力也相對(duì)較好,且價(jià)格低廉,是實(shí)現(xiàn)鎂合金工業(yè)應(yīng)用最具潛力的鎂合金系之一[5]。然而,目前國(guó)內(nèi)外關(guān)于鎂合金動(dòng)態(tài)變形行為及其變形機(jī)制的研究仍較少且不系統(tǒng),主要集中在Mg-Al-Zn系合金上,尤其是AZ31鎂合金[2-4, 6-9]。對(duì)于AM系鎂合金,主要是研究其不同溫度下的靜態(tài)壓縮變形[10-12],關(guān)于動(dòng)態(tài)沖擊載荷下的力學(xué)行為及變形機(jī)制的研究很少。本文對(duì)AM80鎂合金進(jìn)行室溫下不同應(yīng)變速率時(shí)的單向壓縮,研究合金高應(yīng)變速率下的變形行為、力學(xué)本構(gòu)關(guān)系及其仿真模擬。
配制AM80鎂合金的原料化學(xué)成分(質(zhì)量分?jǐn)?shù))分別為:鎂99.98%、鋁99.6%、電解錳95%。保護(hù)氣體采用四氟乙稀(C2F4)+二氧化碳(CO2)的混合氣體。熔煉合金時(shí)對(duì)溶劑錳進(jìn)行阻燃處理,合金經(jīng)精煉之后,在720 ℃下澆鑄到預(yù)熱至250 ℃的金屬模(φ90 mm×500 mm)中,冷卻凝固后獲得實(shí)驗(yàn)所需的AM80鎂合金。其主要化學(xué)成分(質(zhì)量分?jǐn)?shù))為:Al 8.0%,Mn 0.1%,Mg余量。鑄錠在450 ℃下均勻化處理16 h后空冷,最后切割成φ10 mm×10 mm的圓柱形壓縮試樣,取樣位置和晶粒尺寸如圖1所示。
(a)取樣位置
(b)金相顯微組織圖1 取樣位置和金相顯微組織Fig.1 Schematic of samples cutting from the as-received alloys and optical microstructure
準(zhǔn)靜態(tài)壓縮實(shí)驗(yàn)在INSTRON-4206電液伺服實(shí)驗(yàn)機(jī)上進(jìn)行,應(yīng)變速率為0.0001 s-1;動(dòng)態(tài)壓縮實(shí)驗(yàn)采用分離式Hopkinson壓桿(split Hopkinson pressure bar, SHPB)進(jìn)行,應(yīng)變速率分別約為800 s-1、1050 s-1、1600 s-1、1850 s-1、2100 s-1。實(shí)驗(yàn)環(huán)境均為室溫,壓縮前在試樣表面和壓頭間涂潤(rùn)滑劑,以減小試樣和壓頭間的摩擦。壓縮變形后采用火花放電線切割機(jī)沿試樣軸向切取金相試樣。金相試樣經(jīng)打磨、拋光后用5 g苦味酸+5 g冰醋酸+10 mL蒸餾水+80 mL無水乙醇的侵蝕劑腐蝕,通過光學(xué)顯微鏡(optical microscope,OM)觀察其金相組織。
2.1 不同應(yīng)變速率下的力學(xué)響應(yīng)
圖2 鑄態(tài)AM80鎂合金的真應(yīng)力-真應(yīng)變曲線Fig.2 Compressive true stress-strain curves of the casting AM80 alloy
2.2 材料本構(gòu)模型
Johnson-Cook (J-C)模型主要考慮應(yīng)變速率效應(yīng)和溫度效應(yīng),是金屬材料在高應(yīng)變速率下最通用的本構(gòu)模型。本文采用J-C模型來描述鑄造固溶態(tài)AM80鎂合金在高應(yīng)變速率下的力學(xué)行為,并忽略溫度的影響,其本構(gòu)模型表達(dá)式如下:
(1)
(2)
圖3所示為擬合數(shù)據(jù)與實(shí)驗(yàn)結(jié)果的對(duì)比,可以看出,修正參數(shù)C后的J-C本構(gòu)方程能夠很好地描述鑄造固溶態(tài)AM80鎂合金的動(dòng)態(tài)壓縮應(yīng)力-應(yīng)變行為。
圖3 J-C本構(gòu)擬合結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比Fig.3 Experimental and computational results of dynamic compression tests on the studied AM80 alloys
2.3 變形顯微組織
圖4 以0.0001 s-1的應(yīng)變速率壓裂后的顯微組織Fig.4 OM images for the fractured specimen at strain rate of 0.0001 s-1
圖4是以0.0001 s-1的應(yīng)變速率壓裂后試樣中心位置的金相顯微圖,其應(yīng)變約0.246。由圖4可知,此應(yīng)變下晶粒變形嚴(yán)重,部分晶粒已分辨不出晶界。晶粒內(nèi)產(chǎn)生了大量的形變孿晶,大部分形變孿晶與壓縮軸的夾角約60°。一次孿晶間產(chǎn)生了大量的平行或相交于一次孿晶的細(xì)小二次形變孿晶,這表明二次孿晶傾向于在一次孿晶附近形核生長(zhǎng)。此外,在晶界和孿晶界還產(chǎn)生了微小裂紋,且均沿其晶界或?qū)\晶界擴(kuò)展,這表明實(shí)驗(yàn)用鎂合金傾向于在晶界和孿晶界處形核并擴(kuò)展。
圖5是不同應(yīng)變速率下動(dòng)態(tài)壓縮后試樣中心位置的金相顯微圖。由圖5可知,動(dòng)態(tài)壓縮后晶粒內(nèi)產(chǎn)生了大量相互平行的形變孿晶。應(yīng)變速率為1850 s-1時(shí),其壓縮應(yīng)變約0.220,此時(shí)晶粒變形極為嚴(yán)重,晶粒已完全分辨不出晶界,初生形變孿晶間還產(chǎn)生了大量相交于初生孿晶的二次孿晶。與準(zhǔn)靜態(tài)壓縮(應(yīng)變約0.246)相比,該應(yīng)變速率下所產(chǎn)生的形變孿晶密度明顯更高,孿晶間距相對(duì)較小,表明實(shí)驗(yàn)用AM80鎂合金的孿生具有較強(qiáng)的正應(yīng)變速率敏感性。這些形變孿晶將晶粒分割成許多細(xì)小的網(wǎng)格,起到了一定的晶粒細(xì)化作用,從而導(dǎo)致相同應(yīng)變下應(yīng)力水平較高。正是由于孿生的正應(yīng)變速率敏感性,導(dǎo)致隨應(yīng)變速率增大,流變應(yīng)力增大。然而,當(dāng)應(yīng)變速率增大至2100 s-1時(shí),雖然此時(shí)應(yīng)變最大(約0.25),但晶粒內(nèi)的孿晶密度反而降低(圖5d),孿生的正應(yīng)變速率敏感性消失。研究結(jié)果表明,鎂合金的應(yīng)變硬化率與孿晶密度密切相關(guān)。以2100 s-1的應(yīng)變速率壓縮時(shí),由于局部溫升效應(yīng)增加,促進(jìn)了材料的滑移,從而導(dǎo)致孿生密度降低,故其流變應(yīng)力不再增大反而略有減小。
s-1
s-1圖5 不同應(yīng)變速率下動(dòng)態(tài)壓縮變形后的顯微組織Fig.5 OM images of the deformed specimen with strain rates
為了驗(yàn)證本文J-C本構(gòu)模型參數(shù)的有效性,采用顯式動(dòng)力學(xué)分析軟件ABAQUS對(duì)實(shí)驗(yàn)用鎂合金的SHPB實(shí)驗(yàn)進(jìn)行仿真模擬?;赟HPB實(shí)驗(yàn)裝置,將子彈、入射桿、透射桿和壓縮試樣簡(jiǎn)化成共軸的圓柱體。子彈長(zhǎng)度為200mm,入射桿與透射桿長(zhǎng)度為1400m,子彈和壓桿直徑均為14mm,試樣尺寸為φ10 mm×10 mm。計(jì)算單元使用8節(jié)點(diǎn)六面體單元3Dsolid。為兼顧仿真精度與效率,子彈、壓桿和試樣的單元尺寸分別選為1.5 mm、1.5 mm、0.1 mm。此外,為了縮短計(jì)算時(shí)間,仿真采用1/4有限元模型。仿真有限元網(wǎng)格模型如圖6所示,其單元與節(jié)點(diǎn)數(shù)見表1。
圖6 SHPB幾何與有限元模型Fig.6 Finite element mesh model of SHPB
表1 SHPB有限元模型單元與節(jié)點(diǎn)數(shù)
子彈和壓桿采用線彈性鋼材料模型,其密度為7.85×10-6kg/mm3,彈性模量為210 GPa,泊松比為0.3。實(shí)驗(yàn)用AM80鎂合金的動(dòng)態(tài)力學(xué)響應(yīng)行為采用J-C材料本構(gòu)模型來描述,其密度為1.82×10-6kg/mm3,彈性模量為45 GPa,泊松比為0.35,J-C本構(gòu)參數(shù)見式(2)。采用將實(shí)驗(yàn)所采集到的子彈入射速度直接加載至子彈的方式進(jìn)行仿真模擬。不同應(yīng)變速率下仿真與實(shí)驗(yàn)的真應(yīng)力-真應(yīng)變?nèi)鐖D7所示。
圖7 有限元仿真結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比Fig.7 Comparison among simulation and experimental results
由圖7可知,中低應(yīng)變下,數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合得較好;但在高應(yīng)變時(shí),仿真分析結(jié)果與實(shí)驗(yàn)結(jié)果差異比低應(yīng)變時(shí)明顯要大。這主要是因?yàn)橛糜诜抡娴谋緲?gòu)模型未考慮溫升軟化效應(yīng),而實(shí)驗(yàn)過程中材料的應(yīng)變率效應(yīng)和溫度效應(yīng)耦合在一起。低應(yīng)變時(shí),溫升效應(yīng)較小,即形變溫升軟化基本可以忽略,故仿真與實(shí)驗(yàn)的吻合度較高;而在高應(yīng)變時(shí),由于形變溫升軟化加劇,大部分應(yīng)變硬化和應(yīng)變率硬化被絕熱溫升軟所抵消,故在變形中后期仿真與實(shí)驗(yàn)的流變應(yīng)力差異較大,且流變應(yīng)力隨應(yīng)變與變形速率的增大而增大。
(2)基于準(zhǔn)靜態(tài)及動(dòng)態(tài)壓縮,確定了鑄造固溶態(tài)AM80鎂合金的室溫J-C本構(gòu)方程的各參數(shù),且該本構(gòu)方程擬合結(jié)果與實(shí)驗(yàn)基本吻合。
(3)建立了實(shí)驗(yàn)用鎂合金試樣的SHPB有限元模型。中低應(yīng)變下,數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合得較好;高應(yīng)變時(shí),仿真分析結(jié)果與實(shí)驗(yàn)結(jié)果的差異比低應(yīng)變時(shí)明顯要大。研究結(jié)果為金屬高應(yīng)變率力學(xué)行為的研究提供了重要的手段和理論指導(dǎo)。
[1] MILLER W S, ZHUANG L, BOTTENMA J, et al. Experimental Study on Tensile Property of AZ31B Magnesium Alloy at Different High Strain Rates and Temperatures [J]. Materials Science and Engineering A, 2000, 280: 37-49.
[2] FENG Fei, HUANG Shangyu, MENG Zhenghua, et al.Experimental Study on Tensile Property of AZ31B Magnesium Alloy at Different High Strain Rates and Temperatures [J]. Materials & Design, 2014, 57(5): 10-20.
[3] YANG Suyuan. The Microstructure Features and the Deformation Mechanism of a Fine Grained Magnesium Alloy under Dynamic Loading [J]. Reviews on Advanced Materials Science, 2010, 25(2): 122-127.
[4] MAO Pingli, LIU Zheng, WANG Changyi. Texture Effect on High Strain Rates Tension and Compression Deformationbehavior of Extruded AM30 Alloy [J]. Materials Science and Engineering A, 2012, 539(2): 13-21.
[5] 許春香, 鞠輝, 張志偉. Sr含量對(duì)AM80鎂合金顯微組織和力學(xué)性能的影響[J]. 中國(guó)有色金屬學(xué)報(bào), 2013, 23(2): 349-355. XU Chunxiang, JU Hui, ZHANG Zhiwei. Effect of Sr Content on Microstrure and Mechanical Properties AM80 Magnesium Alloy [J]. The Chinese Journal of Nonferrous Metals, 2013, 23(2): 349-355.
[6] WAN G, WU B, ZHANG Y, et al. Anisotropy of Dynamic Behavior of Extruded AZ31 Magnesium Alloy [J]. Materials Science and Engineering A, 2010, 527(12): 2915-2924.
[7] ULACIA I, DUDAMELL N V, GALVEZ F, et al. Mechanical Behavior and Microstructural Evolution of a Mg AZ31 Sheet at Dynamic Strain Rates [J]. Acta Materialia, 2010, 58(8): 2988-2998.
[8] WU B, ZHANG Y, WAN G, et al. Primary Twinning Selection with Respect to Orientation of Deformed Grains Inultra-rapidly Compressed AZ31 Alloy [J]. Materials Science and Engineering A, 2012, 541: 120-127.
[9] DUDAMELL N V, ULACIA I, GALVEZ F, et al. Influence of Texture on the Recrystallization Mechanisms in an AZ31 Mg Sheet Alloy at Dynamic Rates[J].Materials Science and Engineering A, 2012, 532(1): 528-535.
[10] MA Q, KADIRII H E, OPPEDAL A L, et al. Twinning and Double Twinning upon Compression of Prismatic Textures in an AM30 Magnesium Alloy[J].Scripta Materialia, 2011, 64(9): 813-816.
[11] MA Q, OPPEDAL A L, BAIRD J C, et al. Twinning Effects in a Rod-textured AM30 Magnesium Alloy [J]. International Journal of Plasticity, 2012, 29(1): 60-76.
[12] BISWAS S, SUWAS S, SIKAND R, et al. Analysis of Texture Evolution in Pure Magnesium and the Magnesium Alloy AM30 during Rod and Tube Extrusion [J]. Materials Science and Engineering A, 2011, 528(10): 3722-3729.
(編輯 陳 勇)
Deformation Behavior and Numerical Simulation of a Casting AM80 Magnesium Alloy under High Strain Rate Loading
GUO Pengcheng1,2YE Tuo1,2CAO Shufen1XU Congchang1,2LI Luoxing1,2
1.State Key Laboratory of Advanced Design and Manufacture for Vehicle Body,Hunan University,Changsha,410082 2.College of Mechanical and Vehicle Engineering,Hunan University,Changsha,410082
AM80 magnesium; dynamic impact; strain rate sensitivity; Johnson-Cook constitutive equation
2016-05-27
國(guó)家科技重大專項(xiàng)(2014ZX0400207);國(guó)家自然科學(xué)基金資助項(xiàng)目(51475156)
TG156
10.3969/j.issn.1004-132X.2017.06.017
郭鵬程,男,1985生。湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院博士研究生。主要研究方向?yàn)檎嚺鲎舶踩败嚿碛娩X鎂合金的動(dòng)態(tài)沖擊變形行為。E-mail:gpch860429@163.com。葉 拓,男,1989年生。湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院博士研究生。曹淑芬,女,1989年生。湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院碩士研究生。徐從昌,男,1989年生。湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院博士研究生。李落星(通信作者),男,1968年生。湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院教授、博士研究生導(dǎo)師。E-mail:llxly2000@163.com。