閱讀文章前辛苦您點下“關注”，友善讨論和分享，為了回饋您的支援，我将每日更新優質内容。
前言
脈沖發動機利用隔離裝置将固體火箭發動機的燃燒室分成若幹部分，可進行多次的關機和啟動，合理配置設定推力及各脈沖間隔時間，實作飛彈飛行彈道的最優控制和發動機能量的最優管理。
采用脈沖發動機的飛彈具有高速、遠射程、高命中率和較輕的系統結構品質等特點，在實際應用中，雙脈沖固體火箭發動機不僅能明顯提升性能，又不會過度提升設計複雜性和大幅降低系統可靠性。
雙脈沖發動機燃燒室内2個藥柱由隔離裝置間隔開。在一脈沖工作時，隔離裝置起到隔離高溫瓦斯的作用，避免二脈沖藥柱受熱燃燒，此外隔離裝置還要承受燃燒室内高壓瓦斯作用，保持自身結構完整性。
對于軟隔層結構的雙脈沖發動機，一脈沖工作時，軟隔層受到高溫高壓瓦斯作用，會産生較大的變形，與二脈沖藥柱接觸貼合。
若隔層結構設計不合理會導緻結構破壞，隔離失效，造成二脈沖藥柱提前燃燒工作，燃燒室内壓強快速升高，最終導緻發動機殼體超壓破壞，工作失利。
本文主要研究對象為雙脈沖發動機軟隔層，根據雙脈沖發動機工作過程中徑向軟隔層大變形、高彈性、小應變的特點，以實測的橡膠材料、推進劑性能資料為基礎，采用超彈性本構模型對隔層三元乙丙橡膠材料以及推進劑材料性能進行表征。
對一脈沖工作時二脈沖藥柱、隔層間隙工況下的軟隔層受力情況，應用Abaqus軟體進行數值模拟，獲得隔層的實際應變情況，結合實際發動機研制基礎，為後續雙脈沖發動機軟隔層設計及優化指出方向。
材料參數确定
本文中雙脈沖發動機所用隔離裝置為橡膠軟隔層，在加載過程中具有大變形、高彈性的性質。對于橡膠材料，常用的本構模型有Mooney-Rivlin模型和Odgen模型。
Mooney-Rivlin模型能很好地展現橡膠材料在中等應變情況下方位應變能，由于結構簡單，在超彈性本構計算中被廣泛應用，其應變能函數如下：
式中，c、d為材料常數，由材料試驗确定；I1、I2為變形張量；J為體積率，橡膠是不可壓縮的，故取J=1。通過應變能函數對主伸長比λ求導，并根據不可壓縮材料單軸拉伸條件可得：
式中，t為單軸拉伸應力；λ為主伸長比：
Odgen模型由Odgen于1972年提出，是将主伸長比λi作為應變能函數自變量的新的應變能函數。Odgen應變能函數定義為：
式中，μn和αn為材料常數；Ｋ為初始體積模量。過應變能函數對主伸長比求導，并根據不可壓縮材料單軸拉伸條件可得：
進一步可得基于Odgen模型的單軸拉伸應力表達式為：
根據上述方式，從業人員分别采用M-R模型和Odgen模型對橡膠标準試件500mm/min拉伸速率單項拉伸試驗資料進行拟合，超彈性本構模型力與位移曲線如圖1所示，本構特性參數見表1和表2。
由圖1可看出，相比Mooney-Rivlin模型，Odgen模型與材料拉伸曲線貼合度更好，能更準确地反應隔層彈性材料性能，且随着參數數量的提升，4階Odgen模型比3階模型具有更高的精度，是以本文中采用4階Odgen模型進行隔層結構受力分析。
材料變形的應變能（或變形能）是在外載荷的作用下發生變形，研究表明，儲存在材料中的能量，應變能函數能較好表征物體的變形特征。
在超彈性材料中，常用的應變能模型有描述橡膠變形行為的Mooney-Rivlin模型、Odgen模型等及描述生物薄膜變形行為的Fung模型、GPR模型等。
固體推進劑是典型的粘超彈性材料，在許多實際應用中，其力學行為表現出明顯的非線性。為了準确表征推進劑在有限變形條件下的力學性能，對推進劑的本構模型進行研究。
由于推進劑在單向拉伸條件下應力應變曲線與橡膠類材料拉伸響應曲線形狀相近，本文考慮以橡膠類材料的應變能表達式作為函數模型來描述推進劑的拉伸變形。固體推進劑單向拉伸應力應變曲線如圖2所示，拟合推進劑本構特性參數見表3。
計算執行個體
為了簡化計算，提高計算速度與精度，本次計算采用二維軸對稱模型，僅對二脈沖側結構進行分析，整個模型由燃燒室殼體、二脈沖藥柱、絕熱層、隔層和頂蓋組成，重點研究隔層在一脈沖工作過程中承載變形的應力應變。
為減輕沙漏效應，隔層與藥柱采用軸對稱雜交縮減積分網格CAX4RH，其餘結構采用軸對稱縮減積分網格CAX4R。
發動機直徑Ｒ，二脈沖藥柱長度3.5R，隔層厚度取0.06R，對隔層與藥柱取2、4和6mm的不同徑向間隙進行分析，發動機軸對稱模型如圖3所示。
在雙脈沖發動機一脈沖工作時，隔層承受一脈沖工作高壓瓦斯作用，會産生較大變形，與二脈沖藥柱貼合。由于隔層頭部與頂蓋連接配接、大端與殼體粘接，變形受到限制，理論分析隔層承載會在筒段産生較大的環向和徑向應變，軸向應變相對較小，是以分析主要針對隔層的筒段部位。
隔層變形後與藥柱接觸，存在一定的摩擦因數，本次計算中對隔層與藥柱接觸摩擦因數取0.15。發動機前裙施加強定限制，對隔層與殼體粘接處采用綁定限制。
在分析過程中對隔層施加法向的壓力，載荷采用3個分析步進行施加，以利于隔層和藥柱建立平穩接觸，加載曲線如圖4所示。
圖5所示為内壓載荷下隔層整體徑向應變，由圖5可知，隔層大端承受拉應力，應變最大值為2.02%，位于錐段與筒段轉折位置；隔層筒段整體承受壓應力，應變最大值為15.44%，位于隔層頭部，筒段中部位置應變值為12.62%。
圖6所示為内壓載荷下隔層整體軸向應變，由圖6可知，隔層整體應變水準較低，在隔層錐段與筒段轉折位置應變水準較大，既存在拉應力也存在壓應力，拉應力最大應變為10.61%，壓應力最大應變為7.64%，筒段中部位置應變值為0.63%。
圖7所示為内壓載荷下隔層整體環向應變，由圖7可知，隔層筒段承受拉應力，應變最大值為13.48%，位于筒段靠近頭部位置，筒段中部位置應變值為10.87%。
圖8所示為内壓載荷下隔層整體徑向應變，由圖8可知，隔層大端承受拉應力，應變最大值為5.07%，位于錐段與筒段轉折位置；隔層筒段整體承受壓應力，應變最大值為23.39%，位于隔層頭部，筒段中部位置應變值為17.21%。
圖9所示為内壓載荷下隔層整體軸向應變，由圖9可知，隔層整體應變水準較低，在隔層錐段與筒段轉折位置應變水準較大，既存在拉應力也存在壓應力，拉應力最大應變為17.16%，壓應力最大應變為13.5%，筒段中部位置應變值為1.12%。
圖10所示為内壓載荷下隔層整體環向應變，由圖10可知，隔層筒段承受拉應力，應變最大值為19.16%，位于筒段靠近頭部位置，筒段中部位置應變值為16.94%。
圖11所示為内壓載荷下隔層整體徑向應變，由圖11可知，隔層大端承受拉應力，應變最大值為8.24%，位于錐段與筒段轉折位置；隔層筒段整體承受壓應力，應變最大值為32.0%，位于隔層頭部，筒段中部位置應變值為22.4%。
圖12所示為内壓載荷下隔層整體軸向應變，由圖12可知，隔層整體應變水準較低，在隔層錐段與筒段轉折位置應變水準較大，既存在拉應力也存在壓應力，拉應力最大應變為20.81%，壓應力最大應變為18.60%，筒段中部位置應變值為1.85%。
圖13所示為内壓載荷下隔層整體環向應變，由圖13可知，隔層筒段承受拉應力，應變最大值為24.72%，位于筒段靠近頭部位置，筒段中部位置應變值為21.88%。
将隔層與藥柱徑向3種間隙距離獲得應變結果進行對比分析，3個方向應變值見表4。通過對隔層承載應力應變進行數值模拟并進行結果對比，可以得出如下結論。
隔層錐段與筒段轉折處、隔層筒段頭部應變水準較高，對于徑向軟隔層雙脈沖發動機，該位置是發動機設計關鍵位置，需要優化結構以降低整體應變水準。
徑向軟隔層軸向整體應變水準較低；徑向方向上隔層與藥柱互相擠壓，主要承受壓應力，對應為壓縮應變。隔層整體受力水準與隔層、藥柱間隙呈正相關，間隙越小，隔層受力水準越小。
隔層在一脈沖工作時承載應變水準遠小于材料的極限延伸率（≥300%），與實際試驗中隔層破壞現象存在沖突，說明隔層在達到強度極限前已經失效。
結語
通過上述研究可以得出如下結論。通過對一脈沖工作時隔層受力情況進行有限元分析，認為隔層筒段頭部以及錐段與筒段轉折處是結構的薄弱點，需要進行優化設計以降低隔層應變水準。
通過有限元分析，并結合實際試驗中出現問題，說明隔層在應變遠小于材料的極限時已經失效，表明現有的以最大延伸率作為材料承載極限的安全判據不适用。
分析不同二脈沖藥柱、隔層間隙下隔層的應變變化情況，發現隔層受力水準與間隙呈正相關，減小藥柱與隔層間隙可降低隔層應變水準。