退火結(jié)束后，在Forcite模塊下對無定形混合體系模型進(jìn)行弛豫。在COMPASS II力場下，首先進(jìn)行100 p、的正則系綜(NVT)動力學(xué)計算穩(wěn)定體系，然后進(jìn)行250 p、的恒溫恒壓系綜(NPT)動力學(xué)計算去調(diào)節(jié)密度。模擬過程均采用Andersen方法控溫器，BeTendsen方法控壓器，壓力為0.000 1 GPa。各分子起始速度按Maxwell分布取樣，VelocityVPT'let算法!la，進(jìn)行求解。最后在相同條件下進(jìn)行1 000 ps的N PT動力學(xué)計算，時間步長為。保存全運(yùn)動軌跡，分析不同溫度下增塑劑分子在氯丁密封膠體系中的均方位移(MSD)，然后用同樣的方法分析常溫下不同增塑劑及增塑劑不同含量時在氯丁密封膠體系中的MSD。

        擴(kuò)散系數(shù)是衡量物質(zhì)遷移能力強(qiáng)弱的一項重要參數(shù)，它是指在單位濃度梯度下，擴(kuò)散物質(zhì)通過單位面積的物質(zhì)的量。利用分子動力學(xué)計算擴(kuò)散系數(shù)，常用Einstein關(guān)系式，如式(1)所示!式中為t時刻分子的坐標(biāo);r為初始坐標(biāo);m為均方位移vISD隨時間t變化曲線的斜率。

        動力學(xué)計算完成后，選取線性較好的100600 ps均方位移軌跡數(shù)據(jù)進(jìn)行擬合，得到直線斜率，代人式(3)求得擴(kuò)散系數(shù)D。由于增塑劑與氯丁密封膠分子結(jié)構(gòu)的差異，兩者存在相容性的問題。通過溶度參數(shù)來判斷兩種物質(zhì)的相容性是一種較為常用且簡便的方法半經(jīng)驗公式如式(4)所示。式中為溶劑與聚合物的體積分?jǐn)?shù)(%)為混合后總體積(mL);s,s2為兩組分的溶度參數(shù)。

    由溶解度參數(shù)相近原則!22，可知:當(dāng)s，與越為接近或相等時，兩種聚合物越容易自發(fā)混合，相容越好。http://zope3.cn/