華磊 劉雪峰 常冬梅

摘要：將陰陽極反應(yīng)分離的元胞自動機法引入到鋁合金多坑腐蝕過程的模擬當(dāng)中。通過腐蝕模擬得到了較為符合實際特征的鋁合金材料多坑腐蝕形貌圖。通過計算不同腐蝕步數(shù)下的蝕坑形貌參數(shù)，發(fā)現(xiàn)蝕坑深度近似于隨腐蝕時間線性增加，同時腐蝕失重也隨時間增加，但是增加的速度逐漸增加。將腐蝕形貌導(dǎo)入ANSYS軟件中進(jìn)行力學(xué)分析，給出了鋁合金材料腐蝕后的應(yīng)力分布，并且由此計算了腐蝕后的剩余疲勞壽命。發(fā)現(xiàn)隨著腐蝕時間的增加，剩余疲勞壽命大體上呈逐漸減少的趨勢。但是，由于蝕坑底端附近局部幾何形貌造成應(yīng)力集中的影響，這一趨勢并不一定總是單調(diào)的。

關(guān)鍵詞：鋁合金；腐蝕；元胞自動機；疲勞

中圖分類號：TG171文獻(xiàn)標(biāo)識碼：ADOI：10.19452/j.issn1007-5453.2021.09.009

基金項目：航空科學(xué)基金（201902067004）；中央高?；究蒲袠I(yè)務(wù)費項目（3122018C008）

鋁合金材料具有重量（質(zhì)量）輕、比強度大等優(yōu)點，因此在航空航天領(lǐng)域得到廣泛的應(yīng)用。但是鋁合金材料在服役過程中容易發(fā)生腐蝕，從而影響材料的可靠性。因此，對鋁合金腐蝕問題的研究受到廣泛的關(guān)注。過去對于鋁合金腐蝕問題的研究以試驗方法為主，但是試驗方法往往費時費力，成本較高，因此近些年在腐蝕數(shù)值模擬方面的研究越來越多。多數(shù)腐蝕數(shù)值模擬方法是針對一系列有關(guān)材料和環(huán)境參數(shù)方程的求解，因此為確定性方法[1]。確定性方法可以較準(zhǔn)確地模擬出材料的宏觀腐蝕進(jìn)程，但是，確定性方法沒有考慮到材料和環(huán)境在介觀或微觀尺度上的不均勻性，因此不適合模擬由此造成的材料較小尺度上的不規(guī)則腐蝕。而這種較小尺度的不規(guī)則腐蝕往往對材料的疲勞、斷裂等力學(xué)性質(zhì)有較大的影響。因此，在研究鋁合金的不規(guī)則腐蝕以及其造成的應(yīng)力集中等力學(xué)問題時，有必要引入非確定性腐蝕模擬方法。

元胞自動機方法是一種典型的非確定性算法，最早被用來模擬生物的生長問題，之后被應(yīng)用于植被生長過程[2]、交通流[3]和火災(zāi)蔓延過程[4]等很多方面的模擬。近年來，該方法作為一種有效的非確定性方法被引入到金屬材料的腐蝕模擬中，取得了較好的效果。2004年，Vautrin等[5]在元胞自動機模型中允許陰極和陽極反應(yīng)在不同位置同時發(fā)生反應(yīng)，很好地模擬了真實的電化學(xué)腐蝕過程，并且由此驗證了蝕坑內(nèi)的氫離子濃度指數(shù)（pH）不均衡現(xiàn)象。Aarao等[6]在2006年對這一模型的相關(guān)參數(shù)進(jìn)行了分析。之后出現(xiàn)了一系列基于這一模型的研究[7-9]。2016年，Pérez-Brokate等[10]將該模型引入到三維情況下，得到了更符合實際的模擬結(jié)果。2013年，Wang[11]和Han[12]將有限元與元胞自動機腐蝕模型結(jié)合，研究了應(yīng)力作用下的腐蝕模擬問題。隨后將研究引申到了雙蝕坑的問題上。但是，到目前為止，基于元胞自動機模擬得到的腐蝕形貌對材料力學(xué)性質(zhì)的研究仍然很少見到，而在這方面的研究有利于分析介觀到宏觀尺度腐蝕對材料力學(xué)性質(zhì)的影響。因此，本文將基于陰陽極分離的元胞自動機腐蝕模型模擬鋁合金多坑腐蝕形貌，并且在此基礎(chǔ)上研究鋁合金疲勞壽命隨腐蝕進(jìn)程的變化。

1基于元胞自動機的多坑腐蝕形貌

1.1元胞自動機模型

在酸性堿性元胞的擴散方面，假設(shè)A、B元胞在每個反應(yīng)步中會發(fā)生Ndiff次布朗運動，當(dāng)其運動到E元胞的位置，則與之交換位置，而當(dāng)A元胞運動到B元胞的位置或相反時，A、B元胞中和，兩個元胞均變?yōu)镋元胞。每個反應(yīng)步中，每個C元胞都有一定概率發(fā)生失效，該元胞和其下面的相鄰M元胞變?yōu)镋元胞。

以上反應(yīng)的概率見表1。其中，NA-B為鄰居中A、B元胞個數(shù)差，NE為鄰居元胞中液態(tài)元胞（A、B、E）的個數(shù)。在每一步腐蝕模擬過程中，每次隨機選取一對可以發(fā)生反應(yīng)的R元胞和S元胞，按照式（6）～式（9）發(fā)生反應(yīng)，反復(fù)選取，直到?jīng)]有可供反應(yīng)的元胞為止。然后進(jìn)行Ndiff步A、B元胞的擴散模擬，每次隨機選取一個A或B元胞進(jìn)行擴散，直到所有A、B元胞選取結(jié)束算作一個擴散步。Ndiff步擴散步結(jié)束后，算作一步腐蝕模擬結(jié)束，進(jìn)行下一步腐蝕模擬，直到結(jié)束所有模擬步。

1.2模擬結(jié)果

本文采用摩爾鄰居，元胞個數(shù)為1024×128個，設(shè)PB為1/100000，PSSE為0.2，PRE為0，PPE為1/8，PPE為1/800，即為PPE的1/100。因此，在中性條件下，再鈍化產(chǎn)物的水解概率為酸性條件下的1/100。每個反應(yīng)步中的擴散步數(shù)Ndiff為200。在MATLAB中編寫程序，由此得到500、1000和1500反應(yīng)步后的腐蝕形貌如圖1（a）～圖1（c）所示，腐蝕形貌與實際腐蝕形貌特征相符。由圖可以明顯地看到，隨著計算步數(shù)的增加，蝕坑逐漸增大，而且有新的蝕坑逐漸產(chǎn)生。蝕坑邊界具有明顯的不規(guī)則性，接近于真實的蝕坑，這與確定性方法得到的結(jié)果是不同的。在第500步時，橫坐標(biāo)0～100區(qū)域內(nèi)有三個小蝕坑。由圖1（b）可見，到第1000步時，這三個小蝕坑已經(jīng)合并成為一個大蝕坑，而且相比于其他蝕坑，這個合并的大蝕坑更不規(guī)則。在試驗中經(jīng)常會看到類似的蝕坑合并現(xiàn)象。

針對模型的隨機性，每個步數(shù)進(jìn)行5次計算取平均值，得到平均腐蝕失重率隨著步數(shù)的變化曲線如圖2所示。其中，腐蝕失重為A、B、E元胞個數(shù)除以總元胞個數(shù)。腐蝕失重隨時間逐漸增加，而且隨著時間的增加，腐蝕失重增加的速度有所增加。

除了材料的腐蝕失重率以外，蝕坑的深度也是受關(guān)注的焦點。尤其是所有蝕坑中的最大深度，對材料的損傷情況有較大的影響。本文中，將R和P元胞作為蝕坑表面，由此將蝕坑表面的最低點與上表面的距離（元胞個數(shù)）作為最大蝕坑深度。與前面相同，計算不同腐蝕步數(shù)情況下5次模擬的最大蝕坑平均深度如圖3所示。可以看到，最大蝕坑平均深度近似隨著腐蝕步數(shù)線性增加。

2疲勞壽命分析

在本文中，利用有限元軟件ANSYS進(jìn)行鋁合金腐蝕后的疲勞壽命分析。首先需要將腐蝕形貌圖導(dǎo)入到ANSYS中。而ANSYS無法自動識別腐蝕形貌，需要用其他軟件將蝕坑邊界信息提取出來再導(dǎo)入到ANSYS中建立模型。

首先，識別蝕坑的邊界。如圖1可以看到，蝕坑的邊界實際上是由點描繪的圖形。這里采用軟件R2V提取材料的邊界，然后將結(jié)果導(dǎo)入AutoCAD，以圖1（b）中的1000步的腐蝕形貌圖為例，得到邊界圖如圖4所示。

在AutoCAD中調(diào)整邊界尺寸，假設(shè)每個元胞的邊長為2μm，則整個模型的尺寸為2048μm×256μm。將圖形保存后導(dǎo)入ANSYS進(jìn)行受力分析。選取的鋁合金材料楊氏模量為7.17×104MPa，泊松比0.3，左側(cè)固定，右側(cè)施加10N/μm的拉應(yīng)力，進(jìn)行計算得到x方向的應(yīng)力云圖如圖5所示?？梢钥吹?，在蝕坑底端附近出現(xiàn)了明顯的應(yīng)力集中現(xiàn)象，通常在較大的蝕坑底端的應(yīng)力集中更明顯。

應(yīng)力集中會大大影響材料的疲勞壽命。以計算得到的應(yīng)力為基礎(chǔ)，可以進(jìn)一步計算鋁合金材料腐蝕后的剩余疲勞壽命?；贏NSYS疲勞分析模塊，以模型的最大Mises應(yīng)力作為疲勞分析的最大應(yīng)力，以0作為最小應(yīng)力計算疲勞壽命，S—N曲線數(shù)值表見表2。其中參數(shù)由試驗數(shù)據(jù)進(jìn)行擬合得到，由于本文中主要針對疲勞壽命的趨勢問題進(jìn)行定性研究，因此具體的參數(shù)選擇僅取為典型值即可。

在本例中，由此得到其對應(yīng)的疲勞壽命為2.43×105次。

按照以上步驟，重新進(jìn)行一次腐蝕進(jìn)程模擬，分別將步數(shù)為500、600、800、1000、1200、1400、1500時的腐蝕形貌導(dǎo)入ANSYS中，得到材料疲勞壽命隨腐蝕步數(shù)的變化曲線如圖6所示。隨著腐蝕步數(shù)的增加，疲勞壽命大體上是逐漸減少的。但是，這一趨勢并不是線性的，甚至當(dāng)腐蝕步數(shù)為600步時的疲勞壽命還高于500步時。這是因為在這種條件下的疲勞壽命變化主要與應(yīng)力集中有關(guān)，而應(yīng)力集中除取決于蝕坑深度外，也與局部的腐蝕形貌有很大關(guān)系。在有些時候盡管蝕坑深度已經(jīng)較大，但是蝕坑底端的表面形貌比較光滑，因此局部應(yīng)力反而有所減小。例如，在本例中，腐蝕步數(shù)為500和600步時，最大應(yīng)力分別為32.14MPa和30.03MPa，由于在500步時應(yīng)力集中造成的最大應(yīng)力更大，因此在500步時計算得到的疲勞壽命反而更少。但是由于蝕坑深度是單調(diào)增加的，因此疲勞壽命的整體趨勢隨腐蝕步數(shù)的增加而逐漸減少。

3結(jié)論

本文基于元胞自動機法模擬了鋁合金材料的多坑電化學(xué)腐蝕過程，在算法中允許陰陽極反應(yīng)發(fā)生在不同位置。通過模擬，成功得到了較符合實際特征的腐蝕形貌圖。在此基礎(chǔ)上，將腐蝕形貌進(jìn)行處理后導(dǎo)入ANSYS軟件中進(jìn)行力學(xué)分析，得到了鋁合金材料腐蝕后的應(yīng)力分布，并由此分析了鋁合金材料腐蝕后的疲勞壽命，得到了以下結(jié)論：

（1）運用陰陽極反應(yīng)分離的元胞自動機算法可以成功模擬鋁合金的腐蝕進(jìn)程，并且可以觀測到實際應(yīng)用中常見到的蝕坑合并現(xiàn)象。

（2）在本算法模擬的條件下，鋁合金的腐蝕失重隨時間逐漸增加，而且增速逐漸加大，蝕坑的最大深度接近于隨時間線性增加。

（3）鋁合金腐蝕后的剩余疲勞壽命大體上隨腐蝕時間減少，但是并不完全單調(diào)遞減。這是因為疲勞壽命與應(yīng)力集中現(xiàn)象有關(guān)，而應(yīng)力集中除取決于蝕坑深度以外，也受到蝕坑底部局部幾何形貌的影響，而這隨著腐蝕過程而劇烈變化。

參考文獻(xiàn)

[1]Walton J C，Cragnolino G，Kalandros S K. A numerical model of crevice corrosion for passive and active metals[J]. Corrosion Science，1996，38（1）：1-18.

[2]Dunn A G，Majer J D. Simulating weed propagation via hierarchical，patch-based cellular automata[C]//Computational Science-ICCS 2007：7th International Conference Proceedings，2007：762-769.

[3]彭麟，譚惠麗，孔令江，等.開放性邊界條件下雙車道元胞自動機交通流模型耦合效應(yīng)研究[J].物理學(xué)報，2003， 52（12）： 3007-3013. Peng Lin， Tan Huili， Kong Lingjiang， et al. A study of coupling effect in cellular automata model of traffic flow for two-lane with open boundary conditions[J]. Acta Physica Sinica，2003， 52（12）： 3007-3013.（in Chinese）

[4]Encinas A H，Encinas L H，White S H，et al. Simulation of forest fire fronts using cellular automata[J]. Advances in Engineering Software，2007，38：372-378.

[5]Vautrin-Ul C，Chausse A，Stafiej J，et al. Simulations of corrosion processes with spontaneous separation of cathodic and anodic reaction zones[J]. Polish Journal of Chemical Technology，2004，78（9）：1795-1810.

[6]Aarao F D A，Stafiej J，Badiali J P. Scaling theory in a model of corrosion and passivation[J]. The Journal of Physical Chemistry B，2006，110（35）：17554-17562.

[7]Vautrin-Ul C，Taleb A，Stafiej J，et al. Mesoscopic modelling of corrosion phenomena：coupling between electro chemical and mechanical processes，analysis of the deviation from the Faraday law[J]. Electrochimica Acta，2007，52（17）：5368-5376.

[8]Di Caprio D，Stafiej J. Cellular automata approach to corrosion and passivity phenomena[J]. Pure and Applied Chemistry，2010，85（1）：62-67.

[9]Bartosik？，Di Caprio D，Stafiej J. Cellular automata approach to corrosion and passivity phenomena[J]. Pure and Applied Chemistry，2013，85（1）：247-256.

[10]Pérez-Brokate C F，Di Caprio D，F(xiàn)éron D. Three dimensional discrete stochastic model of occluded corrosion cell[J]. Corrosion Science，2016，111：230-241.

[11]Wang H，Han E. Simulation of metastable corrosion pit development under mechanical stress[J]. Corrosion Science，2013，90：128-134.

[12]Wang H，Han E. Computational simulation of corrosion pit interactions under mechanochemical effects using a cellular automaton/finite element model[J]. Corrosion Science，2016，103：305-311.

[13]Liu X F，Hua L，Chang D M. Simulation of pit interactions of multi-pit corrosion under an anticorrosive coating with a threedimensional cellular automata model[J]. Modelling and Simula tion in Materials Science and Engineering，2021，29（6）：065018.

Study on Multi-Pits Corrosion of Aluminum Alloy and Its Fatigue Life Based on Cellular Automata

Hua Lei1，Liu Xuefeng1，Chang Dongmei2

1. Civil Aviation University of China，Tianjin 300300，China

2. Tianjin Key Laboratory of High Speed Cutting and Precision Machining，Tianjin University of Technology and Education，Tianjin 300222，China

Abstract： In this paper， the cellular automata method of reaction separation between anode and cathode is introduced into the simulation of multi-pit corrosion process of aluminum alloy. Through corrosion simulation， the multipit corrosion morphology of aluminum alloy was obtained. It is found that the pit depth increases linearly with the corrosion time approximately， and the corrosion weight loss also increases with the corrosion time， but the increasing speed increases with time gradually. The corrosion morphology is imported into ANSYS software for mechanical analysis， the stress distribution of aluminum alloy after corrosion is obtained， and the residual fatigue life after corrosion is calculated. It is found that the residual fatigue life decreases with corrosion time. However， this trend is not always monotonous due to the effect of local geometry near the bottom of the pit on the stress concentration.

Key Words： aluminum alloys； corrosion； cellular automata； fatigue