CN110793852B - 一种基于复杂应力状态模型的管道应力内检测方法 - Google Patents
一种基于复杂应力状态模型的管道应力内检测方法 Download PDFInfo
- Publication number
- CN110793852B CN110793852B CN201911085775.2A CN201911085775A CN110793852B CN 110793852 B CN110793852 B CN 110793852B CN 201911085775 A CN201911085775 A CN 201911085775A CN 110793852 B CN110793852 B CN 110793852B
- Authority
- CN
- China
- Prior art keywords
- stress
- pipeline
- magnetization
- complex
- expressed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/72—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
- G01N3/06—Special adaptations of indicating or recording means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/0617—Electrical or magnetic indicating, recording or sensing means
- G01N2203/0635—Electrical or magnetic indicating, recording or sensing means using magnetic properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/067—Parameter measured for estimating the property
- G01N2203/0676—Force, weight, load, energy, speed or acceleration
Landscapes
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
一种基于复杂应力状态模型的管道应力内检测方法,该方法步骤如下:第一步:利用检测设备探头采集弱磁信号;第二步:将第一步采集的信号输入至构建的复杂应力状态模型,得出管道应力变化;第三步:根据第二步的管道应力变化识别出管道应力集中区。本发明采用一种复杂应力状态模型,建立磁力学耦合数学模型,提供一种复杂应力状态下的应力检测定量算法,实现了对管道应力集中程度的判定检测。
Description
技术领域
本发明属于长输油气管道应力内检测技术领域,尤其涉及一种基于复杂应力状态模型的管道应力内检测方法。
背景技术
长输油气管道长期在外部载荷作用下工作,即使管道的最大工作应力低于强度极限,甚至低于屈服极限,也可能会在没有明显塑性变形的情况下突然发生泄露、爆管等恶性事故。常规的无损检测技术如磁粉、漏磁、涡流和超声等,在金属缺陷监测、事故预防等方面发挥了重要的作用,然而,这些方法在管道应力内检测领域的应用具有一定局限性,无法预判材料的早期损伤。
弱磁检测方法作为近年来新型的无损检测方法对材料应力集中十分敏感,并且具有设备轻巧、操作简单、可靠性好且不需要专门的磁化设备等特点,是一种应用前景十分广阔的应力无损检测方法。但是,管道在内部载荷作用下,管壁受力情况复杂,现有的理论模型无法满足在役油气管道应力内检测信号识别的要求,限制了弱磁法在管道内检测技术领域的应用。
发明内容
发明目的:
本发明针对上述问题,提供一种基于复杂应力状态模型的管道应力内检测方法。为实现上述目的。
技术方案:
本发明采用如下技术方案:
材料在应力下会通过磁致伸缩改变内部有效场。设整个***受到外磁场H和应力σ共同作用,应力对***磁化的作用可以等效为一个外加磁场Hσ,由磁致伸缩系数λ来表示。当材料在受到压应力作用时,大多数磁畴会向垂直于应力的方向转动,造成局部磁极密度的增加,这就会导致在外场方向上严重的退磁,以此来减弱材料的磁化;当材料在受到拉应力作用时,大多数磁畴会向平行于应力的方向转动,材料的磁化则不会受到较大影响。在拉压应力状态下材料的磁化是不对称的,因此需要引入应力退磁因子Dσ。应力退磁因子随压应力的增加而线性增加,在拉应力下可以忽略不计,则有:
式中,Ms为饱和磁化强度,Bs为饱和磁感应强度,λs为饱和磁致伸缩系数,可表示为:
式中,υ为泊松比,E为弹性模量,b(σ)是与应力相关的函数,可以用一分段函数来表示:
式中,b为常数,σ0为初始应力。
在外磁场作用下,铁磁性材料的磁化包含两部分,一部分是由畴壁弯曲引起的可逆磁化部分,一部分是由畴壁位移引起的不可逆磁化部分。磁化强度M、可逆磁化强度Mrev、不可逆磁化强度Mirr和无磁滞磁化强度Man的关系如下:
M=Mrev+Mirr (4)
Mrev=c(Man-Mirr) (5)
式中,c为可逆系数。
受外磁场和应力作用的无磁滞磁化强度可表示为有效磁场He的函数:
式中,a为材料规划常数。
磁化强度对应变能W的微分可表示为:
式中,ξ为系数,与单位体积的能量有关。
实际工程中,在役油气管道在运行时受到内压产生来自各个方向、复杂不定的拉伸应力,此时管道处于复杂应力状态。
由材料力学可知,复杂应力状态下材料单位体积的应变能可以表示为:
式中,σ1、σ2、σ3表示单元体的三个主应力。
材料的应变能分为因形状改变而存储的畸变能ef和因体积改变而存储的体积应变能ev,则W可表示为:
W=ev+ef (9)
畸变能是磁畴转动和畴壁位移的主要动力,而体积应变能主要是克服磁化过程中材料的内部缺陷和钉扎,以热能和静电能的形式耗散,对磁化强度的改变影响很小。因此在对应变能求导时,可以忽略体积应变能的影响,直接采用畸变能来分析磁化过程。则材料的应变能也可以表示为:
W=ef (10)
材料的体积应变能可表示为:
式中,σm=1/3(σ1+σ2+σ3)为单元体各面上的平均应力。
将式(8)、式(10)带入式(9)中整理化简求得畸变能为:
复杂应力可等效为一个单向应力,即相当应力σr,根据第四强度理论,σr可以表示为:
同时产生一个等效磁场Hσr。则有效磁场可表示为:
式中,H为外磁场,α为内部耦合场参数,μ0为真空磁导率,λ为磁致伸缩系数。将式(13)带入式(12)可以得到:
对复杂应力状态下材料单位体积的应变能求导可得:
将式(16)带入式(7)可得复杂应力状态下的应力磁化微分方程:
式(17)可以用于分析油气管道在地磁场和外部载荷作用下应力集中区的弱磁信号特征。
磁化强度与相对磁导率μr的关系为:
M=(μr-1)H (18)
将式(18)带入式(17)可得相对磁导率与复杂应力的关系为:
由公式(17)、(19)可见,复杂应力磁化模型将管道运行时的应力状态考虑进去,建立了适用于实际工程的力磁耦合模型,进而可以有效分析实际工程中管道运行时的磁力学特性。
作为一种优选方案,本发明还包括用于验证数学模型准确性的仿真建模部分,仿真建模部分以实际工程应用中的X70钢为研究对象,具体步骤如下;利用有限元软件Ansys,建立复杂应力状态下的管道力磁耦合模型;以管道模型中心为原点建立三维笛卡尔坐标系,模型内表面存在应力集中区,深度方向沿Z轴;在模型空间中添加外磁场,磁场沿X轴方向分布。
作为另一种优选方案,本发明所述管道模型的磁导率根据公式(19)进行设置,通过设置不同的磁导率计算不同应力集中程度下的弱磁信号。
作为另一种优选方案,本发明还包括验证模型准确性的仿真计算部分,包括以下步骤:
1)应力分布云图计算
对管道模型进行静力学仿真,在管道内部施加不同程度的压力,计算不同内压下,管道模型的应力分布云图;观察裂纹附近的应力分布情况,根据云图的颜色深浅确定模型的应力集中区;
2)磁特性计算
对管道模型进行磁学仿真,计算不同应力集中程度下,管道应力集中区的弱磁信号特征,查看弱磁信号计算结果,应力集中区是否出现波峰;查看随着应力集中程度的增大,弱磁信号曲线峰值是否增大;
若弱磁信号在某处出现波峰,且峰值随应力集中程度的增大而增强,则该处存在应力集中情况且弱磁信号随应力集中程度的增大而增强。
作为另一种优选方案,本发明所述应力集中区峰值信号大于材料表面20%时,认为可以从曲线中明显识别出应力集中区。
作为另一种优选方案,本发明所述在管道应力集中区附近每隔3mm记录一个点的弱磁信号特征,曲线15-45mm为应力集中区域。
作为另一种优选方案,本发明所述管道模型包括空气层和钢管,钢管位于空气层中心区域,钢管直径1016mm,长2000mm,壁厚14.5mm,表面存在尺寸30mm×1mm×3mm的裂纹,外加激励磁场沿X轴方向,强度为50μT;计算路径取Y轴方向-30mm-30mm。
作为另一种优选方案,本发明还包括验证复杂应力状态下铁磁性金属构件弱磁信号变化特征的钢管打压实验,该实验试件为一段含人工裂纹的X70型钢管和自制高精度磁通门传感器检测设备;检测设备探头安装在管道裂纹尖端应力集中区,检测应力集中区附近的弱磁信号;采用压力泵向管道注水打压,连接水压传感器监测钢管内水压的变化情况,观察弱磁信号特征随着管道内压的变化趋势,若随着管道内压的不断增大,应力集中区弱磁信号峰值特征逐渐明显,说明与理论模型和仿真结果相同。
作为另一种优选方案,本发明所述X70型钢管长度为6000mm,直径1016mm,壁厚14.5mm。
作为另一种优选方案,本发明所述钢管表面预先制作一条35mm×1mm×3mm的裂纹。
作为另一种优选方案,本发明所述自制高精度磁通门传感器采样频率1.5KHz,磁场检测精度1nT。
其次,本发明所述实验过程中管道内压从0MPa加至3MPa。
优点效果:
本发明有益效果。
本发明采用一种复杂应力状态模型,建立磁力学耦合数学模型,提供一种复杂应力状态下的应力检测定量算法,实现了对管道应力集中程度的判定检测。
附图说明
下面结合附图和具体实施方式对本发明做进一步说明。本发明保护范围不仅局限于以下内容的表述。
图1为本发明应力磁化曲线。
图2为本发明应力磁导率曲线。
图3为本发明管道模型示意图。
图4为本发明磁场模型图。
图5为本发明管道仿真模型应力分布云图。
图6为本发明不同复杂应力下仿真模型的弱磁信号特征。
图7为本发明不同复杂应力下仿真模型弱磁信号峰值的变化特征。
图8为本发明理论模型和仿真模型的力磁变化率对比图。
图9为本发明实验钢管。
图10为本发明检测仪器。
图11为本发明钢管打压实验的弱磁信号特征。
图12为本发明理论模型和打压实验的力磁变化率对比图。
图13为本发明管道内压与管壁所受复杂应力的对应数值。
具体实施方式
一种基于复杂应力状态模型的管道应力内检测方法,该方法步骤如下:
第一步:利用检测设备探头采集弱磁信号;(采集的是管道表面实时的磁场信号)
第二步:将第一步采集的信号输入至构建的复杂应力状态模型,得出管道应力变化;
第三步:根据第二步的管道应力变化识别出管道应力集中区。
第二步的复杂应力状态模型为:
所述复杂应力状态模型构建方法为:设整个***受到外磁场H和应力σ共同作用,应力对***磁化的作用等效为一个外加磁场Hσ,由磁致伸缩系数λ来表示,并引入应力退磁因子Dσ;
应力退磁因子随压应力的增加而线性增加,则有:
式中,Ms为饱和磁化强度,Bs为饱和磁感应强度,λs为饱和磁致伸缩系数,表示为:
式中,υ为泊松比,E为弹性模量,b(σ)是与应力相关的函数,用一分段函数来表示:
式中,b为常数,σ0为初始应力;(数学中的e约等于2.71828)
在外磁场作用下,铁磁性材料的磁化包含两部分,一部分是由畴壁弯曲引起的可逆磁化部分,一部分是由畴壁位移引起的不可逆磁化部分;磁化强度M、可逆磁化强度Mrev、不可逆磁化强度Mirr和无磁滞磁化强度Man的关系如下:
M=Mrev+Mirr (4)
Mrev=c(Man-Mirr) (5)
式中,c为可逆系数;
受外磁场和应力作用的无磁滞磁化强度表示为有效磁场He的函数:
式中,a为材料规划常数;
磁化强度对应变能W的微分表示为:
式中,ξ为系数,与单位体积的能量有关;
由材料力学得知,复杂应力状态下材料单位体积的应变能表示为:
式中,σ1、σ2、σ3表示单元体的三个主应力;
材料的应变能分为因形状改变而存储的畸变能ef和因体积改变而存储的体积应变能ev,则W表示为:
W=ev+ef (9)
在对应变能求导时,忽略体积应变能的影响,直接采用畸变能来分析磁化过程,则材料的应变能也表示为:
W=ef (10)
材料的体积应变能表示为:
式中,σm=1/3(σ1+σ2+σ3)为单元体各面上的平均应力;
将式(8)、式(10)带入式(9)中整理化简求得畸变能为:
复杂应力可等效为一个单向应力,即相当应力σr,根据第四强度理论,σr表示为:
同时产生一个等效磁场Hσr;则有效磁场表示为:
式中,H为外磁场,α为内部耦合场参数,μ0为真空磁导率,λ为磁致伸缩系数;
将式(13)带入式(12)得到:
对复杂应力状态下材料单位体积的应变能求导得:
将式(16)带入式(7)可得复杂应力状态下的应力磁化微分方程:
式(17)用于分析油气管道在地磁场和外部载荷作用下应力集中区的弱磁信号特征。(得出的是磁化强度随复杂应力变化的趋势。)
得出管道应力变化的方法为:
磁化强度与相对磁导率μr的关系为:
M=(μr-1)H (18)
将式(18)带入式(17)得相对磁导率与复杂应力的关系为:
第三步中根据第二步的管道应力变化识别出管道应力集中区的方法如下:
绘制弱磁信号识别曲线,若弱磁信号在某处出现波峰,则该处存在应力集中情况且弱磁信号随应力集中程度的增大而增强。
应力集中区峰值信号大于材料表面20%时,认为可以从曲线中明显识别出应力集中区。
在管道应力集中区附近每隔3mm记录一个点的弱磁信号特征,曲线15-45mm为应力集中区域。
一种基于复杂应力状态模型的管道应力内检测***,所述***包括信号采集模块、应力变化获得模块和应力集中区输出模块;
信号采集模块利用检测设备探头采集弱磁信号;
应力变化获得模块将信号采集模块采集的信号输入至构建的复杂应力状态模型,得出管道应力变化;
应力集中区输出模块根据应力变化获得模块的管道应力变化识别出管道应力集中区。
下面结合附图对本发明做进一步的详细说明:
本发明采用如下技术方案:
材料在应力下会通过磁致伸缩改变内部有效场。设整个***受到外磁场H和应力σ共同作用,应力对***磁化的作用可以等效为一个外加磁场Hσ,由磁致伸缩系数λ来表示。当材料在受到压应力作用时,大多数磁畴会向垂直于应力的方向转动,造成局部磁极密度的增加,这就会导致在外场方向上严重的退磁,以此来减弱材料的磁化;当材料在受到拉应力作用时,大多数磁畴会向平行于应力的方向转动,材料的磁化则不会受到较大影响。在拉压应力状态下材料的磁化是不对称的,因此需要引入应力退磁因子Dσ。应力退磁因子随压应力的增加而线性增加,在拉应力下可以忽略不计,则有:
式中,Ms为饱和磁化强度,Bs为饱和磁感应强度,λs为饱和磁致伸缩系数,可表示为:
式中,υ为泊松比,E为弹性模量,b(σ)是与应力相关的函数,可以用一分段函数来表示:
式中,b为常数,σ0为初始应力。
在外磁场作用下,铁磁性材料的磁化包含两部分,一部分是由畴壁弯曲引起的可逆磁化部分,一部分是由畴壁位移引起的不可逆磁化部分。磁化强度M、可逆磁化强度Mrev、不可逆磁化强度Mirr和无磁滞磁化强度Man的关系如下:
M=Mrev+Mirr (4)
Mrev=c(Man-Mirr) (5)
式中,c为可逆系数。
受外磁场和应力作用的无磁滞磁化强度可表示为有效磁场He的函数:
式中,a为材料规划常数。
磁化强度对应变能W的微分可表示为:
式中,ξ为系数,与单位体积的能量有关。
实际工程中,在役油气管道在运行时受到内压产生来自各个方向、复杂不定的拉伸应力,此时管道处于复杂应力状态。
由材料力学可知,复杂应力状态下材料单位体积的应变能可以表示为:
式中,σ1、σ2、σ3表示单元体的三个主应力。
材料的应变能分为因形状改变而存储的畸变能ef和因体积改变而存储的体积应变能ev,则W可表示为:
W=ev+ef (9)
畸变能是磁畴转动和畴壁位移的主要动力,而体积应变能主要是克服磁化过程中材料的内部缺陷和钉扎,以热能和静电能的形式耗散,对磁化强度的改变影响很小。因此在对应变能求导时,可以忽略体积应变能的影响,直接采用畸变能来分析磁化过程。则材料的应变能也可以表示为:
W=ef (10)
材料的体积应变能可表示为:
式中,σm=1/3(σ1+σ2+σ3)为单元体各面上的平均应力。
将式(8)、式(10)带入式(9)中整理化简求得畸变能为:
复杂应力可等效为一个单向应力,即相当应力σr,根据第四强度理论,σr可以表示为:
同时产生一个等效磁场Hσr。则有效磁场可表示为:
式中,H为外磁场,α为内部耦合场参数,μ0为真空磁导率,λ为磁致伸缩系数。
将式(13)带入式(12)可以得到:
对复杂应力状态下材料单位体积的应变能求导可得:
将式(16)带入式(7)可得复杂应力状态下的应力磁化微分方程:
式(17)可以用于分析油气管道在地磁场和外部载荷作用下应力集中区的弱磁信号特征。
磁化强度与相对磁导率μr的关系为:
M=(μr-1)H (18)
将式(18)带入式(17)可得相对磁导率与复杂应力的关系为:
由公式(17)、(19)可见,复杂应力磁化模型将管道运行时的应力状态考虑进去,建立了适用于实际工程的力磁耦合模型,进而可以有效分析实际工程中管道运行时的磁力学特性。
本发明还包括计算数学模型磁力学特性以及用于验证数学模型准确性的仿真建模部分,仿真建模部分以实际工程应用中的X70钢为研究对象,建立磁力学模型,具体步骤如下:
1)数学模型计算
式(7)描述了应变能与磁信号的对应关系,将力学特性与磁学特性联系起来;将复杂应力状态下的应变能带入式(7)得到复杂应力磁力学模型,将式(18)带入式(17)得到复杂应力状态下的相对磁导率,如图1,2所示;
图1中X轴为应力值,其中正半轴表示拉应力,负半轴表示压应力,可以看出:在拉应力和压应力阶段,复杂应力状态模型具有很好的对称性,整体变化平稳。由图1和图2可以看出,在拉应力阶段,磁化曲线与磁导率曲线都呈现线性增加的变化趋势。
2)仿真模型计算
利用有限元软件Ansys,建立复杂应力状态下的管道力磁耦合模型,管道表面存在一条轴向裂纹,模型示意图如图3所示;在模型空间中添加外磁场,磁场沿X轴方向分布,磁场模型图如图4所示;管道模型的磁导率根据公式(19)进行设置,通过设置不同的磁导率计算不同应力集中程度下的弱磁信号。
对管道模型进行静力学仿真,在管道内部施加不同程度的压力,计算不同内压下,管道模型的应力分布云图;观察裂纹附近的应力分布情况,根据云图的颜色深浅确定模型的应力集中区,如图5所示。可以看出,管道裂纹的尖端及底部受应力较大,存在应力集中,取裂纹尖端为研究区域。
钢管的周向拉力与内部压力之间满足其中,FT2为周向拉力,单位为MPa/mm;P为管道承受的内压,单位为MPa;d为管道直径,单位为mm;T为管道壁厚,单位为mm。管壁所受复杂应力可分解为轴向应力和周向应力,其中轴向应力大小为周向应力的1/2。管道内压与管壁所受复杂应力的对应数值如下表所示(下表为本发明管道内压与管壁所受复杂应力的对应数值);
对管道模型进行磁学仿真,沿管道外壁扫描裂纹尖端应力集中区,取Y轴方向-30mm-30mm作为检测路径,分析结果如图6,7所示。
图6中X轴为扫描路径位移,其中15-45mm为应力集中区,Y轴为复杂应力值,Z轴为弱磁信号强度,可以看出:随应力的增大,裂纹尖端应力集中区的弱磁信号增强,弱磁信号峰值线性增加,即弱磁信号与复杂应力具有一一对应的线性关系。
计算图1拉应力阶段和图7的曲线斜率,得到理论模型和仿真模型的力磁变化率k,对比结果如图8所示,理论计算结果与仿真计算结果对比得出,仿真模型的力磁变化率与理论模型的力磁变化率误差相差3%,证明了理论模型的正确性。所述管道模型计算中的内压调节范围为0-3Mpa,间隔为0.5MPa。
所述在管道应力集中区附近每隔3mm记录一个点的弱磁信号特征。
所述管道模型包括空气层和钢管,钢管位于空气层中心区域,钢管直径1016mm,长2000mm,壁厚14.5mm,表面存在尺寸30mm×1mm×3mm的裂纹,
所述模型外加激励磁场沿X轴方向,强度为50μT。
本发明还包括验证复杂应力状态下铁磁性金属构件弱磁信号变化规律的钢管打压实验。该实验试件为一段含人工裂纹的X70型钢管和自制高精度磁通门传感器检测设备,实验钢管如图9所示,检测仪器如图10所示;
检测设备探头安装在管道裂纹尖端应力集中区,检测应力集中区附近的弱磁信号;采用压力泵向管道注水打压,连接水压传感器监测钢管内水压的变化情况,观察弱磁信号特征随着管道内压的变化趋势。
钢管打压实验中,应力集中区的弱磁信号变化特征如图11所示,可以看出,打压过程中,管道应力集中区的弱磁信号与应力呈一一对应的线性关系。计算图1拉力阶段和图11曲线的斜率,得到理论模型和钢管打压实验的力磁变化率k,对比结果如图12所示,与打压实验相比,复杂应力状态理论模型的力磁变化率误差极小,误差率约为0.485%,说明复杂应力状态模型适用于实际管道应力检测。
所述X70型钢管长度为6000mm,直径1016mm,壁厚14.5mm。
所述钢管表面预先制作一条35mm×1mm×3mm的裂纹。
所述自制高精度磁通门传感器采样频率1.5KHz,磁场检测精度1nT。
所述管道内压从0MPa加至3MPa,间隔为0.5MPa。
可以理解的是,以上关于本发明的具体描述,仅用于说明本发明而并非受限于本发明实例所描述的技术方案,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换,以达到相同的技术效果;只要满足使用需要,都在本发明的保护范围之内。
Claims (4)
1.一种基于复杂应力状态模型的管道应力内检测方法,其特征在于:该方法步骤如下:
第一步:利用检测设备探头采集弱磁信号;
第二步:将第一步采集的信号输入至构建的复杂应力状态模型,得出管道应力变化;
第三步:根据第二步的管道应力变化识别出管道应力集中区;
第二步的复杂应力状态模型为:
所述复杂应力状态模型构建方法为:设整个***受到外磁场H和应力σ共同作用,应力对***磁化的作用等效为一个外加磁场Hσ,由磁致伸缩系数λ来表示,并引入应力退磁因子Dσ;
应力退磁因子随压应力的增加而线性增加,则有:
式中,Ms为饱和磁化强度,Bs为饱和磁感应强度,λs为饱和磁致伸缩系数,表示为:
式中,υ为泊松比,E为弹性模量,b(σ)是与应力相关的函数,用一分段函数来表示:
式中,b为常数,σ0为初始应力;
在外磁场作用下,铁磁性材料的磁化包含两部分,一部分是由畴壁弯曲引起的可逆磁化部分,一部分是由畴壁位移引起的不可逆磁化部分;磁化强度M、可逆磁化强度Mrev、不可逆磁化强度Mirr和无磁滞磁化强度Man的关系如下:
M=Mrev+Mirr (4)
Mrev=c(Man-Mirr) (5)
式中,c为可逆系数;
受外磁场和应力作用的无磁滞磁化强度表示为有效磁场He的函数:
式中,a为材料规划常数;
磁化强度对应变能W的微分表示为:
式中,ξ为系数,与单位体积的能量有关;
由材料力学得知,复杂应力状态下材料单位体积的应变能表示为:
式中,σ1、σ2、σ3表示单元体的三个主应力;
材料的应变能分为因形状改变而存储的畸变能ef和因体积改变而存储的体积应变能ev,则W表示为:
W=ev+ef (9)
在对应变能求导时,忽略体积应变能的影响,直接采用畸变能来分析磁化过程,则材料的应变能也表示为:
W=ef (10)
材料的体积应变能表示为:
式中,σm=1/3(σ1+σ2+σ3)为单元体各面上的平均应力;
将式(8)、式(10)带入式(9)中整理化简求得畸变能为:
复杂应力可等效为一个单向应力,即相当应力σr,根据第四强度理论,σr表示为:
同时产生一个等效磁场Hσr;则有效磁场表示为:
式中,H为外磁场,α为内部耦合场参数,μ0为真空磁导率,λ为磁致伸缩系数;
将式(13)带入式(12)得到:
对复杂应力状态下材料单位体积的应变能求导得:
将式(16)带入式(7)可得复杂应力状态下的应力磁化微分方程:
式(17)用于分析油气管道在地磁场和外部载荷作用下应力集中区的弱磁信号特征。
磁化强度与相对磁导率μr的关系为:
M=(μr-1)H (18)
将式(18)带入式(17)得相对磁导率与复杂应力的关系为:
2.根据权利要求1所述的一种基于复杂应力状态模型的管道应力内检测方法,其特征在于:
第三步中根据第二步的管道应力变化识别出管道应力集中区的方法如下:
绘制弱磁信号识别曲线,若弱磁信号在某处出现波峰,则该处存在应力集中情况且弱磁信号随应力集中程度的增大而增强。
3.根据权利要求1所述的一种基于复杂应力状态模型的管道应力内检测方法,其特征在于:应力集中区峰值信号大于材料表面20%时,认为可以从曲线中明显识别出应力集中区。
4.根据权利要求1所述的一种基于复杂应力状态模型的管道应力内检测方法,其特征在于:在管道应力集中区附近每隔3mm记录一个点的弱磁信号特征,曲线15-45mm为应力集中区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911085775.2A CN110793852B (zh) | 2019-11-08 | 2019-11-08 | 一种基于复杂应力状态模型的管道应力内检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911085775.2A CN110793852B (zh) | 2019-11-08 | 2019-11-08 | 一种基于复杂应力状态模型的管道应力内检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110793852A CN110793852A (zh) | 2020-02-14 |
CN110793852B true CN110793852B (zh) | 2022-07-15 |
Family
ID=69443295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911085775.2A Active CN110793852B (zh) | 2019-11-08 | 2019-11-08 | 一种基于复杂应力状态模型的管道应力内检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110793852B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112504112B (zh) * | 2020-12-01 | 2022-04-05 | 西南石油大学 | 一种山区管道应变监测安全管环与方法 |
CN112562799B (zh) * | 2020-12-11 | 2024-04-26 | 西安建筑科技大学 | 一种磁致伸缩材料的力磁本构模型及参数标定方法 |
CN113092572B (zh) * | 2021-04-09 | 2024-01-26 | 中国矿业大学 | 一种管道达到磁饱和所需磁化强度的确定方法 |
CN113375846B (zh) * | 2021-05-17 | 2021-12-03 | 西南石油大学 | 一种管道轴向应力快速检测装置与方法 |
CN114544043B (zh) * | 2022-01-29 | 2023-07-07 | 中国石油大学(北京) | 多钢种应力检测装置及方法 |
CN117571815B (zh) * | 2023-11-15 | 2024-06-21 | 中磁数智(北京)科技有限公司 | 一种供热直埋管道的预留支管位置的弱磁检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009222425A (ja) * | 2008-03-13 | 2009-10-01 | Katsuyuki Kida | 非破壊検査装置及び非破壊検査方法 |
CN103376289A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油天然气股份有限公司 | 一种长输油气管道应力集中区域的弱磁内检测方法 |
CN105203629A (zh) * | 2015-09-24 | 2015-12-30 | 三峡大学 | 一种磁性金属构件表面应力集中区及微裂纹的磁探测方法 |
CN109974792A (zh) * | 2019-03-22 | 2019-07-05 | 西安石油大学 | 基于磁力耦合效应的管线无损应力检测试验***及方法 |
CN110095525A (zh) * | 2019-05-23 | 2019-08-06 | 合肥工业大学 | 一种基于内聚力-磁机械耦合模型的界面裂纹评估方法 |
-
2019
- 2019-11-08 CN CN201911085775.2A patent/CN110793852B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009222425A (ja) * | 2008-03-13 | 2009-10-01 | Katsuyuki Kida | 非破壊検査装置及び非破壊検査方法 |
CN103376289A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油天然气股份有限公司 | 一种长输油气管道应力集中区域的弱磁内检测方法 |
CN105203629A (zh) * | 2015-09-24 | 2015-12-30 | 三峡大学 | 一种磁性金属构件表面应力集中区及微裂纹的磁探测方法 |
CN109974792A (zh) * | 2019-03-22 | 2019-07-05 | 西安石油大学 | 基于磁力耦合效应的管线无损应力检测试验***及方法 |
CN110095525A (zh) * | 2019-05-23 | 2019-08-06 | 合肥工业大学 | 一种基于内聚力-磁机械耦合模型的界面裂纹评估方法 |
Non-Patent Citations (3)
Title |
---|
拉压不同应力对磁记忆信号的影响及机理;冷建成等;《工程科学学报》;20180531;第3节 * |
长输油气管道弱磁应力内检测技术;闵希华等;《机械工程学报》;20170630;第3-4节 * |
闵希华等.长输油气管道弱磁应力内检测技术.《机械工程学报》.2017, * |
Also Published As
Publication number | Publication date |
---|---|
CN110793852A (zh) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110793852B (zh) | 一种基于复杂应力状态模型的管道应力内检测方法 | |
CN108896230B (zh) | 一种基于有限元的螺栓紧固力超声检测及关键检测参数确定方法 | |
Shen et al. | A novel method for determining surface residual stress components and their directions in spherical indentation | |
Liu et al. | Study on internal detection in oil–gas pipelines based on complex stress magnetomechanical modeling | |
CN110672001A (zh) | 一种铁磁性材料表面非铁磁材料厚度测量方法及装置 | |
Tamhane et al. | Feature engineering of time-domain signals based on principal component analysis for rebar corrosion assessment using pulse eddy current | |
CN105352433A (zh) | 一种船体典型焊接结构表面裂纹深度及形状测量装置及其方法 | |
Zhao et al. | Coiled tubing wall thickness evaluation system using pulsed alternating current field measurement technique | |
CN109870258B (zh) | 一种平面任意残余应力的仪器化球形压入检测方法 | |
CN112630023A (zh) | 基于热致磁变原理的铁磁性金属材料轴向应力检测方法 | |
CN116184279A (zh) | 一种弱磁信号定量化的方法 | |
CN108229052B (zh) | 基于flapw算法的磁记忆信号特征研究方法 | |
CN116087319A (zh) | 一种漏磁信号特征的量化方法 | |
Fujimoto et al. | Stress intensity factor measurement of cracks using a piezoelectric element | |
CN114740081B (zh) | 一种应力沿深度分布的微磁检测方法 | |
Zhou et al. | Research on the influence of composite defects on quantitative pattern in MFL detection of tubing | |
Alimey et al. | Gradient feature extract for the quantification of complex defects using topographic primal sketch in magnetic flux leakage | |
Zhang et al. | Quantitative study on cross-section damage of steel wire rope based on magnetic signal characteristics under weak magnetic excitation | |
CN102645484A (zh) | 一种金属材料的电磁声发射无损检测装置 | |
CN108280301B (zh) | 一种磁记忆信号特征研究方法 | |
Zhang et al. | Weak Magnetic Flux Leakage Detection and Numerical Method Under Inclined Pressure Conditions | |
Li et al. | Calculation and experimental verification of force-magnetic coupling model of magnetised rail based on density functional theory | |
Gang et al. | Signal Processing and Defect Characterization of Pulsed Eddy Current Testing Based on Cross-Correlation | |
CN106777865B (zh) | 岩体渗透性Dice相似度量预测分类方法 | |
Tang et al. | The Analysis of Weak Magnetic Inspection of Standard Test Block Defects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |