基于Kriging模型抗冰海洋平台管节点疲劳可靠性研究_宋晓杰

第25卷第1期2010年2月中国海洋平台CHINAOFFSHOREPLATFORMVol.25No.1Feb.,2010文章编号:1001-4500(2010)01-0032-04

基于Kriging模型抗冰海洋平台管节点疲劳可靠性研究

(1.大连理工大学,大连116011;2.大连理工大学船舶工程学院,大连116023)

摘 要:以渤海冰情和冰荷载现场观测数据为基础,建立JZ20-2NW海域锥体模型冰力谱。同时利用谱

分析方法,对该平台进行冰激疲劳分析。并结合冰疲劳环境荷载及冰力谱函数,提出了相对于冰速、冰厚随机

冰载的管节点热点应力幅的近似计算方法)))Kriging插值方法。结果表明,Kriging插值方法及疲劳可靠性

分析方法合理有效。

关键词:导管架平台;冰激振动;谱分析;疲劳可靠性;Kriging模型

中图分类号:O346.5 文献标识码:A宋晓杰1, 黄 一2, 张 崎2

StudyonFatigueReliabilityAnalysisofTubularJointsfor

Ice-ResistantOffshorePlatformBasedonKrigingModel

SONGXiao-jie, HUANGYi, ZHANGQi122

(1.DepartmentofContinuingEducation,DalianUniversityofTechnology,Dalian116011,China;

2.DepartmentofNavalArchitecture,DalianUniversityofTechnology,Dalian116023,China)

Abstract:Theice-forcespectrumofpyramidalmodelforJZ20-2NWplatforminitsseaar-

eaisestablishedbasedoniceconditionsintheBohaiandtheobservationalfielddataofice

load.Atthesametime,theice-inducedfatigueofplatformisanalyzedwithspectrumanalysis

method.Accordingtotheicefatigueloadenvironmentandtheice-forcespectrumfunction,

thispaperputsforwardanapproximatecalculationmethodofthetubularjoints'hotspotstress

rangerelativetoicevelocityandthicknessofstochasticiceload,whichcalledKriginginterpo-

lation.TheresultsshowthattheKriginginterpolationandtheanalysismethodoffatiguerel-i

abilityarereasonableandeffective.

Key words:jacketplatform;ice-inducedvibration;spectrumanalysis;fatiguereliabil-i

ty;Krigingmodel

0 引言

海洋平台是一种海洋工程结构物,它为开发和利用海洋资源提供了海上作业与生活的场所。

高纬度海域海冰对结构作用是结构设计和安全运营时需要考虑的重要因素。目前,渤海海洋平台在交变冰力作用下平台动力响应引起的管节点疲劳等问题越来越突出。有事实表明[2,3],只要有少数几个管节点发生疲劳破坏就可能导致整个平台的破坏。

本文基于疲劳可靠性分析方法,力求考虑更多的不确定影响因素,以对导管架平台在疲劳方面的安全性作出比确定性方法更加合理的评估。

收稿日期:2009-03-23; 修改稿收到日期:2009-09-30,女,,

第1期 宋晓杰等 基于Kriging模型抗冰海洋平台管节点疲劳可靠性研究#33#1 随机冰载作用下海洋平台管节点的热点应力计算

1.1 频域法确定热点应力幅值

频域法亦称谱疲劳分析方法[5,6],考虑荷载随机特性采用统计理论表述工况。假定冰载荷在短时间内是各态历经、窄带平稳的正态过程,根据应力循环转换的观点,结构的应力分布为窄带随机正态过程,其峰值服从瑞利分布。

岳前进等利用多年累积的冰力数据,对冰力进行谱分析工作,得到锥体冰力谱函数[1]:

S(f)=exp(-5.47T-0.64(1))ff

式中:F为冰力幅值,F=BRft()0.34;B为待定系数,实取3.7;Rf为冰的弯曲强度,实取700kPa;t为冰LC

0.25厚;D为冰作用处的锥体直径;LC为冰板的特征长度,LC=(),其中E为冰的弹性模量,实取0.512Qwg32-2.5

GPa;g为重力加速度;Qw为海水密度,实取1.0@10kg/m;k为破碎长厚比,在4~10之间取值;h为冰厚;V为冰速。

1.2 相对冰速、冰厚随机冰载的热点应力幅

由式(1)可看出,冰力谱函数与冰速、冰厚有关[7,8],而冰速、冰厚在冰期内是随机变化的。要进行随机荷载下的海洋平台疲劳可靠性分析,需要计算大量冰况组合下的管节点热点应力幅,但计算量很大。因此,本文提出基于Kriging方法近似计算热点应力幅[10,11]:(1)选择一系列有代表性的冰厚和冰速冰况组合,对海洋平台进行谱分析,计算出相对应的热点应力幅;(2)根据已知的样本点信息建立无偏的Kriging模型;

(3)利用Kiging模型插值计算整个疲劳环境的热点应力幅值。33

2 冰振下海洋平台管节点疲劳可靠性分析[7,8]

疲劳可靠性分析的最终目的是要从概率的角度对结构的疲劳寿命进行预测,从而对结构在疲劳方面的安全性作出评估。海洋平台管节点疲劳安全性的度量是计算所得的疲劳寿命Tf大于等于设计寿命TD的概率,即Pr=P(Tf\TD)。对于船舶及海洋工程结构,设计寿命通常取为20年或25年。

2.1 海洋平台疲劳寿命计算

疲劳寿命分析计算采用累积损伤模型,其以S-N曲线和Miner法则作为基础[9]。在船舶及海洋工程结构疲劳可靠性分析中应用最为广泛的是Wirsching模型[12,13]。

疲劳累积损伤: D=NTE(S)/A

式中NT为环循环次数总和。

设Sa为构件的应力幅,S为计算得到的应力幅,则有

Sa=BS

式中B是应力计算误差因子,属于随机变量。

定义应力循环的平均频率: f0=NT/T

则有 D=TB8/A

这里定义8为应力参数,并且

8=f0E(S)

不同的计算模型,对应着不同的8。

可得到疲劳寿命为 T=$A/(B8)(7)

2.2 结构疲劳可靠度计算的对数正态格式

工程实际中,常对基本随机变量或疲劳寿命本身的分布形式作一定的假设,这样可以使可靠度的计算得[14]mmmm[5](2)(3)(4)(5)(6)

#34#中国海洋平台 第25卷 第1期为广泛[15,16]。

由上(7)式知疲劳寿命:T=$A/(Bm8)

由于$、K、B是随机变量,T也是随机变量,疲劳失效概率为

Pf=P(T[TD)

应用对数正态分布法,即假定每一随机变量服从对数正态分布。根据对数正态分布的特性,则有:

Pf=<(-B)

式中B=ln(T/TD)/RlnT。

式中:RlnT为lnT的标准偏差;T为T的平均值;TD为设计寿命。T=$K/(Bm8)

lnT=ln[(1+C$)(1+CK)(1+CB)R2222m2(8)(9)(10)(11)](12)

式中Ci为变量的变异系数。

3 JZ20-2NW平台管节点疲劳可靠性分析

以渤海石油公司JZ20-2NW平台为例,在疲劳环境中取一个45b的冰力方向,在冰速0~100cm/s范围内,平均选取41种冰速;在冰厚0~45cm范围内,平均选取41种冰厚。S-N曲线采用API-X曲线。相关数据参考相关文献[5~9]。

选取有代表性的300种工况,利用ANSYS软件对海洋平台结构进行谱分析,最大节点应力均出现在两个斜撑的交叉处,计算出对应的热点应力幅值S。

根据上述建立的Kriging模型计算待测点的响应。在此过程中实现了MATLAB软件与ANSYS软件的互相调用,共插值得到1681种组合工况的结果。

程序运行结果如图1所示。

为了验证插值近似热点应力计算的精度,在冰速

为0~100cm/s、冰厚为0~45cm范围内,每间隔3.4

cm取一种冰厚,共13种冰厚;每间隔7.7cm/s取一

种冰速,共13种冰速,共有169种组合工况。

依据得到的结果统计得出:在给定的冰况条件下,

Kriging插值计算得到的近似结果与ANSYS谱分析

得到的计算结果基本吻合,二者误差小于1%,且随着

冰厚的增加,误差远低于1%。

李刚等提出的二次项拟合的近似计算热点应力幅

值的方法得到的误差为4%。

在运算时间方面 …… 此处隐藏：4280字，全部文档内容请下载后查看。喜欢就下载吧 ……