舰船科学技术  2023, Vol. 45 Issue (9): 19-23    DOI: 10.3404/j.issn.1672-7649.2023.09.004   PDF    
大厚度环肋圆柱壳结构应力特性研究
陈沙古1,2, 高原1,2, 吴智睿1,2, 谢晓忠1,2     
1. 深海技术科学太湖实验室,江苏 无锡 214082;
2. 中国船舶科学研究中心,江苏 无锡 214082
摘要: 环肋圆柱壳是水下工程产品中广泛采用的一种典型耐压结构。随着工作深度需求越来越深,环肋圆柱壳结构设计逐渐呈现出大厚度的新结构特征。本文以某深海无人系统耐压结构参数为例,采用理论方法和有限元仿真方法进行计算分析,研究不同厚度特征参数(t/R)下的环肋圆柱壳结构应力特性。结果表明,基于薄壳理论的解析计算方法不适用于大厚度环肋圆柱壳结构应力计算,大厚度环肋圆柱壳结构的应力特性与传统圆柱壳结构存在较大差异,在大厚度环肋圆柱壳结构设计时应综合权衡其壳板应力和肋骨应力的特性。
关键词: 环肋圆柱壳     大厚度     有限元仿真     应力特性    
Research on stress property of ring-stiffened cylindrical shell with large-thickness
CHEN Sha-gu1,2, GAO Yuan1,2, WU Zhi-rui1,2, XIE Xiao-zhong1,2     
1. Taihu Laboratory of Deepsea Technological Science, Wuxi 214082, China;
2. China Ship Scientific Research Center, Wuxi 214082, China
Abstract: Ring-stiffened cylindrical shell is a typical pressure structure widely used in underwater engineering products. A new structural feature of large thickness would be gradually presented in the design of ring-stiffened cylindrical shell with the increasing demand for its working depth. Taking the pressure structure parameters of a deep-sea unmanned system as an example, series of calculation and analysis to study the stress property of ring-stiffened cylindrical shell under different thickness characteristic parameters(t/R) had been carry out using theoretical methods and finite element simulation methods in this paper. The results show that the analytical method based on thin shell theory is not suitable for stress calculation of ring-stiffened cylindrical shells with large thickness. The stress property of the ring-stiffened cylindrical shell with large thickness are quite different from those of the traditional cylindrical shell structure. The property of the shell stress and rib stress should be comprehensively weighed in the design of the ring-stiffened cylindrical shell with large thickness.
Key words: ring-stiffened cylindrical shell     large thickness     finite element simulation     stress property    
0 引 言

环肋圆柱壳是一种最为典型的水下耐压结构形式,其受力特性、空间利用率和加工工艺等综合性能良好,在潜艇、潜水器、深海无人系统等水下工程产品中被广泛采用。传统意义上,耐压圆柱壳的厚度一般远小于其半径和长度方向的尺寸,比如潜艇上运用的耐压壳体,其厚度一般不大于壳体半径的1%[1],是一种薄壳性质的概念。但随着技术的不断发展与进步,水下工程产品的工作深度指标和需求越来越深,为抵抗巨大的深水压力,保证耐压结构的安全可靠,环肋圆柱壳的壳板厚度势必越来越厚,逐渐呈现出大厚度的新结构特征。当然,这里说的大厚度仍是一个相对概念,指的是圆柱壳的厚度与壳体半径之比大于5%,远超过传统意义上的薄壳范畴。环肋圆柱壳结构应力分析的力学模型一般是将其简化为两端刚性固定在弹性支座上的复杂弯曲弹性基础梁,目前在用的潜艇、潜水器等结构设计规范[2,3]均是采用该力学模型,其实质均是基于薄壳理论假设。为进一步考虑肋骨形状和内外布置的影响,朱邦俊[4]提出一种将环肋圆柱壳分解为肋骨腹板(二维环板)、肋骨翼板(一维圆环)和圆柱壳3个部分进行联立求解的理论解析方法。针对具有大厚度新结构特点的环肋圆柱壳结构应力计算,目前还没有成熟统一的规范标准。

本文以某深海无人系统耐压结构参数为例,采用文献[4]提出的解析方法和有限元仿真方法,研究不同厚度特征参数(t/R)下的环肋圆柱壳结构应力特性,研究结果可为深海无人系统耐压结构设计和强度评估提供技术支撑。

1 环肋圆柱壳结构

在圆柱壳外部或内部间隔设置环形肋骨是为了提高圆柱壳的结构强度与稳定性,肋骨一般采用T型结构,即包括肋骨腹板、肋骨面板两部分。对于等间距布置肋骨的环肋圆柱壳,其结构参数主要包括:圆柱壳内半径 $ R $ ,壳板厚度 $ t $ ;肋骨的腹板高度 $ H $ 、厚度 $ \delta $ ,面板宽度 $ W $ 、厚度 $ b $ ;相邻肋骨的间距 $ l $ ,其结构及对应参数符号如图1所示。

图 1 环肋圆柱壳结构及参数符号 Fig. 1 Structure parameter of ring-stiffened cylindrical shell

以某深海无人系统耐压结构参数[5]为例进行研究,其主要结构参数如表1所示。

表 1 环肋圆柱壳算例参数表 Tab.1 Parameters table of ring-stiffened cylindrical shell
2 应力计算方法 2.1 文献[4]计算方法

文献[4]将环肋圆柱壳分解为肋骨腹板、肋骨翼板和圆柱壳3个部分,以外置肋骨为例,其力学模型如图2所示。圆柱壳简化为承受静水压力和纵向压力 $ T $ 联合作用的弹性基础梁,肋骨腹板简化为承受 $ {f_1} $ $ {f_0} $ 作用的平面应力二维环板,肋骨面板简化为承受 $ {f_0} $ 作用的一维圆环。计算分析的思路是将这3部分之间的作用力( $ {f_1} $ $ {f_0} $ )作为载荷输入分别独立进行解析求解,再结合三者之间的变形协调关系求解作用力参数。

图 2 外置环肋圆柱壳结构力学模型 Fig. 2 Mechanical model of external ring-stiffened cylindrical shell

由于圆柱壳壳内最大应力可能发生在相邻肋骨的跨中处(此处径向位移最大),也可能发生在肋骨跨端处[6](该处壳上的弯曲应力为最大值)。文献[4]进一步给出了环肋圆柱壳结构典型部位的应力计算式。

1)肋骨根部的耐压壳板纵向应力

$ {\sigma _1}\left( {\frac{l}{2}} \right) = - \frac{{pR}}{{2t}} \pm \frac{{6D}}{{{t^2}}}w''\left( {\frac{l}{2}} \right)。$ (1)

式中,“±”符号,上符号为内表面,下符号为外表面。

2)肋骨根部的耐压壳板周向应力

$ {\sigma _2}\left( {\frac{l}{2}} \right) = - \frac{{(1 + 0.5\mu )pR}}{t} \pm \frac{{6\mu D}}{{{t^2}}}w''\left( {\frac{l}{2}} \right)。$ (2)

式中,“±”符号,上符号为内表面,下符号为外表面。

3)跨中壳板纵向应力

$ {\sigma _1}\left( 0 \right) = - \frac{{pR}}{{2t}} \pm \frac{{6D}}{{{t^2}}}w''\left( 0 \right)。$ (3)

式中,“±”符号,上符号为内表面,下符号为外表面。

4) 跨中壳板周向应力

$ {\sigma _2}\left( 0 \right) = - \frac{{(1 + 0.5\mu )pR}}{t} \pm \frac{{6\mu D}}{{{t^2}}}w''\left( 0 \right)。$ (4)

式中,“±”符号,上符号为内表面,下符号为外表面。

5)肋骨面板应力

$ {\sigma _f} = - \frac{{{f_0}{R_f}}}{{{F_f}}}。$ (5)

6)肋骨腹板周向应力

$ \left\{ \begin{gathered} {\sigma _\theta }(R) = \frac{{1 + {m^2}}}{{1 - {m^2}}}\frac{{{f_1}}}{\delta } - \frac{{2{m^2}}}{{1 - {m^2}}}\frac{{{f_0}}}{\delta },\\ {\sigma _\theta }({R_f}) = \frac{2}{{1 - {m^2}}}\frac{{{f_1}}}{\delta } - \frac{{1 + {m^2}}}{{1 - {m^2}}}\frac{{{f_0}}}{\delta }。\\ \end{gathered} \right. $ (6)

式中,肋骨腹板上、下端圆弧半径之比 $ m = {R_f}/R $

2.2 有限元仿真方法

环肋圆柱壳为轴对称结构,可采用平面实体单元建立平面轴对称模型进行仿真计算分析,有限元模型如图3所示。

图 3 环肋圆柱壳平面实体单元模型 Fig. 3 Finite element model of ring-stiffened cylindrical shell

在外压载荷作用下,环肋圆柱壳的外表面(包括肋骨部分)直接承受深水压力,施加压力 $ p $ ;轴向上,环肋圆柱壳模型一端施加刚体位移约束,一端施加轴向等效压力:

$ {p_1} = p{(1 + \eta )^2}/\left[ {(2 + \eta )\eta } \right]。$ (7)

$ p $ =45 MPa, $ \eta $ =0.07时, $ {p_1} $ =355.56 MPa,有限元模型边界条件施加如图4所示,其应力仿真云图如图5所示,根据仿真计算结果可以准确提取环肋圆柱壳结构各部位的应力。

图 4 环肋圆柱壳模型边界条件施加方式 Fig. 4 Mode of applying boundary conditions to the ring-stiffened cylindrical shell model

图 5 环肋圆柱壳结构应力仿真云图 Fig. 5 Stress nephogram of ring-stiffened cylindrical shell structure

将某环肋骨圆柱壳结构( $ \eta $ =0.028)的有限元仿真与模型试验结果[7]进行比较,两者关于跨中壳板、跨端壳板、肋骨面板等应力的结果吻合较好(见表2),表明采用平面实体单元模型能较准确仿真计算较厚环肋圆柱壳结构的应力。

表 2 某环肋圆柱壳结构应力的仿真计算值与试验值 Tab.2 Simulation and test results of a ring-stiffened cylindrical shell
3 应用算例与分析

表1某深海无人系统耐压结构参数为应用算例,对不同结构特征参数( $ \eta = t/R $ )下的环肋骨圆柱壳进行系列计算,圆柱壳壳板厚度分别取为15 mm,20 mm,25 mm,30 mm,35 mm,40 mm,45 mm,50 mm,对应的结构特征参数 $ \eta $ 介于0.03~0.1之间。

文献[4]和有限元仿真2种应力计算方法的计算结果如表3表5所示。主要给出了相邻肋骨跨中处壳板的周向应力、纵向应力,肋骨根部跨端壳板的周向应力、纵向应力,以及肋骨面板应力,腹板上、下圆弧处的周向应力。

表 3 跨中壳板应力计算结果 Tab.3 Stress results of mid-span shell plate

表 4 跨端壳板应力计算结果 Tab.4 Stress results of cross-end shell plate

表 5 肋骨应力计算结果 Tab.5 Stress results of rib

分析表中数据,可以看出:

1)与有限元仿真结果相比,理论解析方法关于耐压壳板应力的计算结果偏大,跨中壳板应力解析计算值偏大约30%,跨端壳板应力解析计算值偏大超过30%;关于肋骨应力的计算结果偏小,解析计算值偏小约15%,具体如表6所示。

表 6 两种应力计算方法之间的误差 Tab.6 Error between two methods of stress calculation

2)跨中壳板的周向应力大于其纵向应力,如图6所示。外表面周向应力稍大于内表面周向应力,且随着圆柱壳厚度的增加,两者逐渐趋近,表明大厚度环肋圆柱壳跨中壳板中面周向应力可以反映该处的应力特性。

图 6 跨中壳板应力随厚度特征参数的变化曲线 Fig. 6 Stress curve of mid-span shell plate with thickness characteristic parameter

3)由于肋骨的加强作用引起局部弯曲应力,跨端壳板内表面应力大于外表面应力,如图7所示。随着圆柱壳厚度的增加,内表面纵向应力会出现小于其周向应力的情况,与传统的环肋圆柱壳跨端壳板应力特性存在较大的差异。

图 7 跨端壳板应力随厚度特征参数的变化曲线 Fig. 7 Stress curve of cross-end shell plate with thickness characteristic parameter

4)肋骨面板应力小于肋骨腹板应力。为分析圆柱壳厚度增加对结构典型应力水平的改善效益,对其进行归一化比较,如图8所示。随着圆柱壳厚度的增加,壳板应力将急剧降低,肋骨应力水平也会降低,但改善效率不如壳板应力。

图 8 结构典型应力比随厚度特征参数的变化曲线 Fig. 8 Curve of stress ratio with thickness characteristic parameter
4 结 语

针对大厚度环肋圆柱壳结构,本文以某深海无人系统耐压结构参数为例,采用文献解析方法和有限元仿真方法对不同厚度特征参数下的环肋圆柱壳结构进行计算和比较分析,研究结论如下:

1)基于薄壳理论的解析计算方法不适用于大厚度环肋圆柱壳结构应力计算,其关于壳板应力和肋骨应力计算结果都存在较大的误差;

2)与传统环肋骨圆柱壳类似,大厚度环肋圆柱壳的跨中壳板应力以周向应力为主,跨中壳板中面周向应力能够反映该处的应力特性;

3)大厚度环肋圆柱壳的跨端壳板内表面纵向应力会出现小于其周向应力的情况,与传统环肋圆柱壳跨端壳板应力特性存在较大的差异;

4)随着圆柱壳厚度的增加,圆柱壳壳板应力的改善效益远高于肋骨应力,在大厚度环肋圆柱壳结构设计时,应综合权衡其壳板应力和肋骨应力的特性。

上述研究结果可以为采用大厚度环肋圆柱壳的某深海无人系统耐压结构设计提供基础支撑。

参考文献
[1]
谢祚水, 王自力, 吴剑国. 潜艇结构分析[M]. 武汉: 华中科技大学出版社, 2004.
[2]
中华人民共和国国家军用标准. GJB/Z21A-2001潜艇结构设计计算方法[S]. 北京: 国防科学技术工业委员会, 2001.
[3]
中国船级社. 潜水系统和潜水器入级规范 2013 综合文本[S]. 北京, 2014.
[4]
朱邦俊, 万正权. 环肋圆柱壳应力分析的一种新方法[J]. 船舶力学, 2004, 8(4): 61−67.
ZHU Bangjun, WAN Zhengquan. A method for stress analysis of ring-stiffened cylindrical shell[J]. Journal of Ship Mechanics, 2004, 8(4): 61−67.
[5]
高原, 陈沙古, 吴智睿. 某结构样机设计方案报告[R]. 中国船舶科学研究中心科技报告, 2022.
[6]
黄进浩,邱昌贤,蔡新钢,等. 主要潜水器规范圆柱形耐压壳体设计计算方法比较分析[J]. 船舶力学,2017(z1):351-362.
HUANG Jinhao, QIU Changxian, CAI Xingang,et al. Comparative analysis on cylindrical structure design calculation method by main submersible rules[J]. Journal of Ship Mechanics, 2017(z1): 351−362.
[7]
陈沙古, 高原, 陆波, 等. 结构安全性验证缩比模型试验报告[R]. 中国船舶科学研究中心科技报告, 2020.