一种焊接结构疲劳寿命的动态结构应力计算方法技术

本发明专利技术一种焊接结构疲劳寿命的动态结构应力计算方法，包括以下步骤：首先建立关键部件有限元模型，选定关键部件的焊缝局部模型的边界单元及节点子集，采用改进的Craig

【技术实现步骤摘要】
一种焊接结构疲劳寿命的动态结构应力计算方法


[0001]本专利技术属于全自动化产品领域，涉及一种焊接结构疲劳寿命的动态结构应力计算方法。

技术介绍

[0002]现有的焊接结构动态结构应力计算方法及疲劳评估技术包括以下三种：
[0003]1)名义应力法：提取距离焊缝一定距离处的名义应力来评估焊接结构的疲劳寿命，具体的应用可以基于准静态法和动力学计算结果，然后从IIW、BS7608、JIS等标准种选取与所分析焊接结构对应的结构类型及载荷模式的名义应力S
‑
N曲线，进行疲劳评估[1
‑
4]。
[0004]2)准静态结构应力法：准静态结构应力法主要是基于单位载荷下结构的静力学计算结果，计算对应的结构应力，假设载荷与应力直线存在线性关系，基于动载荷数据将静结构应力转换动态结构应力，对结构应力进行雨流技术利用S
‑
N曲线评估焊缝的疲劳
[5
‑
8]。
[0005]3)模态结构应力法:该方法首先要建立带焊缝细节的有限元模型，计算并截取部分低阶模态，然后将动力学方程转换到模态坐标下，计算求解动力学方程获得模态坐标，计算焊缝焊趾处的模态结构应力，采用结构应力与模态坐标进行叠加获得结构应力时间历程，经雨流计数后可以采用主S
‑
N曲线评估焊缝疲劳寿命[9]。
[0006]名义应力法依赖的这些标准种的焊接结构类型是有限的,很多情况是无法找到与工程结构一致的焊接接头类型及受力方向，因此，评估的结果会因人而异，不具唯一性。
[0007]准静态结构应力法忽略结构的振动响应，无法考虑载荷的频率对结构模态振动及焊缝疲劳寿命的影响。当外载荷的激励频率远远低于结构的最第阶模态频率的情况下，可以获得较好的结果；但是，当外部载荷激励与结构的模态频率接近或者存在交集的情况下，载荷幅值与应力之间存在非线性关系，准静态的评估结果会出现较大偏差。
[0008]模态结构应力法虽可以考虑结构振动的影响，但是由于截取的是部分低阶模态，低阶模态振型一般描述的是结构的整体变形；而焊缝局部几何的变形属于高频的模态振型，即使在模型中建立焊缝细节，但高阶模态截断后的部分低阶模态不能较好的模拟焊缝局部变形，因此焊缝局部应力的计算会存一定误差。

技术实现思路

[0009]为了解决上述问题，本专利技术提供本专利技术采用的技术方案是：主要是为了解决，刚柔耦合动力学仿真模型中柔性体焊缝局部应力集中难以有效模拟，计算精度低，柔性体建模复杂等困难；一种焊接结构疲劳寿命的动态结构应力计算方法，包括以下步骤：
[0010]S1:首先建立关键部件有限元模型，选定关键部件的焊缝局部模型的边界单元及节点子集，采用改进的Craig
‑
Bampton模态综合法计算获得柔性体，定义焊缝边界单元和节点子集并输出边界子集的模态节点力；
[0011]S2:建立刚柔耦合动力学模型，施加相关载荷进行动力学仿真计算，并输出柔性体
模态坐标时间历程；
[0012]S3:基于模态叠加原理，通过关键部件的某条焊缝的模型边界模态节点力与模态坐标的叠加获得焊缝局部边界子集的动态节点力时间历程；
[0013]S4:将某条焊缝的模型进行细化，建立焊缝细节并细化局部模型，保留边界处节点数量和位置与边界子集一一对应；
[0014]S5:将边界子集的动态节点力施加到细化的焊缝模态，采用静态自由度凝聚与惯性释放相结合的方法，计算获得焊缝子模型的变形及节点力；
[0015]S6:基于结构应力法计算公式针对每一个时刻计算对应的焊趾处的结构应力，从而获得动态结构应力随时间的变化；
[0016]S7:将动态结构应力经雨流技术统计应力变化范围及循环次数，基于ASME标准公式计算对应的等效结构应力变化范围。
[0017]进一步地：所述采用改进的Craig
‑
Bampton模态综合法计算获得柔性体，定义焊缝边界单元和节点子集并输出边界子集的模态节点力的过程如下：
[0018]将每阶模态对应的节点力进行叠加处理，需要将同一节点的节点力进行叠加，形成模态边界节点力集合[F]，另外在装配节点力之前，需要将不同单元对同一个节点的节点力进行数据加和处理；
[0019][0020]其中：f
xmn
表示第m阶模态n号节点的x方向节点力，每个节点包含六个自由度方向的节点力，即f
x
，f
y
，f
z
三个平动方向的力和m
x
，m
y
，m
z
三个转动方向的力矩，集合后的节点力矩阵元素的总数＝边界节点数n
×6×
模态数m。如果有存在焊缝边界单元集合中有多个单元对某个相同节点有节点力作用，那么首先需要将节点力在全局坐标下进行合并处理即可。
[0021]进一步地：所述建立刚柔耦合动力学模型，的过程如下：
[0022]所述刚柔耦合动力学模型通过将柔性体文件与系统刚体进行连接组成，模态坐标的数量与柔性体定义选取的模态坐标数量一致；柔性车体与刚体之间通过柔性体外接口处的约束方程与多刚体系统建立联系，并通过拉格朗日乘子法将约束条件引入到系统方程中，通过与多刚体系统方程进行联立可以获得刚柔耦合系统动力学方程：
[0023][0024]式(2)中，Ψ(ξ,t)为约束方程；ξ为广义坐标，包括位移坐标X、欧拉角坐标ω和模态坐标q；Q为包括有势力在内的广义力；L为拉格朗日函数；λ为待定因子；F为如下形式的耗散函数，采用New
‑
mark积分法求解动力学方程，获得柔性体模态坐标时间历程q
i
(t),(i＝0,1,
…
m)。
[0025]进一步地：所述通过关键部件的某条焊缝的模型的边界模态节点力与模态坐标的叠加获得焊缝局部边界子集的动态节点力时间历程表达式如下：
[0026][0027]其中：q
i
(t)带表第i阶柔性体模态坐标时间历程，为模态坐标向量，{f
b
(t)}
6n
×1为子模型边界节点力向量随时间的变化。
[0028]进一步地：将某条焊缝的模型进行细化，建立焊缝细节并细化局部模型，保留边界处节点数量和位置与边界子集一一对应；
[0029]S5:将边界子集的动态节点力施加到细化的焊缝模态，采用静态自由度凝聚与惯性释放相结合的方法，计算获得焊缝子模型全部节点的位移及节点力的过程如下：
[0030]将焊缝子模型的刚度矩阵，对应内部自由度u
i
和边界自由度u
b
进行分块表达如下：
[0031][0032]其中，f
b
(t)为计算出的边界载荷，由于焊缝子模型内部不受外载荷作用，基于式(4)的第一式得出内部自由度和边界自由度之间的关系如下：
[0033]{u
i
(t)}＝...

【技术保护点】

【技术特征摘要】
1.一种焊接结构疲劳寿命的动态结构应力计算方法，其特征在于：包括以下步骤：S1:首先建立关键部件有限元模型，选定关键部件的焊缝局部模型的边界单元及节点子集，采用改进的Craig
‑
Bampton模态综合法计算获得柔性体，定义焊缝边界单元和节点子集并输出边界子集的模态节点力；S2:建立刚柔耦合动力学模型，施加相关载荷进行动力学仿真计算，并输出柔性体模态坐标时间历程；S3:基于模态叠加原理，通过关键部件的某条焊缝的模型边界模态节点力与模态坐标的叠加获得焊缝局部边界子集的动态节点力时间历程；S4:将某条焊缝的模型进行细化，建立焊缝细节并细化局部模型，保留边界处节点数量和位置与边界子集一一对应；S5:将边界子集的动态节点力施加到细化的焊缝模态，采用静态自由度凝聚与惯性释放相结合的方法，计算获得焊缝子模型的变形及节点力；S6:基于结构应力法计算公式针对每一个时刻计算对应的焊趾处的结构应力，从而获得动态结构应力随时间的变化；S7:将动态结构应力经雨流技术统计应力变化范围及循环次数，基于ASME标准公式计算对应的等效结构应力变化范围。2.根据权利要求1所述的一种焊接结构疲劳寿命的动态结构应力计算方法，其特征在于：所述采用改进的Craig
‑
Bampton模态综合法计算获得柔性体，定义焊缝边界单元和节点子集并输出边界子集的模态节点力的过程如下：将每阶模态对应的节点力进行叠加处理，需要将同一节点的节点力进行叠加，形成模态边界节点力集合[F]，另外在装配节点力之前，需要将不同单元对同一个节点的节点力进行数据加和处理；其中：f
xmn
表示第m阶模态n号节点的x方向节点力，每个节点包含六个自由度方向的节点力，即f
x
，f
y
，f
z
三个平动方向的力和m
x
，m
y
，m
z
三个转动方向的力矩，集合后的节点力矩阵元素的总数＝边界节点数n
×6×
模态数m。如果有存在焊缝边界单元集合中有多个单元对某个相同节点有节点力作用，那么首先需要将节点力在全局坐标下进行合并处理即可。3.根据权利要求1所述的一种焊接结构疲劳寿命的动态结构应力计算方法，其特征在于：所述建立刚柔耦合动力学模型，的过程如下：所述刚柔耦合动力学模型通过将柔性体文件与系统刚体进行连接组成，模态坐标的数量与柔性体定义选取的模态坐标数量一致；柔性车体与刚体之间通过柔性体外接口处的约束方程与多刚体系统建立联系，并通过拉格朗日乘子法将约束条件引入到系统方程中，通过与多刚体系统方程进行联立可以获得刚柔耦合系统动力学方程：
式(2)中，Ψ(ξ,t)为约束方程；ξ为广义坐标，包括位移坐标X、欧拉角坐标ω和模态坐标q；Q为包括有势力在内的广义力；L为拉格朗日函数；λ为待定因子；F为如下形式的耗散函数，采用New
‑
mark积分法求解动力学方程，获得柔性体模态坐标时间历程q
i
(t),(i＝0,1,
…
m)。4.根据权利要求1所述的一种焊接结构疲劳寿命的动态结构应力计算方法，其特征在于：所述通过关键部件的某条焊缝的模型的边界模态节点力与模态坐标的叠加获得焊缝局部边界子集的动态节点力时间历程表达式如下：其中：q
i
(t)带表第i阶柔性体模态坐标时间历程，为模态坐标向量，{f
b
(t)}
6n
×1为子模型边界节点力向量随时间的变化。5.根据权利要求1所述的一种焊接结构疲劳寿命的动态结构应力计算方法，其特征在于：所述S4:将某条焊缝的模型进行细化，建立焊缝细节并细化局部模型，保留边界处节点数量和位置与边界子集一一对应；S5:将边界子集的动态节点力施加到细化的焊缝模态，采用静态自由度凝聚与惯性释放相结合的方法，计算获得焊缝子模型全部节点的位移及节点力的过程如下：将焊缝子模型的刚度矩阵，对应内部自由度u
i

【专利技术属性】
技术研发人员：方吉，刘昊霖，彭瑞清，
申请(专利权)人：大连交通大学，
类型：发明
国别省市：