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

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

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

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