HNBR與DBP共混物玻璃化溫度及機械性能分子模擬*

王世杰，孫小涵，楊 斌，唐美玲，孫靜婷

(沈陽工業(yè)大學(xué) 機械工程學(xué)院，沈陽 110870)

定子橡膠在一定程度上影響著潛油螺桿泵采油系統(tǒng)中的使用壽命，而定子橡膠在其加工和使用過程中因受到溫度和化學(xué)介質(zhì)等作用，會逐漸出現(xiàn)粘連、變硬、龜裂等由于物理性能失衡導(dǎo)致的一系列現(xiàn)象.因此，在潛油螺桿泵行業(yè)中使用的橡膠制品應(yīng)具有良好的耐高溫、耐低溫、耐油與耐臭氧老化性能.氫化丁腈橡膠(HNBR)是一種分子主鏈高度飽和的橡膠材料，主要通過丁腈橡膠(NBR)氫化制得[1]，即將NBR鏈段上的丁二烯單元進行選擇性氫化，將C—C不飽和雙鍵加氫生成單鍵，既保持了原有NBR的耐油、耐介質(zhì)和耐磨性，又兼具優(yōu)良的耐熱、耐候、耐氧、耐化學(xué)介質(zhì)腐蝕性、高拉伸強度和長期壓縮永久變形較小的綜合性能，特別是其熱老化性能得以極大提升.HNBR分子鏈中的丙烯腈單元提供了材料優(yōu)異的耐油性和高拉伸強度規(guī)格穩(wěn)定性，可用于生產(chǎn)重復(fù)性優(yōu)異的聚合物.然而，高玻璃化轉(zhuǎn)變溫度(Tg)和不良的機械性能阻礙了其應(yīng)用.為了改善HNBR的Tg和機械性能，常用的方法是添加適量的合適增塑劑.普通增塑劑鄰苯二甲酸二丁酯(DBP)[2]可以顯著降低HNBR的分子間相互作用，從而可以降低Tg并改善機械性能.因此，有必要研究DBP對HNBR性能的影響.物理上穩(wěn)定的分散體需要考慮由溶度參數(shù)(δ)判斷的增塑劑與粘合劑之間的互溶性.聚合物的性能可以通過關(guān)鍵指標(biāo)Tg進行評估，該指標(biāo)的值決定工藝溫度范圍.根據(jù)Fox-Flory的自由體積理論，大多使用自由體積(VF)和密度(ρ)方法來計算Tg，即溫度對應(yīng)于密度隨溫度變化曲線的調(diào)整點.此外，機械性能對于聚合物的制備、加工和應(yīng)用也很重要.

許多研究人員通過MD方法對聚合物的混溶性、玻璃化轉(zhuǎn)變溫度和力學(xué)性能進行了模擬，所得結(jié)果與實驗結(jié)果具有很好的一致性，證實了模擬的可靠性[3].分子模擬可用來計算聚合物的物理和機械性能[4].

1 分子建模過程

構(gòu)建HNBR單鏈[1]，由丙烯腈、順式1，4-丁二烯、反式1，4-丁二烯、1，2-丁二烯和丁烷單體組成的重復(fù)單元結(jié)構(gòu)并將其首尾連接，各基體組分如表1所示.

表1 單鏈各基體結(jié)構(gòu)及配比

單鏈重復(fù)數(shù)暫定50，構(gòu)建后的HNBR單鏈構(gòu)成及DBP分子模型分別如圖1～2所示.對HNBR單鏈和DBP分子分別在Forcite模塊中利用Smart Minimizer功能對能量進行最小化，當(dāng)能量達(dá)到最低時，用Amorphous Cell模塊分別建立了純HNBR，以及DBP體積分?jǐn)?shù)為0.1，0.15，0.2，0.25的HNBR/DBP二元體系的混合物模型.繪制出85HNBR/15DBP的AC模型，如圖3所示.

圖1 HNBR單鏈

圖2 DBP分子模型

圖3 85HNBR/15DBP的AC模型

通過分子力學(xué)軟件Material Studio7.0的Forcite模塊對混合物模型進行優(yōu)化.對于DBP較小的分子結(jié)構(gòu)，在273 K下進行優(yōu)化可以得到最優(yōu)的幾何構(gòu)型.然而，對于HNBR復(fù)雜的高分子結(jié)構(gòu)，在273 K下優(yōu)化得不到最優(yōu)的幾何構(gòu)型，因此，將模擬溫度設(shè)置為298 K.選用NVT系綜、COMPASS力場[5]對高分子體系進行MD模擬，時間為500 ps，在0.1 MPa下進行退火處理，將溫度從100 K升溫至遠(yuǎn)高于HNBR玻璃化溫度的500 K再降溫至100 K，經(jīng)過5次循環(huán)，體系在高溫下經(jīng)歷多次松弛，得到能量最低的最優(yōu)結(jié)構(gòu)[6].

2 溶度參數(shù)的確定

溶度參數(shù)(δ)的定義為內(nèi)聚能密度(CED)的平方根，已廣泛用于研究材料的物理混溶性，它是指每單位體積蒸發(fā)一摩爾液體所需的能量，其數(shù)學(xué)表達(dá)式[7]為

(1)

(2)

式中：EV為汽化過程中內(nèi)部能量的變化；HV為蒸發(fā)的焓；Vm為摩爾數(shù)在蒸發(fā)溫度下液體的體積；R和T分別為氣體常數(shù)和絕對溫度；φ1、φ2分別為組分1、2的體積分?jǐn)?shù).由式(1)、(2)可得，δ1、δ2的值越接近，則ΔHV越??；ΔHV越小，體系相容性越好.因此，溶度參數(shù)差對各組分間的相容性預(yù)測有一定參考價值.

基本溶解度理論表明，分子溶質(zhì)的溶解度與CED有關(guān).可以通過MD模擬計算出的δ來預(yù)測混合材料的相溶性.對HNBR和DBP的平衡AC結(jié)構(gòu)在NPT系綜下進行500 ps的MD模擬，然后在NVT系綜下進行再一次500 ps的MD模擬來計算溶度參數(shù).在模擬過程中，均采用周期性邊界條件.溫度控制選擇Andersen恒溫器，壓力控制則選擇Berendsen恒壓器.選用Verlet速度時間積分方法(步長為1 fs)，用于積分牛頓運動方程.范德華相互作用選取Atom Based，截斷半徑為19.5 ?.靜電相互作用選擇Ewald方法，精度為0.000 1 kcal/mol[8].模擬結(jié)果為不同聚合度下HNBR與DBP的溶度參數(shù)，如圖4所示.表2為HNBR與DBP溶度參數(shù)范圍.

圖4 HNBR與DBP溶度參數(shù)

表2 HNBR與DBP溶度參數(shù)范圍

對于高分子體系，分子間如果不存在強極性基團，若溶度參數(shù)之差的絕對值小于1.3～2.1(J·cm-3)1/2，則混合體系相容[9].顯然，本文中二者的溶度參數(shù)絕對值遠(yuǎn)小于1.3(J·cm-3)1/2，故二者相容.由圖4可知，HNBR的溶度參數(shù)在聚合度達(dá)到30之后保持在20.4(J·cm-3)1/2，即HNBR單體鏈長設(shè)置為50合理.

3 體系相容性預(yù)測

徑向分布函數(shù)g(r)是反映材料微觀結(jié)構(gòu)的特征物理量，它是在距離原子一定距離處找到另一個原子的概率與在完全均勻分布的情況下找到一個原子的概率之間的比[10].g(r)反映出原子排布的有序性問題.如果近程出現(xiàn)一個峰值，峰高且陡，表明有序性較高，屬于晶體體系，原子緊密相連.如果出現(xiàn)小于0.3 nm的峰，則表示分子鏈有序性較低，屬于非晶結(jié)構(gòu)即無定型結(jié)構(gòu).分子間徑向分布函數(shù)表示原子自由散落之間的相互作用方式，氫鍵作用范圍為0.26～0.31 nm，范德華力作用范圍為0.31～0.5 nm[11].

HNBR、DBP物質(zhì)的分子間和分子內(nèi)的徑向分布函數(shù)如圖5～6所示.由圖5可知，HNBR分子在0.3 nm直到平穩(wěn)過程中沒有出現(xiàn)強峰，說明結(jié)構(gòu)為無定型結(jié)構(gòu)，分子間主要作用方式為范德華力.由圖6可知，HNBR與DBP二者分子內(nèi)作用力的徑向分布函數(shù)曲線分布較為接近，依據(jù)相似相容原理，可以證明二者相容性較優(yōu).

圖5 HNBR與DBP分子間作用力

圖6 HNBR與DBP共混物的分子內(nèi)作用力

4 玻璃化轉(zhuǎn)變溫度

玻璃化轉(zhuǎn)變溫度(Tg)是聚合物的重要物理性能，決定了其應(yīng)用溫度范圍.獲取Tg的方法是在NPT系綜下設(shè)置平衡混合系統(tǒng)，然后使用腳本完成計算.本質(zhì)上，玻璃化轉(zhuǎn)變是聚合物材料中的常見現(xiàn)象.當(dāng)聚合物經(jīng)歷玻璃化轉(zhuǎn)變時，其各種物理性質(zhì)都會發(fā)生劇烈變化.目前，還沒有關(guān)于玻璃化轉(zhuǎn)變的完整理論來解釋其實驗事實.在玻璃化轉(zhuǎn)變中，摩爾體積根據(jù)溫度變化而變化.Fox和Flory提出了自由體積理論來解釋玻璃化轉(zhuǎn)變，在Tg以下，隨著溫度的升高，聚合物的自由體積或空體積的變化很小，但卻在玻璃化狀態(tài)時自由體積出現(xiàn)突變[12]；共聚物的密度(ρ)隨溫度的拐點也是求得Tg的一種方法，本文根據(jù)不同物質(zhì)含量的二元混合體系的摩爾體積隨溫度變化的圖像得出其玻璃化轉(zhuǎn)變溫度.在NPT系綜下將模擬溫度設(shè)定在一定溫度范圍內(nèi)(200～600 K).用最小二乘法分別調(diào)整玻璃化轉(zhuǎn)變溫度以上和以下的溫度范圍的數(shù)據(jù)點，得到兩條直線，兩條線的交點為玻璃化轉(zhuǎn)變溫度Tg.

在NPT下進行了500 ps的MD模擬，采用平衡體積分?jǐn)?shù)為0、0.10、0.15、0.20和0.25的HNBR/DBP二元體系的AC結(jié)構(gòu).隨后進行高溫退火處理，使用NPT系綜在大氣壓下以20 K的間隔將溫度從600 K降至200 K.選擇Andersen恒溫器和Berendsen恒壓器分別控制溫度和壓力.通過Ewald方法計算靜電相互作用，并且通過基于原子的方法計算范德華相互作用.重復(fù)該步驟，直到溫度遠(yuǎn)低于HNBR/DBP的Tg.從平均偏摩爾體積隨溫度的變化曲線中觀察到的斜率轉(zhuǎn)折點被認(rèn)為是HNBR/DBP的Tg，計算結(jié)果如圖7所示.

圖7 HNBR與DBP共混物的玻璃化轉(zhuǎn)變溫度

在Tg以下，密度隨溫度升高而降低，因為分子鏈的能量不足以克服在一定溫度范圍內(nèi)的旋轉(zhuǎn)能壘.一旦溫度超過一定范圍，由于系統(tǒng)處于高能量狀態(tài)，平均偏摩爾體積的提升速度會加快.此外，由于DBP降低了HNBR之間的相互作用，因此，Tg隨著DBP體積分?jǐn)?shù)的增加而降低.

5 力學(xué)性能

力學(xué)性能主要是由于聚合物系統(tǒng)在不同的彈性性能測試模式下表現(xiàn)出的性能.定義如下特征模量：楊氏模量(E)、體積模量(K)和剪切模量(G).材料的應(yīng)力σ和應(yīng)變張量ε之間的最一般關(guān)系遵循廣義的胡克定律，其表達(dá)式為

σi=Cijεj

(3)

式中，C為彈性系數(shù)的對稱(6×6)矩陣.矩陣的對稱性將系數(shù)減小到21，以描述各向異性材料的應(yīng)力和應(yīng)變之間的關(guān)系，并且當(dāng)材料中存在彈性對稱表面時，數(shù)量將繼續(xù)減少.基于羅伊斯均值方法，可以通過下式計算有效的體積模量和剪切模量，即

KR=[S11+S11+S11+2(S12+S23+S31)]-1

(4)

GR=5[4(S11+S22+S33)-4(S12+S23+S31)+

(S44+S55+S66)]

(5)

式中：下標(biāo)R表示羅伊斯平均值；柔度矩陣S為彈性系數(shù)矩陣C的逆矩陣，即S=C-1，可以通過MD模擬來計算.將DBP體積分?jǐn)?shù)為0、0.10、0.15、0.20、0.25的平衡HNBR/DBP的AC模型在NPT系綜進行進一步的動力學(xué)計算，持續(xù)時間為500 ps，然后通過NVT系綜的MD模擬執(zhí)行500 ps的平衡結(jié)構(gòu)計算HNBR/DBP的K和G.獲得K和G后，計算拉伸模量(E)和泊松比(v)，即

E=2G(1+v)=3K(1-2v)

(6)

機械性能是聚合物最重要的性能之一，主要由聚合物增塑劑決定.拉伸模量(E)與拉伸彈性有關(guān)；體積模量(K)測量聚合物的不可壓縮性；剪切模量(G)是指材料的剛度，即E、K的值表明HNBR/DBP在受到外力作用時易變形的程度.泊松比(v)反映了材料的橫向變形.此外，K/G和柯西壓力是描述材料延展性的另外兩個主要參數(shù).它們與材料承受應(yīng)力時抵抗裂紋形成的能力有關(guān).K/G表示塑性范圍(拉伸伸長率)，因此K/G的高值與延展性相關(guān)，而K/G的低值與脆性相關(guān).K/G值越高，延展性越好.對于柯西壓力，負(fù)值表示脆性，而正值表示延展性.為了獲得上述參數(shù)，在298 K下進行了HNBR的力學(xué)性能模量計算，結(jié)果如表3所示.為了研究DBP對HNBR力學(xué)性能的影響，在298 K下對DBP體積分?jǐn)?shù)為0、0.10、0.15、0.20、和0.25的HNBR/DBP體系進行了力學(xué)性能模量分析，如表3所示.

表3 不同條件下HNBR的有效各向同性力學(xué)性能

表4列出了298 K下HNBR/DBP系統(tǒng)的彈性系數(shù)和機械性能.從C12～C44的正值可以得出在298 K下的HNBR/DBP共混物是易延展材料的結(jié)論.此外，HNBR/DBP體系的泊松比均約為0.285，表明HNBR/DBP體系具有塑性體的性質(zhì)(0.2

表4 不同條件下HNBR的彈性系數(shù)

圖8 HNBR機械模量

6 結(jié) 論

本文通過分子動力學(xué)(MD)模擬預(yù)測HNBR/DBP二元體系的混溶性、玻璃化轉(zhuǎn)變溫度和機械性能，得出如下結(jié)論：通過MD模擬計算了HNBR和DBP的徑向分布函數(shù)，以預(yù)測其相容性，表明HNBR和DBP是可混溶的；根據(jù)平均偏摩爾體積-溫度圖得出0、0.10、0.15、0.20和0.25的HNBR/DBP玻璃化轉(zhuǎn)變溫度變化曲線，結(jié)果表明，Tg隨DBP體積分?jǐn)?shù)的增加而降低；分析了彈性常數(shù)，例如體積模量、剪切模量、拉伸模量、柯西壓力、K/G及泊松比v，結(jié)果表明，適當(dāng)添加DBP可以提高HNBR的延展性并降低其脆性，可以促進HNBR的應(yīng)用.