柳彬，游世輝，趙樹勛，曾憲任

?

基于復(fù)雜網(wǎng)絡(luò)的磁流變橡膠磁致壓縮力學(xué)性能研究

柳彬，游世輝，趙樹勛，曾憲任

(湘潭大學(xué) 土木工程與力學(xué)學(xué)院，湖南 湘潭 411105)

基于磁流變橡膠在機(jī)車的主動(dòng)減振及降噪領(lǐng)域中的重要作用，引入復(fù)雜網(wǎng)絡(luò)理論和方法來研究磁流變橡膠的磁致壓縮力學(xué)性能。通過建立磁流變橡膠的RVE模型，利用多物理場耦合有限元軟件研究磁流變橡膠的靜態(tài)磁致壓縮力學(xué)性能；計(jì)算磁流變橡膠的磁性顆粒復(fù)雜網(wǎng)絡(luò)特征參數(shù)；探索磁流變橡膠磁性顆粒復(fù)雜網(wǎng)絡(luò)的特征參數(shù)與磁流變橡膠磁致壓縮模量之間的關(guān)系；與實(shí)例研究結(jié)果進(jìn)行對比，驗(yàn)證結(jié)論的準(zhǔn)確性。研究結(jié)果表明：磁流變橡膠的磁致壓縮模量與其磁性顆粒網(wǎng)絡(luò)的特征參數(shù)有關(guān)，其磁致壓縮模量的峰值是由磁性顆粒的復(fù)雜網(wǎng)絡(luò)分形特征參數(shù)決定的。

磁流變橡膠；復(fù)雜網(wǎng)絡(luò)；磁致壓縮模量；數(shù)值模擬；實(shí)驗(yàn)研究

磁流變橡膠(MR rubber)是由鐵磁顆粒嵌入橡膠基體中制成的磁流變材料。相較于傳統(tǒng)的橡膠材料，磁流變橡膠兼具了磁性顆粒和橡膠材料的優(yōu)點(diǎn)，并且其力學(xué)和磁學(xué)性能可隨外加磁場的強(qiáng)度變化而改變，因此廣泛地應(yīng)用于高速列車、鐵路斜拉索橋、汽車和艦船等的主動(dòng)減振及降噪等工程領(lǐng)域方面[1?2]。丁志華等[3]研究設(shè)計(jì)的磁流變橡膠空氣彈簧可實(shí)時(shí)調(diào)節(jié)高度和剛度的問題，提高各類機(jī)車的減振性能，該設(shè)計(jì)已獲得國家發(fā)明專利[4]。段富貴等[5?6]研究的磁流變橡膠的吸聲與隔聲性能可應(yīng)用于吸聲板和隔聲板。近年來，磁流變橡膠的理論實(shí)驗(yàn)研究[7?9]已經(jīng)相當(dāng)成熟，在數(shù)值模擬方面，Lokander等[10]研究磁場強(qiáng)度和顆粒對磁流變橡膠的磁流變效應(yīng)的影響，以及磁場對磁流變橡膠阻尼的影響。段富貴等[5?6]分別研究磁流變橡膠的吸聲和隔聲性能。Watts等[11?12]發(fā)表了研究復(fù)雜網(wǎng)絡(luò)的特性的方法，掀起了復(fù)雜網(wǎng)絡(luò)的研究熱潮。Silva等[13]把網(wǎng)絡(luò)理論運(yùn)用到復(fù)合材料電導(dǎo)系數(shù)與網(wǎng)絡(luò)障礙之間的關(guān)系中，論證了復(fù)雜網(wǎng)絡(luò)可以應(yīng)用于材料科學(xué)。Roberto等[14]研究了顆粒狀材料接觸網(wǎng)絡(luò)的分形維數(shù)與力閾值的函數(shù)關(guān)系。段富貴等[5?6]將復(fù)雜網(wǎng)絡(luò)應(yīng)用到多孔橡膠材料的吸聲與隔聲性能與復(fù)雜網(wǎng)絡(luò)參數(shù)的關(guān)聯(lián)性的研究中。隨著計(jì)算機(jī)技術(shù)的進(jìn)步，現(xiàn)階段可通過核磁共振成像技術(shù)或CT掃描技術(shù)獲得磁流變橡膠細(xì)觀結(jié)構(gòu)。在細(xì)觀尺度上，磁流變橡膠中的磁性顆粒形成磁彈耦合結(jié)構(gòu)網(wǎng)絡(luò)，并與基體的分子鏈網(wǎng)絡(luò)相互作用，在整體上形成互相關(guān)聯(lián)的結(jié)構(gòu)網(wǎng)絡(luò)系統(tǒng)。當(dāng)在結(jié)構(gòu)網(wǎng)絡(luò)系統(tǒng)上施加耦合場時(shí)，系統(tǒng)會(huì)產(chǎn)生細(xì)觀和宏觀的磁彈性耦合和粘彈性效應(yīng)，以及細(xì)觀和宏觀之間的相互作用的結(jié)構(gòu)網(wǎng)絡(luò)系統(tǒng)的演化[15]。在此基礎(chǔ)上，本文提出引入復(fù)雜網(wǎng)絡(luò)理論來描述細(xì)觀尺度上的磁性顆粒結(jié)構(gòu)網(wǎng)絡(luò)系統(tǒng)，進(jìn)而來描述磁流變橡膠的磁性顆粒網(wǎng)絡(luò)與磁致力學(xué)性能之間的關(guān)聯(lián)性。

1 復(fù)雜網(wǎng)絡(luò)理論基礎(chǔ)

復(fù)雜網(wǎng)絡(luò)的研究基礎(chǔ)是圖論，即它是一個(gè)有著復(fù)雜拓?fù)浣Y(jié)構(gòu)特征的圖。復(fù)雜網(wǎng)絡(luò)特征的主要度量參數(shù)包括度分布，集聚系數(shù)、平均路徑長度和分形維數(shù)。網(wǎng)絡(luò)的度分布反映了整個(gè)網(wǎng)絡(luò)拓?fù)涞闹匾畔ⅲ侵妇W(wǎng)絡(luò)中任一個(gè)節(jié)點(diǎn)所具有的與其他節(jié)點(diǎn)相連接的邊的數(shù)目，節(jié)點(diǎn)度分布用函數(shù)()來表示。網(wǎng)絡(luò)的集聚系數(shù)[13]表示節(jié)點(diǎn)與k個(gè)相鄰節(jié)點(diǎn)連接的緊密程度，可用公式表達(dá)為

將2個(gè)節(jié)點(diǎn)連接起來的最短邊數(shù)即為節(jié)點(diǎn)的距離，那么網(wǎng)絡(luò)的平均路徑長度則為任意2個(gè)節(jié)點(diǎn)之間的距離的平均值。分形維數(shù)是對復(fù)雜網(wǎng)絡(luò)分形定量表征的重要參數(shù)，本文的磁性顆粒復(fù)雜網(wǎng)絡(luò)分形維數(shù)采用LIU[16]提出的改進(jìn)沙箱算法計(jì)算。

現(xiàn)階段，能夠描述接近真實(shí)世界特性的復(fù)雜網(wǎng)絡(luò)模型主要為小世界網(wǎng)絡(luò)和無標(biāo)度網(wǎng)絡(luò)模型[17]。在本文中磁性顆粒間作用并沒有方向性，故選取構(gòu)造無向無權(quán)網(wǎng)絡(luò)模型。利用大型復(fù)雜網(wǎng)絡(luò)分析軟件Pajek對每個(gè)磁性顆粒隨機(jī)編號，如圖1所示，將每個(gè)磁性顆粒等效成一個(gè)節(jié)點(diǎn)，根據(jù)3-正則圖的對偶圖和鄰近原則構(gòu)造磁性顆粒復(fù)雜網(wǎng)絡(luò)模型如圖2所示。本文中復(fù)雜網(wǎng)絡(luò)平均度、簇聚系數(shù)和平均最短路徑的值，都是采用Pajek軟件計(jì)算的。

為了便于計(jì)算機(jī)處理，網(wǎng)絡(luò)圖也可用矩陣表示。鄰接矩陣是表示頂點(diǎn)之間相鄰關(guān)系的矩陣，鄰接矩陣可分為有向和無向鄰接矩陣。對于本文中的無向網(wǎng)絡(luò)而言，鄰接矩陣一定是對稱的，且主對角線一定為0，則圖2的鄰接矩陣的表示如式(2)。將鄰接矩陣導(dǎo)入到MATLAB中根據(jù)改進(jìn)的沙箱算法進(jìn)一步計(jì)算出復(fù)雜網(wǎng)絡(luò)的分形維數(shù)。

圖1 磁性顆粒分布示意圖

圖2 磁性顆粒復(fù)雜網(wǎng)絡(luò)模型

2 磁流變橡膠的磁致壓縮力學(xué)性能研究

對于耦合場作用下的磁流變橡膠的數(shù)值模擬，需利用大型多物理場耦合軟件COMSOL構(gòu)建磁流變橡膠RVE模型，并施加周期性邊界條件，設(shè)定耦合場的方向和強(qiáng)度如圖3所示。該模型尺寸為20 μm×20 μm，顆粒直徑介于3~5 μm。在磁性顆粒的邊界用Maxwell應(yīng)力張量施加磁場力，圖4顯示了磁場條件下磁流變橡膠RVE模型在壓縮條件下的應(yīng)力變形。

圖3 磁流變橡膠壓縮RVE模型

圖4 各向同性磁流變橡膠壓縮變形

圖5顯示了壓縮條件下的磁性顆粒組分對磁流變橡膠的磁致力學(xué)性能的影響。磁流變橡膠的磁致壓縮應(yīng)力隨體積分?jǐn)?shù)的增加而增加。而磁性顆粒體積分?jǐn)?shù)20%以上的磁致壓縮應(yīng)力值變化較明顯。圖6所示的磁流變橡膠磁致壓縮模量與應(yīng)變曲線關(guān)系是由將圖5的曲線的斜率轉(zhuǎn)換成的，隨著磁性顆粒體積分?jǐn)?shù)的增加，磁流變橡膠的磁致壓縮模量逐漸增大。當(dāng)磁性顆粒體積分?jǐn)?shù)增至25%時(shí)，磁性顆粒組分對磁流變橡膠的磁致壓縮模量的影響較為 顯著。

圖5 磁流變橡膠的磁致壓縮應(yīng)力-應(yīng)變曲線

圖6 磁流變橡膠磁致壓縮模量-應(yīng)變曲線

磁致壓縮模量的影響因素眾多，本文除了研究磁性顆粒組分對磁流變橡膠的磁致壓縮模量的影響外，還將探討磁場強(qiáng)度對壓縮條件下磁流變橡膠磁致模量的影響。

圖7顯示了磁場強(qiáng)度對磁流變橡膠磁致壓縮力學(xué)性能的影響。結(jié)果表明磁流變橡膠的磁致壓縮應(yīng)力隨外部磁場強(qiáng)度增大而增加。而磁場強(qiáng)度對磁流變橡膠壓縮模量的影響如圖8所示。在小變形條件下，磁流變橡膠磁性顆粒體積分?jǐn)?shù)小于10%時(shí)，外部磁場對磁流變橡膠磁致壓縮模量幾乎沒有影響。當(dāng)磁流變橡膠磁性顆粒體積比大于15%時(shí)，磁場強(qiáng)度對磁流變橡膠磁致壓縮模量的影響逐步增大。以上數(shù)值模擬結(jié)果已在文獻(xiàn)[18]與實(shí)驗(yàn)數(shù)據(jù)對比 研究。

圖7 磁流變橡膠磁致壓縮應(yīng)力-應(yīng)變曲線

圖8 磁流變橡膠的磁致壓縮模量-顆粒體積分?jǐn)?shù)曲線(%)

3 磁流變橡膠磁性顆粒復(fù)雜網(wǎng)絡(luò)特性

在復(fù)雜網(wǎng)絡(luò)系統(tǒng)上施加耦合場時(shí)，網(wǎng)絡(luò)系統(tǒng)會(huì)產(chǎn)生細(xì)觀和宏觀的磁彈性耦合效應(yīng)，而復(fù)雜網(wǎng)絡(luò)就承擔(dān)起構(gòu)建宏觀力學(xué)性能與細(xì)觀拓?fù)浣Y(jié)構(gòu)特性的橋梁。本節(jié)利用復(fù)雜網(wǎng)絡(luò)的特征參數(shù)來構(gòu)建細(xì)觀網(wǎng)絡(luò)結(jié)構(gòu)特征與宏觀磁流變橡膠的磁性顆粒體積分?jǐn)?shù)的基本規(guī)律。

研究結(jié)果如圖9~12所示。圖9顯示了各向同性磁流變橡膠磁性顆粒網(wǎng)絡(luò)的度值分布情況，從圖中可以看出磁性顆粒網(wǎng)絡(luò)的平均度值隨磁性顆粒體積分?jǐn)?shù)的增大而增大，最后趨近于平穩(wěn)。而磁性顆粒網(wǎng)絡(luò)集聚系數(shù)的分布情況圖10所示，磁性顆粒網(wǎng)絡(luò)的集聚系數(shù)隨體積分?jǐn)?shù)的增大而減小，最終趨近于平穩(wěn)。圖11中的磁性顆粒網(wǎng)絡(luò)平均最短路徑值則隨體積分?jǐn)?shù)的增大而增大，曲線接近線性關(guān)系。磁流變橡膠磁性顆粒網(wǎng)絡(luò)的分形維數(shù)隨磁性顆粒體積分?jǐn)?shù)的增長先迅速增加，而后在某個(gè)區(qū)域上下波動(dòng)，有接近平穩(wěn)的趨勢，如圖12所示。

圖9 磁流變橡膠磁性顆粒網(wǎng)絡(luò)的平均度-體積分?jǐn)?shù)(100%)

圖10 磁流變橡膠顆粒網(wǎng)絡(luò)的集聚系數(shù)-體積分?jǐn)?shù)(100%)

圖11 磁流變橡膠顆粒網(wǎng)絡(luò)的平均最短路徑-體積分?jǐn)?shù)(100%)

圖12 磁性顆粒體積分?jǐn)?shù)(%)-磁性顆粒網(wǎng)絡(luò)分形維數(shù)曲線

4 磁流變橡膠的磁致力學(xué)性能與磁性顆粒網(wǎng)絡(luò)特征分析

本文分別研究磁流變橡膠的磁致壓縮模量、磁流變橡膠的磁性顆粒復(fù)雜網(wǎng)絡(luò)參數(shù)與磁性顆粒組分關(guān)系，本節(jié)則研究磁流變橡膠的磁致力學(xué)性能與磁性顆粒復(fù)雜網(wǎng)絡(luò)參數(shù)的關(guān)聯(lián)性。

分別將圖9~12與圖8進(jìn)行轉(zhuǎn)換，得到壓縮條件下磁流變橡膠的磁致力學(xué)參數(shù)與磁性顆粒復(fù)雜網(wǎng)絡(luò)參數(shù)的關(guān)系曲線，結(jié)果如圖13~16所示。圖13中磁流變橡膠的磁致壓縮模量值隨磁性顆粒網(wǎng)絡(luò)平均度的增加而增大；當(dāng)磁性顆粒網(wǎng)絡(luò)的平均度大于3時(shí)，磁致壓縮模量的變化開始明顯。圖14中磁流變橡膠磁致壓縮模量隨磁性顆粒網(wǎng)絡(luò)平均集聚系數(shù)的增加而逐漸減??；當(dāng)平均集聚系數(shù)大于0.8時(shí)，磁致壓縮模量幾乎不隨平均集聚系數(shù)而變化。圖15顯示了磁流變橡膠磁性顆粒網(wǎng)絡(luò)的平均最短路徑長度與磁致壓縮模量之間的關(guān)系，總體上，磁流變橡膠的壓縮模量值隨平均最短路徑長度的增加而增大；但當(dāng)平均最短路徑長度小于1.5時(shí)，磁流變彈性體的磁致壓縮模量幾乎不受外部磁場的影響。在圖16中，磁流變橡膠的磁致壓縮模量隨磁性顆粒網(wǎng)絡(luò)的分形維數(shù)增加而先增大后減?。淮帕髯兿鹉z磁致壓縮模量的峰值處的分形維數(shù)對應(yīng)的磁性顆粒體積分?jǐn)?shù)介于25%~30%之間，與現(xiàn)有實(shí)驗(yàn)研究[9]結(jié)果吻合。

圖13 磁性顆粒網(wǎng)絡(luò)的平均度-磁致壓縮模量曲線

圖14 磁性顆粒網(wǎng)絡(luò)的平均集聚系數(shù)-磁致縮模量曲線

圖15 磁性顆粒網(wǎng)絡(luò)的平均路徑長度-磁致壓縮模量曲線

圖16 磁性顆粒網(wǎng)絡(luò)分形維數(shù)-磁致壓縮模量曲線

5 磁流變橡膠的實(shí)驗(yàn)研究

實(shí)驗(yàn)制備磁流變橡膠所用的球形羥基鐵粉顆粒的平均尺寸為3~5 μm，電鏡照片見圖17。在制備磁流變橡膠時(shí)，先干燥羥基鐵粉，然后逐步將羰基鐵粉與硅油混合確保鐵磁顆粒分布均勻且不粘連，再加入到橡膠中攪拌均勻。將均質(zhì)的鐵粉橡膠混合物倒入平面模具中開模，脫氣。再在液壓機(jī)中固化。以同樣的方法制備球形羥基鐵粉體積比為5%~30%的磁流變橡膠樣品。用掃描電子顯微鏡SEM和光學(xué)顯微鏡檢查實(shí)驗(yàn)制備樣品的SEM表面如圖18所示。

圖17 羰基鐵粉的SEM形貌

圖18 各向同性磁流變橡膠樣本斷裂表面SEM形貌

準(zhǔn)備羥基鐵粉體積比為30%的磁流變橡膠材料樣本進(jìn)行壓縮測試，用于壓縮測試的磁流變橡膠樣本如圖19。磁流變橡膠的靜態(tài)壓縮測試結(jié)果以及實(shí)驗(yàn)與數(shù)值模擬結(jié)果對比可參考文獻(xiàn)[18]，本章直接研究實(shí)驗(yàn)制備磁流變橡膠的復(fù)雜網(wǎng)絡(luò)特性。在計(jì)算實(shí)驗(yàn)制備磁流變橡膠的復(fù)雜特性之前，需將樣本SEM圖像通過圖形圖像轉(zhuǎn)換技術(shù)[19]轉(zhuǎn)換成可直接導(dǎo)入有限元軟件中計(jì)算的矢量圖，在轉(zhuǎn)換中，排除相交的磁性顆粒，只選取最上層的磁性顆粒，得到磁流變橡膠矢量圖如圖20所示。任意選取20 μm× 20 μm細(xì)觀局部結(jié)構(gòu)圖如圖21。

在計(jì)算實(shí)驗(yàn)制備的磁流變橡膠樣本的復(fù)雜網(wǎng)絡(luò)特性時(shí)，需要借助Pajek構(gòu)建接近真實(shí)的復(fù)雜網(wǎng)絡(luò)模型，并對磁性顆粒的復(fù)雜網(wǎng)絡(luò)參數(shù)進(jìn)行處理。將實(shí)驗(yàn)制備的磁流變橡膠的磁性顆粒復(fù)雜網(wǎng)絡(luò)參數(shù)與數(shù)值模擬結(jié)果對比可得到圖22~25。對于復(fù)雜網(wǎng)絡(luò)的4個(gè)參數(shù)：平均度、平均集聚系數(shù)、平均最短路徑和分形維數(shù)，雖然數(shù)值模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)有細(xì)微的差別，但各參數(shù)與體積分?jǐn)?shù)曲線的發(fā)展趨勢是一致的，驗(yàn)證了數(shù)值模擬結(jié)果的正確性。

圖19 磁敏橡膠壓縮測試樣本

圖20 各向同性磁流變橡膠樣本二維矢量圖

圖21 磁流變橡膠樣本局部細(xì)觀圖

圖22 平均度與顆粒體積分?jǐn)?shù)(%)的模擬與實(shí)驗(yàn)數(shù)據(jù)對比

圖23 平均集聚系數(shù)與顆粒體積分?jǐn)?shù)(%)的模擬與實(shí)驗(yàn)數(shù)據(jù)對比

圖24 平均最短路徑與顆粒體積分?jǐn)?shù)(%)的模擬與實(shí)驗(yàn)數(shù)據(jù)對比

圖25 分形維數(shù)與顆粒體積分?jǐn)?shù)(%)的模擬與實(shí)驗(yàn)數(shù)據(jù)對比

6 結(jié)論

1) 磁流變橡膠的磁致壓縮模量可以通過磁性顆粒復(fù)雜網(wǎng)絡(luò)特征參數(shù)來定性描述。通過磁流變橡膠的細(xì)觀結(jié)構(gòu)特征，估算磁流變橡膠的宏觀力學(xué)性能。

2) 磁流變橡膠的磁致壓縮模量的峰值只與磁性顆粒網(wǎng)絡(luò)的分形維數(shù)有關(guān)，而與磁場強(qiáng)度大小無關(guān)。即可通過磁流變橡膠的復(fù)雜網(wǎng)絡(luò)分形維數(shù)確定磁性顆粒的最佳體積比。本文中的磁性顆粒最佳體積比介于25%~30%之間，與已有的研究結(jié)果相 吻合。

[1] 劉少謙. 磁流變橡膠減振器的設(shè)計(jì)與研究[D]. 廈門: 華僑大學(xué), 2008. LIU Shaoqian. Design and research of magnetorheological rubber shock absorber[D]. Xiamen: Huaqiao University, 2008.

[2] 孫航, 高秀云, 陳大鵬. 基于碳納米管加強(qiáng)的磁流變橡膠支座在連續(xù)梁橋減隔震中的應(yīng)用分析[EB/OL]. http:// www.paper.edu.cn/releasepaper/content/201401[2014?01 ?26]. SUN Hang, GAO Xiuyun, CHEN Dapeng. Application analysis of magnetorheological rubber support based on carbon nanotubes in seismic reduction of continuous beam bridges[EB/OL]. http://www.paper.edu.cn/release epaper/content/201401[2014?01?26].

[3] 丁志華, 游世輝. 磁流變彈性體橡膠空氣彈簧的結(jié)構(gòu)設(shè)計(jì)[J]. 機(jī)械設(shè)計(jì)與制造, 2008(10): 17?18. DING Zhihua, YOU Shihui. Structural design of magnetorheological elastomeric rubber air spring[J]. Mechanical Design & Manufacture, 2008(10): 17?18.

[4] 游世輝, 肖靜, 馮云華, 等. 一種磁流變液壓橡膠彈簧: 中國, 201010508586.4[P]. 2014?12?10. YOU Shihui, XIAO Jing, FENG Yunhua, et al. A magnetorheological hydraulic rubber spring: China, 201010508586.4[P]. 2014?12?10.

[5] 段富貴. 磁流變橡膠材料隔聲性能數(shù)值模擬[D]. 湘潭:湘潭大學(xué), 2017. DUAN Fugui. Numerical simulation of sound insulation performance of magnetorheological rubber material[D]. Xiangtan: Xiangtan University, 2017.

[6] 楊罡星. 磁流變橡膠材料吸聲性能數(shù)值模擬[D]. 湘潭:湘潭大學(xué), 2017. YANG Gangxing. Numerical simulation of sound absorption properties of magnetorheological rubber materials[D]. Xiangtan: Xiangtan University, 2017.

[7] 孫友誼. 磁流變橡膠材料的研究現(xiàn)狀[R]. 國家自然基金委員會(huì)、揚(yáng)州大學(xué)、北京市新型高分子材料制備加工成型重點(diǎn)實(shí)驗(yàn)室、中國化工學(xué)會(huì)橡膠專業(yè)委員會(huì), 2016: 1. SUN Youyi. Research status of magnetorheological rubber materials[R]. National Natural Science Foundation of China, Yangzhou University, Beijing Key Laboratory of New Polymer Material Preparation, Processing and Molding, China Rubber Industry Committee, 2016: 1.

[8] 游世輝. 磁敏橡膠力學(xué)性能的研究與應(yīng)用[D]. 廣州: 華南理工大學(xué), 2008. YOU Shihui. Research and application of mechanical properties of magnetic resin[D]. Guangzhou: South China University of Technology, 2008.

[9] Ioan Bicaa, Eugen M Anitas, Madalin Bunoiuc, et al. Hybrid magnetorheological elastomer: Influence of magnetic field and compression pressure on its electrical conductivity[J]. Journal of Industrial and Engineering Chemistry, 2014, 20(6): 3994?3999.

[10] Lokander M, Stenberg B. Improving the magnet- orheological effect in isotropic magnetorheological rubber materials[J]. Polymer Testing, 2003, 22(6): 677? 680.

[11] Watts D J, Strogatz S H. Collective dynamics of small world networks[J]. Nature, 1998(393): 440?442.

[12] Barabási A L, Albert R. Emergence of scaling in random networks[J]. Science, 1999(286): 509?512.

[13] Silva J, Simoes R, Lanceros-Mendez S, et al. Applying complex network theory to the understanding of high-aspect-ratio carbon-filled composites[J]. Euro Physics Letters, 2011, 93(3): 37005.

[14] Roberto Arévalo, Iker Zuriguel, Diego Maza. Topological properties of the contact network of granular materials[J]. International Journal of Bifurcation & Chaos, 2009, 19(2): 695?702.

[15] 張波. 復(fù)雜網(wǎng)絡(luò)的構(gòu)建及演化方式研究[D]. 長春: 吉林大學(xué), 2014. ZHANG Bo. Research on the construction and evolution of complex networks[D]. Changchun: Jilin University, 2014.

[16] LIU J L, YU Z G, Anh V. Determination of multifractal dimensions of complex networks by means of the sandbox algorithm[J]. Chaos, 2015, 25(2): 392?395.

[17] 孫璽菁, 司守奎. 復(fù)雜網(wǎng)絡(luò)算法與應(yīng)用[M]. 北京: 國防工業(yè)出版社, 2015. SUN Xijing, SI Shoukui. Complex network algorithm and application[M]. Beijing: Defense Industry Press, 2015.

[18] 柳彬, 游世輝, 趙樹勛, 等. 基于Digimat中RVE模型的磁流變彈性體的磁致壓縮力學(xué)性能的數(shù)值模擬[J].功能材料, 2018, 49(2): 67?75. LIU Bin, YOU Shihui, ZHAO Shuxun, et al. Simulation of magnetic field induced compressive mechanical properties of magneto-rheological elastomers based on RVE model in Digimat[J]. Functional Materials, 2018, 49(2): 67?75.

[19] 宋歡. 基于圖形圖像轉(zhuǎn)換的磁流變彈性體拉伸力學(xué)性能的數(shù)值模擬[D]. 湘潭: 湘潭大學(xué), 2017. SONG Huan. Numerical simulation of tensile mechanical properties of magnetorheological elastomer based on graphical image transformation[D]. Xiangtan: Xiangtan University, 2017.

Research on magneto-rheological rubber compression mechanical properties based on complex network

LIU Bin, YOU Shihui, ZHAO Shuxun, ZENG Xianren

(College of Civil Engineering and Mechanics, Xiangtan University, Xiangtan 411105, China)

Magneto-rheological rubber plays an important role in the field of active vibration and noise reduction of locomotives. This paper aims to introduce complex network theory and methods to study the magneto-induced compression mechanical properties of magneto-rheological rubber. By establishing the RVE model of magneto-rheological rubber, the static magneto-induced compressive mechanical properties of magneto- rheological rubber were simulated using COMSOL multiphysics software. And the topological structure and fractal characteristics of the complex network of magnetic particles were studied. Finally, the association of the magneto-induced compression mechanical properties and complex network parameters was explored. And the results of the simulation study were compared with the experimental study data to verify the accuracy of the conclusions. The results show that the magneto-induced modulus of magneto-rheological rubber is related to the characteristic parameters of the magnetic particles network, and the peak value of magneto-induced modulus is determined by the fractal characteristic parameters.

magneto-rheological rubber; complex network; magneto-induced compressive modulus; numerical simulation; experimental study

10.19713/j.cnki.43?1423/u.2019.04.027

TB381

A

1672 ? 7029(2019)04 ? 1041 ? 09

2018?05?27

國家自然科學(xué)基金資助項(xiàng)目(51375416)；湖南省教育廳創(chuàng)新平臺開放基金資助項(xiàng)目(15K131)

游世輝(1962?)，男，江西萍鄉(xiāng)人，教授，博士，從事工程中的力學(xué)問題與計(jì)算研究；E?mail：shihuiyou@hotmail.com

(編輯 蔣學(xué)東)