CN106557638A - 焊接混合双向过渡单元网格模型的建立方法 - Google Patents
焊接混合双向过渡单元网格模型的建立方法 Download PDFInfo
- Publication number
- CN106557638A CN106557638A CN201611112718.5A CN201611112718A CN106557638A CN 106557638 A CN106557638 A CN 106557638A CN 201611112718 A CN201611112718 A CN 201611112718A CN 106557638 A CN106557638 A CN 106557638A
- Authority
- CN
- China
- Prior art keywords
- model
- transition
- grid
- welding
- zone
- 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]
Abstract
本发明提供一种焊接混合双向过渡单元网格模型的建立方法,包括以下步骤:步骤一、建立三维焊接结构几何模型并将其划分为焊缝区、近焊区、过渡区和远焊区;步骤二、进一步将几何模型划分为最小几何模型组元;步骤三、将上述最小几何模型组元的过渡区部分进一步切分为长方体有限元胞体;步骤四、将最小几何模型组元从远焊区到焊缝区依次进行实体有限元网格划分;步骤五、将最小几何模型组元的有限元网格模型进行空间复制,与三维几何模型完全重合;步骤六、将整个有限元模型合并重复节点,删除三维几何模型,完成建模。此方法简单方便、过渡效率高、适应性好、单元规整,既减少了有限元计算时间,又提高了计算精度。
Description
技术领域
本发明涉及焊接数值计算有限元分析CAE领域,特别涉及一种焊接有限元分析用混合双向过渡单元网格模型的建立方法。
背景技术
有限元分析技术由于核心算法的不同可以分为显式和隐式两种计算分析类型(流体除外)。有限元分析以单元网格模型为有限元基础,相同单元数量下,显式计算的速度比隐式计算快很多,隐式计算每一步求解整个模型的刚度矩阵的平衡解;并且隐式计算随着网格单元数量的增加,计算时间成本呈指数增加。焊接有限元模拟属于隐式计算,采用热-力耦合有限元分析方法。由于焊接过程的局部高温与快冷,使得焊接模拟的过程中,为了获得较为准确的结果,在焊缝区域必须采用较为密集的网格单元,在远离焊缝区由于影响较小,可以采用较稀的单元分布,这样的疏密过渡以达到减少单元数量,在保证计算精度的同时降低计算成本的作用。
传统的采用等比数列节点分布的由密集到稀疏的过渡过程(只针对六面体),单元数量较多,计算成本太大,在复杂的模型或是多层多道焊中并不实用。由此产生了过渡单元的网格划分方式,现有的单元过渡形式有单向过渡,分为2:1和3:1两种过渡方式;这种单向过渡方式存在以下几个问题:单元的长宽比较大;随焊缝区网格的加密,单元的长宽比非常大,单元数量较多,因此此单元只适用于厚度较小的模型、单元的层数不多、焊缝区不密集的焊接计算,不适用于厚板和焊缝区较为密集的有限元焊接计算。
为了解决上述问题,不完全双向过渡单元由此产生。不完全双向过渡单元解决了在焊缝区网格较密和焊接结构厚度较大时出现的单元奇异问题,单元的过渡效率较高。但是此单元网格模型也存在以下几个问题:同样存在长宽比较大的问题;在计算焊接残余应力和变形时,由于过渡单元网格的异向性,可导致计算结果精度较低;单元的灵活性不够好,无法根据实际的模型调整过渡方式实现灵活的处理。另外,现有的焊接有限元模型有3:1双向过渡单元,在厚度和宽度方向均为六面体3:1过渡,该网格的建立方式是采用C++编程开发,与商业软件相结合,能解决部分上述问题,但是还是存在单元建立的灵活度不够,在厚、薄焊接结构的有限元分析中通用性不强,且建立此单元模型有一定的难度。
发明内容
为了克服上述现有技术的不足,本发明提供了一种焊接混合双向过渡单元网格模型的建立方法。
为解决上述技术问题,本发明采用的一个技术方案是:一种焊接混合双向过渡单元网格模型的建立方法,包括以下步骤:
步骤一、建立三维焊接结构几何模型并将其划分为焊缝区、近焊区、过渡区和远焊区;
步骤二、进一步将几何模型划分为最小几何模型组元;
步骤三、将上述最小几何模型组元的过渡区部分进一步切分为长方体有限元胞体;
步骤四、将最小几何模型组元从远焊区到焊缝区依次进行实体有限元网格划分;
步骤五、将最小几何模型组元的有限元网格模型进行空间复制,与三维几何模型完全重合;
步骤六、将整个有限元模型合并重复节点,删除三维几何模型,完成建模。
其中,由远焊区至焊缝区,网格趋势为逐渐由疏至密;
所述划分的方法为:依托软件,进行手动操作划分单元,画出部分的网格单元,再使用脚本语言进行自动生成;
所述的最小几何模型组元应包含有焊缝区、近焊区、过渡区、远焊区四大区域;
在步骤三中,过渡区网格划分时,对胞体进行网格划分,其形式为在胞体的上、下表面为3:1(或2:1)面网格形式过渡,在两侧表面为2:1(或3:1)面网格过渡;
在步骤四中,过渡区网格过渡前后区域的网格,并且疏密过渡效率为6倍,区域间节点共用,网格节点连续,近焊区与焊缝区网格单元密度相同;
在步骤五中,空间复制有限元网格模型与几何模型网格完全重合,且整体模型为N个正整数倍最小几何模型;
所述有限元网格模型为HEX六面体网格实体单元组成。
本发明的有益效果是:
1、此方法简单方便、过渡效率高、适应性好、单元规整,即减少了有限元计算时间,又提高了计算精度;
2、单元的疏密过渡效率大大提高了,有效的减少了单元的数量,降低了计算时间成本;
3、单元建立的灵活度好,在厚、薄焊接结构的有限元分析中通用性强,且建立此单元模型的难度低;
4、单元模型的灵活性较好,便于适应不同焊缝区单元数量的需要来调整过渡的形式,在不同的厚度的焊接结构中可以灵活的使用混合过渡的方式,且不会由于单元的长宽比过大发生单元奇异的问题,提高了有限元计算的收敛速度;
5、既适用于厚度较小的模型,单元的层数不多、焊缝区不密集的焊接计算,又适用于厚板和焊缝区较为密集的有限元焊接计算。
附图说明
图1是本发明混合双向过渡单元网格模型示意图;
图2为本发明多重双向过渡网格示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1:
一种焊接混合双向过渡单元网格模型的建立方法,包括以下步骤:
步骤一、建立三维焊接结构几何模型并将其划分为焊缝区、近焊区、过渡区和远焊区;
步骤二、进一步将几何模型划分为最小几何模型组元;
步骤三、将上述最小几何模型组元的过渡区部分进一步切分为长方体有限元胞体;
步骤四、将最小几何模型组元从远焊区到焊缝区依次进行实体有限元网格划分;
步骤五、将最小几何模型组元的有限元网格模型进行空间复制,与三维几何模型完全重合;
步骤六、将整个有限元模型合并重复节点,删除三维几何模型,完成建模,如图1所示。
进一步的,由远焊区至焊缝区,网格趋势为逐渐由疏至密;
进一步的,所述划分的方法为:依托软件,进行手动操作划分单元,画出部分的网格单元,再使用脚本语言进行自动生成;
进一步的,所述的最小几何模型组元应包含有焊缝区A、近焊区B、过渡区C、远焊区D四大区域;
进一步的,在步骤三中,过渡区网格划分时,对胞体进行网格划分,其形式为在胞体的上、下表面为3:1(或2:1)面网格形式过渡,在两侧表面为2:1(或3:1)面网格过渡;
进一步的,在步骤四中,过渡区网格过渡前后区域的网格,并且疏密过渡效率为6倍,区域间节点共用,网格节点连续,近焊区与焊缝区网格单元密度相同;
进一步的,在步骤五中,空间复制有限元网格模型与几何模型网格完全重合需保证整体模型为N个正整数倍最小几何模型;
进一步的,所述有限元网格模型均为HEX六面体网格实体单元组成。
本发明提出的混合双向过渡单元的有限元模型,即综合2:1与3:1两种过渡方式,并且实现双向过渡,具体过渡情况由焊接结构的厚度和要求的单元网格尺寸决定。
实施例2:
在焊缝处网格密度很大的情况下还可以采用多重双向过渡单元进行疏密过渡,如图2所示。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (9)
1.一种焊接混合双向过渡单元网格模型的建立方法,其特征在于,包括以下步骤:
步骤一、建立三维焊接结构几何模型并将其划分为焊缝区、近焊区、过渡区和远焊区;
步骤二、进一步将几何模型划分为最小几何模型组元;
步骤三、将上述最小几何模型组元的过渡区部分进一步切分为长方体有限元胞体;
步骤四、将最小几何模型组元从远焊区到焊缝区依次进行实体有限元网格划分;
步骤五、将最小几何模型组元的有限元网格模型进行空间复制,与三维几何模型完全重合;
步骤六、将整个有限元模型合并重复节点,删除三维几何模型,完成建模。
2.根据权利要求1所述的焊接混合双向过渡单元网格模型的建立方法,其特征在于,由远焊区至焊缝区,网格趋势为逐渐由疏至密。
3.根据权利要求1所述的焊接混合双向过渡单元网格模型的建立方法,其特征在于,划分所采用的方法为:依托软件,进行手动操作划分单元,画出部分的网格单元,再使用脚本语言进行自动生成。
4.根据权利要求1所述的焊接混合双向过渡单元网格模型的建立方法,其特征在于,所述的最小几何模型组元包含焊缝区、近焊区、过渡区、远焊区四大区域。
5.根据权利要求1所述的焊接混合双向过渡单元网格模型的建立方法,其特征在于,在步骤三中,过渡区网格划分时,对胞体进行网格划分,其形式为在胞体的上、下表面为3:1面网格形式过渡,在两侧表面为2:1面网格过渡。
6.根据权利要求1所述的焊接混合双向过渡单元网格模型的建立方法,其特征在于,在步骤三中,过渡区网格划分时,对胞体进行网格划分,其形式为在胞体的上、下表面为2:1面网格形式过渡,在两侧表面为3:1面网格过渡。
7.根据权利要求1所述的焊接混合双向过渡单元网格模型的建立方法,其特征在于,在步骤四中,过渡区网格过渡前后区域的网格,并且疏密过渡效率为6倍,区域间节点共用,网格节点连续,近焊区与焊缝区网格单元密度相同。
8.根据权利要求1所述的焊接混合双向过渡单元网格模型的建立方法,其特征在于,在步骤五中,空间复制有限元网格模型与几何模型网格完全重合,整体模型为N个正整数倍最小几何模型。
9.根据权利要求1所述的焊接混合双向过渡单元网格模型的建立方法,其特征在于,所述有限元网格模型为HEX六面体网格实体单元组成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611112718.5A CN106557638B (zh) | 2016-12-07 | 2016-12-07 | 焊接混合双向过渡单元网格模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611112718.5A CN106557638B (zh) | 2016-12-07 | 2016-12-07 | 焊接混合双向过渡单元网格模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106557638A true CN106557638A (zh) | 2017-04-05 |
CN106557638B CN106557638B (zh) | 2019-11-05 |
Family
ID=58445984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611112718.5A Active CN106557638B (zh) | 2016-12-07 | 2016-12-07 | 焊接混合双向过渡单元网格模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106557638B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107016192A (zh) * | 2017-04-07 | 2017-08-04 | 重庆科技学院 | 焊接数值模拟中动态生死单元的有限元建模方法 |
CN107169194A (zh) * | 2017-05-11 | 2017-09-15 | 常州轻工职业技术学院 | 一种汽车排气管消音器的有限元建模方法 |
CN109101671A (zh) * | 2018-06-11 | 2018-12-28 | 北京航空航天大学 | 一种变密度与变构型三维点阵结构建模方法 |
CN109614739A (zh) * | 2018-12-24 | 2019-04-12 | 成都安世亚太科技有限公司 | 一种新型的针对带冠叶片的非线性接触处理方法 |
CN109885909A (zh) * | 2019-01-31 | 2019-06-14 | 中国石油大学(华东) | 超大型压力容器局部热处理加固建模及网格划分方法 |
CN110728080A (zh) * | 2018-06-27 | 2020-01-24 | 株洲中车时代电气股份有限公司 | 焊接有限元模型构建方法及校核方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246772A (zh) * | 2013-05-11 | 2013-08-14 | 天津大学 | 一种基于abaqus的超声喷丸矫正焊接变形的有限元模拟方法 |
CN103530463A (zh) * | 2013-10-05 | 2014-01-22 | 黑龙江工程学院 | 一种有限元三维模型网格快速过渡的方法 |
-
2016
- 2016-12-07 CN CN201611112718.5A patent/CN106557638B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103246772A (zh) * | 2013-05-11 | 2013-08-14 | 天津大学 | 一种基于abaqus的超声喷丸矫正焊接变形的有限元模拟方法 |
CN103530463A (zh) * | 2013-10-05 | 2014-01-22 | 黑龙江工程学院 | 一种有限元三维模型网格快速过渡的方法 |
Non-Patent Citations (3)
Title |
---|
毛卫民等: "《金属材料成形与加工》", 31 December 2008, 清华大学出版社 * |
申刚: "《基于典型焊接构件模型的有限元六面体网格自动划分算法研究》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
蹤雪梅等: "《网格过渡形式对焊接残余应力及变形的影响》", 《焊接技术》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107016192A (zh) * | 2017-04-07 | 2017-08-04 | 重庆科技学院 | 焊接数值模拟中动态生死单元的有限元建模方法 |
CN107016192B (zh) * | 2017-04-07 | 2019-10-01 | 重庆科技学院 | 焊接数值模拟中动态生死单元的有限元建模方法 |
CN107169194A (zh) * | 2017-05-11 | 2017-09-15 | 常州轻工职业技术学院 | 一种汽车排气管消音器的有限元建模方法 |
CN107169194B (zh) * | 2017-05-11 | 2020-07-07 | 常州轻工职业技术学院 | 一种汽车排气管消音器的有限元建模方法 |
CN109101671A (zh) * | 2018-06-11 | 2018-12-28 | 北京航空航天大学 | 一种变密度与变构型三维点阵结构建模方法 |
CN109101671B (zh) * | 2018-06-11 | 2023-01-17 | 北京航空航天大学 | 一种变密度与变构型三维点阵结构建模方法 |
CN110728080A (zh) * | 2018-06-27 | 2020-01-24 | 株洲中车时代电气股份有限公司 | 焊接有限元模型构建方法及校核方法 |
CN109614739A (zh) * | 2018-12-24 | 2019-04-12 | 成都安世亚太科技有限公司 | 一种新型的针对带冠叶片的非线性接触处理方法 |
CN109885909A (zh) * | 2019-01-31 | 2019-06-14 | 中国石油大学(华东) | 超大型压力容器局部热处理加固建模及网格划分方法 |
CN109885909B (zh) * | 2019-01-31 | 2023-07-04 | 中国石油大学(华东) | 超大型压力容器局部热处理加固建模及网格划分方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106557638B (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106557638A (zh) | 焊接混合双向过渡单元网格模型的建立方法 | |
CN104036095B (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
CN103699744B (zh) | 一种基于有限元法的风电主控柜体热管理分析及优化方法 | |
CN107273582A (zh) | 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 | |
CN110188423A (zh) | 一种基于有限元网格划分的线性工程结构快速bim建模方法 | |
CN111259547B (zh) | 一种用于综合能源系统运行控制的天然气气路建模方法 | |
CN111523233B (zh) | 一种针对三维压水堆堆芯的中子输运计算方法 | |
CN110046400A (zh) | 一种基于h-p型有限元法对混凝土重力坝进行静力特性分析的高精度数值模拟计算方法 | |
CN106021720A (zh) | 一种大体积混凝土水管冷却温度场的模拟方法 | |
CN108763647A (zh) | 一种浇注机底座轻量化设计方法 | |
CN103729506B (zh) | 一种复杂模型完全六面体建模及几何体重塑加密方法 | |
CN105138772A (zh) | 变截面构件电子束焊的有限元模拟方法 | |
CN106202841A (zh) | 一种机床床身三点支撑位置的优化设计方法 | |
CN109920036A (zh) | 一种基于bim技术的弧形钢结构施工方法 | |
CN106844963B (zh) | 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法 | |
CN110765695B (zh) | 一种基于高阶有限元法获取混凝土重力坝裂纹扩展路径的模拟计算方法 | |
CN102837435A (zh) | 基于中面模型的非等温树脂传递模塑的流动前沿预测方法 | |
CN111553011A (zh) | 一种基于bim的重力坝参数化建模方法 | |
CN109101671A (zh) | 一种变密度与变构型三维点阵结构建模方法 | |
CN111553010A (zh) | 一种基于bim的重力坝智能化设计方法 | |
CN116484668A (zh) | 一种电子束增材制造工艺仿真方法 | |
CN105302979A (zh) | 两相流体网络模型中阀门组的建模方法和系统 | |
CN111797447B (zh) | 一种复杂地基混凝土坝孔道配筋方法 | |
CN108694299A (zh) | 基于icem-cfd的二维有限元中子学稳态计算方法 | |
CN110009746B (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 |