CN109740232B - 液滴沉积过程模拟的边界条件处理方法 - Google Patents
液滴沉积过程模拟的边界条件处理方法 Download PDFInfo
- Publication number
- CN109740232B CN109740232B CN201811604200.2A CN201811604200A CN109740232B CN 109740232 B CN109740232 B CN 109740232B CN 201811604200 A CN201811604200 A CN 201811604200A CN 109740232 B CN109740232 B CN 109740232B
- Authority
- CN
- China
- Prior art keywords
- activated
- unit
- layer
- units
- boundary
- 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.)
- Active
Links
- 238000005137 deposition process Methods 0.000 title claims abstract description 15
- 238000004088 simulation Methods 0.000 title claims abstract description 13
- 238000003672 processing method Methods 0.000 title claims abstract description 12
- 239000007788 liquid Substances 0.000 title description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 18
- 210000004027 cell Anatomy 0.000 claims description 40
- 210000003888 boundary cell Anatomy 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 8
- 230000008021 deposition Effects 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000000465 moulding Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011960 computer-aided design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了液滴沉积过程模拟的边界条件处理方法,在生死单元法的基础上,根据几何位置和特征对边界单元进行分类,提高边界条件处理过程的运算效率。本发明首先根据单元三维坐标对其进行一级位置判定和二级位置判定,再根据相邻未激活单元的相对位置进行判定,并对单元所属的边界类型进行编号,实现对26种边界条件的快速分类。该方法能够对计算区域内任何单元最多仅需5次判断,即可准确判断出该单元所属的边界类型。在液滴沉积过程的模拟中使用该方法,可以进一步提高液滴沉积过程模拟的计算效率,在液滴沉积制造领域的实际应用中具有中要推动作用。
Description
技术领域
本发明涉及一种液滴沉积过程模拟的边界条件处理方法,可以将边界单元按照几何位置和相邻未激活单元的相对位置进行分类,能够提高边界条件设置的运算效率及准确性。本发明属于有限元数值算法领域。
背景技术
金属增材制造技术是一种直接成型的制造技术,具有成型过程简单、速度快、损耗低等优势,在航空、航天和汽车工业中备受瞩目。该技术能够结合计算机辅助设计(CAD)和计算机辅助制造(CAM)技术,采用逐层沉积的方式,完成复杂三维结构的制造。
为了提高生产过程中的样件成型质量、降低次品率,有必要根据实际工况下的材料参数、样件几何结构和尺寸等条件,对铺展过程中单元间界面温度变化过程进行监测。根据材料的热容、热导率、线膨胀率、相变潜热等热物理性质,判断出生产过程中是否会产生冷裂纹、粗晶核、缩孔等缺陷。通过数值模拟方法***样件成型质量,可以为调整温度条件和沉积频率等参数提供指导。
生死单元法是金属液滴多层沉积过程的数值模拟方法之一,通过逐个单元激活的方法降低计算量、提高计算效率。该算法需要在每次激活新单元时对温度边界条件进行更新,是该算法的难点问题之一。需要设计出更加高效的分类和处理方法提高运算效率,以提高数值模拟预测在实际工程应用中的价值。
发明内容
本发明基于液滴沉积过程模拟,在生死单元法的基础上,根据几何位置和特征对边界单元进行分类,提高边界条件处理过程的运算效率。分别经过一级位置判定、二级位置判定和单元相对位置判定,最后根据判定结果对该边界单元进行编号,实现边界条件的快速分类。
本发明采用的技术方案为液滴沉积过程模拟的边界条件处理方法,该方法主要包括以下步骤:
1)一级位置判定
将所有网格按照单元的三维坐标对其进行一级位置判定,分为顶侧层、底侧层、左侧层、右侧层、上侧层、下侧层和内部区域;
2)二级位置判定
在一级位置判定的基础之上,根据单元的三维坐标对其进行二级位置判定。其中,顶侧层与底侧层的网格被分为前棱、后棱、左棱、右棱、左前顶点、左后顶点、右前顶点、右后顶点和中央区域;左侧层和右侧层网格被分为前棱、后棱和中央区域;前侧层、后侧层和内部区域不需二级位置判定。
3)相邻单元相对位置判定
在二级位置判定的基础之上,根据相邻未激活单元的相对位置进行判定。根据判定结果,对26种边界单元进行编号,并将内部单元编号为0;
本发明可以对计算区域内任一单元的空间坐标、相邻未激活单元的相对位置进行最多5次判断,即可准确判断出该单元所属的边界类型。通过该方法可以进一步提高液滴沉积过程模拟的计算效率,在液滴沉积制造领域的实际应用中,具有重要推动作用。
附图说明
图1是本发明液滴沉积过程模拟的边界条件处理方法的边界类型图。
图2是本发明液滴沉积过程模拟的边界条件处理方法的一、二级位置判定。
图3.1-图3.7是本发明液滴沉积过程模拟的边界条件处理方法的相邻单元相对位置判定;其中,图3.1.上侧层中央区域,图3.2.下侧层中央区域图,3.3.左侧层中央区域,图3.4.右侧层中央区域,图3.5.前侧层中央区域图,3.6.后侧层中央区域,图3.7.中央区域。
图4是本发明液滴沉积模拟的边界条件处理方法的流程框图。
具体实施方式
1)一级位置判定:
将长方体计算区域根据坐标划分为7个一级位置区域。首先对单元的X坐标进行判断,若X坐标为计算区域内的最大值,则该单元处于计算区域底侧层;若X坐标为计算区域内的最小值,则该单元处于计算区域顶侧层。如果单元X坐标位于最大值与最小值之间,则对单元Y坐标进行判断,若Y坐标为计算区域内的最大值,则该单元处于计算区域右侧层;若Y坐标为计算区域内最小值,则该单元位于计算区域左侧层。如果X、Y坐标均位于最大值与最小值之间,则对单元Z坐标进行判断,若Z坐标为计算区域内的最大值,则该单元处于计算区域前侧层;若Z坐标为计算区域内的最小值,则该单元处于计算区域内后侧层;若Z坐标为最大值与最小值之间,则该单元位于中央区域。
(2)二级位置判定:
根据单元坐标对一级位置进一步细分。如果单元处于顶侧层,对单元Y、Z坐标分别进行判断。当Y坐标为顶侧层内的最大值时,若Z坐标为顶层内的最大值,则该单元处于顶侧层左前顶点;若Z坐标为顶侧层内的最小值,则该单元处于顶侧层的左后顶点;若Z坐标位于顶侧层内最大值与最小值之间,则该单元处于顶侧层的左棱。当Y坐标为顶侧层的最小值时,若Z坐标为顶层内的最小值,则该单元处于顶侧层的右后顶点;若Z坐标为顶侧层的最大值,则该单元处于顶侧层的右前顶点;若Z坐标为顶侧层的最大值与最小值之间,则该单元处于顶侧层的右棱。当Y坐标为顶侧层的最大值与最小值之间时,若Z坐标为顶侧层的最大值,则该单元处于顶侧层的前棱;若Z坐标为顶层的最小值,则该单元处于顶层的后棱;若Z坐标为顶层的最大值与最小值之间,则该单元处于顶层的中央。如果单元处于底侧层,通过与上述相同方法进行分类,可以得到底侧层的前棱、后棱、左棱、右棱、左前顶点、左后顶点、右前顶点、右后顶点和中央区域。如果单元处于右侧层,对单元Z坐标进行判断。若Z坐标为右侧层最大值,则该单元处于右侧层的前棱;若Z坐标为右侧层最小值,则该单元处于右侧层的后棱;若Z坐标为右侧层最大值与最小值之间,则该单元处于右侧层的中央区域。如果单元处于左侧层,通过与上述相同方法进行分类,可以得到左侧层的前棱、后棱和中央区域。
(3)单元相对位置判定及边界分类编号:
将上侧层左后顶点编号为19,左前顶点编号为20,右后顶点编号为21,右前顶点编号为22。将下侧层左后顶点编号为23,左前顶点编号为24,右后顶点编号为25,右前顶点编号为26。对于上侧层的前棱,若其左侧单元未被激活编号为20,若其右侧单元未被激活编号为22,若两侧均激活编号为10;对于上侧层的后棱,若其左侧单元未被激活编号为19,若其右侧单元未被激活编号为21,若两侧均激活编号为9;对于下侧层的前棱,若其左侧单元未被激活编号为24,若其右侧单元未被激活编号为26,若两侧均激活编号为14;对于下侧层的后棱,若其左侧单元未被激活编号为23,若其右侧单元未被激活编号为25,若两侧均激活编号为13;对于上侧层的左棱,若其后侧单元未被激活编号为19,若其前侧单元未被激活编号为18,若两侧均激活编号为7;对于上侧层的右棱,若其后侧单元未被激活编号为21,若其前侧单元未被激活编号为22,若两侧均激活编号为8;对于下层的左棱,若其后侧单元未被激活编号为23,若其前侧单元未被激活编号为24,若两侧均激活编号为11;对于下侧层的右棱,若其后侧单元未被激活编号为25,若其前侧单元未被激活编号为26,若两侧均激活编号为12;对于左侧层的前棱,若其上侧单元未被激活编号为20,若其下侧单元未被激活编号为24,若两侧均激活编号为18;对于左侧层的后棱,若其上侧单元未被激活编号为19,若其下侧单元未被激活编号为23,若两侧均激活编号为17;对于右侧层的前棱,若其上侧单元未被激活编号为22,若其下侧单元未被激活编号为26,若两侧均激活编号为16;对于右侧层的后棱,若其上侧单元未被激活编号为21,若其下侧单元未被激活编号为25,若两侧均激活编号为15。
对于上侧层的中央区域(图3.1),若其左侧单元未激活,则编号为7;若其后侧单元未激活,则编号9;若其右侧单元未激活,则编号8;若其前侧单元未激活,则编号10;若其左侧、后侧单元未激活,则编号为19;若其左侧、前侧单元未激活,则编号为20;若其右侧、后侧单元未激活,则编号为21;若其右侧、前侧单元未激活,则编号为22;若其左前方或者左后方单元未激活,则编号为7;若其右前方或者右后方单元未激活,则编号为8。
对于下侧层的中央区域(图3.2),若其左侧单元未激活,则编号为11;若其后侧单元未激活,则编号13;若其右侧单元未激活,则编号12;若其前侧单元未激活,则编号14;若其左侧、后侧单元未激活,则编号为23;若其左侧、前侧单元未激活,则编号为24;若其右侧、后侧单元未激活,则编号为25;若其右侧、前侧单元未激活,则编号为26;若其左前方或者左后方单元未激活,则编号为11;若其右前方或者右后方单元未激活,则编号为12。
对于左侧层的中央区域(如图3.3),若其上侧单元未激活,则编号为7;若其后侧单元未激活,则编号17;若其下侧单元未激活,则编号为11;若其前侧单元未激活,则编号18;若其上侧和后侧单元未激活,则编号为19;若其上侧和前侧单元未激活,则编号为20;若其下侧和后侧单元未激活,则编号为23;若其下侧和前侧单元未激活,则编号为24;若其上前方或者下前方单元未激活,则编号为18;若其上后方或者下后方单元未激活,则编号为17。
对于右侧层的中央区域(如图3.4),若其上侧单元未激活,则编号为8;若其后侧单元未激活,则编号15;若其下侧单元未激活,则编号为12;若其前侧单元未激活,则编号为16;若其上侧和后侧单元未激活,则编号为21;若其上侧和前侧单元未激活,则编号为22;若其下侧和后侧单元未激活,则编号为25;若其下侧和前侧单元未激活,则编号为26;若其上前方或者下前方单元未激活,则编号为16;若其上后方或者下后方单元未激活,则编号为21。
对于前侧层的中央区域(如图3.5),若其上侧单元未激活,则编号为10;若其右侧单元未激活,则编号16;若其下侧单元未激活,则编号为14;若其左侧单元未激活,则编号为18;若其上侧和右侧单元未激活,则编号为22;若其下侧和右侧单元未激活,则编号为26;若其下侧和左侧单元未激活,则编号为24;若其上侧和左侧单元未激活,则编号为20;若其上左方或者下左方单元未激活,则编号为18;若其上右方或者下右方单元未激活,则编号为16。
对于后侧层的中央区域(如图3.6),若其上侧单元未激活,则编号为9;若其右侧单元未激活,则编号15;若其下侧单元未激活,则编号为13;若其左侧单元未激活,则编号为17;若其上侧和右侧单元未激活,则编号为21;若其下侧和右侧单元未激活,则编号为25;若其下侧和左侧单元未激活,则编号为23;若其上侧和左侧单元未激活,则编号为19;若其上左方或者下左方单元未激活,则编号为17;若其上右方或者下右方单元未激活,则编号为15。
对于内部区域(如图3.7),若其下侧单元未激活,则编号为2;若其上侧单元未激活,则编号为1;若其左侧单元未激活,则编号为3;若其右侧单元未激活,则编号为4;若其前侧单元未激活,则编号为5;若其后侧单元未激活,则编号为6;若其左侧和上侧单元未激活,则编号为7;若其右侧和上侧单元未激活,则编号为8;若其上侧和后侧单元未激活,则编号为9;若其上侧和前侧单元未激活,则编号为10;若其下侧和左侧单元未激活,则编号为11;若其下侧和右侧单元未激活,则编号为12;若其下侧和后侧单元未激活,则编号为13;若其下侧和前侧单元未激活,则编号为14;若其后侧和右侧单元未激活,则编号为15;若其前侧和右侧单元未激活,则编号为16;若其后侧和左侧单元未激活,则编号为17;若其前侧和左侧单元未激活,则编号为18;若其左侧、后侧和上侧未激活,则编号为19;若其左侧、前侧和上侧未激活,则编号为20;若其右侧、后侧和上侧未激活,则编号21;若其右侧、前侧和上侧未激活,则编号22;若其左侧、后侧和下侧未激活,则编号23;若其左侧、前侧和下侧未激活,则编号为24;若其右侧、后侧和下侧未激活,则编号为25;若其右侧、前侧和下侧未激活,则编号为26。当该单元在六个方向上的单元均已激活时,若其左上方、左后方、左下方或左前方单元未激活,则均编号为3;若右上方、右后方、右下方、右前方单元未激活,则编号为4;若后上方或后下方单元未激活,则编号为5;若前上方或前下方单元未激活,则编号为6。当该单元在以上18个方向上的单元均已激活时,若其左后上方、左前上方、左前下方或左后下方单元未激活,则均编号为3;若其右后上方、右前上方、右前下方或右后下方单元未激活,则均编号为4。当该单元在以上26个方向上单元均已激活,则该单元不属于边界单元,编号为0。
Claims (2)
1.液滴沉积过程模拟的边界条件处理方法,其特征在于:该方法主要包括以下步骤:
1)一级位置判定
将所有网格按照单元的三维坐标对其进行一级位置判定,分为顶侧层、底侧层、左侧层、右侧层、上侧层、下侧层和内部区域;
2)二级位置判定
在一级位置判定的基础之上,根据单元的三维坐标对其进行二级位置判定;其中,顶侧层与底侧层的网格被分为前棱、后棱、左棱、右棱、左前顶点、左后顶点、右前顶点、右后顶点和中央区域;左侧层和右侧层网格被分为前棱、后棱和中央区域;前侧层、后侧层和内部区域不需二级位置判定;
3)相邻单元相对位置判定
在二级位置判定的基础之上,判断相邻未激活单元的相对位置,再进行至多三轮的位置判定,包括:①判断上、下、左、右、前、后是否存在未激活单元;②判断左上、左下、左前、左后、右上、右下、右前、右后、上前、上后、下前、下后是否存在未激活单元;③判断左上前、左上后、左下前、左下后、右上前、右上后、右下前、右下后是否存在未激活单元;即可对26种边界单元进行编号,而上述位置判断均不存在未激活单元时,则将单元编号为0。
2.根据权利要求1所述的液滴沉积过程模拟的边界条件处理方法,其特征在于:对计算区域内任一单元的空间坐标、相邻未激活单元的相对位置进行最多5次判断,即可准确判断出该单元所属的边界类型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811604200.2A CN109740232B (zh) | 2018-12-26 | 2018-12-26 | 液滴沉积过程模拟的边界条件处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811604200.2A CN109740232B (zh) | 2018-12-26 | 2018-12-26 | 液滴沉积过程模拟的边界条件处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109740232A CN109740232A (zh) | 2019-05-10 |
CN109740232B true CN109740232B (zh) | 2023-12-08 |
Family
ID=66361291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811604200.2A Active CN109740232B (zh) | 2018-12-26 | 2018-12-26 | 液滴沉积过程模拟的边界条件处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109740232B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002288240A (ja) * | 2001-03-26 | 2002-10-04 | Foundation For The Promotion Of Industrial Science | 圧延解析システム及びプログラム |
CN105760588A (zh) * | 2016-02-04 | 2016-07-13 | 国家***第海洋研究所 | 一种基于二层规则网格的sph流体表面重建方法 |
CN107016192A (zh) * | 2017-04-07 | 2017-08-04 | 重庆科技学院 | 焊接数值模拟中动态生死单元的有限元建模方法 |
CN107545577A (zh) * | 2017-08-23 | 2018-01-05 | 电子科技大学 | 基于神经网络的沉积相图像分割方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2317348B1 (en) * | 2009-10-30 | 2014-05-21 | Services Pétroliers Schlumberger | Method for building a depositional space corresponding to a geological domain |
-
2018
- 2018-12-26 CN CN201811604200.2A patent/CN109740232B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002288240A (ja) * | 2001-03-26 | 2002-10-04 | Foundation For The Promotion Of Industrial Science | 圧延解析システム及びプログラム |
CN105760588A (zh) * | 2016-02-04 | 2016-07-13 | 国家***第海洋研究所 | 一种基于二层规则网格的sph流体表面重建方法 |
CN107016192A (zh) * | 2017-04-07 | 2017-08-04 | 重庆科技学院 | 焊接数值模拟中动态生死单元的有限元建模方法 |
CN107545577A (zh) * | 2017-08-23 | 2018-01-05 | 电子科技大学 | 基于神经网络的沉积相图像分割方法 |
Non-Patent Citations (1)
Title |
---|
基于有限差分非规则边界网格剖分技术;张伟;方健儒;熊守美;郭志鹏;;铸造(第05期) * |
Also Published As
Publication number | Publication date |
---|---|
CN109740232A (zh) | 2019-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105488562B (zh) | 一种基于多因子粒子群算法的不规则件排样方法 | |
CN108062427A (zh) | 基于数值计算的梯度控速降低涡轮盘锻造残余应力的方法 | |
CN104008181B (zh) | 一种基于三维模型特征的电子零件相似数控工艺检索方法 | |
CN104933220B (zh) | 复杂曲面汽车注塑模具高精密制造方法及注塑模具 | |
Zhu et al. | A novel decision-making logic for hybrid manufacture of prismatic components based on existing parts | |
Shao et al. | Evolutionary forging preform design optimization using strain-based criterion | |
Zhao et al. | Multiobjective optimization design of porthole extrusion die using Pareto-based genetic algorithm | |
CN110232742A (zh) | 3d打印分层控制算法 | |
CN113591350B (zh) | 一种材料挤出成形3d打印成形质量提升方法 | |
Zhang et al. | Process planning of automatic wire arc additive remanufacturing for hot forging die | |
CN109034226A (zh) | 一种基于图论的车辆轨迹聚类方法 | |
Sedighi et al. | Optimisation of gate location based on weld line in plastic injection moulding using computer-aided engineering, artificial neural network, and genetic algorithm | |
CN110415331A (zh) | 一种基于点云数据的轮廓或孔洞快速检测判别和孔洞修补方法 | |
Samant et al. | Octree data structure for support accessibility and removal analysis in additive manufacturing | |
CN109740232B (zh) | 液滴沉积过程模拟的边界条件处理方法 | |
CN104657565B (zh) | 近型面水道的热作模具设计方法 | |
CN111951399A (zh) | 一种増减材制造中基于体素曲面距离场的轨迹规划方法 | |
Jauregui-Becker et al. | Performance evaluation of a software engineering tool for automated design of cooling systems in injection moulding | |
Li et al. | Complexity analysis and calculation for sculptured surface in multi-axis CNC machining based on surface subdivision | |
CN103658520B (zh) | 一种金属件铸型剖分方法 | |
CN104014633A (zh) | 一种基于有限元分析法的去除棒料芯部缺陷的模拟方法及基于该模拟方法的冲孔方法 | |
CN111036909A (zh) | 金属快速成型过程中晶粒组织数值预测方法 | |
CN108804791B (zh) | 一种适用于埋入式进气道布局的飞行器参数化方法 | |
Ma et al. | Data driven modeling and optimization for energy efficiency in additive manufacturing process with geometric accuracy consideration | |
Nitnara et al. | Simulation-Based Optimization of Injection Molding Process Parameters for Minimizing Warpage by ANN and GA |
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 |