您的當(dāng)前位置: 首頁 > 鍛壓文獻(xiàn)庫 > 沖壓計算機(jī)&自動化&安環(huán) > 正文|
1引言 板料成形是材料加工成形的一個重要分支，它廣泛應(yīng)用于汽車、航天、航空、家電等各個部門。隨著現(xiàn)代工業(yè)的發(fā)展，板料成形件越來越復(fù)雜，人們對板料成形的質(zhì)量和速度的要求也越來越高。傳統(tǒng)的板料成形模具的設(shè)計依賴的是經(jīng)驗(yàn)和直覺，并且通過反復(fù)試驗(yàn)調(diào)試來保證成形的質(zhì)量。這不僅需要消耗大量的人力物力，而且周期長，效率低，不能適應(yīng)社會發(fā)展的需要。上世紀(jì)七十年代以來，人們逐漸以數(shù)值模擬技術(shù)為輔助設(shè)計手段，大大降低了生產(chǎn)制造的成本。然而，由于板料成形是一種復(fù)雜的力學(xué)過程，其中包含幾何非線性、材料非線性、接觸非線性等強(qiáng)非線性問題，影響的參數(shù)非常多，這對數(shù)值模擬技術(shù)造成了極大的挑戰(zhàn)。雖然目前板料成形的數(shù)值模擬軟件已經(jīng)商業(yè)化，但板料成形的模擬技術(shù)還不夠完善，仍然是國內(nèi)國外研究的熱點(diǎn)。本文將主要介紹金屬板料成形數(shù)值模擬的研究現(xiàn)狀。 2板料成形有限元算法 用于板料成形的有限元算法大體可以分為彈-(粘)塑性和剛-(粘)塑性。粘塑性有限元法主要應(yīng)用于熱加工，而剛塑性有限元法在板料成形中應(yīng)用有限；目前，彈塑性有限元法在板料成形數(shù)值模擬中應(yīng)用較廣。用彈塑性有限元法分析板料成形問題，不僅能計算工件的變形和應(yīng)力、應(yīng)變分布；而且還能計算工件的回彈和殘余應(yīng)力、殘余應(yīng)變等。由于板料成形過程中板料與模具具有相對滑動、粘著和脫落，所以必須控制增量步的大小從而盡量反映真實(shí)情況。 根據(jù)對時間積分方法的不同，板料成形有限元算法可以分為靜力隱式、靜力顯式和動力顯式。隱式算法是非條件穩(wěn)定的，它在解決低速接觸問題中更有優(yōu)勢[1]，而在解決復(fù)雜的三維模型時將會遇到困難[2]：當(dāng)時間步長減小時，內(nèi)存消耗會急劇增大，甚至造成收斂問題；局部的不穩(wěn)定性很難達(dá)到力的平衡，這也不符合靜態(tài)隱式的先決條件。顯式算法克服了隱式算法的缺點(diǎn)[3]，然而它的不足之處在于，在解決像板料成形這樣的條件穩(wěn)定問題時，必須盡量消除慣性力的影響。對此一般有兩個辦法可采用：一個就是將運(yùn)動能限制在應(yīng)變能的5%以下；另外就是限制元素類型，一般只采用四節(jié)點(diǎn)的四邊形或者八節(jié)點(diǎn)的實(shí)體型。由于回彈是一個準(zhǔn)靜態(tài)問題宜采用隱式算法。板料成形中常常先用顯式算法模擬成形階段，而用隱式算法模擬回彈。 Belytschco[4]等人提出并行混合時間算法，即將兩種方法——顯式方法和隱式方法綜合起來，并行處理算法的一大難點(diǎn)就是進(jìn)程間的交互。研究表明[5]并行混合時間算法可以用來進(jìn)行板料成形的模擬，但是當(dāng)計算大的工件時，時間步長增大了，計算精度也不夠高；當(dāng)時間步長變小時，并行程度就會下降。 Kwansoo Chung[6，7]等將形變理論與極值功路徑相結(jié)合，提出了理想形變的思想，蘭箭等人[8]編寫了一步法有限元程序，并將其成功地應(yīng)用于板料成形的數(shù)值模擬中。 3板料成形接觸處理 板料成形過程是板料與模具不斷的接觸、作用與分離的過程，所以接觸問題是板料成形數(shù)值模擬中關(guān)鍵技術(shù)之一。 要解決接觸問題首先必須對模具進(jìn)行描述，模具的一般表示方法有：解析函數(shù)法、參數(shù)曲面法、網(wǎng)格法等。解析函數(shù)法只能用來表示類似于圓筒件等簡單的模型，而參數(shù)曲面法雖然能較準(zhǔn)確的表達(dá)模具曲面，但也存在算法復(fù)雜效率低的問題。網(wǎng)格法在一定程度上克服了以上兩者的缺點(diǎn)，加上有限元研究對象也用網(wǎng)格法表示，網(wǎng)格法成為大多數(shù)研究者采用的方法。陳中奎[9]提出了接觸特征空間的概念，并建立了接觸特征空間和模具曲面網(wǎng)格的對應(yīng)關(guān)系，成功地解決了三維板料沖壓成形過程有限元分析中的邊界條件非線性問題。 對于幾何形狀復(fù)雜的模具往往采用曲面單元描述法[10]，也就是將參數(shù)曲面描述法和有限元網(wǎng)格描述法相結(jié)合，即首先利用CAD等造型軟件生成模具的參數(shù)曲面，再通過IGES讀入FEA系統(tǒng)，并沿參數(shù)方向?qū)?shù)曲面劃分為一系列曲面片，從而得到由曲面單元描述的模具。王曉林等[10]在此基礎(chǔ)上采用逐級更新的搜索算法，推導(dǎo)出求迭代初值優(yōu)化選擇方案及滿足無穿透原則的判斷準(zhǔn)則，并建立了板料成形的數(shù)值模擬接觸搜索模型，這種模型不但改善了數(shù)值模擬的收斂性，而且提高了計算效率和穩(wěn)定性。 目前求解接觸問題常用的算法有拉格朗日乘子法和罰函數(shù)法，它們的共同點(diǎn)在于將有約束條件的變分問題轉(zhuǎn)化為無約束條件的變分問題。拉格朗日乘子法模擬出的結(jié)果更為準(zhǔn)確，但是計算效率較差，而且對于變形大的單元容易造成收斂困難。罰函數(shù)控制方程的階數(shù)和帶寬都小于拉格朗日乘子法，但是罰函數(shù)法的因子的取值對計算結(jié)果的精度影響很大。Shimizu和Sano[11，12]將罰函數(shù)發(fā)展為可以同時處理接觸問題和摩擦問題，其主要思想就是將由于接觸壓力造成的罰因子和由于摩擦造成的罰因子加起來作為總的罰因子來進(jìn)行搜索計算。 4摩擦處理 板料成形過程是板料相對于模具運(yùn)動變形的過程，所以板料成形的結(jié)果必定受到板料與模具之間摩擦力的影響。研究表明摩擦力與以下一些參數(shù)有關(guān)：接觸壓力、滑動速度、鋼板和模具的材料特性、表面粗糙度、潤滑劑和同向變形。尤其是板料厚度/板料面積很小時，摩擦力對板料的流動和最后的應(yīng)變分布有重要影響。除此之外[13]，研究還表明摩擦力的大小還與模具的幾何特性有關(guān)。因此對摩擦力描述的正確與否直接影響到數(shù)值模擬的精確性。傳統(tǒng)的用于有限元模擬的摩擦力模型有兩種：一是Amonton準(zhǔn)則，也就是庫侖準(zhǔn)則τ=μq，q是正應(yīng)力；另一個是固定摩擦(constant friction)τ=mk，k是剪應(yīng)力。一般情況下，正應(yīng)力比剪應(yīng)力大得多，所以第二種模型運(yùn)用得比較多，但是它們只有在主應(yīng)力值不是很大時模擬效果比較準(zhǔn)確。1976年Wanheim和Bay[14]提出了一個通用的模型用來描述摩擦力：在正應(yīng)力比較低(q/σ0<1.5，)時摩擦力正比于正應(yīng)力，而當(dāng)正應(yīng)力比較高(q/σ0>3)時摩擦力為一常數(shù)，中間為過渡區(qū)間。后來陸續(xù)有人對此模型進(jìn)行了大量的研究，絕大多數(shù)的研究集中在小變形的情況。S.B.Petersen[15]等人將其運(yùn)用于大塑性變形的情況，綜合運(yùn)用模擬和實(shí)驗(yàn)的方法對三種模型進(jìn)行了比較，研究結(jié)果表明Wanheim-Bay模型能夠更加精確地反映實(shí)際情況。也有研究[16]表明基于局部接觸條件的摩擦模型，特別是在板料高速滑動時，比傳統(tǒng)的庫侖模型更有優(yōu)勢。趙振鐸等人[17]從微觀角度出發(fā)，指出在變形過程中，真實(shí)接觸面積只占名義接觸面積很小的一部分，并對板料與模具相對滑動時的摩擦力學(xué)模型進(jìn)行了初步的探討。也有人[18]運(yùn)用神經(jīng)網(wǎng)絡(luò)的原理對摩擦行為研究。有關(guān)數(shù)值模擬之摩擦模型的研究還有待進(jìn)一步深入。 5成形極限圖FLD(forming limit diagram) 在上世紀(jì)六十年代，Keeler首先提出成形極限圖的思想，成形極限圖的兩個數(shù)據(jù)軸分別代表主應(yīng)變和次應(yīng)變，連接材料發(fā)生縮頸或斷裂時對應(yīng)的應(yīng)變狀態(tài)所得的點(diǎn)，就得到成形極限圖，起初它只局限于次應(yīng)變大于零的情況，后來由Goowin發(fā)展了次應(yīng)變小于零的情況。雖然FLD的原理比較簡單，但是其制作方法卻很復(fù)雜，比較常用的是網(wǎng)格應(yīng)變分析法，即將試樣首先劃分為一定的網(wǎng)格，根據(jù)樣品失效后網(wǎng)格的形變求出應(yīng)變值，從而作出成形極限圖。這種成形極限圖的缺陷在于：它依賴于材料的應(yīng)變軌跡，對工業(yè)生產(chǎn)中材料的非線性應(yīng)變軌跡，則不能用來預(yù)測成形。R.Arrieux[19]提出了基于應(yīng)力而作的成形極限圖即FLSD(Forming limit stress diagram)，它不受材料應(yīng)變軌跡的影響。R.Arrieux分別將其應(yīng)用于研究各向同性材料和各向異性材料，都得到不錯的效果。它的優(yōu)點(diǎn)在于，它可以用來預(yù)測復(fù)雜及多道次工序零件的縮頸的發(fā)生，可以分段表示出應(yīng)變狀況，這是經(jīng)典的成形極限圖所做不到的。Z.Zimniak[20]將FLD與MARC結(jié)合起來，數(shù)值模擬和物理實(shí)驗(yàn)證明基于攝動理論形成的FLSD能夠比較精確的預(yù)測板料成形過程的開裂問題。 FLD可以由實(shí)驗(yàn)得方法測得，但是測量的周期比較長，而且無論哪種極限圖都必須受到實(shí)驗(yàn)條件的影響，人們一方面尋找盡可能方便簡潔的實(shí)驗(yàn)方法來制作成形極限圖，另一方面也在積極尋求從理論上來對FLD進(jìn)行預(yù)測，這就必須知道FLD的影響因素。實(shí)驗(yàn)研究結(jié)果表明，F(xiàn)LD受潤滑劑、板料彎曲率，厚度以及落料方向的影響。材料性質(zhì)也影響著FLD的位置，曲線在坐標(biāo)軸上的截距和它的斜率取決于n和r值；不僅如此，預(yù)應(yīng)變和板料厚度也能改變FLD圖。D.W.A.Rees[21]對塑性鍍鋅鋼的FLD進(jìn)行保守的估算，其中考慮了板料落料方向、預(yù)應(yīng)變、厚度以及n和r的值，研究表明：n值必須相對較高，預(yù)應(yīng)變要盡量消除，落料方向?yàn)?5度較好。為了判別材料成形狀態(tài)，人們又先后發(fā)展了幾種不同的屈服準(zhǔn)則，其中具有代表意義的有von Mises模型、Hill(1948、1979、1993)模型、Gotoh模型以及Barlat-Lian模型等等[22]。 隨著計算機(jī)技術(shù)的發(fā)展，人們寄希望于利用數(shù)值模擬的方法來制作FLD，Marciniak和Kuczynski在上世紀(jì)60年代第一次對FLD進(jìn)行模擬，即M-K模型，實(shí)驗(yàn)證明由這種模型作出的FLD比較精確[23]。為了得到更好的效果Choi等人和Pishbin以及Gillis提出了另外一種數(shù)值模型稱為JG方法，但是它的用戶界面友好性不如M-K模型[23]。建立更好的數(shù)值模型仍然是FLD的研究方向之一。 用數(shù)值模擬的方法結(jié)合實(shí)驗(yàn)數(shù)據(jù)研究FLD的影響因素以及優(yōu)化屈服準(zhǔn)則也是FLD的熱點(diǎn)之一。Antonio F.Ávila[23]等人在前人的基礎(chǔ)上基于M-K模型來預(yù)測FLD的算法，其中用到了五種不同的屈服準(zhǔn)則：von Mises、Hill(1948)、Hill(1979)、Logan和Hosford、Hill(1993)；FLD隨著判據(jù)的不同而有顯著的不同，在此與實(shí)驗(yàn)結(jié)果最為符合的是用Hill(1993)準(zhǔn)則。 6缺陷的研究 板料成形的缺陷有起皺、開裂和回彈。 6.1起皺 在薄板沖壓成形過程中當(dāng)切向壓應(yīng)力達(dá)到或超過板料的臨界應(yīng)力時，就可能發(fā)生起皺。實(shí)際沖壓成形時，一般常采用調(diào)節(jié)壓邊力和設(shè)置拉伸筋的方法控制拉深過程中沖壓件突緣與側(cè)壁起皺，但由此又會導(dǎo)致板材成形極限下降，使沖壓件過早地發(fā)生破裂。特別是對于形狀復(fù)雜或深度大的沖壓件來說，既要保證成形形狀或深度又要抑制起皺發(fā)生，在實(shí)際加工工藝中常常是比較困難的。當(dāng)起皺與破裂兩者發(fā)生矛盾時，起皺的抑制必須以板材不破裂為基本條件，所以防止起皺的難點(diǎn)在于抑制起皺成形規(guī)律的強(qiáng)有力工具。由于起皺經(jīng)常發(fā)生在汽車的覆蓋件上，形成波紋或折疊，極大地影響了表面質(zhì)量，因此近年來對起皺的研究，特別是理論研究越來越多。 起皺的形成可以由Hill的分叉準(zhǔn)則(Hutchinson將其特別用于薄壁薄殼的情況)來預(yù)測，可以假設(shè)一個局部的分歧位移區(qū)，在這個位移區(qū)內(nèi)忽略接觸條件和連續(xù)性條件，就可以推導(dǎo)出只決定于材料性能和幾何性質(zhì)的臨界應(yīng)力。以主應(yīng)力為坐標(biāo)軸的起皺極限圖就是依據(jù)這種方法制得的。研究表明[24]應(yīng)力狀態(tài)的各向異性以及進(jìn)入側(cè)壁的板料的曲率對起皺的影響最大，各向異性指數(shù)增大時，板料發(fā)生起皺時拉伸距離也變大。起皺還與坯料的長寬比有關(guān)，也與沖頭沖壓的距離有關(guān)[25]。 用于起皺的有限元方法有兩種：一種是具有完美結(jié)構(gòu)的分叉分析，另外一種是具有初始非理想度的非分叉分析。對于第一種方法，它能夠非常精確地預(yù)測起皺現(xiàn)象，但是分析過程復(fù)雜，而且需要一定的起皺判據(jù)；第二種方法的精確度受到初始非理想度的影響。Xu Weili[26]采用了一種簡單的模擬起皺現(xiàn)象的方法：在接觸算法中采用節(jié)點(diǎn)直接投影以及使用BHG(blankholder gap)，這種方法避免使用起皺判別準(zhǔn)則，不需要復(fù)雜的塑性分叉理論以及坯料的初始非理想度。謝暉等人[27]從16個單元區(qū)域的起皺分析入手，用能量法計算，初步得到各單元的起皺臨界因子，然后進(jìn)一步進(jìn)行區(qū)域搜索或選定區(qū)域計算，找到包含某單元最可能發(fā)生起皺的區(qū)域，得到能反映各單元起皺發(fā)生難易程度的起皺臨界因子，并用云圖將它描述出來，從而較精確、直觀地預(yù)測起皺的發(fā)生。 對不均勻拉應(yīng)力下的起皺機(jī)理、起皺規(guī)律、影響起皺高度的因素、消皺措施等各國都進(jìn)行了大量研究工作。但對剪應(yīng)力下的起皺的研究還很少。崔令江[28]研究表明：(1)剪應(yīng)力起皺實(shí)驗(yàn)中試樣的起皺區(qū)長度越大，其抗剪應(yīng)力起皺能力就越差；(2)板材的硬化指數(shù)n值和厚向異性系數(shù)r值越大，抗剪應(yīng)力起皺的能力越強(qiáng)。 6.2開裂 開裂是由于材料的強(qiáng)度或塑性不足，當(dāng)拉應(yīng)力超過臨界值時便會發(fā)生。開裂發(fā)生的位置主要在：凸模端部、側(cè)壁、凸模圓角部位、法蘭部分和拉伸筋部分[29]。由于開裂的影響因素也很多，因此到目前為止還沒有十分精確的判斷開裂的準(zhǔn)則。對于開裂的判斷現(xiàn)在用的比較多的主要還是成形極限圖[30]：包括FLCN(forming limit curves at neck)以及FLCF(forming limit curves at fracture)；最大拉深率也是主要判據(jù)之一，Min Wan等人[31]結(jié)合了圓柱和圓錐杯形的內(nèi)在聯(lián)系，從理論上給出了圓錐杯形的極限拉深系數(shù)的判斷公式。 在有限元模擬上，人們一方面在探索新的模型，另一方面將現(xiàn)有的判斷準(zhǔn)則與有限元模擬結(jié)合起來進(jìn)行研究。K.Komori[32]提出了一種節(jié)點(diǎn)分離的模型：將服從開裂準(zhǔn)則的單元體和與它鄰近的單元體共有的節(jié)點(diǎn)分離開來，也就是假設(shè)開裂沿著單元體的邊界擴(kuò)展；而一般的有限元模型中常常將服從開裂準(zhǔn)則的單元直接刪除，在一定程度上改變了材料本來的特性，導(dǎo)致模擬結(jié)果誤差相對較大。Z.H.Chen等人[33]應(yīng)用了混合有限元法(Mixed finite element method)，混合了位移和壓力的有限元法，可以對板料成形中的不可壓縮的塑性變形進(jìn)行比較精確的模擬。而Ridha Hambli等人[34]用了一種反向技術(shù)的模型，可以識別開裂準(zhǔn)則的臨界值，從而將其應(yīng)用于開裂模擬中，對研究開裂準(zhǔn)則和開裂模擬都有很大的幫助。 另外，由于摩擦力在板料成形過程常常充當(dāng)為拉應(yīng)力的一部分，所以摩擦和潤滑也是研究開裂問題的主要方向之一。 6.3回彈 回彈是板料成形后不可避免的現(xiàn)象，回彈現(xiàn)象主要表現(xiàn)為整體卸載回彈、切邊回彈和局部卸載回彈，當(dāng)回彈量超過允許容差后，就成為成形缺陷，影響零件的幾何精度。因此，回彈一直是影響、制約模具和產(chǎn)品質(zhì)量的重要因素。早期的工作主要基于解析法對一些簡單零件純彎曲或拉彎成形的回彈進(jìn)行分析[35]。Gardiner最早基于彎曲工程理論模型對理想彈塑性板彎曲的回彈問題進(jìn)行了研究[36]。Hill[37]首先建立了板彎曲的精確數(shù)學(xué)理論，并對剛塑性寬板純彎曲問題進(jìn)行了研究。現(xiàn)在人們逐漸尋找一些新的研究回彈問題的模型，這樣能夠研究更多種因素。Lumin Geng[38]采用了一種新的各向不等硬化模型，并且將這種硬化模型和四個屈服公式用在了ABQUS中：von Mises、Hill二次方程式、Barlat三參量和Barlat1996。 由于應(yīng)用于板料成形的有限元方法的不斷發(fā)展，人們逐漸對有限元模擬回彈進(jìn)行了研究。由于回彈是成形的最后一步，成形過程模擬中產(chǎn)生的任何誤差都會積累到回彈計算階段；因此，回彈模擬結(jié)果的準(zhǔn)確性很大程度上取決于成形過程的模擬精度。回彈問題的有限元模擬受到許多因素的影響：厚度方向積分點(diǎn)的數(shù)量、元素類型、網(wǎng)格大小、與凹模臺肩處每個元素接觸角度、慣性和接觸算法的影響。不僅如此回彈還受到很多物理參量的影響：材料性質(zhì)、硬化規(guī)則、摩擦系數(shù)、壓邊力(BHF)和卸載過程。K.P.Li[39]等人以及Luc Papeleux等人[40]利用有限元模擬的方法對影響板料回彈的物理參數(shù)和數(shù)值因素進(jìn)行了詳細(xì)的研究。Jenn-Terng Gau等人[41]在研究回彈問題時對包申格效應(yīng)進(jìn)行了不少研究。研究包申格效應(yīng)普遍用單軸拉伸/壓縮實(shí)驗(yàn)和循環(huán)扭曲實(shí)驗(yàn)。Jenn-Terng Gau等人研究發(fā)現(xiàn)兩種相似的樣品有相同的總的應(yīng)變，卻有完全不同的回彈效果。這是因?yàn)樵趹?yīng)變區(qū)它們變形的歷程不一樣，包申格效應(yīng)不僅與材料有關(guān)還與變形的歷程有關(guān)。包申格效應(yīng)對鋼類材料影響不大對鋁材影響較大。P.Xue等人[42]基于膜殼理論和能量方法研究發(fā)現(xiàn)彎曲半徑越大回彈越大。 在有限元算法上普遍采用顯式算法計算成形過程，再用隱式算法計算回彈。由于影響因素較多，過程復(fù)雜，今后回彈仍將是板料成形研究的主攻方向之一。 7展望 板料成形的數(shù)值模擬經(jīng)過20多年的研究已經(jīng)有了很大的發(fā)展。在實(shí)際生產(chǎn)中往往由于零件所要求的變形程度超過了材料的一次成形所允許的最大變形程度，必須采用多道次拉深成形的方法。如何處理多道次拉深變形的分配，以及如何精確確定材料首次拉深后的性能參數(shù)是板料成形方面亟待解決的問題。另外由于板料成形模擬技術(shù)在汽車覆蓋件方面的應(yīng)用，如何精確模擬復(fù)雜曲面的成形過程也是大家正在研究的熱點(diǎn)之一。今后板料成形的數(shù)值模擬研究還要注意以下幾方面： 1)進(jìn)一步優(yōu)化時間積分方案，能夠?qū)追N時間積分方案綜合起來運(yùn)用，各取所長。2)繼續(xù)對摩擦行為進(jìn)行研究，建立更合適的接觸算法。3)優(yōu)化各種本構(gòu)關(guān)系，能夠?qū)⒊尚螛O限圖和有限元分析更好地結(jié)合起來。4)對缺陷問題還要深入研究，優(yōu)化數(shù)值模擬模型。5)將有限元模擬系統(tǒng)和神經(jīng)網(wǎng)絡(luò)、自動控制等結(jié)合起來，形成大的分析系統(tǒng)。 |