2. 深海技术科学太湖实验室,江苏 无锡 214082
2. Taihu Laboratory of Deepsea Technological Science, Wuxi 214082, China
深海耐压结构可以为人员、设备或系统提供常压环境,是深海装备不可或缺的重要组成部分。圆柱壳是一种最为典型的深海耐压结构形式,其受力特性、空间利用率和加工工艺等综合性能良好,在潜艇、潜水器、深海无人系统等深海装备中被广泛采用[1 − 2]。
受静水外压的环肋圆柱壳应力计算一般通过薄壳理论来分析求解,包括2种简化力学模型:一种是将环肋视为圆柱壳梁的弹性支座,假定环肋的形心与圆柱壳壳板形心重合,肋骨的作用以其横剖面面积的方式折算到壳板中曲面上[3];另一种是将环肋圆柱壳分解为肋骨腹板、肋骨翼板和圆柱壳3个部分,分别进行受力分析并通过连接处的变形相等来进行联立求解[4]。这2种力学模型均只考虑了面内二维应力状态,忽略了静水压力载荷以及面内法向应力的影响。在传统薄壳结构范畴,这种处理可以简化工程计算且满足误差要求,但随着载荷的增加和耐压结构壳体厚度设计的增大,这种基于薄壳理论的计算结果会引起较大误差[5]。陈沙古等[6]通过肋骨作用力函数假设,提出一种考虑面内法向载荷和沿厚度方向应力影响的环肋圆柱壳三维应力分析方法,可以用于计算圆柱壳和肋骨任意位置的结构应力。
深海耐压结构尤其是单层耐压结构一般采用内置肋骨进行加强[7],本文根据三维应力分析方法导出了内置环肋圆柱壳的典型应力计算式,并通过文献中的圆柱壳算例进行检验计算。在此基础上,进一步采用三维应力计算方法对不同结构特征参数下的环肋骨圆柱壳结构进行计算,通过与采用实体单元的有限元仿真计算结果进行比较,研究分析导出的内置环肋圆柱壳三维应力计算方法的适应性。
1 三维应力计算方法 1.1 内置环肋圆柱壳力学模型假定肋骨对圆柱壳的支撑等效为内作用力函数
![]() |
图 1 内置环肋圆柱壳结构力学模型 Fig. 1 Mechanical model of inside ring-stiffened cylindrical shell |
首先定义内置环肋圆柱壳结构的一般参数:
参考文献[7]中的推导过程,首先利用弹性力学平面轴对称问题方法导出圆柱壳梁环、肋骨腹板圆环和面板圆环的三维应力和变形等式,并通过各自的力边界条件求出相关系数;再结合圆柱壳与肋骨腹板、腹板与面板连接处的尾翼边界条件,联立方程可求出之间的内作用力。经推导,圆柱壳梁环的三维应力和变形可以表达为
{σr(r,z)=A(z)r2+2C(z)+μ1−μσz(r,z),σθ(r,z)=−A(z)r2+2C(z)+μ1−μσz(r,z),σz(r,z)=m21−m2p+E(r−Reff)f1(1−μ2)[Q3coshα1zcosα2z+Q4sinhα1zsinα2z],ur(r,z)=1E[−(1+μ)A(z)r+2(1−μ)C(z)r]。 | (1) |
式中:
肋骨腹板的应力和变形可以表达为
{σrF(r)=Ri2r211−n2(−f1δ+f2δ)+(n21−n2f1δ−11−n2f2δ),σθF(r)=−Ri2r211−n2(−f1δ+f2δ)+(n21−n2f1δ−11−n2f2δ),urF(r)=1E[Ri2r(1+μ)1−n2(f1δ−f2δ)+r[(1−μ)n21−n2f1δ−1−μ1−n2f2δ]]。 | (2) |
式中:
肋骨面板的应力和变形可以表达为:
{σrM(r)=−Rfo2r211−k2f2W+k21−k2f2W,σθM(r)=Rfo2r211−k2f2W+k21−k2f2W,μrM(r)=1E[Rfo2r1+μ1−k2f2W+r(1−μ)k21−k2f2W]。 | (3) |
式中:
经推导上述方程中的作用力参数的一般解析式为:
{A(z)=−Reff2Q5[Q6p+Etf1(Q7coshα1zcosα2z+Q8sinhα1zsinα2z)],C(z)=12[Q11p−Etf1(Q23coshα1zcosα2z+Q24sinhα1zsinα2z)],f(z)=Q22p+Etf1(Q25coshα1zcosα2z+Q26sinhα1zsinα2z),f1=Q20Q12−Q13Q17−Q21σp;f2=WQ14σQ16+WQ15f1。 | (4) |
式中,相关的符号参数如下:
{α=4√3(1−μ2)√Refft;γ=√3(1−μ2)2pReff2Et2;α1=α√1−γ;α2=α√1+γ;u1=α1l2;u2=α2l2;Q1=−Reff2Etl√1−γ2u2sinhu1cosu2+u1coshu1sinu2cosh2u1−cos2u2;Q2=−Reff2Etl√1−γ2u2coshu1sinu2−u1sinhu1cosu2cosh2u1−cos2u2;Q3=Q1α12+2Q2α1α2−Q1α22;Q4=Q2α12−2Q1α1α2−Q2α22;Q5=(1−μ)Reff2Ro2+(1+μ);Q6=μm21−m2−(1−0.5μ)Refft;Q7=μQ32(1−μ2)−Q1tReff;Q8=μQ42(1−μ2)−Q2tReff;Q9=μ1−μm21−m2;Q10=μ1−μ12(1−μ2);Q11=Q6Q5Reff2Ro2−Q9−1;Q12=(1+μ)+(1−μ)n21−n2;Q13=21−n2;Q14=2n21−n2;Q15=(1−μ)+(1+μ)n21−n2;Q16=(1+μ)+(1−μ)k21−k2;Q17=WQ14δQ16+WQ15;Q18=Etδ(Q7coshu1cosu2+Q8sinhu1zsinu2);Q19=Etδ(Q23coshu1cosα2+Q24sinhu1sinu2);Q20=(1+μ)Reff2Ri2Q6Q5+(1−μ)Q11;Q21=(1+μ)Reff2Ri2Q18Q5−(1−μ)Q19;Q22=m2−1m2Reff2Ri2Q6Q5+1;Q23=Q3Q10−Q7Q5Reff2Ro2;Q24=Q4Q10−Q8Q5Reff2Ro2;Q25=2Q3Q10+m2−1m2Q7Q5Reff2Ri2;Q26=2Q4Q10+m2−1m2Q8Q5Reff2Ri2;A(0)=−Reff2Q5(Q6p+Etf1Q7);A(l/2)=−R2effQ5(Q6p+Q18f1δ);C(0)=12(Q11p−Etf1Q23);C(l/2)=12(Q11p−Q19f1δ)。 | (5) |
根据式(1)~式(4),代入坐标参数
环肋圆柱壳为轴对称结构,静水外压作用下其典型部位的应力包括相邻肋骨跨中壳板内、外表面的纵向应力和周向应力,肋骨根部跨端壳板内、外表面的纵向应力和周向应力,肋骨腹板和面板的周向应力。根据三维应力计算方法可导出内置环肋骨圆柱壳的典型应力计算式,如表1所示。
![]() |
表 1 内置环肋圆柱壳典型应力计算式 Tab.1 Stress calculation formula of inside ring-stiffened cylindrical shell |
为检验上述应力计算式,对文献[5]中的内置环肋圆柱壳进行算例验证,圆柱壳内半径为
![]() |
表 2 内置环肋圆柱壳算例计算结果 Tab.2 Calculation results of inside ring-stiffened cylindrical shell |
在静水外压作用下,跨中壳板的外表面周向应力、中面周向应力,跨端壳板的内表面纵向应力和内表面周向应力,肋骨最大周向应力和最小周向应力是内置环肋圆柱壳结构的典型应力,具体位置如图2所示。
![]() |
图 2 内置环肋圆柱壳典型位置应力示意图 Fig. 2 Typical stress of inside ring-stiffened cylindrical shell |
通过与有限元仿真计算结果比较研究上述典型应力计算式的适应性,内置环肋圆柱壳为轴对称结构,选用二维实体结构单元PLANE42建立平面轴对称模型进行仿真计算,如图3所示,在模型外部施加线分布载荷
![]() |
图 3 内置环肋圆柱壳平面实体单元模型仿真计算 Fig. 3 Plane solid element model of inside ring-stiffened cylindrical shell |
不同圆柱壳厚度特征(
![]() |
图 4 不同圆柱壳厚度特征下的典型应力计算结果比值 Fig. 4 Typical stress results of different thickness of cylindrical shell |
不同肋骨高度特征(
![]() |
图 5 不同肋骨高度特征下的典型应力计算结果比值 Fig. 5 Typical stress results of different height of ring-stiffened |
不同肋骨宽度特征(
![]() |
图 6 不同肋骨宽度特征下的典型应力计算结果比值 Fig. 6 Typical stress results of different width of ring-stiffened |
本文以内置环肋圆柱壳结构为对象,开展基于三维应力分析方法的应力计算方法适应性研究,主要工作和结论如下:
1)推导给出了内置环肋圆柱壳结构典型应力计算式,可以计算典型部位的结构应力。
2)选取不同结构特征进行计算分析,包括圆柱壳厚度特征(覆盖薄壳、中厚壳范畴),肋骨高度特征(矮肋骨、高肋骨)和肋骨宽度特征。
3)完成内置环肋圆柱壳典型应力系列计算,结果表明本文的解析计算结果与采用二维实体结构单元的有限元仿真结果较为一致,除跨端壳板的纵向应力相对误差稍大之外,其他典型应力的计算结果相对误差均很小。
研究表明,内置环肋圆柱壳三维应力计算方法不仅可以适应薄壳和中厚壳结构,且对于不同的肋骨形状也具有较好的适应性,是一种适用于深海耐压结构强度设计与评估的解析计算方法。
[1] |
徐秉汉, 朱邦俊, 欧阳吕伟, 等. 现代潜艇结构强度的理论与试验[M]. 北京: 国防工业出版社, 2007.
|
[2] |
苗兰森, 胡勇, 沈允生. 大深度救生钟耐压结构的优化设计[J]. 海洋技术, 2011, 30(1): 95-99. MIAO L S, HU Y, SHEN Y S. Stress resistant framework optimization design of deep water submarine rescue bell[J]. Ocean Technology, 2011, 30(1): 95-99. |
[3] |
吴梵, 朱锡, 梅志远. 船舶结构力学[M]. 北京: 国防工业出版社, 2010.
|
[4] |
朱邦俊, 万正权. 环肋圆柱壳应力分析的一种新方法[J]. 船舶力学, 2004, 8(4): 61−67. ZHU B J, WAN Z Q. A method of stress analysis of ring-stiffened cylindrical shell[J]. Journal of Ship Mechanics, 2004, 8(4): 61−67. |
[5] |
陈沙古, 高原, 吴智睿, 等. 大厚度环肋圆柱壳结构应力特性研究[J]. 舰船科学技术, 2023, 45(9): 19-23. CHEN S G, GAO Y, WU Z R, et al. Research on stress property of ring-stiffened cylindrical shell with large-thickness[J]. Ship Science And Technology, 2023, 45(9): 19-23. DOI:10.3404/j.issn.1672-7649.2023.09.004 |
[6] |
陈沙古, 谢晓忠, 高原, 等. 大厚度环肋圆柱壳三维应力分析方法研究[J]. 船舶力学, 2023, 27(8): 1188-1197. CHEN S G, XIE X Z, GAO Y, et al. Three dimensional stress analysis method of ring-stiffened cylindrical shell with large thickness[J]. Journal of Ship Mechanics, 2023, 27(8): 1188-1197. DOI:10.3969/j.issn.1007-7294.2023.08.008 |
[7] |
邱昌贤, 黄进浩, 秦天, 等. 内肋骨长舱段塑性总体稳定性计算方法研究[J]. 船舶力学, 2021, 25(9): 1232-1237. QIU C X, HUANG J H, QIN T, et al. Theoretical research on calculation methods for plastic general buckling stability of long cylindrical pressure hull with internal ribs[J]. Journal of Ship Mechanics, 2021, 25(9): 1232-1237. DOI:10.3969/j.issn.1007-7294.2021.09.011 |