发明名称 焊接结构多轴疲劳寿命评估方法
摘要 本发明公开了一种焊接结构多轴疲劳寿命评估方法,具有如下步骤:建立待测定焊接结构的有限元模型;在焊接结构处施加多种实验载荷,求解得到焊趾端面沿板厚方向所有节点的结构应力张量[σ];根据节点的结构应力张量[σ]计算节点对应的载荷比例系数R,进而确定应力峰值临界点的位置;使用插值的方法得到所述应力峰值临界点的结构应力张量[σ<sub>CPSP</sub>];根据应力峰值临界点的结构应力张量[σ<sub>CPSP</sub>]计算得出疲劳寿命临界损伤面的最大剪应力和正应力范围,应用联合多轴疲劳寿命预测方法MWCM进行多轴疲劳寿命预测。
申请公布号 CN103714204B 申请公布日期 2017.04.12
申请号 CN201310712305.0 申请日期 2013.12.18
申请人 大连理工大学 发明人 刘刚;黄一;刘英芳;张崎;陈景杰;李红霞;王文华
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 大连东方专利代理有限责任公司 21212 代理人 李馨;李洪福
主权项 一种焊接结构多轴疲劳寿命评估方法,具有如下步骤:—建立待测定焊接结构的有限元模型;在焊接结构处施加多种实验载荷,求解得到焊趾端面沿板厚方向所有节点的结构应力张量[σ];—根据节点的结构应力张量[σ]计算节点对应的载荷比例系数R,进而确定应力峰值临界点的位置;使用插值的方法得到所述应力峰值临界点的结构应力张量[σ<sub>CPSP</sub>];—根据应力峰值临界点的结构应力张量[σ<sub>CPSP</sub>]计算得出疲劳寿命临界损伤面的最大剪应力和正应力范围,应用联合多轴疲劳寿命预测方法进行多轴疲劳寿命预测;所述应力峰值临界点定义:焊趾处切口应力σ(x)在板厚方向的分布分解为三部分:膜应力σ<sub>m</sub>,弯曲应力σ<sub>b</sub>和非线性应力峰值σ<sub>nlp</sub>;计算公式分别如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>m</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>t</mi></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><mi>&sigma;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mi>d</mi><mi>x</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001123008900000011.GIF" wi="1004" he="103" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>b</mi></msub><mo>=</mo><mfrac><mn>6</mn><msup><mi>t</mi><mn>2</mn></msup></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><mi>&sigma;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><mi>t</mi><mn>2</mn></mfrac><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mi>d</mi><mi>x</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001123008900000012.GIF" wi="1038" he="103" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mrow><mi>n</mi><mi>l</mi><mi>p</mi></mrow></msub><mo>=</mo><mi>&sigma;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>m</mi></msub><mo>-</mo><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>x</mi></mrow><mi>t</mi></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>&sigma;</mi><mi>b</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001123008900000013.GIF" wi="1078" he="106" /></maths>式中,t为板厚,x为当前节点与板上表面的距离;当非线性应力峰值沿着板厚方向衰减为零时,该点的切口只包含所述的膜应力和弯曲应力,该点即为所述的应力峰值临界点;在应力峰值临界点处焊趾处切口应力σ(x)与外加载荷的关系如下:<maths num="0004"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mi>m</mi></msub><mrow><msub><mi>&sigma;</mi><mi>m</mi></msub><mo>+</mo><msub><mi>&sigma;</mi><mi>b</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001123008900000014.GIF" wi="964" he="102" /></maths>所述步骤“—根据节点的结构应力张量[σ]计算节点对应的载荷比例系数R,进而确定应力峰值临界点的位置;”采用临界点拟合公式计算临界点的位置:<img file="FDA0001123008900000015.GIF" wi="1246" he="189" />式中,x<sub>CPSP</sub>为临界点到板表面的距离,a=0.06833,b=0.06374;根据该临界点周边节点的结构应力张量,使用插值的方法计算得出临界点的结构应力张量[σ<sub>CPSP</sub>]。
地址 116024 辽宁省大连市高新园区凌工路2号