基于Voronoi圖理論的微體積成形數(shù)值建模及驗(yàn)證

席道英鄭偉

(山東建筑大學(xué)工程力學(xué)研究所，山東濟(jì)南250101)

0 引言

隨著科學(xué)技術(shù)的快速發(fā)展，人們對微型零件的需求正逐漸增加，微塑性成形技術(shù)成為材料加工領(lǐng)域的研究熱點(diǎn)。微塑性成形技術(shù)是指通過材料發(fā)生一定的形變，制造出二維或三維方向上尺寸＜1 mm零件的技術(shù)[1-4]，具有高生產(chǎn)效率、性能好、誤差小及無污染等優(yōu)點(diǎn)。微塑成形技術(shù)在微機(jī)械及微電子等領(lǐng)域具有十分重要的應(yīng)用前景，近幾十年在市場因素的推動下迅速發(fā)展成為材料加工領(lǐng)域的熱點(diǎn)。

數(shù)值模擬方法是研究微塑性成形的主要手段之一，在生產(chǎn)工藝中，可以減少工藝設(shè)計(jì)的時間和降低生產(chǎn)成本。目前提出的微塑性成形理論中的區(qū)域化模型不包含各種尺度參數(shù)，因此基于成形理論的數(shù)值模擬不能直接使用現(xiàn)有的有限元軟件。很多學(xué)者都基于實(shí)驗(yàn)數(shù)據(jù)推導(dǎo)本構(gòu)關(guān)系及研究微塑性成形過程中的尺度效應(yīng)現(xiàn)象，但缺少一些理論依據(jù)[5-6]。彭林法[7]考慮了表面層模型理論的情況，將表面層晶粒看作單晶，從而改進(jìn)了Hall-Petch公式，建立了一種復(fù)合材料區(qū)域化模型并實(shí)現(xiàn)了微塑性成形過程的數(shù)值模擬。但在微尺度下，變形區(qū)域只有幾個或十幾個晶粒時，不能忽視單個晶粒取向?qū)φ麄€變形的影響，因此需要對整個材料模型進(jìn)行晶?；幚恚员愀逦貙?shí)現(xiàn)對微塑性成形過程的數(shù)值模擬分析。常妍妍[8]在研究微成形工藝時引用一種新的數(shù)值模擬方法，將所用的多晶體材料看作由正六邊形組成的晶粒，每個晶粒代表一個特殊作用的有限單元，其每個晶粒單元由晶粒邊界單元連接，這些晶粒邊界單元用來描述晶粒之間的剪切變形。王廣春等[9]將模型中的每個晶粒單元視為由6個節(jié)點(diǎn)組成，在具體的網(wǎng)格劃分中將其分為2個四邊形單元。將該模型引入到研究激光微拉深和脹形的過程中，雖然模型有局限性，但為在連續(xù)介質(zhì)力學(xué)框架下有限元模擬微塑性成形提供了有效的建模思路。Zheng等[10]利用二維Voronoi多晶簡化模型數(shù)值研究了微成形過程中試樣晶粒分布對摩擦行為的影響，為后來建立基于Voronoi圖的三維多晶模型奠定了基礎(chǔ)。由于二維的多晶模型不能精確地反映晶粒取向分布不同對微塑性成形過程中摩擦因子和流動應(yīng)力的影響，因此建立三維的多晶幾何模型并進(jìn)行數(shù)值模擬，可以更清晰地反映變形過程中的尺寸效應(yīng)現(xiàn)象。Jiang等[11]基于Voronoi圖理論建立了一種適用于微尺度下塑性成形的材料多晶模型，在此模型中每個晶粒的材料流動應(yīng)力與晶粒的尺寸和成形溫度有關(guān)。此外，考慮晶界強(qiáng)化作用并簡化晶粒為立方體，提出了一種多晶幾何模型，從而實(shí)現(xiàn)了微尺度下鐓粗過程的數(shù)值模擬[12]。

文章研究尺寸為微尺度，不能忽視坯料中單個晶粒取向?qū)φ麄€變形的影響，因此基于Voronoi圖理論，建立考慮表面層理論的三維幾何模型，根據(jù)各個晶粒的取向不同，將區(qū)域分為3類，賦予其不同的本構(gòu)關(guān)系，得到三維材料本構(gòu)模型。利用有限元軟件ABAQUS實(shí)現(xiàn)微鐓粗過程的數(shù)值模擬，得到的流動應(yīng)力結(jié)果與實(shí)驗(yàn)數(shù)據(jù)相對比，以驗(yàn)證該模型的可行性。

1 微尺度下基于Voronoi圖的數(shù)值模型建立

1.1 多晶幾何模型的建立

在三維空間中，過去很多研究學(xué)者一般在多晶體模型中以六面體單元代表晶粒，但模擬結(jié)果不能準(zhǔn)確地反映晶粒的形狀不規(guī)則性和對材料變形過程的影響，因此并不是最好的模型選擇。研究發(fā)現(xiàn)，Voronoi圖從種子點(diǎn)生成六邊形的過程與真實(shí)多晶體材料的結(jié)晶過程具有一定的相似性[13-14]?；赩oronoi圖理論在有限元軟件ABAUQS中建立三維多晶幾何模型的具體過程如下:

(1)通過第三方軟件MATLAB，在Mpt工具箱中生成多晶體內(nèi)所有晶粒表示點(diǎn)的幾何信息。

(2)將多晶體中的每個晶粒所包含點(diǎn)的幾何信息與周圍其他晶粒相交，得到共同的幾何信息集合。若結(jié)果≥3個點(diǎn)，則得到了晶粒每個面所包含點(diǎn)的幾何信息；若結(jié)果為2個點(diǎn)，則得到了晶粒每個面所包含邊的信息。再將多晶體中晶粒所包含面和邊的全部幾何信息按照格式一并輸出為txt文本并保存。

(3)Python作為ABAQUS/CAE的腳本語言，可以實(shí)現(xiàn)模型的前處理過程。通過ABAQUS的腳本接口讀取 MATLAB中生成的數(shù)據(jù)，通過WirePloyLine命令自動生成晶粒每個面的邊，用AddFaces命令將晶粒的邊生成面。

(4)在Assembly模塊中將晶粒的面通過布爾運(yùn)算生成晶粒的殼，然后選中生成的殼和實(shí)體進(jìn)行移動到重合，再利用Assembly模塊中cut geometry功能，通過生成晶粒的殼和目標(biāo)實(shí)體可得到任意形狀的多晶幾何模型。

文章主要以圓柱體來建立材料多晶模型，如圖1所示。

圖1 基于Voronoi圖理論的三維多晶幾何模型圖

1.2 區(qū)域化本構(gòu)模型的建立

與單晶體相比，受位錯塞積影響的多晶體的晶粒在塑性變形過程中會在晶界附近產(chǎn)生很大的應(yīng)力場[15]。其中，多晶體中的晶粒遵循Hall-Petch公式，由式(1)表示為

式中:σpoly為多晶體變形中的流動應(yīng)力，MPa；M為多晶體晶粒的取向因子；τ為塑性變形中的總應(yīng)力，MPa；τi為沿滑移面滑動時產(chǎn)生的摩擦切應(yīng)力，MPa；ε為塑性變形中的應(yīng)變；為晶界附近由于位錯塞積對位錯運(yùn)動所產(chǎn)生的抵抗切應(yīng)力，MPa；d為平均晶粒尺寸，μm?？紤]塑性變形過程中的硬化作用，K′和τi都是應(yīng)變的函數(shù)。

由于傳統(tǒng)的多晶體包含晶粒數(shù)量非常大，單個晶粒對整個變形的影響可忽略，取其各個晶粒取向因子的平均值，對于面心立方晶格的多晶體，內(nèi)部晶粒取向因子的平均值為3.06。然而在微尺度下，變形區(qū)域只有幾個或十幾個晶粒，單個晶粒取向?qū)φ麄€變形的影響作用不能忽視。Huang[15]對多晶銅拉伸試樣的顯微組織進(jìn)行了分析。在純銅試樣的模擬中，每個三維模型均由上百個晶粒組成，可以根據(jù)二維的圖像反映三維空間晶粒的分布。根據(jù)在三維空間的不同取向分布，將晶粒大致分為3大類，如圖2所示。第1類晶粒主要分布在反極圖三角形的中間部分，以[411]為代表；第 2類晶粒大致分布在[100]周圍15°內(nèi)，并以[100]為代表；第3類晶粒分布在[111]周圍20°內(nèi)，作為[111]類型。課題組前期對純銅試樣晶粒位向分布進(jìn)行了測試[16]，電子背散射衍射技術(shù)測試的結(jié)果如圖3所示。

從知識產(chǎn)權(quán)評議在我國十余年的實(shí)踐來看，其所承載的功能不斷豐富和拓展：一方面，知識產(chǎn)權(quán)評議的目的從單純的防范重大經(jīng)濟(jì)科技活動的知識產(chǎn)權(quán)風(fēng)險(xiǎn)，逐漸延伸至優(yōu)化政府經(jīng)濟(jì)科技決策和政府管制；另一方面，知識產(chǎn)權(quán)評議的適用范圍從具體的經(jīng)濟(jì)科技活動（初期主要針對政府投資項(xiàng)目），逐漸擴(kuò)展至產(chǎn)業(yè)發(fā)展、科技進(jìn)步、國際投資和國際貿(mào)易等領(lǐng)域的政策制定、政策實(shí)施與宏觀層面的經(jīng)濟(jì)科技決策。

圖2 變形方向上純銅晶粒取向的分布圖[15]

圖3 前期測得純銅試樣電子背散射衍射反極圖

為了研究多晶與單晶純銅塑性變形行為的相關(guān)性，Huang等[15]分別將取向平行于[411]、[100]和[111]的晶粒作為面心立方多晶中出現(xiàn)的3種類型晶粒的代表，其取向因子M和體積分?jǐn)?shù)f見表1。

表1 3種不同類型晶粒的取向因子M及體積分?jǐn)?shù)f表[15]

在連續(xù)介質(zhì)力學(xué)理論的框架下，考慮到微塑性成形的變形機(jī)理，需要對建立的有限元多晶模型做以下假設(shè):

(1)表面層晶粒的位錯可以在外應(yīng)力作用下沿最有利滑移系統(tǒng)滑移，其取向因子M為2。

(2)在變形方向上，只有[411]和[100]及[111]3種具有代表性取向的晶粒屬于不存在自由表面的內(nèi)部晶粒，建模過程中3種晶粒按照1∶2∶7的比例進(jìn)行分配。

(3)在變形過程中不考慮晶粒轉(zhuǎn)動。

簡化流動應(yīng)力僅為應(yīng)變的冪指數(shù)，對文獻(xiàn)[1]給出的CuZn15單向壓縮實(shí)驗(yàn)數(shù)據(jù)進(jìn)行回歸分析，求解得到Hall-Petch公式中的待定系數(shù)，由式(2)表示為在變形方向上，不同取向晶粒的流動應(yīng)力應(yīng)變曲線如圖4所示。3種類型的代表晶粒的流動應(yīng)力都大于表面層晶粒，[111]方向上的晶粒的流動應(yīng)力最大?？梢杂帽砻鎸幽Ｐ屠碚摻忉?在試樣微小的情況下，表面層的晶粒相對較多，在試樣變形時，與內(nèi)部的各種類型晶粒比較，表面層晶粒受到的約束較小，在較低的流動應(yīng)力下就能發(fā)生變形，由此可得，微塑性成形中試樣的流動應(yīng)力隨尺寸的減小而減小。根據(jù)這些應(yīng)力—應(yīng)變數(shù)據(jù)可以定義在三維多晶模型的屬性上，使其進(jìn)行數(shù)值模擬分析。

圖4 不同取向晶粒的流動應(yīng)力—應(yīng)變曲線圖

2 微體積成型的模擬及驗(yàn)證

根據(jù)晶粒取向分布的體積分?jǐn)?shù)，隨機(jī)生成了3種不同晶粒取向分布的有限單元法模型(分別為模型1、模型2和模型3)。壓縮后3個模型剖面變形結(jié)果如圖6所示。在3種不同晶粒取向的多晶模型中，2側(cè)均出現(xiàn)了不同程度的凹凸不平，這與實(shí)驗(yàn)結(jié)果相吻合。其產(chǎn)生原因?yàn)?由于隨著試樣尺寸的減小，試樣中包含的晶粒個數(shù)也隨之減小，晶粒與晶粒之間由于晶粒取向的差異引起的力學(xué)性能的差異也愈加明顯，最終造成了試樣在變形過程中的不均勻性。

通過模擬結(jié)果及流動應(yīng)力—應(yīng)變關(guān)系可以得到3種不同晶粒取向隨機(jī)生成的多晶模型的應(yīng)力—應(yīng)變曲線，如圖7所示。在微鐓粗的過程中，由于尺寸發(fā)生微小的變化，材料產(chǎn)生了明顯的尺寸效應(yīng)。由于晶粒的分布和取向不同，最終由3種有限單元法模型得到的流動應(yīng)力也有所不同。相同應(yīng)變下，模型1的流動應(yīng)力較大，但三者相差不大，最大相差為4.33%。

圖5 微尺度下試樣單向壓縮的模擬過程圖

圖6 不同晶粒取向分布的有限單元法模型的剖面變形圖

圖7 由3種有限單元法模型得到的應(yīng)力—應(yīng)變關(guān)系圖

將3種FEM模型模擬得到的流動應(yīng)力的平均值與文獻(xiàn)[1]中的實(shí)驗(yàn)結(jié)果進(jìn)行對比，如圖8所示。模擬得到的流動應(yīng)力要略高于實(shí)驗(yàn)中得到的流動應(yīng)力，分析認(rèn)為是受模擬過程中選取的[111]取向的晶粒的影響，其加工硬化率要高于其他取向的晶粒，造成模擬結(jié)果中內(nèi)部多晶的流動應(yīng)力要高于實(shí)際情況。因此，所建立的三維有限元模型具有一定的可行性，能夠準(zhǔn)確地反映微鐓粗過程中的尺寸效應(yīng)現(xiàn)象。

圖8 模擬結(jié)果與實(shí)驗(yàn)結(jié)果的對比圖

3 結(jié)論

通過研究，得到以下結(jié)論:

(1)基于Voronoi圖理論，在ABAQUS中生成了圓柱體和立方體的三維多晶幾何模型，能夠更真實(shí)地反映試樣的晶粒形狀，修正了宏觀的Hall-Petch公式，得到了微尺度下的有限元多晶模型，證實(shí)了尺寸效應(yīng)現(xiàn)象及表面層晶粒的流動應(yīng)力小于內(nèi)部晶粒。

(2)通過對三維有限元模型進(jìn)行微鐓粗過程的數(shù)值模擬，晶粒分布不同對變形后的形貌產(chǎn)生不同的結(jié)果。數(shù)值模擬的數(shù)據(jù)與實(shí)驗(yàn)的結(jié)果相吻合，實(shí)現(xiàn)了實(shí)驗(yàn)與模擬的一致性，驗(yàn)證了三維多晶幾何模型和本構(gòu)模型的可行性。并且，此建模方法能更廣泛地分析微塑性成形中尺寸效應(yīng)現(xiàn)象。