潛艇發射筒保護蓋力學性能分析

焦冀光，楊春，劉文一

?

潛艇發射筒保護蓋力學性能分析

焦冀光1，楊春2，劉文一1

（1.91550部隊91分隊，遼寧 大連 116023；2.北京宇航系統工程研究所，北京 100076）

目的研究艇載發射筒保護蓋在深水壓力和溫度載荷下的力學性能。方法建立保護蓋力學載荷和溫度載荷模型，采用有限元技術仿真計算保護蓋在這些載荷下的形變、應力、溫度分布和熱結構情況。結果深水壓力載荷作用下，保護蓋出現的形變大小在0.977~1.47 cm范圍內，應力大小在1.22~7.91 MPa范圍內。溫度載荷下出現的形變大小在0.92~1.38 cm范圍內，應力大小在0.272~1.26 MPa范圍內。結論兩種工況下保護蓋最大形變均出現在中部，因此這個部位結構強度需加強，以防其被破壞。

發射筒；潛艇；保護蓋

某潛艇發射筒保護蓋由非金屬材料制成，可以對發射筒內的雷彈起到保護作用。潛艇在水下執行巡航任務時，水深產生的壓力會對保護蓋結構產生靜壓影響；在水面上時，筒內外溫度差會對保護蓋的結構產生影響[1]。為了研究這些載荷對非金屬保護蓋的影響，文中建立了溫度變化模型和熱結構模型，采用有限元技術計算了保護蓋在這些載荷下的應力、變形和溫度分布情況。

1 保護蓋模型

保護蓋為圓柱形殼體，整體制造，用夾緊環夾緊保護蓋下端面法蘭，固定在發射筒上，且法蘭上下端面有密封圈以保證發射筒的氣、水密。圖1為保護蓋的三維模型和有限元模型，使用IsoMesh劃分器，劃分為Quad4型單元。這樣劃分的原因是便于控制網格疏密，共劃分了25 994個單元[2]。

圖1 保護蓋的及其有限元模型

保護蓋為非金屬材料，其性能參數：密度為1730 kg/m3，泊松比為0.495，彈性模量為290 MPa，導熱系數為0.57 W/(m·K)，熱膨脹系數為3.2×10－6K－1。

2 保護蓋內外壓差對保護蓋的影響

為了滿足裝載在發射筒內的雷彈貯存需要，發射筒內充有一定壓力的氣體，發射筒位于水下時，水深會在發射筒保護蓋內外產生一定的壓差，若保護蓋承壓能力差，則壓差有可能使保護蓋失穩或破壞保護蓋，使水進入發射筒，進而使雷彈失效[3]。

2.1 保護蓋穩定性分析模型

在有限元方法中，線形穩定性理論通過提取使線性系統剛度矩陣奇異的特征值來獲得結構的失穩載荷，即失穩模態。線性穩定性分析的特點是變形為小變形，單元應力必須是彈性的，載荷引起的內部單元力分布是不變的[4]。離散后有限元模型的屈曲分析控制方程為：

式中：為剛度矩陣；d為剛度微分矩陣；為特征值；為特征向量。

這個方程可簡化為一種標準形式求解。采用Lanczos 法對其實特征值求解，將式（1）改寫為[5—6]：

選擇初始迭代矢量：

(3)

依次取= 1，2，3，…，直到出現收斂的迭代根，迭代函數為：

其中：

(5)

將α，β合并為一個變換矩陣，如式（6）：

設為特征方程=的特征值根，則式（2）的特征值為：

(7)

因此為使線性系統剛度矩陣奇異的特征值，即為所求保護蓋的的臨界壓力[7]。

2.2 保護蓋靜力分析模型

對于結構材料，可以由Kirchoff應力張量t和Green應變張量γ來描述其應力應變狀態，它們滿足以下關系式：

式中：1，2，3為變形張量不變量；為應變能密度函數[6]。根據計算所用C/C材料模型，可以推導出其主應力1和主伸長比1之間的關系為[8]：

(9)

2.3 仿真計算及結果分析

計算了發射筒在水下時，保護蓋在水靜壓下的受力情況，得到了形變和應力，如圖2所示。

圖2 保護蓋在水壓載荷下的形變與應力

從計算結果來看，在深水靜壓的作用下，保護蓋出現了較大的形變和應力，形變大小在0.977~1.47 cm范圍內，形變最大的地方在保護蓋幾何中心。應力大小在1.22~7.91 MPa范圍內，應力最大的地方在保護蓋中部和邊緣位置，在邊緣位置出現了應力集中現象。

3 溫度對保護蓋的影響

3.1 熱結構分析模型

當結構的內部存在溫度梯度的時候，溫度的變化必然會導致結構產生膨脹或者收縮的現象。當結構的外部存在某些約束的時候，這些變形必然會受到抑制，此時結構內部就會存在某種應力來抵抗這些變形，這種由于結構的溫度場引起的應力叫做熱應力[9]。

當結構的內部存在熱應力時，結構的總應力分為兩部分：結構機械應力以及熱應力[10]。由于結構的小變形假設，因此滿足線性疊加原理，那么此時結構的應力可以表示為[11]：

對于結構的應變可以寫成：

(11)

式中：σ為結構的應力分量；ε為結構的應變分量；c為結構的剛度系數；s為結構的柔度系數；β為結構的熱模量；α為結構的熱膨脹系數[12]。

3.2 仿真計算及結果分析

潛艇在巡航時，在某些海域發射筒保護蓋外表面受到陽光照射而有較高的溫度，保護蓋內表面則由于雷彈貯存溫度需要而保持一定的溫度。這樣會在保護蓋外表面產生溫度差，計算得到了保護蓋的溫度分布和溫度梯度分布，如圖3所示。

從計算結果來看，在保護蓋下圓弧處，溫度變化較大，因而出現了較大的溫度梯度。保護蓋在溫度載荷作用下，會出現變形，計算得到了保護蓋在溫度載荷下的形變和應力，如圖4所示。

圖3 保護蓋溫度與溫度梯度

從計算結果來看，保護蓋在溫度載荷下出現了較大的形變，形變大小在0.92~1.38 cm范圍內，形變最大的地方在保護蓋幾何中心。應力大小在2.72 kPa~1.26 MPa范圍內，應力最大的地方在保護蓋根部位置，在根部位置出現了應力集中現象。

4 結論

1）深水靜壓載荷作用下，保護蓋出現了較大的形變和應力，形變大小在0.977~1.47 cm范圍內，形變最大的地方在保護蓋幾何中心。應力大小在1.22~7.91 MPa范圍內，應力最大的地方在保護蓋中部和邊緣位置，在邊緣位置出現了應力集中現象。

2）太陽直射時，保護蓋內外溫差較大。在保護蓋下圓弧處，溫度變化較大，因而出現了較大的溫度梯度。

3）保護蓋在溫度載荷下出現了較大的形變，形變大小在0.92~1.38 cm范圍內，形變最大的地方在保護蓋幾何中心。應力大小在272 kPa~1.26 MPa范圍內，應力最大的地方在保護蓋根部位置，在根部位置出現了應力集中現象。

4）水壓和溫度載荷作用下，保護蓋最大形變均在保護蓋中部，因此這個部位需加強，以防其被破壞。

[1] 劉文一, 焦冀光. 固體發動機裝藥熱安全性分析[J]. 裝備環境工程, 2016, 13(2): 129—131.

[2] MIL-STD-810G, Department of Defense Test Method Standard—Environmental Engineering Considerations and Laboratory Tests[S].

[3] 黃樹濤, 周旺, 宋博, 等. 扭轉載荷作用下熔接涂層與基體結合界面處應力特性仿真研究[J]. 表面技術, 2017, 46(1): 57—58.

[4] 高帥, 曹磊, 張泉, 等. 金屬鋁表面超疏水薄膜的構筑及減摩特性[J]. 表面技術, 2016, 45(1):76—78.

[5] 胡建朋, 劉智勇, 胡山山, 等. 304 不銹鋼在模擬深海和淺海環境中的應力腐蝕行為[J]. 表面技術, 2015, 44(3): 9—11.

[6] 趙建忠, 葉文, 田建海. 艦載環境對導彈武器裝備可靠性的影響分析及對策[J]. 質量與可靠性, 2014(2): 5—9.

[7] 劉文一, 焦冀光. 艦載環境對固體發動機裝藥影響分析[J]. 彈箭與制導學報, 2015, 35(6): 80—82.

[8] 林琳, 張熙川, 葉濤. MIL-STD-810F 低溫試驗方法研究[J]. 電子產品可靠性與環境試驗, 2010, 28(2): 5—8.

[9] 劉文一, 焦冀光. 某飛行器復合材料薄壁加筋結構艙段穩定性分析[J]. 艦船電子境工程, 2016, 36(2): 110—113.

[10] DEF STAN 00-35, General Specification for Aircraft Gas Turbine Engines[S].

[11] 唐荻, 趙愛民, 武會賓, 等. 板帶鋼軋制新技術及品種研發進展[J]. 鋼鐵, 2012, 47(11): 1—47.

[12] 石文天, 劉玉德, 張永安, 等. 芳綸纖維復合材料切削加工研究進展[J]. 表面技術, 2016, 45(1): 31—32.

Mechanical Property of Submarine Launch Canister Protection Cover

JIAO Ji-guang1, YANG Chun2, LIU Wen-yi2

(1.Unit 91of PLA 91550, Dalian 116023, China；2.Beijing Institute of Astronautical System Engineering, Beijing 100076, China)

Objective To study mechanical properties of submarine launch canister protection cover under deepwater pressure and temperature load. Methods Model for mechanical and temperature load of protection cover was established. Finite-element method was used to analyze the stress, deformation, temperature distribution and thermal structure under those loads. Results Under deep water pressure, the deformation was 0.977 cm to 1.47 cm, the stress was 1.22 MPa ~ 7.91 MPa. Under temperature load, the deformation of the protection cover is 0.92~1.38 cm and the stress was 0.272 kPa ~ 1.26 MPa. Conclusion The maximum deformation of protection cover was on centre portion in two working conditions, so the structural strength of centre portion should be reinforced, otherwise it would be damaged.

launch canister; submarine; protection cover

10.7643/ issn.1672-9242.2017.04.018

TJ07；TG174

A

1672-9242(2017)04-0087-04

2017-01-20；

2017-02-23

焦冀光（1984—），女，湖北武漢人，主要研究方向為武器控制及優化設計。