CN110646184B - 一种船舶动力管道疲劳应力的监测方法 - Google Patents
一种船舶动力管道疲劳应力的监测方法 Download PDFInfo
- Publication number
- CN110646184B CN110646184B CN201910901092.3A CN201910901092A CN110646184B CN 110646184 B CN110646184 B CN 110646184B CN 201910901092 A CN201910901092 A CN 201910901092A CN 110646184 B CN110646184 B CN 110646184B
- Authority
- CN
- China
- Prior art keywords
- stress
- cycle
- pipeline
- calculating
- monitoring point
- 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
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种船舶动力管道疲劳应力的监测方法,涉及动力管道的安全控制技术领域,包括以下步骤:建立动力管道的应力分析模型,生成管道在单位载荷情况下,所有监测点的应力对照表;获取船舶航行过程中,管道受到的低周循环载荷和高周循环载荷的时间历程曲线;计算并绘制监测点的低周应力循环的累积频数表;计算管道的等效重力载荷因子;计算监测点的高周应力循环的累积频数表;合并监测点的低周应力循环的累积频数表和高周应力循环的累积频数表,得到完整的应力循环的累积频数表,以完成对船舶管道的监测点进行疲劳分析与校核,本发明能够减少传感器的设置,并避免传感器直接贴在危险点测量,间接完成对船舶动力管道任意区域的应力监测。
Description
技术领域
本发明涉及动力管道的安全控制技术领域,具体涉及一种船舶动力管道疲劳应力的监测方法。
背景技术
动力管道是能源系统中重要组成部分,承担着汽、水、油等工质输送和能量传递的任务,对能源系统来说至关重要。一般的,动力管道内部的工质存在较高压力,在交变载荷的作用下,一旦结构发生疲劳破坏,会导致大量工质泄漏,可能产生爆炸、火灾、中毒等严重的安全事故,造成巨大的经济损失。特别对于船舶能源系统,由于远离陆地,外部的救援力量无法快速到达事故区域,一旦船舶能源系统发生动力管道泄漏事故,产生的事故后果更为严重。为此,对船舶能源系统复杂动力管道上可能发生破坏的区域进行疲劳应力监测,对疲劳破裂危险进行预报显得尤其重要。
目前,现有的应力测量主要采用在测量点布置应力传感器的方式,测量比较直接。但是对于船舶管道系统,这种方法存在以下问题:一是船舶能源系统的动力管道复杂庞大,存在的可能发生疲劳破坏的区域较多,实现整个管路的应力监测需要设置大量传感器以及采集数据的导线,投资与维护成本高,布置空间也不允许;二是船舶舱室环境恶劣,高温、高湿、振动等环境并不适合应力传感器长期工作。
发明内容
针对现有技术中存在的缺陷,本发明的目的在于提供一种船舶动力管道疲劳应力的监测方法,能够减少传感器的设置,并避免传感器直接贴在危险点测量,间接完成对船舶动力管道任意区域的应力监测。
为达到以上目的,本发明采取的技术方案是:
提供一种船舶动力管道疲劳应力的监测方法,包括以下步骤:
建立动力管道的ANSYS应力分析模型,分别计算在各载荷形式的单位载荷下,监测点产生的应力值,并完成对所有监测点的计算,生成管道在单位载荷情况下,所有监测点的应力对照表,所述载荷形式包括压力、温度、摇摆和振荡;
在船舶不同用途的动力管道上均设置一压力传感器和温度传感器,并在船舶的船体上设置加速度传感器;
通过压力传感器、温度传感器和加速度传感器获取船舶航行过程中,管道受到的低周循环载荷和高周循环载荷的时间历程曲线,所述低周循环载荷包括压力和温度,所述高周循环载荷包括摇摆和振荡;
根据应力对照表和低周循环载荷的时间历程曲线,计算并绘制监测点的低周应力循环的累积频数表;
根据高周循环载荷的时间历程曲线和管道中心距船舶中心位置的距离,计算管道的等效重力载荷因子;
根据应力对照表和管道的等效重力载荷因子,计算监测点的高周应力循环的累积频数表;
合并监测点的低周应力循环的累积频数表和高周应力循环的累积频数表,得到完整的应力循环的累积频数表,以完成对船舶管道的监测点进行疲劳分析与校核。
在上述技术方案的基础上,分别计算在各载荷形式的单位载荷下,监测点产生的应力值时,具体包括以下步骤:
计算管道温度相对初始状态升高1℃时,监测点产生的等效应力值σ0t;
计算管道内压相对初始状态升高0.1MPa时,监测点产生的等效应力值σ0y;
计算管道横向方向受到1g加速度时,监测点产生的等效应力值σ0r;
计算管道垂向方向受到1g加速度时,监测点产生的等效应力值σ0e;
计算管道纵向方向受到1g加速度时,监测点产生的等效应力值σ0p,其中,g为重力加速度。
在上述技术方案的基础上,计算并绘制监测点的低周应力循环的累积频数表时,具体包括以下步骤:
根据应力对照表,得到监测点实时温度的σ0t和实时压力的σ0y数据;
计算监测点在低周循环载荷下的实时载荷应力值:
σ(t)=P(t)×σ0y+T(t)×σ0t
其中,P(t)为实时压力,T(t)为实时温度;
根据实时载荷应力值绘制低周应力循环历程曲线;
通过雨流计数法对低周应力循环历程曲线进行处理,得到监测点的低周应力循环累积频数表。
在上述技术方案的基础上,计算管道的等效重力载荷因子时,具体包括以下步骤:
计算管道中心距船舶中心位置的距离;
根据以上数据,计算管道的等效重力载荷因子。
在上述技术方案的基础上,管道中心距船舶中心位置的距离包括管道中心距船舶中心位置的纵向距离X、横向距离Y和垂向距离Z;
计算管道的等效重力载荷因子时,具有包括以下步骤:
在上述技术方案的基础上,计算监测点的高周应力循环的累积频数表时,具体包括以下步骤,
根据应力对照表,得到监测点的σ0r、σ0e和σ0p数据;
根据监测点的σ0r、σ0e和σ0p数据,计算监测点在高周循环载荷下的最大应力值,
根据振荡的时间历程曲线,提取垂向振荡频数Ne,得到监测点最大应力值和频数对应的高周应力循环的累积频数表。
在上述技术方案的基础上,所述高周循环载荷下的最大应力值包括纵向最大应力值、横向最大应力值和垂向最大应力值,监测点的纵向最大应力值σXmax的计算公式为:
监测点的横向最大应力值σYmax的计算公式为:
监测点的垂向最大应力值σZmax的计算公式为:
与现有技术相比,本发明的优点在于:
本发明的一种船舶动力管道疲劳应力的监测方法,由于具有同样用途的管道内的工质压力和工质温度是相同的,且管道是固定于船舶的船体上的,因此,在进行测量时,只需要设置少量的传感器进行测量不同用途的管道内的工质压力和工质温度,以及管道的摇摆和振荡载荷,相对于在每个监测点设置应力传感器来说,减少了传感器和导线的设置,进而降低了投资与维护的成本;同时,通过计算低周应力循环和高周应力循环的累积频数表,得到完整的管道应力循环的累积频数表,进而完成对船舶动力管道的疲劳分析,间接完成了对船舶动力管道任意区域的应力监测,避免传感器直接贴在危险点测量,提高船舶海上航行的安全性。
附图说明
图1为本发明实施例中船舶动力管道疲劳应力的监测方法的流程图。
具体实施方式
以下结合附图及实施例对本发明作进一步详细说明。
参见图1所示,本发明实施例提供一种船舶动力管道疲劳应力的监测方法,其特征在于,包括以下步骤:
建立动力管道的ANSYS应力分析模型,分别计算在各载荷形式的单位载荷下,监测点产生的应力值,并完成对所有监测点的计算,生成管道在单位载荷情况下,所有监测点的应力对照表,载荷形式包括压力、温度、摇摆和振荡,分别计算在各载荷形式的单位载荷下,监测点产生的应力值时,具体包括以下步骤:
计算管道温度相对初始状态升高1℃时,监测点产生的等效应力值σ0t;
计算管道内压相对初始状态升高0.1MPa时,监测点产生的等效应力值σ0y;
计算管道横向方向受到1g加速度时,监测点产生的等效应力值σ0r;
计算管道垂向方向受到1g加速度时,监测点产生的等效应力值σ0e;
计算管道纵向方向受到1g加速度时,监测点产生的等效应力值σ0p,其中,g为重力加速度;
生成的每个监测点的应力对照表如下:
在船舶不同用途的动力管道上均设置一压力传感器和温度传感器,并在船舶的船体上设置加速度传感器;
通过压力传感器、温度传感器和加速度传感器获取提取船舶航行过程中,管道受到的低周循环载荷和高周循环载荷的时间历程曲线,所述低周循环载荷包括压力和温度,其中,测量的压力和温度为管道工质压力和工质温度,所述高周循环载荷包括摇摆和振荡;
根据应力对照表和低周循环载荷的时间历程曲线,计算并绘制监测点的低周应力循环的累积频数表,计算并绘制监测点的低周应力循环的累积频数表时,具体包括以下步骤:
根据应力对照表,得到监测点实时温度的σ0t和实时压力的σ0y数据;
计算监测点在低周循环载荷下的实时载荷应力值:
σ(t)=P(t)×σ0y+T(t)×σ0t
其中,P(t)为实时压力,T(t)为实时温度;
根据实时载荷应力值绘制低周应力循环历程曲线;
通过雨流计数法对低周应力循环历程曲线进行处理,得到监测点的低周应力循环的累积频数表;
根据高周循环载荷的时间历程曲线和管道中心距船舶中心位置的距离,计算管道的等效重力载荷因子,计算管道的等效重力载荷因子时,具体包括以下步骤:
计算管道中心距船舶中心位置的距离;
根据以上数据,计算管道的等效重力载荷因子;
其中,管道中心距船舶中心位置的距离包括管道中心距船舶中心位置的纵向距离X、横向距离Y和垂向距离Z;
计算管道的等效重力载荷因子时,具有包括以下步骤:
根据应力对照表和管道的等效重力载荷因子,计算监测点的高周应力循环的累积频数表,计算监测点的高周应力循环的累积频数表时,具体包括以下步骤,
根据应力对照表,得到监测点的σ0r、σ0e和σ0p数据;
根据监测点的σ0r、σ0e和σ0p数据,计算监测点在高周循环载荷下的最大应力值,
根据振荡的时间历程曲线,提取垂向振荡频数Ne,得到监测点最大应力值和频数对应的高周应力循环的累积频数表,其中,最大应力值对应为幅值;
其中,所述高周循环载荷下的最大应力值包括纵向最大应力值、横向最大应力值和垂向最大应力值,监测点的纵向最大应力值σXmax满足以下公式:
监测点的横向最大应力值σYmax满足以下公式:
监测点的垂向最大应力值σZmax满足以下公式:
进一步的,高周应力循环的累积频数表为:
幅值 | σ<sub>Xmax</sub> | σ<sub>Ymax</sub> | σ<sub>Zmax</sub> |
频数 | N<sub>e</sub> | N<sub>e</sub> | N<sub>e</sub> |
合并监测点的低周应力循环的累积频数表和高周应力循环的累积频数表,得到完整的应力循环的累积频数表,并通过疲劳与断裂理论根据完整的应力循环的累积频数表完成对船舶管道的监测点进行疲劳分析与校核。
建立动力管道的ANSYS应力分析模型,生成所有监测点的应力对照表,并提取低周循环载荷和高周循环载荷的时间历程曲线,通过计算低周应力循环和高周应力循环的累积频数表,得到完整的管道应力循环的累积频数表,能够实现船舶动力管道任意位置(任意监测点)的应力疲劳应力监测,并用于完成对船舶动力管道的疲劳分析,提高船舶海上航行的安全性,具有较高的经济效益;同时,本实施中监测温度、压力、摇摆等参数,通过间接计算实现管道上任意位置点的应力载荷的测量,相比于在危险点直接贴应力传感器的方法,提高船舶海上航行的安全性,具有测量数据量小,装置成本低,运行可靠性高,耐受性强,以及实施便捷的特点。
同时,由于具有同样用途的管道内的工质压力和工质温度是相同的,且管道是固定于船舶的船体上的,因此,各部位的管道的摇摆和振荡是一致的,则在进行测量时,只需要设置少量的传感器进行测量不同用途的管道内的工质压力和工质温度,相对于在每个监测点设置应力传感器来说,减少了传感器和导线的设置,进而降低了投资与维护的成本,缩小了布置的空间,其中,测量温度和压力的传感器,为船舶本身的设置,在进行检测管道疲劳应力时,不需要额外增设传感器,进一步的减小了投资成本。
本发明不局限于上述实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (5)
1.一种船舶动力管道疲劳应力的监测方法,其特征在于,包括以下步骤:
建立动力管道的ANSYS应力分析模型,分别计算在各载荷形式的单位载荷下,监测点产生的应力值,并完成对所有监测点的计算,生成管道在单位载荷情况下,所有监测点的应力对照表,所述载荷形式包括压力、温度、摇摆和振荡;
在船舶不同用途的动力管道上均设置一压力传感器和温度传感器,并在船舶的船体上设置加速度传感器;
通过压力传感器、温度传感器和加速度传感器获取船舶航行过程中,管道受到的低周循环载荷和高周循环载荷的时间历程曲线,所述低周循环载荷包括压力和温度,所述高周循环载荷包括摇摆和振荡;
根据应力对照表和低周循环载荷的时间历程曲线,计算并绘制监测点的低周应力循环的累积频数表;
根据高周循环载荷的时间历程曲线和管道中心距船舶中心位置的距离,计算管道的等效重力载荷因子;
根据应力对照表和管道的等效重力载荷因子,计算监测点的高周应力循环的累积频数表;
合并监测点的低周应力循环的累积频数表和高周应力循环的累积频数表,得到完整的应力循环的累积频数表,以完成对船舶管道的监测点进行疲劳分析与校核;
分别计算在各载荷形式的单位载荷下,监测点产生的应力值时,具体包括以下步骤:
计算管道温度相对初始状态升高1℃时,监测点产生的等效应力值σ0t;
计算管道内压相对初始状态升高0.1MPa时,监测点产生的等效应力值σ0y;
计算管道横向方向受到1g加速度时,监测点产生的等效应力值σ0r;
计算管道垂向方向受到1g加速度时,监测点产生的等效应力值σ0e;
计算管道纵向方向受到1g加速度时,监测点产生的等效应力值σ0p,其中,g为重力加速度;
计算并绘制监测点的低周应力循环的累积频数表时,具体包括以下步骤:
根据应力对照表,得到监测点实时温度的σ0t和实时压力的σ0y数据;
计算监测点在低周循环载荷下的实时载荷应力值:
σ(t)=P(t)×σ0y+T(t)×σ0t
其中,P(t)为实时压力,T(t)为实时温度;
根据实时载荷应力值绘制低周应力循环历程曲线;
通过雨流计数法对低周应力循环历程曲线进行处理,得到监测点的低周应力循环累积频数表。
4.如权利要求3所述的一种船舶动力管道疲劳应力的监测方法,其特征在于,计算监测点的高周应力循环的累积频数表时,具体包括以下步骤,
根据应力对照表,得到监测点的σ0r、σ0e和σ0p数据;
根据监测点的σ0r、σ0e和σ0p数据,计算监测点在高周循环载荷下的最大应力值,
根据振荡的时间历程曲线,提取垂向振荡频数Ne,得到监测点最大应力值和频数对应的高周应力循环的累积频数表。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910901092.3A CN110646184B (zh) | 2019-09-23 | 2019-09-23 | 一种船舶动力管道疲劳应力的监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910901092.3A CN110646184B (zh) | 2019-09-23 | 2019-09-23 | 一种船舶动力管道疲劳应力的监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110646184A CN110646184A (zh) | 2020-01-03 |
CN110646184B true CN110646184B (zh) | 2021-08-31 |
Family
ID=68992355
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910901092.3A Active CN110646184B (zh) | 2019-09-23 | 2019-09-23 | 一种船舶动力管道疲劳应力的监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110646184B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111625980A (zh) * | 2020-05-28 | 2020-09-04 | 中国船舶工业集团公司第七0八研究所 | 一种fpso结构疲劳寿命计算方法 |
CN114199438B (zh) * | 2021-12-10 | 2023-07-21 | 哈尔滨工程大学 | 一种基于船舶结构监测数据分析的载荷识别方法 |
CN114216599B (zh) * | 2021-12-10 | 2023-07-21 | 哈尔滨工程大学 | 一种基于传感器数据分析的船舶监测结构载荷识别方法 |
CN114199440B (zh) * | 2021-12-10 | 2023-07-21 | 哈尔滨工程大学 | 一种船舶加筋板结构应力监测数据的转换处理方法 |
CN114199436B (zh) * | 2021-12-10 | 2023-07-21 | 哈尔滨工程大学 | 一种基于传感器数据的船舶监测结构整体应力状态快速获取方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007086049A (ja) * | 2005-09-23 | 2007-04-05 | Jiaotong Univ | 地面の変位をモニターするファイバブラッググレーティングにより検知するセグメント化された歪み計 |
CN102567567A (zh) * | 2011-11-15 | 2012-07-11 | 北京宇航系统工程研究所 | 一种基于有限元分析的管路随机振动疲劳寿命分析方法 |
CN107559595A (zh) * | 2017-08-04 | 2018-01-09 | 武汉理工大学 | 一种船舶管道振动应急处理系统 |
CN107742039A (zh) * | 2017-10-31 | 2018-02-27 | 合肥通用机械研究院 | 一种基于悬跨涡激振动疲劳寿命的海底管道设计方法 |
CN109885920A (zh) * | 2019-02-13 | 2019-06-14 | 西安电子科技大学 | 一种航空发动机涡轮叶片的高低周复合疲劳寿命预测方法 |
CN110108455A (zh) * | 2019-04-09 | 2019-08-09 | 中广核工程有限公司 | 一种核电厂小型管道振动应力测量及疲劳寿命评估方法 |
-
2019
- 2019-09-23 CN CN201910901092.3A patent/CN110646184B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007086049A (ja) * | 2005-09-23 | 2007-04-05 | Jiaotong Univ | 地面の変位をモニターするファイバブラッググレーティングにより検知するセグメント化された歪み計 |
CN102567567A (zh) * | 2011-11-15 | 2012-07-11 | 北京宇航系统工程研究所 | 一种基于有限元分析的管路随机振动疲劳寿命分析方法 |
CN107559595A (zh) * | 2017-08-04 | 2018-01-09 | 武汉理工大学 | 一种船舶管道振动应急处理系统 |
CN107742039A (zh) * | 2017-10-31 | 2018-02-27 | 合肥通用机械研究院 | 一种基于悬跨涡激振动疲劳寿命的海底管道设计方法 |
CN109885920A (zh) * | 2019-02-13 | 2019-06-14 | 西安电子科技大学 | 一种航空发动机涡轮叶片的高低周复合疲劳寿命预测方法 |
CN110108455A (zh) * | 2019-04-09 | 2019-08-09 | 中广核工程有限公司 | 一种核电厂小型管道振动应力测量及疲劳寿命评估方法 |
Non-Patent Citations (2)
Title |
---|
海上浮动核电站管道系统摇摆载荷分析;何风 等;《机械与电子》;20160229;第34卷(第2期);正文第1-4节 * |
船舶动力装置高温高压管道破前漏分析;白凡 等;《工程力学》;20180331;第35卷(第3期);正文第1-3节 * |
Also Published As
Publication number | Publication date |
---|---|
CN110646184A (zh) | 2020-01-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110646184B (zh) | 一种船舶动力管道疲劳应力的监测方法 | |
CN102914427B (zh) | 一种多轴随机载荷下疲劳损伤评估方法 | |
CN107324214B (zh) | 海洋平台起重机智能状态监测方法 | |
CN110738753B (zh) | 一种船舶动力管道疲劳损伤预测系统及预测方法 | |
CN107064315B (zh) | 基于小波分析的叶片裂纹位置及深度识别方法 | |
CN106441851A (zh) | 一种机械零件疲劳寿命的检测方法 | |
CN106656669A (zh) | 一种基于阈值自适应设置的设备参数异常检测系统及方法 | |
CN111625980A (zh) | 一种fpso结构疲劳寿命计算方法 | |
CN104807586A (zh) | 海上浮式风机系泊系统断裂失效预报方法与系统 | |
CN113374652A (zh) | 一种风力发电机组寿命评估方法 | |
CN105241589A (zh) | 一种机器人手臂应变测试数据处理方法 | |
CN103471840B (zh) | 箱型梁垂向弯矩试验机构 | |
CN107607282A (zh) | 液货船舶碰撞实验装置及其实验方法 | |
CN109614689B (zh) | 一种确定船舶结构在组合载荷作用下长期疲劳寿命的方法 | |
CN110390173B (zh) | 考虑剩余强度退化的千米深井提升机时变可靠性评估方法 | |
CN107808046A (zh) | 一种船体梁爆炸冲击动弯矩确定方法 | |
CN104880331A (zh) | 起重机风载模拟测试使用方法 | |
CN204831675U (zh) | 海上浮式风机系泊系统断裂失效预报系统 | |
CN113591272B (zh) | 基于预应力修正的复杂管路结构疲劳损伤评估方法和系统 | |
CN103678775B (zh) | 一种输送管动强度分析方法 | |
CN104880293A (zh) | 起重机风载模型建立方法 | |
CN110287619B (zh) | 一种多轴随机载荷疲劳寿命预测方法 | |
CN105527114A (zh) | 起重机风载实测方法 | |
CN112761818B (zh) | 一种可重复使用冲压发动机的状态监视系统 | |
CN117782647B (zh) | 一种地铁行走部轴承健康状态监测系统 |
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 |