CN103714204A - 焊接结构多轴疲劳寿命评估方法 - Google Patents
焊接结构多轴疲劳寿命评估方法 Download PDFInfo
- Publication number
- CN103714204A CN103714204A CN201310712305.0A CN201310712305A CN103714204A CN 103714204 A CN103714204 A CN 103714204A CN 201310712305 A CN201310712305 A CN 201310712305A CN 103714204 A CN103714204 A CN 103714204A
- Authority
- CN
- China
- Prior art keywords
- stress
- sigma
- critical point
- fatigue life
- cpsp
- 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.)
- Granted
Links
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种焊接结构多轴疲劳寿命评估方法,具有如下步骤:建立待测定焊接结构的有限元模型;在焊接结构处施加多种实验载荷,求解得到焊趾端面沿板厚方向所有节点的结构应力张量[σ];根据节点的结构应力张量[σ]计算节点对应的载荷比例系数R,进而确定应力峰值临界点的位置;使用插值的方法得到所述应力峰值临界点的结构应力张量[σCPSP];根据应力峰值临界点的结构应力张量[σCPSP]计算得出疲劳寿命临界损伤面的最大剪应力和正应力范围,应用联合多轴疲劳寿命预测方法MWCM进行多轴疲劳寿命预测。
Description
技术领域
本发明涉及一种焊接疲劳结构分析方法,尤其涉及一种焊接结构多轴疲劳寿命的评估方法。
背景技术
海洋平台焊接结构处于多轴疲劳损伤状态。一是焊接结构本身受到不同类型载荷的联合作用;二是焊接部位复杂的几何外形,导致局部区域处于多向应力状态。
目前,规范上主要基于单轴疲劳方法来评估海洋平台的疲劳寿命,是导致疲劳断裂事故仍未得到完全控制的主要原因之一。国际上,虽然提出了多种多轴疲劳寿命评估方法,但是到目前为止,还没有方法可以应用于实际复杂的海洋平台焊接结构。
热点应力法已成为工程上最受欢迎的焊接结构疲劳寿命评估方法。从工程应用角度来看,国际上最新的成果是基于热点应力法并结合修正沃勒曲线法MWCM评估焊接接头的多轴疲劳寿命,但是该方法只能处理载荷类型和结构形式都非常简单的焊接接头,其主要原因是利用结构表面应力插值获得的应力难以处理实际海洋平台焊接结构焊趾处复杂且不断变化的应力状态。
发明内容
本发明针对以上问题的提出,而研制的一种焊接结构多轴疲劳寿命评估方法,具有如下步骤:
—建立待测定焊接结构的有限元模型;在焊接结构处施加多种实验载荷,求解得到焊趾端面沿板厚方向所有节点的结构应力张量[σ];
—根据节点的结构应力张量[σ]计算节点对应的载荷比例系数R,进而确定应力峰值临界点的位置;使用插值的方法得到所述应力峰值临界点的结构应力张量[σCPSP];
—根据应力峰值临界点的结构应力张量[σCPSP]计算得出疲劳寿命临界损伤面的最大剪应力和正应力范围,应用联合多轴疲劳寿命预测方法MWCM进行多轴疲劳寿命预测。
所述应力峰值临界点定义:
焊趾处切口应力σ(x)在板厚方向的分布分解为三部分:膜应力σm,弯曲应力σb和非线性应力峰值σnlp;计算公式分别如下:
当非线性应力峰值沿着板厚方向衰减为零时,该点的切口只包含所述的膜应力和弯曲应力,该点即为所述的应力峰值临界点。
在应力峰值临界点处焊趾处切口应力σ(x)与外加载荷的关系如下:
所述步骤“—根据节点的结构应力张量[σ]计算节点对应的载荷比例系数R,进而确定应力峰值临界点的位置;”采用临界点拟合公式计算临界点的位置:
式中,xCPSP为临界点到板表面的距离,a=0.06833,b=0.06374;根据该临界点周边节点的结构应力张量,使用插值的方法计算得出临界点的结构应力张量[σCPSP]。
所述建立待测定焊接结构的有限元模型步骤中:在焊趾位置进行局部网格细化。
所述联合多轴疲劳寿命预测方法为修正沃勒曲线法。
由于采用了上述技术方案,本发明提供的焊接结构多轴疲劳寿命评估方法,实现了海洋平台焊接结构从单轴疲劳寿命评估到多轴疲劳寿命评估的跨越,对提高海洋平台焊接结构疲劳寿命评估精度具有重要的理论意义和工程应用价值。其优越性具体体现在以下四个方面:1本发明方法不需要考虑海洋平台焊接结构的焊趾形状及尺寸大小,继承了传统热点应力法的优点,便于应用于工程实际;2本发明方法从临界点位置直接提取结构应力作为评估应力,比传统热点应力法的插值取法具有更明确的物理意义;3本发明方法能够考虑厚度效应对疲劳寿命的影响,通过对疲劳试验数据的计算分析,证明比传统热点应力法具有更高的精度;4传统的热点应力法是通过表面插值计算焊趾处的应力状态,不能计算复杂焊接结构的应力状态,因此无法应用于工程实际。而本发明方法通过有限元计算,可以直接在临界点处提取多轴应力状态,能够应用于实际复杂的海洋平台焊接结构。
附图说明
为了更清楚的说明本发明的实施例或现有技术的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的流程图
图2为本发明实施例1的焊接结构示意图:图中具有经过热点A的纵截面ABEF和横截面ABCD;
图3为图2中纵截面ABEF上正应力分解;
图4为图2中横截面ABCD上剪应力分解;
图5为本发明实施例1中采用的焊接模型几何尺寸示意图;
图6为发明方法的多轴疲劳寿命预测精度对比示意图;
具体实施方式
为使本发明的实施例的目的、技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述:
如图1所示:主要包括如下步骤:一种焊接结构多轴疲劳寿命评估方法,具有如下步骤:
首先建立待测定焊接结构的有限元模型,在焊趾位置进行局部细化,在本发明的实施例中焊趾细化网格尺寸为0.1t×0.1t,其余位置粗网格尺寸0.5t×0.5t。
在焊接结构处施加多种实验载荷,根据试验数据对所建立的有限元模型施加相应的疲劳载荷:纯弯曲、纯扭转、同相加载、非比例加载,具体载荷见附表1:
其中,非比例加载的相位差是90°。非比例载荷加载时采用正弦函数加载,将一个循环周期内的载荷离散为36个载荷工况(间隔10°)。
提取焊趾端面沿板厚方向所有节点的结构应力张量[σ]:
应用膜应力和弯曲应力计算公式,计算对应的载荷比例系数R,然后代入本发明方法通过大量有限元分析计算得到的临界点拟合公式(所述公式5),确定临界点位置以及相应的临界点结构应力。
在确定临界点位置后,使用插值的方法得到临界点的结构应力张量[σCPSP]
根据临界点位置的结构应力张量计算疲劳临界损伤平面——最大剪应力面。
对于给定应力峰值临界点xCPSP,x1≤xCPSP≤x2,在节点的结构应力张量[σ]中分别对应[σ1]、[σCPSP]、[σ2],使用线性插值的方法计算[σCPSP],插值公式如下:
式中,x1和x2为所述的应力峰值临界点xcpsp在板厚度方向的两个相邻节点与板上表面的距离,[σCPSP]即为应力峰值临界点的结构应力张量。
本发明所采用的修正沃勒曲线法MWCM(modifiedcurve method),认为对于焊接结构复杂应力状态可以用最大剪应力范围Δτ与疲劳临界面上的法向正应力范围Δσn的比值来评估,因此使用MWCM进行多轴疲劳寿命预测则必须计算疲劳临界损伤平面以确定最大剪应力范围Δτ与疲劳临界面上的法向正应力范围Δσn。
在本实施例中,采用如图5所示的焊接模型:以1°为间隔,遍历经过临界点的任意空间平面共361×181=65341个平面,计算每个空间平面经历一个加载循环后的最大剪应力范围Δτ及相应的正应力范围Δσn,然后比较所有平面,得到最大剪应力范围所在的平面,即疲劳临界损伤平面。确定最终最大剪应力范围Δτ及相应的正应力范围Δσn。
选用国际公开发表的实验数据(Yousefi F,Witt M,Zenner H.Fatigue strengthof welded joints under multiaxial loading:experiments and calculations,2001)对本方法进行验证。
该方法由Luca Susmel提出,认为对于焊接结构复杂应力状态可以用最大剪应力范围Δτ与疲劳临界面上的法向正应力范围Δσn的比值来评估,如下式所示:
(6)式中:ρω为应力状态;具体而言,在给定循环寿命NA下,MWCM定义焊接结构的多轴疲劳循环次数,如下式所示:预测寿命Nf为:
式中,反向斜率kτ(ρω)和参考剪应力幅值ΔτA,Ref(ρw)由下式确定。
MWCM方法的优点是可以与各种疲劳强度计算方法联合,如名义应力法、热点应力法及切口应力法。其中热点应力法的MWCM控制方程是根据标准疲劳曲线FAT90(单轴拉压曲线)和FAT100(纯扭转曲线)计算得到,在假定失效寿命NA=2×106,FAT90标准差为0.25(IIW建议),生存概率为50%时,其控制方程为:
在双对数坐标系下,以横轴为预测寿命Nf,e、纵轴为真实试验寿命Nf,做精度验证图,并将校准曲线FAT90和FAT100的预测分布带在附图6中画出。
附图6可以看出应用本发明方法得到的寿命预测结果基本落于单轴疲劳寿命预测分布带内部,可知本发明方法可以成功与修正沃勒曲线法MWCM联合进行多轴疲劳寿命预测,并且可以得到满意的预测结果。
附表1
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (5)
1.一种焊接结构多轴疲劳寿命评估方法,具有如下步骤:
—建立待测定焊接结构的有限元模型;在焊接结构处施加多种实验载荷,求解得到焊趾端面沿板厚方向所有节点的结构应力张量[σ];
—根据节点的结构应力张量[σ]计算节点对应的载荷比例系数R,进而确定应力峰值临界点的位置;使用插值的方法得到所述应力峰值临界点的结构应力张量[σCPSP];
—根据应力峰值临界点的结构应力张量[σCPSP]计算得出疲劳寿命临界损伤面的最大剪应力和正应力范围,应用联合多轴疲劳寿命预测方法进行多轴疲劳寿命预测。
2.根据权利要求1所述的焊接结构多轴疲劳寿命评估方法,其特征还在于:所述应力峰值临界点定义:
焊趾处切口应力σ(x)在板厚方向的分布分解为三部分:膜应力σm,弯曲应力σb和非线性应力峰值σnlp;计算公式分别如下:
当非线性应力峰值沿着板厚方向衰减为零时,该点的切口只包含所述的膜应力和弯曲应力,该点即为所述的应力峰值临界点。
3.根据权利要求2所述的焊接结构多轴疲劳寿命评估方法,其特征还在于:在应力峰值临界点处焊趾处切口应力σ(x)与外加载荷的关系如下:
所述步骤“—根据节点的结构应力张量[σ]计算节点对应的载荷比例系数R,进而确定应力峰值临界点的位置;”采用临界点拟合公式计算临界点的位置:
式中,xCPSP为临界点到板表面的距离,a=0.06833,b=0.06374;根据该临界点周边节点的结构应力张量,使用插值的方法计算得出临界点的结构应力张量[σCPSP]。
4.根据权利要求1所述的焊接结构多轴疲劳寿命评估方法,其特征还在于:所述建立待测定焊接结构的有限元模型步骤中:在焊趾位置进行局部网格细化。
5.根据权利要求1所述的焊接结构多轴疲劳寿命评估方法,其特征还在于:所述联合多轴疲劳寿命预测方法为修正沃勒曲线法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310712305.0A CN103714204B (zh) | 2013-12-18 | 2013-12-18 | 焊接结构多轴疲劳寿命评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310712305.0A CN103714204B (zh) | 2013-12-18 | 2013-12-18 | 焊接结构多轴疲劳寿命评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103714204A true CN103714204A (zh) | 2014-04-09 |
CN103714204B CN103714204B (zh) | 2017-04-12 |
Family
ID=50407175
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310712305.0A Active CN103714204B (zh) | 2013-12-18 | 2013-12-18 | 焊接结构多轴疲劳寿命评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103714204B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104200122A (zh) * | 2014-09-22 | 2014-12-10 | 大连交通大学 | 复杂焊接结构随机振动疲劳寿命预测方法 |
CN104408253A (zh) * | 2014-11-26 | 2015-03-11 | 沈阳华创风能有限公司 | 一种风机焊接件疲劳寿命校核的计算方法 |
CN104462790A (zh) * | 2014-11-21 | 2015-03-25 | 南京衍达软件科技有限公司 | 疲劳耐久分析的自由表面方法 |
CN104462834A (zh) * | 2014-12-16 | 2015-03-25 | 中国汽车工程研究院股份有限公司 | 包括焊接模拟的车架复杂工况多轴疲劳计算方法 |
CN105953954A (zh) * | 2016-04-29 | 2016-09-21 | 河海大学 | 一种组合钢桥面板界面母材焊趾疲劳热点应力测试方法 |
CN110232243A (zh) * | 2019-06-12 | 2019-09-13 | 四川大学 | 基于晶体塑性焊接工艺模型的损伤与疲劳寿命评估方法 |
CN110991104A (zh) * | 2019-11-21 | 2020-04-10 | 武汉理工大学 | 焊接空间网架节点焊缝风致多轴高周疲劳损伤评定方法 |
CN112836402A (zh) * | 2021-01-06 | 2021-05-25 | 海洋石油工程股份有限公司 | 海洋工程结构物多轴疲劳应力监测方法 |
CN112834093A (zh) * | 2021-01-08 | 2021-05-25 | 中南大学 | 多点表面应力监测的应力场重构与损伤评估方法及系统 |
CN113111310A (zh) * | 2021-03-29 | 2021-07-13 | 北京理工大学重庆创新中心 | 多受力工况下焊点疲劳寿命测试数据归一化处理方法 |
CN114297893A (zh) * | 2021-12-28 | 2022-04-08 | 北京工业大学 | 一种焊接结构多轴疲劳失效寿命评估方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101592552A (zh) * | 2009-07-01 | 2009-12-02 | 同济大学 | 预测汽车后悬架多轴疲劳的方法 |
US20100235110A1 (en) * | 2009-03-12 | 2010-09-16 | Gm Global Technology Operations, Inc. | Systems and methods to predict fatigue lives of aluminum alloys under multiaxial loading |
-
2013
- 2013-12-18 CN CN201310712305.0A patent/CN103714204B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100235110A1 (en) * | 2009-03-12 | 2010-09-16 | Gm Global Technology Operations, Inc. | Systems and methods to predict fatigue lives of aluminum alloys under multiaxial loading |
CN101592552A (zh) * | 2009-07-01 | 2009-12-02 | 同济大学 | 预测汽车后悬架多轴疲劳的方法 |
Non-Patent Citations (2)
Title |
---|
刘刚等: "基于临界面理论的焊接结构多轴疲劳寿命评估方法", 《船舶力学》 * |
刘刚等: "复杂焊接接头多轴疲劳强度评估的等效热点应力法", 《焊接学报》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104200122A (zh) * | 2014-09-22 | 2014-12-10 | 大连交通大学 | 复杂焊接结构随机振动疲劳寿命预测方法 |
CN104200122B (zh) * | 2014-09-22 | 2017-02-15 | 大连交通大学 | 复杂焊接结构随机振动疲劳寿命预测方法 |
CN104462790A (zh) * | 2014-11-21 | 2015-03-25 | 南京衍达软件科技有限公司 | 疲劳耐久分析的自由表面方法 |
CN104408253A (zh) * | 2014-11-26 | 2015-03-11 | 沈阳华创风能有限公司 | 一种风机焊接件疲劳寿命校核的计算方法 |
CN104462834A (zh) * | 2014-12-16 | 2015-03-25 | 中国汽车工程研究院股份有限公司 | 包括焊接模拟的车架复杂工况多轴疲劳计算方法 |
CN104462834B (zh) * | 2014-12-16 | 2018-01-09 | 中国汽车工程研究院股份有限公司 | 包括焊接模拟的车架复杂工况多轴疲劳计算方法 |
CN105953954A (zh) * | 2016-04-29 | 2016-09-21 | 河海大学 | 一种组合钢桥面板界面母材焊趾疲劳热点应力测试方法 |
CN105953954B (zh) * | 2016-04-29 | 2018-09-04 | 河海大学 | 一种组合钢桥面板界面母材焊趾疲劳热点应力测试方法 |
CN110232243A (zh) * | 2019-06-12 | 2019-09-13 | 四川大学 | 基于晶体塑性焊接工艺模型的损伤与疲劳寿命评估方法 |
CN110232243B (zh) * | 2019-06-12 | 2020-03-10 | 四川大学 | 基于晶体塑性焊接工艺模型的损伤与疲劳寿命评估方法 |
CN110991104A (zh) * | 2019-11-21 | 2020-04-10 | 武汉理工大学 | 焊接空间网架节点焊缝风致多轴高周疲劳损伤评定方法 |
CN110991104B (zh) * | 2019-11-21 | 2023-07-25 | 武汉理工大学 | 焊接空间网架节点焊缝风致多轴高周疲劳损伤评定方法 |
CN112836402A (zh) * | 2021-01-06 | 2021-05-25 | 海洋石油工程股份有限公司 | 海洋工程结构物多轴疲劳应力监测方法 |
CN112834093A (zh) * | 2021-01-08 | 2021-05-25 | 中南大学 | 多点表面应力监测的应力场重构与损伤评估方法及系统 |
CN112834093B (zh) * | 2021-01-08 | 2022-02-11 | 中南大学 | 多点表面应力监测的应力场重构与损伤评估方法及系统 |
CN113111310A (zh) * | 2021-03-29 | 2021-07-13 | 北京理工大学重庆创新中心 | 多受力工况下焊点疲劳寿命测试数据归一化处理方法 |
CN113111310B (zh) * | 2021-03-29 | 2022-09-06 | 北京理工大学重庆创新中心 | 多受力工况下焊点疲劳寿命测试数据归一化处理方法 |
CN114297893A (zh) * | 2021-12-28 | 2022-04-08 | 北京工业大学 | 一种焊接结构多轴疲劳失效寿命评估方法 |
CN114297893B (zh) * | 2021-12-28 | 2023-03-28 | 北京工业大学 | 一种焊接结构多轴疲劳失效寿命评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103714204B (zh) | 2017-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103714204A (zh) | 焊接结构多轴疲劳寿命评估方法 | |
US20210012046A1 (en) | Meshless method for solid mechanics simulation, electronic device, and storage medium | |
CN106886663B (zh) | 齿轮弯曲疲劳寿命预测方法及装置 | |
CN104573392B (zh) | 一种焊点疲劳寿命预测方法 | |
CN103712865A (zh) | 高温部件的蠕变损伤评价方法以及损伤评价系统 | |
Khoei et al. | Crack growth modeling via 3D automatic adaptive mesh refinement based on modified-SPR technique | |
CN104200122A (zh) | 复杂焊接结构随机振动疲劳寿命预测方法 | |
CN106223201A (zh) | 梁段桥梁线形监控的纠偏方法 | |
CN103324782A (zh) | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 | |
CN105260574A (zh) | 一种基于临界面法疲劳破坏准则的高周多轴疲劳寿命预测方法 | |
CN110263401B (zh) | 钢箱梁顶板-纵肋焊接细节残余应力松弛效应的评估方法 | |
CN103559390A (zh) | 一种基于平均失效指数的复合材料π形胶接连接结构拉伸强度预测方法 | |
Khoei et al. | Three-dimensional cohesive fracture modeling of non-planar crack growth using adaptive FE technique | |
CN106021709A (zh) | 一种混凝土早期开裂风险的评估和控制方法 | |
Liang et al. | Estimating inherent deformation in thin-plate Al-alloy joint by means of inverse analysis with the help of cutting technique | |
CN103559361B (zh) | 一种构件强度的优化方法及其应力分析方法 | |
Khodaie et al. | Parametric analyses on the initial stiffness of the SHS column base plate connections using FEM | |
Luca et al. | Estimation of multi-pass welds deformations with Virtual Weld Bead method | |
Chiocca et al. | Rapid and accurate fatigue assessment by an efficient critical plane algorithm: application to a FSAE car rear upright | |
CN105574307A (zh) | 一种展向连接结构dfr双向修正临界点的判断方法 | |
CN102620694B (zh) | 一种机翼精加工姿态计算方法 | |
Sáenz et al. | Analysis and prediction of welding distortion in complex structures using elastic finite element method | |
Amorim | On the lumped damage mechanics for nonlinear structural analyses: new developments and applications | |
CN105136558A (zh) | 一种获取高强钢筋母材及其连接结构的疲劳强度曲线的方法 | |
Lackner et al. | Adaptivity in computational mechanics of concrete structures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |