李亞莎, 胡豁然, 夏 宇, 李光竹, 劉清東

(三峽大學(xué) 電氣與新能源學(xué)院，宜昌 443002)

1 引 言

聚乙烯電纜因其具有優(yōu)異的導(dǎo)電性能、機(jī)械性能以及耐環(huán)境應(yīng)力而備受關(guān)注，自問世以來就廣泛應(yīng)用于高壓輸電線路和配電網(wǎng)中.然而隨著現(xiàn)代電力工業(yè)的急速發(fā)展，日益增加的用電需求迫使電力電纜的傳動性和安全性需要有高度保障，聚乙烯材料易變形、不耐熱的缺點(diǎn)也隨之放大[1].長時間的熱老化作用以及水分子的滲透會使聚乙烯內(nèi)部產(chǎn)生缺陷或結(jié)構(gòu)性破壞，降低聚乙烯材料的絕緣性能[2-4].

為了解決這一系列問題，人們提出聚合物無機(jī)納米復(fù)合材料，借助納米顆粒的某些優(yōu)良特性來改善聚合物的性能[5].納米顆粒本身具有獨(dú)特的小尺寸效應(yīng)、表面界面效應(yīng)以及量子效應(yīng)而使得其不僅能提高復(fù)合材料的耐直流擊穿能力，同時也能提高復(fù)合材料的力學(xué)性能等特性[6-8].趙等人[9]將硅烷偶聯(lián)處理后的MgO與LDPE共混制得MgO/LDPE納米復(fù)合材料，通過測試發(fā)現(xiàn)質(zhì)量分?jǐn)?shù)為4%的納米MgO填料能有效抑制空間電荷的注入，且復(fù)合介質(zhì)可以抑制電樹枝的引發(fā)和生長.王有元等采用MgO、ZnO、SiO2三種納米氧化物顆粒為填料對LDPE進(jìn)行改性[10].結(jié)果表明，摻雜MgO使得LDPE具有最佳的抗氧化性和熱氧老化防護(hù)效果.然而因?yàn)榧{米顆粒尺寸過于微小，且絕緣材料的諸多優(yōu)異性能由納米復(fù)合材料的界面結(jié)構(gòu)[11]決定，復(fù)合材料的微觀理化過程無法通過實(shí)驗(yàn)從分子層面獲取，分子模擬技術(shù)相比而言在微觀研究方面有得天獨(dú)厚的優(yōu)勢與潛力[12].

在高電壓與絕緣技術(shù)領(lǐng)域內(nèi)，常利用分子模擬技術(shù)仿真計算絕緣體系在特定溫度或電場等狀況下的反應(yīng)過程，從微觀機(jī)制角度去揭示絕緣材料的劣化與損壞過程.以分子動力學(xué)(MD)為依據(jù)對聚乙烯/蒙脫土納米復(fù)合材料的研究發(fā)現(xiàn)，蒙脫土改性聚乙烯會產(chǎn)生氫鍵作用，使材料結(jié)構(gòu)更加穩(wěn)定，同時力學(xué)性能也有所改善[8].王等人[10]同樣在實(shí)驗(yàn)的基礎(chǔ)上進(jìn)一步借助分子動力學(xué)模擬計算納米顆粒和材料分子鏈之間的相互作用，并研究納米顆粒對材料內(nèi)部氧氣擴(kuò)散速率的影響，發(fā)現(xiàn)納米顆粒通過抑制重結(jié)晶和氧氣擴(kuò)散的作用延后了熱氧老化反應(yīng)發(fā)生的時間，起到了熱氧老化防護(hù)的效果.凡等人[13]研究了LDPE、納米SiO2以及水分的復(fù)合模擬模型，結(jié)果表明水分會降低材料的玻璃轉(zhuǎn)化溫度，納米SiO2的加入則會提升其玻璃轉(zhuǎn)化溫度，增強(qiáng)復(fù)合材料熱穩(wěn)定性.然而對于低密度聚乙烯摻雜MgO后，其熱穩(wěn)定性以及水分子滲透能力變化的研究鮮有報道.

本文基于MD模擬方法對比研究了不同納米顆粒半徑的MgO團(tuán)簇改性LDPE復(fù)合體系的玻璃化溫度、分子鏈運(yùn)動以及水分子的滲透系數(shù)，研究結(jié)果可為LDPE的熱老化實(shí)驗(yàn)和抑制水樹枝等方面提供參考.

2 計算模型和方法

高壓電力電纜的主要絕緣材料為聚乙烯，聚乙烯則是由結(jié)晶區(qū)和無定形區(qū)兩種存在形式構(gòu)成的兩相結(jié)構(gòu).結(jié)晶區(qū)的微晶因排列緊密，相互之間結(jié)合力大，從而具有較高的機(jī)械強(qiáng)度，相對于無定形區(qū)更為穩(wěn)定不易被破壞，故建立納米氧化鎂摻雜聚乙烯模型時優(yōu)先考慮無定形區(qū).一般而言，聚乙烯因?yàn)槲孕〉奶匦裕肿硬灰诐B透到聚乙烯絕緣層中，然而無定形區(qū)是最容易吸水和融入雜質(zhì)的地方，在潮濕的環(huán)境中，水分在電場的作用下會經(jīng)材料的無定形區(qū)滲入到材料內(nèi)部的自由體積空穴中逐漸形成水樹枝，導(dǎo)致電力設(shè)備尤其是高壓電力電纜絕緣破壞.因此本文在構(gòu)建模型時，為了對比研究MgO納米團(tuán)簇和水分對聚乙烯的性能影響，分別建立了不同種類的復(fù)合體系.

2.1 模型的構(gòu)建

為了對比研究MgO納米團(tuán)簇和水分對聚乙烯的影響，綜合考慮實(shí)際情況和模擬需要耗費(fèi)的時間，本文采用Materials Studio 17.1構(gòu)建和計算了不同類型的聚乙烯無定形區(qū)模型，分別為純聚乙烯(PE)模型、MgO/PE模型、含有相同水分含量的PE模型和MgO/PE模型.不失一般性，為了計算方便，基于一個聚合度為120的聚乙烯分子單鏈利用Amorphous Cell模塊構(gòu)建無定形區(qū)復(fù)合模型.綜合考慮構(gòu)建模型的尺寸以及MgO納米團(tuán)簇在模型中所占的質(zhì)量分?jǐn)?shù)和自由體積，同時為了對比分析不同納米顆粒半徑的MgO對復(fù)合模型微觀特性的影響，納米MgO的顆粒半徑分別設(shè)置為3 ?和5 ?.為使模擬結(jié)果不僅能反映真實(shí)聚乙烯的性質(zhì)，并盡量減少鏈端作用，遂以5條聚乙烯單鏈構(gòu)建的基體進(jìn)行摻雜.在實(shí)際中，初始制備的納米MgO表面存在不飽和的化學(xué)鍵，置于空氣中會與O2和H2O發(fā)生反應(yīng)，從而處于穩(wěn)定狀態(tài).為了接近實(shí)際情況，需要對納米團(tuán)簇進(jìn)行氧化處理，即對納米團(tuán)簇表面不飽和的O原子進(jìn)行H處理，不飽和的MgO原子進(jìn)行-OH處理.有相關(guān)研究表明，聚乙烯材料內(nèi)的水分含量可達(dá)1%以上[14].因此在構(gòu)建模型時考慮水分子在模型中所占的自由體積和質(zhì)量分?jǐn)?shù)，并且為了避免因水分含量不同而造成的影響，需要確保每個含水的復(fù)合模型中水分所占的質(zhì)量分?jǐn)?shù)均約為1%左右，即在復(fù)合模型中分別加入10個水分子.則構(gòu)建PE、PE-10H2O、PE-MgO(3 ?)、PE-MgO(3 ?)-10H2O、PE-MgO(5 ?)、PE-MgO(5 ?)-10H2O六個復(fù)合體系模型，如圖1所示.

2.2 MD模擬

為了使建立的六個模型達(dá)到最穩(wěn)定的初始狀態(tài)，首先，使用Forcite模塊中的Geometry Optimization對復(fù)合模型進(jìn)行結(jié)構(gòu)優(yōu)化以及能量最小化計算，選擇正則系綜(NVT)和等溫等壓系綜(NPT)依次對體系進(jìn)行一次分子動力學(xué)弛豫使其處于平衡狀態(tài)，2次模擬的時長為200 ps，溫度為300 K，壓強(qiáng)為0.1 MPa.在此工況下對復(fù)合體系進(jìn)行從100 K到600 K的5次循環(huán)退火處理，溫度間隔為50 K，每個循環(huán)周期的模擬時間為100 ps，系綜為NVT，并在每一步退火處理中都進(jìn)行模型能量最小化.經(jīng)過以上處理，使得構(gòu)建的無定形模型中產(chǎn)生的局部最小值的勢阱得以消除，使材料構(gòu)型更加合理，并為進(jìn)行下一步分子動力學(xué)模擬提供了比較合理的平衡幾何構(gòu)象.然后提取退火處理后能量最小的結(jié)構(gòu)構(gòu)象，先進(jìn)行200 ps的NVT系綜下的動力學(xué)模擬，將模型升溫至600 K，再進(jìn)行200 ps的NPT系綜下的動力學(xué)模擬獲取此溫度下的參數(shù)，隨后以50 K為梯度使用NVT系綜按100 ps的模擬時長將模型的熱力學(xué)溫度依次降溫至100 K，并在相同的目標(biāo)溫度下又進(jìn)行200 ps的NPT模擬.在以上過程中，積分步長均設(shè)置為1 fs，壓強(qiáng)設(shè)置為0.1 MPa，體系內(nèi)各原子的動力學(xué)信息每隔200 fs收集一次.考慮體系在新溫度點(diǎn)下需要重新平衡的問題，遂提取NPT系綜下后100 ps的模擬數(shù)據(jù)用于分析.

整個模擬采用COMPASS力場[15]，此力場是把以往分別處理的有機(jī)分子體系的力場與無機(jī)分子體系的力場統(tǒng)一的分子力場，已被證實(shí)非常適合有機(jī)分子和無機(jī)分子構(gòu)成的體系.溫度控制采用Anderse熱浴法[16]，壓強(qiáng)控制采用Berendsen壓浴法[17]，牛頓運(yùn)動方程的求解采用Velocity Verlet跳蛙算法.對于非鍵相互作用，范德華作用采用Atom Based方法[18]，靜電作用采用Ewald法[19].

3 仿真結(jié)果及分析

3.1 玻璃化溫度分析

很多高分子材料都存在著結(jié)晶區(qū)和無定形區(qū)，其無定形區(qū)就屬于玻璃態(tài)物質(zhì).玻璃化轉(zhuǎn)變溫度(Tg)是指非晶區(qū)中聚合物由“玻璃態(tài)”轉(zhuǎn)變?yōu)椤案邚棏B(tài)”的一個轉(zhuǎn)折溫度，也是無定形聚合物大分子鏈段自由運(yùn)動的最低溫度.一旦聚合物開始由平衡轉(zhuǎn)變?yōu)殒湺芜\(yùn)動，其多項(xiàng)性質(zhì)也會發(fā)生顯著變化，例如自由體積、膨脹系數(shù)、彈性模量、介電常數(shù)等，因此玻璃化轉(zhuǎn)變溫度是高聚物材料穩(wěn)定性的重要評判指標(biāo).微觀分子動力學(xué)模擬中確立玻璃化溫度Tg最常用的可靠方法是比體積-溫度曲線法，該方法采用NPT模擬取得各溫度點(diǎn)下體系的比體積對溫度數(shù)據(jù)作散點(diǎn)圖，然后以拐點(diǎn)兩側(cè)的數(shù)據(jù)作分段線性擬合分別得到兩條直線，直線交點(diǎn)的橫坐標(biāo)即為Tg，分析得到模型中的玻璃化溫度如圖2所示.

由圖2可知，純聚乙烯的Tg值比文獻(xiàn)[20]低8 K，這主要是模擬時間和力場不同造成的，且與聚乙烯建模時初始構(gòu)象的隨機(jī)性有關(guān)，但并不影響對所有體系的對比分析.隨著溫度的升高，模型中比體積均在增大，且在玻璃轉(zhuǎn)化溫度附近出現(xiàn)跳變現(xiàn)象，即玻璃態(tài)向高彈態(tài)轉(zhuǎn)化，分子鏈段開始運(yùn)動的過程.其中不含水分的模型，即PE、PE-MgO(3 nm)、PE-MgO(5 nm)的Tg值分別為247 K、292 K、301 K，往模型中添加質(zhì)量分?jǐn)?shù)為1%的水分后，其Tg值分別變?yōu)?32 K、274 K、282 K.對比各體系的Tg值可知水分會降低復(fù)合模型的玻璃轉(zhuǎn)化溫度，而摻雜MgO納米顆粒則會提高PE的玻璃轉(zhuǎn)化溫度.玻璃轉(zhuǎn)化溫度的變化對應(yīng)材料熱穩(wěn)定工作溫度的變化，因此從這個意義上講，水分會降低聚乙烯體系的熱穩(wěn)定性，而MgO則很好的提升了聚乙烯的熱穩(wěn)定性.

高分子材料內(nèi)部的體積分為自由體積(VFV)和占有體積(VOV)，被高分子材料自身原子或分子占據(jù)的體積是占有體積，具有原子尺寸的缺陷和由于原子無規(guī)則堆砌而產(chǎn)生的空穴的體積是自由體積.正是由于自由體積的存在，分子鏈才可能通過轉(zhuǎn)動和位移而調(diào)整構(gòu)象.根據(jù)Fox和Flory的自由體積理論[21]和分子運(yùn)動理論可知，空穴是材料中微觀粒子運(yùn)動的重要原因之一，因此可以采用自由體積理論來解釋絕緣材料內(nèi)部分子的熱穩(wěn)定性和機(jī)械性能的變化.本文通過Material Studio軟件中atom volumes & surfaces工具采用硬球探針法計算復(fù)合體系中介質(zhì)即聚乙烯的自由體積分?jǐn)?shù)(FFV)，以PE+H2O為例，如圖3所示.由于原子、分子之間具有相互作用，在此利用范德華半徑為硬球探針，其具體數(shù)據(jù)如表1所示.

表1 300 K下各體系模型的VFV、VOV和FFV

表1為不同復(fù)合模型介質(zhì)的占有體積VOV、自由體積VFV和自由體積分?jǐn)?shù)FFV.從自由體積分?jǐn)?shù)可以看出，水分子的加入使得聚乙烯的自由體積和自由體積分?jǐn)?shù)均增大.自由體積的增大，必然會給聚乙烯的的鏈運(yùn)動增加空間，水分子體積很小，對聚乙烯分子的鏈運(yùn)動造成的空間阻礙也很小，并且因?yàn)榫垡蚁┓肿渔溄Y(jié)構(gòu)無羥基等親水性集團(tuán)，水分子與聚乙烯之間無法產(chǎn)生氫鍵，因此水分的加入會降低聚乙烯的玻璃化溫度，起到增塑劑的作用.而MgO納米團(tuán)簇?fù)诫s到聚乙烯中后，會填補(bǔ)聚乙烯原有的自由體積空隙，從而減小復(fù)合體系的自由體積，同時MgO與聚乙烯分子之間存在著一定的相互作用，都會減弱聚乙烯分子鏈的運(yùn)動，并且氧化鎂表面羥基與水分子之間形成氫鍵，水分子的運(yùn)動強(qiáng)度得到抑制，使得復(fù)合體系的結(jié)構(gòu)更加穩(wěn)定，從而提高了體系的玻璃轉(zhuǎn)化溫度.

3.2 聚乙烯鏈運(yùn)動

聚乙烯在玻璃化轉(zhuǎn)變之前呈玻璃態(tài)，玻璃轉(zhuǎn)變后呈高彈態(tài)，在高彈態(tài)階段，鏈與鏈之間的阻礙作用減小，高分子鏈開始發(fā)生整體運(yùn)動.聚乙烯的機(jī)械性能及熱穩(wěn)定性與聚乙烯的鏈運(yùn)動有很大關(guān)系.聚乙烯鏈運(yùn)動越強(qiáng)烈，聚乙烯的機(jī)械性能及熱穩(wěn)定性則越差.均方位移(MSD)能夠描述鏈的整體運(yùn)動情況，其計算公式如下所示：

MSD=〈|ri(t)-ri(0)|2〉

(1)

式中，ri(t)和ri(0)分別代表t時刻和初始時刻i分子或原子的位置向量，其中〈〉表示系綜平均.

圖4所示為不同溫度下PE模型和PE-MgO(3 nm)模型中聚乙烯鏈的均方位移曲線.由圖4可知隨著溫度的升高，兩個體系聚乙烯鏈的MSD都在逐漸增大，說明溫度的升高使得PE分子鏈的整體運(yùn)動開始活躍起來.同時由圖4可知，兩個體系的MSD均出現(xiàn)了兩次較為明顯的跳變.對于純PE模型第一次跳變出現(xiàn)在300～350 K之間，對于PE-MgO(3 nm)模型第一次跳變出現(xiàn)在350～400 K之間，均高于前述各體系的玻璃化轉(zhuǎn)變溫度，這正好說明模型從玻璃態(tài)向高彈態(tài)轉(zhuǎn)變后聚乙烯分子鏈之間的相互阻礙作用減小，自由體積增大，鏈運(yùn)動的自由度增加，整體移動加劇，從而導(dǎo)致模型的穩(wěn)定性變?nèi)?此外，兩個體系分別在350～400 K和500～550 K出現(xiàn)第二次跳變，這可能對應(yīng)著聚乙烯從高彈態(tài)到粘流態(tài)的轉(zhuǎn)變.從整體來看，PE-MgO(3 nm)模型的均方位移數(shù)值明顯小于PE模型中的均方位移的值，而且前者的均方位移出現(xiàn)一次跳變和二次跳變的溫度范圍均大于后者，進(jìn)一步說明摻雜MgO納米團(tuán)簇能減弱聚乙烯的鏈運(yùn)動強(qiáng)度，從而提升聚乙烯材料的熱穩(wěn)定性.模擬體系中分子鏈的運(yùn)動和玻璃化轉(zhuǎn)變過程之間的密切聯(lián)系正好反映兩者所具有的一致性.

3.3 水分子滲透

3.3.1基礎(chǔ)理論

水分子在復(fù)合模型中的滲透能力可以通過計算水分子的溶解系數(shù)及擴(kuò)散系數(shù)來進(jìn)行表征.為了定量表征水分子在聚合物中的滲透系數(shù)，并且將溶解性及擴(kuò)散速度定量化，滲透系數(shù)為

P=D×S

(2)

式中，P為滲透系數(shù)，表征水分子通過高分子膜難易度的一個重要參數(shù)，D為分子的擴(kuò)散系數(shù)，反映滲透分子與聚合物體系的動力學(xué)特性，S為溶解度系數(shù)，反映滲透分子與聚合物相互作用的熱力學(xué)特性.

(1)溶解度系數(shù)S是一個熱力學(xué)參數(shù)，表示在滲透達(dá)平衡時，溶質(zhì)在聚合物膜中滲透的物質(zhì)量.分子在聚合物膜中的滲透過程一般采用溶解擴(kuò)散理論來描述，聚合物膜的分離性能則由擴(kuò)散系數(shù)和溶解度系數(shù)共同表征.在計算機(jī)模擬中，溶解系數(shù)通常通過蒙特卡洛統(tǒng)計力學(xué)方法[22]求得，有兩種等效的方法進(jìn)行計算.文中采用吸附等溫線法，吸附等溫線法是通過在一定的壓力范圍內(nèi)進(jìn)行一系列的巨正則蒙特卡洛(GCMC)模擬，得到吸附分子的濃度隨壓力變化的曲線，通常表現(xiàn)出非線性行為，可以通過dual-mode模型來描述[23]：

(3)

式中，C是偏壓p下高分子中的平衡滲透物的濃度，其單位是cm3(STP)/cm3(聚合物)，kD是Henry常數(shù)，CH是Langmuir容量參數(shù)，b是Langmuir親和參數(shù)，表示氣體分子吸附/解吸的能力.通過求出壓力為零時的極限斜率，得到溶解系數(shù)S，即：

(4)

(2)分子動力學(xué)中水分子的擴(kuò)散系數(shù)可以通過Einstein關(guān)系式[24]求解得到，公式為

(5)

由于模擬過程中已經(jīng)對N個擴(kuò)散原子數(shù)的MSD值作了平均，因此當(dāng)模擬時間足夠長時，擴(kuò)散系數(shù)可以簡化為如下形式：

D=m/6

(6)

其中m為MSD對t的曲線通過最小二乘擬合得到的直線的斜率.

3.3.2滲透機(jī)理分析

(1)圖5是298 K時水分子在各個復(fù)合膜中的吸附等溫曲線，由圖可以看出，水分子在3種復(fù)合膜中的吸附等溫線趨勢相同.壓力較低時，水分子以單分子形式優(yōu)先在聚合物微孔區(qū)域被吸附，吸附速度較快；隨著壓力升高，水分子開始聚集成簇，且微孔區(qū)趨于飽和，水分子吸附在聚合物稠密區(qū)進(jìn)行，速度變慢，并逐漸達(dá)到吸附平衡.分析圖中各曲線可以看出，當(dāng)壓力相同時，純聚乙烯膜內(nèi)吸附的水分子個數(shù)最多，摻雜3 nm MgO的聚乙烯膜次之，5 nm MgO的最少.壓力不同時，吸附在各個復(fù)合膜上的水分子個數(shù)隨著壓力的增大而增多，當(dāng)壓力增大到一定程度時，吸附分子的個數(shù)隨壓強(qiáng)的變化幾乎不變，逐漸達(dá)到吸附平衡，且達(dá)到吸附平衡時在不同膜中的濃度值也有很大差別.

對其所得的數(shù)據(jù)采用最小二乘法進(jìn)行擬合，結(jié)果如表2所示.表2給出了水分子在3種復(fù)合膜中的雙重吸附模型參數(shù)及計算所得的溶解度系數(shù).水分子在膜中吸附溶解過程受高分子鏈自由體積、水分子與高分子鏈相互作用等多種因素的共同影響.水分子在不同復(fù)合膜中的擴(kuò)散系數(shù)和溶解度系數(shù)值差別較大，且溶解度系數(shù)大小順序與復(fù)合膜自由體積分?jǐn)?shù)大小順序基本相同.

表2 300 K時水分子的雙重吸附模型參數(shù)及溶解度系數(shù)

(2)以PE+H2O在不同溫度下的水分子擴(kuò)散為例，作出均方位移如圖6所示，從圖6中可以發(fā)現(xiàn)水分子在體系中的運(yùn)動比較穩(wěn)定，均方位移曲線與時間呈很好的線性關(guān)系，且隨著溫度的升高，水分子的均方位移和絕對值斜率均逐步上升.

對各體系水分子MSD曲線進(jìn)行擬合求出不同溫度下水分子在各模型中的擴(kuò)散系數(shù)，如表3所示.

表3 不同溫度下各模型中水分子的擴(kuò)散系數(shù)(10-6cm2·s-1)

由表3中的數(shù)據(jù)可知，隨著溫度的升高，水分子熱運(yùn)動加劇，分子間的碰撞越強(qiáng)烈，所以擴(kuò)散系數(shù)變大.水分子在聚合物體系中是不斷隨機(jī)運(yùn)動的，但在大多數(shù)情況下作小幅度的振動，只有在合適的時候從一個自由體積空穴跳躍到另一自由體積空穴.由此可見體系內(nèi)自由體積的大小、分布都對滲透分子的擴(kuò)散有影響.自由體積分?jǐn)?shù)越大，體系內(nèi)的自由體積空穴越多，則滲透分子在體系內(nèi)空穴間的跳躍所需的距離越小，擴(kuò)散越快.這恰好是添加MgO納米團(tuán)簇會使水分子擴(kuò)散系數(shù)降低的原因之一，另一個原因則是MgO納米團(tuán)簇能與水分子形成氫鍵從而抑制其擴(kuò)散運(yùn)動.由圖6還可看出PE體系中水分子的均方位移在250 K～300 K之間有一個跳變，這與其玻璃化溫度有關(guān)，也從側(cè)面證明了玻璃化溫度結(jié)果的準(zhǔn)確性.

(3)對所得的水分子的擴(kuò)散系數(shù)及溶解度系數(shù)進(jìn)行運(yùn)算，可得相應(yīng)的滲透系數(shù)，如表4所示.表4給出了300 K時水分子在各體系中的溶解度系數(shù)、擴(kuò)散系數(shù)及滲透系數(shù).由表中可以看出，在相同溫度下，水分子在三個體系中滲透系數(shù)的大小依次為PE、PE/MgO(3 nm)和PE/MgO(5 nm)，滲透系數(shù)的變化趨勢與溶解度系數(shù)和擴(kuò)散系數(shù)的變化趨勢是一致的，說明在水分子滲透過程中，滲透系數(shù)的大小主要由水分子對聚合物膜的溶解度系數(shù)和擴(kuò)散系數(shù)值共同決定.水分子在摻雜納米MgO的聚乙烯中的滲透能力小于純聚乙烯，這是由于納米MgO的加入占據(jù)了聚合物體系中的部分自由體積，變相減弱了水分子在聚乙烯體系內(nèi)的溶解性與擴(kuò)散能力，從而導(dǎo)致其滲透能力的減弱.

表4 300 K時水分子在3種模型中的溶解度系數(shù)、擴(kuò)散系數(shù)及滲透系數(shù)

4 結(jié) 論

(1)利用經(jīng)典的比體積-溫度曲線法確定了各個聚乙烯復(fù)合模型的玻璃化溫度.對于PE體系而言，摻雜水分使得其玻璃化溫度下降了15 K，而MgO納米團(tuán)簇的加入則使得其玻璃化溫度有所提升，且納米MgO的顆粒半徑越大，則復(fù)合模型的玻璃化溫度提升的也越高，反映出水分會減弱聚乙烯的熱穩(wěn)定性，降低聚乙烯的熱穩(wěn)定溫度工作點(diǎn)，納米MgO則可以增強(qiáng)聚乙烯體系的熱穩(wěn)定性.

(2)聚乙烯分子鏈的均方位移隨著溫度的升高逐漸增大，并且隨著溫度的增加，分子鏈的均方位移增加的幅度也越大，并在不同的溫度范圍內(nèi)均有兩次明顯的跳變，說明溫度的升高不僅能使聚乙烯分子鏈的運(yùn)動加劇，并且能引起聚乙烯分子結(jié)構(gòu)從玻璃態(tài)向高彈態(tài)轉(zhuǎn)變，進(jìn)一步破壞聚乙烯結(jié)構(gòu)的穩(wěn)定性.添加MgO納米團(tuán)簇則能占據(jù)聚乙烯體系中的部分自由體積，有效抑制聚乙烯分子鏈的運(yùn)動，從而提高聚乙烯體系的熱穩(wěn)定性.

(3)隨著溫度的升高，復(fù)合體系中水分子的擴(kuò)散系數(shù)也逐漸增大，復(fù)合體系的自由體積越大，越有利于水分子的擴(kuò)散.水分子在復(fù)合體系中的滲透能力由模擬得到的擴(kuò)散系數(shù)和溶解度系數(shù)共同表征，納米MgO的加入使得復(fù)合體系的自由體積分?jǐn)?shù)變小，從而使其擴(kuò)散系數(shù)和溶解度系數(shù)減小，因此導(dǎo)致水分子對聚乙烯的滲透能力減弱.