職世君，曹付齊，申志彬，韓 波，李記威

(1.中國(guó)空空導(dǎo)彈研究院，洛陽(yáng) 471000；2.國(guó)防科技大學(xué) 航天科學(xué)與工程學(xué)院，長(zhǎng)沙 410073)

復(fù)合固體推進(jìn)劑雙折線脫濕損傷模型參數(shù)影響分析

職世君1，曹付齊1，申志彬2，韓 波1，李記威1

(1.中國(guó)空空導(dǎo)彈研究院，洛陽(yáng) 471000；2.國(guó)防科技大學(xué) 航天科學(xué)與工程學(xué)院，長(zhǎng)沙 410073)

為分析雙折線損傷模型參數(shù)對(duì)復(fù)合固體推進(jìn)劑細(xì)觀損傷及宏觀非線性力學(xué)性能的影響，采用分子動(dòng)力學(xué)方法建立復(fù)合固體推進(jìn)劑顆粒夾雜模型，根據(jù)Surface-based cohesive方法，在高氯酸銨(AP)顆粒與基體之間的界面處設(shè)置接觸損傷。利用有限元方法對(duì)具有不同損傷參數(shù)的顆粒夾雜模型進(jìn)行計(jì)算，并對(duì)比數(shù)值仿真結(jié)果。結(jié)果表明，損傷起始應(yīng)力對(duì)復(fù)合固體推進(jìn)劑抗拉強(qiáng)度、最大延伸率有較大的影響；界面初始剛度在一定程度內(nèi)的變化，對(duì)復(fù)合固體推進(jìn)劑宏觀力學(xué)性能及細(xì)觀損傷形貌影響較小；界面失效距離主要影響復(fù)合固體推進(jìn)劑的最大延伸率。

固體推進(jìn)劑；細(xì)觀力學(xué)；界面損傷；數(shù)值模擬

0 引言

復(fù)合固體推進(jìn)劑是一種多相體(基體、增強(qiáng)相、界面相等)復(fù)合材料。該類型推進(jìn)劑具有較良好的力學(xué)性能，燃燒穩(wěn)定。因此，在導(dǎo)彈和宇航火箭發(fā)動(dòng)機(jī)中被廣泛應(yīng)用。復(fù)合固體推進(jìn)劑的力學(xué)性能和損傷破壞規(guī)律受增強(qiáng)相的體積分?jǐn)?shù)及組分材料性質(zhì)影響較大，但同時(shí)也取決于其細(xì)觀結(jié)構(gòu)特征，這些特征包括增強(qiáng)相的分布規(guī)律、形狀以及界面相的性質(zhì)等。早期對(duì)復(fù)合固體推進(jìn)劑力學(xué)性能的研究，主要基于連續(xù)介質(zhì)力學(xué)，不考慮材料的內(nèi)部細(xì)觀結(jié)構(gòu)變化，通過(guò)大量的固體推進(jìn)劑力學(xué)性能試驗(yàn)確定固體推進(jìn)劑的非線性本構(gòu)關(guān)系[1-3]。隨著對(duì)固體推進(jìn)劑力學(xué)性能研究的深入，微裂紋和損傷的影響逐漸引起了許多學(xué)者的重視。根據(jù)固體推進(jìn)劑在拉伸載荷下的體積膨脹現(xiàn)象及體積膨脹和材料損傷之間的關(guān)系，引進(jìn)一些數(shù)學(xué)函數(shù)來(lái)描述固體推進(jìn)劑的損傷[4-6]。近年來(lái)，隨著計(jì)算機(jī)性能的大幅度提高，計(jì)算細(xì)觀力學(xué)得到了迅速發(fā)展。Matous等[7-8]利用他們的自主開發(fā)軟件Rocpack生成了固體推進(jìn)劑顆粒夾雜模型，并在顆粒與基體之間的界面層設(shè)置了cohesive單元模擬了固體顆粒和基體之間損傷的產(chǎn)生及發(fā)展。Tan等[9-10]通過(guò)數(shù)字圖象等相關(guān)技術(shù)獲得了高能炸藥PBX9501緊湊拉伸試樣裂尖周圍的應(yīng)力場(chǎng)及位移場(chǎng)，利用擴(kuò)展的Mori-Tanaka方法對(duì)試驗(yàn)結(jié)果做了均勻化處理，從而提取了PBX9501炸藥中固體顆粒與基體之間界面損傷的相關(guān)參數(shù)。隨后，又結(jié)合試驗(yàn)研究得到顆粒與基體之間的非線性界面損傷模型，對(duì)含不同粒徑顆粒的細(xì)觀模型進(jìn)行了數(shù)值計(jì)算。國(guó)內(nèi)一些學(xué)者也展開了復(fù)合固體推進(jìn)劑的細(xì)觀力學(xué)性能研究。有的學(xué)者采用有限元法對(duì)復(fù)合固體推進(jìn)劑進(jìn)行直接數(shù)值模擬，對(duì)推進(jìn)劑內(nèi)部界面脫粘過(guò)程進(jìn)行了有限元分析[11]。有的學(xué)者根據(jù)數(shù)值仿真結(jié)果，結(jié)合細(xì)觀力學(xué)方法，如：Mori-Tanaka方法或改進(jìn)的Mori-Tanaka方法，研究了固體推進(jìn)劑的模量、界面脫粘對(duì)固體推進(jìn)劑力學(xué)性能的影響等[12-14]，或采用多步法，將基體與部分顆粒均質(zhì)化為一種混合物，計(jì)算出較復(fù)雜的固體推進(jìn)劑等效模量[15]。

在之前的工作中，已對(duì)復(fù)合固體推進(jìn)劑顆粒夾雜模型的幾何建模及細(xì)觀尺度上復(fù)合固體推進(jìn)劑損傷的產(chǎn)生、演化和聚合及其對(duì)宏觀非線性力學(xué)性能的影響進(jìn)行了研究[16-17]。本文主要在之前的研究基礎(chǔ)上，討論復(fù)合固體推進(jìn)劑雙折線脫濕損傷模型參數(shù)對(duì)其細(xì)觀損傷形貌及宏觀力學(xué)性能的影響。

1 物理模型和計(jì)算方法

1.1 復(fù)合固體推進(jìn)劑組分及材料屬性

根據(jù)復(fù)合固體推進(jìn)劑的各組分含量，確定AP顆粒的填充體積分?jǐn)?shù)為65.3%。采用分子動(dòng)力學(xué)方法[17]，建立復(fù)合固體推進(jìn)劑顆粒夾雜模型，如圖1所示。復(fù)合固體推進(jìn)劑基體松弛模量如圖2所示。

假設(shè)AP顆粒為彈性體，取其彈性模量和泊松比分別為E=32 450 MPa，ν=0.143 3[7]。

1.2 Surface-based cohesive方法

Surface-based cohesive方法是基于面與面的接觸行為而定義的損傷模型。接觸點(diǎn)對(duì)的確定采用的是從屬表面上的離散點(diǎn)對(duì)主控表面進(jìn)行投影，如圖3所示。

從圖3可看出，從屬表面上的As、Bs和Cs的3個(gè)節(jié)點(diǎn)分別向主控表面進(jìn)行投影，取與其距離最短的點(diǎn)為接觸點(diǎn)對(duì)。在計(jì)算時(shí)，若對(duì)應(yīng)點(diǎn)不在主控表面的網(wǎng)格節(jié)點(diǎn)處，其具體值如位移和接觸壓力等，則可根據(jù)鄰近節(jié)點(diǎn)的計(jì)算結(jié)果進(jìn)行插值計(jì)算得到，從而建立了主控表面和從屬面的一一對(duì)應(yīng)關(guān)系。在接觸點(diǎn)對(duì)之間引入cohesive損傷模型，便是本文所采用的Surface-based cohesive方法。

界面的力學(xué)響應(yīng)可描述為

(1)

式中t為界面的名義應(yīng)力矩陣，包含2個(gè)含量，分別代表法向和切向的名義應(yīng)力；K為界面剛度矩陣。

損傷起始準(zhǔn)則取最大名義應(yīng)力準(zhǔn)則，即

(2)

式中tn為法向應(yīng)力；ts為切向應(yīng)力；tn0、ts0分別為對(duì)應(yīng)的臨界應(yīng)力；〈〉為麥考利符號(hào)。

當(dāng)界面位移δ達(dá)到δ0時(shí)，界面開始損傷，界面應(yīng)力按式(3)、式(4)計(jì)算：

(3)

(4)

式中tn和ts分別為無(wú)損傷時(shí)計(jì)算所得到的應(yīng)力分量；D為損傷因子。

定義損傷因子為

(5)

從圖4可看出：

(6)

(7)

故可求得損傷因子為

(8)

式中δ為計(jì)算過(guò)程中界面張開位移。

2 計(jì)算結(jié)果與討論

復(fù)合固體推進(jìn)劑典型的損傷過(guò)程一般可分為3個(gè)階段。第1階段：復(fù)合固體推進(jìn)劑無(wú)損傷，顆粒周圍真空孔數(shù)量很低，或者為零。第2階段：隨著載荷的繼續(xù)增大，發(fā)生脫濕損傷的顆粒逐漸增多，顆粒周圍真空孔的數(shù)量增加，并在某一點(diǎn)處達(dá)到最大值。第3階段：隨著顆粒脫濕損傷的繼續(xù)演化，前一階段形成的真空孔不斷增大，損傷不斷聚合并伴隨基體材料的損傷、斷裂，最后導(dǎo)致整個(gè)復(fù)合固體推進(jìn)劑斷裂失效[18]?？梢姡趶?fù)合固體推進(jìn)劑斷裂失效的過(guò)程中，顆粒脫濕損傷的影響非常大。

根據(jù)圖1所示的復(fù)合固體推進(jìn)劑顆粒夾雜模型，對(duì)有限元程序進(jìn)行二次開發(fā)，完成顆粒與基體之間接觸損傷的自定義，實(shí)現(xiàn)復(fù)合固體推進(jìn)劑細(xì)觀含損傷模型的自動(dòng)化建模。通過(guò)對(duì)界面損傷參數(shù)反演識(shí)別，當(dāng)取界面損傷參數(shù)如表1所示時(shí)，采用二維顆粒夾雜模型計(jì)算固體推進(jìn)劑的損傷過(guò)程，其表現(xiàn)出的宏觀應(yīng)力-應(yīng)變曲線與試驗(yàn)結(jié)果吻合良好，如圖5所示。在損傷的第2階段，仿真計(jì)算的應(yīng)力-應(yīng)變曲線沒(méi)有形成較明顯的過(guò)渡段，即曲線的第2階段表現(xiàn)的不明顯。這說(shuō)明采用的界面損傷模型仍需要進(jìn)一步的完善，但從整體來(lái)看，仿真計(jì)算的應(yīng)力-應(yīng)變曲線與試驗(yàn)曲線較接近，說(shuō)明該損傷模型可近似表征，并模擬復(fù)合固體推進(jìn)劑的損傷過(guò)程。

表1 界面損傷參數(shù)

圖6為復(fù)合固體推進(jìn)劑顆粒夾雜模型在拉伸載荷下失效時(shí)的損傷形貌。從圖6可看出，當(dāng)推進(jìn)劑在拉伸載荷下失效時(shí)，部分顆粒脫濕較嚴(yán)重，顆粒與基體之間形成了較大的空隙。在受載過(guò)程中，該區(qū)域內(nèi)的基體材料逐漸被拉成絲狀結(jié)構(gòu)，多個(gè)顆粒脫濕損傷形成的低應(yīng)力區(qū)域互相影響，從而產(chǎn)生聚合效應(yīng)，最終導(dǎo)致了固體推進(jìn)劑失效。雖然文中沒(méi)有考慮基體材料的開裂及裂紋匯聚過(guò)程，無(wú)法展現(xiàn)固體推進(jìn)劑在失效時(shí)的斷裂，但從圖6可看出裂紋產(chǎn)生的位置及開裂時(shí)裂紋可能的擴(kuò)展路徑，如圖7所示。

為考慮復(fù)合固體推進(jìn)劑顆粒與基體之間界面損傷參數(shù)的影響，分別對(duì)表2中6組參數(shù)進(jìn)行計(jì)算。

表2 不同的界面損傷參數(shù)

2.1 界面損傷起始應(yīng)力影響分析

表2中，前2組參數(shù)考慮的是界面損傷起始應(yīng)力的影響。圖8為界面損傷起始應(yīng)力分別為0.35、0.5 MPa時(shí)顆粒夾雜模型的損傷結(jié)果。

對(duì)比圖6和圖8中的結(jié)果可看出，顆粒夾雜模型主要的損傷區(qū)域沒(méi)有發(fā)生明顯的變化，但損傷起始應(yīng)力的大小對(duì)局部區(qū)域的損傷演化速率有一定的影響。當(dāng)損傷起始應(yīng)力為0.35 MPa時(shí)，使得顆粒夾雜模型裂紋可能的擴(kuò)展路徑發(fā)生變化，如圖9所示。這主要是因?yàn)楫?dāng)損傷起始應(yīng)力為0.35 MPa時(shí)，區(qū)域A的界面損傷演化速率比區(qū)域B快，相應(yīng)的區(qū)域A內(nèi)的損傷程度較大。隨著損傷起始應(yīng)力的增大，區(qū)域B損傷演化速率逐漸增大，當(dāng)損傷起始應(yīng)力為0.5 MPa時(shí)，區(qū)域A和B內(nèi)的損傷程度基本相當(dāng)，當(dāng)損傷起始應(yīng)力為0.665 MPa時(shí)，區(qū)域B內(nèi)的損傷程度比區(qū)域A大。因此，形成了圖7和圖9所示2種可能的裂紋擴(kuò)展路徑。

在全省宣傳思想工作會(huì)議上，省委書記駱惠寧強(qiáng)調(diào)，要以習(xí)近平新時(shí)代中國(guó)特色社會(huì)主義思想為指導(dǎo)，立足新方位，擔(dān)負(fù)新使命，奮力推動(dòng)全省宣傳思想工作開創(chuàng)新局面。

圖10為復(fù)合固體推進(jìn)劑界面損傷起始應(yīng)力不同時(shí)的應(yīng)力-應(yīng)變曲線。從圖10可看出，界面損傷起始應(yīng)力越小，顆粒夾雜模型的損傷起始時(shí)間越早，損傷速率越快。因此，相應(yīng)的抗拉強(qiáng)度和最大延伸率越小。

以上現(xiàn)象可通過(guò)雙折線損傷模型來(lái)解釋，如圖11所示。從圖11可看出，損傷起始應(yīng)力較小時(shí)，對(duì)應(yīng)的界面損傷起始位移也較小，因此容易較早形成損傷，且在損傷產(chǎn)生后，相同的界面張開位移所對(duì)應(yīng)的界面剛度較小，使得在同一階段，顆粒夾雜模型的整體剛度較低。

2.2 界面初始剛度影響分析

圖12(a)、(b)分別是界面剛度為500、10 000 MPa/mm時(shí)顆粒夾雜模型的損傷形貌。對(duì)比圖6可知，界面剛度在500～15 000 MPa/mm范圍內(nèi)的改變，對(duì)顆粒夾雜模型最終的損傷形貌并沒(méi)有太大影響。

圖13為界面初始剛度不同時(shí)模型的應(yīng)力-應(yīng)變曲線。從圖13可看出，界面初始剛度為15 000、10 000 MPa/mm時(shí)的2條應(yīng)力-應(yīng)變曲線較接近，與界面初始剛度為500 MPa/mm的曲線有一定的差異。在起始階段，該曲線與另2條曲線的差異隨拉伸載荷的增大而增大，當(dāng)拉伸載荷達(dá)到一定程度時(shí)，該曲線與另2條曲線的差異又隨拉伸載荷的增大而減小，到最后階段，3條曲線逐步趨于重合。

從界面初始剛度不同時(shí)的界面損傷模型來(lái)解釋以上現(xiàn)象，如圖14所示。當(dāng)損傷起始應(yīng)力一定時(shí)，界面初始剛度的不同，使得界面損傷起始位移不同，界面初始剛度越小，界面損傷起始位移越大。

2.3 界面失效距離影響分析

圖15(a)、(b)分別是界面失效距離為0.25、0.3 mm時(shí)顆粒夾雜模型的損傷形貌。從圖15可看出，界面失效距離不同時(shí)，模型損傷的區(qū)別主要在于形成裂紋時(shí)裂紋的最大張開位移不同。界面失效距離越大，顆粒夾雜模型失效時(shí)，裂紋的張開位移越大。

圖16為顆粒夾雜模型界面失效距離不同時(shí)的應(yīng)力-應(yīng)變曲線。從圖16可看出，3條曲線在起始階段完全重合，隨著拉伸載荷的增大，區(qū)別逐漸明顯并不斷增大，界面失效距離較小的模型首先形成裂紋，界面失效距離大的模型，相應(yīng)的最大延伸率較大。

圖17為界面失效距離不同時(shí)的損傷模型示意圖。從圖17可看出，界面失效距離不同時(shí)的區(qū)別主要在于剛度衰減階段。因此，在拉伸的初始階段，曲線完全重合。隨載荷逐漸增大界面開始損傷，損傷程度越大，界面剛度的差異越大。因此，界面失效距離越小，顆粒夾雜模型裂紋產(chǎn)生時(shí)的最大張開位移及模型的最大延伸率越小。

3 結(jié)論

(1)界面損傷起始應(yīng)力的大小，在宏觀尺度上直接影響復(fù)合固體推進(jìn)劑的抗拉強(qiáng)度和最大延伸率。在細(xì)觀尺度上，影響局部區(qū)域損傷速率，導(dǎo)致裂紋可能的擴(kuò)展路徑產(chǎn)生變化。當(dāng)模型中形成的裂紋較單一時(shí)，界面損傷起始應(yīng)力越大，復(fù)合固體推進(jìn)劑的抗拉強(qiáng)度和最大延伸率越大。因此，提高顆粒與基體之間的界面強(qiáng)度，可有效改善復(fù)合固體推進(jìn)劑的力學(xué)性能。

(2)界面初始剛度在一定范圍內(nèi)的變化，對(duì)顆粒夾雜模型的損傷形貌影響較小。在宏觀尺度上，對(duì)復(fù)合固體推進(jìn)劑的抗拉強(qiáng)度和最大延伸率影響不大。

(3)界面失效距離主要影響復(fù)合固體推進(jìn)劑裂紋形成時(shí)的最大張開位移。界面失效距離越大，裂紋形成時(shí)的張開位移越大，相應(yīng)的在宏觀尺度上，復(fù)合固體推進(jìn)劑最大延伸率越大。因此，在保證界面損傷起始應(yīng)力的同時(shí)，增大固體推進(jìn)劑顆粒與基體之間界面的失效距離，可有效提升復(fù)合固體推進(jìn)劑的最大延伸率。

[1] Swanson S R,Christensen L W.A constitutive formulation for high-elongation propellants[J].Journal of Spacecraft,1983,20(6):559-566.

[2] Francis E C,Thompson R E.Nonlinear structural modeling of solid propellants[R].AIAA 84-1290.

[3] Ozupek S,Becker E B.Constitutive modeling of high-elongation solid propellants[J].Journal of Engineering Material and Technology,1992,114(1):111-115.

[4] Gazonas G A.A uniaxial nonlinear viscoelastic constitutive model with damage for M30 gun propellant[J].Mechanics of Materials,1993,15(4):323-335.

[5] Ozupek S,Becker E B.Constitutive equations for solid propellants[J].Journal of Engineering Materials and Technology,1997,119(2):125-132.

[6] Park S W,Schapery R A.A viscoelastic constitutive model for particulate composites with growing damage[J].International Journal of Solids and Structures,1997,34(8):931-947.

[7] Matous K,Inglis H M,Gu X F,et al.Multiscale damage modeling of solid propellants:theory and computational framework[R].AIAA 2005-4347.

[8] Matous K,Geubelle P H.Multiscale modeling of particle debonding in reinforced elastomers subjected to finite deformation[J].International Journal for Numerical Methods in Engineering,2006,65(2):190-223.

[9] Tan H,Liu C,Huang Y,et al.The cohesive law for the particle/matrix interfaces in high explosives[J].Journal of the Mechanics and Physics of Solids,2005,53(8):1892-1917.

[10] Tan H,Huang Y,Liu C,et al.The uniaxial tension of particulate composite materials with nonlinear interface debonding[J].International Journal of Solids and Structures,2007,44(6):1809-1822.

[11] 曲凱,張旭東,李高春.基于內(nèi)聚力界面脫粘的復(fù)合固體推進(jìn)劑力學(xué)性能研究[J].火炸藥學(xué)報(bào),2008,31(6):77-81.

[12] 李高春,邢耀國(guó),王玉峰.基于細(xì)觀力學(xué)的復(fù)合固體推進(jìn)劑模量預(yù)估方法[J].推進(jìn)技術(shù),2007,28(4):441-444.

[13] 李高春,邢耀國(guó),戢治洪,等.復(fù)合固體推進(jìn)劑細(xì)觀界面脫粘有限元分析[J].復(fù)合材料學(xué)報(bào),2011,28(3):229-235.

[14] 劉承武,陽(yáng)建紅,陳飛.改進(jìn)的Mori-Tanaka法在復(fù)合推進(jìn)劑非線界面脫粘中的應(yīng)用[J].固體火箭技術(shù),2011,34(1):67-70.

[15] 馬昌兵,強(qiáng)洪夫,武文明,等.顆粒增強(qiáng)復(fù)合材料有效彈性模量預(yù)測(cè)的多步法[J].固體力學(xué)學(xué)報(bào),2010,31(S1):12-16.

[16] Zhi Shi-jun,Sun Bing,Zhang Jian-wei.Multiscale modeling of heterogeneous propellants from particle packing to grain failure using a surface-based cohesive approach[J].Acta Mechanica Sinica,2012,28(3):746-759.

[17] 職世君,孫冰,張建偉.基于表面粘結(jié)損傷的復(fù)合固體推進(jìn)劑細(xì)觀損傷數(shù)值模擬[J].推進(jìn)技術(shù),2013,34(2):273-279.

[18] 阿蘭,達(dá)文納斯.固體火箭推進(jìn)劑技術(shù)[M].巴黎:巴黎馬森出版社,1989.

(編輯：劉紅利)

Parameter analysis of the conduplicate-line dewetting damage model of composite solid propellants

ZHI Shi-jun1，CAO Fu-qi1，SHEN Zhi-bin2，HAN Bo1，LI Ji-wei1

(1.China Airborne Missile Academy，Luoyang 471000，China；2.College of Aerospace Science and Engineering，National University of Defense Technology，Changsha 410073，China)

In order to analyze the influence of parameters of the conduplicate-line damage model on the meso-damage and macro-mechanical properties of composite solid propellants,a particle packing model of the solid propellant was created by molecular dynamics method.The surface-based cohesive approach was adopted to set the interfacial damage between particles and the bind.The packing model with different damage parameters was computed by FEM,and the numerical simulation results were compared.The results show that damage initial stress has a great effect on the tensile strength and maximum elongation of solid propellants.Interfacial initial stiffness changes in a certain extent,which affects the macro-mechanical properties and meso-damage morphology of solid propellants slightly.Interfacial failure distance has a great influence on the maximum elongation of solid propellants.

solid propellant；mesomechanics；interfacial damage；numerical simulation

2015-01-10；

2016-09-30。

國(guó)家自然科學(xué)基金(U1404106)；“十二五”總裝預(yù)研項(xiàng)目(51328050101)；國(guó)防科技大學(xué)科研計(jì)劃資助項(xiàng)目(JC13-01-03)。

職世君(1983—)，男，博士，研究方向?yàn)榘l(fā)動(dòng)機(jī)仿真設(shè)計(jì)。E-mail： zhishijun@buaa.edu.cn

V512

A

1006-2793(2017)02-0183-06

10.7673/j.issn.1006-2793.2017.02.009