朱建勇,何兆益

(重慶交通大學(xué)　材料科學(xué)與工程學(xué)院,重慶　400074)

?

抗剝落劑與瀝青相容性的分子動(dòng)力學(xué)研究

朱建勇,何兆益

(重慶交通大學(xué)材料科學(xué)與工程學(xué)院,重慶400074)

采用Materials Studio研究抗剝落劑與瀝青的相容性，通過(guò)構(gòu)建基質(zhì)瀝青分子平均結(jié)構(gòu)模型，利用分子動(dòng)力學(xué)得到穩(wěn)定構(gòu)象，通過(guò)對(duì)分子運(yùn)動(dòng)軌跡數(shù)據(jù)進(jìn)行分析得出基質(zhì)瀝青與瀝青各組分的內(nèi)聚能密度和溶解度參數(shù)。通過(guò)分子動(dòng)力學(xué)方法考察了5種胺類抗剝落劑與瀝青分子的內(nèi)聚能密度和溶解度參數(shù)。結(jié)果表明，各組分溶解度參數(shù)差值都在2.0以內(nèi)，相容性較好;星形抗剝落劑與基質(zhì)瀝青的相容性較好，而線形抗剝落劑與瀝青的相容性較差，與試驗(yàn)結(jié)果較為符合。

道路工程;抗剝落劑;分子動(dòng)力學(xué);相容性;溶解度參數(shù)

ZHU Jian-yong,HE Zhao-yi

(School of Material Science and Engineering,Chongqing Jiaotong University，Chongqing 400074，China)

0　引言

近年來(lái)，由于交通運(yùn)輸?shù)难杆侔l(fā)展，交通量和汽車軸載快速增加，對(duì)于瀝青和瀝青混合料的性能提出了更高的要求，普通基質(zhì)瀝青已難以滿足路用要求。因此通過(guò)在瀝青中添加其他材料，改善瀝青性能成為了一種趨勢(shì)。瀝青改性效果是否良好主要依賴于改性劑與瀝青的相容性?，F(xiàn)有的改性瀝青相容性評(píng)價(jià)體系主要是以相容性試驗(yàn)結(jié)果作為依據(jù)，沒(méi)有足夠的理論加以指導(dǎo)，因此改性劑的配方具有一定的盲目性和局限性。

分子模擬以原子水平的分子模型來(lái)模擬分子的結(jié)構(gòu)與行為，進(jìn)而模擬分子體系的各種物理化學(xué)性質(zhì)。由于其高效性、實(shí)用性、通用性、穩(wěn)定性等特點(diǎn)，成為了當(dāng)今材料設(shè)計(jì)的主要手段?，F(xiàn)有研究表明[1-7]采用分子模擬技術(shù)分析瀝青材料的結(jié)構(gòu)及化學(xué)變化，是一種有效的技術(shù)手段。

本文利用商業(yè)分子模擬軟件Material Studio構(gòu)建了基質(zhì)瀝青的三維(3D)結(jié)構(gòu)模型，并采用分子動(dòng)力學(xué)(MD)方法得到最穩(wěn)定的三維分子結(jié)構(gòu)，對(duì)基質(zhì)瀝青和5種抗剝落劑之間的相容性進(jìn)行了研究。試驗(yàn)結(jié)果表明通過(guò)分子動(dòng)力學(xué)模擬方法，研究抗剝落劑和基質(zhì)瀝青的相容性，對(duì)抗剝落劑的初選及優(yōu)選提供理論指導(dǎo)。

1　相關(guān)理論基礎(chǔ)

相容性是道路瀝青的重要性質(zhì)，只有瀝青具有良好的相容性 ,才能保證瀝青具有優(yōu)異性能，所以有必要研究該模型各組分的相容性情況，從而判斷該瀝青分子平均結(jié)構(gòu)模型是否合理。研究高分子的相容性方法較多,內(nèi)聚能密度和溶度參數(shù)是其中較為簡(jiǎn)便的表征參數(shù)。溶度參數(shù)對(duì)了解、判斷高分子體系的溶解性、相容性時(shí)具有十分重要的參考價(jià)值。為了預(yù)測(cè)瀝青各組分的相容性,采用分子動(dòng)力學(xué)模擬方法對(duì)瀝青和瀝青各組分的內(nèi)聚能密度(CED)和溶解度參數(shù)(δ)進(jìn)行模擬計(jì)算。

由于高分子的溶解過(guò)程是溶質(zhì)分子和溶劑分子互相作用的過(guò)程,在等溫等壓下溶解過(guò)程的自由能變化可寫為：

(1)

式中，ΔGM,ΔHM和ΔSM分別為高分子混合時(shí)的混合自由能、混合焓和混合熵；T為溶解溫度。在高分子和溶劑混合時(shí)，只有當(dāng)ΔGM<0時(shí)，高分子才可能溶解。對(duì)于非極性高分子，溶解過(guò)程一般是吸熱的，即ΔHM>0，所以,要使高分子溶解，必須滿足以下條件，即ΔHM

(2)

式中，φ1，φ2分別為溶質(zhì)和溶劑的體積分?jǐn)?shù)(φ1+φ2= 1)；δ1，δ2分別為溶質(zhì)和溶劑的溶度參數(shù)。由式(2)可見(jiàn)ΔHM總是正值，要保證ΔGM<0 ，必然是ΔHM越小越好，也就是說(shuō)δ1與δ2必須接近或相等，這就是著名的溶度參數(shù)相近原則。

溶解度參數(shù)δ定義為：

(3)

式中，E為凝聚態(tài)物質(zhì)的內(nèi)聚能；V為摩爾體積；E/V為內(nèi)聚能密度。當(dāng)δ1和δ2的差值Δδ增大時(shí)，溶解的傾向降低。一般認(rèn)為，Δδ≤4.10(J/cm3)1/2時(shí)，才可能有較好的相容性。

2　瀝青分子模型

瀝青組分的組成和結(jié)構(gòu)十分復(fù)雜，本文擬將瀝青分成4組分，用平均分子結(jié)構(gòu)表征各組分?，F(xiàn)有研究表明[8-16]，瀝青質(zhì)的基本結(jié)構(gòu)是以多個(gè)芳香環(huán)組成的稠合芳香環(huán)系為核心，周圍連接有若干的環(huán)烷環(huán)，芳香環(huán)和環(huán)烷環(huán)上都還帶有若干個(gè)長(zhǎng)度不一的正構(gòu)烷基側(cè)鏈，其中還含有各種含硫、氮、氧的基團(tuán)；石油膠質(zhì)分子的基本結(jié)構(gòu)是具有稠環(huán)芳烴和脂環(huán)烴,并帶有碳鏈長(zhǎng)度不等的正構(gòu)或異構(gòu)烷基的層狀結(jié)構(gòu)分子,含有硫、氮和氧等各種雜原子和相應(yīng)的基團(tuán)；飽和分主要化學(xué)結(jié)構(gòu)是烷烴、環(huán)烷烴，而芳香分主要化學(xué)結(jié)構(gòu)是芳烴、含硫衍生物。

董喜貴[8]等研究結(jié)果表明瀝青質(zhì)的基本結(jié)構(gòu)單元可以用稠環(huán)芳烴連接環(huán)烷烴和烷基側(cè)鏈并含雜原子的單元來(lái)表示，結(jié)構(gòu)單元之間形成締合體，締合數(shù)為4～6。王大喜等[9]采用量子化學(xué)方法對(duì)石油膠質(zhì)進(jìn)行了優(yōu)化計(jì)算，得到石油膠質(zhì)單層結(jié)構(gòu)SG的優(yōu)化構(gòu)型；齊邦峰[12]等以紫外吸收光譜理論為基礎(chǔ)，通過(guò)剖析勝利渣油中膠質(zhì)的紫外吸收光譜及其二階導(dǎo)數(shù)光譜數(shù)據(jù),與模型化合物的紫外吸收特征相比較,研究了膠質(zhì)化學(xué)結(jié)構(gòu)。兩種膠質(zhì)分子結(jié)構(gòu)模型相比較而言，前者的分子量大且縮合度小，其分子結(jié)構(gòu)中不含雜原子，因此，可以推測(cè)前者膠質(zhì)分子結(jié)構(gòu)模型比后者膠質(zhì)分子結(jié)構(gòu)模型更具合理性。芳香分分子結(jié)構(gòu)可用1，7-二甲基萘代替，該分子模型比較接近實(shí)際路用瀝青的飽和分、芳香分分子結(jié)構(gòu)，且與瀝青質(zhì)形成早期分子結(jié)構(gòu)類似，芳香環(huán)數(shù)量與支鏈數(shù)量使其較適合作為油分和瀝青質(zhì)的中間物。 D.A.Storm[15]等人的研究表明，瀝青烷烴分子鏈長(zhǎng)度分布在C16與C36之間，C22位于這個(gè)區(qū)間中間部分。I.Kowalewski[16]等人的研究指出在瀝青中所含烷烴C22H46的含量最高，而且C22H46軟化點(diǎn)和沸點(diǎn)與多數(shù)瀝青油分軟化點(diǎn)、沸點(diǎn)較為一致。因此，選取二十二烷C22H46作為飽和分的代表分子模型。

綜上所述，本文所選道路瀝青分子各組分的結(jié)構(gòu)模型如圖1～圖4所示。4組分試驗(yàn)結(jié)果如表1所示。

確定物質(zhì)溶解度參數(shù)的可由物質(zhì)的摩爾汽化熱求得，而分子動(dòng)力學(xué)則是通過(guò)計(jì)算分子之間的氣化熱求得分子的溶解度參數(shù)。利用Accelrys公司開發(fā)的Materials Studio 4.0軟件包中的Amorphous Cell對(duì)基質(zhì)瀝青各組分進(jìn)行建模，利用Compass力場(chǎng)可以在很大的溫度、壓力范圍內(nèi)精確地預(yù)測(cè)孤立體系或凝聚態(tài)體系中各種分子的結(jié)構(gòu)、構(gòu)象、振動(dòng)以及熱物理性質(zhì)。以瀝青質(zhì)為例，設(shè)置10個(gè)不同構(gòu)象的Amorphous cell，溫度為298 K，選用Compass力場(chǎng)優(yōu)化模型，由于瀝青質(zhì)、膠質(zhì)分子結(jié)構(gòu)復(fù)雜，進(jìn)行分子動(dòng)力學(xué)計(jì)算平衡態(tài)時(shí)分子所需的馳豫時(shí)間較長(zhǎng)，為了讓Amorphous Cell內(nèi)的分子有足夠大的能量進(jìn)行調(diào)整，得到最低能量構(gòu)象，最大幅度地逼近實(shí)際體系，應(yīng)先在Compass力場(chǎng)下采用非周期性邊界條件進(jìn)行退火模擬。第一次進(jìn)行動(dòng)力學(xué)模擬所設(shè)溫度為1 500 K，動(dòng)力學(xué)過(guò)程采用NPT系綜，步數(shù)設(shè)置為100 000步，進(jìn)行0.10 ns的分子動(dòng)力學(xué)運(yùn)算，每1 000 步輸出一個(gè)構(gòu)象，得到時(shí)間與體系溫度、能量波動(dòng)圖，如圖5、圖6所示。

圖1　瀝青質(zhì)分子結(jié)構(gòu)[8]Fig.1　Molecule structure of asphaltene

圖2　膠質(zhì)分子結(jié)構(gòu)[9]Fig.2　Molecule structure of colloid

圖3　芳香分分子結(jié)構(gòu)[13]Fig.3　Molecule structure of aromatic

圖4　飽和分分子結(jié)構(gòu)[16]Fig.4　Molecule structure of saturate

瀝青種類化學(xué)組分/%AsRArS70#基質(zhì)瀝青11.3226.3743.1519.16

圖5　1 500 K下模擬時(shí)間與溫度波動(dòng)Fig.5　Relation between simulation time and temperature fluctuation at 1 500 K

圖6　1 500 K下模擬時(shí)間與能量波動(dòng)Fig.6　Relation between simulation time and energy f luctuation at 1 500 K

從圖5、圖6可以看到，經(jīng)過(guò)100 000步的分子動(dòng)力學(xué)運(yùn)算后，體系的溫度和能量波動(dòng)很小，體系達(dá)到平衡狀態(tài)，選擇最后一個(gè)構(gòu)象進(jìn)行第二步的動(dòng)力學(xué)運(yùn)算。

第二次進(jìn)行動(dòng)力學(xué)計(jì)算所設(shè)溫度為298 K，采用NVT系綜，步數(shù)設(shè)置為200 000步，進(jìn)行0.2 ns的分子動(dòng)力學(xué)運(yùn)算。每1 000步輸出一個(gè)構(gòu)象，以下兩圖是體系的溫度和能量隨模擬時(shí)間的變化情況，如圖7，圖8所示。

圖7　298 K下模擬時(shí)間與溫度波動(dòng)Fig.7　Relation between simulation time and temperature fluctuation at 298 K

圖8　298 K下模擬時(shí)間與能量波動(dòng)Fig.8　Relation between simulation time and energy fluctuation at 298 K

從以上兩圖可以看到，經(jīng)過(guò)200 000步的分子動(dòng)力學(xué)運(yùn)算后，體系的溫度和能量波動(dòng)很小，體系達(dá)到平衡狀態(tài)，最后一個(gè)構(gòu)象即為最合理的體系構(gòu)象，利用Amorphous Cell的分析模塊可以直接對(duì)分子運(yùn)動(dòng)軌跡數(shù)據(jù)進(jìn)行分析，得到瀝青質(zhì)內(nèi)聚能密度和溶度參數(shù)。根據(jù)表1結(jié)果，按照相應(yīng)比例，得到瀝青分子模型如圖9所示。

圖9　基質(zhì)瀝青的Amorphous Cell 3D模型Fig.9　Amorphous Cell 3D model of base asphalt

依據(jù)以上步驟，可得瀝青與瀝青各組分的內(nèi)聚能密度和溶度參數(shù)的計(jì)算結(jié)果如表2所示。

表2　瀝青與瀝青各組分的內(nèi)聚能密度和溶度參數(shù)

按照溶度參數(shù)相近原則，即Δδ≤4.10(J/cm3)1/2，對(duì)于瀝青和瀝青各組分，不僅要求兩者的δ相近，而且要求兩者的靜電溶度參數(shù)δe和范德華溶度數(shù)δv也要分別相近，這樣各組分間才能較好地相容。瀝青和瀝青各組分間的最大溶度參數(shù)差值如表3所示。

瀝青和瀝青各組分的溶度參數(shù)曲線都比較平緩，而且從表3可以看到瀝青和瀝青各組分間的最大溶度參數(shù)差值均在4.10(J/cm3)1/2以內(nèi)，表明瀝青各組分之間的相容性較好，證明了適用于道路瀝青分子平均結(jié)構(gòu)模型的合理性。

表3　瀝青和瀝青各組分間的最大溶度參數(shù)差值

3　抗剝落劑模型及其與瀝青相容性研究

本文所選用的5種抗剝落劑描述如下：

抗剝落劑A為環(huán)氧化松香二取代二乙烯三胺；抗剝落劑B為環(huán)氧化松香四取代二乙烯三胺；抗剝落劑C為環(huán)氧化十八酸二取代二乙烯三胺；抗剝落劑D為環(huán)氧化十八酸四取代二乙烯三胺；抗剝落劑E為環(huán)氧化松香十八酸二取代二乙烯三胺。

抗剝落劑A，B，C，D，E的分子結(jié)構(gòu)式和標(biāo)準(zhǔn)命名如圖10～圖14所示。

圖10　抗剝落劑A分子結(jié)構(gòu)Fig.10　Molecule structure of anti-stripping agent A

圖11　抗剝落劑B分子結(jié)構(gòu)Fig.11　Molecule structure of anti-stripping agent B

圖12　抗剝落劑C分子結(jié)構(gòu)Fig.12　Molecule structure of anti-stripping agent C

圖13　抗剝落劑D分子結(jié)構(gòu)Fig.13　Molecule structure of anti-stripping agent D

圖14　抗剝落劑E分子結(jié)構(gòu)Fig.14　Molecule structure of anti-stripping agent E

按照溶度參數(shù)相近原則，即要求Δδ≤4.10 (J/cm3)1/2，考察上述5種抗剝落劑的溶度參數(shù)δ、靜電溶度參數(shù)δe、范德華溶度數(shù)δv和瀝青3項(xiàng)溶度參數(shù)的差值情況，

表4　瀝青與5種抗剝落劑的內(nèi)聚能密度和溶度參數(shù)

結(jié)果如表5所示。

胺類試劑的抗剝落性能的根本原因在于胺類試劑是一種陽(yáng)離子型表面活性劑，其極性的陽(yáng)離子端頭與酸性石料表面存在著離子鍵的化學(xué)吸附作用，具有強(qiáng)烈的抗剝落效果。因此，雜原子或基團(tuán)是胺類試劑抗剝落性的關(guān)鍵因素。

根據(jù)5種抗剝落劑的分子結(jié)構(gòu)可以看出，該5種試劑的基本結(jié)構(gòu)類似，雜原子和基團(tuán)都是氮、酯基、羥基，這些極性原子和基團(tuán)與石料的相互作用是產(chǎn)生黏附效果的主要原因。因此，5種抗剝落劑的改性機(jī)理都是相同的。

表5　5種抗剝落劑的各溶度參數(shù)和瀝青溶度參數(shù)的差值

總結(jié)以上觀點(diǎn)，得到如下推論:上述5種抗剝落劑與基質(zhì)瀝青的相容性是影響5種試劑改性效果的最主要因素。

(1)抗剝落劑C(環(huán)氧化十八酸二取代二乙烯三胺)與瀝青的范德華溶度數(shù)差值Δδv控制在4.10(J/cm3)1/2以內(nèi)，但是它的溶度參數(shù)差值Δδ和靜電溶度參數(shù)差值Δδe都已超過(guò)4.10(J/cm3)1/2，這可能是由于抗剝落劑C分子極性較弱，分子鏈?zhǔn)情L(zhǎng)線形結(jié)構(gòu)而難以分散在瀝青中，與瀝青的裹覆能力較差，所以與瀝青的相容性很差，可以預(yù)測(cè)抗剝落劑C的改性效果較差。

(2)對(duì)于抗剝落劑A(環(huán)氧化松香二取代二乙烯三胺)，雖然抗剝落劑A與瀝青的溶度參數(shù)差值Δδ和范德華溶度數(shù)差值Δδv都在4.10(J/cm3)1/2以內(nèi)，但是靜電溶度參數(shù)差值Δδe為4.80(J/cm3)1/2，已經(jīng)超過(guò)了4.10(J/cm3)1/2，所以抗剝落劑A與瀝青的相容性較差，可以預(yù)測(cè)改性效果較差。

阿強(qiáng)趕緊跑下了樓，到樓下一看，哪里見(jiàn)到人影。他心想不好，自己一定是中了調(diào)虎離山計(jì)。他趕緊回家，見(jiàn)桌子上還留了個(gè)條子：“阿強(qiáng)先生，我已經(jīng)贏了一個(gè)回合，這次先拿一個(gè)手機(jī)作紀(jì)念，下次就是你書桌上的那支金筆。江城神偷?！?/p>

(3)抗剝落劑D(環(huán)氧化十八酸四取代二乙烯三胺)和抗剝落劑E(環(huán)氧化松香十八酸二取代二乙烯三胺)的3項(xiàng)溶度參數(shù)與瀝青的3項(xiàng)溶度參數(shù)的差值都在1～2之間，所以與瀝青的相容性良好，這可能是由于抗剝落劑D的分子鏈較長(zhǎng)，星形結(jié)構(gòu)裹覆瀝青的能力較強(qiáng)；對(duì)于抗剝落劑E，可能是由于其極性較強(qiáng)，它在瀝青中溶解時(shí)由于極性作用所以與瀝青的相容性良好。可以預(yù)測(cè)，抗剝落劑D和抗剝落劑E的改性效果良好。

(4)抗剝落劑B(環(huán)氧化松香四取代二乙烯三胺)的3項(xiàng)溶度參數(shù)與瀝青的3項(xiàng)溶度參數(shù)的差值是最小的，即抗剝落劑B與瀝青的相容性是最好的，這可能是由于抗剝落劑B極性較強(qiáng)，而星形結(jié)構(gòu)與瀝青的裹覆能力也很強(qiáng)，所以相容性很好，可以預(yù)測(cè)抗剝落劑B的抗剝落性應(yīng)該是較好的。

(5)線形抗剝落劑A和星形抗剝落劑B都為環(huán)氧化十八酸取代產(chǎn)物，線形抗剝落劑A相容性較差而星形抗剝落劑B相容性較好；線形抗剝落劑C和星形抗剝落劑D都為環(huán)氧化松香取代產(chǎn)物，線形抗剝落劑C相容性較差而星形抗剝落劑D相容性較好。這是由于星形結(jié)構(gòu)裹覆瀝青的能力比線形結(jié)構(gòu)更強(qiáng)，使得星形抗剝落劑與瀝青的相容性比線形抗剝落劑更好。

(6)抗剝落劑A和抗剝落劑C都為線形抗剝落劑，抗剝落劑A的靜電溶度參數(shù)差值超過(guò)4.10(J/cm3)1/2，而抗剝落劑C的總?cè)芏葏?shù)差值和靜電溶度參數(shù)差值均超過(guò)4.10(J/cm3)1/2，說(shuō)明環(huán)氧化松香取代比十八酸取代相容性更好；抗剝落劑D和抗剝落劑E都為星形抗剝落劑，兩者的總?cè)芏葏?shù)差值和各部分溶度參數(shù)差值都在4.10(J/cm3)1/2以內(nèi)，但環(huán)氧化松香取代比十八酸取代相容性更好。這可能是由于十八酸分子鏈過(guò)長(zhǎng),減弱了抗剝落劑在瀝青中的分散程度，造成相容性相對(duì)減弱。這說(shuō)明不管是線形抗剝落劑或者星形抗剝落劑，環(huán)氧化松香取代均比十八酸取代具有更好的相容性。

為了驗(yàn)證分子模擬的正確性，按瀝青質(zhì)量的0.4%加入合成的抗剝落劑，通過(guò)離析試驗(yàn)和水煮法試驗(yàn)進(jìn)行評(píng)價(jià)，試驗(yàn)結(jié)果如表6所示，表明分子模擬結(jié)果與試驗(yàn)結(jié)果相關(guān)性較好。

表6　試驗(yàn)結(jié)果

4　結(jié)論

(1)采用分子動(dòng)力學(xué)的方法考察了5種抗剝落劑的內(nèi)聚能密度和溶度參數(shù)，發(fā)現(xiàn)瀝青與抗剝落劑B，D，E的相容性良好，而與抗剝落劑A，C的相容性較差，預(yù)測(cè)了抗剝落劑B，D，E的改性效果較好而抗剝落劑A，C的改性效果較差。這對(duì)抗剝落劑的初選提供了理論依據(jù)，具有較大的指導(dǎo)意義。

(2)對(duì)環(huán)氧化松香或十八酸取代二乙三胺，總結(jié)出如下3個(gè)結(jié)論：

①瀝青與抗剝落劑B，D，E的相容性良好，而與抗剝落劑A，C的相容性較差，具體相容性順序?yàn)椋嚎箘兟鋭〣>抗剝落劑D>抗剝落劑E>抗剝落劑A>抗剝落劑C；

②星形抗剝落劑與瀝青的相容性較好，而線形抗剝落劑與瀝青的相容性較差；

③不管是線形抗剝落劑或者星形抗剝落劑，環(huán)氧化松香取代均比環(huán)氧化十八酸取代具有更好的相容性。

(3)試驗(yàn)結(jié)果與預(yù)測(cè)結(jié)果能較好符合，趨勢(shì)比較明確，說(shuō)明了分子動(dòng)力學(xué)方法在瀝青材料研究中具有重要的理論支持和現(xiàn)實(shí)指導(dǎo)意義，能夠較好地預(yù)測(cè)瀝青材料的宏觀性能，是一種研究瀝青材料物理化學(xué)性質(zhì)的重要方法。

(4)現(xiàn)有文獻(xiàn)報(bào)道中瀝青的平均分子結(jié)構(gòu)模型有許多，本文根據(jù)四分法從中篩選了各組分的分子結(jié)構(gòu)，構(gòu)建了瀝青分子平均結(jié)構(gòu)模型，其模型選用的合理性有待進(jìn)一步進(jìn)行論證。

[1]楊小震.分子模擬與高分子材料[M].北京：科學(xué)出版社，2002.

YANG Xiao-zhen.Molecular Simulations and Polymers[M].Beijing：Science Press，2002.

[2]唐伯明，丁勇杰，朱洪洲，等.瀝青分子聚集狀態(tài)變化特征研究[J].中國(guó)公路學(xué)報(bào)，2013，26(3)：50-56，76.

TANG Bo-ming，DING Yong-jie，ZHU Hong-zhou，et al.Study on Agglomeration Variation Pattern of Asphalt Molecules[J].China Journal of Highway and Transport,2013，26(3)：50-56，76.

[3]王淑娟，鄭永昌，丁勇杰.基于分子模擬技術(shù)的單體接枝SBS與瀝青相容性研究[J].公路與汽運(yùn)，2013，1(1)：100-102.

WANG Shu-juan，ZHENG Yong-chang，DING Yong-jie.Study on Compatiblity of Monomer Grafting SBS and Asphalt Using Molecular Simulation [J].Highway &Automotive Application，2013，1(1)：100-102.

[4]丁勇杰，周超.分子模擬技術(shù)在復(fù)合材料研究中的應(yīng)用[J].玻璃鋼/復(fù)合材料，2012(增1):297-300.DING Yong-jie，ZHOU Chao.Application of Molecular Simulation in Composite[J].Fiber Reinforced Plastic/Composites，2012(S1):297-300.

[5]曹斌，高金森，徐春明.分子模擬技術(shù)在石油相關(guān)領(lǐng)域的應(yīng)用[J].化學(xué)進(jìn)展，2004，16(2)：291-298.

CAO Bin，GAO Jin-sen，XU Chun-ming.The Applications of Molecular Simulation Technology in the Fields of Petroleum[J].Progress in Chemistry，2004，16(2)：291-298.

[6]周涵，任強(qiáng).分子模擬技術(shù)在石油化工領(lǐng)域的應(yīng)用進(jìn)展[J].計(jì)算機(jī)與應(yīng)用化學(xué)，2006，23(1)：15-19.

ZHOU Han，REN Qiang.The Applications of Molecular Simulation in Petrochemical Field [J].Computers and Applied Chemistry，2006，23(1)：15-19.

[7]任強(qiáng)，代振宇，周涵．重油膠體結(jié)構(gòu)的介觀模擬[J].石油學(xué)報(bào)：石油加工，2013，29(1)：86-94.

REN Qiang，DAI Zhen-yu，ZHOU Han.Mesoscale Simulation on Colloidal Structures of Heavy Oil[J].Acta Petrolei Sinica：Petroleum Process Section，2013，29(1)：86-94.

[8]董喜貴，雷群芳，俞慶森.石油瀝青質(zhì)的NMR測(cè)定及其模型分子推測(cè)[J].燃料化學(xué)學(xué)報(bào)，2004，32(6)：668-672.

DONG Xi-gui，LEI Qun-fang，YU Qing-sen.NMR Determination of Petroleum Asphaltenes and Their Model Molecules Evaluation [J].Journal of Fuel Chemistry and Technology，2004，32(6)：668-672.

[9]王大喜，趙玉玲，潘月秋，等.石油膠質(zhì)結(jié)構(gòu)性質(zhì)的量子化學(xué)研究[J].燃料化學(xué)學(xué)報(bào)，2006，34(6)：690-694.

WANG Da-xi，ZHAO Yu-ling，PAN Yue-qiu，et al.Quantum Chemical Study on the Structure of Petroleum Gums [J].Journal of Fuel Chemistry and Technology，2006，34(6)：690-694.

[8]LESUEUR D.The Colloidal Structure of Bitumen：Consequences on the Rheology and on the Mechanisms of Bitumen Modification[J].Advances in Colloid and Interface Science，2009，145(1)：42-82.

[9]ZHANG L，GREENFIELD M L.Molecular Orientation in Model Asphalts using Molecular Simulation [J].Energy Fuels,2007，21(11)：2-11.

[10]HEADEN T F，BOEK E S.Molecular Dynamics Simulations of Asphaltene Aggregation in Supercritical Carbon Dioxide with and Without Limonene [J].Energy Fuels,2011，25(5)：3-8.

[11]TAREFDER R A，ARISA I.Molecular Dynamic Simulations for Determining Change in Thermodynamic Properties of Asphaltene and Resin Because of Aging [J].Energy Fuels,2011，25(22) ：11-22.

[12]齊邦峰，曹祖賓，陳立仁，等.紫外吸收光譜研究勝利渣油膠質(zhì)、瀝青質(zhì)結(jié)構(gòu)特性[J].石油化工高等學(xué)校學(xué)報(bào)，2001，14(3)：14-17.

QI Bang-feng，CAO Zu-bin，CHEN Li-ren，et al.Study on Structure of Resins and Asphaltenes With U.V.Absorption Spectrum[J].Journal of Petrochemical Universities，2001，14(3)：14-17.

[13]ZHANG L Q，GREENFIELD M L.Relaxation Time，Diffusion，and Viscosity Analysis of Model Asphalt Systems Using Molecular Simulation[J].Journal of Chemical Physics，2007，127(19)：199-215.

[14]BOEK E S，YAKOVLEV D S，HEADEN T F.Quantitative Molecular Representation of Asphaltenes and Molecular Dynamics Simulation of Their Aggregation [J].EnergyFuels，2009，23(3)：1209-1219.

[15]STORM D A，EDWARDS J C，DECANIO S J，et al.Molecular Representations of Ratawi and Alaska North Slope Asphaltenes based on Liquid-and Solid-state NMR[J].Energy Fuels，1994，8(3)：561-566.

[16]KOWALEWSKI I，VAMDEMBROUCKE M，HUC A Y，et al.Preliminary Results on Molecular Modeling of Asphaltenes Using Structure Elucidation Programs in Conjunction with Molecular Simulation Programs[J].Energy Fuels，1996，10(1)：97-107.

Research of Compatiblity of Asphalt and Anti-stripping Agent Using Molecular Dynamics

The compatibility of anti-stripping agent and asphalt is studied using Materials Studio.By using the stable conformation obtained by molecular dynamics,the cohesive energy density and the solubility parameters of the base asphalt and each component in the established base asphalt molecular average structure model are obtained by analysing the data of molecular trajectories.The cohesive energy density and the solubility parameters of 5 amine anti-stripping agents are studied by molecular dynamics.The result shows that (1) the difference values of solubility parameters of the asphalt and each component are all less than 2.0,which demonstrates that their compatiblity satisfactory;(2) the compatiblity of the star-shaped anti-stripping agents in base asphalt is satisfactory while that of the line-shaped anti-stripping agents is poor,which agrees with the test result.

road engineering；anti-stripping agent；molecular dynamics；compatiblity；solubility parameter

2014-09-26

朱建勇(1977- )，男，湖北武漢人，碩士.(751332737@qq.com)

10.3969/j.issn.1002-0268.2016.01.006

U416.217

A

1002-0268(2016)01-0034-07