張中水，許進升，陳　雄，蔣　晶

(南京理工大學(xué) 機械工程學(xué)院，南京　210094)

?

EPDM包覆層熱解狀態(tài)下的力學(xué)行為及本構(gòu)模型

張中水，許進升，陳雄，蔣晶

(南京理工大學(xué) 機械工程學(xué)院，南京210094)

為了準確描述固體火箭發(fā)動機三元乙丙(EPDM)包覆層在不同熱解狀態(tài)下的力學(xué)特性，通過熱重分析的方法得到EPDM的熱解溫度范圍，并采用準靜態(tài)單軸拉伸實驗的方法獲得EPDM在不同熱解溫度下的應(yīng)力-應(yīng)變曲線。以O(shè)gden(n=2)模型為基礎(chǔ)，分別對不同熱解狀態(tài)下的力學(xué)曲線進行建模，并通過對實驗數(shù)據(jù)的擬合獲得模型參數(shù)。結(jié)果表明，在初始熱解狀態(tài)下，材料只發(fā)生了少量的失水和氣體的逃逸，其力學(xué)行為可通過粘超彈本構(gòu)模型進行描述；隨著材料熱解程度增加，基體材料和填充纖維均遭到破壞，其力學(xué)行為表現(xiàn)為超彈脆性材料特性。

三元乙丙；包覆層；熱解；粘超彈本構(gòu)模型；超彈脆性本構(gòu)模型

0　引言

三元乙丙(EPDM)作為固體火箭發(fā)動機常用的包覆材料，位于發(fā)動機殼體與裝藥之間。發(fā)動機工作時，包覆層可有效地控制裝藥燃面的變化規(guī)律，得到所需要的內(nèi)彈道性能。此外，包覆層還對殼體與藥柱之間應(yīng)力、應(yīng)變的傳遞也起到了重要的緩沖作用，緩解了殼體的震動和沖擊對裝藥燃燒的影響[1-2]。

為了保證固體推進劑在生產(chǎn)、運輸和使用的過程中在受到復(fù)雜應(yīng)力狀態(tài)下能保持結(jié)構(gòu)完整性，要求包覆層必須具有良好的力學(xué)性能。李冬[3]等通過單軸拉伸試驗的方法，得到了EPDM包覆層的力學(xué)行為與試驗溫度和加載速率有關(guān)，并從微觀上根據(jù)分子鏈段運動的原理，解釋了EPDM在準靜態(tài)拉伸過程中所表現(xiàn)出的高彈性和率相關(guān)現(xiàn)象。Cheng和Chen[4]對EPDM在準靜態(tài)條件下的Mullins效應(yīng)進行研究，研究結(jié)果表明，在常溫下EPDM是一種超粘彈性材料，其率相關(guān)性可采用應(yīng)變率參數(shù)來進行描述。作為固體火箭發(fā)動機的包覆層材料，要求EPDM必須具有良好的燒蝕性能。隨著EPDM橡膠應(yīng)用范圍越來越廣，對其熱學(xué)性能和燒蝕性能的研究越來越多[5-9]。通過模擬發(fā)動機工作過程，在高溫高壓及高速氣流和顆粒的沖刷條件下，研究EPDM的燒蝕性能，并建立EPDM在熱化學(xué)燒蝕、氣流剝蝕和顆粒侵蝕作用下的燒蝕模型已成為當代學(xué)者的主要研究課題。然而，對于EPDM作為固體火箭發(fā)動機的包覆層，在使用過程中，其熱解后力學(xué)性能的研究，國內(nèi)外鮮有報道。

為了進一步探索EPDM包覆層在發(fā)動機工作過程中的力學(xué)行為，本文對EPDM熱解后的力學(xué)性能展開研究，分析在發(fā)動機工作過程中，不同熱解程度下，EPDM的應(yīng)力響應(yīng)特性，并根據(jù)其材料力學(xué)性質(zhì)的范疇建立本構(gòu)模型。

1　實驗

本研究采用西安近代化學(xué)研究所提供的固體火箭發(fā)動機端面包覆層材料。該包覆層由傳統(tǒng)的三元乙丙復(fù)合材料制成，其中含一定比例的增強劑、耦合劑和固化劑，以滿足其生產(chǎn)和使用性能的要求。

1.1熱重分析

取少量EPDM橡膠樣品置于鋁制坩堝內(nèi)，通過天平測得樣品質(zhì)量為4.3 mg。材料在氮氣環(huán)境下以10 ℃/min的速率升溫，控制溫度由25 ℃升至700 ℃。

1.2試件制備

參照標準QJ924—85，通過使用該標準的啞鈴型刀具沖壓材料的方法，制備了厚度為6 mm的EPDM啞鈴型試件。

1.3熱解實驗

在熱解溫度范圍內(nèi)有規(guī)律地選取若干個點，將試件置于保溫箱內(nèi)加熱至不同的熱解溫度，保溫2 h后取出，在室溫環(huán)境中靜置2 h以上，使其恢復(fù)至常溫。

1.4力學(xué)實驗

實驗在微機控制電子萬能材料試驗機上完成。實驗的環(huán)境溫度為20 ℃，相對濕度為40%。在100 mm/min的加載速率下，分別對原始材料和不同熱解溫度下的熱解材料進行等速單軸拉伸實驗，獲得應(yīng)力-應(yīng)變曲線。

2　結(jié)果與討論

2.1熱解溫度范圍測定

國內(nèi)外大量學(xué)者在對于非均勻體系橡膠材料的研究過程中，都涉及到采用熱重分析的方法獲得材料的熱解溫度范圍[5,10-11]。對于含纖維填充的非均勻體材料，熱重曲線會因試件中包含填充纖維的比例的不同而有所差異，但對于讀取材料發(fā)生熱解的起止溫度，該實驗方法還是十分有效的。如圖1所示，當溫度上升至120 ℃時，曲線開始緩慢下降，EPDM開始熱解；當溫度超過480 ℃時，曲線趨于平緩，熱解反應(yīng)基本完成。

2.2原始材料和熱解材料的力學(xué)曲線

分別對原始材料和不同熱解溫度下的熱解材料進行等速單軸拉伸實驗，結(jié)合推進劑的測試標準，選用100 mm/min的加載速率。得到應(yīng)力-應(yīng)變曲線，如圖2所示。

圖1　熱重分析曲線

(a) 初始熱解狀態(tài)

(b) 深度熱解狀態(tài)

實驗結(jié)果表明，當加熱溫度低于453 K時，材料雖然已開始發(fā)生熱解，但力學(xué)行為與原始材料趨勢相近，依然屬于超粘彈材料范疇。隨著溫度繼續(xù)增加，材料基體與填充纖維均發(fā)生質(zhì)變，在該溫度范圍內(nèi)熱解的材料表現(xiàn)出彈脆性力學(xué)特性。當溫度達到673 K以后，材料質(zhì)地變得疏松易碎，已無法進行單軸拉伸實驗。

3　本構(gòu)模型的建立

3.1超粘彈本構(gòu)模型建立

根據(jù)實驗結(jié)果，如圖2(a) 所示，熱解溫度介于393～453 K的熱解材料與原始材料的單軸拉伸曲線趨勢相近，可參照超粘彈本構(gòu)模型以及時溫等效模型的建立方法[12-13]，在常溫下的超粘彈性本構(gòu)模型基礎(chǔ)上，引入應(yīng)變率參數(shù)函數(shù)及溫度參數(shù)函數(shù)，從而得到熱解材料的本構(gòu)方程：

(1)

3.1.1常溫下的超粘彈本構(gòu)

參照在常溫下關(guān)于EPDM超粘彈本構(gòu)模型的建立[14]，其超粘彈本構(gòu)方程的表達式如下：

基中

對常溫下EPDM在加載速率為100 mm/min的單軸拉伸曲線進行擬合，結(jié)果如圖3、表1、表2所示。

圖3　常溫下擬合結(jié)果與實驗結(jié)果的對比

表1　超彈本構(gòu)擬合參數(shù)

表2　率相關(guān)參數(shù)擬合結(jié)果

3.1.2熱解溫度影響因子

根據(jù)式(1)提出的本構(gòu)模型，溫度參數(shù)函數(shù)f(T,ε)的求取應(yīng)采用不同熱解溫度下的應(yīng)力與常溫下的應(yīng)力之比得到，即

(3)

式中σtest(ε,T)為單軸拉伸實驗中不同熱解溫度下的應(yīng)力值；T0在文中默認為常溫，T0=273 K。

由于實驗初始階段誤差較大，曲線的變化規(guī)律不明顯，因此只選取ε>0.1的實驗數(shù)據(jù)進行分析計算。計算結(jié)果如圖4所示，不同熱解溫度下的參數(shù)函數(shù)f(ε,T)與應(yīng)變的關(guān)系曲線可由OriginPro中的遞減指數(shù)函數(shù)ExpDecay表示，其表達式為

(4)

設(shè)定初始應(yīng)變ε0=0.1，采用式(4)對熱解溫度參數(shù)函數(shù)進行擬合，擬合結(jié)果如圖4所示，得到初始溫度參數(shù)f0=0.85，得到不同熱解溫度下參數(shù)Ai、ti的擬合結(jié)果如表3所示。

圖4　計算結(jié)果的擬合曲線

表3　計算結(jié)果參數(shù)擬合結(jié)果

對于參數(shù)xi、yi關(guān)于熱解溫度的變化規(guī)律經(jīng)統(tǒng)計，可采用OriginPro中的遞增指數(shù)函數(shù)ExpGrow擬合，其表達式為

(5)

(6)

擬合結(jié)果如圖5、表4所示。

表4　熱解溫度參數(shù)擬合結(jié)果

圖5　熱解溫度參數(shù)擬合曲線

將參數(shù)函數(shù)的表達式代入式(1)，得到熱解材料的超粘彈本構(gòu)方程表達式如下：

(7)

將表1～表4中的參數(shù)代入式(7)，得到不同熱解溫度下的應(yīng)力-應(yīng)變曲線，將其與實驗結(jié)果進行對比，如圖6所示。可見，實驗結(jié)果與擬合結(jié)果吻合較好。因此，可選用該本構(gòu)模型對EPDM在393～453 K的熱解溫度下的力學(xué)行為進行描述。

圖6　EPDM 不同熱解溫度下擬合曲線與實驗結(jié)果的對比

3.2超彈脆性模型建立

參照超彈性材料本構(gòu)模型的建立[15]，根據(jù)熱力學(xué)基本定律，在均熱與小變形條件下，其本構(gòu)關(guān)系表達形式如下：

(8)

應(yīng)變能函數(shù)選取用于擬合原始材料超彈性部分的Ogden模型。Ogden模型采用主伸長來表征應(yīng)變能函數(shù)，從而克服了采用變形張量不變量導(dǎo)致的關(guān)系式復(fù)雜[16]，其應(yīng)變能函數(shù)表達式如下：

(9)

對于基于伸長率的應(yīng)變能函數(shù)，名義應(yīng)力可直接由應(yīng)變能函數(shù)對伸長率求偏導(dǎo)得到。

由σ2=σ3=0，可得到：

(10)

(11)

如圖2(b)所示，當熱解溫度大于300 ℃時，材料的最大拉伸應(yīng)力應(yīng)變均顯著降低，試驗機所造成的實驗誤差較大，無法獲取較準確的模型參數(shù)。因此，只對熱解溫度介于200～300 ℃的實驗曲線進行擬合，擬合結(jié)果如圖7、表5所示。

圖7　彈脆性力學(xué)曲線擬合結(jié)果和試驗結(jié)果的對比

表5　參數(shù)擬合結(jié)果

4　結(jié)論

(1)通過熱重分析，得到EPDM包覆材料的熱解溫度范圍為120～480 ℃。對不同熱解溫度下的EPDM試件進行等速單軸拉伸試驗，得到在熱解溫度為120～180 ℃的范圍內(nèi)，熱解材料與原始材料力學(xué)曲線相近，屬于超粘彈材料范疇；在熱解溫度為180～400 ℃的范圍內(nèi)，材料的延伸率顯著降低，表現(xiàn)出彈脆性特性。

(2)在超粘彈本構(gòu)模型的基礎(chǔ)上，引入熱解溫度參數(shù)函數(shù)，得到EPDM在初始熱解溫度下的本構(gòu)方程，經(jīng)擬合獲得模型參數(shù)的變化規(guī)律，從而可預(yù)測EPDM在初始熱解溫度范圍內(nèi)任意溫度下的力學(xué)響應(yīng)特性；推導(dǎo)出以O(shè)gden模型為基礎(chǔ)的彈脆性本構(gòu)方程，對深度熱解程度下的EPDM力學(xué)曲線進行擬合，獲得模型參數(shù)。擬合結(jié)果表明，含熱解溫度項的超粘彈模型和彈脆性模型均可很好地描述EPDM在不同熱解程度下的力學(xué)行為，為EPDM絕熱層在固體火箭發(fā)動機工作過程中的應(yīng)力、應(yīng)變傳遞提供了一定的理論依據(jù)。

[1]張平偉. 三元乙丙內(nèi)絕熱材料及其性能研究[D]. 國防科學(xué)技術(shù)大學(xué), 2009.

[2]徐海平. 絕熱材料炭化層的性能研究[D]. 內(nèi)蒙古大學(xué), 2010.

[3]李冬，陳竚，任黎，等. 三元乙丙基固體推進劑包覆層拉伸性能影響因素研究 [J]. 化工新型材料，2011，39(9)：98-99.

[4]Cheng M，Chen W. Experimental investigation of the stress-stretch behavior of EPDM rubber with loading rate effects [J]. International Journal of Solids and Structures，2003，40(18)：4749-4768.

[5]Natali M, Rallini M, Puglia D, et al. EPDM based heat shielding materials for solid rocket motors: A comparative study of different fibrous reinforcements[J]. Polymer Degradation and Stability, 2013, 98(11):2131-2139.

[6]Komalan C, George K E, Varughese K T, et al. Thermogravimetric and wide angle X-ray diffraction analysis of thermoplastic elastomers from nylon copolymer and EPDM rubber[J]. Polymer Degradation and Stability, 2008, 93(12):2104-2112.

[7]Singh S, Guchhait P K, Chaki T K. Development of polyimide-nanosilica filled EPDM based light rocket motor insulator compound: Influence of polyimide-nanosilica loading on thermal, ablation, and mechanical properties[J]. Composites Part A, 2013, 44(1):8-15.

[8]楊颯，何國強，李江，等. 碳化層疏松/致密結(jié)構(gòu)的三元乙丙燒蝕模型 [J]. 航空動力學(xué)報, 2012, 27(5), 1172-1178.

[9]陳劍，李江，李強，等. EPDM絕熱材料碳化層結(jié)構(gòu)特征及其對燒蝕的影響 [J]. 固體火箭技術(shù)，2011, 34(1), 122-130.

[10]王天波，鞠玉濤. Mechanical properties with high temperature and the microstructure of carbon/phenolic ablative composites[J]. 武漢理工大學(xué)學(xué)報：材料科學(xué)英文版, 2012, 27(5):967-972.

[11]Komalan C, George K E, Mathew V S, et al. Thermogravimetric and wide angle X-ray diffraction analysis of thermoplastic elastomers from nylon copolymer and EPDM rubber[J]. Polymer Degradation & Stability, 2008, 93:2104-2112.

[12]Zhang J F, Zheng J, Chen X, et al. A thermovisco-hyperelastic constitutive model of NEPE propellant over a large range of strain rates[J]. Journal of Engineering Materials and Technology, 2014, 136(3):337-346.

[13]Yu P, Yao X, Han Q, et al. A visco-elastoplastic constitutive model for large deformation response of polycarbonate over a wide range of strain rates and temperatures[J]. Polymer, 2014, 55:6577-6593.

[14]張中水，陳雄，周清春，等. EPDM絕熱層的粘超彈本構(gòu)模型[J]. 固體火箭技術(shù)，2015,38(2): 273-277.

[15]常武軍, 鞠玉濤, 胡少青. HTPB固化膠片的超彈性本構(gòu)模型[J]. 推進技術(shù), 2012, 33(5): 795-798.

[16]Ogden R. Large deformation isotropic elasticity on the correlation of theory and experiment for incompressible rubberlike solids[J]. Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences，1972，326(1567)：565-584.

(編輯：薛永利)

Mechanical behavior and constitutive model of EPDM inhibitor

ZHANG Zhong-shui, XU Jin-sheng, CHEN Xiong, JIANG Jing

(School of Mechanical Engineering，NUST，Nanjing210094，China)

To describe the mechanical behavior of EPDM inhibitor of solid rocket in different thermal degradation, the range of thermal degradation temperature of EPDM was determined by thermogravimetric analysis. Stress-strain curves of EPDM under different degree of pyrolysis were obtained by quasi-static mechanics experiments. Based on Ogden(n=2) model, constitutive model were established respectively to obtain the model parameters by fitting mechanics characteristic curves of EPDM under different degrees of pyrolysis. The results show that primarily water and gas are released at the initial stage of pyrolysis, and mechanical behaviors of the composite could be described by visco-hyperelastic constitutive model. With the increase of degree of pyrolysis, rubber matrix and filler fibers were degraded and the composite presents superelastic-brittle characteristics.

EPDM；inhibitor；pyrolysis；visco-hyperelastic constitutive model；superelastic-brittle constitutive model

2015-06-05；

2015-08-11。

高能聚合物沖擊載荷下熱-力耦合模型研究(BK20140772)。

張中水(1990—)，女，碩士生，研究方向為航空宇航推進理論與工程。E-mail:zhongshuizhang@126.com

V258

A

1006-2793(2016)04-0560-05

10.7673/j.issn.1006-2793.2016.04.020