CN113849923A - 一种基于abaqus的钢结构整体焊接变形的数值模拟方法 - Google Patents

一种基于abaqus的钢结构整体焊接变形的数值模拟方法 Download PDF

Info

Publication number
CN113849923A
CN113849923A CN202110952622.4A CN202110952622A CN113849923A CN 113849923 A CN113849923 A CN 113849923A CN 202110952622 A CN202110952622 A CN 202110952622A CN 113849923 A CN113849923 A CN 113849923A
Authority
CN
China
Prior art keywords
welding
steel structure
node
parameters
determining
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
CN202110952622.4A
Other languages
English (en)
Other versions
CN113849923B (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.)
Beijing Construction Industry Refco Group Ltd
Beijing Machinery Construction Group Co ltd
Original Assignee
Beijing Machinery Construction Group Co ltd
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 Beijing Machinery Construction Group Co ltd filed Critical Beijing Machinery Construction Group Co ltd
Priority to CN202110952622.4A priority Critical patent/CN113849923B/zh
Publication of CN113849923A publication Critical patent/CN113849923A/zh
Application granted granted Critical
Publication of CN113849923B publication Critical patent/CN113849923B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,包括:利用ABAQUS软件,建立钢结构的焊接结构模型,确定所述焊接结构模型中的焊接节点和焊接类型;设置焊接工艺参数,并按照所述焊接工艺参数完成模拟焊接;完成对所述钢结构的整体焊接变形的数值计算,保证了对钢结构整体焊接变形模拟计算的检测精度,并提高了检测效率,根据焊接变形模拟计算的检测结果,保证了钢结构的稳定性。

Description

一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法
技术领域
本发明涉及建筑施工技术领域,特别涉及一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法。
背景技术
目前,随着科学技术的不断发展,钢结构造型和结构形式变得也越来越复杂,同时必然会给钢结构工程带来一系列挑战,在钢结构焊接时,由于局部高温加热而造成焊件上温度分布不均匀,最终导致在结构内部产生了焊接应力与变形,以及钢结构的整体焊接变形,降低其承载能力和机械性能。
现有的钢结构整体焊接变形的检测方法大多为人工检测,存在检测误差大、精度低,检测效率低的问题,无法保证后续钢结构的稳定性。
因此,本发明提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法。
发明内容
本发明提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,基于ABAQUS软件,保证了对钢结构整体焊接变形模拟计算的检测精度,并提高了检测效率,根据焊接变形模拟计算的检测结果,保证了钢结构的稳定性。
本发明提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,包括:
步骤1:利用ABAQUS软件,建立钢结构的焊接结构模型;
步骤2,对所述焊接结构模型进行分析,确定所述焊接结构模型中的焊接节点和焊接类型;
步骤3:根据所述焊接节点和焊接类型,设置焊接工艺参数,并按照所述焊接工艺参数完成模拟焊接;
步骤4:根据模拟焊接结果,完成对所述钢结构的整体焊接变形的数值计算。
在一种可能实现的方式中,
步骤1中,利用ABAQUS软件,建立钢结构的焊接结构模型包括:
步骤101:获取所述钢结构的材质特征和部件特征,并基于所述ABAQUS软件,建立模板库;
步骤102:对所述钢结构进行三维激光扫描,获取点云数据,并对所述点云数据进行遍历,删除与所述钢结构无关的数据,得到目标点云数据;
步骤103:根据所述目标点云数据构建钢结构框架图;
步骤104:基于所述模板库,对所述钢结构框架图进行填充,得到焊接结构模型。
在一种可能实现的方式中,
步骤2中,对所述焊接结构模型进行分析,确定所述焊接结构模型中的焊接节点和焊接类型包括:
基于边折叠算法对所述焊接结构模型进行特征约束和几何形状约束,得到所述焊接模型的简化模型;
对所述简化模型进行表面分解得到多个曲面,并对所述多个曲面进行网格划分,生成多个网格;
利用3D技术将所有网格进行排序预览,并根据预览结果将所有网格加载在一起,生成六面体网格;
对所述六面体网格进行分析,得到焊接节点;
对焊接节点进行分析,确定所述焊接节点对应的焊接类型。
在一种可能实现的方式中,
对所述六面体网格进行分析,得到焊接节点包括:
获取所述六面体网格中每个网格单元的第一图像,和相邻的两个网格单元之间连接节点的第二图像;
对所述第一图像进行高斯滤波处理后,确定所述第一图像中像素值在预设范围内的目标区域,并利用纹理直方图计算方法获取所述目标区域的纹理特征,并计算所述纹理特征与预设纹理特征之间的相似性;
选取相似性大于预设相似值的目标区域作为第一目标区域;
获取所述第一目标区域的基本信息,并基于所述基本信息确定第一目标区域对应的钢结构表面的阻抗数据,且根据所述阻抗数据确定第一目标区域的光滑度;
选择所述光滑度大于预设光滑度的第一目标区域,作为第二目标区域;
判断所述第二目标区域是否存在所述连接节点的部分区域;
若是,对所述第二图像进行检测,确定所述第二图像中所包含的部分区域,并对所述部分区域的边界进行拓展,得到拓展区域,选择所述拓展区域的灰度值和像素值在预设阈值范围内的目标子区域;
基于所述目标子区域,对所述第二目标区域进行修正,并将修正后的第二目标区域作为焊接区域;
否则,将所述第二目标区域作为焊接区域;
基于所述六面体网格与所述焊接结构模型的映射关系和缩放比例,确定所述焊接区域在所述焊接结构模型中的第一位置;
判断所述第一位置是否在所述焊接结构模型的预设焊接位置范围内;
若是,对所述第一位置进行标记,得到焊接节点;
否则,基于误差范围对所述第一位置进行误差补偿,得到第二位置,并对所述第二位置进行标记,得到焊接节点。
在一种可能实现的方式中,
对所述焊接节点进行分析,确定所述焊接节点对应的焊接类型包括:
建立所述焊接节点的坐标系,确定每个焊接节点在所述坐标系中的坐标位置,并根据所述坐标位置对所述焊接节点进行分组;
根据每个分组中焊接节点的坡口类型确定所述焊接节点的焊接类型。
在一种可能实现的方式中,
步骤3中,根据所述焊接节点和焊接类型,设置焊接工艺参数包括:
获取所述焊接节点和焊接类型在历史焊接中的历史焊接参数;
对所述历史焊接参数进行预分析,确定所述历史焊接参数基于焊接结果的评定数据,并根据所述评定数据对所述历史焊接参数进行等级划分,得到所述历史焊接参数的等级;
对等级小于预设等级的历史焊接参数设置第一权重值,对等级大于等于所述预设等级的历史焊接参数设置第二权重值,得到初始焊接参数;
基于所述初始焊接参数,为所述焊接节点设置初始焊接工艺参数;
遍历每个焊接节点,确定每个焊接节点的节点属性,并根据所述节点属性,筛选出第一焊接节点,其余作为第二焊接节点;
根据所述第一焊接节点的节点属性,为所述第一焊接节点确定第一焊接顺序,根据所述第二焊接节点的节点属性,在所述第一焊接顺序的基础上,为所述第二焊接节点设置第二焊接顺序;
基于所述第一焊接顺序和第二焊接顺序,得到焊接工艺顺序;
基于所述焊接节点的节点属性,确定每个焊接节点的焊接约束条件;
判断所述焊接节点的初始焊接工艺参数和焊接工艺顺序是否满足所述约束条件;
若是,则将所述初始焊接工艺参数和焊接工艺顺序作为焊接工艺参数;
否则,基于所述焊接约束条件对所述初始焊接工艺参数和到焊接工艺顺序进行优化,得到焊接工艺参数。
在一种可能实现的方式中,
步骤3中,按照所述焊接工艺参数完成模拟焊接包括:
根据所述焊接工艺参数确定模拟焊接的材料参数和热边界条件,并将所述材料参数和热边界条件载入所述ABAQUS软件中,构建焊接条件数据;
基于所述焊接条件数据对所述焊接结构模型进行模拟焊接。
在一种可能实现的方式中,
在步骤4之前,还包括,获取模拟结果,其过程如下:
对所述模拟焊接过程的温度进行监测,并根据所述焊接结构模型中的材料参数建立所述焊接结构模型对应的热源模型,并基于所述热源模型确定在所述焊接工艺参数下的温度场;
获取所述温度场中各个焊接节点的温度值,并将所述温度值,结合材料参数转化为载荷加入热-结构转化单元模型中,得到所述热源模型确定在所述焊接工艺参数下的应力场;
所述温度场和应力场作为所述模拟结果。
在一种可能实现的方式中,
步骤4中,根据模拟焊接结果,完成对所述钢结构的整体焊接变形的数值计算包括:
步骤401:根据所述模拟结果的温度场确定所述钢结构在焊接过程中的热流密度;
步骤402:根据所述模拟结果的应力场确定所述钢结构在焊接过程中的残余应力;
步骤403,基于所述钢结构在焊接过程中的热流密度和残余应力,确定所述钢结构的整体变形系数。
在一种可能实现的方式中,
获取所述钢结构的整体变形系数后还包括:
判断所述整体变形系数是否小于预设变形系数;
若是,按照设置的焊接工艺参数对所述钢结构进行实际焊接;
否则,基于所述整体变形系数,对所述焊接工艺参数进行调整。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1为本发明实施例中一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法的流程图;
图2为本发明实施例中所述焊接结构模型建立的流程图;
图3为本发明实施例中对所述钢结构的整体焊接变形的数值计算的及流程图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
实施例1
本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,如图1所示,包括:
步骤1:利用ABAQUS软件,建立钢结构的焊接结构模型;
步骤2:对所述焊接结构模型进行分析,确定所述焊接结构模型中的焊接节点和焊接类型;
步骤3:根据所述焊接节点和焊接类型,设置焊接工艺参数,并按照所述焊接工艺参数完成模拟焊接;
步骤4:根据模拟焊接结果,完成所述钢结构的整体焊接变形的数值计算。
在该实施例中,所述ABAQUS软件是功能强大的有限元软件,可以对应力、位置等进行静态和动态分析。
在该实施例中,所述焊接结构模型为所述钢结构在ABAQUS软件上的三维仿真模型,用来反应所述钢结构的外形以及材质等特征。
在该实施例中,所述焊接工艺参数包括:焊接顺序、焊接电流电压、焊接速度、多层多道焊的数量、环境温度、预热和保温措施、降温速度等参数的设定。
在该实施例中,所述焊接类型包括熔焊、压焊、钎焊等。
上述设计方案的有益效果是:通过ABAQUS软件,建立钢结构的焊接结构模型和焊接模拟过程,实现对钢结构的整体焊接变形的数值计算,基于ABAQUS软件,保证焊接结构模型和模拟焊接过程的准确性,从而保证了对钢结构整体焊接变形模拟计算的检测精度,并提高了检测效率,根据焊接变形模拟计算的检测结果,方便及时进行对焊接工艺参数进行调整,保证了钢结构的稳定性。
实施例2
基于实施例1的基础上,本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,如图2所示,步骤1中,利用ABAQUS软件,建立钢结构的焊接结构模型包括:
步骤101:获取所述钢结构的材质特征和部件特征,并基于所述ABAQUS软件,建立模板库;
步骤102:对所述钢结构进行三维激光扫描,获取点云数据,并对所述点云数据进行遍历,删除与所述钢结构无关的数据,得到目标点云数据;
步骤103:根据所述目标点云数据构建钢结构框架图;
步骤104:基于所述模板库,对所述钢结构框架图进行填充,得到焊接结构模型。
在该实施例中,建立基于ABAQUS软件,建立模板库后,可以方便相似的钢结构进行重复使用,提高利用率和类似工程的构建效率。
在该实施例中,删除与所述钢结构无关的数据为在三维激光扫描中因为光线、障碍物等的干扰造成的无关数据,删除与所述钢结构无关的数据可以提高对钢结构框架图构建的准确性,更好的反应钢结构的外形。
在该实施例中,对所述钢结构框架图进行填充具体为根据所述钢结构的材质特征和部件特征为所述钢结构框架图设置钢结构材质参数,并根据结构材质参数用不同的颜色进行填充,来表示不同的材质参数,例如硬度、张力等。
上述设计方案的有益效果是:通过根据钢结构的三维点云数据、材质特征和部件特征来构建焊接结构模型,保证了钢结构的外形以及所述钢结构的材质参数的准确性,使得焊接结构模型更好地反应钢结构的实际情况,为保证对钢结构整体焊接变形模拟计算的检测精度提供了基础。
实施例3
基于实施例1的基础上,本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,步骤2中,对所述焊接结构模型进行分析,确定所述焊接结构模型中的焊接节点和焊接类型包括:
基于边折叠算法对所述焊接结构模型进行特征约束和几何形状约束,得到所述焊接模型的简化模型;
对所述简化模型进行表面分解得到多个曲面,并对所述多个曲面进行网格划分,生成多个网格;
利用3D技术将所有网格进行排序预览,并根据预览结果将所有网格加载在一起,生成六面体网格;
对所述六面体网格进行节点分析,得到焊接节点;
对在所述六面体网格中的焊接节点进行分析,确定所述焊接节点对应的焊接类型。
在该实施例中,所述边折叠算法属于几何元素删除法的一种,每次简化时,通过算法选定一条有向边e以及相关的2个点(u,v),将其中一个点u“折叠”至v,然后修改拓扑关系,将与u相关的边映射到v,最后完成简化操作。一次简化可以减少所述焊接结构模型的1条边和2个面。
在该实施例中,所述特征约束为通过边折叠算法对所述焊接结构模型进行简化,生成连续的细节层次,保持对应的纹理不变,更好地反应所述焊接结构模型的特征,所述几何约束为通过边折叠算法对所述焊接结构模型进行简化,保持对应的顶点不变,更好地反应所述焊接结构模型的几何形状。
在该实施例中,所述特征识别结果为对应曲面的参数特征,例如曲面的凹凸程度、表面光滑度、曲面面积大小、曲面形状。
在该实施例中,对所述多个曲面进行网格划分的划分尺寸与所述曲面的凹凸程度、表面光滑度、曲面面积大小、曲面形状有关,所述曲面的凹凸程度越大、表面光滑度越小、曲面面积越大、曲面形状越不规则,对应的划分尺寸越精细,提高了对所述焊接结构模型的分析精度,从而提高了焊接节点和焊接类型的确定精度。
上述设计方案的有益效果是:通过对焊接结构模型进行分析,生成六面体网格,从而得到焊接节点和焊接类型,保证了获取焊接节点和焊接类型的精度。
实施例4
基于实施例3的基础上,本发明实施例提供一种一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,对所述六面体网格进行分析,得到焊接节点包括:
获取所述六面体网格中每个网格单元的第一图像,和相邻的两个网格单元之间连接节点的第二图像;
对所述第一图像进行高斯滤波处理后,确定所述第一图像中像素值在预设范围内的目标区域,并利用纹理直方图计算方法获取所述目标区域的纹理特征,并计算所述纹理特征与预设纹理特征之间的相似性;
选取相似性大于预设相似值的目标区域作为第一目标区域;
获取所述第一目标区域的基本信息,并基于所述基本信息确定第一目标区域对应的钢结构表面的阻抗数据,且根据所述阻抗数据确定第一目标区域的光滑度;
选择所述光滑度大于预设光滑度的第一目标区域,作为第二目标区域;
判断所述第二目标区域是否存在所述连接节点的部分区域;
若是,对所述第二图像进行检测,确定所述第二图像中所包含的部分区域,并对所述部分区域的边界进行拓展,得到拓展区域,选择所述拓展区域的灰度值和像素值在预设阈值范围内的目标子区域;
基于所述目标子区域,对所述第二目标区域进行修正,并将修正后的第二目标区域作为焊接区域;
否则,将所述第二目标区域作为焊接区域;
基于所述六面体网格与所述焊接结构模型的映射关系和缩放比例,确定所述焊接区域在所述焊接结构模型中的第一位置;
判断所述第一位置是否在所述焊接结构模型的预设焊接位置范围内;
若是,对所述第一位置进行标记,得到焊接节点;
否则,基于误差范围对所述第一位置进行误差补偿,得到第二位置,并对所述第二位置进行标记,得到焊接节点。
在该实施例中,所述基本信息包括所述第一目标区域的灰度值、像素值信息,不同的灰度值、亮度值取值对应不同的阻抗数据,所述阻抗数据越小,对应的光滑度越大。
在该实施例中,焊接节点的表面应该光滑、洁净,有利于焊接。
在该实施例中,判断所述第二目标区域是否存在所述连接节点的部分区域,若存在,则通过对所述第二图像进行检测来对所述第二目标区域进行修正,可以避免因为对相邻的两个网格单元之间的连接节点缺少检测,而造成对最终确定的焊接区域的误差,通过对相邻的两个网格单元之间的连接节点的第二图像进行检测,可以对焊接区域在所述第二图像中进行仔细检测,保证焊接区域的准确性。
在该实施例中,所述预设纹理特征为预先确定的标准焊接节点的纹理特征。
在该实施例中,所述焊接结构模型的预设焊接位置范围根据钢结构特点确定的大致焊接位置范围,通过进一步检测分析,确定在所述大致焊接位置范围中准确的焊接位置。
在该实施例中,基于误差范围对所述第一位置进行误差补偿具体为根据误差范围将所述第一位置进行整体移动,所述误差补偿为移动的位置量。
上述设计方案的有益效果是:通过对六面体网格的每个网格单元集相邻网格单元的节点进行检测分析,保证了确定的焊接节点的精确性,从而保证了模拟焊接和数值模拟计算的准确性。
实施例5
基于实施例3的基础上,本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,对所述焊接节点进行分析,确定所述焊接节点对应的焊接类型包括:
建立所述焊接节点的坐标系,确定每个焊接节点在所述坐标系中的坐标位置,并根据所述坐标位置对所述焊接节点进行分组;
根据每个分组中焊接节点的坡口类型确定所述焊接节点的焊接类型。
在该实施例中,根据所述坐标位置对所述焊接节点进行分组具体为横坐标或纵坐标相同的焊接节点为一组。
在该实施例中,所述焊接节点的坡口类型包括I形坡口、V形坡口、X形坡口、U形坡口。
在该实施例中,根据所述坐标位置对所述焊接节点进行分组,方便了对焊接节点的焊接类型的判断,提高了效率。
上述设计方案的有益效果是:通过根据坐标位置对所述焊接节点进行分组,以分组结果为所述焊接节点确定焊接类型,提高了对焊接类型确定的效率,通过根据焊接节点的坡口类型确定焊接类型,为模拟焊接提供基础。
实施例6
基于实施例1的基础上,本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,步骤3中,根据所述焊接节点和焊接类型,设置焊接工艺参数包括:
获取所述焊接节点和焊接类型在历史焊接中的历史焊接参数;
对所述历史焊接参数进行预分析,确定所述历史焊接参数基于焊接结果的评定数据,并根据所述评定数据对所述历史焊接参数进行等级划分,得到所述历史焊接参数的等级;
对等级小于预设等级的历史焊接参数设置第一权重值,对等级大于等于所述预设等级的历史焊接参数设置第二权重值,得到初始焊接参数;
基于所述初始焊接参数,为所述焊接节点设置初始焊接工艺参数;
遍历每个焊接节点,确定每个焊接节点的节点属性,并根据所述节点属性,筛选出第一焊接节点,其余作为第二焊接节点;
根据所述第一焊接节点的节点属性,为所述第一焊接节点确定第一焊接顺序,根据所述第二焊接节点的节点属性,在所述第一焊接顺序的基础上,为所述第二焊接节点设置第二焊接顺序;
基于所述第一焊接顺序和第二焊接顺序,得到焊接工艺顺序;
基于所述焊接节点的节点属性,确定每个焊接节点的焊接约束条件;
判断所述焊接节点的初始焊接工艺参数和焊接工艺顺序是否满足所述约束条件;
若是,则将所述初始焊接工艺参数和焊接工艺顺序作为焊接工艺参数;
否则,基于所述焊接约束条件对所述初始焊接工艺参数和到焊接工艺顺序进行优化,得到焊接工艺参数。
在该实施例中,所述设置初始焊接工艺参数具体是是为所述焊接节点匹配相似度最高的基于所述初始焊接参数中的历史焊接节点,将所述历史焊接节点对应的初始焊接参数作为所述焊接节点的初始焊接工艺参数。
在该实施例中,所述评定数据为以所述历史焊接参数进行焊接后对结构的焊接残余应力、局部变形和整体变形的评定,所述焊接残余应力、局部变形和整体变形均越小,对应的等级越高。
在该实施例中,所述第一权重值大于所述第二权重值,等级越低,需要的调整越大,故对应的第一权重越大。利用第一权重值和第二权重值对所述历史焊接参数进行调整,可以使焊接参数加强对结构的稳定性。
在该实施例中,所述焊接节点的节点属性包括节点位置、焊接类型、焊点厚度。
在该实施例中,所述第一焊接节点为主焊点,第二焊接点为次焊点,于焊点的节点位置、焊接类型相关,第一焊接顺序为所述第一焊接节点的焊接顺序,在确定焊接顺序是根据所述节点位置和焊接类型确定,例如节点位置在边上的焊接节点比节点位置在中间的焊接顺序靠前,焊接类型难度较大的比焊接类型难度较小的焊接顺序靠前。
在该实施例中,所述焊接约束条件包括对每个焊接节点的焊接工艺参数例如焊接电流电压、焊接速度、多层多道焊的数量、环境温度等的范围设定,以及对焊接顺序的约束,例如第三个焊接节点的焊接顺序必须在第二个焊接节点的前面。
上述设计方案的有益效果是:通过根据历史焊接参数和每个焊接节点的节点属性,为所述焊接节点设置更优的焊接工艺参数,减小对钢结构的变形,保证了钢结构的稳定性,为模拟焊接提供优秀的工艺参数。
实施例7
基于实施例1的基础上,本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,步骤3中,按照所述焊接工艺参数完成模拟焊接包括:
根据所述焊接工艺参数确定模拟焊接的材料参数和热边界条件,并将所述材料参数和热边界条件载入所述ABAQUS软件中,构建焊接条件数据;
基于所述焊接条件数据对所述焊接结构模型进行模拟焊接。
在该实施例中,所述材料参数包括焊料类型、焊料用量等。
在该实施例中,所述热边界条件包括焊接环境温度场、焊接温度等。
上述设计方案的有益效果是:根据焊接工艺参数确定材料参数和热边界条件来完成模拟焊接,基于ABAQUS软件提高模拟焊接的焊接精度和焊接效率。
实施例8
基于实施例1的基础上,本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,在步骤4之前,还包括,获取模拟结果,其过程如下:
对所述模拟焊接过程的温度进行监测,并根据所述焊接结构模型中的材料参数建立所述焊接结构模型对应的热源模型,并基于所述热源模型确定在所述焊接工艺参数下的温度场;
获取所述温度场中各个焊接节点的温度值,并将所述温度值,结合材料参数转化为载荷加入热-结构转化单元模型中,得到所述热源模型确定在所述焊接工艺参数下的应力场;
所述温度场和应力场作为所述模拟结果。
上述设计方案的有益效果是:通过获取所述钢结构焊接过程中的温度场和应力场,为数值模拟计算提供基础。
实施例9
基于实施例1的基础上,本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,如图3所示,步骤4中,根据模拟焊接结果,完成对所述钢结构的整体焊接变形的数值计算包括:
步骤401:根据所述模拟结果的温度场确定所述钢结构在焊接过程中的热流密度,其计算公式如下:
Figure BDA0003219047170000151
其中,δ表示所述热流密度,V表示所述温度场的体积,U表示焊接电压,I表示焊接电流,e表示自然常数取值为2.72,a表示所述温度场的最小长度,b表示所述温度场的最大长度;
步骤402:根据所述模拟结果的应力场确定所述钢结构在焊接过程中的残余应力;
步骤403:基于所述钢结构在焊接过程中的热流密度和残余应力,确定所述钢结构的整体变形系数,其计算公式如下:
Figure BDA0003219047170000152
其中,ω表示所述钢结构的整体变形系数,R表示所述钢结构的表面传热系数,τ表示所述钢结构在焊接过程中的残余应力,m表示所述钢结构的质量,g表示重力加速度,取值为9.8N/Kg,ΔT表示所述钢结构在焊接过程中的温差变化。
在该实施例中,所述温度场的最小长度与最大长度的平方的比值越小,表明所述温度场的可流动性越好,对应的热流密度就越小。
在该实施例中,所述热流密度表示钢结构在单位体积内的热功率,单位为w/m3,w表示瓦,m3表示立方米。
在该实施例中,对于公式
Figure BDA0003219047170000153
例如可以是:U=200v,V=50m3,I=20A,a=4m,b=5m,则对应的δ大约为84w/m3
在该实施例中,对于公式
Figure BDA0003219047170000161
例如可以是:δ=84w/m3,R=1000w/m·℃,τ=50MPa=50000N/m2,m=485kg,mg约等于4850N,ΔT=5℃,得到ω大约等于2.24。
上述设计方案的有益效果是:基于ABAQUS,根据模拟焊接结果,完成对所述钢结构的整体焊接变形的数值计算,保证了对钢结构整体焊接变形模拟计算的检测精度,并提高了检测效率。
实施例10
基于实施例9的基础上,本发明实施例提供一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,获取所述钢结构的整体变形系数后还包括:
判断所述整体变形系数是否小于预设变形系数;
若是,按照设置的焊接工艺参数对所述钢结构进行实际焊接;
否则,基于所述整体变形系数,对所述焊接工艺参数进行调整。
在该实施例中,所述预设变形系数的取值为2.5,从实施例9得到整体变形系数为2.24,则说明所述钢结构的整体变形在可接受范围内。
上述设计方案的有益效果是:通过比较整体变形系数是否小于预设变形系数,判断钢结构的整体变形是否在可接受范围内,若没有,则重新设置接工艺参数来减小钢结构的整体变形系数,保证了钢结构的稳定性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (10)

1.一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,包括:
步骤1:利用ABAQUS软件,建立钢结构的焊接结构模型;
步骤2:对所述焊接结构模型进行分析,确定所述焊接结构模型中的焊接节点和焊接类型;
步骤3:根据所述焊接节点和焊接类型,设置焊接工艺参数,并按照所述焊接工艺参数完成模拟焊接;
步骤4:根据模拟焊接结果,完成对所述钢结构的整体焊接变形的数值计算。
2.根据权利要求1所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,步骤1中,利用ABAQUS软件,建立钢结构的焊接结构模型包括:
步骤101:获取所述钢结构的材质特征和部件特征,并基于所述ABAQUS软件,建立模板库;
步骤102:对所述钢结构进行三维激光扫描,获取点云数据,并对所述点云数据进行遍历,删除与所述钢结构无关的数据,得到目标点云数据;
步骤103:根据所述目标点云数据构建钢结构框架图;
步骤104:基于所述模板库,对所述钢结构框架图进行填充,得到焊接结构模型。
3.根据权利要求1所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,步骤2中,对所述焊接结构模型进行分析,确定所述焊接结构模型中的焊接节点和焊接类型包括:
基于边折叠算法对所述焊接结构模型进行特征约束和几何形状约束,得到所述焊接模型的简化模型;
对所述简化模型进行表面分解得到多个曲面,并对所述多个曲面进行网格划分,生成多个网格;
利用3D技术将所有网格进行排序预览,并根据预览结果将所有网格加载在一起,生成六面体网格;
对所述六面体网格进行分析,得到焊接节点;
对焊接节点进行分析,确定所述焊接节点对应的焊接类型。
4.根据权利要求3所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,对所述六面体网格进行分析,得到焊接节点包括:
获取所述六面体网格中每个网格单元的第一图像,和相邻的两个网格单元之间连接节点的第二图像;
对所述第一图像进行高斯滤波处理后,确定所述第一图像中像素值在预设范围内的目标区域,并利用纹理直方图计算方法获取所述目标区域的纹理特征,并计算所述纹理特征与预设纹理特征之间的相似性;
选取相似性大于预设相似值的目标区域作为第一目标区域;
获取所述第一目标区域的基本信息,并基于所述基本信息确定第一目标区域对应的钢结构表面的阻抗数据,且根据所述阻抗数据确定第一目标区域的光滑度;
选择所述光滑度大于预设光滑度的第一目标区域,作为第二目标区域;
判断所述第二目标区域是否存在所述连接节点的部分区域;
若是,对所述第二图像进行检测,确定所述第二图像中所包含的部分区域,并对所述部分区域的边界进行拓展,得到拓展区域,选择所述拓展区域的灰度值和像素值在预设阈值范围内的目标子区域;
基于所述目标子区域,对所述第二目标区域进行修正,并将修正后的第二目标区域作为焊接区域;
否则,将所述第二目标区域作为焊接区域;
基于所述六面体网格与所述焊接结构模型的映射关系和缩放比例,确定所述焊接区域在所述焊接结构模型中的第一位置;
判断所述第一位置是否在所述焊接结构模型的预设焊接位置范围内;
若是,对所述第一位置进行标记,得到焊接节点;
否则,基于误差范围对所述第一位置进行误差补偿,得到第二位置,并对所述第二位置进行标记,得到焊接节点。
5.根据权利要求3所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,对所述焊接节点进行分析,确定所述焊接节点对应的焊接类型包括:
建立所述焊接节点的坐标系,确定每个焊接节点在所述坐标系中的坐标位置,并根据所述坐标位置对所述焊接节点进行分组;
根据每个分组中焊接节点的坡口类型确定所述焊接节点的焊接类型。
6.根据权利要求1所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,步骤3中,根据所述焊接节点和焊接类型,设置焊接工艺参数包括:
获取所述焊接节点和焊接类型在历史焊接中的历史焊接参数;
对所述历史焊接参数进行预分析,确定所述历史焊接参数基于焊接结果的评定数据,并根据所述评定数据对所述历史焊接参数进行等级划分,得到所述历史焊接参数的等级;
对等级小于预设等级的历史焊接参数设置第一权重值,对等级大于等于所述预设等级的历史焊接参数设置第二权重值,得到初始焊接参数;
基于所述初始焊接参数,为所述焊接节点设置初始焊接工艺参数;
遍历每个焊接节点,确定每个焊接节点的节点属性,并根据所述节点属性,筛选出第一焊接节点,其余作为第二焊接节点;
根据所述第一焊接节点的节点属性,为所述第一焊接节点确定第一焊接顺序,根据所述第二焊接节点的节点属性,在所述第一焊接顺序的基础上,为所述第二焊接节点设置第二焊接顺序;
基于所述第一焊接顺序和第二焊接顺序,得到焊接工艺顺序;
基于所述焊接节点的节点属性,确定每个焊接节点的焊接约束条件;
判断所述焊接节点的初始焊接工艺参数和焊接工艺顺序是否满足所述约束条件;
若是,则将所述初始焊接工艺参数和焊接工艺顺序作为焊接工艺参数;
否则,基于所述焊接约束条件对所述初始焊接工艺参数和到焊接工艺顺序进行优化,得到焊接工艺参数。
7.根据权利要求1所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,步骤3中,按照所述焊接工艺参数完成模拟焊接包括:
根据所述焊接工艺参数确定模拟焊接的材料参数和热边界条件,并将所述材料参数和热边界条件载入所述ABAQUS软件中,构建焊接条件数据;
基于所述焊接条件数据对所述焊接结构模型进行模拟焊接。
8.根据权利要求1所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,在步骤4之前,还包括,获取模拟结果,其过程如下:
对所述模拟焊接过程的温度进行监测,并根据所述焊接结构模型中的材料参数建立所述焊接结构模型对应的热源模型,并基于所述热源模型确定在所述焊接工艺参数下的温度场;
获取所述温度场中各个焊接节点的温度值,并将所述温度值,结合材料参数转化为载荷加入热-结构转化单元模型中,得到所述热源模型确定在所述焊接工艺参数下的应力场;
所述温度场和应力场作为所述模拟结果。
9.根据权利要求1所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,步骤4中,根据模拟焊接结果,完成对所述钢结构的整体焊接变形的数值计算包括:
步骤401:根据所述模拟结果的温度场确定所述钢结构在焊接过程中的热流密度;
步骤402:根据所述模拟结果的应力场确定所述钢结构在焊接过程中的残余应力;
步骤403:基于所述钢结构在焊接过程中的热流密度和残余应力,确定所述钢结构的整体变形系数。
10.根据权利要求9所述的一种基于ABAQUS的钢结构整体焊接变形的数值模拟方法,其特征在于,获取所述钢结构的整体变形系数后还包括:
判断所述整体变形系数是否小于预设变形系数;
若是,按照设置的焊接工艺参数对所述钢结构进行实际焊接;
否则,基于所述整体变形系数,对所述焊接工艺参数进行调整。
CN202110952622.4A 2021-08-19 2021-08-19 一种基于abaqus的钢结构整体焊接变形的数值模拟方法 Active CN113849923B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110952622.4A CN113849923B (zh) 2021-08-19 2021-08-19 一种基于abaqus的钢结构整体焊接变形的数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110952622.4A CN113849923B (zh) 2021-08-19 2021-08-19 一种基于abaqus的钢结构整体焊接变形的数值模拟方法

Publications (2)

Publication Number Publication Date
CN113849923A true CN113849923A (zh) 2021-12-28
CN113849923B CN113849923B (zh) 2022-05-24

Family

ID=78976025

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110952622.4A Active CN113849923B (zh) 2021-08-19 2021-08-19 一种基于abaqus的钢结构整体焊接变形的数值模拟方法

Country Status (1)

Country Link
CN (1) CN113849923B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103246772A (zh) * 2013-05-11 2013-08-14 天津大学 一种基于abaqus的超声喷丸矫正焊接变形的有限元模拟方法
CN106326574A (zh) * 2016-08-27 2017-01-11 哈尔滨工业大学(威海) 一种计及焊接变形和残余应力的极限强度数值计算方法
CN110227876A (zh) * 2019-07-15 2019-09-13 西华大学 基于3d点云数据的机器人焊接路径自主规划方法
CN111733969A (zh) * 2020-06-30 2020-10-02 浙江精工钢结构集团有限公司 一种焊接球节点外包锥管加强构造及其制作方法
CN113065272A (zh) * 2021-04-30 2021-07-02 中国核工业华兴建设有限公司 一种焊接变形数值模拟快速实现方法
CN113158380A (zh) * 2021-05-26 2021-07-23 石河子大学 集箱管座多层多道焊焊接过程数值模拟计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103246772A (zh) * 2013-05-11 2013-08-14 天津大学 一种基于abaqus的超声喷丸矫正焊接变形的有限元模拟方法
CN106326574A (zh) * 2016-08-27 2017-01-11 哈尔滨工业大学(威海) 一种计及焊接变形和残余应力的极限强度数值计算方法
CN110227876A (zh) * 2019-07-15 2019-09-13 西华大学 基于3d点云数据的机器人焊接路径自主规划方法
CN111733969A (zh) * 2020-06-30 2020-10-02 浙江精工钢结构集团有限公司 一种焊接球节点外包锥管加强构造及其制作方法
CN113065272A (zh) * 2021-04-30 2021-07-02 中国核工业华兴建设有限公司 一种焊接变形数值模拟快速实现方法
CN113158380A (zh) * 2021-05-26 2021-07-23 石河子大学 集箱管座多层多道焊焊接过程数值模拟计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
王能庆: "616装甲钢焊接残余应力与变形有限元分析", 《中国优秀硕士学位论文全文数据库工程科技I辑》 *
贾坤荣等: "厚板焊接温度场和残余应力场的数值模拟", 《热加工工艺》 *
赵炳桢: "《现代刀具设计与应用》", 30 September 2014 *

Also Published As

Publication number Publication date
CN113849923B (zh) 2022-05-24

Similar Documents

Publication Publication Date Title
CN112581457A (zh) 一种基于三维点云的管道内表面检测方法及装置
CN107767457A (zh) 一种基于点云快速重建的stl数模生成方法
JP2013507679A (ja) 三次元物体モデルの3dプリントが可能な方法及びシステム
CN109872394B (zh) 基于最小二乘支持向量机的狭长三角形网格优化方法
CN110060342B (zh) 一种三维曲面拟合方法
US6573892B1 (en) Method of constructing surface element layers of hexahedral mesh for finite element analysis and method of constructing hexahedral mesh using the same
CN112258521B (zh) 基于单张椭圆成像射线照片的管道焊缝三维重建方法
AU769960B2 (en) Mesh generator for and method of generating meshes in an extrusion process
CN111581776A9 (zh) 一种基于几何重建模型的等几何分析方法
CN117047286B (zh) 一种激光加工工件表面的方法、加工系统、处理器和存储介质
CN108335354A (zh) 一种基于移动最小二乘法的贪婪投影三角化算法单木重建方法
Gargallo-Peiró et al. Mesh generation for atmospheric boundary layer simulation in wind farm design and management
CN113849923B (zh) 一种基于abaqus的钢结构整体焊接变形的数值模拟方法
CN103679814A (zh) 一种三角形网格模型的处理方法及装置
CN111445569A (zh) 一种沉积地质演化动态模拟方法
KR101462861B1 (ko) 선박 곡외판 열간 가공 정보 산출 방법 및 시스템
CN113536617A (zh) 一种复杂结构的精细有限元模型快速生成方法
CN109885946B (zh) 一种确定复合热源的能量分配的方法及焊接模拟方法
CN115906528A (zh) 一种焊接结构模型的自动预处理方法、系统、设备和存储介质
CN115795971A (zh) 一种适用于复杂焊接结构疲劳裂纹扩展分析的网格生成方法、装置及设备
Stepanenko et al. Technique of increasing the accuracy of GTE parts manufactured by selective laser melting
CN111950182A (zh) 一种提高大型复杂构件有限元计算效率的模型简化方法
Patel et al. Stitching and filling: Creating conformal faceted geometry
Barazzetti et al. Vault modeling with neural networks
CN113849924B (zh) 一种基于abaqus的钢结构焊接残余应力和变形的方法和系统

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20220413

Address after: No.15, South Lishi Road, Fuwai, Xicheng District, Beijing

Applicant after: Beijing Machinery Construction Group Co.,Ltd.

Applicant after: Beijing construction industry Refco Group Ltd

Address before: No.15, South Lishi Road, Fuwai, Xicheng District, Beijing

Applicant before: Beijing Machinery Construction Group Co.,Ltd.

GR01 Patent grant
GR01 Patent grant