夏海明，趙 晶，夏麗麗

(1.沈陽工業(yè)大學(xué)機械工程學(xué)院，遼寧 沈陽 110870) (2.特變電工沈陽變壓器股份有限公司，遼寧 沈陽 110144)

電力變壓器各法蘭盤聯(lián)接處常常會產(chǎn)生振幅極小的微動，當(dāng)微動作用在電力變壓器各聯(lián)接處的橡膠密封件表面時，會導(dǎo)致其接觸界面磨損直至密封失效。本文通過有限元軟件ABAQUS進(jìn)行仿真計算，根據(jù)不同壓縮率和位移幅值計算得到的橡膠密封件應(yīng)力和應(yīng)變結(jié)果，并對橡膠密封件在變壓器實際工況中的微動損傷行為[1-2]進(jìn)行探討。

1 有限元仿真過程

圖1所示為變壓器法蘭盤實際工況圖，可以看出電力變壓器橡膠密封件兩側(cè)始終承受著法蘭盤施加的徑向載荷，在電力變壓器橡膠密封件凸起的接觸表面發(fā)生微動磨損。為了有效地研究變壓器橡膠密封件的實際損傷狀況，同時使仿真結(jié)果呈現(xiàn)的變化規(guī)律與實際工況更加吻合，本文對橡膠密封件的夾具進(jìn)行設(shè)計，同時在密封件上放置磨塊，在橡膠密封件與磨塊之間的接觸界面設(shè)置振幅極小的微動，建立仿真模型的優(yōu)點是在仿真分析過程中可以準(zhǔn)確地施加壓縮量和位移幅值這兩個運動參數(shù)，變壓器橡膠密封件分析模型如圖2所示。

圖1 變壓器法蘭盤實際工況圖

圖2 變壓器橡膠密封件分析模型

1.1 橡膠材料參數(shù)的確定

橡膠材料是超彈性材料，通常基于應(yīng)變勢能函數(shù)來表征超彈性材料的應(yīng)力與應(yīng)變關(guān)系，例如Arruda-Boyce、Marlow、Mooney-Rivlin以及Neo-hooke等應(yīng)變勢能函數(shù)。根據(jù)大量實驗取證，Mooney-Rivlin函數(shù)[3]可以準(zhǔn)確地貼合橡膠特性，故本文采用該函數(shù)來設(shè)置仿真過程的材料參數(shù)。Mooney-Rivlin函數(shù)具體公式如下:

(1)

式中：i,j為多項式階次；W為應(yīng)變能密度；Cij為材料常數(shù)；I1和I2分別為第1和第2應(yīng)變不變量。材料常數(shù)C01和C10通過經(jīng)驗公式計算:

(2)

C01=0.25C10

(3)

E=6(C10+0.25C01)

(4)

(5)

式中：E為彈性模量；HA為材料常數(shù)。由于橡膠材料具有不可壓縮性，因此在ABAQUS軟件中建立材料模型時，壓縮比設(shè)置為0，選取的橡膠密封件硬度為70，計算出彈性模量后，根據(jù)公式(2)～(5)計算得Mooney-Rivlin函數(shù)的兩項材料常數(shù)C10和C01分別為0.738 9和0.184 7，夾具和磨塊設(shè)置為剛體，楊氏模量為206 000 Pa，泊松比為0.3。

1.2 雙鼓型橡膠壓縮量計算

根據(jù)目前國內(nèi)對于大型電力變壓器橡膠密封領(lǐng)域相關(guān)研究現(xiàn)狀，橡膠材料的壓縮率通常情況下應(yīng)控制在10%～20%，常見的計算壓縮率ε的公式為：

(6)

式中：d為壓縮前橡膠密封件的初始長度，mm；l為壓縮后橡膠密封件的長度，mm。如圖3所示，在橡膠密封件兩側(cè)由夾具施加夾緊力并將橡膠材料壓縮率控制在10%～20%。為了更好地表現(xiàn)出橡膠材料壓縮率對于橡膠材料損傷行為的影響，分別設(shè)置壓縮率為10%、12%、14%以及16%，通過式(6)計算得出壓縮量分別為3.0,3.6,4.2以及4.8 mm。

2 設(shè)定運動條件并創(chuàng)建網(wǎng)格劃分

圖3所示為施加運動條件示意圖，為了使不銹鋼磨塊作穩(wěn)定的往復(fù)運動，將往復(fù)運動設(shè)置為正弦型往復(fù)運動，材料接觸界面之間的往復(fù)運動頻率設(shè)定為100 Hz,位移幅值分別設(shè)定為0.2,0.5,1.0以及2.0 mm，通過橡膠專用夾具對橡膠密封件施加沿Y軸方向由外向里的位移，例如為了實現(xiàn)總壓縮量達(dá)到3.0 mm，使兩端夾具分別壓縮1.5 mm，即可達(dá)到10%的壓縮率，其余橡膠密封件壓縮率的設(shè)置以此類推。

圖3 施加運動條件示意圖

模型網(wǎng)格采用二次減縮積分六面體單元C3D20R進(jìn)行劃分，而不規(guī)則形狀的夾具部分則采用二次四面體單元進(jìn)行劃分，劃分后通過命令Verify Mesh檢查網(wǎng)格，最終結(jié)果如圖4所示。

圖4 仿真模型劃分網(wǎng)格示意圖

3 有限元模型計算結(jié)果分析

3.1 不同橡膠材料壓縮率有限元分析結(jié)果

圖5所示為不同壓縮率條件下橡膠密封件應(yīng)力和應(yīng)變的變化情況，可以發(fā)現(xiàn)雙鼓型橡膠密封件的最大應(yīng)力值以及最大應(yīng)變值均在正弦型往復(fù)運動的最大極值位置。壓縮率越高，雙鼓型橡膠密封件的最大應(yīng)力值以及最大應(yīng)變值越大。

圖5 不同材料壓縮率條件下應(yīng)力、應(yīng)變變化曲線

圖6所示為不同壓縮率條件下橡膠密封件應(yīng)力變化情況，可以觀察到在不同壓縮率條件下，應(yīng)力集中即橡膠材料易出現(xiàn)疲勞失效的位置大致相同，均處在橡膠凸起表面與橡膠平面交界處；同時根據(jù)仿真情況可以看出，雙鼓型橡膠密封件底部的平面部分在運行過程中也比較容易出現(xiàn)疲勞損傷現(xiàn)象。根據(jù)仿真結(jié)果可知，橡膠密封件的變形往往都是從外側(cè)表面先開始，且變形較為劇烈，而橡膠密封件中心部位則依靠橡膠的彈塑性變形抵消材料變形，所以在橡膠密封件中心部位并未發(fā)現(xiàn)較明顯的變形。

圖6 不同材料壓縮率條件下的應(yīng)力圖

當(dāng)橡膠密封件壓縮率分別為12%、14%和16%時，雙鼓型橡膠密封件的最大應(yīng)力值即凸起部位與平面交界處的應(yīng)力值分別為2.018,2.411和2.744 MPa。從仿真結(jié)果可以發(fā)現(xiàn)，隨著橡膠壓縮率增大，在往復(fù)運動過程中雙鼓型橡膠密封件的應(yīng)力值也逐漸增大，表明壓縮率增大會導(dǎo)致橡膠密封件的疲勞損傷更加嚴(yán)重[4]。

圖7所示為不同壓縮率條件下橡膠密封件應(yīng)變變化情況。當(dāng)橡膠材料壓縮率分別為12%、14%和16%時，雙鼓型橡膠密封件的最大應(yīng)變值即凸起部位與平面交界處的應(yīng)變值分別為29.61%、34.60%和38.60%。從應(yīng)變結(jié)果中可以發(fā)現(xiàn)，隨著雙鼓型橡膠密封件壓縮率增大，在高頻微動下橡膠的變形也越來越大，因此壓縮率成為影響變壓器橡膠密封件微動磨損的主要因素。

圖7 不同材料壓縮率條件下的應(yīng)變圖

3.2 不同位移幅值有限元分析結(jié)果

圖8所示為不同位移幅值條件下橡膠密封件應(yīng)力及應(yīng)變的變化情況，隨著位移幅值增大，雙鼓型橡膠密封件的最大應(yīng)力值和應(yīng)變值同時增大，表明隨著位移幅值的增大會加速雙鼓型橡膠密封件的磨損[5-6]。

圖8 不同位移幅值條件下應(yīng)力、應(yīng)變變化曲線

圖9所示為不同位移幅值條件下橡膠密封件應(yīng)力變化情況，當(dāng)施加位移幅值u分別達(dá)到0.2，0.5，1.0和2.0 mm時，雙鼓型橡膠密封件的凸起部位與平面交界處的應(yīng)力值分別為1.797，1.841，2.030和2.336 MPa。

圖9 不同位移幅值條件下的應(yīng)力圖

圖10所示為不同位移幅值條件下橡膠密封件應(yīng)變變化情況。當(dāng)施加位移幅值u分別達(dá)到0.2，0.5，1.0和2.0 mm時，橡膠密封件凸起部位與平面交界處的應(yīng)變值分別為26.67%、27.27%、29.78%和33.67%。在同一個周期內(nèi)，當(dāng)運動到往復(fù)運動極值點位置時，應(yīng)力值和應(yīng)變值達(dá)到了最大值，雙鼓型橡膠易產(chǎn)生疲勞損傷，而且隨著位移幅值的增大，損傷情況也更加嚴(yán)重。

圖10 不同位移幅值條件下的應(yīng)變圖

4 結(jié)束語

在變壓器實際工況中，位移幅值和壓縮率都會成為影響橡膠密封件密封能力的重要因素[7-8]。橡膠密封件壓縮量越大，橡膠密封件的微動損傷越嚴(yán)重；如果加載較大的橡膠材料壓縮率，在一定程度上會加速橡膠密封件的變形，從而導(dǎo)致橡膠密封件接觸表面發(fā)生疲勞損傷甚至萌生裂紋；在高頻極小位移幅值的條件下，隨著位移幅值的增加會加劇橡膠密封件損傷，但如果位移幅值達(dá)到一定值時，繼續(xù)增大位移幅值則材料會達(dá)到穩(wěn)定損傷狀態(tài)。本文運用ABAQUS仿真軟件來模擬橡膠密封件實際工況，為處在各類實際工況中的橡膠類密封件仿真計算提供了一些計算思路和問題解決方法。