圓柱殼變形計算方法研究

周海波，彭云飛，吳鴻敏，計　方(. 武漢第二船舶設計研究所，湖北　武漢 430064；. 中國艦船研究院，北京 009)

?

圓柱殼變形計算方法研究

周海波1，彭云飛1，吳鴻敏1，計方2
(1. 武漢第二船舶設計研究所，湖北武漢 430064；2. 中國艦船研究院，北京 100192)

摘要:圓柱殼的變形特性計算一直受到廣泛的重視，解決這一問題的經典板殼理論解法主要有解析法、數值法和半解析法。由于環肋圓柱殼結構的復雜性，采用純解析的方法求解非常困難。隨著計算機技術的發展，有限元法在船體結構分析中已占主導地位，并已成為數值解法的典型代表。本文通過對底部簡支和懸臂薄壁圓柱殼的受力狀態進行分析，建立有限元模型進行仿真計算并與理論計算結果加以比對，分析不同計算方法對圓柱殼變形結果的影響。

關鍵詞：圓柱殼；變形；解析法；數值法

0　引　言

水下結構物的主要載荷是承受深水靜壓力，這一受力特性決定了其耐壓結構采用橫剖面為圓形的薄殼結構，同時為了有效地提高耐壓結構的穩定性，還必須在殼板上設置一系列橫向加強筋（肋骨）。因此，從結構力學的角度來看，水下結構物的耐壓船體主要是以一系列橫向肋骨加強的環肋圓柱殼。雖然耐壓船體鋼板厚度較大，但其直徑較大，屬于薄壁結構。解決這一問題的經典板殼理論解法主要有解析法、數值法和半解析法。由于環肋圓柱殼結構的復雜性，采用純解析的方法求解非常困難。隨著計算機技術的發展，有限元法在船體結構分析中已占主導地位，并已成為數值解法的典型代表。由于有限元法把連續的結構離散成若干個幾何形狀規則的理想單元，所以能夠解決復雜的高次靜不定問題。但有限元法計算結果的精確與否，直接依賴于軟件使用的合理性和網格劃分的準確性，所以又必須以經典的板殼理論作支撐。

殼體結構（即使是薄殼結構）的穩定性問題就其本質來說，應該歸屬于非線性理論范疇，但是一些重要的基本結果仍然是從線性理論中獲得。薄殼結構的分支點問題可以用線性理論判別其特征。目前工程計算中大多仍是以線性理論解為基礎，依靠大量試驗的結果得出某種規律，給出各類近似的修正公式。當前對于線性理論解的研究，采用的主要方法有能量法、傳遞矩陣法等。

隨著現代力學、計算數學和計算機技術的發展，數值計算無論在理論，還是在計算技術方面都取得了巨大的進步。常用的數值計算方法包括有限元法、有限差分法和邊界元法，而有限元法基于全離散的原理，可分析復雜結構屈曲和結構的復雜非線性屈曲問題，因此成為殼體結構屈曲分析（尤其是非線性屈曲分析）的首要手段。

本文依據船舶結構力學基礎理論，通過對理想底部簡支和懸臂薄壁圓柱殼的受力狀態進行分析，同時建立有限元模型進行仿真計算并與理論計算結果加以比對，分析不同計算方法對于圓柱殼力學性能結果的影響。

1　圓柱殼的無矩理論

工程中，最大厚度遠小于中面曲率半徑和另外2個方向尺寸的殼體稱為薄殼，薄殼理論是 19 世紀末在基爾霍夫?樂甫假設的基礎上建立起來的。在薄殼理論中，如果殼體的幾何形狀和表面載荷都是連續可微函數，則殼體處于無彎矩的應力狀態，此稱之為柱殼的無矩理論[1]。

如圖1 所示，q1，q2，q3分別為柱殼所受載荷分別在縱向、環向及法向的分量；FT1，FT2及FT12= FT21分別為縱向拉壓力、環向拉壓力及平錯力，則柱殼的無矩理論平衡方程由下式給出[2]：

彈性方程由下式給出：

式中：u，v，w分別為柱殼中面內各點的縱向、環向及法向位移。

現假設一各向同性材料的懸臂圓柱殼，其彈性模量為E，泊松比μ，密度為ρ，長為l，R-r = δ，一端固支，一端自由，如圖2 所示。

圖1　柱殼的無矩理論Fig. 1　Non-distance theoretical of cylindrical shell

圖2　懸臂圓柱殼Fig. 2　Cantilevered cylindrical shell

取圓柱殼中點為坐標原點，因為左右兩邊對稱，所以只對 α 正方向進行計算，每單位面積的載荷為q0= ρδg，q1= 0，q2= q0sinφ，q3= q0cosφ，由式（1）得：

因此，對于懸臂圓柱殼在自身重力作用下，柱殼任一點縱向、環向及法向變形量u，v，w由下式給出[3]：

式中：α 為圓柱殼上任意一點距固定端的距離，自由端 α = l[4]；

2　圓柱殼的有矩理論

如圖3 所示，圓柱殼臥置在平臺上，A 為圓柱殼上任意位置的截面，M0，M1為作用于圓柱殼任意截面A 的彎矩。取縱向右半圓的圓柱殼進行分析，薄壁圓柱殼的截面厚度遠小于曲率半徑，作用于筒體任意截面 A 的垂直力對彎矩M的影響甚微[5]，故省略計算。

A0截面：X0= R，Z0= 0；A2截面：X2= 0，Z2=R；A 截面：X = Rcosφ，Z = Rsinφ。M2為第Ⅰ象限 A2至 A 的圓柱殼自重作用于 A 截面的彎矩，P/2 是平臺對右半圓圓柱殼的反作用力傳遞到 A0截面的值。設第Ⅰ象限的圓柱殼于 A0截面處固定之，則在 M2作用下，A2分別向下和向左位移[6]。

圖3　底部簡支圓柱殼Fig. 3　Bottom simply supported cylindrical shell

式中：F為圓柱殼的截面積；γ為材料比重。

圓柱殼下半圓的第Ⅳ象限，設第Ⅳ象限的圓柱殼于 A0截面處固定之，則在 M1作用下，A1分別向上和向左位移。P/2 是第Ⅰ象限圓柱殼的自重作用于 A0截面的值，M1為平臺對右半圓圓柱殼向上的反作用力 P和第Ⅳ象限圓柱殼向下的自重，共同作用于任意截面A 的彎矩。

則分別作用于第Ⅰ和第Ⅳ象限的彎矩相等M1= M2= M。即作用于圓柱殼任意截面 A 的彎矩相等[7]。

以第Ⅳ象限為例計算彎矩M值：令ε0= 0，。其中，ε0為線應變，ω為角應變。由于理想圓柱殼，其截面 A0和 A1之間的夾角，自重變形后仍為直角，故 A0和 A1之間的角度改變Δφ1= 0。

垂直方向的自重變形量 U 是指垂直位置的圓柱殼中徑在垂直方向的減少量。第Ⅳ象限 A1點在垂直方向往上的位移 U1：

正號表示圓柱殼的最低點 A1垂直向上位移，垂直位置的圓柱殼中徑在垂直方向的減少量 U 由下式給出[8]：

其中：I 為圓柱殼的截面慣性矩；E 為材料彈性模量。

3　有限元仿真

隨著計算機技術的發展，有限元法在船體結構分析中已占主導地位[9]，目前主要有MSC/Ansys/ Abaqus等大型通用有限元軟件，本文通過Ansys建立典型圓柱殼有限元模型，計算底部簡支及懸臂圓柱殼變形問題。

1）理想圓柱殼底部簡支狀態的自重變形

選取圓柱殼的直徑分別為 4 000 mm，5 000 mm，6 000 mm，圓柱殼壁厚為 20 mm，30 mm，40 mm 共 9種計算工況。材料屬性為 ρ = 7.85 × 10-9t/mm3，E = 2.1 × 105MPa，μ = 0.3。模型邊界條件為底部簡支，對整個模型施加 g = 9 800 mm/s2的重力載荷。9種計算工況下圓柱殼變形仿真計算結果如表 1 和圖4 所示，有限元仿真解與有矩理論值對比如表 2 所示。

仿真計算結果表明底部簡支的理想圓柱殼在自身重力作用下垂直變形量隨內徑的增大而增大，隨厚度的增大而減小。同時，仿真解與理論解的差值控制在10% 左右，且圓柱殼直徑越小，厚度越大，即圓柱殼變形值越小，計算誤差越小。

表1　底部簡支圓柱殼仿真計算結果Tab. 1　Simulation calculation results of bottom simply supported cylindrical shell

圖4　底部簡支的理想圓柱殼垂直變形量隨直徑和壁厚變化曲線Fig. 4　Vertical deformation of bottom simply supported cylindrical shell under different diameter and thickness

表2　簡支圓柱殼垂直變形有限元仿真解與有矩理論值對比Tab. 2　The constrast of analytical method and numerical method

2）懸臂圓柱殼的自重變形

為研究懸臂圓柱殼的自重變形，依然選取圓柱殼的直徑分別為 4 000 mm，5 000 mm，6 000 mm，圓柱殼壁厚為 20 mm，30 mm，40 mm 共 9 種計算工況，自由端變形量的仿真結果如表 3 和圖5 所示，有限元仿真解與無矩理論值對比如表 4 所示。

表3　懸臂圓柱殼仿真計算結果Tab. 3　Simulation calculation results of cantilevered cylindrical shell

圖5　懸臂圓柱殼端面變形量隨直徑和壁厚變化曲線Fig. 5　Surface deformation of cantilevered cylindrical shell under different diameter and thickness

表4　懸臂圓柱殼端面變形有限元仿真解與無矩理論值對比Tab. 4　The constrast of analytical method and numerical method

從仿真結果可知，對于自身重力作用下的懸臂圓柱殼，仿真解與理論解誤差較小，且自由端的變形量隨圓柱殼內徑的增大而減小，與圓柱殼厚度無關，同一端面各點的變形量大致相同，即對于懸臂圓柱殼，重力載荷對自由端面圓度影響不大。由式（4）可知，不同圓柱殼厚度δ下，無矩理論給出的變形量結果一致，這是因為無矩理論將此類圓柱殼假設為薄壁圓柱殼，忽略了圓柱殼厚度對自由端變形的影響，仿真結果與理論解具有極高的一致性，表明該假設合理。

4　結　語

本文采用無矩理論及有矩理論對理想圓柱殼的變形進行分析，通過有限元方法對典型圓柱殼變形仿真計算并與理論值進行對比，得出以下結論：

1）無矩理論適用于一端固支，一端自由的圓柱殼自重變形計算，有矩理論僅適用于圓柱殼為理想圓柱殼，底部簡支的情況；

2）對于一端固支，一端自由的圓柱殼，仿真結果表明，自由端的變形量隨圓柱殼直徑的增大而減小，與圓柱殼厚度無關，同一端面各點的變形量大致相同，即重力載荷對自由端面圓度影響不大，這與經典無矩理論規律一致；

3）對于理想圓柱殼底部簡支的情況，仿真結果表明，圓柱殼在自身重力作用下，垂直位置的圓柱殼中徑在垂直方向的減少量隨直徑的增大而增大，隨厚度的增大而減小，變形量基本上與內徑的四次方成正比，與厚度的二次方成反比，這與經典有矩理論規律一致；

4）對于復雜圓柱殼結構，圓柱殼的幾何形狀和所受載荷不再是連續可微函數，理論計算具有較大的局限性，限元方法較為合理準確。

參考文獻：

[1]LIU J H, FRANCIS A. Theoretical analysis of local indentation on pressured pipes[J]. International Journal of Pressure Vessels and Piping, 2004, 81(12): 931–939.

[2]JAROSLAV MACKERLE. Fipressured pipes[J]. International Journal of Pressure Vessels and Piping, 2004, 81: 931–939

[3]張鐘華. 簡支圓管在重力作用下的應力狀態[J]. 力學與實踐, 1982, 4(1): 32–35.

[4]王敏中. 簡支圓管在自重作用下的彎曲[J]. 力學與實踐, 1983, 5(6): 52–54.

[5]PLANCQ D, BERTON M N, PIERRE G, et al. Complete analytic elastic study of a 90° branch connection of a cylinder on a sphere subjected to a bending load[J]. International Journal of Pressure Vessels and Piping, 1997, 71(1): 61–69.

[6]ZINGONI A, PAVLOVI? M N. Discontinuity phenomena around the supports of stepwise-thickened spherical steel tanks. Part 1: theoretical considerations and parametric results[J]. International Journal of Pressure Vessels and Piping, 1993, 53(3): 405–435.

[7]呂延茂. 薄壁筒體臥置狀態圓度的探索[J]. 化學工業與工程技術, 2004, 25(1): 49–54.

[8]呂延茂. 薄壁筒體臥置狀態圓度的測量和計算[J]. 壓力容器, 2004, 21(8): 16–21.

[9]MACKERLE J. Finite element analysis of fastening and joining: A bibliography (1990–2002)[J]. The International Journal of Pressure Vessels and Piping, 2003, 80(4): 253–271.

Research on calculation methods of cylindrical shell Deformation

ZHOU Hai-bo1, PENG Yun-fei1, WU Hong-min1, JI Fang2
(1. Wuhan Second Ship Design and Research Institute, Wuhan 430064, China; 2. China Ship Research and Development Academy, Beijing 100192, China)

Abstract:Calculation on the deformation of cylindrical shell is a key problem, analytical method, numerical method and semi-analytic method become effective method to solve these problems. Because of the complexity of ring-stiffened cylindrical shell,it is difficult to use analytical method, FEM becomes the representation of numerical method by the developing of computer technology. The stress condition of typical cylindrical shell was anslysed in this paper, and the FEM model was established. According to the contrast of analytical method and numerical method, The influence of different methods to the cylindrical shell deformation was analysed.

Key words:cylindrical shelldeformation；analytical method；numerical method

作者簡介:周海波(1987–)，男，工程師，研究方向為船舶總體設計。

基金項目:國家自然科學基金資助項目（51409139）

收稿日期:2015–09–06； 修回日期: 2015–11–05

文章編號：1672–7619(2016)03–0001–04

doi：10.3404/j.issn.1672–7619.2016.03.001

中圖分類號：U661.44

文獻標識碼：A