焊趾处热点应力的计算方法
【技术领域】
[0001] 本发明涉及计算领域，特别地涉及焊趾处热点应力的计算方法。
【背景技术】
[0002] 焊接结构中，焊接接头的固有疲劳强度往往较低。这是由于，焊接特殊的制造工艺 包含了太多的对结构疲劳性能不利的附加因素：连接附件、加强件以及焊缝几何等不连续 结构。这些几何结构上的不连续性导致了结构刚度的不连续性。最终，很容易导致在焊接 接头处产生高度的应力集中。而应力集中的严重程度，决定了焊接结构的疲劳强度。而实 际上，在高应力集中区域，由于应力梯度的高度非线性，很难截取一个统一的具有力学意义 的应力值来考核的结构的疲劳寿命。
[0003] 针对这种问题，国内外学者提出了两种基于S-N曲线的焊接结构的应力截取方 法：名义应力法和基于国际焊接协会的热点应力法。名义应力法是焊接结构疲劳分析中研 究最广、应用最广的传统方法。应用这种方法时，不需要特别考虑焊接接头自身的应力集 中。但是，不同的焊接接头需要通过不同的S-N曲线来区分。从应力截取的角度看，名义应 力法显然比较简单。但是，名义应力法也存在严重的局限性：1)焊接结构的疲劳破坏是焊 接接头处的应力集中而导致，而名义应力法忽略了焊接接头的应力集中。比如，相同的焊接 接头，如果对焊缝进行后处理以降低应力集中，其疲劳寿命肯定是不一样的。2)由于焊接接 头的多样性，焊接接头与S-N曲线的一一对应是很难做到的。3)焊接结构处的名义应力是 高度网格敏感的，应力值会有较大的离散。
[0004] 基于如上的考虑，国际焊接协会推出了热点应力法。热点应力法考虑的焊接接头 的应力集中，它是一种基于几何结构的外推法。但是，热点应力法也不能解决接头与S-N曲 线对应及应力与网格敏感的问题。后来，有学者提出了在有限元中基于节点力的热点应力 计算方法。这种方法一方面不需要考虑单元类型、网格尺寸等敏感因素对热点应力的影响， 另一方面，由于较准确的考虑了焊接接头的应力集中，也突破了焊接接头与S-N曲线对应 的限制。但是此法在实施前需要将焊接结构3D几何模型转化成板壳2D几何模型，对于复 杂钢结构很难准确的用2D板壳模型来表达3D实体模型，且花费时间、人力成本巨大，工程 应用中存在一定的局限性。
【发明内容】
[0005] 本发明提供了一种焊趾处热点应力的计算方法，以解决现有技术中的问题。
[0006] 为实现上述目的，本发明提供了一种焊趾处热点应力的计算方法，包括以下步 骤：S10,有限元模型前处理，包括以下步骤：S11，在焊趾处建立虚拟应变片的几何模型； S12,对所述虚拟应变片的几何模型进行网格划分；S13,将所述虚拟应变片的几何模型分为 应变片外围和应变片核心块，将所述应变片核心块相应的面均分为N个子面，在每个所述 子面的侧边生成1个线特征，所述N个子面共生成（N+1)个所述线特征；S14,为所述虚拟应 变片定义材料、单元属性、施加载荷和边界条件；S2,根据所述有限元模型进行计算；S3,提 取每个节点的线性力和线性力矩；S4,根据S3中提取每个所述节点的所述线性力和所述线 性力矩，计算出所述焊趾处的拉应力、弯曲应力和真实应力；S5,将所述真实应力叠加得到 所述焊趾处每个所述节点的热点应力。
[0007] 优选地，N大于或者等于8。
[0008] 优选地，为所述虚拟应变片建立的网格为3D网格。
[0009] 优选地，当步骤S2完成后，在所述有限元模型的计算结果文件中，提取所述虚拟 应变片中选取的节点在板厚方向的所述线性力与所述线性力矩，其中：
[0010] 真实应力由计算公式
确定，其中：
[0011] f为节点对于相邻一个子面的线性力，m为相应的线性力矩，t是板厚。
[0012] 优选地，将每个所述节点的所述真实应力叠加得出每个所述节点的热点应力：
[0015] 其中，Fi是每个节点的热点应力，1代表单元沿焊缝方向的长度。
[0016] 优选地，根据每个所述节点的所述真实应力和所述热点应力，推广至η个线性力/ 力矩，形成矩阵方程：
[0018] 优选地，将所述矩阵方程中的Fn替换成Mn，fn替换成mn即得到所述线性力矩的 相关矩阵。
[0019] 本发明涉及的焊趾处热点应力的计算方法基于有限元3D模型来计算焊趾处的热 点应力，这种计算方法步骤简单，效率较高，并且计算精度高，极大地降低了人力物力，提高 了焊接结构的后期设计验证的效率和焊接结构的设计质量。
【附图说明】
[0020] 图1是本发明涉及的一种焊趾处热点应力的计算方法中的焊趾位置的示意图；
[0021] 图2是本发明涉及的一种焊趾处热点应力的计算方法的步骤流程图；
[0022] 图3是本发明涉及的一种焊趾处热点应力的计算方法中的虚拟应变片网格划分 示意图；
[0023] 图4是本发明涉及的一种焊趾处热点应力的计算方法中的虚拟应变片示意图；
[0024] 图5是本发明涉及的一种焊趾处热点应力的计算方法中的虚拟应变片示意图；
[0025] 图6是本发明涉及的一种焊趾处热点应力的计算方法中的虚拟应变片3D示意 图；
[0026] 图7是本发明涉及的一种焊趾处热点应力的计算方法中的虚拟应变片局部示意 图；
[0027] 图8是本发明涉及的一种焊趾处热点应力的计算方法中线性力和线性力矩示意 图；
[0028] 图9是本发明涉及的一种焊趾处热点应力的计算方法中线性力和线性力矩示意 图。
【具体实施方式】
[0029] 为了更好地说明本发明的意图，下面结合附图对本
【发明内容】
做进一步说明。
[0030] 本实施例涉及一种焊趾处热点应力的计算方法，其中，如图1所示，焊趾处是指焊 接时，焊缝1与图1中母材的工字钢梁2两个实体相交的区域，在进行计算时，在相交区域 选取虚拟应变片区域3,如图2所示，在对焊趾处进行热点应力的计算方法包括以下步骤：
[0031] S10,有限元模型前处理，包括以下步骤：
[0032] S11，利用有限元分析软件或3D软件在焊趾处建立虚拟应变片的几何模型；
[0033] S12,如图3和图6所示，对所述虚拟应变片的几何模型进行网格划分，具体地，将 几何模型划分为3D网格；
[0034] S13,如图4和图5所示，将所述虚拟应变片的几何模型分为应变片外围4和应变 片核心块5,将所述应变片核心块5相应的面均分为N个子面，在每个所述子面的侧边生成 1个线特征，所述N个子面共生成（N