李高春，姜愛民，黃衛(wèi)東，王朝轟，王玉峰，劉 鐵

(1. 海軍航空大學岸防兵學院，山東 煙臺 264001；2. 海軍航空大學航空基礎(chǔ)學院，山東 煙臺 264001；3. 91115部隊，浙江 舟山 316000)

引 言

對于固體火箭發(fā)動機粘接界面，由于界面兩側(cè)材料屬性差異，容易在界面附近產(chǎn)生應(yīng)力集中導致微裂紋的萌生、擴展，使得界面成為固體火箭發(fā)動機的薄弱環(huán)節(jié)之一[1-2]。在服役過程中發(fā)動機結(jié)構(gòu)完整性破壞往往發(fā)生在粘接界面上，界面失效嚴重影響發(fā)動機工作安全[3]。當前，國內(nèi)外對粘接界面開展了大量研究，包括界面粘接機理[4]、試驗表征技術(shù)[5-6]、粘接界面強度測試[7]以及界面變形過程數(shù)值模擬[8]等。李高春等[9]對粘接界面試件在拉伸過程的變形破壞過程進行了觀察，對拉伸過程圖像采用數(shù)字散斑相關(guān)方法進行處理，獲得了界面拉伸過程中的位移場，進而分析其變形演化過程及其破壞機理。史宏斌[10]考慮到絕熱層、推進劑和襯層多層材料的性能，采用有限元方法計算了含預制脫粘層藥柱的應(yīng)力應(yīng)變場，探討了脫粘層前緣附近有內(nèi)聚空洞時對界面應(yīng)力應(yīng)變的影響。Kulkarni[11]根據(jù)粘接界面內(nèi)部的組成，建立了表征界面細觀形態(tài)結(jié)構(gòu)的數(shù)值模型，分析了推進劑界面內(nèi)部顆粒與基體的粘接強度以及基體內(nèi)微裂紋等細觀結(jié)構(gòu)參數(shù)對宏觀力學性能的影響。王廣等[12]對復合固體推進劑及推進劑/襯層粘接試件微CT掃描后重構(gòu)的圖像進行了統(tǒng)計，得到了其細觀形貌特征及推進劑顆粒典型尺寸，采用粘接界面單元表征顆粒/基體和推進劑/襯層界面的脫粘過程。

盡管有不少學者對粘接界面進行了詳細的試驗和數(shù)值模擬，由于固體火箭發(fā)動機粘接界面組成的復雜性，對粘接界面細觀損傷和擴展過程以及失效機理尚不能全面把握，有待于從細觀試驗和數(shù)值模擬進一步分析其破壞過程。細觀破壞分析的最大優(yōu)點是能在細觀尺度上研究其力學性能，即從粘接界面實際組成層次研究其破壞機理，有利于掌握界面破壞的本質(zhì)機理。

本研究以固體火箭發(fā)動機粘接界面試件為對象，采用SEM對粘接界面試件在拉伸過程的變形和破壞過程進行細觀觀察，分析粘接界面試件的失效機理，根據(jù)SEM結(jié)果，建立描述界面細觀形態(tài)結(jié)構(gòu)模型，對界面細觀脫粘過程進行了數(shù)值模擬，分析了粘接界面在拉伸作用下的細觀損傷機理。研究結(jié)果有利于理解粘接界面變形破壞過程，為改善粘接界面宏觀力學性能提供必要的依據(jù)。

1 試 驗

1.1 設(shè)備及過程

為了能更加真實地再現(xiàn)拉伸過程中粘接界面變形破壞過程，采用SEM觀察了粘接界面的拉伸過程，試驗在日本島津JSM-5410LV型SEM試驗系統(tǒng)上完成，試驗裝置由高精度的掃描電鏡和全數(shù)字電液伺服加載系統(tǒng)組成，能實時在線觀測載荷作用下材料的力學行為及表面細觀結(jié)構(gòu)的變化。由于掃描電鏡試驗無法采用標準粘接界面試件，需要設(shè)計微形非標準試件?？紤]到加載系統(tǒng)的夾持方式，微形試件的尺寸如圖1所示，厚度為3.0mm。

但制作上述粘接界面試件難度較大，為了能將粘接界面觀察區(qū)移至試件中心，將微形試件兩端采用有機玻璃加工，中間平行段10mm為微形粘接界面試件。微形粘接界面試件從標準矩形試件中切取，采用302膠將兩端與有機玻璃板粘接，自然固化兩天以上，制作好的微形粘接界面試件如圖2所示。

試驗時先將加載樣品臺、試件夾具及試件表面等吹除干凈；將加載樣品臺推入電鏡真空室并按要求抽真空，加5kV電壓。在需要拍攝的倍數(shù)下，調(diào)節(jié)對比度和亮度，獲得所需圖片。試驗過程中，由伺服控制器進行自動加載直至粘接界面試件被拉斷。由于掃描電鏡成像速度較慢，為了得到各應(yīng)變條件下的清晰圖像，加載速率選定為0.12mm/min。

1.2 界面粘接破壞過程

圖3為粘接界面試件從開始拉伸到破壞過程的掃描電鏡圖像。

由圖3(a)可看出，拉伸前推進劑內(nèi)部存在一些初始缺陷，如微孔洞?；w與顆粒之間的界面粘接狀態(tài)不均勻；未施加載荷時，顆粒與基體之間存在粘接不完好的區(qū)域；由圖3(b)可看出，拉伸初期，襯層和推進劑由于彈性模量較小，在載荷作用下發(fā)生較大變形；由圖3(c)可看出，隨拉伸位移載荷的增加，在推進劑/襯層界面處的顆粒發(fā)生脫濕；界面脫濕處不再承受載荷，基體被拉伸，基體內(nèi)部產(chǎn)生應(yīng)力集中，而使其附近區(qū)域產(chǎn)生微孔洞；由圖3(d)～(e)可看出，隨著位移載荷的繼續(xù)增加，推進劑內(nèi)部顆粒與基體出現(xiàn)脫濕，并隨載荷的增加，微孔洞逐漸增大。隨載荷進一步增加，推進劑內(nèi)部微裂紋出現(xiàn)匯合，導致宏觀裂紋的形成，進而發(fā)生破壞，裂紋的擴展方向與拉伸方向近乎垂直。

2 粘接界面細觀破壞數(shù)值模擬

2.1 計算模型

為了簡便，將粘接界面看作絕熱層、襯層、推進劑組成的三相非均質(zhì)復合材料。根據(jù)試驗結(jié)果，在拉伸過程中，推進劑內(nèi)部顆粒與基體的界面以及推進劑/襯層界面是其薄弱環(huán)節(jié)，在界面處會發(fā)生脫粘，為此引入了界面元[13]。粘接界面細觀填充模型以及有限元網(wǎng)格如圖4所示，圖中只顯示粘接界面中顆粒的有限元網(wǎng)格。計算邊界條件為：絕熱層邊(即AB邊)在x、y方向均受約束作用，其位置保持不動；推進劑邊(即CD邊)沿y方向受均布位移載荷作用。

推進劑基體膠片、絕熱層和襯層視為超彈性材料，對單向拉伸數(shù)據(jù)用N=3的Ogden模型進行擬合，求出Ogden模型的系數(shù)，其結(jié)果如表1所示。

表1 基體、襯層和絕熱層的Ogden模型參數(shù)Table 1 Ogden model parameters of binder, liner and insulator

推進劑中的顆粒，由于其強度比基體大幾個數(shù)量級，可認為是彈性材料，彈性模量E和泊松比υ分別為32.4GPa和 0.14。

2.2 界面元模型

界面元模型采用一種簡單的雙線性界面元模型，界面力和位移分離量關(guān)系如圖5所示，包括材料的損傷起始(0<Δ<Δ0)、損傷擴展(Δ0<Δ<Δf)以及完全脫粘(Δf<Δ)3個階段。將n表示界面的法向，t表示切向，Δn、Δt分別表示界面的法向和切向位移分離量，L0為界面元指定的初始厚度。界面的應(yīng)變?yōu)?/p>

(1)

界面元的初始界面力為

(2)

Tn和Tt分別表示界面法向、切向力，Knn和Ktt分別表示界面的初始線性模量。

采用二次應(yīng)變率準則表示在混合模式作用下界面損傷的起始，其形式為

(3)

(4)

式中：D表示界面損傷程度，其起始值為0；當其值為1時，表示界面完全損傷。D的表達式如下：

(5)

表2 不同界面的粘接強度Table 2 Different interfacial adhesion strength

由于基體等本構(gòu)關(guān)系的材料非線性、材料大變形的幾何非線性以及界面脫粘過程的軟化現(xiàn)象，使得在計算過程中獲得收斂解變得非常困難。為了解決上述問題，采用位移控制的Newton-Raphson方法，并采用線性搜尋算法解決非線性問題。

為避免界面完全脫粘時界面元產(chǎn)生相互滲透，在界面元的上、下表面引入接觸單元來抵抗界面相互滲透。

2.3 計算結(jié)果及分析

圖6為計算模型在不同階段下的細觀損傷情況。

由圖6(a)可以看出，由于粘接界面組分材料屬性不同導致界面內(nèi)的應(yīng)力分布不均勻，并且由于顆粒與顆粒周圍的應(yīng)力場發(fā)生相互作用，產(chǎn)生較大的應(yīng)力集中。由圖6(b)可看出，隨著應(yīng)變的增加，推進劑內(nèi)部顆粒與基體的分離程度不斷增加，并在大顆粒附近形成較大的裂隙，表明上述區(qū)域容易脫濕。由圖6(c)可看出，推進劑內(nèi)部顆粒界面的脫濕將導致推進劑內(nèi)應(yīng)力的重新分布，在脫濕區(qū)域產(chǎn)生應(yīng)力集中，使其附近的顆粒產(chǎn)生脫濕。另外，由于復合推進劑顆粒填充分數(shù)比較高，顆粒之間的相互作用比較明顯，使顆粒與基體之間因脫濕形成的微孔洞大小各不相同，脫濕形成的裂隙方向與所受的應(yīng)力方向也不重合。由圖6(d)、(e)可知，隨著外界應(yīng)變的增加，推進劑內(nèi)部顆粒的脫濕分數(shù)進一步增加，造成進一步損傷。

本研究從粘接界面的細觀結(jié)構(gòu)出發(fā)，采用界面元對推進劑/襯層之間脫粘以及推進劑內(nèi)部顆粒與基體之間脫濕進行了模擬，粘接界面細觀破壞分析結(jié)果與試驗吻合，正確反映了粘接界面在拉伸過程中細觀損傷萌生與擴展規(guī)律。

在分析粘接界面的細觀破壞時，只考慮了推進劑內(nèi)部顆粒與基體之間脫濕以及推進劑/襯層之間脫粘的影響，然而粘接界面的破壞是一個復雜的過程，還存在其他不同的損傷模式，如基體的損傷。為了更加精確地模擬其細觀破壞，需要考慮這些損傷模式，從而有效地模擬其破壞的整個過程。

3 結(jié) 論

(1)在外界應(yīng)變?yōu)?%時，粘接界面表現(xiàn)為非均質(zhì)材料內(nèi)部應(yīng)力分布不均，隨應(yīng)變的增加，推進劑內(nèi)部脫濕形成的微孔洞不斷擴展，最終導致界面破壞，界面拉伸失效過程表現(xiàn)為損傷的起裂和擴展，是推進劑內(nèi)部脫濕和粘接界面脫粘共同作用的結(jié)果。

(2)針對界面的細觀破壞特點，采用界面元對推進劑/襯層之間脫粘以及推進劑內(nèi)部顆粒與基體之間脫濕進行了模擬，得到了粘接界面的脫粘過程以及其內(nèi)部應(yīng)力分布。粘接界面細觀破壞分析結(jié)果與試驗結(jié)果吻合，正確反映了粘接界面在拉伸過程中細觀損傷萌生與擴展規(guī)律。

參考文獻:

[1] Propulsion and Energetics Panel Working Group. Structural assessment of solid propellant grains , AGARD AR-350 [R].Neuilly-Sur-Seine： Advisory Group for Aerospace Research and Development, 1997.

[2] 尹華麗, 王清和. 界面粘接性能的影響因素[J]. 固體火箭技術(shù), 1998, 21(3): 40-46.

YIN Hua-li WANG Qing-he. Factors of influencing the bond characteristics at interface[J].Journal of Solid Rocket Technology, 1998, 21(3): 40-46.

[3] 邢耀國, 金廣文, 許學春, 等. 某型固體火箭發(fā)動機綜合性能試驗與壽命評估[J]. 推進技術(shù), 2004, 25(2): 176-179.

XING Yao-guo, JIN Guang-wen, XU Xue-chun, et al. Comprehensive property tests and service life evaluation for solid rocket motors [J]. Journal of Propulsion Technology, 2004, 25(2): 176-179.

[4] Grythe K F, Hansen F K. Diffusion rates and the role of diffusion in solid propellant rocket motor adhesion [J]. Journal of Applied Polymer Science, 2007, 103: 1529-1538.

[5] 邱欣, 李高春, 邢耀國. HTPB推進劑與襯層界面破壞過程試驗研究[J].兵工學報,2013,34(1):66-71.

QIU Xin, LI Gao-chun, XING Yao-guo. Experimental study on failure behavior of HTPB propellant and liner interface [J]. Acta Armamentarii,2013,34(1):66-71.

[6] 陳剛, 王小英, 張雪. 未老化NEPE推進劑/襯層粘接試件拉伸失效模式研究[J]. 固體火箭技術(shù), 2012, 35(6): 787-791.

CHEN Gang, WANG Xiao-ying, ZHANG Xue. Study on the failure mode of not aging NEPE propellant/liner bondline sample [J]. Journal of Solid Rocket Technology, 2012, 35(6): 787-791.

[7] Choupani N. Interfacial mixed-mode fracture characterization of adhesively bonded joints [J]. International Journal of Adhesion and Adhesives, 2008, 28(6): 267-282.

[8] 姜愛民, 李高春, 黃衛(wèi)東. 固體火箭發(fā)動機粘接界面力學性能的有限元計算及參數(shù)分析[J]. 火炸藥學報, 2012, 35 (4): 54-57.

JIANG Ai-min, LI Gao-chun, HUANG Wei-Dong. Finite element computation and parametric analysis on mechanical property of solid rocket motor bondline[J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2012, 35 (4): 54-57.

[9] 李高春, 劉著卿, 唐波，等. 基于數(shù)字圖像相關(guān)方法的粘接界面細觀破壞分析[J]. 宇航材料工藝, 2012, 21(3): 82-85.

LI Gao-chun, LIU Zhu-qing, TANG Bo,et al. Investigation on adhesive interface failure using digital image correlation method [J]. Aerospace Materials and Technology, 2012,21(3): 82-85.

[10] 史宏斌, 侯曉, 朱祖念, 等. 含缺陷藥柱預制脫粘層前緣應(yīng)力分析[J]. 固體火箭技術(shù), 1999, 22(3): 46-49.

SHI Hong-bin,HOU Xiao, ZHU Zu-nian,et al. Stress analysis of flap of imperfect grain in consideration of various materials [J]. Journal of Solid Rocket Technology, 1999, 22(3): 46-49.

[11] Kulkarni M G. Multiscale cohesive failure modeling of heterogeneous adhesives [D]. Urbana： University of Illinois at Urbana-Champaign, 2009.

[12] 王廣, 趙奇國, 武文明. 復合固體推進劑/襯層粘接界面細觀結(jié)構(gòu)數(shù)值建模及脫粘過程模擬[J]. 科學技術(shù)與工程, 2012, 10: 7972-7977.

WANG Guang, ZHAO Qi-guo, WU Wen-ming. Mesostructure numerical modeling and debonding procedure simulation of composite solid propellant/liner bonding interface [J]. Science Technology and Engineering, 2012, 10: 7972-7977.

[13] 李高春, 邢耀國, 戢治洪，等. 復合固體推進劑細觀界面脫粘有限元分析[J]. 復合材料學報, 2011, 28(3): 229-235.

LI Gao-chun, XING Yao-guo, ZHE Zhi-hong,et al. Finite element analysis of microscale interfacial debonding in composite solid propellants[J]. Acta Materiae Compositae Sinica, 2011, 28(3):229-235.