崔培培， 高 才， 楊 鎖， 葉 斌， 唐景春， 倪宜華

（合肥工業(yè)大學(xué) 機(jī)械與汽車(chē)工程學(xué)院，安徽 合肥 230009）

玻璃態(tài)是熱力學(xué)上的非平衡態(tài)，處于該狀態(tài)的物質(zhì)會(huì)經(jīng)歷一個(gè)自發(fā)的向同溫度下平衡態(tài)趨近的連續(xù)過(guò)程，稱(chēng)之為結(jié)構(gòu)松弛［1］。該過(guò)程中伴有過(guò)剩焓的釋放，故也稱(chēng)焓松弛，這是一個(gè)非常緩慢的過(guò)程，有許多實(shí)驗(yàn)手段可用于研究這種現(xiàn)象，如DSC法、動(dòng)態(tài)機(jī)械分析法（DMA）、介電法和熱膨脹測(cè)定法等。DSC具有很高的測(cè)量精度且重復(fù)性好，是研究物質(zhì)玻璃化轉(zhuǎn)變和焓松弛行為的有力工具［1－7］。文獻(xiàn)［1－2，8］結(jié)合該方法，采用Tool－Narayanaswamy－Moynihan （TNM ） 或AGV現(xiàn)象學(xué)模型對(duì)一些物質(zhì)玻璃轉(zhuǎn)變和焓松弛行為進(jìn)行了研究。這2種模型均能很好地描述物質(zhì)玻璃體系的焓松弛現(xiàn)象，但文獻(xiàn)［5－8］發(fā)現(xiàn)，不同現(xiàn)象學(xué)模型對(duì)物質(zhì)焓松弛行為預(yù)測(cè)效果往往存在差異。一些學(xué)者在原有模型的基礎(chǔ)上作了改進(jìn)，針對(duì)聚合物玻璃體系分析了不同現(xiàn)象學(xué)模型的預(yù)測(cè)力。文獻(xiàn)［3－4］在AG協(xié)同松弛理論的基礎(chǔ)上引入松弛極限態(tài)，建立了GR焓松弛模型，繼而運(yùn)用這2種模型擬合的參數(shù)與介電實(shí)驗(yàn)分析的結(jié)果作了對(duì)比［5］；文獻(xiàn)［6］在AG理論的基礎(chǔ)上引入非線性參數(shù)xs，建立了H模型，并比較了H模型與TNM和AGV模型預(yù)測(cè)力大?。?］；針對(duì)聚合物玻璃形成體系，文獻(xiàn)［2］將AGV與改進(jìn)的TNM多參數(shù)模型作了比較，并對(duì)GR與H模型的預(yù)測(cè)力進(jìn)行了對(duì)照分析［7］。目前所有此類(lèi)對(duì)比分析的文獻(xiàn)都無(wú)一例外地針對(duì)聚合物系統(tǒng)，這些模型對(duì)除聚合物以外的玻璃形成系的預(yù)測(cè)能力的差異報(bào)道不多。

本文選擇木糖醇作為模型系統(tǒng)，擬采用AGV、GR和H3種焓松弛現(xiàn)象學(xué)模型，結(jié)合DSC技術(shù)，比較3種模型描述小分子氫鍵玻璃形成系焓松弛行為時(shí)的預(yù)測(cè)能力。

1 數(shù)學(xué)模型

焓松弛是一個(gè)非線性和非指數(shù)性的過(guò)程，其非指數(shù)性質(zhì)由Kohlrausch－Williams－Watts（KWW）延展指數(shù)函數(shù)描述，即

其中，φ（t）為松弛函數(shù)；τ為特征松弛時(shí)間；β（0＜β≤1）為焓松弛過(guò)程非指數(shù)性質(zhì)的無(wú)量綱量，用于表征松弛時(shí)間分布寬度。

焓松弛不僅依賴于溫度變化，還與體系結(jié)構(gòu)變化有關(guān)，Tool引入Tf作為玻璃化物質(zhì)的結(jié)構(gòu)參數(shù)。Scherer等基于AG理論，給出了AGV模型焓松弛時(shí)間的表達(dá)式為：

其中，A為指前因子，可看成物質(zhì)在溫度無(wú)限高時(shí)的松弛時(shí)間；B＝ScΔμ／k，Sc為最小協(xié)同重排域的位形熵，Δμ為阻止單體（鏈段或分子）協(xié)同重排的能量勢(shì)壘，是一個(gè)溫度無(wú)關(guān)的量，k表示Boltzmann常數(shù)；R為氣體常數(shù)；T2為Adam－Gibbs溫度，表示位形熵消失時(shí)的溫度。

文獻(xiàn)［3］將位形熵Sc作為結(jié)構(gòu)參數(shù)取代Tf，提出位形熵模型，其松弛時(shí)間表達(dá)式為：文獻(xiàn)［4］引入松弛極限態(tài)假設(shè)，即GR模型。

文獻(xiàn)［6］在玻璃化轉(zhuǎn)變過(guò)程引入極限情況（即H模型），將玻璃化轉(zhuǎn)變過(guò)程假設(shè)成2個(gè)相互獨(dú)立的過(guò)程：①以玻璃化轉(zhuǎn)變溫度為T(mén)f′時(shí)的極慢降溫過(guò)程（xs＝0）；②以玻璃化轉(zhuǎn)變溫度為T(mén)f″時(shí)的極快降溫過(guò)程（xs＝1）。通過(guò)引入非線性參數(shù)xs（0＜xs＜1），將Sc表示成依賴于熱力學(xué)溫度T和極限假想溫度Tf（Tf′與Tf″范圍內(nèi)某一值）2部分，即

2 實(shí)驗(yàn)部分

2.1 實(shí)驗(yàn)儀器及試劑

采用配備G－M二級(jí)制冷系統(tǒng)的DSC儀器（Pyris－Diamond，美國(guó) Perkin－Elmer公司），標(biāo)稱(chēng)最低溫度193K，實(shí)測(cè)DSC內(nèi)樣品溫度最低可達(dá)約203K（隨室溫變化1K左右）。DSC以金屬銦（99.999%，Perkin－Elmer）和純水（實(shí)驗(yàn)室自制二次蒸餾去離子水）做兩點(diǎn)溫度標(biāo)定，標(biāo)定范圍為210～320K；以冰的融化潛熱（333.88kJ／kg）進(jìn)行單點(diǎn)焓尺度標(biāo)定，標(biāo)定速率10K／min。氮?dú)猓?9.999%）氣氛，流量約20mL／min并保持恒定。分析純結(jié)晶態(tài)木糖醇（Sigma，上海君創(chuàng)生物科技有限公司進(jìn)口分裝，99%）未經(jīng)進(jìn)一步純化直接使用。取樣品8～10mg（精確到±0.01mg），置于標(biāo)準(zhǔn)液體鋁皿中央并壓蓋密封。

2.2 實(shí)驗(yàn)過(guò)程

先將結(jié)晶態(tài)樣品從室溫以10K／min升溫至410K充分熔融，然后降溫至Tg（247K左右）以下。將玻璃態(tài)樣品升溫至加樣溫度（室溫）后，以0.5K／min的降溫速率降低至215K，待熱流穩(wěn)定后按如下溫度程序采集數(shù)據(jù)：① 等溫1min；② 以10K／min升溫到310K；③ 等溫1min。將降溫速率由0.5K／min 分別改為1、2、5、10、20K／min，重復(fù)以上實(shí)驗(yàn)。取與樣品皿質(zhì)量相同的空皿（±0.01mg），驟冷至215K，熱流穩(wěn)定后，按與分析樣品相同的溫度程序獲得基線熱流，用于比定壓熱容計(jì)算。為了保證熱接觸狀況的一致性，針對(duì)樣品的所有分析結(jié)束后再獲得一條用于計(jì)算比定壓熱容的基線。

2.3 曲線擬合算法

本文對(duì)3個(gè)模型分別利用曲線擬合技術(shù)，采用曲線擬合程序以Matlab（附優(yōu)化工具箱）語(yǔ)言編制。由于模型中各個(gè)參數(shù)之間存在很強(qiáng)的關(guān)聯(lián)性，在對(duì)AGV模型進(jìn)行曲線擬合時(shí)，不能對(duì)模型中4個(gè)參數(shù)同時(shí)進(jìn)行擬合，故采用嵌套循環(huán)搜索，將T2置于外循環(huán)固定，搜索其他3個(gè)參數(shù)（A，B，β）的最優(yōu)值。對(duì)GR和H模型采取同樣的方法，將T2置于外循環(huán)分別搜索A，B，β，δ和A，B，β，xs。

現(xiàn)象學(xué)模型中參數(shù)值依賴于DSC掃描時(shí)的熱歷史，玻璃體系經(jīng)歷不同的熱歷史，模型的各參數(shù)值可能會(huì)有所不同［2］。但各現(xiàn)象學(xué)模型提出的初衷是尋求材料常數(shù)，即和玻璃體系所經(jīng)歷的熱歷史無(wú)關(guān)［1］，模型參數(shù)擬合值與實(shí)驗(yàn)數(shù)據(jù)的一致程度是本文討論模型預(yù)測(cè)力大小的重要標(biāo)準(zhǔn)之一。

文獻(xiàn)［1，3］中用現(xiàn)象學(xué)模型對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合時(shí)，只考慮單純的評(píng)估模型參數(shù)與玻璃體系所經(jīng)歷的熱歷史之間的關(guān)系，沒(méi)有考慮模型與實(shí)驗(yàn)數(shù)據(jù)之間擬合的均方差，為了更加精確地比較3種模型對(duì)實(shí)驗(yàn)數(shù)據(jù)的預(yù)測(cè)力，本文針對(duì)6次實(shí)驗(yàn)數(shù)據(jù)利用Nedler－Mead搜索方法尋求最優(yōu)參數(shù)，使模型與實(shí)驗(yàn)數(shù)據(jù)之間擬合的均方差σa最小，即

其中，i為實(shí)驗(yàn)次數(shù)；j為實(shí)驗(yàn)掃描時(shí)的溫度點(diǎn)；w（i）為由峰高確定的權(quán)重因子（0＜w（i）＜1）。

3 結(jié)果與討論

3.1 模型參數(shù)

文獻(xiàn)［1－9］表明，利用現(xiàn)象學(xué)模型對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行非線性曲線擬合獲得模型參數(shù)時(shí)，模型參數(shù)的初值選取尤其重要，不同的初始值可能會(huì)導(dǎo)致不同的搜索結(jié)果。

本文在對(duì)木糖醇實(shí)驗(yàn)數(shù)據(jù)擬合之前，首先取升降溫速率為10K／min，根據(jù)松弛時(shí)間τ≈100s時(shí)的玻璃化轉(zhuǎn)變情況來(lái)估算指前因子A，采用文獻(xiàn)［1］中DSC法測(cè)定的木糖醇玻璃化轉(zhuǎn)變溫度Tg（取247.5K），AGV與TNM模型參數(shù)之間的關(guān)系式［1］為：

從文獻(xiàn)［2］中獲得（6）式中非線性參數(shù)x值，可估算AGV模型參數(shù)T2≈149K。Williams等采用黏彈性實(shí)驗(yàn)，結(jié)合 WLF方程分析得出大多數(shù)玻璃化物質(zhì)Tg－T2≈50K。結(jié)合這2種方法，設(shè)定可變參數(shù) T2為127.5 ～197.5K，搜索AGV模型中其他3個(gè)參數(shù)的值，見(jiàn)表1所列。

表1 設(shè)定不同的T2值A(chǔ)GV模型對(duì)木糖醇實(shí)驗(yàn)的擬合參數(shù)

由表1可以看出，隨著T2的增加，指前因子A呈增加趨勢(shì)，而非指數(shù)性質(zhì)β與B均呈現(xiàn)出總體減小的趨勢(shì)，和已報(bào)道的聚合物［2，6］變化趨勢(shì)相同，但擬合的非指數(shù)性質(zhì)β值較TNM模型［1］擬合值小。木糖醇在不同降溫速率下實(shí)驗(yàn)比定壓熱容（散點(diǎn)）與GR模型擬合（實(shí)線）最好時(shí)的結(jié)果對(duì)照（T2＝137.5K）如圖1所示，模型各參數(shù)擬合結(jié)果見(jiàn)表2所列。

圖1 木糖醇實(shí)驗(yàn)cp（T）與GR模型擬合結(jié)果

表2 設(shè)定不同的T2值GR模型對(duì)木糖醇實(shí)驗(yàn)的擬合參數(shù)

從圖1可以看出，GR模型可以很好地重現(xiàn)6種不同降溫速率的DSC實(shí)驗(yàn)數(shù)據(jù)。

由表2可以看出，隨著T2取值的增加，GR模型加入新參數(shù)δ的值并沒(méi)有呈現(xiàn)規(guī)律性變化，而是在一定范圍內(nèi)（0.1～0.44J／（g·K））波動(dòng)，是木糖醇玻璃化轉(zhuǎn)變過(guò)程中比定壓熱容增量（Δcp（Tg）≈0.96J／（g·K））［1］的10%～45%。δ在木糖醇中所表現(xiàn)的性質(zhì)與已報(bào)道聚合物結(jié)果有一定的差異，文獻(xiàn)［3－4］發(fā)現(xiàn)一些聚合物δ的值是其玻璃化轉(zhuǎn)變過(guò)程中比熱容增量的45%左右，說(shuō)明木糖醇玻璃體系存在的假設(shè)亞穩(wěn)極限態(tài)比聚合物更加趨近它們各自對(duì)應(yīng)溫度下的平衡態(tài)。木糖醇在不同降溫速率下實(shí)驗(yàn)數(shù)據(jù)無(wú)量綱比定壓熱容（散點(diǎn)）與H模型擬合（實(shí)線）最好時(shí)的結(jié)果對(duì)照（T2＝157.5K）如圖2所示，模型各參數(shù)擬合值見(jiàn)表3所列。

圖2 木糖醇實(shí)驗(yàn)cp（T）與H模型擬合結(jié)果

表3 設(shè)定不同的T2值H模型對(duì)木糖醇實(shí)驗(yàn)的擬合參數(shù)

H模型對(duì)木糖醇實(shí)驗(yàn)數(shù)據(jù)擬合結(jié)果顯示，隨著T2值的增加，新參數(shù)xs值增加明顯，和聚甲基丙烯酸甲酯［7］xs值表現(xiàn)趨勢(shì)相同，但木糖醇玻璃體系的xs值（0.035＜xs＜0.18）變化幅度遠(yuǎn)大于PAMM中xs值（0.02＜xs＜0.08）。同時(shí)，由TNM模型參數(shù)［1］和H模型參數(shù)的近似關(guān)系式［6］獲得xs值為0.03～0.34。

3個(gè)模型同時(shí)對(duì)木糖醇實(shí)驗(yàn)數(shù)據(jù)進(jìn)行曲線擬合，對(duì)比表1～表3發(fā)現(xiàn)，隨著T2值的增加，3個(gè)模型得出的指前因子A均呈增加趨勢(shì)，而非指數(shù)性質(zhì)β與B均呈現(xiàn)出總體減小的趨勢(shì)。但GR模型得出的指前因子A值較其他2個(gè)模型變化幅度要小，且由GR和H模型擬合的非指數(shù)性質(zhì)β明顯高于AGV模型擬合結(jié)果，同時(shí)也高于TNM模型［1］擬合的β值，這表明GR和H模型對(duì)位形熵模型改進(jìn)時(shí)，在模型中引入新參數(shù)對(duì)其他模型參數(shù)有一定的影響。

3.2 均方差σa

由表1～表3可看出，AGV和GR 2種模型擬合的T2值一致表現(xiàn)在137.5K時(shí)（Tg－T2≈110K）σa最小，而H模型擬合的T2值在157.5K時(shí)（Tg－T2≈90K）σa達(dá)到最小，表明模型擬合與實(shí)驗(yàn)數(shù)據(jù)的一致性最好，但這不切實(shí)際，與 WLF方程分析得出大多數(shù)玻璃化物質(zhì)Tg－T2≈50K相差甚遠(yuǎn)，類(lèi)似的情況在聚合物［1，7－8］的研究中也有顯現(xiàn)。

對(duì)AGV模型而言，隨著T2取值的繼續(xù)增大，σa呈現(xiàn)明顯增大的趨勢(shì)，對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合效果表現(xiàn)不理想。GR和H模型表現(xiàn)出相同的變化趨勢(shì)，但這2種模型對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合的σa與AGV模型相比總體上有所減小，在T2取值較小時(shí)，GR模型對(duì)AGV模型改善較大，隨著T2的增大，H模型對(duì)AGV改善效果較GR模型明顯。

同時(shí)可以看出，3種模型在對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合結(jié)果最好時(shí)，σa值變化不大，3種模型對(duì)實(shí)驗(yàn)數(shù)據(jù)的擬合效果相當(dāng)。如以σa最小時(shí)的擬合結(jié)果作為模型參數(shù)的預(yù)測(cè)值，H模型預(yù)測(cè)的T2（T2＝157.5K）明顯高于前2種模型對(duì)T2值（T2＝137.5K）的預(yù)測(cè)，且隨著T2的繼續(xù)增大，H模型對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合的σa仍然較前2種模型表現(xiàn)理想。從模型與實(shí)驗(yàn)數(shù)據(jù)擬合的σa總體改善情況討論模型預(yù)測(cè)力，H模型對(duì)實(shí)驗(yàn)數(shù)據(jù)的預(yù)測(cè)力較AGV模型和GR模型有了一定提高。

由表3可以發(fā)現(xiàn)，T2值從157.5K擬合效果最好結(jié)果開(kāi)始繼續(xù)增大到187.5K，σa仍然表現(xiàn)理想，以此值作為H模型對(duì)T2的預(yù)測(cè)，可以發(fā)現(xiàn)這一結(jié)果和介電松弛實(shí)驗(yàn)通過(guò)Vogel－Tammann－Fulcher（VTF）方程計(jì)算得到的木糖醇T2值［9］（190±5）K結(jié)果一致。

結(jié)合介電法測(cè)得的結(jié)果，選擇T2＝187.5K，由表1～表3中的σa可知，H模型和GR模型對(duì)實(shí)驗(yàn)數(shù)據(jù)的預(yù)測(cè)力較AGV模型分別提高了約36%和80%。

3.3 Adam－Gibbs溫度T2

通過(guò)現(xiàn)象學(xué)模型對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行曲線擬合，擬合出不同溫度下焓松弛時(shí)間，對(duì)合理選擇貯藏條件及延長(zhǎng)貨架壽命有一定的理論指導(dǎo)價(jià)值。在比較3種現(xiàn)象學(xué)模型對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合情況時(shí)，有必要考慮不同現(xiàn)象學(xué)模型在選擇相同T2值時(shí)，它們對(duì)松弛時(shí)間計(jì)算結(jié)果有無(wú)影響，3種現(xiàn)象學(xué)模型在T2＝197.5K時(shí)，Tg附近焓松弛時(shí)間計(jì)算結(jié)果如圖3所示。

由圖3可以看出，在靠近玻璃化轉(zhuǎn)變溫度Tg時(shí)，3種模型對(duì)松弛時(shí)間的擬合結(jié)果比較一致，隨著T相對(duì)于Tg偏離距離的增大，3種模型計(jì)算出的松弛時(shí)間出現(xiàn)一定差異，在T為225K時(shí)，H模型計(jì)算出的松弛時(shí)間是AGV模型和GR模型計(jì)算出松弛時(shí)間的3倍。另外，在計(jì)算過(guò)程中發(fā)現(xiàn)，對(duì)于同一種模型而言，T2選取在一定范圍內(nèi)計(jì)算出體系的松弛時(shí)間表現(xiàn)出很好的一致性（其他T2值計(jì)算結(jié)果和圖示幾乎重合，為清晰起見(jiàn)未在圖中標(biāo)出），和苯乙烯－丙烯腈共聚物表現(xiàn)出的結(jié)果類(lèi)似［5］。所以，在木糖醇玻璃體系中，這3種模型對(duì)T2值預(yù)測(cè)的準(zhǔn)確性關(guān)系到模型對(duì)實(shí)驗(yàn)數(shù)據(jù)預(yù)測(cè)能力。

圖3 3種模型計(jì)算木糖醇的松弛時(shí)間

AGV、GR和H模型均為基于AG理論的位形熵焓松弛現(xiàn)象學(xué)模型，GR和H模型又同時(shí)是對(duì)位形熵模型的非線性部分的改進(jìn)。已有文獻(xiàn)［6］表明一些模型對(duì)應(yīng)參數(shù)之間具有一定的近似關(guān)系，由（6）式發(fā)現(xiàn)，T2值的大小和現(xiàn)象學(xué)模型的非線性程度有關(guān)，即非線性程度越大，T2值反而越小，如以擬合σa最小時(shí)結(jié)果作為模型參數(shù)，則發(fā)現(xiàn)3種模型預(yù)測(cè)的T2普遍小于由其他實(shí)驗(yàn)方法（介電松弛法、WLF方程分析等）分析得到的T2，AGV和GR模型擬合的Tg－T2的差值甚至高達(dá)110K左右，文獻(xiàn)［6－8］在利用AG理論模型對(duì)聚合物玻璃體系焓松弛分析時(shí)也發(fā)現(xiàn)和木糖醇相同的結(jié)果。

采用位形熵模型對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合過(guò)程中，Sc不完全取決于影響非線性性質(zhì)的結(jié)構(gòu)參數(shù)Tf［9］，為了補(bǔ)償位形熵理論中的非線性，H模型改進(jìn)過(guò)程中添加更多的表征非線性特征的參數(shù)，對(duì)木糖醇玻璃系而言，將Tg－T2的差值降至90K，對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合有了一定的提高。

用現(xiàn)象學(xué)模型對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行曲線擬合時(shí)，比定壓熱容增量Δcp（T）的處理方式對(duì)T2的結(jié)果也有一定的影響，Δcp（T）取常數(shù)或曲線表達(dá)形式時(shí)，基于AG理論的現(xiàn)象學(xué)模型松弛時(shí)間趨近于無(wú)窮大得到的T2值和Kauzmann溫度（TK）相等［6］。文獻(xiàn)［10］對(duì)T2和TK之間的關(guān)系也作了研究，并通過(guò)介電數(shù)據(jù)和比熱頻譜TK和T2之間的關(guān)系，發(fā)現(xiàn)在小分子有機(jī)物T2≈TK－20K。以擬合σa最小時(shí)結(jié)果作為模型的參數(shù)，在木糖醇中運(yùn)用Dixon方法，可將H模型預(yù)測(cè)的T2值升至177.5K，與其他實(shí)驗(yàn)方法分析得到的結(jié)果接近。

對(duì)木糖醇而言，H模型引入xs將物質(zhì)的玻璃化過(guò)程設(shè)想成2種極限過(guò)程，從熱力學(xué)方面看有些不切實(shí)際。但物質(zhì)在降溫過(guò)程中，玻璃體系在玻璃化轉(zhuǎn)變溫度Tg附近會(huì)發(fā)生一個(gè)純動(dòng)力學(xué)過(guò)程，文獻(xiàn)［11］發(fā)現(xiàn)在玻璃形成體系中有動(dòng)力學(xué)非均質(zhì)特性存在。另外，Glotzer通過(guò)數(shù)值分析研究發(fā)現(xiàn)在Lennard－Jones系統(tǒng)中快速域和緩慢域是可以共存的。將木糖醇玻璃體系焓松弛過(guò)程描述成快速和緩慢2個(gè)過(guò)程的綜合，可使模型對(duì)實(shí)驗(yàn)數(shù)據(jù)的預(yù)測(cè)有一定的提高，表明H模型對(duì)木糖醇體系的焓松弛行為的描述更加成功。

4 結(jié) 論

運(yùn)用AGV、GR和H 3種模型，結(jié)合DSC實(shí)驗(yàn)考察木糖醇玻璃體系的焓松弛行為，進(jìn)行了模型預(yù)測(cè)力的分析。結(jié)果表明：

（1）3種模型均能很好地重現(xiàn)木糖醇經(jīng)歷不同熱歷史玻璃體系的升溫cp曲線。

（2）3種模型預(yù)測(cè)的T2值普遍小于其他實(shí)驗(yàn)方法分析（如 WLF方程分析、介電松弛等）得到的結(jié)果。

（3）從現(xiàn)象學(xué)模型對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合的均方差和模型對(duì)T2預(yù)測(cè)能力分析發(fā)現(xiàn)，H模型對(duì)實(shí)驗(yàn)數(shù)據(jù)的預(yù)測(cè)力優(yōu)于AGV和GR模型。

［1］高 才，楊 鎖，劉向農(nóng)，等.木糖醇玻璃焓松弛的量熱分析［J］.物理化學(xué)學(xué)報(bào)，2010，26（1）：7－12.

［2］Andreozzi L，F(xiàn)aetti M，Giordano M，et al.Enthalpy relaxation in polymers：a comparison among different multiparameter approaches extending the TNM／AGV model［J］.Mac－romolecules，2003，36：7379－7387.

［3］Ribelles G J L，Pradas M M.Structural relaxation of glassforming polymers based on an equation for configurational entropy：1，DSC experiments on polycarbonate［J］.Macromolecules，1995，28：5867－5877.

［4］Ribelles G J L，Pradas M M，Vidaurre Garayo A，et al.Structural relaxation of glass－forming polymers based on an equation for configurational entropy：3，on attained at infinite time in the relaxation process results on imide［J］.Polymer，1997，38：963－969.

［5］Due?nas J M M.Structural relaxation of glass－forming polymers based on an equation for configurational entropy：4，structural relaxation in styrene－acrylonitrile copolymer［J］.Polymer，1997，35：2201－2217.

［6］Hutchinson J M.Application of the Adam－Gibbs equation to the non－equilibrium glassy state ［J］.Macromolecules，2000，33：5252－5262.

［7］Andreozzi L，F(xiàn)aetti M，Zulli F，et al.Enthalpy relaxation of polymers：comparing the predictive power of two configura－tional entropy models extending the AGV approach［J］.Eur Phys J B，2004，41：383－393.

［8］Hutchinson J M，Montserrat S，Calventus Y，et al.On the application of the Adam－Gibbs equation to the non－equilibrium glassy state ［J］.Non－Cryst Solids，2002，307（310）：412－416.

［9］Wehn R，Lunkenheimer P，Loidl A.Broadband dielectric spectroscopy and aging of glass formers［J］.Non－Cryst Solids，2007，353：3862－3870.

［10］Dixon P.Specific－h(huán)eat spectroscopy and dielectric susceptibility measurements of salol at the glass transition ［J］.Phys Rev B，1990，42：8179－8186.

［11］Ferguson A，Chakraborty B.Spatially heterogenous dynamics in dense，driven granular flows［J］.EPL，2007，78（2）：28003.doi：10.1209／0295－5075／78／28003.