張晉輝,牛 婷,溫 凱,謝瑞珍

(晉中學(xué)院 機(jī)械學(xué)院,山西 晉中 030619)

0 引言

LF2M屬于典型的Al-Mg系防銹鋁合金，由于具有比強(qiáng)度高、耐磨性高、抗蝕性好、疲勞性能良好等優(yōu)點(diǎn)，已廣泛應(yīng)用于航空航天、武器裝備、石油化工等領(lǐng)域的重要部件中[1]。錐形件強(qiáng)旋是一個(gè)多因素交互作用下的復(fù)雜的局部連續(xù)塑性成形過(guò)程[2]，在實(shí)際生產(chǎn)中工藝參數(shù)的選擇主要依據(jù)經(jīng)驗(yàn)而定，同時(shí)該成形工藝的復(fù)雜性加深了仿真模擬的難度。本文根據(jù)具體的錐形件強(qiáng)旋變形特點(diǎn)[3-4]，在對(duì)建模過(guò)程涉及到的關(guān)鍵技術(shù)進(jìn)行探討與處理的基礎(chǔ)上，建立了可靠的符合實(shí)際的三維有限元模型，進(jìn)而為錐形件旋壓成形的工藝試驗(yàn)研究與參數(shù)確定及優(yōu)化提供了有效的參考。

1 錐形件強(qiáng)旋成形有限元建模的關(guān)鍵技術(shù)

1.1 錐形件強(qiáng)旋成形幾何模型

本文研究的錐形件坯料外徑D0=100 mm，內(nèi)徑D0′=20 mm，厚度t0=3 mm；旋輪圓角半徑r=6 mm，直徑Dρ=80 mm，錐角β=45°，寬度B=25 mm，進(jìn)給比f(wàn)=1.8 mm/r；芯模小端直徑D1=40 mm，大端直徑D2=160 mm，圓角半徑R=6 mm，錐角α=45°，轉(zhuǎn)速n=195 r/min。

為了使建立的三維模型簡(jiǎn)單適用，采用Tie方式將毛坯與芯模固定在一起，實(shí)現(xiàn)芯模與毛坯的完全同速運(yùn)轉(zhuǎn)，從而避免了尾頂?shù)慕９ぷ?；采用雙旋輪裝置，旋輪對(duì)稱(chēng)分布，旋輪沿芯模母線(xiàn)方向進(jìn)給。建立的錐形件強(qiáng)旋成形三維模型如圖1所示。

圖1 錐形件強(qiáng)旋成形三維模型

1.2 材料屬性設(shè)置

1.3 單元選取

芯模采用離散剛體，旋輪采用解析剛體進(jìn)行處理，以減少模擬的仿真運(yùn)算量和計(jì)算模擬時(shí)間。由于線(xiàn)性減縮積分單元能夠很好地承受扭曲變形，故毛坯采用六面體顯式線(xiàn)性減縮積分實(shí)體單元C3D8R進(jìn)行離散。為了便于模擬過(guò)程中有關(guān)應(yīng)力應(yīng)變等場(chǎng)變量的輸出和處理，在模型中定義了坯料的局部材料方向，即局部材料方向是隨著板料的變形而轉(zhuǎn)動(dòng)的，因此定義變形后1方向?yàn)榕髁夏妇€(xiàn)方向，2方向?yàn)榕髁现芟颍?方向?yàn)榕髁戏ㄏ?，如圖2所示。

圖2 坯料局部材料方向定義

1.4 網(wǎng)格劃分

本文采用各向同性材料，考慮到數(shù)值模擬的精確性和模擬的效率問(wèn)題，對(duì)坯料、芯模以及旋輪采用不同的單元進(jìn)行不同密度的網(wǎng)格劃分。在ABAQUS軟件中，坯料采用實(shí)體單元C3D8R進(jìn)行離散，為了模擬的精確性，對(duì)坯料進(jìn)行網(wǎng)格局部細(xì)化，在徑向、軸向壁厚和同向進(jìn)行單元網(wǎng)格劃分，整個(gè)工件共劃分為6 600個(gè)單元、9 200個(gè)節(jié)點(diǎn)。對(duì)于芯模和旋輪，由于在旋壓過(guò)程中理想的情況被認(rèn)為是不發(fā)生變形的，因此假設(shè)成剛體，網(wǎng)格的密度也相應(yīng)地比坯料要小得多，以提高計(jì)算效率。其中芯模視為R3D4殼單元?jiǎng)傮w進(jìn)行離散，共802個(gè)單元、824個(gè)節(jié)點(diǎn)；旋輪視為解析剛體進(jìn)行離散，無(wú)需網(wǎng)格的劃分。坯料和芯模的單元網(wǎng)格劃分分別如圖3、圖4所示。

圖3 坯料網(wǎng)格劃分

圖4 芯模網(wǎng)格劃分

1.5 邊界條件處理

1.5.1 接觸邊界條件處理

在ABAQUS中分別定義了坯料頂面接觸面(Blank Top)、坯料底面外部接觸面(Blank Bot Out)、坯料底面中心接觸面(Blank Bot Core)、芯模頂面接觸面(Mandrel Top Core)、芯模錐體接觸面(Mandrel Body)、右旋輪接觸面(Roller Right)、左旋輪接觸面(Roller Left)以及接觸對(duì)Roller Right-Blank Top、Roller Left-Blank Top、Blank Bot Out-Mandrel Body。由于模型中旋輪和芯模都視為剛體處理，故在以上的接觸對(duì)定義中將右旋輪接觸面、左旋輪接觸面、芯模接觸面都定義為主控表面，而坯料的所有接觸面都定義為從屬表面。

1.5.2 摩擦邊界條件處理

錐形件強(qiáng)旋成形過(guò)程中需要采用簡(jiǎn)化的摩擦力學(xué)模型，旋輪與坯料之間的摩擦當(dāng)成定摩擦來(lái)處理，而坯料與芯模之間的接觸摩擦相對(duì)要簡(jiǎn)單得多，理想的狀態(tài)是相對(duì)不發(fā)生摩擦，模型中視為無(wú)摩擦處理。因此，在錐形件強(qiáng)旋模型中定義了兩種摩擦：Fric和Nofric。Fric定義坯料與旋輪之間的摩擦系數(shù)為0.04。Nofric指忽略表面間的摩擦，表征坯料與芯模之間沒(méi)有摩擦。

2 模型可靠性驗(yàn)證

2.1 理論評(píng)估

在建模過(guò)程中，為了提高計(jì)算精度和效率，對(duì)不同單元類(lèi)型、網(wǎng)格密度以及質(zhì)量放大因子情況下的旋壓過(guò)程進(jìn)行了模擬分析。分析結(jié)果表明：當(dāng)坯料采用六面體減縮積分實(shí)體單元C3D8R、網(wǎng)格劃分情況如圖3所示(單元數(shù)6 600個(gè)單元，節(jié)點(diǎn)數(shù)9 200個(gè))、質(zhì)量放大因子為225時(shí)，模擬所獲得的偽應(yīng)變能、動(dòng)能與內(nèi)能的比很小(如圖5所示)，說(shuō)明本文所建立的模型是可靠的。

圖5 錐形件強(qiáng)旋過(guò)程能量變化

2.2 實(shí)驗(yàn)對(duì)比

為了進(jìn)一步驗(yàn)證本文所建仿真模型的可靠性，采用文獻(xiàn)[5]實(shí)驗(yàn)參數(shù)為模擬條件，對(duì)材料為L(zhǎng)F2M的坯料采用不同的偏離率進(jìn)行了模擬計(jì)算，以壁厚差為評(píng)價(jià)指標(biāo)，將計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比，如圖6所示。

圖6 計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

從圖6中可以看出：計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的壁厚差變化趨勢(shì)基本一致，且最大誤差不超過(guò)15%，表明本文所建立的三維有限元仿真模型是可行和可靠的。

3 結(jié)論

本文針對(duì)錐形件強(qiáng)旋變形特點(diǎn)，基于ABAQUS軟件平臺(tái)，在對(duì)幾何模型、單元選擇、材料屬性設(shè)置、網(wǎng)格劃分、邊界條件等關(guān)鍵技術(shù)探討與合理處理的基礎(chǔ)上建立了三維有限元模型，理論評(píng)估和實(shí)驗(yàn)驗(yàn)證證明該模型與實(shí)際成形過(guò)程十分吻合。