CN114564868B - 一种锚链疲劳寿命预测方法 - Google Patents

一种锚链疲劳寿命预测方法 Download PDF

Info

Publication number
CN114564868B
CN114564868B CN202210214483.XA CN202210214483A CN114564868B CN 114564868 B CN114564868 B CN 114564868B CN 202210214483 A CN202210214483 A CN 202210214483A CN 114564868 B CN114564868 B CN 114564868B
Authority
CN
China
Prior art keywords
crack
anchor chain
mooring
model
fatigue
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
Application number
CN202210214483.XA
Other languages
English (en)
Other versions
CN114564868A (zh
Inventor
何文涛
曹诗卉
谢凌俊
王昌梓
王树青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ocean University of China
Original Assignee
Ocean University of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ocean University of China filed Critical Ocean University of China
Priority to CN202210214483.XA priority Critical patent/CN114564868B/zh
Publication of CN114564868A publication Critical patent/CN114564868A/zh
Application granted granted Critical
Publication of CN114564868B publication Critical patent/CN114564868B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04—Ageing analysis or optimisation against ageing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Software Systems (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Mathematical Physics (AREA)
  • Computer Graphics (AREA)
  • Computational Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种锚链疲劳寿命预测方法，属于疲劳寿命的技术领域。该方法是通过疲劳裂纹扩展方法实现的，包括：前处理、迭代求解以及裂纹扩展，可分别借助SESAM，ABAQUS，FRANC 3D平台实现，具体步骤包括：S1模型建模、S2水动力分析、S3系泊分析；S4锚链有限元模型的建模、S5裂纹插入以及网格更新、S6增量计算判断、S7终止条件检验、S8疲劳寿命统计。本发明以裂纹深度以及裂纹长度双参数描述裂纹，利用椭圆公式准确计算更新后裂纹形状；实现了从系泊张力时程到裂纹扩展的发展，可有效表征真实海况下锚链的疲劳性能，并对其剩余寿命予以准确预测。

Description

一种锚链疲劳寿命预测方法
技术领域
本发明属于疲劳寿命的技术领域，具体涉及一种锚链疲劳寿命预测方法。
背景技术
对于海上浮式结构的系泊系统，结构物的往复运动、持续不断的波浪载荷将会导致锚链在服役期间始终处于交变载荷的作用下，并导致疲劳损伤的出现以及恶化。随着海上开发向深远海的发展、浮式结构的大量应用，锚链的疲劳破坏日益突出，已被大量的船级社、机构、学者公认为最重要的锚链失效模式之一。疲劳裂纹在锚链中的萌生以及扩展会明显减低锚链的疲劳强度以及承载能力，甚至会导致整根锚链的疲劳断裂失效并进而威胁整个海上结构物的安全。另外，海上浮式结构的锚链连续服役时间一般长达20年，其定期检查以及维护更换工作并不好进行，因此，基于断裂力学的方法对锚链的疲劳性能进行合理的评估，利用裂纹扩展的方法对锚链的剩余疲劳寿命进行准确预测是很有必要的。
针对锚链疲劳问题，S-N/T-N曲线法由于其便捷性、工程上可接受的精度，在锚链的疲劳寿命评估中得到了广泛的认可以及应用。但是，S-N/T-N曲线法理论本身的缺陷也限制了该方法的更深入的应用，比如：
(1)S-N/T-N曲线法无法考虑锚链在制造、运输、安装时引入的初始缺陷，严重高估锚链疲劳性能；
(2)S-N/T-N无法考虑载荷次序、材料特性等问题；
(3)无法基于裂纹扩展路径以及裂纹形貌对锚链寿命进行准确预测；
由于S-N/T-N曲线法的诸多不足，断裂力学在锚链疲劳性能、寿命预测中的应用也越来越迫切，目前，限制断裂力学在锚链疲劳问题中应用的主要因素如下：
(1)锚链服役期内遭受的实际海况极其复杂且极不规则，将实际海况进行简化并与断裂力学进行结合存在着较大问题；
(2)基于断裂力学开展疲劳寿命预测，将疲劳问题视为裂纹扩展问题，即裂纹与有限元模型的更新与迭代计算。锚链在服役期内能够遭受超过百万量级的循环载荷，如果直接将如此大量的循环载荷应用于迭代计算，工作量巨大且计算效率低下，不利于基于断裂力学实现锚链的疲劳寿命快速预测；
发明内容
为将断裂力学成功应用到锚链的疲劳断裂问题中，本发明提出了一种锚链疲劳寿命预测方法。本发明利用线性插值、椭圆方程等成功地将实际海况与裂纹扩展的精确计算结合了起来。推进了断裂力学在锚链疲劳性能研究、剩余寿命预测中应用。
本发明是由以下计算方案实现的：
一种锚链疲劳寿命预测方法，该方法是通过疲劳裂纹扩展方法实现的，包括：前处理、迭代求解以及裂纹扩展，可分别借助SESAM，ABAQUS，FRANC 3D平台实现，其中，SESAM平台作为前处理器，其工况包括：S1模型建模、S2水动力分析、S3系泊分析；ABAQUS负责S4锚链有限元模型的建模，具体包括以下几个步骤：S41简化锚链模型的建模、S42裂纹扩展区域的定义、S43边界载荷以及系泊张力的施加、S44有限元分析；FRANC 3D负责S5裂纹插入以及网格更新，整个过程包括初始裂纹的插入以及裂纹自动扩展，其中，裂纹自动扩展的实现包括以下主要环节：S6增量计算判断、S7终止条件检验、S8疲劳寿命统计；
所述S1模型包括S11几何模型、S12水动力模型，其中，S11几何模型包含结构物尺寸及结构信息；S12水动力模型是通过对几何模型的水动力性质的定义和赋值来实现，是赋予几何模型水动力属性后完成，包含湿表面的定义以及压载舱水动力属性的定义；
所述S2水动力分析为求解结构自身运动特性和波浪载荷，其工作的开展是基于势流理论完成，内容包括：S21波浪环境、S22水动力属性赋予以及S23浮态调整步骤；其中，S21波浪环境是指定义势流理论分析用一系列规则波，包括浪向、波频以及工作水深等信息定义；S22水动力属性赋予是指莫里森结构阻尼系数等水动力相关参数的定义；S23浮态调整则是针对重浮力不相等或首尾吃水与实际存在差异的情况进行的模型姿态调整，通过压载水的方式进行处理；
所述S3系泊分析为求解锚链系泊张力时程，其主要工作包括：S31系泊系统定义、S32环境条件定义、S33系泊张力幅值与循环次数提取；其中，S31系泊系统定义指系泊缆方位定义、长度定义、材料属性及水动力属性定义；S32环境条件定义指根据散布图进行定义，包括波高、周期、浪向、风速、风向、流速流向等；S33系泊张力幅值与循环次数提取指利用雨流计数法对锚链系泊张力时程进行张力幅值序列的提取；
所述S41简化锚链模型的建模：建模工作以规范中规定的标准无档/有档锚链为参考，取至少一个完整链环与两段半环建模；
所述S42裂纹扩展区域的定义：根据疲劳热点所在位置或者常见疲劳裂纹源对模型进行分割，进行网格划分，并在裂纹扩展区域划分较为精细的网格、而其他区域保持相对较粗的网格以平衡计算精度以及计算效率；
所述S43边界载荷以及系泊张力的施加：为准确模拟链环之间的相互作用，链环接触面之间施加面面接触，在整段锚链的一端加对称边界条件，而另一端耦合到参考点并施加系泊张力；
所述S44有限元分析：对插入裂纹的锚链模型进行有限元分析以求得求解裂纹前缘应力强度因子所需的应力、位移等；
所述S5裂纹插入及网格更新：将初始裂纹以及更新后裂纹插入到适当的位置，在FRANC 3D中完成裂纹的插入以及裂纹扩展区网格的更新；
所述S6增量计算判断，是指对当前海况下的裂纹扩展增量进行判断，以确定其值是否需要记入裂纹累积扩展量；
所述S7终止条件校验，是指对程序的终止条件进行判断；
所述S8服役年限统计根据波浪散布图之遍历次数确定，散布图遍历一遍则计数一年，总共服役年限视为裂纹尺寸扩展至临界值时散布图遍历次数。
进一步的，该方法还包括自适应裂纹增量优化控制，疲劳裂纹扩展过程中裂纹前缘应力强度因子的计算需要多次插入裂纹并更新有限元模型，因此扩展增量阈值采用自适应步长设置，仿真过程中，将其视为一个与裂纹大小相关的参数，只有当裂纹扩展增量超过阈值时，裂纹的插入和有限元模型的更新才会被执行。如果对每一个系泊张力作用下的裂纹尺寸都予以更新无疑可以有效的提高仿真精度，但频繁的裂纹插入和模型更新势必会导致有限元模型的反复计算，这无疑会大大降低断裂力学方法在实际工程结构疲劳寿命预测中的可行性，因此需要对裂纹扩展增量的阈值予以合理控制。考虑到疲劳裂纹扩展是一个逐渐加速扩展的过程，固定的阈值通常不能适用于整个裂纹扩展的过程。在裂纹萌生阶段阈值可能过大，导致预测精度降低，而在裂纹快速扩展阶段固定阈值又可能过小，不利于计算效率的提高。
进一步的，所述S6增量计算判断包括如下步骤：
S61、根据应力强度因子、线性插值、线性累加计算当前海况下的扩展增量da,dL；
S62、对比当前裂纹扩展量与预设阈值th_Δa，如果扩展增量小于阈值，则循环将自动进入下一个海况的计算，如果裂纹扩展增量大于预设阈值，则当前扩展增量将累积至上一个循环结果，并更新裂纹尺寸和对应的有限元模型；
S63、利用椭圆公式计算更新后裂纹尺寸形状。
进一步的，所述S61包括如下步骤：
S611、应力强度因子确定：对于施加的系泊张力幅值T0，可得到等效应力强度因子
Figure BDA0003533829350000041
其中，KI,KII和KIII分别代表三种典型的裂纹的应力强度因子，即张开型(Mode I)、滑开型(Mode II)和撕开型(Mode III)，v是材料的泊松比；
S612、线性插值：对于任一海况对应的系泊张力时程，经雨流计数均可得到一系列张力幅值，因此，利用线性插值可以得到该海况对应的一系列张力幅值对应的应力强度因子
Figure BDA0003533829350000042
其中Keqi为当前海况张力幅值Ti所对应的等效应力强度因子，Keq0为所施加张力幅值T0所对应的等效应力强度因子；
S613、裂纹增量确定：运用Paris公式，
Figure BDA0003533829350000043
可以求得裂纹上任一点扩展量，其中
Figure BDA0003533829350000051
为裂纹前沿上某点的扩展速率，ΔKeqs为该点所对应的等效应力强度因子范围，C,m分别为材料常数；
对于本发明，主要关注裂纹前沿最深点以及裂纹前沿表面点，因此利用上述S612得到的张力幅值序列对应的等效应力强度因子分别代入S613的Paris公式便可得到每一个张力幅值所对应的裂纹扩展量，然后进行线性累加便可得到该海况对应的裂纹增量dai,dLi；
上述S612所得张力幅值序列对应的等效应力强度因子
对于上述S62对比当前裂纹扩展量与预设阈值th_Δa，通过上述S61计算的裂纹扩展量与th_Δa对比，当裂纹扩展量小于阈值时，循环将自动进入下一个海况重复S612、S613并累加裂纹扩展量，如果裂纹扩展量大于阈值，则当前裂纹扩展量累加至上一循环结果。
进一步的，所述S63包括如下步骤：
裂纹尺寸ai+1,ci+1计算：为准确计算裂纹尺寸，本步骤引入椭圆公式计算出更新后裂纹尺寸，
Figure BDA0003533829350000052
其中，ai+1,ci+1分别为更新后裂纹深度以及裂纹半长，xi+1,yi+1分别为更新后裂纹前沿与链环表面交点的坐标，其中，ai+1通过S61以及S62计算可直接得到，而ci+1为椭圆公式求解对象，另外，xi+1,yi+1可分别表达如下：
xi+1＝xi+dL·cosθ
yi+1＝yi+dL·sinθ
其中dL为当前裂纹前沿在其与链环表面交点处的裂纹扩展量，θ为dL对应的平面内扩展角，cosθ,sinθ可根据下式求解，
Figure BDA0003533829350000053
Figure BDA0003533829350000054
Figure BDA0003533829350000061
Figure BDA0003533829350000062
其中
Figure BDA0003533829350000063
为裂纹面的法向量，
Figure BDA0003533829350000064
为裂纹面内由裂纹前沿表面点B指向裂纹前沿表面点A的向量，
Figure BDA0003533829350000065
和
Figure BDA0003533829350000066
分别为更新前裂纹前沿表面点的切向量以及法向量，且点C与点A同侧。
进一步的，所述S7终止条件校验，包括如下步骤：
S71、如果更新后的裂纹尺寸达到临界值，终止计算并跳出循环；
S72、如果当前裂纹尺寸尚未触及临界条件，则循环继续，并进入下一个海况的循环。
有益效果
(1)本发明提出的方法可有效表征真实海况下锚链的疲劳性能，并对其剩余寿命予以准确预测；
(2)通过线性插值、线性累加的应用，将每个海况对应的总扩展量转化为对应的一系列张力幅值分别对应的裂纹扩展量的线性叠加。此法既保证了裂纹扩展量的精确计算，也保证了计算效率。
(3)通过椭圆公式的应用，精确计算了更新后的裂纹形。将裂纹简化为由裂纹深度以及裂纹长度双参数描述的裂纹，通过深度以及长度的计算来描述更新后的裂纹，既能保证计算精度，也能保证计算效率。
(4)通过二次开发所建立的锚链疲劳寿命评估平台，实现了学科间协同作业以及数据共享和信息传递，推进了断裂力学在锚链疲劳性能研究、剩余寿命预测的应用。
本发明为解决断裂力学在锚链疲劳问题研究中的应用，利用线性插值、线性叠加将每个海况的裂纹总扩展量转化为对应的一系列张力幅值的裂纹扩展量的叠加；同时以裂纹深度以及裂纹长度双参数描述裂纹，利用椭圆公式准确计算更新后裂纹形状；通过对现有商业软件的二次开发，建立跨学科耦合分析平台，实现了从系泊张力时程到裂纹扩展的发展。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本发明的疲劳裂纹扩展寿命评估流程图；
图2为疲劳裂纹自动扩展实现方法；
图3为验证模型系泊方案，其中(a)为系泊系统立体示意图，(b)为平面示意图；
图4为系泊张力时程曲线与雨流计数结果，其中，(a)为海况Hs＝1.5m、Ts＝6.5s系泊张力时程对应的结果图，(b)为海况Hs＝1.5m、Ts＝6.5s张力幅值序列对应的结果图，(c)为海况Hs＝6.5m、Ts＝8.5s系泊张力时程对应的结果图，(d)为海况Hs＝6.5m、Ts＝8.5s张力幅值序列对应的结果图，(e)为海况Hs＝10.5m、Ts＝11.5s系泊张力时程对应的结果图，(f)为海况Hs＝10.5m、Ts＝11.5张力幅值序列对应的结果图；
图5为锚链有限元模型，其中(a)为无档锚链尺寸标准，(b)为锚链的有限元模型、边界条件、相互作用，(c)为常见裂纹起源位置所对应的裂纹扩展区划分；
图6为裂纹插入及网格划分；
图7为应力强度因子在裂纹扩展过程中的变化，其中(a)为等效应力强度因子，(b)为I型应力强度因子KI，(c)为II型应力强度因子KII，(d)为III型应力强度因子KIII；
图8为线性插值线性累加示意图；
图9裂纹尺寸计算示意图；
图10为典型的裂纹扩展形貌，其中(a)为锚链裂纹扩展应力图，(b)为锚链链环中的裂纹形貌，(c)为裂纹扩展过程中裂纹每年的形貌；
图11为疲劳寿命示意图，其中(a)为不同锚链的疲劳寿命曲线，(b)为锚链上位置S1、位置S2、位置S3、位置S4在Crown截面上的疲劳寿命曲线(c)为锚链上位置S1、位置S2、位置S3、位置S4在Crown截面上的寿命预测，(d)为链环上Crown截面、Bend截面、Straight截面的疲劳寿命曲线，(e)为链环上Crown截面、Bend截面、Straight截面的疲劳寿命预测值。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
为验证本发明所提出的基于断裂力学的锚链疲劳寿命预测的方法，本发明以浮式生产储油船(FPSO)的系泊锚链的疲劳寿命预测为例做进一步阐述。
如图1～11所示，本实施例具体包括如下步骤：
S1：模型建模
本实施例模型包括S11几何模型以及S12水动力模型，其中，
S11：几何模型的创建应包含结构物尺寸及结构信息，根据结构特点，选择不同建模方式。
S12：水动力模型是通过对几何模型的水动力性质的定义和赋值来实现，是赋予几何模型水动力属性后完成，包含湿表面的定义以及压载舱水动力属性的定义，其结构可只包含湿表面区域，而忽略其内部结构细节，但应包含舱室信息，以保证外海压力和舱室内部压力均可得以施加。
S2：水动力分析
水动力分析为求解结构自身运动特性和波浪载荷，其工作的开展是基于势流理论完成，内容包括：S21波浪环境定义、S22水动力属性赋予以及S23浮态调整步骤；
S21：波浪环境是指定义势流理论分析用一系列规则波，包括浪向、波频以及工作水深等信息定义；
S22：水动力属性赋予是指类似莫里森结构阻尼系数等水动力相关参数的定义；本实施例水动力分析计算主要基于三维势流理论完成；
S23：浮态调整则是针对重浮力不相等或首尾吃水与实际存在差异的情况进行的模型姿态调整，多通过压载水的方式进行处理；
S3：系泊分析
系泊分析为求解锚链系泊张力时程，内容包括：S31系泊系统定义、S32环境条件定义、S33系泊张力幅值与循环次数提取；
S31系泊系统定义指系泊缆方位定义、长度定义、材料属性及水动力属性定义，如图3所示；
S32环境条件定义指根据散布图进行定义，包括波高、周期、浪向、风速、风向、流速流向等；
S33系泊张力幅值与循环次数提取指利用雨流计数法对锚链系泊张力时程进行对应的张力幅值序列的提取。系泊张力时程与雨流计数结果如图4所示。
S4：锚链有限元模型的建模
锚链有限元模型的建模，具体包括以下几个步骤：S41简化锚链模型的建模、S42裂纹扩展区域的定义、S43边界载荷以及系泊张力的施加、S44有限元分析。模型如图5所示。
S41简化锚链模型的建模：建模工作以规范中规定的标准无档/有档锚链为参考，取至少一个完整链环与两段半环建模；
S42裂纹扩展区域的定义：根据疲劳热点所在位置或者常见疲劳裂纹源对模型进行分割，进行网格划分，并在裂纹扩展区域划分较为精细的网格、而剩余区域保持相对较粗的网格以平衡计算精度以及计算效率；
S43边界载荷以及系泊张力的施加：为准确模拟链环之间的相互作用，链环接触面之间施加面面接触，在整段锚链的一端施加对称边界条件，而另一端耦合到参考点并施加系泊张力；
S44有限元分析：对插入裂纹的锚链模型进行有限元分析以求得求解裂纹前缘应力强度因子所需的应力、位移等；
S5：裂纹插入及网格更新：
将初始裂纹以及更新后裂纹插入到适当的位置，在FRANC 3D中完成裂纹的插入以及裂纹扩展区网格的更新；裂纹插入及网格更新流程如图6所示；
本实施例以初始裂纹a0＝0.5mm,c0＝0.5mm进行研究，根据插入的裂纹(包括初始裂纹以及更新后的裂纹)重新划分S42所定义裂纹扩展区域网格；
裂纹自动扩展方法如图2所示，具体包括步骤S6、S7、S8；详述如下：
S6：裂纹增量计算
计算模型在初始裂纹和更新后裂纹前缘处的应力强度因子值，并通过线性插值、线性累加计算当前海况下的裂纹增量，通过当前裂纹扩展量与预设阈值th_Δa的对比，判断下一步是更新裂纹抑或进入下一个海况的循环；若更新裂纹，则通过椭圆公式进行更新后裂纹形状的计算。
S61、如图8所示，根据应力强度因子、线性插值、线性累加计算当前海况下的扩展增量da,dL，具体步骤如下：
S611、应力强度因子确定，以张力幅值T0计算裂纹前沿应力强度因子并求得裂纹深处以及裂纹表面点的等效应力强度因子
Figure BDA0003533829350000101
Figure BDA0003533829350000102
其中，Keqa,KeqL分别为裂纹深处以及裂纹表面点的等效应力强度因子，v为材料泊松比。
S612、通过线性插值、线性累加计算当前海况裂纹扩展量，利用等效应力强度因子插值公式
Figure BDA0003533829350000103
Figure BDA0003533829350000104
Ti为张力时程对应的一系列张力幅值的遍历，Keqai,KeqLi为Ti对应的裂纹深度、裂纹表面点的等效应力强度因子。
S613、裂纹增量确定：运用Paris公式，本实施例针对裂纹最深处、表面点运用如下公式：
Figure BDA0003533829350000105
Figure BDA0003533829350000106
可得到所有张力幅值所对应的裂纹增量，而后，进行线性叠加即可得到该海况所对应的裂纹扩展量。
Figure BDA0003533829350000111
Figure BDA0003533829350000112
Figure BDA0003533829350000113
S62、对比当前裂纹扩展量da与预设阈值th_Δa，如果扩展增量小于阈值，则循环将自动进入下一个海况的计算，如果裂纹扩展增量大于预设阈值，则当前扩展增量将累积至上一个循环结果，并更新裂纹尺寸和对应的有限元模型；在本实施例中，th_Δa设置为10％当前裂纹的深度即10％ai。
S63、裂纹尺寸ai+1,ci+1计算：如图9所示，为准确计算裂纹尺寸，本发明引入椭圆公式计算出更新后裂纹尺寸。
Figure BDA0003533829350000114
其中，ai+1,ci+1分别为更新后裂纹深度以及裂纹半长，xi+1,yi+1分别为更新后裂纹前沿与链环表面交点的坐标，ai+1通过S61以及S62计算可直接得到；而ci+1为椭圆公式求解对象，另外，xi+1,yi+1可分别表达如下：
xi+1＝xi+dL·cosθ
yi+1＝yi+dL·sinθ
其中dL为当前裂纹前沿在其与链环表面交点处的裂纹扩展量，θ为dL对应的平面内扩展角，cosθ,sinθ可根据下式求解，
Figure BDA0003533829350000115
Figure BDA0003533829350000116
Figure BDA0003533829350000121
Figure BDA0003533829350000122
其中
Figure BDA0003533829350000123
为裂纹面的法向量，
Figure BDA0003533829350000124
为裂纹面内由裂纹前沿表面点B指向裂纹前沿表面点A的向量，
Figure BDA0003533829350000125
和
Figure BDA0003533829350000126
分别为更新前裂纹前沿表面点的切向量以及法向量，且点C与点A同侧。典型的裂纹扩展形貌如图10所示。
S7终止条件校验，包括如下步骤:
S71、如果更新后的裂纹尺寸达到临界值，终止计算并跳出循环；
S72、如果当前裂纹尺寸尚未触及临界条件，则循环继续，并进入下一个海况的循环。
S8服役年限统计根据散布图的遍历次数确定，散布图遍历一遍计数一年，依据裂纹尺寸扩展至临界值时的散布图遍历计数便可统计出服役年限。
如图11所示，(a)为本实施例所选锚链2、5、8、11位置为S3的Crown截面计算所得的疲劳寿命曲线；
(b)、(c)为本实施例分别计算了锚链位置S1、位置S2、位置S3、位置S4在Crown截面的寿命预测，认为当裂纹深度a达到15％锚链直径时锚链失效；
(d)、(e)为本实施例分别计算了在Crown截面、Bend截面、Straight截面进行的寿命预测，认为裂纹深度a分别达到了15％、30％、12％的锚链直径时锚链失效。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (6)

1.一种锚链疲劳寿命预测方法，该方法是通过疲劳裂纹扩展方法实现的，其特征在于，该方法包括：前处理、迭代求解以及裂纹扩展，可分别借助SESAM，ABAQUS，FRANC 3D平台实现；其中，SESAM平台作为前处理器，其工况包括：S1模型建模、S2水动力分析、S3系泊分析；ABAQUS负责S4锚链有限元模型的建模，具体包括以下几个步骤：S41简化锚链模型的建模、S42裂纹扩展区域的定义、S43边界载荷以及系泊张力的施加、S44有限元分析；FRANC 3D负责S5裂纹插入以及网格更新，整个过程包括初始裂纹的插入以及裂纹自动扩展，其中，裂纹自动扩展的实现包括以下主要环节：S6增量计算判断、S7终止条件检验、S8疲劳寿命统计；
所述S1模型包括S11几何模型、S12水动力模型，其中，S11几何模型包含结构物尺寸及结构信息；S12水动力模型是通过对几何模型的水动力性质的定义和赋值来实现，是赋予几何模型水动力属性后完成，包含湿表面的定义以及压载舱水动力属性的定义；
所述S2水动力分析为求解结构自身运动特性和波浪载荷，其工作的开展是基于势流理论完成，内容包括：S21波浪环境、S22水动力属性赋予以及S23浮态调整步骤；其中，S21波浪环境是指定义势流理论分析用一系列规则波，包括浪向、波频以及工作水深等信息定义；S22水动力属性赋予是指莫里森结构阻尼系数等水动力相关参数的定义；S23浮态调整则是针对重浮力不相等或首尾吃水与实际存在差异的情况进行的模型姿态调整，通过压载水的方式进行处理；
所述S3系泊分析为求解锚链系泊张力时程，其主要工作包括：S31系泊系统定义、S32环境条件定义、S33系泊张力幅值与循环次数提取；其中，S31系泊系统定义指系泊缆方位定义、长度定义、材料属性及水动力属性定义；S32环境条件定义指根据散布图进行定义，包括波高、周期、浪向、风速、风向、流速流向等；S33系泊张力幅值与循环次数提取指利用雨流计数法对锚链系泊张力时程进行张力幅值序列的提取；
所述S41简化锚链模型的建模：建模工作以规范中规定的标准无档/有档锚链为参考，取至少一个完整链环与两段半环建模；
所述S42裂纹扩展区域的定义：根据疲劳热点所在位置或者常见疲劳裂纹源对模型进行分割，进行网格划分，并在裂纹扩展区域划分较为精细的网格、而其他区域保持相对较粗的网格以平衡计算精度以及计算效率；
所述S43边界载荷以及系泊张力的施加：为准确模拟链环之间的相互作用，链环接触面之间施加面面接触，在整段锚链的一端施加对称边界条件，而另一端耦合到参考点并施加系泊张力；
所述S44有限元分析：对插入裂纹的锚链模型进行有限元分析以求得求解裂纹前缘应力强度因子所需的应力、位移等；
所述S5裂纹插入及网格更新：将初始裂纹以及更新后裂纹插入到适当的位置，在FRANC3D中完成裂纹的插入以及裂纹扩展区网格的更新；
所述S6增量计算判断，是指对当前海况下的裂纹扩展增量进行判断，以确定其值是否需要记入裂纹累积扩展量；
所述S7终止条件校验，是指对程序的终止条件进行判断；
所述S8服役年限统计根据波浪散布图之遍历次数确定，散布图遍历一遍则计数一年，总共服役年限视为裂纹尺寸扩展至临界值时散布图遍历次数。
2.根据权利要求1所述疲劳寿命预测方法，其特征在于，该方法还包括自适应裂纹增量优化控制，疲劳裂纹扩展过程中裂纹前缘应力强度因子的计算需要多次插入裂纹并更新有限元模型，因此扩展增量阈值采用自适应步长设置，仿真过程中，将其视为一个与裂纹大小相关的参数，只有当裂纹扩展增量超过阈值时，裂纹的插入和有限元模型的更新才会被执行。
3.根据权利要求1或2任一项所述的疲劳寿命预测方法，其特征在于，所述S6增量计算判断包括如下步骤：
S61、根据应力强度因子、线性插值、线性累加计算当前海况下的扩展增量da,dL；
S62、对比当前裂纹扩展量与预设阈值th_△a，如果扩展增量小于阈值，则循环将自动进入下一个海况的计算，如果裂纹扩展增量大于预设阈值，则当前扩展增量将累积至上一个循环结果，并更新裂纹尺寸和对应的有限元模型；
S63、利用椭圆公式计算更新后裂纹尺寸形状。
4.根据权利要求3所述的疲劳寿命预测方法，其特征在于，所述S61包括如下步骤：
S611、应力强度因子确定：对于施加的系泊张力幅值T0，可得到等效应力强度因子
Figure FDA0003533829340000021
其中，KI,KII和KIII分别代表三种典型的裂纹的应力强度因子，即张开型(Mode I)、滑开型(Mode II)和撕开型(Mode III)，v是材料的泊松比；
S612、线性插值：对于任一海况对应的系泊张力时程，经雨流计数均可得到一系列张力幅值，因此，利用线性插值可以得到该海况对应的一系列张力幅值对应的应力强度因子
Figure FDA0003533829340000022
其中Keqi为当前海况张力幅值Ti所对应的等效应力强度因子，Keq0为所施加张力幅值T0所对应的等效应力强度因子；
S613、裂纹增量确定：运用Paris公式，
Figure FDA0003533829340000031
可以求得裂纹上任一点扩展量，其中
Figure FDA0003533829340000032
为裂纹前沿上某点的扩展速率，△Keqs为该点所对应的等效应力强度因子范围，C,m分别为材料常数。
5.根据权利要求3所述的疲劳寿命预测方法，其特征在于，所述S63包括如下步骤：
裂纹尺寸ai+1,ci+1计算：为准确计算裂纹尺寸，本步骤引入椭圆公式计算出更新后裂纹尺寸，
Figure FDA0003533829340000033
其中，ai+1,ci+1分别为更新后裂纹深度以及裂纹半长，xi+1,yi+1分别为更新后裂纹前沿与链环表面交点的坐标，ai+1通过S61以及S62计算可直接得到；而ci+1为椭圆公式求解对象，另外，xi+1,yi+1可分别表达如下：
xi+1＝xi+dL·cosθ
yi+1＝yi+dL·sinθ
其中dL为更新后裂纹前沿在其与链环表面交点处的裂纹扩展量，θ为dL对应的平面内扩展角，cosθ,sinθ可根据下式求解，
Figure FDA0003533829340000034
Figure FDA0003533829340000035
Figure FDA0003533829340000036
Figure FDA0003533829340000037
其中
Figure FDA0003533829340000038
为裂纹面的法向量，
Figure FDA0003533829340000039
为裂纹面内由裂纹前沿表面点B指向裂纹前沿表面点A的向量，
Figure FDA00035338293400000310
和
Figure FDA00035338293400000311
分别为更新前裂纹前沿表面点的切向量以及法向量，且点C与点A同侧。
6.根据权利要求1或2任一项所述的疲劳寿命预测方法，其特征在于，所述S7终止条件校验，包括如下步骤：
S71、如果更新后的裂纹尺寸达到临界值，终止计算并跳出循环；
S72、如果当前裂纹尺寸尚未触及临界条件，则循环继续，并进入下一个海况的循环。
CN202210214483.XA 2022-03-07 2022-03-07 一种锚链疲劳寿命预测方法 Active CN114564868B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210214483.XA CN114564868B (zh) 2022-03-07 2022-03-07 一种锚链疲劳寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210214483.XA CN114564868B (zh) 2022-03-07 2022-03-07 一种锚链疲劳寿命预测方法

Publications (2)

Publication Number Publication Date
CN114564868A CN114564868A (zh) 2022-05-31
CN114564868B true CN114564868B (zh) 2023-05-12

Family

ID=81718071

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210214483.XA Active CN114564868B (zh) 2022-03-07 2022-03-07 一种锚链疲劳寿命预测方法

Country Status (1)

Country Link
CN (1) CN114564868B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115238583B (zh) * 2022-07-27 2024-02-13 山东理工大学 一种支持增量日志的业务流程剩余时间预测方法与系统
CN116840356B (zh) * 2023-09-01 2023-11-17 南京安盛电子有限公司 一种灌封变压器裂纹的监测方法
CN118114384B (zh) * 2024-04-30 2024-07-12 北京航空航天大学 保证裂纹前缘形状一致的涡轮盘裂纹扩展模拟件设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002047313A2 (en) * 2000-10-26 2002-06-13 Vextec Corporation Method and apparatus for predicting the failure of a component
JP2006010378A (ja) * 2004-06-23 2006-01-12 Bridgestone Corp タイヤ経時変化予測方法、装置、プログラム及び媒体
CN105005695A (zh) * 2015-07-14 2015-10-28 中国海洋大学 一种用于时域疲劳分析的波浪散布图组块等效方法
WO2020042753A1 (zh) * 2018-08-30 2020-03-05 长沙理工大学 季节性腐蚀与疲劳耦合作用下钢筋混凝土桥梁寿命预测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210383035A1 (en) * 2018-11-01 2021-12-09 Siemens Energy Global GmbH & Co. KG Computer-implemented method for the probabilistic estimation of a probability of failure of a component, a data processing system, a computer program product and a computer-readable storage medium

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002047313A2 (en) * 2000-10-26 2002-06-13 Vextec Corporation Method and apparatus for predicting the failure of a component
JP2006010378A (ja) * 2004-06-23 2006-01-12 Bridgestone Corp タイヤ経時変化予測方法、装置、プログラム及び媒体
CN105005695A (zh) * 2015-07-14 2015-10-28 中国海洋大学 一种用于时域疲劳分析的波浪散布图组块等效方法
WO2020042753A1 (zh) * 2018-08-30 2020-03-05 长沙理工大学 季节性腐蚀与疲劳耦合作用下钢筋混凝土桥梁寿命预测方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
Effect of single tensile overload on fatigue crack growth behavior based on plastically dissipated energy and critical distance theory;he wentao and etc,.;Engineering Fracture Mechanics;全文 *
Numerical study on fatigue crack growth at a web-stiffener of ship structural details by an objected-oriented approach in conjunction with ABAQUS;Wentao He and etc,.;Marine Structures;全文 *
Probabilistic life assessment on fatigue crack growth in mixed-mode by coupling of Kriging model and finite element analysis;he wentao and etc,.;Engineering Fracture Mechanics;全文 *
一种新型半潜式平台概念设计研究;王天英;王树青;;海洋工程(05);全文 *
二阶差频力对半潜式平台系泊锚链疲劳损伤影响;宋宪仓，王树青，杜君峰;第十七届中国海洋（岸）工程学术讨论会论文集;全文 *
单点系泊系统将军柱疲劳分析;孟轲;董海防;朱刚;;舰船科学技术(05);全文 *
深水半潜式钻井平台节点疲劳寿命谱分析研究;马网扣;王志青;张海彬;;海洋工程(03);全文 *
疲劳谱分析方法及其在立柱式平台结构分析中的应用;李华祥;;海洋工程装备与技术(06);全文 *

Also Published As

Publication number Publication date
CN114564868A (zh) 2022-05-31

Similar Documents

Publication Publication Date Title
CN114564868B (zh) 一种锚链疲劳寿命预测方法
Shittu et al. Comparative study between SN and fracture mechanics approach on reliability assessment of offshore wind turbine jacket foundations
CN114218661B (zh) 一种基于疲劳裂纹扩展的疲劳寿命预测方法
Park et al. Numerical investigation and development of design formula for cylindrically curved plates on ships and offshore structures
CN116205109A (zh) 多裂纹耦合扩展轨迹计算方法、设备及存储介质
CN105005695A (zh) 一种用于时域疲劳分析的波浪散布图组块等效方法
CN107180123B (zh) 一种高强度钢潜水器耐压球壳极限承载力估算方法
CN115860197A (zh) 一种基于数据驱动的煤层气产量预测方法及系统
CN112836401B (zh) 提高复杂环境条件下系泊系统安全性的分析方法
CN107944100B (zh) 综合生产脐带缆的线型设计方法
CN116579266B (zh) 一种波动水位影响下的河岸边坡崩塌预警方法及系统
CN113609555A (zh) 一种基于大数据技术的水工金属结构设计方法
Lacasta et al. Simulation of PID control applied to irrigation channels
KR101627771B1 (ko) 석션 케이슨 앵커의 인발력을 예측하기 위한 데이터 지향 모델링 방법
CN107563013A (zh) 三维编织复合材料的车辆配件及其制作方法
CN108563846B (zh) 一种基于可靠性的海洋立管波致疲劳安全系数的确定方法
Vanem et al. Environmental contours based on a direct sampling approach and the IFORM approach: Contribution to a benchmark study
CN113836696A (zh) 一种持续锈蚀状态下闸门极限承载力校核方法
Birhane et al. A Computational Framework for the Aerodynamic Shape Optimization of Long-Span Bridge Decks
Dermatis et al. Stochastic Prediction of the Extreme Responses of a Floating Offshore Wind Turbine Using Response-Conditioned Waves
Muys et al. Cycle-by-cycle simulation of variable amplitude fatigue crack propagation
CN118395739B (zh) 一种浮式平台运动响应计算方法、系统和存储介质
Crouch The effect of crack front curvature and side-grooving on three point bend specimen fracture toughness measurements
Zhang et al. Stochastic damage modelling of mixed mode fatigue delamination of composite wind turbine blades
CN118296933A (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