天然橡膠/順丁橡膠共混物相容性的分子模擬研究

袁曉欽，曹獻(xiàn)超，馬艷彬，于人同，趙艷芳，陳永平，廖建和

天然橡膠/順丁橡膠共混物相容性的分子模擬研究

袁曉欽，曹獻(xiàn)超，馬艷彬，于人同，趙艷芳，陳永平，廖建和*

海南大學(xué)材料科學(xué)與工程學(xué)院，海南海口 570228

天然橡膠（NR）具有優(yōu)異的綜合性能而廣泛應(yīng)用于各個(gè)領(lǐng)域，但較差的耐熱、耐老化性能影響了NR制品的使用壽命。天然橡膠與順丁橡膠（BR）共混可改善制品的拉伸強(qiáng)度、耐磨性、抗撕裂性能，提高制品的使用壽命和耐低溫性能，二者共混常用于輪胎制造。NR與BR共混，二者的相容性對(duì)共混物的相結(jié)構(gòu)和性能具有重要影響。本研究利用Material Studio軟件構(gòu)建模型，采用分子動(dòng)力學(xué)（MD）和耗散粒子動(dòng)力學(xué)（DPD）模擬方法在分子尺度和介觀尺度上對(duì)室溫下NR/BR共混的相容性進(jìn)行探討。通過(guò)MD模擬，確定了NR和BR在模擬過(guò)程的最小聚合度均為20；構(gòu)建10/90、30/70、50/50、70/30、90/10的NR/BR共混物模型，并計(jì)算每個(gè)NR/BR共混模型的相互作用參數(shù)，發(fā)現(xiàn)在NR/BR共混體系中NR與BR分子間的相互作用參數(shù)NR-BR始終比臨界相互作用參數(shù)C小；分析了不同質(zhì)量比NR/BR共混物分子間的徑向分布函數(shù)，發(fā)現(xiàn)NR與BR分子間的徑向分布函數(shù)ginter()始終高于NR與NR、BR與BR同種分子間的徑向分布函數(shù)，表明共混物中不同種類分子間（NR與BR）的相互作用比同種類分子間（NR與NR、BR與BR）的相互作用要強(qiáng)；DPD模擬得到不同質(zhì)量比NR/BR共混物的介觀形貌，發(fā)現(xiàn)共混物分散相的相疇尺寸會(huì)隨著分散相含量的增加而增加，但均未出現(xiàn)宏觀相分離的情況。上述模擬結(jié)果表明，NR和BR在室溫下以不同質(zhì)量比混合時(shí)，始終表現(xiàn)出良好的相容性。為驗(yàn)證模擬結(jié)論是否準(zhǔn)確，制備了NR/BR共混物試樣，并采用動(dòng)態(tài)熱機(jī)械分析儀對(duì)樣品進(jìn)行測(cè)試，結(jié)果證實(shí)了NR/BR共混體系的混溶性。模擬結(jié)論與實(shí)驗(yàn)結(jié)果相匹配，證明分子模擬研究和預(yù)測(cè)共混物相容性是可行的，可為共混實(shí)驗(yàn)提供可靠的參考和指導(dǎo)。

天然橡膠；順丁橡膠；分子動(dòng)力學(xué)；耗散粒子動(dòng)力學(xué)；相容性

天然橡膠（NR），由于其優(yōu)異的彈性、絕緣性、隔水性而廣泛應(yīng)用在輪胎、鞋業(yè)、航空航天等領(lǐng)域。但由于其耐老化、耐候性能較差，通常需要與其他橡膠并以提高其使用性能。順丁橡膠（BR）與天然橡膠相比具有優(yōu)良的耐低溫、耐磨性能，并且生熱低，動(dòng)態(tài)機(jī)械性能好[1-3]。天然橡膠與順丁橡膠共混不僅能有效改善制品的老化和耐磨性能，并且能提升產(chǎn)品的耐低溫效果，當(dāng)前許多工程輪胎配方都用到NR與BR共混[4]。NR、BR都是聚合物，當(dāng)2種聚合物共混時(shí)，涉及到兩組分的相容性問題，完全不相容導(dǎo)致的宏觀相分離以及完全相容導(dǎo)致的均相材料都不是共混所希望的，材料的性能也會(huì)受到影響，因此，研究聚合物共混體系的相容性十分必要[5-6]。

通過(guò)傳統(tǒng)差示掃描量熱法（DSC）、動(dòng)態(tài)熱機(jī)械分析法（DMA）、掃描電鏡（SEM）等方法可研究聚合物共混的相容性，但過(guò)程比較繁瑣且耗時(shí)也較長(zhǎng)，經(jīng)濟(jì)效應(yīng)比較低。分子模擬依賴于計(jì)算機(jī)和模擬軟件，以牛頓運(yùn)動(dòng)方程為基礎(chǔ)描述體系粒子的運(yùn)動(dòng)，從而得到模擬體系的一些基本性質(zhì)和性能。得益于當(dāng)代計(jì)算機(jī)技術(shù)的飛速發(fā)展，分子模擬不斷改進(jìn)與發(fā)展，如今已經(jīng)廣泛應(yīng)用于物理、化學(xué)、生物、材料等領(lǐng)域。使用分子模擬的方法對(duì)聚合物共混進(jìn)行模擬不僅節(jié)能經(jīng)濟(jì)，而且還可以給實(shí)驗(yàn)及配方設(shè)計(jì)提供可靠的參考，有助于共混材料改性和新型功能材料的制備[7-8]。李飛舟等[9-10]用分子動(dòng)力學(xué)和介觀動(dòng)力學(xué)模擬的方法，從Flory-Huggins參數(shù)、相圖、混合能等研究了杜仲膠（EUG）與多種橡膠的相容性，結(jié)果表明，EUG在350 K以上可與氯丁橡膠（CR）、順丁橡膠（BR）、丁苯橡膠（SBR）、天然橡膠（NR）、丁腈橡膠（NBR）混溶，模擬預(yù)測(cè)與實(shí)驗(yàn)結(jié)果基本一致。FU等[11]通過(guò)原子分子動(dòng)力學(xué)（MD）和介觀動(dòng)力學(xué)（MesoDyn）模擬研究了不同重量比的聚丙烯（PP）/聚酰胺-11（PA11）共混物的混溶性，從徑向分布函數(shù)、共混物的介觀結(jié)構(gòu)等方面準(zhǔn)確預(yù)測(cè)了PP/PA11的相容性。劉倩等[12]采用分子動(dòng)力學(xué)模擬和耗散力子動(dòng)力學(xué)模擬的方法，從共混物徑向分布函數(shù)，共混物等密度圖，內(nèi)聚能密度等方面研究了杜仲膠和天然橡膠的相容性，結(jié)果表明天然橡膠和杜仲膠具有良好的相容性。

本研究利用分子動(dòng)力學(xué)（MD）和耗散粒子動(dòng)力學(xué)（DPD）模擬，采用Material Studio 4.0軟件構(gòu)建NR/BR（質(zhì)量比為90/10，70/30，50/50，30/70，10/90）共混物模型，從NR/BR共混物的Flory-Huggins參數(shù)，徑向分布函數(shù)，共混物介觀形貌來(lái)分析不同配比下NR/BR共混體系的相容性。

1 材料與方法

1.1 材料

天然橡膠，泰3#煙片膠；順丁橡膠，BR9000，購(gòu)于中國(guó)石油天然氣股份有限公司；硬脂酸（AR），購(gòu)于阿拉丁試劑公司；氧化鋅（AR），購(gòu)于阿拉丁試劑公司；促進(jìn)劑CZ（AR），購(gòu)于羅恩試劑公司；硫黃（AR），購(gòu)于山東齊魯石化工程有限公司。

1.2 方法

1.2.1 NR/BR共混物試樣制備 試樣配方：硬脂酸1 phr，氧化鋅5 phr，促進(jìn)劑CZ 0.5 phr，硫黃3 phr，NR、BR為變量。

按照試樣配方，制備不同質(zhì)量比（0/100，10/90，30/70，50/50，70/30，90/10，100/0）的NR/BR共混物試樣。首先將天然橡膠置于開放式煉膠機(jī)中塑煉3 min，再將順丁橡膠及各種助劑加入開煉機(jī)中與天然橡膠混煉，得到NR/BR共混膠。根據(jù)測(cè)得的硫化時(shí)間，在平板硫化機(jī)上將NR/BR共混膠在145℃條件下進(jìn)行硫化，制備出NR/BR共混物試樣。

1.2.2 性能測(cè)試 動(dòng)態(tài)熱機(jī)械分析測(cè)試：采用動(dòng)態(tài)熱機(jī)械分析儀（Q850，美國(guó)TA公司），在薄膜拉伸模式下，對(duì)樣品（寬度×厚度為4 mm×1 mm，長(zhǎng)度適宜的長(zhǎng)條矩形式樣）進(jìn)行了DMA測(cè)試。掃描溫度范圍–130～60℃，升溫速率為3℃/min，頻率1 Hz，應(yīng)變幅度為25 μm。

1.2.3 分子動(dòng)力學(xué)（MD）模擬 用分子動(dòng)力學(xué)模擬計(jì)算NR/BR共混物的相互作用參數(shù)，徑向分布函數(shù)。利用Masterial Studio軟件的Visualizer模塊構(gòu)建聚合度為5～50（每5聚合度構(gòu)建一個(gè)模型）NR、BR的分子鏈，用Smart Minimization對(duì)各分子鏈進(jìn)行結(jié)構(gòu)優(yōu)化。然后利用Amorphous Cell模塊構(gòu)建NR、BR三維周期性的無(wú)定形單分子鏈模型，通過(guò)對(duì)模型的結(jié)構(gòu)和能量進(jìn)行優(yōu)化之后，分析不同模型的溶度參數(shù)，建立溶度參數(shù)與聚合重復(fù)單元數(shù)的關(guān)系曲線圖，從而確定模擬過(guò)程N(yùn)R與BR的最小聚合度。

在Materials Studio中對(duì)建好的不同分子模型先用Forcite模塊的Geometry Optimization，方法選Smart Minimization進(jìn)性結(jié)構(gòu)優(yōu)化，然后對(duì)優(yōu)化后的模型進(jìn)行退火處理，溫度梯度50 K，從300 K升到500 K，再?gòu)?00 K下降到300 K，循環(huán)5次，整個(gè)模擬共50 ps，以消除模型內(nèi)局部不合理結(jié)構(gòu)。再將退火后能量最小的結(jié)構(gòu)在298 K條件下進(jìn)行100 ps的NVT、200 ps的NPT平衡，每1 ps保存一次軌跡，后50 ps軌跡用于計(jì)算參數(shù)、分析性能。

整個(gè)分子動(dòng)力學(xué)模擬過(guò)程采用COMPASS[13]力場(chǎng)，Nose控溫法，BERENDSEN等[14]研究控壓法，靜電力（electrostatic）和范德華力（van der Waals）分別采用Ewald[15]、Atom based[16]方法計(jì)算。模型的非鍵截?cái)喟霃剑╟utoff distance）設(shè)置為0.95 nm，樣條寬度、緩沖寬度設(shè)置為0.1 nm、0.05 nm，時(shí)間步伐為1 fs（圖1）。

圖1 NR/BR共混物模型構(gòu)建流程

1.2.4 耗散粒子動(dòng)力學(xué)（DPD）模擬 DPD模擬是用一系列珠子代表體系中的原子簇，利用柔性勢(shì)函數(shù)計(jì)算體系能量，并通過(guò)運(yùn)動(dòng)方程和3種作用力來(lái)描述珠子的運(yùn)動(dòng)[17]。因此，在做DPD模擬之前需要對(duì)聚合物分子鏈按下式進(jìn)行粗粒化處理。

式中，NDPD為DPD模擬中分子鏈中的珠子數(shù)，是聚合物分子量，n是聚合物的聚合度，為極限特征比，是分子鏈的一個(gè)固有屬性，利用MS的Synthia模塊求出NR、BR的極限特征比分別為5.45和5.60近似于5，即在DPD模擬中NR、BR分子鏈中5個(gè)重復(fù)單元可用一個(gè)珠子表示，分子鏈的粗粒化示意圖如圖2所示。

DPD屬于介觀尺度上的模擬，模擬的體系與分子動(dòng)力學(xué)模擬相比要大得多，故在模擬過(guò)程中設(shè)定的分子量要大一些，具體的粗粒化參數(shù)如表1所示。

表1 NR、BR粗粒化參數(shù)

模擬體系兩組分珠子間的作用參數(shù)見表2。

表2 DPD模擬參數(shù)

2 結(jié)果與分析

2.1 聚合度與溶度參數(shù)的關(guān)系分析

分子模擬中計(jì)算分子數(shù)量和分子鏈長(zhǎng)度越大，計(jì)算時(shí)間就越長(zhǎng)，對(duì)計(jì)算機(jī)或服務(wù)器要求就越高，選擇合適而相對(duì)準(zhǔn)確的聚合度進(jìn)性模擬可達(dá)到事半功倍的效果。構(gòu)建不同重復(fù)單元的NR、BR分子鏈模型，對(duì)不同的模型進(jìn)行結(jié)構(gòu)和能量?jī)?yōu)化后計(jì)算出不同重復(fù)單元模型的溶度參數(shù)，最終得到了NR、BR聚合度和溶度參數(shù)的關(guān)系（圖3）。

從圖3可看到，NR、BR的溶度參數(shù)一開始比較大，隨著聚合度的逐步增加而下降。當(dāng)聚合度超過(guò)20以后，2種聚合物的溶度參數(shù)變化幅度減小并趨于穩(wěn)定，并且浮動(dòng)的范圍在NR、BR的溶度參數(shù)的實(shí)際范圍之內(nèi)，NR溶度參數(shù)范圍為16.2～17.0(J/cm3)1/2，BR溶度參數(shù)范圍為16.6～ 17.2(J/cm3)1/2[18]。因此，20聚合度能代表NR、BR最小鏈長(zhǎng)，故模擬實(shí)驗(yàn)均以20聚合度來(lái)構(gòu)建NR、BR共混物模型。

圖3 NR、BR聚合度和溶度參數(shù)的關(guān)系

2.2 NR/BR共混物Flory-Huggins相互作用參數(shù)χ分析

二元聚合物共混體系的相容性可通過(guò)Flory- Huggins相互作用參數(shù)來(lái)預(yù)測(cè)[19]。根據(jù)Flory-Huggins理論，可求解出二元聚合物共混的相互作用參數(shù)AB和臨界相互作用參數(shù)C。

比較兩組分間的相互作用參數(shù)AB與臨界相互作用參數(shù)C的大小可判斷二者之間的相容性。當(dāng)AB>C時(shí)，共混體系中兩組分不相容，發(fā)生相分離；當(dāng)AB≤C時(shí)，共混體系的兩組分具有相容性。

前面的實(shí)驗(yàn)已經(jīng)確定了能代表NR、BR的最小聚合度為A和B均為20，構(gòu)建NR/BR（0/100、10/90、70/30、50/50、30/70、90/10、100/0）共混模型，并進(jìn)行結(jié)構(gòu)優(yōu)化，對(duì)優(yōu)化模型進(jìn)行分析求出體積分?jǐn)?shù)、混合能等參數(shù)（表3）。

表3 共混模型體積分?jǐn)?shù)、混合能參數(shù)

構(gòu)建好的模型進(jìn)行結(jié)構(gòu)和能量?jī)?yōu)化（298 K條件下）后，對(duì)其進(jìn)行分析。將得到的相關(guān)數(shù)據(jù)代入公式（4）、（5）、（6）中，求出不同組分NR/BR共混物的相互作用參數(shù)及臨界相互作用參數(shù)，并建立相互作用參數(shù)與BR質(zhì)量分?jǐn)?shù)的關(guān)系（圖4）。

圖4 NR/BR共混物的相互作用參數(shù)

由圖4可知，NR/BR二者相互作用參數(shù)的臨界值C=0.1，BR的質(zhì)量分?jǐn)?shù)為10%時(shí)，二者相互作用參數(shù)NR-BR最小，隨著BR質(zhì)量分?jǐn)?shù)的增加NR-BR逐步增加；當(dāng)BR升高到30%之后，NR-BR起伏趨于穩(wěn)定。雖然NR-BR隨著BR質(zhì)量分?jǐn)?shù)變化而變化，但其值始終比C小，即不同質(zhì)量比的NR/BR共混體系，都有NR-BR

2.3 NR/BR共混物的徑向分布函數(shù)分析

徑向分布函數(shù)（radial distribution function）是反應(yīng)材料微觀結(jié)構(gòu)的物理量，它表示一個(gè)原子或分子周圍距離的地方出現(xiàn)另一個(gè)原子或分子的概率密度與隨機(jī)分布概率密度的比值。可通過(guò)下式求解：

為了進(jìn)一步探究NR與BR的相容性關(guān)系，對(duì)NR/BR共混體系分子間的徑向分布函數(shù)進(jìn)行了分析。根據(jù)文獻(xiàn)資料的報(bào)道[22]，在共混體系中如果不同分子之間的徑向分布函數(shù)大于同種分子之間的徑向分布函數(shù)，表明不同分子之間的相互作用較大，共混體系具有相容性。反之，同種分子間的徑向分布函數(shù)大于不同分子間的徑向分布函數(shù)，表明通用分子間相互作用大，發(fā)生聚集現(xiàn)象，共混體系不相容。

圖5是不同NR/BR共混體系分子間的徑向分布函數(shù)。由圖5可見，不同質(zhì)量比共混的NR/BR體系都有NR-BR分子間的徑向分布函數(shù)大于同種分子NR-NR、BR-BR間的徑向分布函數(shù)，表明不同共混體系都有NR與BR的相互作用大于NR與NR、BR與BR的相互作用，說(shuō)明NR與BR在不同質(zhì)量比共混條件下均具有良好的相容性。

2.4 DPD模擬結(jié)果分析

圖6是DPD模擬得到的NR/BR共混物粒子的介觀分散圖（A）和等密度圖（B）。為了便于區(qū)分，分散圖中紅色珠子代表NR、綠色代表BR。等密度圖中紅色代表NR等值面，灰色代表BR等值面。觀察圖6A，NR（90/10，70/30）或BR（10/90、30/70）為分散相時(shí)，分散相粒子分散均勻；NR/BR為50/50時(shí)，BR粒子存在輕微的聚集現(xiàn)象。而由圖6B可見，當(dāng)NR為分散相時(shí)，隨著NR含量的升高，NR等值面（灰色面積增大）增加，即分散相相疇變大；同理，BR為分散相時(shí)，分散相相疇尺寸也隨著BR含量增加而增加（紅色面積增大），相疇增大表明粒子聚集程度升高，相容性下降。但是，不同比例的NR/BR共混物的相界面模糊，均未出現(xiàn)明顯相分離現(xiàn)象，表明不同質(zhì)量比的NR/BR共混體系是相容的。

2.5 NR/BR共混物DMA曲線分析

眾所周知，動(dòng)態(tài)熱機(jī)械分析是研究聚合物共混相容性的一種重要手段，借助動(dòng)態(tài)熱機(jī)械分析儀可得到聚合物共混物的玻璃化轉(zhuǎn)變溫度（g）。當(dāng)2種聚合物共混時(shí)，如果能夠以任意比例互溶，則共混物只出現(xiàn)一個(gè)g；若2種聚合物的相容性很差，共混時(shí)發(fā)生相分離，共混物表現(xiàn)出2個(gè)g，分別處于純物質(zhì)g處；而如果兩聚合物具有良好的相容性，但又不能任意比例互溶，共混物則出現(xiàn)2個(gè)g，這2個(gè)g相互靠近且與純物質(zhì)的g有所偏離。

圖5 NR/BR共混物分子間徑向分布函數(shù)

圖6 NR/BR共混物介觀形貌

圖7是NR/BR共混物的DMA曲線圖，tan δ峰值對(duì)應(yīng)的溫度為NR、BR的玻璃化轉(zhuǎn)變溫度g。由圖7可見，NR、BR的g分別為43、98℃。當(dāng)二者以不同質(zhì)量比共混時(shí)，共混物出現(xiàn)2個(gè)相互靠近的g，大小介于NR與BR的g之間，表明不同質(zhì)量比的NR/BR共混體系具有良好相容性，但不能以任意比例互溶，只能部分相容。這一實(shí)驗(yàn)結(jié)果與MD和DPD模擬的結(jié)論相一致，驗(yàn)證了分子模擬研究和預(yù)測(cè)聚合物共混相容性是可行、有效的，可為實(shí)驗(yàn)提供可靠性的參照和指導(dǎo)。

3 討論

聚合物共混改性是開發(fā)具有多種優(yōu)良性能聚合物新材料的一種相對(duì)經(jīng)濟(jì)的、重要的手段。2種聚合物共混，組分間的相容性對(duì)共混物的相結(jié)構(gòu)和性能起決定作用[23-24]。因此，研究聚合物共混物的相容性，對(duì)聚合物共混材料合成和制備具有重要意義。

圖7 NR/BR共混物DMA曲線

天然橡膠常與順丁橡膠共混，用于輪胎胎面膠的制造，二者的相容性對(duì)產(chǎn)品的性能也有一定影響。可通過(guò)溶脹法測(cè)量溶度參數(shù)或者差示掃描量熱法測(cè)量玻璃化轉(zhuǎn)變溫度等實(shí)驗(yàn)方法對(duì)NR/BR共混物相容性進(jìn)行判斷，但是耗時(shí)耗材。分子模擬可實(shí)現(xiàn)對(duì)聚合共混體系粒子分布、形態(tài)結(jié)構(gòu)的變化、相分離等準(zhǔn)確把控，還能從分子、原子尺度上分析共混物的結(jié)構(gòu)和性能，不僅對(duì)聚合物共混原理、相分離機(jī)制及共混物結(jié)構(gòu)性能的研究起到關(guān)鍵的指導(dǎo)作用，還減少實(shí)驗(yàn)上人力、財(cái)力和時(shí)間的消耗。目前分子模擬在聚合共混體系的應(yīng)用非常廣泛，許多科研工作者利用分子模擬研究聚合物共混體系取得了良好的效果[25-27]。但是大部分集中在塑料共混體系或其他復(fù)合材料體系，目前有關(guān)橡膠共混體系相容性的分子模擬研究較少。因此，利用分子模擬的方法構(gòu)建合理的分子模型、準(zhǔn)確研究橡膠共混體系的相容性仍是一項(xiàng)較大的挑戰(zhàn)。

本研究利用Materials Studio構(gòu)建NR/BR共混物模型，采用分子動(dòng)力學(xué)模擬方法計(jì)算NR/BR共混物模型分子間的相互作用參數(shù)、徑向分布函數(shù)，實(shí)現(xiàn)在分子尺度上對(duì)NR/BR共混物的相容性進(jìn)行探討。分子動(dòng)力學(xué)模擬結(jié)果顯示，NR/BR共混體系中NR-BR分子間的相互作用參數(shù)始終比臨界相互作用參數(shù)c小，并且NR與BR分子間的徑向分布函數(shù)比同種類分子間的徑向分布函數(shù)高，表明NR/BR共混體系具有相容性。

為進(jìn)一步研究NR/BR共混體系的相容性，采用耗散粒子動(dòng)力學(xué)模擬方法對(duì)NR/BR共混物的粒子分散情況及介觀相貌進(jìn)行研究。結(jié)果發(fā)現(xiàn)，NR/BR共混體系粒子分散比較均勻，隨著分散相含量增加，分散相粒子出現(xiàn)輕微的聚集現(xiàn)象，共混物分散相相疇尺寸增加，體系相容性下降，但共混物始終未出現(xiàn)明顯的相分離現(xiàn)象，表明NR/BR共混體系存在相容性，但不能互溶，只能部分相容。

為驗(yàn)證模擬方法和模擬結(jié)論是否可行，制備了不同質(zhì)量比NR/BR共混物試樣，并對(duì)試樣進(jìn)行動(dòng)態(tài)熱機(jī)械分析。結(jié)果顯示NR/BR共混物出現(xiàn)2個(gè)相互靠近的玻璃化轉(zhuǎn)變溫度，且偏離純NR、BR的玻璃化轉(zhuǎn)變溫度，說(shuō)明不同質(zhì)量比的NR/BR共混體系具有相容性，存在部分相容的現(xiàn)象。實(shí)驗(yàn)結(jié)果與模擬的結(jié)論相一致，證明本研究構(gòu)建的模型及模擬方法是合理且可行的。綜上所述，分子模擬可用于研究和預(yù)測(cè)橡膠共混體系的相容性，對(duì)橡膠共混配方的設(shè)計(jì)和新型共混材料的開發(fā)具有重要指導(dǎo)意義。

理論、實(shí)驗(yàn)、模擬是科學(xué)發(fā)現(xiàn)的三大支柱，隨著計(jì)算機(jī)技術(shù)的不斷變革和創(chuàng)新，分子模擬的準(zhǔn)確度越來(lái)越高，方法層出不窮。并且與傳統(tǒng)的實(shí)驗(yàn)相比，分子模擬只需在計(jì)算機(jī)上進(jìn)行操作，具有成本低、安全性高等優(yōu)點(diǎn)，還可研究化學(xué)反應(yīng)、相結(jié)構(gòu)演變等快速發(fā)生的現(xiàn)象，模擬得到的結(jié)果準(zhǔn)確度也比較高，可進(jìn)一步對(duì)現(xiàn)象或問題的了解，有利于理論的解釋和新材料的開發(fā)，未來(lái)分子模擬在材料科學(xué)領(lǐng)域的應(yīng)用將越來(lái)越廣泛。

[1] 趙艷芳, 鄭詩(shī)選, 劉 丹, 林升博, 廖雙泉. 天然橡膠共混改性最新研究綜述[J]. 熱帶作物學(xué)報(bào), 2014, 35(2): 413-418.

ZHAO Y F, ZHENG S X, LIU D, LIN S B, LIAO S Q. New advance in study on NR blend system[J]. Chinese Journal of Tropical Crops, 2014, 35(2): 413-418. (in Chinese)

[2] 龐 松, 徐新建, 陳家輝, 吳友平. 硫化時(shí)間對(duì)NR/BR并用膠中交聯(lián)密度的分布以及膠料耐磨性的影響[J]. 彈性體, 2019, 29(5): 11-16.

PANG S, XU X J, CHEN J H, WU Y P. Prepration and properties of water-swellable rubber based on hydrogenated nitrile rubber covalently grafted with polyether-PU polymer[J]. China Elastomerics, 2019, 29(5): 11-16. (in Chinese)

[3] 鮑中村, 宗成中. NR/BR/TPI并用膠性能之研究[J]. 世界橡膠工業(yè), 2013, 40(12): 3-6.

BAO Z C, ZONG C Z. Study on the properties of NR/BR/ TPI composites[J]. World Rubber Industry, 2013, 40(12): 3-6. (in Chinese)

[4] 都昌澤, 李玲麗, 陳 濤,肖建斌. 天然橡膠/順丁橡膠并用膠抗撕裂性能[J]. 彈性體, 2019, 29(5): 17- 20.

DU C Z, LI L L, CHEN T, XIAO J B. Improvement of tear resistance of NR/BR blends[J]. China Elastomerics, 2019, 29(5): 17-20. (in Chinese)

[5] 鄭 強(qiáng), 趙鐵軍. 多相/多組分聚合物動(dòng)態(tài)流變行為與相分離的關(guān)系[J]. 材料研究學(xué)報(bào), 1998(3): 225-232.

ZHENG Q, ZHAO T J. Dynamic rheologucal behavior phase-separation of polymers[J]. Chinese Journal of Materials Research, 1998(3): 225-232. (in Chinese)

[6] 徐鴻升, 李忠明, 楊鳴波. 多組分聚合物的動(dòng)態(tài)流變特性[J]. 高分子材料科學(xué)與工程, 2004(6): 24- 28.

XU H S, LI Z M, YANG M B. Dynamic rheology for multi-component polymer blends[J]. Polymeric Materials Science and Engineering, 2004(6): 24-28. (in Chinese)

[7] 付一政. 固體推進(jìn)劑和高分子共混物的微觀、介觀和宏觀多尺度模擬研究[D]. 太原: 中北大學(xué), 2012.

FU Y Z. Multiscale modeling of solid propellants and polymer blends: from the atomistictothe mesoscopic and macroscopic scales[D]. Taiyuan: North University of China，2012. (in Chinese)

[8] SETHI S K, SONI L, SHANKAR U, CHAUHAN R P, MANIK G. A molecular dynamics simulation study to investigate poly (vinyl acetate)-poly (dimethyl siloxane) based easy-clean coating: an insight into the surface behavior and substrate interaction[J]. Journal of Molecular Structure, 2020, 1202: 127342.

[9] 李飛舟. 杜仲膠與常用橡膠的相容機(jī)制及性能研究[D]. 西安: 西安理工大學(xué), 2019.

LI F Z. Study on compatible mechanism and properties of blend systems from eucommia ulmoides gum and common rubber[D]. Xi’an: Xi’an University of Technology, 2019. (in Chinese)

[10] LI F Z, LU Z L, TIAN D P. A combined experimental and molecular dynamics simulation study on the miscibility of eucommia ulmoides gum with several rubbers[J]. Polymers and Polymer Composites, 2017, 25(1): 87-91.

[11] FU Y Z, LIAO L Q, LAN Y H, YANG L X, MEI L Y, LIU Y Q, HU S Q. Molecular dynamics and mesoscopic dynamics simulations for prediction of miscibility in polypropylene/polyamide-11 blends[J]. Journal of Molecular Structure, 2012, 1012: 113-118.

[12] 劉 倩, 岳 紅, 江 浩, 陳 沖. 杜仲膠/天然橡膠共混物的分子動(dòng)力學(xué)模擬和耗散粒子動(dòng)力學(xué)模擬[J]. 材料導(dǎo)報(bào), 2012, 26(6): 141-145.

LIU Q, YUE H, JIANG H, CHEN C. Molecular dynamics simulation and dissipative particle dynamics simulation of eucommia gum/natural rubber blends[J]. Materials Reports, 2012, 26(6): 141-145. (in Chinese)

[13] SUN H. COMPASS: an ab initio force-field optimized for condensed-phase applications overview with details on alkane and benzene compounds[J]. Journal of Physical Chemistry B, 1998, 102(38): 7338-7364.

[14] BERENDSEN H J C, POSTMA J P M, VAN GUNSTEREN W F, DINOLA A, HEAAK J R. Molecular dynamics with coupling to an external bath[J]. The Journal of Chemical Physics, 1984, 81(8): 3684-3690.

[15] EWARD P. Evaluation of optical and electrostatic lattice potentials[J].Annalen Der Physik, 1921, 369: 253-287.

[16] KARASAWA N, GODDARD III W A. Force fields, structures, and properties of poly (vinylidene fluoride) crystals[J]. Macromolecules, 2002, 25(26): 7268-7281.

[17] 苑世領(lǐng), 張 恒, 張冬菊. 分子模擬: 理論與實(shí)驗(yàn)[M]. 北京: 化學(xué)工業(yè)出版社, 2016: 155-166.

YUAN S L, ZHANG H, ZHANG D J. Molecular simulation theory and experiment[M]. Beijing: Chemical Industry Press, 2016: 155-166. (in Chinese)

[18] 吳其曄, 張 萍, 楊文君, 劉潤(rùn)雄. 高分子物理學(xué)[M]. 北京: 高等教育出版社, 2011: 272-277.

WU Q Y, ZHANG P, YANG W J, LIU R X. Polymer physics[M]. Beijing: Higher Education Press, 2011: 272-277. (in Chinese)

[19] DE GENNES P G. Scaling concepts in polymer physics[M]. Ithaca: Cornell University Press, 1979: 69-91.

[20] 游長(zhǎng)江. 橡膠并用與橡塑共混[M]. 北京: 化學(xué)工業(yè)出版社, 2014: 63-65.

YOU C J. Combined use of rubber and blending of rubber and plastic[M]. Beijing: Chemical Industry Press, 2014: 63-65. (in Chinese)

[21] 賈紅兵, 王經(jīng)逸. 橡膠材料學(xué)[M]. 南京: 南京大學(xué)出版社, 2018: 99-100.

JIA H B, WANG J Y. Rubber material science[M]. Nanjing: Nanjing University Press, 2018: 99-100. (in Chinese)

[22] GESTOSO P, BRISSON J. Towards the simulation of poly (vinyl phenol) / poly (vinyl methyl ether) blends by atomistic molecular modelling[J]. Polymer, 2003, 44(8): 2321-2329.

[23] HOSOYA R, ITO M, NAKAJIMA K, MORITA H. Coarse- grained molecular dynamics study of styrene-block-isoprene- block-styrene thermoplastic elastomer blends[J]. ACS Applied Polymer Materials, 2022, 4(4): 2401-2413.

[24] LIU D M, LIN Y, GONG K, BO H F, LI D Y, ZHANG Z X, CHEN W D. Phase behavior and interfacial tension of ternary polymer mixtures with block copolymers[J]. RSC Advances, 2021, 11(61): 38316-38324.

[25] PANNUZZO M, HORTA B A C, LA ROSA C, ROSA C L, DECUZZI P. Predicting the miscibility and rigidity of poly (lactic-co-glycolic acid)/polyethylene glycol blends via molecular dynamics simulations[J]. Macromolecules, 2020, 53(10): 3643-3654.

[26] LIU D M, GONG K, LIN Y, LIU T, LIU Y, DUAN X Z. Dissipative particle dynamics study on interfacial properties of symmetric ternary polymeric blends[J]. Polymers, 2021, 13(9): 1516.

[27] LEMOS T, ABREU C, PINTO J C. DPD simulation of homopolymer-copolymer-homopolymer mixtures: effects of copolymer structure and concentration[J]. Macromolecular Theory and Simulations, 2020, 29(4): 2000014.

Molecular Simulation on Miscibility of NR/BR Blends

YUAN Xiaoqin, CAO Xianchao, MA Yanbin, YU Rentong, ZHAO Yanfang, CHEN Yongping, LIAO Jianhe*

School of Material Science and Engineering, Hannan University, Haikou, Hainan 570228, China

Natural rubber (NR) is widely used in various fields for its excellent comprehensive properties, however its poor heat resistance and aging resistance have great influence on the service life of NR products. The blending of natural rubber and-butadiene rubber (BR) could not only improve the tensile strength, wear resistance and tear resistance of the product, but also improve the service life and cold resistance of the product, and the NR/BR blends are widely used in tire manufacturing. When NR and BR are blended, the miscibility between NR and BR has an important influence on the phase structure and properties of the blends. In this study, Material Studio software was used to build the models required for the simulation, and molecular dynamics (MD) and dissipative particle dynamics (DPD) simulation methods were performed to investigate the miscibility of NR/BR blends from the molecular and mesoscopic scales at room temperature. The MD simulation results showed that the minimum repeating units of NR and BR in the NR/BR blend models were both 20; Models of NR/BR blends with mass ratios of 10/90, 30/70, 50/50, 70/30, 90/10 were constructed by Material Studio, by calculating the interaction parameterof each NR/BR blends, it was found that the interaction parameters of NR-BR in NR/BR blends (NR/BR) were always smaller than the critical interaction parameterC; The radial distribution function between molecules of NR/BR blends with different mass ratios were analyzed, and it was found that the radial distribution functioninter() between NR and BR molecules were always higher than that between NR and NR, BR and BR molecules of the same species, which showed that the interaction between different species (NR-BR) in the blends is stronger than the interaction between the same species (NR-NR, BR-BR). In addition, DPD simulations results calculated from the NR/BR blends with different mass ratios showed that the domain size of the dispersed phase of the blends increased with the increase of the dispersed phase content while macroscopic phase separation had always not been detected. Furthermore, the samples of NR/BR blends were prepared, and the miscibility of each blending system was confirmed with the results of dynamic mechanical analysis (DMA). The simulation results show that when NR and BR are blended with different mass ratios at room temperature, they always show good compatibility, and the experimental results are consistent with the results of molecular simulation.For this reason, molecular simulation will be of significance to study and predict the miscibility of polymer blend, and it will provide much more reliable guidance for experiment.

natural rubber;butadiene rubber; molecular dynamics; dissipative particle dynamics; miscibility

S794.1；TQ330.1

A

10.3969/j.issn.1000-2561.2022.11.005

2022-04-11；

2022-05-05

廣東省重點(diǎn)領(lǐng)域研發(fā)計(jì)劃項(xiàng)目（No. 2020B020217001）；海南省重大科技項(xiàng)目（No. ZDKJ2016020）。

袁曉欽（1998—），男，碩士研究生，研究方向：天然橡膠基礎(chǔ)與改性研究。*通信作者（Corresponding author）：廖建和（LIAO Jianhe），E-mail：990359@hainanu.edu.cn。