常應(yīng)變率下某固體推進(jìn)劑非線性粘彈性本構(gòu)關(guān)系研究*

陳 鑫，賈有軍，郜 婕，李 偉，李錄賢

(1.西安交通大學(xué) 航天航空學(xué)院，機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室/陜西省先進(jìn)飛行器服役環(huán)境與控制重點(diǎn)實(shí)驗(yàn)室，西安 710049；2.中國(guó)航天科技集團(tuán)四院四十一所，西安 710025)

0 引言

固體推進(jìn)劑發(fā)動(dòng)機(jī)具有大推力、高可靠性、方便使用維護(hù)等優(yōu)點(diǎn)，已成為航天運(yùn)載火箭的主要?jiǎng)恿ν七M(jìn)裝置，在世界各國(guó)航天運(yùn)載技術(shù)發(fā)展中發(fā)揮著重要作用。固體推進(jìn)劑是固體火箭發(fā)動(dòng)機(jī)的能源材料，也是結(jié)構(gòu)材料，其力學(xué)性能直接影響著整個(gè)發(fā)動(dòng)機(jī)的性能。

固體推進(jìn)劑力學(xué)性能的主要研究?jī)?nèi)容之一是固體推進(jìn)劑的本構(gòu)關(guān)系。固體推進(jìn)劑在制造、貯存、使用等過(guò)程中，其力學(xué)性能受到溫度、應(yīng)變率的影響[1]。目前，國(guó)內(nèi)外學(xué)者已經(jīng)根據(jù)各種實(shí)驗(yàn)對(duì)固體推進(jìn)劑的力學(xué)性能進(jìn)行了廣泛深入的研究，并建立相應(yīng)的本構(gòu)模型，但由于固體推進(jìn)劑的力學(xué)性能十分復(fù)雜，已有的相關(guān)研究與實(shí)際需求仍然存在較大差距。如Ho[2]建立了考慮損傷、溫度和應(yīng)變率因素的非線性本構(gòu)模型，能夠預(yù)測(cè)固體推進(jìn)劑高應(yīng)變率下的非線性粘彈性力學(xué)響應(yīng)，但溫度較低時(shí)，預(yù)測(cè)效果不是很好。朱-王-唐本構(gòu)模型[3-4]雖較好地考慮了準(zhǔn)靜態(tài)和動(dòng)態(tài)條件，被大量運(yùn)用于高聚物、炸藥等脆性材料中，但只適用于變形小于20%的情況，對(duì)于高達(dá)80%以上的應(yīng)變，仍需進(jìn)一步完善。Schapery課題組[5-7]、孟紅磊等[8]、姚東等[9]在這方面也做了大量工作，其采用宏觀物理量的變化來(lái)描述固體推進(jìn)劑變形時(shí)所表現(xiàn)出的非線性行為，回避了對(duì)微觀結(jié)構(gòu)的復(fù)雜分析。但是，這些非線性本構(gòu)模型中一般都含有多個(gè)物理意義含糊的材料參數(shù)，且所采用的是與線性本構(gòu)模型基本無(wú)關(guān)的另一套體系。

固體推進(jìn)劑是典型的粘彈性材料，其粘彈性行為因應(yīng)變率不同而顯著不同。實(shí)驗(yàn)表明，不同應(yīng)變率下固體推進(jìn)劑的拉伸曲線差異非常明顯。基于該實(shí)驗(yàn)現(xiàn)象，本文從線性粘彈性本構(gòu)關(guān)系出發(fā)，研究應(yīng)變率對(duì)固體推進(jìn)劑材料力學(xué)性能的影響，擬發(fā)展一種參數(shù)數(shù)量少、物理意義明確的固體推進(jìn)劑材料的唯象非線性粘彈性本構(gòu)關(guān)系，以描述固體推進(jìn)劑的實(shí)驗(yàn)結(jié)果，方便應(yīng)用于藥柱結(jié)構(gòu)的實(shí)際分析，為高裝填比固體火箭發(fā)動(dòng)機(jī)裝藥結(jié)構(gòu)設(shè)計(jì)提供支持。

1 線性粘彈性基本理論

單軸應(yīng)力狀態(tài)下材料的線性粘彈性本構(gòu)關(guān)系用積分形式可表示為

(1)

式中σ為應(yīng)力；ε為應(yīng)變；t為時(shí)間；E(t)為材料的松弛模量。

松弛模量E(t)可通過(guò)將式(1)應(yīng)用于松弛實(shí)驗(yàn)來(lái)獲得。松弛實(shí)驗(yàn)是在初始時(shí)刻給定一個(gè)突變的初始應(yīng)變并予以保持、然后測(cè)取應(yīng)力隨時(shí)間變化的實(shí)驗(yàn)。也就是說(shuō)，在松弛實(shí)驗(yàn)中，其應(yīng)變變化為

ε(t)=ε0H(t-0)

(2)

式中H(t-t0)為Heaviside函數(shù)。

這樣，松弛實(shí)驗(yàn)的應(yīng)變率為

(3)

式中δ(t-t0)為Dirac delta函數(shù)。

將式(3)代入式(1)，對(duì)于松弛實(shí)驗(yàn)，得

(4)

一般地，松弛模量可用Prony級(jí)數(shù)表示：

在給定環(huán)境溫度下，對(duì)于ε0=5%，松弛實(shí)驗(yàn)得到的松弛模量隨時(shí)間的變化如圖2中實(shí)線所示。

圖1 廣義Maxwell模型Fig.1 Generalized Maxwell model

圖2 松弛模量的實(shí)驗(yàn)值與理論計(jì)算值Fig.2 Relaxation modulus from test and theory

(6)

將式(6)中的參數(shù)值代入式(5)中，即得到松弛模量的表達(dá)式，預(yù)測(cè)的松弛模量如圖2中各點(diǎn)所示?？梢钥闯觯c實(shí)驗(yàn)值吻合良好。

得到松弛模量后，就建立起了如式(1)的線性粘彈性本構(gòu)關(guān)系。下面將該關(guān)系用于常應(yīng)變率拉伸實(shí)驗(yàn)。

(7)

式(7)即為依據(jù)線性粘彈性本構(gòu)關(guān)系得到的常應(yīng)變率時(shí)的應(yīng)力隨時(shí)間變化關(guān)系。

考察應(yīng)變率為0.02 s-1的情形，實(shí)驗(yàn)獲得的應(yīng)力-時(shí)間曲線如圖3中的實(shí)線所示，根據(jù)式(7)預(yù)測(cè)的應(yīng)力-時(shí)間曲線如圖3中的虛線所示?？梢钥闯?，二者只在6%以內(nèi)的小變形階段具有較好的一致性；在其余的較大應(yīng)變范圍內(nèi)，二者出現(xiàn)較大偏差。因而，需要研究非線性粘彈性本構(gòu)關(guān)系，以期在更寬應(yīng)變范圍或整個(gè)過(guò)程與實(shí)驗(yàn)結(jié)果一致。

圖3 常應(yīng)變率實(shí)驗(yàn)結(jié)果與線性粘彈性理論計(jì)算結(jié)果對(duì)比Fig.3 Comparison of the results from the constant- strain-rate test and the theory

2 非線性粘彈性基本理論

2.1 推進(jìn)劑材料非線性性能的影響因素

在第1節(jié)中，利用松弛實(shí)驗(yàn)獲得松弛模量，進(jìn)而建立了線性粘彈性本構(gòu)關(guān)系。但是，如此得到的線性粘彈性本構(gòu)關(guān)系對(duì)常應(yīng)變率拉伸實(shí)驗(yàn)預(yù)測(cè)的效果并不理想，只在小變形階段適用。隨著應(yīng)變?cè)龃?，線性本構(gòu)關(guān)系與實(shí)驗(yàn)曲線出現(xiàn)明顯差異，而且越來(lái)越大，說(shuō)明線性粘彈性本構(gòu)關(guān)系不再適用。

雖然根據(jù)式(6)的松弛模量參數(shù)獲得的常應(yīng)變率應(yīng)力時(shí)間關(guān)系式(7)不能很好描述實(shí)驗(yàn)結(jié)果，但分析式(7)函數(shù)的基本形態(tài)是一件饒有興趣的事。直接采用式(7)對(duì)0.02 s-1的實(shí)驗(yàn)結(jié)果進(jìn)行擬合，得到的參數(shù)值為

(8)

式(7)表示的應(yīng)力隨時(shí)間變化規(guī)律將如圖4所示?？梢钥闯觯允?8)代替式(6)、式(7)與實(shí)驗(yàn)曲線的吻合程度在整個(gè)變形過(guò)程都得到了顯著改觀。

圖4 改進(jìn)預(yù)測(cè)與實(shí)驗(yàn)對(duì)比Fig.4 Comparison of the results from the improved theory and the test

由此看來(lái)，在廣義Maxwell模型中采用恰當(dāng)?shù)膹椥阅Ａ亢驼承韵禂?shù)，就可很好預(yù)測(cè)常應(yīng)變率實(shí)驗(yàn)結(jié)果，進(jìn)而獲得相應(yīng)的本構(gòu)關(guān)系。由于這種本構(gòu)關(guān)系中材料參數(shù)是根據(jù)常應(yīng)變率實(shí)驗(yàn)結(jié)果擬合而得，與松弛實(shí)驗(yàn)獲得的相應(yīng)參數(shù)并不相同，其機(jī)理可解釋為粘彈性元件中的彈簧和粘壺性能因應(yīng)變率效應(yīng)發(fā)生了改變，因而，將基于常應(yīng)變率實(shí)驗(yàn)擬合得到的參數(shù)(如式(8))所對(duì)應(yīng)的粘彈性本構(gòu)關(guān)系稱為應(yīng)變率相關(guān)的非線性粘彈性本構(gòu)關(guān)系。

2.2 非線性粘彈性本構(gòu)關(guān)系

由2.1節(jié)分析可知，應(yīng)變率效應(yīng)使得粘彈性元件的材料參數(shù)發(fā)生變化，引起了材料的非線性粘彈性行為。據(jù)此，仿照積分型線性粘彈性本構(gòu)關(guān)系，提出積分型非線性粘彈性本構(gòu)關(guān)系為

(9)

式中ER(t)為非線性松弛模量，由常應(yīng)變率實(shí)驗(yàn)確定。

仿照式(5)，ER(t)也可表示為

(10)

為便于實(shí)際應(yīng)用并保證精度，對(duì)式(10)仍取二項(xiàng)三參數(shù)形式，即

(11)

其物理意義可表示為如圖5所示的非線性三參量固體。

σ(t)=at+b(1-e-ct)

(12)

根據(jù)式(7)，不難得到

(13)

根據(jù)式(13)，可反演得到材料參數(shù)為

(14)

圖5 非線性三參量固體材料參數(shù)的物理解釋Fig.5 Physical interpretation of the three nonlinear parameters

式(12)的函數(shù)形式可描述多種不同形態(tài)的變化關(guān)系。例如，當(dāng)b=5、c=0.1及a=-0.01時(shí)，式(12)可描述先增加后減小的形態(tài)，如圖6(a)所示；當(dāng)b=5、c=0.1及a=0時(shí)，式(12)可描述先增加、然后基本保持不變的形態(tài)，如圖6(b)所示；當(dāng)b=5、c=0.1及a=+0.01時(shí)，式(12)可描述單調(diào)增加的形態(tài)，如圖6(c)所示。需要特別指出的是，這三種基本形態(tài)幾乎涵蓋了常應(yīng)變率實(shí)驗(yàn)結(jié)果中應(yīng)力隨時(shí)間變化的各種規(guī)律，因而可用來(lái)刻畫應(yīng)變率對(duì)材料性能的影響。

3 非線性粘彈性本構(gòu)關(guān)系的應(yīng)用

3.1 本構(gòu)參數(shù)識(shí)別方法

由于具有應(yīng)變率相關(guān)性，雖然與線性粘彈性本構(gòu)關(guān)系具有相同的函數(shù)形式，但非線性粘彈性本構(gòu)關(guān)系的材料參數(shù)不能再借助于松弛實(shí)驗(yàn)獲得，而是通過(guò)常應(yīng)變率實(shí)驗(yàn)結(jié)果擬合得到，以實(shí)現(xiàn)更寬應(yīng)變范圍或整個(gè)變形過(guò)程中與實(shí)驗(yàn)結(jié)果一致。

Step-1：將給定應(yīng)變率下的應(yīng)力-時(shí)間實(shí)驗(yàn)結(jié)果準(zhǔn)備成Excel格式；

Step-2：在MATLAB軟件的Curve Fitting中自定義形如式(12)的函數(shù)；

Step-3：將Step-1中Excel格式的實(shí)驗(yàn)結(jié)果通過(guò)Import Data接口輸入到MATLAB中；

Step-4：?jiǎn)?dòng)MATLAB中Curve Fitting的Fit功能，獲得式(12)中參數(shù)a、b和c的值；

(a)先增加后減小的形態(tài)

(b)基本保持不變的形態(tài)

(c)單調(diào)增加的形態(tài)

3.2 常應(yīng)變率拉伸實(shí)驗(yàn)預(yù)測(cè)

某環(huán)境溫度下不同常應(yīng)變率時(shí)的實(shí)驗(yàn)結(jié)果如圖7中實(shí)線所示。

按照3.1節(jié)中的參數(shù)識(shí)別方法，根據(jù)圖7的實(shí)驗(yàn)曲線，對(duì)式(12)中的3個(gè)參數(shù)a、b和c分別進(jìn)行識(shí)別，再利用式(12)預(yù)測(cè)得到圖7中虛線所示的理論值?？梢钥闯?，非線性粘彈性本構(gòu)關(guān)系在0.000 1～0.1 s-1的5個(gè)不同常應(yīng)變率條件下，100%應(yīng)變范圍內(nèi)，預(yù)測(cè)的應(yīng)力-時(shí)間關(guān)系都與實(shí)驗(yàn)結(jié)果吻合良好。

(a) 應(yīng)變率為0.000 1 s-1 (b) 應(yīng)變率為0.000 4 s-1 (c) 應(yīng)變率為0.01 s-1

(d) 應(yīng)變率為0.02 s-1 (e) 應(yīng)變率為0.1 s-1

3.3 參數(shù)值統(tǒng)計(jì)

按照3.1節(jié)方法識(shí)別得到的函數(shù)參數(shù)和根據(jù)式(14)反演得到的非線性材料參數(shù)如表1所示。

(15)

表1 識(shí)別得到的3個(gè)非線性材料參數(shù)值Table1 Identified values of the three nonlinear constants

將式(15)通過(guò)式(11)代入式(9)，最終得到某溫度下考慮應(yīng)變率效應(yīng)的非線性粘彈性本構(gòu)關(guān)系為

(16)

隨應(yīng)變率變化隨應(yīng)變率變化隨應(yīng)變率變化

3.4 討論

應(yīng)力-時(shí)間的關(guān)系(式(12))具有如圖6所示的基本變化形態(tài)，但是，圖7(e)所示的實(shí)驗(yàn)曲線在應(yīng)變較大時(shí)呈現(xiàn)出了略微的起伏變化，因而，與其他相比，式(12)對(duì)圖7(e)變化的描述尚待進(jìn)一步改進(jìn)。

實(shí)際上，在式(12)基礎(chǔ)上再增加一項(xiàng)，可使應(yīng)力-時(shí)間關(guān)系變?yōu)?/p>

σ(t)=at+b(1-e-ct)+nt2

(17)

利用3.1節(jié)的識(shí)別方法和圖7(e)的實(shí)驗(yàn)數(shù)據(jù)，得到函數(shù)參數(shù)值分別為

(18)

圖9 式(17)的改進(jìn)預(yù)測(cè)結(jié)果Fig.9 Improved prediction by Eq.(17)

式(17)的預(yù)測(cè)與實(shí)驗(yàn)結(jié)果比較如圖9所示，與式(12)預(yù)測(cè)值(圖7(e)中的虛線)相比，式(17)的預(yù)測(cè)(圖9中的虛線)與實(shí)驗(yàn)結(jié)果符合更好。

此時(shí)，令

(19)

對(duì)式(19)兩邊求導(dǎo)得材料的非線性松弛模量:

(20)

但式(20)再不能像式(14)那樣反演得到具有明確物理意義的材料參數(shù)。但從實(shí)際使用角度，數(shù)量很有限的唯象參數(shù)a、b、c及n仍可方便地用于實(shí)際的工程結(jié)構(gòu)分析。

4 結(jié)論

(1) 對(duì)因應(yīng)變率引起的非線性粘彈性行為，可通過(guò)引入材料的非線性松弛模量概念，由線性粘彈性本構(gòu)關(guān)系發(fā)展成為非線性粘彈性本構(gòu)關(guān)系。

(2) 非線性本構(gòu)關(guān)系的材料參數(shù)雖然由實(shí)驗(yàn)數(shù)據(jù)唯象地?cái)M合得到，但在物理上可解釋為三參數(shù)固體中的非線性彈簧系數(shù)和非線性粘壺系數(shù)。

(3) 應(yīng)變率效應(yīng)可用較少(例如3個(gè))的材料參數(shù)對(duì)整個(gè)變形過(guò)程予以描述，符合工程領(lǐng)域設(shè)計(jì)的簡(jiǎn)單性要求。

由于實(shí)驗(yàn)數(shù)據(jù)量的限制及較大的應(yīng)變率跨度，本文非線性粘彈性本構(gòu)關(guān)系材料參數(shù)隨應(yīng)變率的變化規(guī)律尚不明晰，需要更多實(shí)驗(yàn)的支持。實(shí)際上，假以豐富的大溫度范圍多個(gè)常應(yīng)變率拉伸實(shí)驗(yàn)數(shù)據(jù)，本文非線性粘彈性本構(gòu)關(guān)系的研究思想也可進(jìn)一步延伸，以考慮對(duì)推進(jìn)劑材料性能顯著影響的溫度因素。