全文快速搜索:   高级搜索

  中国石油大学学报(自然科学版)  2019, Vol. 43 Issue (6): 112-122  DOI:10.3969/j.issn.1673-5005.2019.06.013
0

引用本文 [复制中英文]

曹宇光, 魏云港, 田宏军, 等. 海底管道中泡沫清管器接触应力有限元分析[J]. 中国石油大学学报(自然科学版), 2019, 43(6): 112-122. DOI: 10.3969/j.issn.1673-5005.2019.06.013.
[复制中文]
CAO Yuguang, WEI Yungang, TIAN Hongjun, et al. Finite element analysis of contact stress of foam pig in submarine pipeline[J]. Journal of China University of Petroleum (Edition of Natural Science), 2019, 43(6): 112-122. DOI: 10.3969/j.issn.1673-5005.2019.06.013.
[复制英文]

基金项目

国家自然科学基金项目(11472309);国家重点研发计划项目(2016YFC0802306);中国石油大学(华东)研究生创新工程项目(YCX2018069)

作者简介

曹宇光(1979-),男,教授,博士,博士生导师,研究方向为海洋工程技术、管道完整性评价、断裂失效与断口分析等。E-mail:cao_yuguang@qq.com

文章历史

收稿日期:2019-01-22
海底管道中泡沫清管器接触应力有限元分析
曹宇光, 魏云港, 田宏军, 刘畅     
中国石油大学(华东)山东省油气储运安全重点实验室, 山东青岛 266580
摘要: 为分析清管器在海底管道中的清洁能力, 以泡沫清管器为研究对象, 基于有限元软件ANSYS建立清管器与管道的3维模型, 进行数值模拟, 并通过与理论解析解对比以验证模型准确性。在此基础上, 从清管器接触应力分布出发, 对泡沫密度、管道椭圆度、清管器过盈量、管壁摩擦等影响因素进行研究。结果表明:清管器在无损的均匀直管内接触应力均匀分布、在椭圆管道内则呈“8”字型不均匀分布; 接触应力随清管器过盈量的增大而增大, 但过盈量对清管器表面接触应力分布形式无影响; 管壁摩擦使得清管器表面轴向接触应力分布不均匀, 且接触应力随摩擦系数增大而指数增大; 清管器密度只影响接触应力大小, 对接触应力的分布规律无影响。
关键词: 海底管道    泡沫清管器    接触应力    数值模拟    
Finite element analysis of contact stress of foam pig in submarine pipeline
CAO Yuguang , WEI Yungang , TIAN Hongjun , LIU Chang     
Shandong Province Key Laboratory of Safety of Oil & Gas Storage and Transportation in China University of Petroleum (East China), Qingdao 266580, China
Abstract: In order to analyze the cleaning ability of the pig in the submarine pipeline, the foam pig was taken as the research object, and the 3D models of the pig and pipeline were established based on the finite element software ANSYS. The accuracy of the model was verified by comparing with the theoretical analytical solution. Based on this analysis, starting from the contact stress distribution of the pig, the factors such as the density of the foam, the ovality of the pipe, the magnitude of the interference and the friction of the pipe wall were studied. The results show that the contact stress is uniformly distributed in the nondestructive homogeneous straight pipe, but it is not uniformly distributed in the elliptical pipe shaped by the character "8". The contact stress increases with the increase of the magnitude of the interference, but the magnitude of the interference has no effect on the distribution form of the contact stress on the pig. The friction of the pipe wall makes the distribution of the axial contact stress on the surface of pig nonuniform, and the contact stress increases exponentially with the increase of the friction coefficient. Moreover, the density of the pig only affects the magnitude of the contact stress, but it has no effect on the distribution law of the contact stress.
Keywords: submarine pipeline    foam pig    contact stress    numerical simulation    

随着海上油气资源的大量开采, 海底管道作为海上油气生命线在社会生活中所起到的作用越来越大。海底管道在长久运行后会出现诸如管道结蜡、黏附杂质等影响管输效率甚至堵塞管道的情况, 因此, 需要对海底管道进行定期、有效地清洁。清管器作为清洁管道的主要工具, 在海底管道内的应用不可或缺, 其中泡沫清管器结构简单、形变量大、耐磨性好、通过能力强、密封效果也较好[1-3], 能避免清洁过程中造成的卡堵, 并能较好地清除管道内杂质、很好的完成海底管道清洁任务。清管器与管道之间过盈配合, 通过密封形成压差推动清管器前进, 两者相互作用的关键是其间的接触应力[4], 分析接触应力的分布规律能够为研究清管器在海底管道中的清洁能力起到一定的指导意义。同时考虑到海底管道铺设和所处环境的特殊性, 海底管道会因弯曲而导致管道截面椭圆化[5]。国内外学者针对海底管道椭圆化进行过大量研究[6-12]。椭圆化后的海底管道可能会减小与清管器的接触而影响密封效果, 进而影响清洁效果。此外, 泡沫清管器规格多样, 不同密度、不同过盈量的清管器清洁效果可能不同, 接触应力分布也会有所不同。鉴于此, 笔者以泡沫清管器清洁海底管道为背景, 利用有限元软件ANSYS建立3维数值模型并分析泡沫密度、管道椭圆度、清管器过盈量、管壁摩擦等因素对清管器接触应力分布的影响, 研究泡沫清管器接触应力分布规律。

1 有限元模型验证

为保证ANSYS中建立的有限元模型的准确性, 利用过盈配合的解析解同有限元模拟解进行对比, 在保证有限元模型正确的前提下对清管器接触应力分布进行研究。

1.1 过盈配合的解析解

清管器与管道之间过盈配合(即清管器直径大于需要清洁的管道直径), 过盈量计算公式为

$ \delta = \frac{{{D_1} - {D_2}}}{{{D_2}}}. $

式中, δ为清管器过盈量; D1D2分别为清管器直径和管道内径, m。

清管器因过盈而在清洁管道过程中受管道径向挤压进而发生变形。

泡沫清管器是由泡沫橡胶或泡沫塑料制成的圆柱体[13], 长度约为管道内径的1.3~2倍, 其整体形状如图 1所示。为验证有限元模型的正确性, 将清管器和管道简化为轴对称平面问题, 通过解析解求得清管器过盈配合下的轴向变形, 与有限元模拟解对比进行验证。

图 1 聚氨酯泡沫清管器 Fig.1 Polyurethane foam pig

以过盈量2%的泡沫清管器清洁DN450(名义管道外径457 mm, 壁厚14.3 mm)均匀直管为例, 清管器与管道的过盈配合如图 2所示。

图 2 泡沫清管器与管道过盈配合示意图 Fig.2 Schematic diagram of interference fit between foam pig and pipeline

图 2中, Rr为泡沫清管器半径, b为管道外半径, Rp为管道内半径。

根据超弹性材料过盈配合的轴对称平面应力解答[14], 结合Cao等[15]求解聚氨酯皮碗过盈配合解析解时所做的工作, 通过拉梅公式、平面应力下的本构方程、几何方程, 可求得管道的应力、应变分量以及径向位移分别为

$ \left\{ \begin{array}{l} {\sigma _{rr}} = \frac{{{{\left( {{R_r} - {\delta _1}} \right)}^2}}}{{{b^2} - {{\left( {{R_r} - {\delta _1}} \right)}^2}}}\left( {1 - \frac{{{b^2}}}{{{r^2}}}} \right){p_0},\\ {\sigma _{\theta \theta }} = \frac{{{{\left( {{R_{\rm{r}}} - {\delta _1}} \right)}^2}}}{{{b^2} - {{\left( {{R_{\rm{r}}} - {\delta _1}} \right)}^2}}}\left( {1 + \frac{{{b^2}}}{{{r^2}}}} \right){p_0}. \end{array} \right. $ (1)
$ {\varepsilon _{\theta \theta }} = \frac{1}{E}\frac{{{{\left( {{R_{\rm{r}}} - {\delta _1}} \right)}^2}}}{{E{b^2} - {{\left( {{R_{\rm{r}}} - {\delta _1}} \right)}^2}}}\left[ {1 - \mu + \left( {1 + \mu } \right)\frac{{{b^2}}}{{{r^2}}}} \right]{p_0}, $ (2)
$ u = R{\varepsilon _{\theta \theta }}. $ (3)

泡沫清管器由泡沫橡胶或泡沫塑料制成, 其中弹性聚氨酯泡沫塑料的生产规模大, 用途广泛[16], 选用材料为聚氨酯泡沫的清管器进行研究。聚氨酯泡沫为超弹性材料, 超弹性材料的本构关系为

$ S' = \frac{2}{J}DEV\left[ {\left( {\frac{{\partial U}}{{\partial {{\bar I}_1}}} + {{\bar I}_1}\frac{{\partial U}}{{\partial {{\bar I}_2}}}} \right)\bar B - \frac{{\partial U}}{{\partial {{\bar I}_2}}}\bar B \cdot \bar B} \right]. $ (4)

式中, S′为Cauchy应力的偏量; U为应变能密度。

B的第一和第二不变量${\bar I_1} = {\bar I_2} = 1 + {r^2}/{R^2} + {R^2}/{r^2}$, 记做I

利用超弹性材料的本构关系得

$ {{\sigma '}_{\theta \theta }} - {{\sigma '}_{rr}} = 2\left[ {\left( {\frac{{\partial U}}{{\partial {{\bar I}_1}}} + {{\bar I}_1}\frac{{\partial U}}{{\partial {{\bar I}_2}}}} \right)\left( {\frac{{{r^2}}}{{{R^2}}} - \frac{{{R^2}}}{{{r^2}}}} \right) - \frac{{\partial U}}{{\partial {{\bar I}_2}}}\left( {\frac{{{r^4}}}{{{R^4}}} - \frac{{{R^4}}}{{{r^4}}}} \right)} \right]. $ (5)

式中, R为变形前泡沫清管器内某点的半径, m; r为过盈配合后该点的半径, m。

利用轴对称平面问题的平衡方程

$ \frac{{{\rm{d}}{\sigma _{rr}}}}{{{\rm{d}}r}} + \frac{1}{r}\left( {{\sigma _{rr}} - {\sigma _{\theta \theta }}} \right) = 0. $ (6)

可推出泡沫清管器的应力计算公式:

$ \left\{ {\begin{array}{*{20}{l}} {{\sigma _{rr}} = \int_0^{{R_{\rm{r}}}} {\frac{1}{r}} \left( {{{\sigma '}_{\theta \theta }} - {{\sigma '}_{rr}}} \right){\rm{d}}r,}\\ {{\sigma _{\theta \theta }} = {\sigma _{rr}} + \left( {{{\sigma '}_{\theta \theta }} - {{\sigma '}_{rr}}} \right).} \end{array}} \right. $ (7)

式中, σ′θθσ′rr分别为应力分量σθθσrr的应力偏量。

假设泡沫清管器的长度为1 m, 变形后轴向伸长量为z, 材料取中高密度聚氨酯泡沫, 忽略材料的体积可变性, 体积不变条件可表示为

$ {R^2} \times 1 = {r^2}\left( {1 + z} \right), $ (8)

$ r = R/\sqrt {1 + z} . $ (9)

将式(9)代入式(7)中第一式, 得

$ \begin{array}{l} {\left. {{\sigma _{rr}}} \right|_{R = {R_{\rm{r}}}}} = - {P_0} = \int_0^R {\frac{{\sqrt {1 + z} }}{R}} \left( {{{\sigma '}_{\theta \theta }} - {{\sigma '}_{rr}}} \right){\rm{d}}\frac{R}{{\sqrt {1 + z} }} = \\ \int_0^{{R_{\rm{r}}}} {\frac{1}{R}} \left( {{{\sigma '}_{\theta \theta }} - {{\sigma '}_{rr}}} \right){\rm{d}}R = \int_0^{{R_{\rm{r}}}} 2 \left[ {\left( {\frac{{\partial U}}{{\partial {{\bar I}_1}}} + {{\bar I}_1}\frac{{\partial U}}{{\partial {{\bar I}_2}}}} \right) \times } \right.\\ \left. {\left( {\frac{{{r^2}}}{{{R^2}}} - \frac{{{R^2}}}{{{r^2}}}} \right) - \frac{{\partial U}}{{\partial {{\bar I}_2}}}\left( {\frac{{{r^4}}}{{{R^4}}} - \frac{{{R^4}}}{{{r^4}}}} \right)} \right]\frac{1}{R}{\rm{d}}R. \end{array} $ (10)

由接触面上的位移关系rr=Rr-δ+u|R=Rr-δ, 其中以rr表示Rr对应的过盈后的半径, u|R=Rr-δ表示外层钢半径为Rr-δ处的位移, 可得

$ \frac{{{R_{\rm{r}}}}}{{\sqrt {1 + z} }} - \left( {{R_{\rm{r}}} - {\delta _1}} \right) = \frac{{\left( {{R_{\rm{r}}} - {\delta _1}} \right)}}{E}\left[ {\mu + \frac{{{b^2} + {{\left( {{R_{\rm{r}}} - {\delta _1}} \right)}^2}}}{{{b^2} - {{\left( {{R_{\rm{r}}} - {\delta _1}} \right)}^2}}}} \right]{p_0}. $ (11)

泡沫清管器种类繁多, 其中裸体泡沫清管器由高强度聚氨酯材料直接模具发泡成型, 泡沫的密度可以通过添加或多或少的吹气剂来改变[17], 密度变化为16~160 kg/m3[18], 选用常见的80 kg/m3(中密度), 同时为了研究泡沫密度对清管器接触应力的影响, 将选用160 kg/m3(高密度)进行对比。两种密度的聚氨酯泡沫应力应变曲线如图 3所示。可将试验得到的应力应变曲线代入有限元材料库[19], 并在ANSYS中选用5参数的Mooney-Rivlin模型对应力应变曲线进行拟合, 从而得到材料参数[20], 得到的具体材料参数如表 1所示。

图 3 两种密度泡沫的应力应变曲线 Fig.3 Stress-strain curves of two density foam
表 1 材料参数 Table 1 Material parameters

以密度为80 kg/m3的聚氨酯泡沫为例, 其超弹性本构关系为C10=-16.439, C01=18.804, C20=33.9, C11=-99.902, C02=81.086。本构关系Mooney-Rivlin应变能函数为

$ \begin{array}{l} U = - 16.439\left( {{{\bar I}_1} - 3} \right) + 18.804\left( {{{\bar I}_2} - 3} \right) + 33.9{\left( {{{\bar I}_1} - 3} \right)^2} - \\ 99.902\left( {{{\bar I}_1} - 3} \right)\left( {{{\bar I}_2} - 3} \right) + 81.086{\left( {{{\bar I}_2} - 3} \right)^2}. \end{array} $

代入式(10)得

$ \begin{array}{l} - {p_0} = \int_0^{{R_{\rm{r}}}} {\frac{2}{R}} \left( {\frac{{ - 57.971}}{{1 + z}} + \frac{{30.168}}{{{{\left( {1 + z} \right)}^2}}} + 27.803 - } \right.\\ \left. {2.365z - 30.168{z^2}} \right){\rm{d}}R. \end{array} $ (12)

经计算, 过盈2%时清管器过盈量为4.284 mm, 清管器半径为218.484 mm; 将Rr=0.218 484 m, δ=0.004 284 m, b=0.228 5 m, E=201 GPa, μ=0.3代入式(11)得

$ \frac{{3.383}}{{2.01 \times {{10}^5}}}{p_0} = \frac{{0.218\;484}}{{\sqrt {1 + z} }} - 0.2142. $ (13)

Rr=0.218 484(m)代入式(12)得

$ \begin{array}{l} {p_0} = \frac{{ - 176.23184}}{{{{\left( {1 + z} \right)}^2}}} + \frac{{91.71072}}{{{{\left( {1 + z} \right)}^2}}} + 84.52112 - 7.1896z - \\ 91.710\;72{z^2}. \end{array} $ (14)

将式(14)代入式(13)得

$ \begin{array}{l} \frac{{0.218\;484}}{{\sqrt {1 + z} }} - 0.2142 = \frac{{3.383}}{{2.01 \times {{10}^5}}} \times \\ \left[ {\frac{{ - 176.231\;84}}{{{{\left( {1 + z} \right)}^2}}} + \frac{{91.710\;72}}{{{{\left( {1 + z} \right)}^2}}} + 84.521\;12 - 7.189\;6z - } \right.\\ \left. {91.710\;72{z^2}} \right]. \end{array} $

最终利用MATLAB求得轴向变形z=40.5 mm, p0=-0.582 MPa。

1.2 有限元模型数值解

截取长度3 m的DN450管道, 将泡沫清管器简化为长1 m的圆柱, 通过ANSYS并根据清管器与管道结构的对称性建立1/4实体模型。划分网格:管道沿环向划分12份, 沿轴向划分120份; 清管器沿环向划分12份, 沿轴向划分35份。设置管道内壁和清管器外壁为接触对并设置接触对参数。管道和清管器均设置对称位移约束, 并在管道外壁设置全约束, 同时限制清管器某一端的位移。有限元建模与接触对如图 45所示。

图 4 建模与网格划分 Fig.4 Modeling and meshing
图 5 接触对设置 Fig.5 Contact pair setting

采用Newton-Raphson迭代算法对模型进行求解, 得到清管器在过盈配合下的轴向变形与径向变形如图 6所示。其中有限元解的径向压缩为4.231 64 mm, 轴向伸长39.884 6 mm, 与解析解的误差分别为径向误差1.2%, 轴向误差1.5%。误差较小, 验证了有限元模型的正确性。

图 6 清管器径向和轴向位移云图 Fig.6 Radial and axial displacement nephogram of pig
2 管道椭圆度对接触应力分布影响

工程中常以椭圆度定义椭圆管道外形, 根据DNV-OS-F101规范的定义, 管道椭圆度为${f_0} = \frac{{{D_{\max }} - {D_{\min }}}}{D}$, 其中f0为管道椭圆度, DmaxDmin分别为最大、最小管道外径, m; D为名义管道外径[21]。椭圆管道横截面如图 7所示。

图 7 椭圆管道横截面 Fig.7 Cross section of elliptical pipeline

规范规定, 安全的海底管道椭圆度不得超过3%, 现选取椭圆度f0分别为0、1%、2%、3%的管道进行数值模拟分析, 研究泡沫清管器在较小椭圆度管道内的接触应力变化。选用不同椭圆度的DN450管道尺寸如表 2所示。由表 2可知随着管道椭圆度增大, 最小管道外径逐渐减小, 而最大管道外径逐渐增大。

表 2 管道尺寸 Table 2 Pipe size

利用ANSYS前处理器构建管道与清管器模型, 管道长度为3 m。聚氨酯泡沫清管器一般为子弹状, 长径比为1.3~2, 选用清管器总长度为850 mm, 其中有效密封长度为725 mm; 考虑模型特征并为了提高效率, 建立三维1/4实体模型, 选用Solid185单元进行扫掠网格划分(图 8), 模型在设置合理的对称约束后可进行扩展, 扩展后的清管器有限元计算模型如图 9所示。

图 8 管道与清管器1/4模型 Fig.8 1/4 model of pipe and pig
图 9 扩展后的泡沫清管器模型 Fig.9 Expanded model of foam pig

模型中管道相对于清管器可视为刚性体, 其弹性模量为201 GPa, 泊松比为0.3, 根据实际情况, 对管道外壁施加固定约束。所研究问题属于非线性大变形接触问题, 分别用conta174和targe170单元定义接触单元与目标单元, 并选取清管器外表面为接触面、管道内壁为目标面。在清管器尾部施加轴向位移载荷, 并定义左侧为前方, 使其进入管道并沿管道轴线由右向左运动。

以椭圆度2%的管道为例, 过盈量2%清管器在椭圆管道内运行时, 因管道外形不规则而引发清管器的不规则变形。如图 10所示, 清管器主体在管道最小内径处发生最大径向压缩, 在管道最大内径处压缩量最小, 清管器头部子弹头处几乎不压缩反而因主体的径向压缩而轻微地膨胀。

图 10 清管器径向位移云图 Fig.10 Radial displacement nephogram of pig

清管器不均匀变形会导致表面接触应力分布不均匀。为研究接触应力的这种不均匀分布, 现选取清管器与管道最小内径接触处为0°并沿顺时针依次规定角度(规定后的角度分布见图 11), 提取过盈量2%的清管器主体在不同椭圆管道内的环向接触应力进行分析。提取后的中高密度泡沫清管器主体部位的接触应力环向分布如图 12所示。

图 11 横截面及角度分布 Fig.11 Cross section and angle distribution
图 12 不同椭圆管道内接触应力分布 Fig.12 Contact stress distribution in different elliptical pipes

图 12可知, 同一过盈量、不同密度的泡沫清管器在同一椭圆管道内的环向接触应力分布形式类似, 但在不同椭圆管道内的接触应力分布规律差别较大。当清管器在无椭圆度的均匀管道内运行时, 清管器表面的接触应力均匀分布; 而在椭圆管道内运行时, 清管器接触面上应力分布不均匀, 且管道椭圆度越大, 接触应力分布越不均匀。这种不均匀性体现在接触应力在与最小管道内径接触处取得最大值, 在最大管道内径处取得最小值; 且管道椭圆度越大, 最大最小值间的差值越大。这是因为清管器在清洁过程中因过盈受到管道挤压, 且管道内径越小挤压越明显, 产生的接触应力也就越大。形状规则的泡沫清管器在椭圆管道最小内径处受挤压最严重, 因此在此处产生最大的接触应力。

图 12可知, 泡沫清管器在椭圆管道内的接触应力分布类似8字型曲线。泡沫清管器在管道内环向的接触应力分布表示为

$ {\sigma _\theta } = \sqrt {\sigma _{\max }^2{{\cos }^2}\theta + \sigma _{\min }^2{{\sin }^2}\theta } . $ (15)

式中, σmaxσmin分别为清管器在椭圆管道内的最大、最小接触应力; θ为清管器的环向角度。

通过origin软件利用式(15)对环向接触应力进行拟合(图 13), 可见拟合效果较为理想。

图 13 不同椭圆管道内接触应力分布拟合 Fig.13 Fitting of contact stress distribution in different elliptical pipes

为完善泡沫清管器的接触应力分布公式, 提取接触面上最大、最小接触应力并绘制成曲线(图 14), 研究接触应力随管道椭圆度变化规律。

图 14 接触应力随椭圆度变化 Fig.14 Variation in contact stress with ovality

图 14可知, 中高密度泡沫清管器接触应力随椭圆度变化规律相近, 且高密度泡沫的接触应力约为中密度泡沫的4倍。清管器表面最大接触应力随椭圆度增大线性增大, 最小接触应力随椭圆度增大线性减小。这与管道最大、最小内径随椭圆度的变化相关:随管道椭圆度增大, 最小管道内径减小, 对清管器的挤压作用更加明显, 清管器与最小内径接触处产生的最大接触应力也随之增大; 最大管道内径随椭圆度增大而增大, 此处与清管器的接触效果越来越差, 产生的接触应力也越来越小。因为接触应力随椭圆度变化的线性相关性, 利用origin对其拟合可求得σmaxσmin关于椭圆度f0的关系:

$ {\sigma _{\max }} = \gamma + \eta {f_0},{\sigma _{\min }} = \gamma + \lambda {f_0}. $ (16)

式中, γ为泡沫清管器在无椭圆度管道内的接触应力(可由第一节中的解析模型进行求解); ηλ分别为最大、最小接触应力随椭圆度变化曲线的斜率。从图 14中还可发现, 高密度泡沫清管器的曲线斜率约为中密度泡沫清管器的4倍。

3 清管器过盈量对接触应力分布影响

同一清管器在不同椭圆管道内的应力分布不同, 不同过盈量清管器在同一管道内的应力分布也会有所不同。泡沫清管器过盈量一般为2%~4%[22], 选取过盈量分别为2%、2.5%、3%和3.5%的清管器, 对应半径分别为218.484、219.555、220.626和221.697 mm, 分析其在管道内接触应力的不同。

为了控制变量, 不考虑管道对清管器接触应力的影响, 管道椭圆度均设置为1%, 分析过盈量对清管器接触应力的影响。

提取中高密度、不同过盈量清管器在椭圆度1%管道内运行时表面的环向接触应力分布(图 15)。

图 15 不同过盈量清管器接触应力分布 Fig.15 Contact stress distribution with different magnitude of interference

图 15可知, 不同过盈量的泡沫清管器在同一管道内运行时的接触应力分布形式类似, 在椭圆管道内均呈8字型, 且过盈量越小8字型越明显。这说明清管器过盈量越小, 其在椭圆管道内的接触应力分布越不均匀(清管器与椭圆管道最大内径、最小内径接触处的应力差值明显增大)。

利用式(15)同样能对不同过盈量下清管器环向接触应力分布进行拟合, 拟合效果如图 16所示。

图 16 不同过盈量清管器接触应力分布拟合 Fig.16 Fitting of contact stress distribution with different magnitude of interference

图 16可知, 拟合效果较为理想, 可间接证明接触应力分布公式的可靠性。为更清晰地分析接触应力随过盈量变化趋势, 提取不同过盈量清管器表面的最大、最小接触应力并绘制成曲线如图 17所示。

图 17 接触应力随过盈量变化曲线 Fig.17 Variation in contact stress with magnitude of interference

当泡沫清管器在椭圆度1%管道内运行时, 清管器与管道内壁接触产生应力, 且清管器过盈量越大, 接触面上的接触应力越大。这是因为清管器过盈量越大, 清管器清洁管道时受挤压越严重, 因此产生更大的应力。根据接触应力随过盈量变化的线性相关性, 利用Origin对其拟合可求得σmax、σmin关于过盈量δ的关系:

$ {\sigma _{\max }} = \gamma ' + \zeta \delta ,{\sigma _{\min }} = \gamma ' + \nu \delta . $ (17)

式中, γ′为过盈量2%的清管器在管道内的接触应力; ζν分别为最大、最小接触应力随过盈量变化曲线的斜率。

同时还可发现高密度泡沫清管器的曲线斜率约为中密度泡沫清管器的5倍。此外, 结合图 1617可知, 泡沫密度只改变清管器接触应力的大小, 对清管器表面接触应力分布规律无影响, 对接触应力随过盈量变化的形式同样无影响。

4 摩擦系数对接触应力分布影响

清管器在清洁管道时, 因清管目的和管道内杂质的不同, 其与管壁间会产生不同的摩擦, 这可能会对清管器的接触应力分布产生影响。

不考虑管道椭圆度和清管器过盈量的影响, 选取管道椭圆度为0, 清管器过盈量为2%的基本工况, 依次设置摩擦系数为0、0.05、0.1、0.15、0.2进行有限元模拟。模拟后首先提取高密度泡沫清管器在摩擦系数分别为0和0.1时的接触应力云图(图 18)进行对比分析。

图 18 泡沫清管器接触应力云图 Fig.18 Contact stress nephograms of foam pig

当泡沫清管器在无摩擦的管道内运行时, 清管器头部与管道接触处存在接触应力变化区, 清管器主体部位接触应力均匀分布; 在受到摩擦的影响时, 清管器主体表面轴向接触应力不再均匀分布, 且越靠近尾部接触应力越大。为更形象研究这一规律, 提取不同摩擦系数时清管器头部到尾部轴向的接触应力绘制成曲线(图 19), 分析接触应力随摩擦系数变化规律。

图 19 不同摩擦系数下接触应力沿轴向变化 Fig.19 Variation in contact stress along axial direction with different friction coefficient

图 19可知, 中高密度下清管器接触应力沿轴向分布形式相同。清管器“子弹头”部位在管道内运行时未接触管壁, 因而无接触应力。管壁无摩擦时, 随着清管器在管道内运行, 清管器与管壁最先接触的部位接触应力取得最大值, 后接触应力逐渐减小并趋于稳定; 当管壁有摩擦时, 接触应力先减小后增大, 且摩擦系数越大, 接触应力增大的越明显。

为进一步分析接触应力随摩擦系数变化规律, 提取清管器整个表面的平均接触应力并绘制曲线(图 20)。由图 20可知, 摩擦系数对清管器接触应力的变化影响较大, 随摩擦系数增大接触应力非线性增大, 且摩擦越大, 接触应力增大的越明显。这是因为摩擦系数越大, 清管器行进越困难, 其在管道内运行时产生的应力便随之增大。

图 20 平均接触应力随摩擦系数变化 Fig.20 Variation in average contact stress with friction coefficient
5 结论

(1) 清管器清洁带椭圆度的管道时, 因椭圆管道对清管器的不规则挤压, 使得清管器表面接触应力呈8字型分布, 最大、最小接触应力σmaxσmin与管道椭圆度、泡沫密度、清管器过盈量等因素相关。

(2) 清管器接触应力随管道椭圆度线性变化, 且最大接触应力随椭圆度增大而增大、最小接触应力随椭圆度增大而减小。这与管道最大、最小外径随管道椭圆度变化呈正相关关系。

(3) 接触应力随清管器过盈量增大呈线性增大关系, 而泡沫密度只影响接触应力值, 对清管器表面接触应力的环向分布形式和接触应力随过盈量变化的趋势均无影响。

(4) 摩擦使接触应力轴向分布不均匀, 且摩擦系数越大应力分布越不均匀; 接触面上的平均接触应力随管壁摩擦增大而增大。

参考文献
[1]
季伟, 杨全蔚, 冯维, 等. 泡沫型清管器在子洲气田的应用效果评价[J]. 石油化工应用, 2016, 35(11): 69-71.
JI Wei, YANG Quanwei, FENG Wei, et al. Evaluation of the application of foam pigging in Zizhou gas field[J]. Petrochemical Industry Application, 2016, 35(11): 69-71. DOI:10.3969/j.issn.1673-5285.2016.11.018
[2]
刘鹏飞, 王雪, 熊波. 浅谈清管器在天然气管线中的运行特点[J]. 广西轻工业, 2011, 27(3): 15-17, 25.
LIU Pengfei, WANG Xue, XIONG Bo. Operation characteristics of pig in natural gas pipeline[J]. Guangxi Journal of Light Industry, 2011, 27(3): 15-17, 25. DOI:10.3969/j.issn.1003-2673.2011.03.007
[3]
刘淼, 曹学文, 杨文, 等. 清管器运行过程中的旁通特性[J]. 油气田地面工程, 2016, 35(10): 32-34.
LIU Miao, CAO Xuewen, YANG Wen, et al. Study of pig's bypass character in running process[J]. Oil-Gas field Surface Engineering, 2016, 35(10): 32-34. DOI:10.3969/j.issn.1006-6896.2016.10.010
[4]
陈浩, 刘小明. 清管器皮碗在不同过盈量下的接触应力分析[J]. 管道技术与设备, 2015, 136(6): 55-58.
CHEN Hao, LIU Xiaoming. Analysis on contact stress of pig's cup under different interference[J]. Pipeline Technique and Equipment, 2015, 136(6): 55-58.
[5]
KEMP D, GOSSARD J, FINNERAN S, et al. Evaluating pipeline ovality acceptability criteria for straight pipe sections[R]. V001T03A078, 2016.
[6]
PARK T D, KYRIAKIDES S. On the collapse of dented cylinders under external pressure[J]. International Journal of Mechanical Sciences, 1996, 38(5): 557-578. DOI:10.1016/0020-7403(95)00065-8
[7]
GUARRACINO F, MALLARDO V. A refined analytical analysis of submerged pipelines in seabed laying[J]. Applied Ocean Research, 1999, 21(6): 281-293. DOI:10.1016/S0141-1187(99)00020-6
[8]
BAEK J H. Effect of ovality and eccentricity on collapse pressure of subsea pipeline[R]. IGU2011-P2-34, 2011.
[9]
YAN S, SHEN X, YE H, et al. On the collapse failure of dented pipes under bending moment and external pressure: OCEANS 2016 MTS/IEEE[C]. Monterey: IEEE, 2016: 1-5.
[10]
王彦辉.深水卷管铺设过程中管道的关键力学问题研究[D].大连: 大连理工大学, 2016.
WANG Yanhui. Research on key mechanical problems of the pipe during the reel-lay process in deep water[D].Dalian: Dalian University of Technology, 2016.
[11]
FAN Z, YU J, SUN Z, et al. Effect of axial length parameters of ovality on the collapse pressure of offshore pipelines[J]. Thin-Walled Structures, 2017, 116: 19-25. DOI:10.1016/j.tws.2017.02.032
[12]
薛陆丰, 余杨, 余建星, 等. 椭圆度-凹坑双缺陷海底管道局部屈曲特性研究[J]. 海洋工程, 2018, 36(2): 92-99.
XUE Lufeng, YU Yang, YU Jianxing, et al. Local buckling study on subsea pipeline with ellipticity-denting defect[J]. The Ocean Engineering, 2018, 36(2): 92-99.
[13]
den HEIJER A. Frictional behaviour of pigs in motion[D]. Holland: Delft University of Technology, 2016. https://repository.tudelft.nl/islandora/object/uuid:2b426d43-e934-425f-8a2b-799aa8b1634e?collection=research
[14]
邹雨, 庄茁, 黄克智. 超弹性材料过盈配合的轴对称平面应力解答[J]. 工程力学, 2004, 21(6): 72-75, 83.
ZOU Yu, ZHUANG Zhuo, HUANG Kezhi. The solutions of axisymmetric plane stress for a hyperelastic material interference problem[J]. Engineering Mechanics, 2004, 21(6): 72-75, 83. DOI:10.3969/j.issn.1000-4750.2004.06.011
[15]
CAO Y G, ZHEN Y, SHI Y J, et al. Stress distribution of the power section cup of pipeline inspection gauges by finite element method[J]. Journal of Pipeline Systems Engineering & Practice, 2018, 9(2): 04017039.
[16]
ZHANG X D, MACOSKO C W, DAVIS H T, et al. Role of silicone surfactant in flexible polyurethane foam[J]. Journal of Colloid & Interface Science, 1999, 215(2): 270-279.
[17]
MACKAY N. Available methodologies concerning the treatment and removal of sand from pipelines and associated case studies[J]. Oil & Gas Facilities, 2013, 2(2): 61-71.
[18]
张磊, 李宜震, 王冲, 等. 清管器应用及分类简析[J]. 化学工程与装备, 2013(7): 173-176.
ZHANG Lei, LI Yizhen, WANG Chong, et al. Application and classification of pig[J]. Chemical Engineering & Equipment, 2013(7): 173-176.
[19]
OHWOKA A O. Polyurethane sealing discs on cleaning pigs-characterization and dynamic behavior[D]. Norway: University of Stavanger, 2015.
[20]
冯敏亮, 刘罡, 汪珺, 等. 基于ANSYS的超弹性模型参数确定[J]. 合肥学院学报(自然科学版), 2012, 22(3): 40-45.
FENG Minliang, LIU Gang, WANG Jun, et al. Determination of superelesticity models parameters based on ANSYS[J]. Journal of Hefei University(Natural Sciences), 2012, 22(3): 40-45. DOI:10.3969/j.issn.1673-162X.2012.03.010
[21]
DNV (Det Norske Veritas). Submarine pipeline systems: DNV-OS-F101[S]. Høvik: DNV, 2012.
[22]
吕惠建, 于达, 郑荣荣. 油气管道清管器的研发进展[J]. 管道技术与设备, 2012(4): 43-46.
LV Huijian, YU Da, ZHENG Rongrong. Research progress of pigs in oil and gas pipeline[J]. Pipeline Technique and Equipment, 2012(4): 43-46. DOI:10.3969/j.issn.1004-9614.2012.04.016