徐學(xué)文,牟俊林,辛慶偉
(海軍航空工程學(xué)院 接改裝訓(xùn)練大隊,煙臺 264001)
?
固體火箭發(fā)動機藥柱裂紋擴展判據(jù)的統(tǒng)計分析
徐學(xué)文,牟俊林,辛慶偉
(海軍航空工程學(xué)院 接改裝訓(xùn)練大隊,煙臺 264001)
針對固體火箭發(fā)動機藥柱裂紋擴展機理及其危險性研究的復(fù)雜性,選擇發(fā)動機工作過程中裂紋失穩(wěn)擴展作為發(fā)動機判廢的標準,采用實驗測量和數(shù)值仿真的方法,確定出了不同貯存期某型固體火箭發(fā)動機在點火發(fā)射過程中,其藥柱星角處橫向貫穿楔形裂紋發(fā)生失穩(wěn)擴展的臨界深度,應(yīng)用數(shù)理統(tǒng)計分析,給出了含該型裂紋的發(fā)動機安全工作的裂紋深度閾值曲線,從而為制訂發(fā)動機的判廢標準提供理論依據(jù)。
固體火箭發(fā)動機;裂紋擴展判據(jù);J積分;統(tǒng)計分析
固體火箭發(fā)動機在制造、運輸、貯存和勤務(wù)等過程中,因受各種載荷的作用和自身老化的影響,發(fā)動機的裝藥中可能會產(chǎn)生各種裂紋、夾雜或內(nèi)聚空洞等缺陷[1]。在這些缺陷中,藥柱裂紋最為危險。在這樣情況下,使用含有裂紋等缺陷的發(fā)動機是不可避免的。在發(fā)動機工作過程中,這些缺陷對發(fā)動機工作性能產(chǎn)生或大或小的影響,一旦裂紋發(fā)生失穩(wěn)擴展,很有可能引起發(fā)動機爆炸等災(zāi)難性后果[2-3]。因此,國內(nèi)外采用理論分析、計算機仿真和實驗測量等手段,對發(fā)動機藥柱上裂紋腔內(nèi)對流燃燒過程、裂紋擴展機理及其影響因素做了大量的研究工作,認識到裂紋擴展的過程是一個涉及到燃氣流動、傳熱、燃燒、結(jié)構(gòu)變形的強烈耦合過程,弄清了影響裂紋擴展的幾個主要因素和它們對裂紋腔內(nèi)火焰?zhèn)鞑ァ⒘鸭y擴展的定性影響[4-6]。但到目前為止,還沒有完全揭示出裂紋腔內(nèi)燃氣流動、裂紋擴展的真正規(guī)律,無法給出裂紋擴展速度與各影響因素之間的精確函數(shù)關(guān)系,不能對裂紋的危險性及發(fā)動機壽命作出準確判斷[7-8]。
根據(jù)目前國內(nèi)外研究狀況,本文將發(fā)動機工作過程中藥柱上裂紋穩(wěn)定性作為判據(jù),選擇五角星型內(nèi)孔燃燒的固體火箭發(fā)動機藥柱星角處的一條典型裂紋——橫向貫穿楔形裂紋為研究對象,根據(jù)J積分判據(jù),通過實驗測量、計算機仿真,確定出不同貯存期發(fā)動機藥柱上裂紋在點火發(fā)射過程中擴展的臨界深度;應(yīng)用數(shù)理統(tǒng)計分析,給出含該型裂紋的發(fā)動機安全工作的裂紋深度閾值曲線,可為制訂該型發(fā)動機的判廢標準提供理論依據(jù)。
該型固體火箭發(fā)動機的推進劑是貼壁澆注的高燃速HTPB復(fù)合推進劑,由于老化因素影響,發(fā)動機貯存期不同,推進劑的力學(xué)性能也不同,根據(jù)J積分判據(jù),首先通過三點彎曲斷裂實驗[9],測量出不同貯存期固體火箭發(fā)動機推進劑裂紋擴展的阻力JR值,來確定裂紋擴展量Δa=0.2 mm時的推進劑斷裂韌性,作為裂紋擴展的臨界J積分值。
本文選擇貯存2、5、8、10 a 4個時間點發(fā)動機的推進劑作為試驗樣本,通過實驗測量和計算,確定出不同貯存期推進劑不同裂紋擴展量Δa下的擴展阻力JR值,如表1所示。然后,通過線性回歸分析,確定出不同貯存年限下推進劑的斷裂韌度J0.2,將其作為裂紋擴展的臨界J積分值JC,如表2所示。
表1 不同貯存年限下的推進劑裂紋擴展阻力JRTable1 Crack propagation resistance JR of propellant with different storage years
表2 不同貯存年限下的推進劑的斷裂韌度J0.2Table2 Fraction toughness J0.2 of propellant with different storage years
本文以五角星型內(nèi)孔燃燒的固體火箭發(fā)動機為研究對象,根據(jù)固體火箭發(fā)動機復(fù)合推進劑裂紋形成的機理和無損檢測,服役期內(nèi)發(fā)動機藥柱上形成的裂紋形狀、大小和數(shù)量等千差萬別,且裂紋可能出現(xiàn)在發(fā)動機藥柱的任何位置上。這里選擇在藥柱星角處存在的一條典型裂紋——橫向貫穿楔形裂紋為研究對象,如圖1所示。裂紋縫線為直線,這主要是為了方便以后采用J積分來研究裂紋的擴展而做的假設(shè)。
首先,以發(fā)動機燃燒室和噴管一體化流場模型為研究對象,采用三維流場控制方程,應(yīng)用有限體積法計算出發(fā)動機點火發(fā)射過程中燃燒室內(nèi)流場分布[10];然后,以裂紋腔體為研究對象,將裂紋腔進口處燃燒室的壓力和溫度作為裂紋腔流場的入口邊界條件,計算出裂紋腔內(nèi)的燃氣流場分布[11];再以固體火箭發(fā)動機藥柱為研究對象,將流場計算出的燃氣壓力施加到固相計算模型中的流固耦合界面上,采用三維更新Lagrangian格式固相控制方程,應(yīng)用有限元法計算出藥柱的應(yīng)力應(yīng)變分布;最后,采用三維J積分理論計算出裂紋縫線上各積分點上的J積分值,得出了該型裂紋縫線上J積分值分布呈現(xiàn)中間高、兩端低的分布規(guī)律,裂紋擴展從縫線中間部位開始的結(jié)論[12]。根據(jù)J積分判據(jù),將裂紋縫線上最大J積分值與裂紋擴展的臨界J積分值JC相比較,從而判斷出該裂紋是否會發(fā)生擴展。本文確定出貯存2、5、8、10 a的發(fā)動機在點火發(fā)射過程中藥柱星角處橫向貫穿楔形裂紋擴展的臨界深度,如表3所示。
圖1 藥柱星角上橫向貫穿楔形裂紋的位置和形狀Fig.1 Geometry shape and position of the traversed wedge crack in the grain projection
表3 不同貯存年限下發(fā)動機藥柱星角處橫向貫穿裂紋臨界深度Table3 Critical depth of crack at the star root of motor-grain with different storage years
將不同貯存期發(fā)動機裂紋擴展的臨界深度應(yīng)用最小二乘法做線性回歸分析,得到的擬合曲線見圖2。
圖2 發(fā)動機藥柱上臨界裂紋深度回歸直線方程Fig.2 Regression equation on the critical depth of crack in motor grain
隨著發(fā)動機貯存年限的增長,藥柱上燃燒室星根處橫向貫穿楔形裂紋擴展的臨界深度也越來越小。該方程為
LC=5.882-0.209 1ta
(1)
式中LC為裂紋擴展的臨界深度,mm;ta為發(fā)動機的貯存年限,a。
采用相關(guān)系數(shù)法對回歸方程做顯著性檢驗,可知該方程在顯著水平α=0.05下,線性回歸是顯著的。根據(jù)數(shù)理統(tǒng)計分析,取置信概率為85%,則該發(fā)動機藥柱上裂紋擴展的臨界深度預(yù)測區(qū)間下限方程為
(2)
該方程即為裝藥含裂紋的發(fā)動機安全工作的裂紋深度閾值曲線方程。當發(fā)動機藥柱星角處含有的橫向貫穿裂紋深度位于該閾值曲線以下時,這種裂紋不會擴展,屬于安全型裂紋,該發(fā)動機可安全使用,且這種裂紋腔內(nèi)壁面最大面積僅占發(fā)動機燃燒室燃面面積的0.046%,裂紋腔內(nèi)的對流燃燒也不會對發(fā)動機內(nèi)彈道曲線產(chǎn)生影響。根據(jù)安全閾值曲線方程,得到該發(fā)動機不同貯存年限下藥柱上燃燒室星根處橫向貫穿裂紋安全深度,如表4所示。
表4 不同貯存年限下發(fā)動機藥柱星角處橫向貫穿裂紋安全深度Table4 Safety depth of crack at the projection of motor-grain with different storage years
(1)該型固體火箭發(fā)動機隨著貯存年限的增長,推進劑的斷裂韌性逐漸降低,發(fā)動機藥柱上的裂紋越容易失穩(wěn)擴展,在藥柱星角處的橫向貫穿楔形裂紋發(fā)生擴展的臨界深度也越來越小,發(fā)動機的工作安全性降低。
[1] 陽建紅, 俞茂宏. HTPB復(fù)合固體推進劑損傷和老化本構(gòu)研究[J]. 推進技術(shù), 2002, 23(6): 509-512.
[2] Zhang Tai-hua,L Bai Yi-long . Damage of a high-energy solid propellant and its deflagration-to-detonation transition[J]. Propellants Explosives Pyrotechnics, 2003,56(1): 37-43.
[3] Philippe H G, Hwang C Y, Robert Fiedler, et al. Simulation of dynamic fracture events in solid propellant rockets[R]. AIAA 2002-4342.
[4] Lu Y C. Combustion-induced crack propagation process in a solid-propellant crack cavity[D]. The Pennsylvania State University, 1992.
[5] 何國強, 蔡體敏, 杜磊. 含缺陷固體裝藥燃燒異常實驗分析[J]. 推進技術(shù), 1999, 20(3): 30-35.
[6] 邢耀國, 王立波, 董可海, 等. 燃燒條件下影響推進劑脫粘面擴展的因素[J]. 推進技術(shù),2001, 22(1): 77-80.
[7] Liu C T, Yang J N. Probabilistic crack growth model for application to composite solid propellants[J]. Journal of Spacecraft and Rocks, 1994, 31(1): 45-49.
[8] 袁端才, 唐國金, 李海陽, 等. 基于加速老化與三維粘彈性有限元分析的固體導(dǎo)彈發(fā)動機壽命預(yù)估[J]. 兵工學(xué)報, 2006, 27(4): 685-689.
[9] GB 2038—80. 利用JR阻力曲線確定金屬材料延性斷裂韌度的試驗方法[S]. 中華人民共和國國家標準, 1980.
[10] 徐學(xué)文, 王連生, 牟俊林. 某固體火箭發(fā)動機點火啟動過程三維流場一體化仿真[J]. 固體火箭技術(shù), 2008,31(1): 8-13.
[11] 徐學(xué)文, 王連生, 牟俊林, 等. 固體火箭發(fā)動機藥柱裂紋腔內(nèi)三維流場瞬態(tài)特性分析[J]. 固體火箭技術(shù), 2007, 30(1): 39-43.
[12] 徐學(xué)文, 孫建國, 牟俊林, 等. 某固體火箭發(fā)動機藥柱上三維裂紋擴展的判定[J]. 固體火箭技術(shù), 2008, 31(4): 331-335.
(編輯:崔賢彬)
Statistic analysis of the propagation criterion about grain crack in solid rocket motor
XU Xue-wen, MU Jun-lin, XIN Qing-wei
(Training Brigade of Equipment Acceptance and Modification, Naval Aeronautical and Astronautical University, Yantai 264001, China)
Aiming at the complexity of the propagation mechanism of the crack in solid-rocket-motor grain and its risk research, the unsTablecrack propagation was selected as the criterion of motor failure-judgment during its working. With experiment measurement and numerical simulation methods, the critical depths of unsTablepropagation about the wedge cracks which traversed through its grain projection were determined during the ignition of the solid rocket motor variation with different storage period. The threshold values of crack depths in grains during the safety working of motor were determined by use of mathematical statistics analysis, which offered the theoretical basis for the institution of motor failure-judgment criterion.
solid rocket motor;crack propagation criterion;J-integral;statistics analysis
2015-01-04;
:2015-06-16。
徐學(xué)文(1971—),男,副教授/博士,主要研究固體火箭發(fā)動機可靠性與推進劑性能。E-mail:xuxuewen02@126.com
V435
A
1006-2793(2015)05-0657-03
10.7673/j.issn.1006-2793.2015.05.010