鄭成志，趙曙輝，楊曉輝

(1.東華大學 化學化工與生物工程學院，上海201620; 2.上海諾科生物科技有限公司，上海 201612)

蒙脫土納米復合材料作為成碳劑的熱氧化分解動力學研究

鄭成志1，趙曙輝1，楊曉輝2

(1.東華大學 化學化工與生物工程學院，上海201620; 2.上海諾科生物科技有限公司，上海 201612)

采用蒙脫土納米復合材料與高聚磷酸銨(APP)組成的阻燃體系對棉織物進行阻燃整理，并用非等溫熱重分析方法對復合阻燃體系整理后棉織物的熱氧化分解特性進行研究.在空氣氣氛下，以不同的升溫速率將樣品從室溫加熱至900 ℃.通過基于動力學補償效應的恒定動力學參數(IKP)法計算整理后織物的非等溫氧化分解動力學參數，并由IKP法求得常用固相機理函數概率分布，確定復合體系的動力學模型.最終確定整理后織物的熱氧化分解過程可由成核和核增長以及冪級反應機理函數進行描述.

蒙脫土納米復合材料；成碳劑；熱氧化分解動力學；反應機理

膨脹阻燃體系主要由酸源、碳源和氣源組成，在高溫條件下，該體系能夠在基體表面形成膨脹的碳層，減緩氣相與凝固相之間的熱量及質量轉移，阻止基體繼續(xù)燃燒，達到提高基體阻燃性能的目的.基體表面膨脹型碳層的產生主要通過以下反應過程實現：酸源高溫分解產生無機酸，無機酸與多元醇發(fā)生反式羰基化反應促進碳層的形成，氣源分解產生的不燃性氣體使碳層進一步膨脹，從而實現阻燃[1].

蒙脫土納米復合材料可以代替多元醇類作為膨脹體系的成碳劑.文獻[2]使用尼龍6(PA6)及PA6/蒙脫土納米復合材料代替季戊四醇制備阻燃膨脹體系，克服了多元醇類在聚合物基體中可遷移及相容性較差的缺點.文獻[3]用PA6/蒙脫土納米復合材料與高聚磷酸銨(APP)混合制備膨脹阻燃體系，證實了納米分散的蒙脫土能夠加速碳層的產生，增加碳層的密度，實現聚合物基體較好的阻燃性能.文獻[4] 用聚氨酯/蒙脫土納米復合材料與三聚氰胺聚磷酸鹽(MPP)制備的膨脹阻燃體系，能夠較好地抑制煙霧以及毒性氣體的釋放，該體系中MPP分解產生的聚磷酸能夠促進聚氨酯交聯成碳，并且聚磷酸與納米分散的蒙脫土形成在高溫下穩(wěn)定存在的磷酸鋁的陶瓷結構，阻止聚合物基體燃燒.

聚合物材料的燃燒主要為氣相與凝固相之間的熱及質量的轉移過程.其中發(fā)生的燃燒降解機理包含以下動力學模型：相界面反應模型、氣體擴散模型、化學反應級數模型、冪級反應級數模型以及成核與生長模型[5].準確應用相應的模型解釋聚合物熱降解過程是非常重要的.然而從非等溫的熱氧化分解模型中確定熱降解動力學參數比較困難，并且計算量很大.為了克服上述缺陷，本文采用由文獻[6]首次提出并由文獻[7]改進的恒定動力學參數(invariant kinetic parameter,IKP)法計算阻燃體系的動力學參數.由于在使用動力學方程計算時，不需要假設相應的前提條件，因此，計算結果可嚴格地反映體系的動力學性能.同時該方法能夠計算體系中所應用的18種降解動力學方程的概率分布[8],從而推斷出整個體系的分解模型.

本文采用合成的硬脂酸酰胺/蒙脫土納米復合材料(EBS/MMT)與APP組成膨脹阻燃體系對純棉織物進行阻燃整理，并對整理后的織物應用恒定動力學參數法研究該體系的分解模型與分解機理.該方法有助于更好地理解整理后織物的熱氧化分解行為，并對防止火災的發(fā)生具有理論指導意義.

1　試驗部分

1.1試劑和織物

織物采用全棉府綢漂白半制品，經緯密度為510根/10 cm×280根/10 cm,面密度為120 g/m2，紗線線密度為25 tex×25 tex.

試劑：二乙烯三胺、冰醋酸，分析純，國藥集團化學試劑有限公司；硬脂酸，分析純，泰科棕化有限公司；鈉基蒙脫土，上海試四赫維化工有限公司；APP，鎮(zhèn)江星星阻燃劑有限公司；乳化劑(T05)，上海諾科生物科技有限公司；改性2D樹脂交聯劑，上海纖化生物科技有限公司.

1.2試驗儀器

德國NRTZSCH公司的NRTZSCH TG 209F1型熱重分析儀，MATHIS公司的Mathis烘蒸兩用機，英國ROACHES公司的EHP 350MM 型微調水平軋車.

1.3硬脂酸酰胺蒙脫土納米復合材料的制備

取0.17 mol的二乙烯三胺與0.17 mol冰醋酸，以及5.6 g的蒙脫土于燒杯中，放置在磁力攪拌器上攪拌分散3 h，將0.34 mol的硬脂酸放入帶有分水器及N2保護裝置的圓底四口燒瓶中，升溫至100 ℃，保溫30 min.待硬脂酸完全融化后，將二乙烯三胺與蒙脫土的分散液滴入燒瓶中，升溫至160～180 ℃反應5 h，冷卻后即得硬脂酸酰胺蒙脫土納米復合材料.

1.4阻燃整理液的制備

稱取定量的復合材料于燒瓶中，加熱融化，滴加冰醋酸和乳化劑，加入去離子水，用高速分散機分散30 min，制得復合材料水分散液.在分散液中加入APP和樹脂交聯劑制備出阻燃整理液APP/EBS/MMT體系，其配方如表1所示.

表1　阻燃整理工藝配方Table 1　The process formula of flame retardant finishing

1.5棉織物的阻燃整理

采用浸軋法，浴比為1∶20,二浸二軋(軋液率為100%)，預烘(80 ℃，5 min),焙烘(150 ℃，3 min).

1.6熱氧化分解動力學分析

分別采用4種不同的升溫速率 (βv=10,20,30,40 ℃/min) 對整理后的棉織物進行熱氧化分解動力學測試分析.測試條件：空氣和N2氛圍，流速為25 mL/min.通過熱重(TG)曲線計算恒定動力學參數以及熱分解方程，研究阻燃體系的熱分解行為.

1.7恒定動力學參數法

在多數情況下，特定升溫速率下的TG試驗曲線可被多種機理函數描述，不同機理函數得到的動力學參數往往差別很大[7].為了減少試驗條件對動力學參數的影響，本文采用恒定動力學參數法.

恒定動力學參數主要是根據降解速率(dα/dt)方程推導的，其表達式為

(1)

式中：α為轉換程度；t為時間；k為速率常數；f(α)為機理函數.

通過Coats-Redfern方法[9]可以計算出表觀活化能及表觀前驅因子.恒定動力學參數法主要依據為補償效應，即對于每一個熱分解方程fj(α),在相應的升溫速率βv下，用logAiv對Eiv作圖，如果函數關系為一條直線，則認為具有補償效應.具體表達式為

logAiv=Bv+lvEiv

(2)

其中：Aiv和Eiv分別為表觀前驅因子和表觀活化能，下標i為機理函數個數，i=1,2,…，18,v為升溫速率，在升溫速率βv下，通過Coats-Redfern方法可計算Aiv與Eiv；Bv和lv為變量函數，其表達式為

Bv=logkv

(3)

lv=(2.3RTv)-1

(4)

其中：kv為體系在溫度為Tv條件下的速率常數；R為氣體常數.Bv與lv隨著試驗條件的變化而變化.用logkv對1/Tv作圖，可以計算得到式(5)的截距及斜率，從而計算得出改變升溫速率的動力學參數Ainv以及Einv.

logkv=logAinv-Einv/2.3RTv

(5)

一旦動力學補償效應成立，就可以計算得到Ainv和Einv,并根據文獻[10]的統(tǒng)計學處理方法，得到每一個分解方程的概率分布.

2　結果與討論

2.1熱氧化分解分析

經過APP/EBS/MMT體系整理后的織物，在空氣與N2氛圍下，升溫速率為10 ℃/min時的TG曲線如圖1所示.

圖1　空氣和N2氛圍下APP/EBS/MMT整理后織物的TG曲線圖Fig.1　TG curves of the APP/EBS/MMT treated fabrics at air and nitrogen condition

由圖1可知，在空氣氛圍下，APP/EBS/MMT膨脹體系的熱氧化分解過程分為3個階段：第一階段主要發(fā)生在溫度為240 ℃左右，該階段聚合物開始分解；第二階段在溫度為260～290 ℃，為膨脹體系的分解成碳過程；第三階段在溫度為290～780 ℃，該階段氧氣參與殘留物的氧化分解過程，分解速率變小，最終殘留物質量分數為5.4%.在N2氛圍下，體系的熱氧化分解過程也可以分為3個階段，其第一、二階段的降解和空氣氛圍下基本一致，但在300 ℃以上，膨脹體系分解完成，分解速率緩慢，說明在此過程中織物表面已形成了一層致密的碳層，阻止織物的進一步降解.

空氣中的O2在膨脹阻燃體系的分解過程中起到了非常重要的作用.由于O2的存在，APP/EBS/MMT膨脹體系在300 ℃以上并未形成穩(wěn)定的碳層，而是進一步與O2作用，發(fā)生氧化分解[10]，使棉織物的殘余質量分數由N2氛圍下的25% 降至5%.

N2氛圍下APP/EBS/MMT體系在不同升溫速率下的TG曲線如圖2所示.空氣氛圍下APP/EBS/MMT體系在不同升溫速率的TG和DTG曲線如圖3所示.由圖2和3可以看出，N2氛圍下，當升溫速率為30和40 ℃/min時，升溫速率對TG曲線的影響不大，第二階段熱分解曲線基本重合，而在空氣氛圍下整個體系的初始熱分解溫度隨著升溫速率的升高而增大.

圖2　氮氣氛圍不同升溫速率下APP/EBS/MMT體系的TG曲線Fig.2　TG curves of APP/EBS/MMT system under different heating rates at nitrogen condition

(a) TG曲線

(b) DTG曲線圖3　空氣氛圍不同升溫速率下APP/EBS/MMT體系的TG與DTG曲線Fig.3　TG and DTG curves of APP/EBS/MMT system under different heating rates at air condition

本文中恒定動力學參數法主要應用于熱氧化分解的第二階段.事實上，在N2氛圍下，第一階段并不能觀察到線性補償效應，即恒定動力學參數法并不能進行第一階段的研究[5].而在空氣氛圍下，用恒定動力學參數進行第一階段研究，線性擬合度較差，誤差較大[11].通過對熱分解數據進行處理發(fā)現,APP/EBS/MMT膨脹體系的TG曲線在第二分解階段較好地符合方程式(5)，存在線性關系.因此，本文對APP/EBS/MMT在空氣氛圍下第二階段(260～290 ℃)熱氧化分解動力學的分解機理和分解模型進行研究.

2.2恒定動力學參數法的計算與分析

對于每一個降解方程fj(α),在相應的升溫速率βv下，用logAj對Ej作圖，如果函數關系為一條直線，則認為具有補償效應，可以用IKP法進行熱氧化分解分析.

基于IKP法,共選取了18種固相反應機理函數(S1～S18)[7]用于動力學參數計算.通過Coats-Redfern積分法[9]計算不同升溫速率下的18個表觀活化能及前驅因子，結果如表2所示.將表2中的數值代入式(2)，可驗證E與logA之間的動力學補償效應,其中l(wèi)ogA=loge×lnA.

表2　不同機理函數和升溫速率下的動力學參數Table 2　The kinetic parameters of different mechanism functions and the heating rates

圖4給出了4種升溫速率下APP/EBS/MMT整理后織物的熱氧化分解動力學參數的線性補償關系.由圖4可知，各直線基本重合，表明APP/EBS/MMT非等溫熱氧化分解可以用IKP方法進行模擬分析.從每一條直線的斜率與截距可以得到每種升溫速率下的動力學補償參數Bv與lv.

圖4　不同升溫速率下APP/EBS/MMT整理后織物的熱氧化分解參數的線性補償關系Fig.4　A liner relation between the thermal oxidation decomposition parameters of fabrics treated with APP/EBS/MMT at different heating rates

本文所列機理函數在4種升溫速率下的動力學補償參數計算結果如表3所示.根據式(5)及Bv與lv的線性關系(如圖5)，由斜率和截距可以求得恒定動力學參數：Einv=137.4 kJ/mol,ln(Ainv/(m3·mol-1)n-1·s-1))=28.3.

表3　不同升溫速率的動力學補償參數Table 3　Compensation parameters at different heating rates

圖5　不同升溫速率下動力學補償參數的線性關系Fig.5　A liner relation between the kinetic compensation parameters at different heating rates

表4給出了棉織物恒定動力學參數[11]及經過APP/EBS/MMT處理后棉織物的恒定動力學參數.由表4可知，經膨脹阻燃體系處理后，棉織物的主要降解溫度由320～380 ℃下降至260～290 ℃，其活化能較純棉織物減少近一半.即織物經過APP/EBS/MMT整理后，阻燃體系在較低的溫度下就能相互作用，在織物表面形成一層致密的碳層，從而實現較好的阻燃性能和熱穩(wěn)定性.

表4　熱氧化分解的恒定動力學參數Table 4　Invariant kinetic parameters of thermo-oxidative degradation

文獻[7]中將聚合物的降解機理分為5個動力學模型：成核與核增長模型，包含機理函數S1～S4；相界面反應模型，包含機理函數S6～S8；擴散模型，包含機理函數S9～S11及S18；冪級反應模型，包含機理函數S12～S14及S17；化學反應級數模型，包含機理函數S5及S15～S16.

通過文獻[7]給出的降解動力學機理函數S1～S18，根據文獻[10]計算出每個機理函數的概率，結果如圖6所示，則可以得到每一種動力學模型的概率分布(機理函數概率之和)，相應的結果總結于表5中.由圖6可知，復合體系降解函數的概率分布主要集中在函數S1～ S4,S13和S14中.因此，APP/EBS/MMT的熱氧化分解主要階段可以用成核與核增長模型(S1～S4)以及冪級反應模型(S13和S14)進行模擬分析.

圖6　IKP法計算得到的機理函數的概率分布Fig.6　Probabilities distribution of the mechanism function evaluated by IKP

表5　每一種動力學模型的概率分布Table 5　Probabilities distribution associated with the kinetic models　%

APP/EBS/MMT膨脹體系的熱氧化降解并不能用單一的降解函數進行表征，而是通過一系列函數進行組合模擬表征.APP/EBS/MMT阻燃膨脹體系最有可能的降解方程f(α)可以用式(6)進行表征.

(6)

將該降解方程的f(α)對α作圖，可以得到圖7.由圖7可知，熱氧化降解方程的值隨著α值的增加而緩慢增大，至α為0.4時達到最大.這可能是因為APP/EBS/MMT不同組分之間相互作用(協同作用)，使得降解增加，形成碳化物，而隨著轉化程度的繼續(xù)增加，函數值減小，降解變得很緩慢，說明此時體系外圍已經形成了一層致密的保護層，阻止了體系進一步氧化降解.

圖7　空氣氛圍下復合體系降解函數曲線圖Fig.7　The curve of degradation function of hybrid system at air condition

由前文分析可知，經APP/EBS/MMT復合體系整理后織物的熱氧化分解發(fā)生了明顯的變化，主要分解溫度由未整理的320～380 ℃下降至260～290 ℃，且在260～290 ℃空氣與氮氣氛圍下的TG曲線重合(圖2)，說明在該階段空氣中的氧氣并沒有參與復合體系的氧化反應，并且由于APP的開始分解溫度在300 ℃左右[3],所以該階段發(fā)生的分解主要為復合阻燃體系中APP與EBS/MMT間的化學反應以及EBS/MMT自身的分解反應，分解機理為成核與核增長及冪級反應模型，這與文獻[5]結論相符.另外，升溫速率的改變不影響分解機理及發(fā)生的化學反應，僅影響復合體系的初始分解溫度與最大失重速率對應的溫度.

3　結　語

(1) APP/EBS/MMT整理后織物的熱重分析結果表明，整理前后織物的主要降解溫度范圍發(fā)生變化.空氣中的O2不影響第二階段(260～290 ℃)的降解過程，但對在第二階段形成的不飽和碳層有很明顯的氧化分解作用.

(2) 經APP/EBS/MMT整理的織物，其主要降解階段(260～290 ℃) 可以用IKP方法計算熱氧化降解動力學參數，得到Einv=137.4 kJ/mol,ln (Ainv/(m3·mol-1)n-1·s-1))=28.3，其參數約為純棉織物的一半，說明經過APP/EBS/MMT整理后，阻燃體系在較低的溫度下就能相互作用，在織物表面形成一層致密的碳層，從而實現較好的阻燃性能和熱穩(wěn)定性.

(3) 經APP/EBS/MMT整理后，織物的主要降解模型為成核與核增長及冪級反應模型.

[1] ANTONOV A,RESHETNIKOV I S,KHALTURINSKIJ N A.Combustion of char forming polymeric systems[J].Russ Chemical Review,1999,68(7):605-614.

[2] BOURBIGOT S,LE BRAS M,DABROWSKI F,et al.PA-6/claynano-composite hybrid as char forming agent in intumescent formulations[J].Fire Material,2000,24(4):201-208.

[3] DABROWSKI F,LE BRAS M,CARTIER L,et al.The use of clay in an EVA-based intumescent formulation: Comparison with the intumescent formulation using polyamide-6 clay nanocomposite as carbonisation agent[J].Fire Sciences，2001,19(3):219- 241.

[4] SONG L,HU Y,ZHANG R,et al.Study on the properties of flame retardant polyurethane/organoclay nanocomposite[J].Polymer Degradation Stability,2005,87(1):111-116.

[5] ALMERAS X,DABROWSKI F,LE BRAS M,et al.Using polyamide 6 as charring agent in intumescent polypropylene formulations[J].Polymer Degradation Stability，2002,77(2):315- 323.

[6] LEVCHIK S V,LEVCHIK G P,LESNIKOVICH A I.Analysis and development of effective invariant kinetic parameters finding method based on the non-isothermal data[J].Thermochim Acta,1985,92:157-160.

[7] BOURBIGOT S,DELOBEL R,LE BRAS M,et al.Comparative study of the integral TG methods used in the invariant kinetic parameter method application to the fire-retardant polypropylene[J].Chim Phys,1993,90(10):1909-1928.

[8] BUDRUGEAC P,SEGAL E,PEREZ-MAQUEDA L A,et al.The use of the IKP method for discerning the kinetic parameters and the conversion function of the thermal dehydroclorination of PVC from non-isothermal data[J].Polymer Degradation Stability,2004,84(2):311-320.

[9] COATS A W,REDFERN J P.Kinetic parameters from thermogravimetric data[J].Journal of Polymer Science Part B: Polymer Letters,1965,3(11):917-920.

[10] TANG Y,HU Y,XIAO J,et al.PA-6 and EVA alloy/clay nano-composites as char forming agents in poly(propylene) intumescent formulations[J].Polymer Advance Technology,2005,16(4):338-343.

[11] 廖艷芬,王樹榮,駱仲泱，等.纖維素熱裂解過程動力學的實驗分析研究[J].浙江大學學報：工學版，2002,36(2): 173-176.

Kinetics of the Thermo-oxidative Decomposition of Montmorillonite Nano-composites as Charring Agent

ZHENGCheng-zhi1，ZHAOShu-hui1，YANGXiao-hui2

(1.College of Chemistry,Chemical Engineering and Biotechnology,Donghua University,Shanghai 201620,China;2.Shanghai RUCO Bio-Tech Co.Ltd.,Shanghai 201612,China)

A kind of intumescent flame retardant finishing solution was prepared by the montmorillonite nano-composites and poly-phosphate ammonium (APP),used for cotton fabric flame-retardant finishing.The thermo-oxidative decomposition of treated fabrics was studied using non-isothermal gravimetric analysis at different heating rates under the air atmosphere within the temperature from room to 900 ℃,respectively.The kinetic parameters were calculated by the invariant kinetic parameter (IKP) method which was based on the kinetic compensation effect.With the aim of determining the kinetic model,the distribution of probabilities of the solid-state conversion functions was evaluated by IKP method.Finally,the nucleation and nuclei growth as well as potential law were obtained to describe the thermo-oxidative decomposition of the treated cotton fabrics.

montmorillonite nano-composites; charring agent; thermo-oxidative decomposition kinetics；reaction mechanism

1671-0444(2015)02-0221-06

2014-01-08

鄭成志(1988—)，男，江蘇徐州人，碩士，研究方向為紡織品的功能整理.E-mail: dazhi007@163.com

趙曙輝(聯系人)，女，教授，E-mail:shzhao@dhu.edu.cn

TS 195.5

A