舰船科学技术  2024, Vol. 46 Issue (24): 15-19    DOI: 10.3404/j.issn.1672-7649.2024.24.003   PDF    
内置环肋圆柱壳三维应力计算方法适应性分析
陈沙古1,2, 高原1,2, 张平平1,2, 周平1,2     
1. 中国船舶科学研究中心,江苏 无锡 214082;
2. 深海技术科学太湖实验室,江苏 无锡 214082
摘要: 环肋圆柱壳三维应力分析方法计及面内法向载荷和应力的影响,可计算圆柱壳和肋骨任意位置的结构应力。本文根据三维应力分析方法导出内置环肋圆柱壳结构的典型应力计算式,并通过文献中的算例检验应力计算式的正确性。采用典型应力计算式对不同结构特征的内置环肋圆柱壳进行系列计算,通过与有限元仿真计算结果比较,研究分析典型应力计算式对不同圆柱壳厚度特征(t/R)、肋骨高度特征(H/R)、肋骨宽度特征(W/l)等参数的适应性。结果表明,基于三维应力分析方法的内置环肋圆柱壳典型应力计算式不仅可以适应薄壳和中厚壳结构应力计算,且对于不同的肋骨形状也具有较好的适应性。
关键词: 三维应力     内置环肋圆柱壳     结构特征    
Adaptability analysis of three dimensional stress calculation method for inside ring-stiffened cylindrical shell
CHEN Shagu1,2, GAO Yuan1,2, ZHANG Pingping1,2, ZHOU Ping1,2     
1. China Ship Scientific Research Center, Wuxi 214082, China;
2. Taihu Laboratory of Deepsea Technological Science, Wuxi 214082, China
Abstract: The three dimensional stress analysis method accounting the influence of in-plane normal loads and stress can calculate structural stress at any position of ring-stiffened cylindrical shell. This paper derived typical stress calculation formula for inside ring-stiffened cylindrical shell based on three dimensional stress analysis method, and verified correctness of the calculation formula through examples in literature. A series of calculation are conducted on inside ring-stiffened cylindrical shell using the formulas. Comparing with finite element simulation, the adaptability of typical stress calculation formula to different structural characteristic parameters such as shell thickness, ring-stiffened height, and ring-stiffened width was analyzed. The results show that the formulas in this paper not only can adapt to the stress calculation of thin shell and medium-thick shell, but also have good adaptability to different ring-stiffened shapes.
Key words: three dimensional stress     inside ring-stiffened cylindrical shell     structural characteristic    
0 引 言

深海耐压结构可以为人员、设备或系统提供常压环境,是深海装备不可或缺的重要组成部分。圆柱壳是一种最为典型的深海耐压结构形式,其受力特性、空间利用率和加工工艺等综合性能良好,在潜艇、潜水器、深海无人系统等深海装备中被广泛采用[12]

受静水外压的环肋圆柱壳应力计算一般通过薄壳理论来分析求解,包括2种简化力学模型:一种是将环肋视为圆柱壳梁的弹性支座,假定环肋的形心与圆柱壳壳板形心重合,肋骨的作用以其横剖面面积的方式折算到壳板中曲面上[3];另一种是将环肋圆柱壳分解为肋骨腹板、肋骨翼板和圆柱壳3个部分,分别进行受力分析并通过连接处的变形相等来进行联立求解[4]。这2种力学模型均只考虑了面内二维应力状态,忽略了静水压力载荷以及面内法向应力的影响。在传统薄壳结构范畴,这种处理可以简化工程计算且满足误差要求,但随着载荷的增加和耐压结构壳体厚度设计的增大,这种基于薄壳理论的计算结果会引起较大误差[5]。陈沙古等[6]通过肋骨作用力函数假设,提出一种考虑面内法向载荷和沿厚度方向应力影响的环肋圆柱壳三维应力分析方法,可以用于计算圆柱壳和肋骨任意位置的结构应力。

深海耐压结构尤其是单层耐压结构一般采用内置肋骨进行加强[7],本文根据三维应力分析方法导出了内置环肋圆柱壳的典型应力计算式,并通过文献中的圆柱壳算例进行检验计算。在此基础上,进一步采用三维应力计算方法对不同结构特征参数下的环肋骨圆柱壳结构进行计算,通过与采用实体单元的有限元仿真计算结果进行比较,研究分析导出的内置环肋圆柱壳三维应力计算方法的适应性。

1 三维应力计算方法 1.1 内置环肋圆柱壳力学模型

假定肋骨对圆柱壳的支撑等效为内作用力函数f(z);肋骨腹板与圆柱壳之间的内作用力f1;肋骨腹板与面板之间的内作用力f2。由于内置肋骨不直接承受水压,因此环肋骨圆柱壳各部分的径向作用力可以简化为如图1所示的力学模型。其中,p为圆柱壳梁环外表面承受面载荷;f(z)为内表面承受等效面载荷;f1为肋骨腹板圆环外圆承受作用载荷;f2为内圆承受作用载荷;f3为肋骨面板圆环外圆承受作用载荷。

图 1 内置环肋圆柱壳结构力学模型 Fig. 1 Mechanical model of inside ring-stiffened cylindrical shell
1.2 三维应力计算方法

首先定义内置环肋圆柱壳结构的一般参数:Ri为圆柱壳内半径,Ro为外半径,t为厚度;H为内肋骨腹板高度,δ为厚度,W为面板宽度,b为厚度,l为相邻肋骨之间的间距;Rfo=RiH为肋骨面板外圆半径,Rfi=RiHb为内圆半径;E为材料弹性模量,μ材料泊松比。

参考文献[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)=m21m2p+E(rReff)f1(1μ2)[Q3coshα1zcosα2z+Q4sinhα1zsinα2z]ur(r,z)=1E[(1+μ)A(z)r+2(1μ)C(z)r] (1)

式中:σr(r,z)σθ(r,z)σz(r,z)ur(r,z)分别为圆柱壳梁环的径向应力、周向应力、纵向应力和径向位移,均为双变量函数;变量r为沿极坐标下圆柱壳的径向方向,z为沿圆柱壳的母线方向;A(z)C(z)为与z相关的参数,m=Ro/Ri为圆柱壳内外半径之比。Reff为圆柱壳等效弯曲面对应的曲率半径,满足关系式要求ur(Reff,z)=1tRoRiur(r,z)dr,可导出Reff=t/lnm,圆柱壳的中面半径Rm相比误差为Δ=1Rmtlnm,对于薄壳和一般中厚壳来说误差很小,为简化可取ReffRm

肋骨腹板的应力和变形可以表达为

{σrF(r)=Ri2r211n2(f1δ+f2δ)+(n21n2f1δ11n2f2δ)σθF(r)=Ri2r211n2(f1δ+f2δ)+(n21n2f1δ11n2f2δ)urF(r)=1E[Ri2r(1+μ)1n2(f1δf2δ)+r[(1μ)n21n2f1δ1μ1n2f2δ]] (2)

式中:σrF(r)σθF(r)urF(r)分别为肋骨腹板圆环的径向应力、周向应力和径向位移,均为单变量函数;n=Ri/Rfo为肋骨腹板圆环的内外圆半径之比。

肋骨面板的应力和变形可以表达为:

{σrM(r)=Rfo2r211k2f2W+k21k2f2WσθM(r)=Rfo2r211k2f2W+k21k2f2WμrM(r)=1E[Rfo2r1+μ1k2f2W+r(1μ)k21k2f2W] (3)

式中:σrM(r)σθM(r)urM(r)分别为肋骨面板圆环的径向应力、周向应力和径向位移,均为单变量函数;k=Rfo/Rfi为肋骨面板圆环的内外圆半径之比。

经推导上述方程中的作用力参数的一般解析式为:

{A(z)=Reff2Q5[Q6p+Etf1(Q7coshα1zcosα2z+Q8sinhα1zsinα2z)]C(z)=12[Q11pEtf1(Q23coshα1zcosα2z+Q24sinhα1zsinα2z)]f(z)=Q22p+Etf1(Q25coshα1zcosα2z+Q26sinhα1zsinα2z)f1=Q20Q12Q13Q17Q21σpf2=WQ14σQ16+WQ15f1 (4)

式中,相关的符号参数如下:

{α=43(1μ2)Refft;γ=3(1μ2)2pReff2Et2;α1=α1γ;α2=α1+γ;u1=α1l2;u2=α2l2Q1=Reff2Etl1γ2u2sinhu1cosu2+u1coshu1sinu2cosh2u1cos2u2Q2=Reff2Etl1γ2u2coshu1sinu2u1sinhu1cosu2cosh2u1cos2u2Q3=Q1α12+2Q2α1α2Q1α22;Q4=Q2α122Q1α1α2Q2α22Q5=(1μ)Reff2Ro2+(1+μ)Q6=μm21m2(10.5μ)RefftQ7=μQ32(1μ2)Q1tReffQ8=μQ42(1μ2)Q2tReffQ9=μ1μm21m2Q10=μ1μ12(1μ2)Q11=Q6Q5Reff2Ro2Q91Q12=(1+μ)+(1μ)n21n2Q13=21n2Q14=2n21n2Q15=(1μ)+(1+μ)n21n2Q16=(1+μ)+(1μ)k21k2Q17=WQ14δQ16+WQ15Q18=Etδ(Q7coshu1cosu2+Q8sinhu1zsinu2)Q19=Etδ(Q23coshu1cosα2+Q24sinhu1sinu2)Q20=(1+μ)Reff2Ri2Q6Q5+(1μ)Q11Q21=(1+μ)Reff2Ri2Q18Q5(1μ)Q19Q22=m21m2Reff2Ri2Q6Q5+1Q23=Q3Q10Q7Q5Reff2Ro2Q24=Q4Q10Q8Q5Reff2Ro2Q25=2Q3Q10+m21m2Q7Q5Reff2Ri2Q26=2Q4Q10+m21m2Q8Q5Reff2Ri2A(0)=Reff2Q5(Q6p+Etf1Q7)A(l/2)=R2effQ5(Q6p+Q18f1δ)C(0)=12(Q11pEtf1Q23)C(l/2)=12(Q11pQ19f1δ) (5)

根据式(1)~式(4),代入坐标参数rz可以计算内置环肋圆柱壳任意部位的结构应力。显然,当W=δ时,则退化为矩形肋骨的特例情形。

2 计算方法适应性分析 2.1 典型应力计算式

环肋圆柱壳为轴对称结构,静水外压作用下其典型部位的应力包括相邻肋骨跨中壳板内、外表面的纵向应力和周向应力,肋骨根部跨端壳板内、外表面的纵向应力和周向应力,肋骨腹板和面板的周向应力。根据三维应力计算方法可导出内置环肋骨圆柱壳的典型应力计算式,如表1所示。

表 1 内置环肋圆柱壳典型应力计算式 Tab.1 Stress calculation formula of inside ring-stiffened cylindrical shell

为检验上述应力计算式,对文献[5]中的内置环肋圆柱壳进行算例验证,圆柱壳内半径为1200 mm;壳板厚度为10 mm;肋骨间距为189 mm;计算压力为6.4 MPa,材料弹性模量为1.96×105 MPa、泊松比为0.3。内置肋骨比较2种结构形式:1)T形肋骨,面板宽度为35 mm、厚度为7 mm,腹板高度为70 mm、厚度为7 mm;2)矩形肋骨,高度为49 mm、厚度为15 mm。表2表明,本文方法计算值与文献的计算值比较一致,证明本文基于三维应力分析方法导出的内置环肋圆柱壳典型应力计算式正确可靠。

表 2 内置环肋圆柱壳算例计算结果 Tab.2 Calculation results of inside ring-stiffened cylindrical shell
2.2 计算式适应性分析

在静水外压作用下,跨中壳板的外表面周向应力、中面周向应力,跨端壳板的内表面纵向应力和内表面周向应力,肋骨最大周向应力和最小周向应力是内置环肋圆柱壳结构的典型应力,具体位置如图2所示。

图 2 内置环肋圆柱壳典型位置应力示意图 Fig. 2 Typical stress of inside ring-stiffened cylindrical shell

通过与有限元仿真计算结果比较研究上述典型应力计算式的适应性,内置环肋圆柱壳为轴对称结构,选用二维实体结构单元PLANE42建立平面轴对称模型进行仿真计算,如图3所示,在模型外部施加线分布载荷p(模拟静水外压),模型一端施加线分布载荷pm2/(m21)(模拟轴向水压),模型一端施加轴向位移约束(模拟约束轴向刚体位移)。

图 3 内置环肋圆柱壳平面实体单元模型仿真计算 Fig. 3 Plane solid element model of inside ring-stiffened cylindrical shell

不同圆柱壳厚度特征(t/Ri)下的圆柱壳计算参数为:Ri=1000 mm,l=180 mm,肋骨高度为90 mm、厚度为20 mm,计算压力p随厚度t变化,有限元网格划分尺寸ESIZE=min(δ/2t/4),材料弹性模量为1.96×105 MPa、泊松比为0.3。改变圆柱壳厚度t进行系列计算,t/Ri取值范围为0.01~0.1,结构参数特征覆盖典型的薄壳和中厚壳结构范畴。将基于三维应力计算式的解析计算值与有限元仿真结果进行无量纲化比值,结果表明两者计算结果具有较好的一致性,其中跨中壳板应力和肋骨应力的相对误差均很小,跨端壳板的纵向应力的相对误差稍大(跨端处局部结构突变,应力集中的仿真计算结果受单元网格尺寸影响较大),如图4所示。

图 4 不同圆柱壳厚度特征下的典型应力计算结果比值 Fig. 4 Typical stress results of different thickness of cylindrical shell

不同肋骨高度特征(H/Ri)下的圆柱壳计算参数为:Ri=1000 mm,t=20 mm,l=180 mm,肋骨厚度为20 mm,计算压力p=20 MPa,有限元网格划分尺寸ESIZE=min(δ/2, t/4),材料弹性模量为1.96×105 MPa、泊松比为0.3。改变肋骨高度H进行系列计算,H/Ri取值范围为0.03~0.3。将两者计算结果进行无量纲化比值,结果表明两者计算结果具有较好的一致性,除跨端壳板的纵向应力相对误差稍大之外,跨中壳板应力和肋骨应力的相对误差均很小,如图5所示。

图 5 不同肋骨高度特征下的典型应力计算结果比值 Fig. 5 Typical stress results of different height of ring-stiffened

不同肋骨宽度特征(W/l)下的圆柱壳计算参数为:Ri=1000 mm,t=20 mm,l=180 mm,肋骨面板厚度为20 mm,腹板高度为200 mm、厚度为20 mm,计算p=20 MPa,有限元网格划分尺寸ESIZE=min(δ/2, t/4),材料弹性模量为1.96×105 MPa、泊松比为0.3。改变肋骨宽度W进行系列计算,W取值范围为20~120 mm。将两者计算结果进行无量纲化比值,结果表明两者计算结果具有较好的一致性,除跨端壳板的纵向应力相对误差稍大之外,跨中壳板应力和肋骨应力的相对误差均很小,如图6所示。

图 6 不同肋骨宽度特征下的典型应力计算结果比值 Fig. 6 Typical stress results of different width of ring-stiffened
3 结 语

本文以内置环肋圆柱壳结构为对象,开展基于三维应力分析方法的应力计算方法适应性研究,主要工作和结论如下:

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
内置环肋圆柱壳三维应力计算方法适应性分析
陈沙古, 高原, 张平平, 周平  &nbs...