CN107016192A - 焊接数值模拟中动态生死单元的有限元建模方法 - Google Patents
焊接数值模拟中动态生死单元的有限元建模方法 Download PDFInfo
- Publication number
- CN107016192A CN107016192A CN201710223624.3A CN201710223624A CN107016192A CN 107016192 A CN107016192 A CN 107016192A CN 201710223624 A CN201710223624 A CN 201710223624A CN 107016192 A CN107016192 A CN 107016192A
- Authority
- CN
- China
- Prior art keywords
- welding
- weld metal
- heat
- unit
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提供了一种焊接数值模拟中动态生死单元的有限元建模方法,包括第一步:建立焊接结构的三维实体模型,并根据焊接材料将三维实体模型分为母材区和至少一个焊缝区;第二步:对三维实体模型进行几何清理,并将母材区和焊缝区创建为新的层或集合;第三步:将母材区划分成体网格,并保存到母材区所在的层或集合中;以及,在焊缝区建立面网格,并保存到焊缝区所在层或集合中;第四步:定义焊缝单元,生成焊接有限元分析模型;第五步:定义材料属性、单元的物理属性、分析步数;第六步:定义焊缝区金属的“生死单元”,完成建模。利用本发明可以准确模拟焊缝金属随焊接过程动态填充,新旧换热面动态变换的过程,提高焊接数值计算的精度。
Description
技术领域
本发明涉及焊接数值模拟有限元分析技术领域,更为具体地,涉及一种焊接数值模拟中动态生死单元的有限元建模方法。
背景技术
“生死单元”是现代焊接数值仿真领域常用的模拟方法,其作用是在焊接有限元计算过程中焊缝单元随焊接过程动态的填充焊缝,与实际的焊接过程相符以获得与实际焊接过程相吻合的温度场,提高数值仿真的精度。Python为大型通用ABAQUS软件的二次开发语言,可以对ABAQUS进行建模、计算一体化,尤其在焊接领域中要使用Python二次开发进行焊接过程“生死单元”的定义;同时Tcl是高级网格前处理HyperMesh软件的二次开发语言,可进行复杂模型网格前处理,可以实现不同行业用户的特殊需求而进行二次开发。
如何实现焊接“生死单元”的自动化建模与定义相关的边界条件,在保证温度场计算精度的同时能简化程序,是焊接数值模拟能在工程上进行应用的前提。基于大型通用软件ABAQUS的传统“生死单元”技术是在ABAQUS里采用Python进行编程、建模、仿真与计算,由于ABAQUS自带的建模与网格划分灵活性不强,复杂网格处理较为困难,使得编写的程序十分复杂,尤其在多层多道焊接中实用性更加不强,而工程焊接过程大多采用多层多道焊,这就限制了焊接数值模拟的应用范围。另一方面,从外部导入ABAQUS的复杂的焊接网格模型由于没有实体模型,传热边界条件无法准确的定义到每一个相关的单元面上,一般的处理方法是采用等效的简化,如采用一个总的传热系数代替热对流与热辐射,这所带来的影响就是由于边界条件没有准确定义使得数值模拟结果精度不高。
发明内容
鉴于上述问题,本发明的目的是提供一种焊接数值模拟中动态生死单元的有限元建模方法,以解决现有焊接有限元分析中编程复杂,温度场计算结果精度不高的问题。
本发明提供的焊接数值模拟中动态生死单元的有限元建模方法,包括:
步骤1:建立焊接结构的三维实体模型,并根据焊接材料将三维实体模型分为母材区和至少一个焊缝区;
步骤2:对三维实体模型进行几何清理,并将母材区和焊缝区创建为新的层或集合;
步骤3:将母材区划分成体网格,并保存到母材区所在的层或集合中;以及,在焊缝区建立面网格,并保存到焊缝区所在层或集合中;
步骤4:将焊缝区的面网格拉伸生成体网格,当每次填充焊缝金属单元时新建一个新的层或集合,提取当前焊接位置所有单元的外表面作为换热面并生成一个集合或层的描述,循环操作,直到全部焊缝区完成新的层或集合的建立及单元面的提取,生成焊接有限元分析模型;
步骤5:定义焊接材料的材料属性并赋予给输入焊接有限元分析模型的单元,以及,定义单元类型为八节点六面体;
步骤6:建立第一分析步骤;其中,第一分析步骤的分析时间为1*10^6秒,不激活每个焊缝区的焊缝金属单元;
步骤7:将第一次提取的单元的换热面所在的集合或层的描述赋予对流换热系数与辐射换热系数;
步骤8:建立第二分析步骤;其中,第二分析步骤的分析时间的秒数为总时间除以换热面的总数,激活第二次提取的单元的换热面所在的集合或层的描述,不激活第一次提取的单元的换热面所在的集合或层的描述,循环操作,直到完成所有的焊接过程的定义,建立生死单元的有限元模型。
利用上述根据本发明提供的焊接数值模拟中动态生死单元的有限元建模方法,在焊接过程中,每一步都将整个换热单元面进行更替,而不是基于焊缝金属单元的单元面的激活与杀死,精准的定义了换热单元面的同时避免了复杂的程序设计。
为了实现上述以及相关目的,本发明的一个或多个方面包括后面将详细说明的特征。下面的说明以及附图详细说明了本发明的某些示例性方面。然而,这些方面指示的仅仅是可使用本发明的原理的各种方式中的一些方式。此外,本发明旨在包括所有这些方面以及它们的等同物。
附图说明
通过参考以下结合附图的说明,并且随着对本发明的更全面理解,本发明的其它目的及结果将更加明白及易于理解。在附图中:
图1为根据本发明实施示例的焊缝金属单元的填充过程图;
图2为图1的局部放大图。
具体实施方式
本发明提供的焊接数值模拟中动态生死单元的有限元建模方法,包括如下步骤:
步骤1:建立焊接结构的三维实体模型,并根据焊接材料将三维实体模型分为母材区和至少一个焊缝区。
通过三维软件建立焊接结构的三维实体模型。
在将三维实体模型分为母材区和至少一个焊缝区后,将母材区与焊缝区存储为STP通用数据格式。
步骤2:对三维实体模型进行几何清理,并将母材区和所述焊缝区创建为新的层或集合。
通过HyperMesh软件读入STP通用数据格式的文件,并对三维实体模型进行几何清理。
对三维实体模型进行几何清理的过程为:删除合并整个三维实体模型中多余的点线面体。
步骤3:将母材区划分成体网格,并保存到母材区所在的层或集合中;以及,在焊缝区建立面网格,并保存到焊缝区所在层或集合中。
分别对母材区和焊缝区的体网格进行疏密过渡以保证焊缝区的网格质量,母材区和所述焊缝区的体网格均为完全六面体网格。
步骤4:将焊缝区的面网格拉伸生成体网格,当每次填充焊缝金属单元时新建一个新的层或集合,提取当前焊接位置所有单元的外表面作为换热面并生成一个集合或层的描述,循环操作,直到全部焊缝区完成新的层或集合的建立及单元面的提取,生成焊接有限元分析模型。
换热面为与单元与空气接触的面,本发明以单元的外表面作为换热面。
在每次填充焊缝金属单元时,可以一次填充一个或多个焊缝金属单元,对应地,新建一个集合或层,在每新建一个集合或层时,对该集合或层的描述进行命名,例如:当填充一个焊缝金属单元新建一个新的层或集合时,将生成的集合或层的描述命名为Surface-0,当再次填充一个焊缝金属单元新建一个新的层或集合时,将生成的集合或层的描述命名为Surface-1,以此类推,直到全部焊缝区完成新的层或集合的建立及单元面的提取。对于一次填充多个焊缝金属单元时,同理适用。
可以用Surface-0表示提取第一个单元的换热面,Surface-1表示提取第二个单元的换热面,每一次填充焊缝金属单元时,均会提取一次换热面,采用新提取的换热面替换旧的换热面,例如用Surface-1表示的换热面代替Surface-0表示的换热面。
步骤4通过TCL程序语言编程。
步骤5:利用合并命令将步骤4中所生成的所有重复的节点和单元合并统一。
步骤6:将步骤2、步骤3、步骤4、步骤5过程中所产生的面网格删除。
在建立焊接有限元分析模型后需要删除面网格,删除面网格的目的在于确保焊接有限元分析模型中只有实体网格。
步骤7:检查步骤6中网格数据的正确性,并输出为Inp后缀的ABAQUS文件。
在输出的Inp文件时,必须进行网格数据有效性的检查,确保无重复节点、单元,以及命名错误等问题,如果发现重复或命名错误问题,返回HyperMesh软件重新运行调试TCL程序。
步骤8:将步骤7所产生的Inp文件读入ABAQUS软件。
在读入步骤7所产生的Inp文件时,检查数据交接时步骤4所建立的层或集合和换热面数据的正确性。
步骤9:定义焊接材料的材料属性并赋予给输入焊接有限元分析模型的单元,以及,定义单元类型为八节点六面体。
材料属性包括材料密度、热膨胀系数、导热系数、对流换热系数、辐射换热系数。
步骤10:定义初始常数Constant,如玻尔兹曼常数、绝对零度。
步骤11:定义Dflux焊接热源集中热流载荷。
Dflux为焊接体热源或面热源的Fortran编写的用户子程序。
需要说明的是,步骤9-11不分先后顺序。
步骤12:建立第一分析步骤;其中,第一分析步骤的分析时间为1*10^6秒,不激活每个焊缝区的焊缝金属单元。
由于最开始时没有焊接,因此,需要将所有单元定义为“死”,即不激活所有单元与单元的换热面。同理,在激活单元时,激活单元及单元的换热面。
步骤13:将第一次提取的单元的换热面所在的集合或层的描述赋予对流换热系数与辐射换热系数。
根据对集合或层的描述的命令来对集合或层的描述赋予对流换热系数与辐射换热系数。
将第一次提取的单元的换热面所在的集合或层的描述赋予对流换热系数与辐射换热系数,即将Surface-0所在的集合或层的描述赋予对流换热系数与辐射换热系数。
步骤14:建立第二分析步骤;其中,第二分析步骤的分析时间的秒数为总时间除以换热面的总数,激活第二次提取的单元的换热面所在的集合或层的描述,不激活第一次提取的单元的换热面所在的集合或层的描述,循环操作,直到完成所有的焊接过程的定义,建立生死单元的有限元模型。
第二次提取的单元的换热面所在的集合或层的描述,即为Surface-1所在的集合或层的描述。
可以理解为Surface-1相对于Surface-0是新的换热面,而Surface-1相对于Surface-2是旧的换热面,每次填充金属焊缝单元时,用新的换热面替换旧的换热面,即在激活新的焊缝金属单元的同时激活新的换热面,同时杀死旧的换热面。
步骤15:检查步骤14中定义有无Bug,完成整个焊接过程“生死单元”的FEA建模过程,新建Job,写入INP文件,等待计算。
步骤8-14通过Python程序实现。
在实际的焊接的冷却过程中,焊缝在高温下的散热方式主要是由热辐射,其次为与空气的对流换热,以及焊接结构中的热传导,在采用“生死单元”有限元数值计算模拟焊接的过程中,焊缝金属为了与真实的焊接过程相符往往采用逐步填充的方式进行,以得到比较精确的温度场。如图1所示,在焊缝金属单元填充的过程中,伴随着新旧换热面的交替过程,在数值模拟中,每一步需要定义新的换热面的换热系数,同时抑制旧的换热面的换热,这个循环的过程采用TCL/Python程序进行编程设计,与一般的方法不同的是每一步都将整个换热面进行更替,而不是基于焊缝单元的换热面的激活与杀死,精准的定义了换热单元面的同时避免了复杂的程序设计。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (4)
1.一种焊接数值模拟中动态生死单元的有限元建模方法,包括:
步骤1:建立焊接结构的三维实体模型,并根据焊接材料将所述三维实体模型分为母材区和至少一个焊缝区;
步骤2:对所述三维实体模型进行几何清理,并将所述母材区和所述焊缝区创建为新的层或集合;
步骤3:将所述母材区划分成体网格,并保存到所述母材区所在的层或集合中;以及,在所述焊缝区建立面网格,并保存到所述焊缝区所在层或集合中;
步骤4:将所述焊缝区的面网格拉伸生成体网格,当每次填充焊缝金属单元时新建一个新的层或集合,提取当前焊接位置所有单元的外表面作为换热面并生成一个集合或层的描述,循环操作,直到全部焊缝区完成新的层或集合的建立及换热面的提取,生成焊接有限元分析模型;
步骤5:定义所述焊接材料的材料属性并赋予给输入所述焊接有限元分析模型的单元,以及,定义模型单元的类型为八节点六面体;
步骤6:建立第一分析步骤;其中,所述第一分析步骤的分析时间为1*10^6秒,不激活每个焊缝区的焊缝金属单元;
步骤7:将第一次提取的单元的换热面所在的集合或层的描述赋予对流换热系数与辐射换热系数;
步骤8:建立第二分析步骤;其中,所述第二分析步骤的分析时间的秒数为总时间除以所述换热面的总数,激活第二次提取的单元的换热面所在的集合或层的描述,不激活第一次提取的单元的换热面所在的集合或层的描述,循环操作,直到完成所有焊接过程的定义,建立生死单元的有限元模型。
2.如权利要求1所述的焊接数值模拟中动态生死单元的有限元建模方法,其中,分别对所述母材区和所述焊缝区的体网格进行疏密过渡,所述母材区和所述焊缝区的体网格均为八节点六面体单元。
3.如权利要求2所述的焊接数值模拟中动态生死单元的有限元建模方法,其中,将步骤4中所生成的所有重复的节点和单元合并统一。
4.如权利要求2所述的焊接数值模拟中动态生死单元的有限元建模方法,其中,在执行步骤4之后,将步骤1-步骤4中所产生的面网格删除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710223624.3A CN107016192B (zh) | 2017-04-07 | 2017-04-07 | 焊接数值模拟中动态生死单元的有限元建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710223624.3A CN107016192B (zh) | 2017-04-07 | 2017-04-07 | 焊接数值模拟中动态生死单元的有限元建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107016192A true CN107016192A (zh) | 2017-08-04 |
CN107016192B CN107016192B (zh) | 2019-10-01 |
Family
ID=59446623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710223624.3A Active CN107016192B (zh) | 2017-04-07 | 2017-04-07 | 焊接数值模拟中动态生死单元的有限元建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016192B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107862157A (zh) * | 2017-12-01 | 2018-03-30 | 辽宁三三工业有限公司 | 基于ansys的生死单元法模拟隧道的盾构开挖方法 |
CN109740232A (zh) * | 2018-12-26 | 2019-05-10 | 北京工业大学 | 液滴沉积过程模拟的边界条件处理方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103366044A (zh) * | 2013-06-04 | 2013-10-23 | 广东电网公司电力科学研究院 | 一种基于ansys预测cmt焊缝形状尺寸的方法 |
CN104217078A (zh) * | 2014-09-03 | 2014-12-17 | 杭州联源重工机械有限公司 | 一种不规则曲面优化下料的方法 |
CN104899381A (zh) * | 2015-06-12 | 2015-09-09 | 卢伟 | 一种焊接箱型截面钢节点多尺度有限元建模方法 |
CN105740577A (zh) * | 2016-03-08 | 2016-07-06 | 南京航空航天大学 | 一种Invar钢PMIG摆动焊温度场及变形模拟方法 |
WO2016173313A1 (zh) * | 2015-04-27 | 2016-11-03 | 江苏金通灵流体机械科技股份有限公司 | 一种基于ansys的双相不锈钢与异种钢焊接变形预测方法 |
CN106557638A (zh) * | 2016-12-07 | 2017-04-05 | 重庆科技学院 | 焊接混合双向过渡单元网格模型的建立方法 |
-
2017
- 2017-04-07 CN CN201710223624.3A patent/CN107016192B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103366044A (zh) * | 2013-06-04 | 2013-10-23 | 广东电网公司电力科学研究院 | 一种基于ansys预测cmt焊缝形状尺寸的方法 |
CN104217078A (zh) * | 2014-09-03 | 2014-12-17 | 杭州联源重工机械有限公司 | 一种不规则曲面优化下料的方法 |
WO2016173313A1 (zh) * | 2015-04-27 | 2016-11-03 | 江苏金通灵流体机械科技股份有限公司 | 一种基于ansys的双相不锈钢与异种钢焊接变形预测方法 |
CN104899381A (zh) * | 2015-06-12 | 2015-09-09 | 卢伟 | 一种焊接箱型截面钢节点多尺度有限元建模方法 |
CN105740577A (zh) * | 2016-03-08 | 2016-07-06 | 南京航空航天大学 | 一种Invar钢PMIG摆动焊温度场及变形模拟方法 |
CN106557638A (zh) * | 2016-12-07 | 2017-04-05 | 重庆科技学院 | 焊接混合双向过渡单元网格模型的建立方法 |
Non-Patent Citations (5)
Title |
---|
BIN LI: "FE Simulation of Size Effects on Interface Fracture Characteristics of Microscale Lead-Free Solder Interconnects", 《2009 INTERNATIONAL CONFERENCE ON ELECTRONIC PACKAGING TECHNOLOGY & HIGH DENSITY PACKAGING》 * |
周杰: "焊接残余应力对构架疲劳强度的影响", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
姜泽东: "304不锈钢薄板脉冲激光焊焊接热过程数值分析", 《金属铸锻焊技术》 * |
申刚: "基于典型焊接构件模型的有限元六面体网格自动划分算法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
高耀东: "基于ANSYS单元生死技术的焊接模拟", 《热处理技术与装备》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107862157A (zh) * | 2017-12-01 | 2018-03-30 | 辽宁三三工业有限公司 | 基于ansys的生死单元法模拟隧道的盾构开挖方法 |
CN109740232A (zh) * | 2018-12-26 | 2019-05-10 | 北京工业大学 | 液滴沉积过程模拟的边界条件处理方法 |
CN109740232B (zh) * | 2018-12-26 | 2023-12-08 | 北京工业大学 | 液滴沉积过程模拟的边界条件处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107016192B (zh) | 2019-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Navangul et al. | Error minimization in layered manufacturing parts by stereolithography file modification using a vertex translation algorithm | |
CN108062432A (zh) | 一种激光选区熔化过程的数值模拟方法 | |
Wassermann et al. | Integrating CAD and numerical analysis:‘Dirty geometry’handling using the Finite Cell Method | |
Shapiro et al. | Geometric issues in computer aided design/computer aided engineering integration | |
CN109766663B (zh) | 一种超大型压力容器超长焊缝焊接残余应力与变形高效计算处理方法 | |
Martha | Topological and geometrical modeling approach to numerical discretization and arbitrary fracture simulation in three dimensions | |
CN107016192B (zh) | 焊接数值模拟中动态生死单元的有限元建模方法 | |
Duan et al. | High-order overset flux reconstruction method for dynamic moving grids | |
Xia et al. | Highly efficient wall-distance-based parallel unstructured overset grid assembly | |
von Danwitz et al. | Four‐dimensional elastically deformed simplex space‐time meshes for domains with time‐variant topology | |
Wang | Deformation analysis in horizontal stabilizer assembly using FEA modeling and multilevel analysis | |
Lai et al. | Error compensation of free-form surface with critical area based on T-spline surface reconstruction | |
TW202333039A (zh) | 使用超級葉單元的記憶體個體重新配置 | |
CN104657557A (zh) | 采用电磁仿真软件获取芯片级电路电磁辐射特性的方法及装置 | |
Slak et al. | Adaptive RBF-FD method for Poisson’s equation | |
Fan et al. | Arbitrary surface data patching method based on geometric convolutional neural network | |
V Nekhoroshev et al. | Computer simulation of high-speed anodic dissolution processes of geometrically-complex surfaces of GTE details | |
Wang et al. | Stream surface-supported fundamental sheets insertion toward high-quality hex meshing | |
Wang et al. | Dual-output LUT merging during FPGA technology mapping | |
Kulkarni et al. | Computation of midsurface by feature-based simplification–abstraction–decomposition | |
CN112927350B (zh) | 多域地质体模型构建方法、装置、设备及存储介质 | |
Marsala et al. | From cad to representations suitable for isogeometric analysis: a complete pipeline | |
Chen et al. | One-parameter optimal flux reconstruction schemes for adaptive mesh refinement | |
Heibrok et al. | Open-source MBSE workflow for automated spacecraft thermal analyses from a system model IAC-22, D1, 4B, 2, x70017 | |
CN116186860B (zh) | 一种基于bim的大跨度结构出图方法 |
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 |