一种含不规则腐蚀缺陷钢制管道临界弹塑性屈曲压力计算方法

本发明涉及一种含不规则腐蚀缺陷钢制管道临界弹塑性屈曲压力计算方法，其属于深海油气开发、生产和运输的。

背景技术：

1、海底管道是海上油气开发项目的重要组成部分。海底管道的外压随着水深的增加而逐渐增大。较大的外部压力可能导致管道屈曲、断裂甚至倒塌，其中外部压力屈曲是海底管道事故的重要原因。腐蚀缺陷的存在导致管道壁厚变薄，大大降低了管道的抗压承载能力。因此，腐蚀缺陷的研究对管道的稳定性评价具有重要意义。

2、针对钢质管道在外压作用下临界屈曲荷载计算公式，目前已有完好管道、双腐蚀管道弹塑性屈曲计算公式。但腐蚀是一个复杂而随机的电化学反应过程。随着管道使用时间的增加，相邻的腐蚀缺陷会重叠形成新的不规则腐蚀缺陷。不规则腐蚀缺陷是在缺陷深度的圆周或轴向分布中具有一个或多个主峰的腐蚀缺陷，其可以简化为小的底层缺陷的大的顶层缺陷组合而成的腐蚀缺陷。目前，国内外还没有含不规则腐蚀缺陷管道的弹塑性屈曲压力计算公式。

技术实现思路

1、本发明为弥补深海管道腐蚀缺陷形式方面研究技术的不足，在考虑了相邻腐蚀缺陷重叠成一个不规则腐蚀缺陷的情况下，提出了一种含不规则腐蚀缺陷管道外压临界弹塑性屈曲压力计算方法，其具有较高的准确性以及实用性。

2、本发明为解决以上技术问题所采用的方案为：

3、基于经典的timoshenko壳体屈曲理论，引入切线模量理论，考虑不规则腐蚀缺陷尺寸对管道屈曲响应的影响，推导出含不规则腐蚀缺陷管道的临界弹塑性屈曲压力计算方程。并考虑管道在发生正对称或反对称两种屈曲模式时的屈曲响应，推导出腐蚀管道临界外压计算方法；并通过与以往实验以及有限元结果的对比，验证计算方法的准确性。

4、一种含不规则腐蚀缺陷钢制管道临界弹塑性屈曲压力计算方法，包括以下步骤：

5、s1、假设一条足够长的不规则缺陷管道，其平均半径为r和壁厚t，并受到外部压力p的影响。当管道和缺陷足够长时，沿管道轴线的每个截面上任何点位移都是相同的，管道的屈曲可以简化为二维平面问题。因此，腐蚀缺陷管道的变形问题可以基于平面应变假设进行分析，位移仅与极角θ有关，管段任意点的径向位移为w。此外，由中心线延伸引起的周向位移v非常小，因此可以忽略。管道的不规则缺陷由底层缺陷(大缺陷)和顶层缺陷(小缺陷)组成。根据对称性，取管道截面圆环的一半进行研究，半管段的环可分为三部分，极角范围为：

6、

7、θ为角度，θ1为半底层缺陷夹角角度，θ2为半底层和半顶层腐蚀缺陷角度之和；

8、底层缺陷、顶层缺陷和完整区域的壁厚和平均半径不同，可以表示为：

9、

10、

11、其中，t为管道厚度，t1为底层缺陷处管道壁厚，t2为顶层缺陷处管道壁厚，t3为管道完好部分壁厚，r为管道中曲面半径，r1为底层缺陷处管道中曲面半径，r2为顶层缺陷处管道中曲面半径，r3为管道完好部分中曲面半径。

12、s2、应力应变曲线采用ramberg-osgood方程：

13、

14、其中，ε为应变，σ为应力，e为弹性模量，εy为应变，σy为屈服应力，β为应变硬化指数。

15、环的每一段的切线模量为：

16、

17、

18、

19、et1，et2，et3分别为底层、顶层以及完好管道的切线模量。

20、s3、中曲面的合力矩m可以通过平衡分析获得：

21、m ＝ m0 + pr(w-w0) (8)

22、其中，w0和m0是由微小的初始位移引起的微小的径向位移和弯矩。然后，通过将方程(2)-(9)与timoshenko(1961)的经典理论相结合，可以得到对称和反对称屈曲模式下具有不规则缺陷的管道临界屈曲压力(pb)的微分方程：

23、1)对称屈曲模态

24、

25、

26、

27、2)反对称屈曲模态

28、

29、

30、

31、s4、两种屈曲模态的边界条件分别为：

32、1)对称屈曲模态

33、

34、

35、2)反对称屈曲模态:

36、

37、通过求解特征方程(10)-(12)，可以获得以下解

38、

39、

40、

41、将方程(16)-(18)与对称屈曲模式的边界条件相结合，得到以下方程：

42、

43、其中：

44、i11＝-sin(k2θ1)；

45、i12＝0；

46、i13＝cos(k1θ1)；

47、i14＝-cos(k2θ1)；

48、

49、

50、i21＝sin[k2θ2]；

51、i22＝-sec(k3π)cos[k3(π-θ2)]；

52、i23＝0；

53、i24＝cos[k2θ2]；

54、

55、

56、i31＝-k2cos(k2θ1)；

57、i32＝0；

58、i33＝-k1sin(k1θ1)；

59、i34＝k2sin(k2θ1)；

60、i35＝0；

61、i36＝0；

62、i41＝cos[k2θ2]；

63、i42＝-k3sec(k3π)sin[k3(π-θ2)]；

64、i43＝0；

65、i44＝-k2sin[k2θ2]；

66、i45＝0；

67、i46＝0；

68、

69、

70、

71、

72、

73、

74、i61＝i62＝0；

75、i63＝1；

76、i64＝0；

77、

78、

79、[i]是包含临界屈曲压力的矩阵。显然，a12、a21、a22、a32、w0和m0都是非零参数。因此，矩阵[i]的行列式的值为零。即

80、|i|＝0(20)对称屈曲模式下不规则腐蚀缺陷管道的临界屈曲压力可以通过求解方程(20)获得。

81、s5、对于反对称屈曲模式，通过求解方程(13)-(15)，我们可以得到

82、

83、w2＝ c21sink2θ+ c22cosk2θ, for θ1 ≤ θ ≤ θ2 (22)

84、w3＝ c31sink3θ +c32cosk3θ, for θ2 ≤ θ ≤ π (23)

85、将方程(21)-(23)与反对称屈曲模式的边界条件相结合，可以得到以下方程：

86、

87、其中[j]中的元素为：

88、j11＝sin(k1θ1)；

89、j12＝-sin(k2θ1)；

90、j13＝-cos(k2θ1)；

91、j14＝0；

92、j21＝0；

93、j22＝sin[k2(θ1+θ2)]；

94、j23＝cos[k2θ2]；

95、j24＝(cos[k3θ2]-sin[k3(π-θ2)]/sink3π；

96、j31＝k1cos(k1θ1)；

97、j32＝-k2cos(k2θ1)；

98、j33＝k2sin(k2θ1)；

99、j34＝0；

100、j41＝0；

101、j42＝k2cos[k2θ2]；

102、j43＝-k2sin[k2θ2]；

103、j44＝-sin[k3θ2]+k3cos[k3(π+θ2)]/sink3π.

104、显然，[j]的行列式也是零：

105、|j| ＝ 0 (25)

106、不规则腐蚀缺陷管道在反对称屈曲模式下的临界屈曲压力可以通过求解方程(25)获得。

107、本发明与现有技术相比具有以下优点：

108、1、在考虑腐蚀缺陷间相互影响的情况下在单腐蚀缺陷管道弹塑性屈曲基础上推导出不规则腐蚀缺陷管道弹塑性屈曲压力计算方法，较双腐蚀缺陷管道情况计算结果更为准确。

109、2、考虑了不规则腐蚀缺陷的复杂形状对临界屈曲压力的影响，对管道的屈曲响应公式进行了进一步改进。