CN110717291A - 一种焊接结构件变形仿真方法、装置、设备及存储介质 - Google Patents

一种焊接结构件变形仿真方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN110717291A
CN110717291A CN201910918981.0A CN201910918981A CN110717291A CN 110717291 A CN110717291 A CN 110717291A CN 201910918981 A CN201910918981 A CN 201910918981A CN 110717291 A CN110717291 A CN 110717291A
Authority
CN
China
Prior art keywords
welding
finite element
plastic strain
coordinate system
element model
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
Application number
CN201910918981.0A
Other languages
English (en)
Other versions
CN110717291B (zh
Inventor
庞盛永
梁吕捷
王靖升
李文
黄安国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201910918981.0A priority Critical patent/CN110717291B/zh
Publication of CN110717291A publication Critical patent/CN110717291A/zh
Application granted granted Critical
Publication of CN110717291B publication Critical patent/CN110717291B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Resistance Welding (AREA)

Abstract

本发明公开了一种焊接结构件变形仿真方法、装置、设备及存储介质，方法包括：获取结构件的焊接接头有限元模型及结构件的有限元模型；对焊接接头有限元模型进行热弹塑性有限元分析，以获得在焊接接头局部坐标系下的焊接残余塑性应变；将焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变；其焊接塑性应变分布于在全局坐标系下的预定平面；获取焊接接头有限元模型及结构件的有限元模型在预定平面的投影信息，并根据投影信息在预定平面上形成计算区域；根据在预定平面上的焊接塑性应变在计算区域的分布，得到结构件对应的塑性应变，能够准确的预测大型复杂结构件的变形，为实际中的结构件提供设计和优化参考。

Description

一种焊接结构件变形仿真方法、装置、设备及存储介质
技术领域
本发明涉及焊接数值模拟应用领域，尤其涉及一种焊接结构件变形仿真方法、装置、设备以及存储介质。
背景技术
焊接是制造业中非常重要而关键的领域，大型复杂结构件的焊接是工程中的前沿课题之一。在大型复杂结构件的焊接中，往往存在焊接变形难预测，难控制，难矫正等问题。大型复杂结构件焊接变形仿真技术可以从设计端到优化端解决上述问题，并提供数字化的解决方案。然而，大型复杂焊接结构件变形仿真计算量巨大，计算精度难以保证。其计算的精度和效率完全取决于核心的算法和求解技术。目前，以固有应变算法为核心的求解技术精度有限，模拟变形形式过于单一，仅能计算焊接收缩变形和角变形，基于热-弹塑性理论的求解方法计算量巨大，计算时间难以忍受，使得超大型超复杂焊接结构变形仿真仍缺乏核心算法。
发明内容
针对上述问题，本发明的目的在于提供一种焊接结构件变形仿真方法、装置、设备以及存储介质，使得整体结构件大量的计算化整为零，避免了对大量的网格进行复杂的非线性计算，极大地提升了大型复杂结构件焊接变形仿真的速度，并准确的预测大型复杂结构件的变形，为实际中的结构件提供设计和优化参考。
第一方面，本发明实施例提供了一种焊接结构件变形仿真方法，包括：
获取结构件的焊接接头有限元模型以及结构件的有限元模型；
对所述焊接接头有限元模型进行热弹塑性有限元分析，以获得在焊接接头的局部坐标系下的焊接残余塑性应变；
将处于局部坐标系下的焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变；其中，所述焊接塑性应变分布于在所述全局坐标系下的预定平面；所述预定平面为结构件焊缝垂面；
获取焊接接头有限元模型以及结构件的有限元模型在所述预定平面的投影信息，并根据所述投影信息在所述预定平面上形成计算区域；
根据在所述预定平面上的焊接塑性应变在所述计算区域的分布，得到结构件对应的塑性应变，以预测结构件的变形情况。
优选地，在所述获取结构件的焊接接头有限元模型以及结构件有限元模型的步骤之前，还包括：
获取预先建立的结构件的三维几何模型；
根据所述结构件的三维几何模型，获取结构件的焊接接头几何模型；
根据所述焊接接头几何模型进行网络划分，生成焊接接头有限元模型；
根据所述结构件的三维几何模型进行网络划分，生成结构件有限元模型；。
优选地，将处于局部坐标系下的焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变，具体为：
将处于局部坐标系下的焊接接头塑性应变从焊接接头焊缝坐标系变换至结构件焊缝坐标系，以获得在结构件全局坐标系下的焊接塑性应变；其中，焊缝坐标系是基于焊缝标志线及与所述焊缝标志线对应的焊枪方向所建立的笛卡尔坐标系。
优选地，所述焊接塑性应变分布于在所述全局坐标系下的预定平面，则获取焊接接头有限元模型以及结构件的有限元模型在所述预定平面投影信息，并根据所述投影信息在所述预定平面上形成计算区域，具体为：
基于焊接接头有限元模型在所述预定平面的第一投影信息，对所述预定平面进行划分，以在所述预定平面上获得母区域；
基于结构件的有限元模型在所述预定平面的第二投影信息，对所述母区域进行划分，以在所述预定平面上形成计算区域；其中，计算区域包括多个子区域。
优选地，根据在所述预定平面上的焊接塑性应变在所述计算区域的分布，得到结构件对应的塑性应变，以预测结构件的变形情况，具体为；
对每一子区域内的焊接塑性应变进行积分求和，以得到每个子区域的积分点塑性应变值；
根据所述子区域的积分点塑性应变值，对母区域进行积分反求，得到结构件母区域对应的焊缝塑性应变，以预测结构件的变形情况。
优选地，还包括：
采用Inteweld软件或者开源可视化程序对映射完成后所获得的数据进行可视化处理。
本发明实施例还提供了一种焊接结构件变形仿真方法装置，包括：
有限元模型获取单元，用于获取结构件的焊接接头有限元模型以及结构件的有限元模型；
焊接残余塑性应变获取单元，用于对所述焊接接头有限元模型进行热弹塑性有限元分析，以获得在焊接接头的局部坐标系下的焊接残余塑性应变；
坐标变换单元，用于将处于局部坐标系下的焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变；其中，所述焊接塑性应变分布于在所述全局坐标系下的预定平面；所述预定平面为结构件焊缝垂面；
计算区域获取单元，用于获取焊接接头有限元模型以及结构件的有限元模型在所述预定平面的投影信息，并根据所述投影信息在所述预定平面上形成计算区域；
变形情况预测单元，用于根据在所述预定平面上的焊接塑性应变在所述计算区域的分布，得到结构件对应的塑性应变，以预测结构件的变形情况。
优选地，还包括：
三维几何模型获取单元，用于获取预先建立的结构件的三维几何模型；
焊接接头几何模型获取单元，用于根据所述结构件的三维几何模型，获取结构件的焊接接头几何模型；
第一生成单元，用于根据所述焊接接头几何模型进行网络划分，生成焊接接头有限元模型；
第二生成单元，用于根据所述结构件的三维几何模型进行网络划分，生成结构件有限元模型。
优选地，坐标变换单元，具体用于：
将处于局部坐标系下的焊接接头塑性应变从焊接接头焊缝坐标系变换至结构件焊缝坐标系，以获得在结构件全局坐标系下的焊接塑性应变；其中，焊缝坐标系是基于焊缝标志线及与所述焊缝标志线对应的焊枪方向所建立的笛卡尔坐标系。
优选地，计算区域获得单元，具体包括：
母区域获取模块，用于基于焊接接头有限元模型在所述预定平面的第一投影信息，对所述预定平面进行划分，以在所述预定平面上获得母区域；
计算区域形成模块，用于基于结构件的有限元模型在所述预定平面的第二投影信息，对所述母区域进行划分，以在所述预定平面上形成计算区域；其中，计算区域包括多个子区域。
优选地，变形情况预测单元，具体包括；
积分求和模块，用于对每一子区域内的焊接塑性应变进行积分求和，以得到每个子区域的积分点塑性应变值；
积分反求模块，用于根据所述子区域的积分点塑性应变值，对母区域进行积分反求，得到结构件母区域对应的焊缝塑性应变，以预测结构件的变形情况。
优选地，还包括：
采用Inteweld软件或者开源可视化程序对映射完成后所获得的数据进行可视化处理。
第三方面，本发明实施例还提供了一种焊接结构件变形仿真方法设备，包括处理器、存储器以及存储在所述存储器内的计算机程序，所述计算机程序能够被所述处理器执行以实现如第一方面所述的焊接结构件变形仿真方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如上述所述的焊接结构件变形仿真方法。
实施本发明实施例，具有如下有益效果：
在上述一实施例中，采用Local-Global的思想，将局部焊接接头在热-弹塑性有限元分析后获得的焊接残余塑性应变快速映射至大型结构件的全局坐标系中，使得整体结构件大量的计算化整为零，避免了对大量的网格进行复杂的非线性计算，极大地提升了大型复杂结构件焊接变形仿真的速度，并最大程度的还原了原始塑性应变的信息，为准确的预测大型复杂焊接结构件变形提供了最重要的保证，为设计和优化工艺参数提供参考。
附图说明
为了更清楚地说明本发明的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明第一实施例提供的焊接结构件变形仿真方法的流程示意图。
图2是本发明实施例提供的焊接接头的焊缝标志线及与所述焊缝标志线对应的焊枪方向所建立的焊缝坐标系的结构示意图。
图3本发明实施例的焊接接头焊缝坐标系，结构件焊缝坐标系和全局坐标系的关系示意图。
图4为本发明实施例基于焊接接头有限元模型的第一投影信息的投影示意图。
图5为本发明实施例的结构件有限元模型的第二投影信息的投影示意图。
图6是本发明实施例的投影区域求交示意图。
图7是本发明实施例的对子区域塑性应变积分求和示意图。
图8是本发明第二实施例提供的焊接结构件变形仿真方法装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为了更好的理解本发明的技术方案，下面结合附图对本发明实施例进行详细描述。
应当明确，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。
应当理解，本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”或“响应于检测”。类似地，取决于语境，短语“如果确定”或“如果检测(陈述的条件或事件)”可以被解释成为“当确定时”或“响应于确定”或“当检测(陈述的条件或事件)时”或“响应于检测(陈述的条件或事件)”。
实施例中提及的“第一\第二”仅仅是是区别类似的对象，不代表针对对象的特定排序，可以理解地，“第一\第二”在允许的情况下可以互换特定的顺序或先后次序。应该理解“第一\第二”区分的对象在适当情况下可以互换，以使这里描述的实施例能够以除了在这里图示或描述的那些以外的顺序实施。
实施例一：
请参阅图1至图5，本发明第一实施例提供了一种焊接结构件变形仿真方法，其可由焊接结构件变形仿真方法设备(以下简称仿真设备)来执行，特别的，由仿真设备内的一个或多个处理器来执行，并至少包括如下步骤：
S101，获取结构件的焊接接头有限元模型以及结构件的有限元模型。
在本实施例中，仿真设备首先获取预先建立的结构件的三维几何模型；根据所述结构件的三维几何模型，获取结构件的焊接接头几何模型，然后对焊接接头几何模型进行网格划分后生成焊接接头有限元模型，以及对三维几何模型进行网络划分后生成结构件的有限元模型。需要说明的是，对焊接接头几何模型进行网格划分得到的焊接接头有限元模型，所采用的单元形式为实体单元；对三维几何模型进行网格划分得到的壳体结构件的有限元模型，所采用的单元形式也为实体单元。
S102，对所述焊接接头有限元模型进行热弹塑性有限元分析，以获得在焊接接头的局部坐标系下的焊接残余塑性应变。
在本实施例中，仿真设备采用Inteweld焊接模拟软件对焊接接头有限元模型进行热弹塑性有限元分析，以获得在焊接接头的局部坐标系下的焊接残余塑性应变，其中，所述Inteweld是一款焊接结构应力变形工艺仿真软件。
S103，将处于局部坐标系下的焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变；其中，所述焊接塑性应变分布于在所述全局坐标系下的预定平面；所述预定平面为结构件焊缝垂面。
在本实施例中，将处于局部坐标系下的焊接接头塑性应变从焊接接头焊缝坐标系变换至结构件焊缝坐标系，以获得在结构件全局坐标系下的焊接塑性应变；其中，如图2所示，局部坐标系下的焊缝坐标系为基于焊缝标志线及与所述焊缝标志线对应的焊枪方向所建立的笛卡尔坐标系，具体地，所述焊缝坐标系是以焊缝上某点为原点，以焊缝和与之垂直的焊枪方向为两轴，用右手螺旋法则确定第三轴，所建立的笛卡尔坐标系。当然，可以理解的是，结构件的焊缝坐标系与焊接接头焊缝坐标系的建立方法一致。则焊接接头焊缝坐标系，结构件焊缝坐标系和全局坐标系的关系如图3所示。
S104，获取焊接接头有限元模型以及结构件的有限元模型在所述预定平面的投影信息，并根据所述投影信息在所述预定平面上形成计算区域。
在本实施例中，基于焊接接头有限元模型在所述预定平面的第一投影信息为焊接接头有限元模型在三维空间中以六面体为单元的应变分布转变为二维平面中以四边形为单元的应变分布，如图4所示。基于结构件的有限元模型在所述预定平面的第二投影信息为遍历大型结构件焊缝区域的单元，对于每一个单元，将其信息投影至结构件结构件焊缝垂面，如图5所示。
其中，首先基于焊接接头有限元模型在所述预定平面的第一投影信息，对所述预定平面进行划分，以在所述预定平面上获得母区域。例如，如图6所示的点A1-A4所围成的区域为母区域。然后基于结构件的有限元模型在所述预定平面的第二投影信息，对所述母区域进行划分，以在所述预定平面上形成计算区域；其中，计算区域包括多个子区域，例如，如图6所示的用点b1-b10的连线对母区域进行划分，得到图6中若干个S或者S’的子区域。
S105，根据在所述预定平面上的焊接塑性应变在所述计算区域的分布，得到结构件对应的塑性应变，以预测结构件的变形情况。
在本实施例中，首先对每一子区域内的焊接塑性应变进行积分求和，以得到每个子区域的积分点塑性应变值，例如如图7中每一子区域Φi(i＝1,2…8)内的接头塑性应变
Figure BDA0002216954360000081
进行积分并求和，以得到ESum 然后根据所述子区域的积分点塑性应变值，对母区域进行积分反求，得到结构件母区域对应的焊缝塑性应变εglobal；以预测结构件的变形情况；其中，Esum＝∫Φ{εglobal}dV。
综上，采用Local-Global的思想，将局部焊接接头在热-弹塑性有限元分析后获得的焊接残余塑性应变快速映射至大型结构件的全局坐标系中，使得整体结构件大量的计算化整为零，避免了对大量的网格进行复杂的非线性计算，极大地提升了大型复杂结构件焊接变形仿真的速度，并最大程度的还原了原始塑性应变的信息，为准确的预测大型复杂焊接结构件变形提供了最重要的保证，为设计和优化工艺参数提供参考。
在上述实施例的基础上，本发明一优选实施例中，还包括：采用Inteweld软件或者开源可视化程序对映射完成后所获得的数据进行可视化处理。在本实施例中，采用Inteweld软件或者开源可视化Paraview对模拟过程中所获得的模拟数据进行可视化处理，对计算得到塑性应变场进行可视化分析，研究判断计算结果的好坏，提取各项指标和曲线如塑性应变数据等，通过采用C++语言编写了相应的物理场求解器，采用Inteweld软件或者开源可视化Paraview对模拟过程中所获得的模拟数据进行可视化处理，并基于VTK和QT平台进行了可视化界面开发，使得分析人员利用计算机便可获取整体结构件变形的结果，为设计和优化工艺参数提供参考，为焊接工艺制定提供了重要工具。
本发明第二实施例：
参见图8，本发明实施例还提供了一种焊接结构件变形仿真方法装置，包括：
有限元模型获取单元100，用于获取结构件的焊接接头有限元模型以及结构件的有限元模型；
焊接残余塑性应变获取单元200，用于对所述焊接接头有限元模型进行热弹塑性有限元分析，以获得在焊接接头的局部坐标系下的焊接残余塑性应变；
坐标变换单元300，用于将处于局部坐标系下的焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变；其中，所述焊接塑性应变分布于在所述全局坐标系下的预定平面；所述预定平面为结构件焊缝垂面；
计算区域获取单元400，用于获取焊接接头有限元模型以及结构件的有限元模型在所述预定平面的投影信息，并根据所述投影信息在所述预定平面上形成计算区域；
变形情况预测单元500，用于根据在所述预定平面上的焊接塑性应变在所述计算区域的分布，得到结构件对应的塑性应变，以预测结构件的变形情况。
优选地，还包括：
三维几何模型获取单元，用于获取预先建立的结构件的三维几何模型；
焊接接头几何模型获取单元，用于根据所述结构件的三维几何模型，获取结构件的焊接接头几何模型；
第一生成单元，用于根据所述焊接接头几何模型进行网络划分，生成焊接接头有限元模型；
第二生成单元，用于根据所述结构件的三维几何模型进行网络划分，生成结构件有限元模型。
在上述实施例的基础上，本发明一优选实施例中，坐标变换单元300，具体用于：
将处于局部坐标系下的焊接接头塑性应变从焊接接头焊缝坐标系变换至结构件焊缝坐标系，以获得在结构件全局坐标系下的焊接塑性应变；其中，焊缝坐标系是基于焊缝标志线及与所述焊缝标志线对应的焊枪方向所建立的笛卡尔坐标系。
在上述实施例的基础上，本发明一优选实施例中，计算区域获得单元400，具体包括：
母区域获取模块，用于基于焊接接头有限元模型在所述预定平面的第一投影信息，对所述预定平面进行划分，以在所述预定平面上获得母区域；
计算区域形成模块，用于基于结构件的有限元模型在所述预定平面的第二投影信息，对所述母区域进行划分，以在所述预定平面上形成计算区域；其中，计算区域包括多个子区域。
在上述实施例的基础上，本发明一优选实施例中，变形情况预测单元500，具体包括；
积分求和模块，用于对每一子区域内的焊接塑性应变进行积分求和，以得到每个子区域的积分点塑性应变值；
积分反求模块，用于根据所述子区域的积分点塑性应变值，对母区域进行积分反求，得到结构件母区域对应的焊缝塑性应变，以预测结构件的变形情况。
在上述实施例的基础上，本发明一优选实施例中，还包括：
采用Inteweld软件或者开源可视化程序对映射完成后所获得的数据进行可视化处理。
本发明第三实施例：
本发明第三实施例提供了一种焊接结构件变形仿真方法设备，包括处理器、存储器以及存储在所述存储器内的计算机程序，所述计算机程序能够被所述处理器执行以实现如上述所述的焊接结构件变形仿真方法。
本发明第四实施例：
本发明第四实施例提供了一种计算机可读存储介质，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如上述所述的焊接结构件变形仿真方法。
示例性的，所述计算机程序可以被分割成一个或多个单元，所述一个或者多个单元被存储在所述存储器中，并由所述处理器执行，以完成本发明。所述一个或多个单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序在焊接结构件变形仿真方法设备中的执行过程。
所述焊接结构件变形仿真方法设备可包括但不仅限于处理器、存储器。本领域技术人员可以理解，所述示意图仅仅是焊接结构件变形仿真方法设备的示例，并不构成对焊接结构件变形仿真方法设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述焊接结构件变形仿真方法设备还可以包括输入输出设备、网络接入设备、总线等。
所称处理器可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述焊接结构件变形仿真方法设备的控制中心，利用各种接口和线路连接整个焊接结构件变形仿真方法设备的各个部分。
所述存储器可用于存储所述计算机程序和/或模块，所述处理器通过运行或执行存储在所述存储器内的计算机程序和/或模块，以及调用存储在存储器内的数据，实现所述焊接结构件变形仿真方法设备的各种功能。所述存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等；存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(SecureDigital,SD)卡，闪存卡(Flash Card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
其中，所述焊接结构件变形仿真方法设备集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。
需说明的是，以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外，本发明提供的装置实施例附图中，模块之间的连接关系表示它们之间具有通信连接，具体可以实现为一条或多条通信总线或信号线。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。

Claims (10)

1.一种焊接结构件变形仿真方法，其特征在于，包括：
获取结构件的焊接接头有限元模型以及结构件的有限元模型；
对所述焊接接头有限元模型进行热弹塑性有限元分析，以获得在焊接接头局部坐标系下的焊接残余塑性应变；
将处于局部坐标系下的焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变；其中，所述焊接塑性应变分布于在所述全局坐标系下的预定平面；所述预定平面为结构件焊缝垂面；
获取焊接接头有限元模型以及结构件的有限元模型在所述预定平面的投影信息，并根据所述投影信息在所述预定平面上形成计算区域；
根据在所述预定平面上的焊接塑性应变在所述计算区域的分布，得到结构件对应的塑性应变，以预测结构件的变形情况。
2.根据权利要求1所述的焊接结构件变形仿真方法，其特征在于，在所述获取结构件的焊接接头有限元模型以及结构件有限元模型的步骤之前，还包括：
获取预先建立的结构件的三维几何模型；
根据所述结构件的三维几何模型，获取结构件的焊接接头几何模型；
根据所述焊接接头几何模型进行网络划分，生成焊接接头有限元模型；
根据所述结构件的三维几何模型进行网络划分，生成结构件有限元模型。
3.根据权利要求1所述的焊接结构件变形仿真方法，其特征在于，将处于局部坐标系下的焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变，具体为：
将处于局部坐标系下的焊接接头塑性应变从焊接接头焊缝坐标系变换至结构件焊缝坐标系，以获得在结构件全局坐标系下的焊接塑性应变；其中，焊缝坐标系是基于焊缝标志线及与所述焊缝标志线对应的焊枪方向所建立的笛卡尔坐标系。
4.根据权利要求1所述的焊接结构件变形仿真方法，其特征在于，所述焊接塑性应变分布于在所述全局坐标系下的预定平面，则获取焊接接头有限元模型以及结构件的有限元模型在所述预定平面投影信息，并根据所述投影信息在所述预定平面上形成计算区域，具体为：
基于焊接接头有限元模型在所述预定平面的第一投影信息，对所述预定平面进行划分，以在所述预定平面上获得母区域；
基于结构件的有限元模型在所述预定平面的第二投影信息，对所述母区域进行划分，以在所述预定平面上形成计算区域；其中，计算区域包括多个子区域。
5.根据权利要求4所述的焊接结构件变形仿真方法，其特征在于，根据在所述预定平面上的焊接塑性应变在所述计算区域的分布，得到结构件对应的塑性应变，以预测结构件的变形情况，具体为；
对每一子区域内的焊接塑性应变进行积分求和，以得到每个子区域的积分点塑性应变值；
根据所述子区域的积分点塑性应变值，对母区域进行积分反求，得到结构件母区域对应的焊缝塑性应变，以预测结构件的变形情况。
6.根据权利要求1所述的焊接结构件变形仿真方法，其特征在于，还包括：
采用Inteweld软件或者开源可视化程序对映射完成后所获得的数据进行可视化处理。
7.一种焊接结构件变形仿真方法装置，其特征在于，包括：
有限元模型获取单元，用于获取结构件的焊接接头有限元模型以及结构件的有限元模型；
焊接残余塑性应变获取单元，用于对所述焊接接头有限元模型进行热弹塑性有限元分析，以获得在焊接接头的局部坐标系下的焊接残余塑性应变；
坐标变换单元，用于将处于局部坐标系下的焊接接头塑性应变通过坐标变换，获得在结构件全局坐标系下的焊接塑性应变；其中，所述焊接塑性应变分布于在所述全局坐标系下的预定平面；所述预定平面为结构件焊缝垂面；
计算区域获取单元，用于获取焊接接头有限元模型以及结构件的有限元模型在所述预定平面的投影信息，并根据所述投影信息在所述预定平面上形成计算区域；
变形情况预测单元，用于根据在所述预定平面上的焊接塑性应变在所述计算区域的分布，得到结构件对应的塑性应变，以预测结构件的变形情况。
8.根据权利要求7所述的焊接结构件变形仿真方法装置，其特征在于，还包括：
三维几何模型获取单元，用于获取预先建立的结构件的三维几何模型；
焊接接头几何模型获取单元，用于根据所述结构件的三维几何模型，获取结构件的焊接接头几何模型；
第一生成单元，用于根据所述焊接接头几何模型进行网络划分，生成焊接接头有限元模型；
第二生成单元，用于根据所述结构件的三维几何模型进行网络划分，生成结构件有限元模型。
9.一种焊接结构件变形仿真方法设备，其特征在于，包括处理器、存储器以及存储在所述存储器内的计算机程序，所述计算机程序能够被所述处理器执行以实现如权利要求1至6任意一项所述的焊接结构件变形仿真方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如权利要求1至6任意一项所述的焊接结构件变形仿真方法。
CN201910918981.0A 2019-09-26 2019-09-26 一种焊接结构件变形仿真方法、装置、设备及存储介质 Active CN110717291B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910918981.0A CN110717291B (zh) 2019-09-26 2019-09-26 一种焊接结构件变形仿真方法、装置、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910918981.0A CN110717291B (zh) 2019-09-26 2019-09-26 一种焊接结构件变形仿真方法、装置、设备及存储介质

Publications (2)

Publication Number Publication Date
CN110717291A true CN110717291A (zh) 2020-01-21
CN110717291B CN110717291B (zh) 2021-05-04

Family

ID=69210968

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910918981.0A Active CN110717291B (zh) 2019-09-26 2019-09-26 一种焊接结构件变形仿真方法、装置、设备及存储介质

Country Status (1)

Country Link
CN (1) CN110717291B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117521452A (zh) * 2023-11-07 2024-02-06 之江实验室 一种基于时空有限元模型的机器人结构件优化方法及装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102152016A (zh) * 2010-02-03 2011-08-17 株式会社日立制作所 焊接变形分析方法
JP2013036902A (ja) * 2011-08-09 2013-02-21 Toshiba Corp 解析装置、評価装置、解析方法および評価方法
KR101745918B1 (ko) * 2016-01-11 2017-06-12 대우조선해양 주식회사 용접 또는 열간가공 시 가공부의 변형량 산출 방법
CN107239591A (zh) * 2017-04-17 2017-10-10 中国核工业华兴建设有限公司 一种核电大型结构焊接变形预测分析方法及其系统
CN107871029A (zh) * 2016-09-26 2018-04-03 首都航天机械公司 预测时效强化铝合金tig焊接头断裂过程的有限元模拟方法
CN108549786A (zh) * 2018-06-21 2018-09-18 中车青岛四方机车车辆股份有限公司 一种大型结构快速焊接数值模拟方法及系统
CN109033496A (zh) * 2018-06-04 2018-12-18 中车青岛四方机车车辆股份有限公司 一种焊接结构的分析方法及装置
CN109543333A (zh) * 2018-12-02 2019-03-29 湖南大学 核电工程用六角管的焊接残余应力与变形的预测方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102152016A (zh) * 2010-02-03 2011-08-17 株式会社日立制作所 焊接变形分析方法
JP2013036902A (ja) * 2011-08-09 2013-02-21 Toshiba Corp 解析装置、評価装置、解析方法および評価方法
KR101745918B1 (ko) * 2016-01-11 2017-06-12 대우조선해양 주식회사 용접 또는 열간가공 시 가공부의 변형량 산출 방법
CN107871029A (zh) * 2016-09-26 2018-04-03 首都航天机械公司 预测时效强化铝合金tig焊接头断裂过程的有限元模拟方法
CN107239591A (zh) * 2017-04-17 2017-10-10 中国核工业华兴建设有限公司 一种核电大型结构焊接变形预测分析方法及其系统
CN109033496A (zh) * 2018-06-04 2018-12-18 中车青岛四方机车车辆股份有限公司 一种焊接结构的分析方法及装置
CN108549786A (zh) * 2018-06-21 2018-09-18 中车青岛四方机车车辆股份有限公司 一种大型结构快速焊接数值模拟方法及系统
CN109543333A (zh) * 2018-12-02 2019-03-29 湖南大学 核电工程用六角管的焊接残余应力与变形的预测方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
L.ZHANG 等: "Evaluation of Applied Plastic Strain Methods for Welding Distortion Prediction", 《JOURNAL OF MANUFACTURING SCIENCE AND ENGINEERING》 *
LVJIE LIANG 等: "A CFD-FEM Model of Residual Stress for Electron Beam Welding Including the Weld Imperfection Effect", 《METALLURGICAL AND MATERIALS TRANSACTIONS A》 *
YOUMIN RONG 等: "Numerical analysis of welding deformation and residual stress in marine propeller nozzle with hybrid laser-arc girth welds", 《INTERNATIONAL JOURNAL OF PRESSURE VESSELS AND PIPING》 *
常勇: "基于局部-整体映射的大型结构件焊接变形与残余应力有限元分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
张忠 等: "基于"局部-全局"映射法的构架装配焊接变形数值仿真", 《国防制造技术》 *
汤一博 等: "夹具拘束下大型薄壁铝合金结构搅拌摩擦焊接变形特性", 《中国机械工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117521452A (zh) * 2023-11-07 2024-02-06 之江实验室 一种基于时空有限元模型的机器人结构件优化方法及装置
CN117521452B (zh) * 2023-11-07 2024-05-07 之江实验室 一种基于时空有限元模型的机器人结构件优化方法及装置

Also Published As

Publication number Publication date
CN110717291B (zh) 2021-05-04

Similar Documents

Publication Publication Date Title
JP5830004B2 (ja) ３次元モデル生成装置、３次元モデル生成方法及び３次元モデル生成プログラム
EP2887241A1 (en) Cad-based initial surface geometry correction
CN112070908B (zh) 三维模型的自动拆分方法、计算机设备以及存储介质
CN112733318A (zh) 一种自适应网格细分方法、装置、设备及存储介质
CN110717291B (zh) 一种焊接结构件变形仿真方法、装置、设备及存储介质
CN110688798B (zh) 一种壳体结构件的变形预测方法、装置、设备及存储介质
CN114861500A (zh) 基于三维点云自动生成隧道结构有限元模型的方法及系统
CN111291450A (zh) 一种确定发动机罩回弹三维模型的方法、装置及设备
CN108744520B (zh) 确定游戏模型摆放位置的方法、装置和电子设备
CN113378435A (zh) 粒子生成方法、装置、设备及存储介质
JP2003263464A (ja) フィレット作成方法、及び３次元ｃａｄプログラム
CN114998381A (zh) 管板焊接中焊接轨迹拟合方法、装置、设备及存储介质
CN114329868A (zh) 变电设备安全距离校验方法、装置和计算机设备
JP2002207777A (ja) 中立面モデルの生成方法
CN109840945B (zh) 有限元预处理方法及装置
CN112329165B (zh) 一种轮罩骨架方管弯曲回弹的建模方法、装置及设备
KR102156336B1 (ko) 3차원 형상 재현 방법
CN116011049B (zh) 一种结构化网格过渡拓扑结构的参数化生成方法及装置
CN112395658B (zh) 模型细节判断方法、终端设备及计算机可读存储介质
US20230196672A1 (en) 3d model generation
US20150106064A1 (en) Computing device and method for processing point clouds
CN116861023B (zh) 一种在三维空间中确定相同几何对象的方法及系统
CN111652957B (zh) 一种线性工程设计图生成方法
CN110335353B (zh) 街区自动细分的方法、装置、设备及存储介质
CN112685871B (zh) 防雷保护范围模拟方法、装置、设备和存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant