CN107066700B - 一种电弧增材制造有限元建模方法 - Google Patents

一种电弧增材制造有限元建模方法 Download PDF

Info

Publication number
CN107066700B
CN107066700B CN201710167587.9A CN201710167587A CN107066700B CN 107066700 B CN107066700 B CN 107066700B CN 201710167587 A CN201710167587 A CN 201710167587A CN 107066700 B CN107066700 B CN 107066700B
Authority
CN
China
Prior art keywords
arc
additive
equation
living
finite
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
Application number
CN201710167587.9A
Other languages
English (en)
Other versions
CN107066700A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710167587.9A priority Critical patent/CN107066700B/zh
Publication of CN107066700A publication Critical patent/CN107066700A/zh
Application granted granted Critical
Publication of CN107066700B publication Critical patent/CN107066700B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; 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; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING; 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

Abstract

本发明公开了一种电弧增材制造有限元建模方法,包括:计算电弧形态,确定增材区域生死单元宽度;计算熔滴形态,确定增材区域生死单元高度;根据增材区域生死单元宽度和增材区域生死单元高度建立电弧增材制造有限元几何模型;通过验证电弧增材制造有限元几何模型的温度分布,进而修正电弧形态和熔滴形态,从而修正电弧增材制造有限元几何模型;本发明采用电弧增材制造模拟过程结合电弧形态模拟技术、熔滴形态模拟技术和有限元模拟技术,并加以实验验证,能够在电弧增材制造有限元建模之前确定增材区域生死单元几何尺寸,实现电弧增材制造有限元模拟准确计算。

Description

一种电弧增材制造有限元建模方法
技术领域
本发明涉及电弧增材制造领域,具体涉及一种电弧增材制造有限元建模方法。
背景技术
近年来,增材制造技术作为新兴技术发展迅速,不断地被应用于工业生产中。成形材料由最初的光固化树脂不断地向金属材料发展,成形结构也向着大型化整体化发展。采用电弧作为热源的电弧增材制造,具有热输入高,成形速度快,成本低,设备简单等特点,非常适用于大型化整体化结构件,包括大型模具、大型航空发动机、大型飞机壁板等。但是,电弧增材制造过程还存在许多成形缺陷,同时成形过程的物理机理还尚未完成探明。依靠传统的工业试验方法研究成形过程物理机理,不仅存在一定盲目性,而且效率低,耗费生产成本,更棘手的是有些物理参数试验难以获得或者不能获得。因此,采用有限元数值模拟方法研究电弧增材制造自然而然成为大多数学者的选择。
电弧增材制造有限元数值模拟方法的一般步骤包括几何建模、网格划分、材料参数加载,增材区域及增材路径定义、边界条件及热源模型定义、计算求解。因此,在电弧增材制造有限元建模之前需要明确增材区域尺寸、材料参数以及电弧增材工艺参数,但是增材区域尺寸在电弧增材制造实验之前是个未知数,传统的建模方法必须在实验之后才能给出增材区域的尺寸参数,以这种方式建模则违背了数值模拟的原则。虽然,现有的有限元软件在建模之前无能为力,但是其它的流体动力学模拟软件的模拟结果可以为有限元建模提供必要参数。
综上所述,利用流体动力学模拟软件率先进行电弧形态和熔滴形态模拟,得出增材区域宽度和高度,实现电弧增材制造有限元建模,对电弧增材制造数值模拟研究以及物理机理研究具有重要的意义。特别是在目前国内增材制造的工业背景下,研究电弧增材制造物理机理以及实现电弧增材制造建模具有一定的迫切性。
发明内容
本发明所要解决的技术问题是针对上述现有技术的不足提供一种电弧增材制造有限元建模方法,本电弧增材制造有限元建模方法以期准确预测增材区域生死单元的宽度和深度,使电弧增材制造能够顺利建模。同时,还采用实验方法,确定测温点位置,验证温度场模型,然后根据相对误差修正模型,直至获得准确的有限元模型,提高计算机模拟的计算精度。
为实现上述技术目的,本发明采取的技术方案为:
一种电弧增材制造有限元建模方法,包括以下步骤:
步骤1:计算电弧形态,确定增材区域生死单元宽度;
步骤2:计算熔滴形态,确定增材区域生死单元高度;
步骤3:根据增材区域生死单元宽度和增材区域生死单元高度建立电弧增材制造有限元几何模型;
步骤4:通过验证电弧增材制造有限元几何模型的温度分布,进而修正电弧形态和熔滴形态,从而修正电弧增材制造有限元几何模型。
作为本发明进一步改进的技术方案,所述步骤1包括以下步骤:
(1)定义电弧增材制造初始工艺参数,材料热物性参数,电势边界条件;
(2)根据电流连续方程及欧姆定律求解轴向和径向的电流密度:
其中σ是电弧的材料电导率,z是轴向坐标,r是径向坐标,是电势,Jz是轴向的电流密度,Jr是径向的电流密度;
根据安培环流定律求解电弧自感应磁场强度:
其中μ0是真空磁导率,B0是电弧自感应磁场强度;
根据轴向的电流密度Jz、径向的电流密度Jr以及公式(5)的能量方程,求解能量方程中的热量源项Q,热量源项Q为电阻热和电子流传递热之和:
其中KB是玻尔兹曼常数,e是电子电荷,T为电弧形态的温度场,Q是热量源项;
根据电弧自感应磁场强度B0、轴向的电流密度Jz、径向的电流密度Jr以及动量守恒方程,求解动量守恒方程中的动量源项:
FZ=Jr×B0+ρg (6),
Fr=-Jz×B0 (7),
其中ρ是电弧的材料密度,g是重力加速度,FZ是指电弧的轴向电磁力,Fr是指电弧的径向电磁力,FZ和Fr均属于电弧的动量源项;
根据动量方程和连续性方程,求解电弧的轴向速度场和电弧的径向速度场:
其中u是电弧的轴向速度场,v是电弧的径向速度场,t是时间,P是压强,μ是材料黏性系数,其中公式(8)为轴向动量方程,公式(9)为径向动量方程,公式(10)是质量连续性方程,公式(10)是对压强P的修正;
根据热量源项Q、电弧的轴向速度场u和电弧的径向速度场v并结合能量方程确定电弧形态的温度的分布T:
其中CP是电弧的材料等压比热容,k是电弧的材料热导率,T为电弧形态的温度场;
(3)根据电弧形态的温度场T获得加热斑点半径rH,则增材区域生死单元宽度为2rH。作为本发明进一步改进的技术方案,所述步骤2包括以下步骤:
(1)确定熔滴的表面张力、重力、电磁力和等离子体流力:
Fγ=2πRγ (12),
其中R是焊丝半径,γ是熔滴表面张力系数,Fγ为熔滴的表面张力;
其中g是重力加速度,ρ是熔滴的材料密度,Fg为熔滴的重力;
其中表示电流密度,表示磁感应强度的矢量,Fm为熔滴的电磁力;
其中CD是等离子流系数,AP是等离子流力作用面积,ρf是等离子流密度,vf是等离子流的速度,FP为熔滴的等离子体流力;
(2)将熔滴的表面张力、重力、电磁力和等离子体流力在径向的分力表示为Fr1,熔滴的表面张力、重力、电磁力和等离子体流力在轴向的分力表示为FZ1,Fr1和FZ1为熔滴的动量源项;
(3)将Fr1代替步骤1中的Fr,FZ1代替步骤1中的FZ,带入步骤1中的公式(8)、公式(9)和公式(10)求解熔滴的轴向速度场u1和熔滴的径向速度场v1
(4)使用VOF模型,对熔滴运动过程自由界面行为进行追踪:
其中F是单元格流体体积分数;
(5)根据求解的熔滴形态即单元格流体体积分数F,获得熔滴半径rD,则增材区域生死单元高度
作为本发明进一步改进的技术方案,所述步骤3包括以下步骤:
(1)根据实际基板几何尺寸,构建基板几何模型;
(2)根据增材区域生死单元宽度增材区域生死单元高度以及增材区域长度L构建增材区域生死单元几何模型;
(3)对基板几何模型和增材区域生死单元几何模型进行网格划分。
作为本发明进一步改进的技术方案,所述步骤4包括以下步骤:
(1)加载电弧增材制造的初始条件;
(2)加载双椭球热源模型,定义热边界条件和换热边界条件,热边界条件包括环境初始温度以及装夹条件,换热边界条件包括环境初始温度、换热单元以及装夹条件:
公式(17)中qf(x,y,z)为前半椭球的热流密度值,公式(17)中的a和b分别表示前半椭球的长度和深度,cf表示前半椭球的宽度,ff是前半椭球的热量分配系数,公式(18)中qr(x,y,z)为后半椭球的热流密度值,公式(18)中的a和b分别表示后半椭球的长度和深度,cr表示后半椭球的宽度,fr是后半椭球的热量分配系数,公式(17)和公式(18)中x,y,z是三维方向的坐标,Q2是增材制造热输入,e是自然数,fr和ff满足条件:
fr+ff=2 (19);
(3)根据导热偏微分方程,得到电弧增材制造过程的温度场Tc
其中CP是增材材料的材料等压比热容,ρ是增材材料的材料密度,k是增材材料的材料热导率,为内热源强度;
(4)在基板上选取5个红外测温点,沿增材区域的起弧处的垂直方向10mm处取第一红外测温点,沿增材区域的收弧处的垂直方向10mm处取第二红外测温点,再在第一红外测温点和第二红外测温点之间选取3个红外测温点,每个相邻的红外测温点之间的距离均为L/4,L为增材区域长度;
(5)将测量的温度分布表示为Tt并记录测温时间,将测温时间作为冷却时间,进而确定模拟的冷却时间的依据,电弧增材制造冷却过程的模拟时间大于记录测温时间;
(6)提取模拟温度场结果与实验测量结果,计算相对误差:
其中Tc是步骤(3)模拟的电弧增材制造过程的温度场结果、Tt是步骤(5)实验测温结果;
(7)预先设定相对误差范围,若步骤(6)获得的相对误差结果不在预先设定的相对误差范围之内,则修正增材区域生死单元宽度和增材区域生死单元高度进而修正前半椭球的长度、深度和宽度以及后半椭球的长度、深度和宽度,直至获得的模拟温度场结果与实验测量结果在预先设定的相对误差范围内,最终确定增材区域生死单元几何模型。
本发明与现有技术相比,有益效果为:本发明针对电弧增材制造有限元建模前无法预先确定增材区域几何尺寸建模问题,采用电弧形态模拟技术以及熔滴形态模拟技术,可以准确获得增材区域生死单元宽度和增材区域生死单元高度的几何尺寸;本发明可以在电弧增材制造实验前,确定增材区域生死单元宽度与高度,减少电弧增材制造数值模拟过程不必要的实验量;本发明通过数值模拟与实验验证相结合的方法,验证模拟获得的温度场,不断修正模型,直至获得准确的有限元模型,提高计算机模拟的计算精度。
附图说明
图1为本发明的流程示意图。
图2为本发明的电弧增材制造过程电弧形态与加热斑点半径的结构示意图。
图3为本发明的电弧增材制造过程熔滴形态与熔滴半径的结构示意图。
图4为本发明的电弧增材制造有限元几何模型的结构示意图。
图5为本发明的温度场实验验证测温点位置的结构示意图。
图6为本发明的电弧增材制造几何尺寸的结构示意图。
图7为本发明的电弧增材制造温度场模拟结果的结构示意图。
具体实施方式
下面根据图1对本发明的具体实施方式作出进一步说明:
本发明提出一种电弧增材制造有限元建模方法,电弧增材制造模拟过程结合电弧形态模拟技术、熔滴形态模拟技术和有限元模拟技术,并加以实验验证,能够在电弧增材制造有限元建模之前确定增材区域生死单元几何尺寸,实现电弧增材制造有限元模拟准确计算,主要流程如图1所示。本发明通过以下步骤实现:
步骤1:计算电弧形态,确定增材区域生死单元宽度;
步骤2:计算熔滴4形态,确定增材区域生死单元高度;
步骤3:根据增材区域生死单元宽度和增材区域生死单元高度建立电弧增材制造有限元几何模型;
步骤4:通过验证电弧增材制造有限元几何模型的温度分布,进而修正电弧形态和熔滴形态,从而修正电弧增材制造有限元几何模型。
所述步骤1包括以下步骤:
(1)定义电弧增材制造初始工艺参数,材料热物性参数,电势边界条件;
(2)根据电流连续方程及欧姆定律求解电弧的轴向和径向的电流密度:
其中σ是电弧的材料电导率,z是轴向坐标,r是径向坐标,是电势,Jz是轴向的电流密度,Jr是径向的电流密度;
根据安培环流定律求解电弧自感应磁场强度:
其中μ0是真空磁导率,r是径向坐标,B0是电弧自感应磁场强度;
根据轴向的电流密度Jz、径向的电流密度Jr以及公式(5)的能量方程,求解能量方程中的热量源项Q,热量源项Q为电阻热和电子流传递热之和:
其中KB是玻尔兹曼常数,e是电子电荷,T为电弧形态的温度场;
根据电弧自感应磁场强度B0、轴向的电流密度Jz、径向的电流密度Jr以及动量守恒方程,求解动量守恒方程中的动量源项:
FZ=Jr×B0+ρg (6),
Fr=-Jz×B0 (7),
其中ρ是电弧中气体的材料密度,g是重力加速度,FZ是指电弧的轴向电磁力,Fr是指电弧的径向电磁力,FZ和Fr均属于电弧的动量源项;
根据动量方程和连续性方程,求解电弧的轴向速度场u和电弧的径向速度场v:
其中u是电弧的轴向速度场,v是电弧的径向速度场,t是时间,P是压强,μ是材料黏性系数,其中公式(8)为轴向动量方程,公式(9)为径向动量方程,公式(10)是质量连续性方程,公式(10)是对压强P的修正;
根据热量源项Q、电弧的轴向速度场u和电弧的径向速度场v并结合能量方程确定电弧形态的温度的分布T:
其中CP是电弧的材料等压比热容,k是电弧的材料热导率,T为电弧形态的温度场;
(3)参见图2,根据电弧形态的温度场T获得加热斑点半径rH,则增材区域生死单元宽度为
所述步骤2包括以下步骤:
(1)确定熔滴4的表面张力、重力、电磁力和等离子体流力:
Fγ=2πRγ (12),
其中R是焊丝半径,γ是熔滴4表面张力系数,Fγ为熔滴4的表面张力;
其中g是重力加速度,ρ是熔滴4的材料密度,Fg为熔滴4的重力;
其中表示电流密度,表示磁感应强度的矢量,Fm为熔滴4的电磁力;
其中CD是等离子流系数,AP是等离子流力作用面积,ρf是等离子流密度,vf是等离子流的速度,FP为熔滴4的等离子体流力;
(2)将熔滴4的表面张力、重力、电磁力和等离子体流力在径向的分力表示为Fr1,熔滴4的表面张力、重力、电磁力和等离子体流力在轴向的分力表示为FZ1,Fr1和FZ1为熔滴4的动量源项;
(3)将Fr1代替步骤1中的Fr,FZ1代替步骤1中的FZ,带入步骤1中的公式(8)、公式(9)和公式(10)求解熔滴4的轴向速度场u1和熔滴4的径向速度场v1
(4)使用VOF模型,对熔滴4运动过程自由界面行为进行追踪:
其中F是单元格流体体积分数;
(5)根据求解的熔滴4形态即界面的单元格流体体积分数F,获得熔滴4半径rD,则增材区域生死单元高度参见图3;其中熔滴4界面处流体体积分数F的范围为0~1,熔滴4中F为1,其他区域F为0,所以根据熔滴4界面处流体体积分数F的值可计算熔滴4界面的曲率,即可获得熔滴4半径rD
所述步骤3包括以下步骤:
(1)根据实际基板2几何尺寸,构建基板2几何模型;
(2)参见图4,根据增材区域生死单元宽度增材区域生死单元高度以及增材区域长度L构建增材区域生死单元几何模型1;
(3)参见图6,对基板2几何模型和增材区域生死单元几何模型1进行网格划分。
所述步骤4包括以下步骤:
(1)加载电弧增材制造的初始条件;
(2)加载双椭球热源模型,定义热边界条件和换热边界条件,热边界条件包括环境初始温度以及装夹条件,换热边界条件包括环境初始温度、换热单元以及装夹条件:
公式(17)中qf(x,f,z)为前半椭球的热流密度值,公式(17)中的a和b分别表示前半椭球的长度和深度,cf表示前半椭球的宽度,ff是前半椭球的热量分配系数,公式(18)中qr(x,y,z)为后半椭球的热流密度值,公式(18)中的a和b分别表示后半椭球的长度和深度,cr表示后半椭球的宽度,fr是后半椭球的热量分配系数,公式(17)和公式(18)中x,y,z是三维方向的坐标,Q2是增材制造热输入,Q=kUI,是焊接效率、电压和电流的乘积,k是焊接效率,e是自然数,fr和ff满足条件:
fr+ff=2 (19);
(3)根据电弧增材制造的初始条件、热边界条件、换热边界条件和导热偏微分方程,得到电弧增材制造过程的温度场Tc
其中CP是增材材料的材料等压比热容,ρ是增材材料的材料密度,k是增材材料的材料热导率,为内热源强度;
(4)参见图5,在基板2上选取5个红外测温点3,沿增材区域的起弧处的垂直方向10mm处取第一个红外测温点3,沿增材区域的收弧处的垂直方向10mm处取第二个红外测温点3,再在第一个红外测温点3和第二个红外测温点3之间选取3个红外测温点3,每个相邻的红外测温点3之间的距离均为L/4,L为增材区域长度;
(5)将测量的温度分布表示为Tt并记录测温时间,将测温时间作为冷却时间,进而确定模拟的冷却时间的依据,电弧增材制造冷却过程的模拟时间大于记录测温时间;
(6)提取模拟温度场结果与实验测量结果,计算相对误差:
其中Tc是步骤(3)模拟的电弧增材制造过程的温度场结果,Tt是步骤(5)实验测温结果;
(7)预先设定相对误差范围,若步骤(6)获得的相对误差结果不在预先设定的相对误差范围之内,则修正增材区域生死单元宽度和增材区域生死单元高度进而修正前半椭球的长度、深度和宽度以及后半椭球的长度、深度和宽度,直至获得的模拟温度场结果与实验测量结果在预先设定的相对误差范围内,最终确定增材区域生死单元几何模型1。
下面以基板2为Q235钢,电弧增材H13热作模具钢的模拟进行分析。
对于建立电弧增材制造有限元模型,现以H13热作模具钢为例,基板2采用Q235钢,根据实际的电弧增材制造过程,确定电弧增材制造初始工艺参数,如表1所示。
表1 H13电弧增材制造工艺参数
电弧形态及熔滴4形态模拟采用的材料热物性参数如表2所示,由于处于电弧状态,氩气的材料参数都是以10000K下的材料参数为基准。参考电弧增材制造过程的实际情况,对计算机模拟的时间、空间进行离散。主要指定义模拟的时间步长(10-6s),网格尺寸(0.1mm*0.1mm)。
表2电弧形态及熔滴4形态模型中采用的材料热物性参数
步骤a:针对氩气保护气体以及相应的工艺参数,求解步骤1中的偏微分方程组,包括能量方程、动量方程以及质量连续方程,求解获得电弧温度分布状态,根据电弧温度分布状态确定电弧能量范围rH,这个范围根据电弧温度来判断,一般要求电弧温度高于12000,进而确定增材区域生死单元宽度(如图2所示),其中步骤1求解的电弧温度分布状态采用的是氩气的材料热物性参数;
步骤b:针对H13热作模具钢焊丝相应的工艺参数,求解步骤2中的偏微分方程组,特别注意自由界面追踪VOF模型求解,求解获得熔滴4形态,根据熔滴4进入熔池5形态,确定熔滴4半径rD,进而确定增材区域生死单元高度参见图3,其中步骤2求解的熔滴4形态采用的是H13的材料热物性参数;
步骤c:建立电弧增材制造几何模型,其中基板2尺寸为250mm*150mm*4mm,根据增材区域生死单元宽度2rH、增材区域生死单元高度2rD以及增材长度L构建增材区域生死单元几何模型1,如图4和图6所示。
步骤d:针对Q235基板和H13增材材料,加载电弧增材制造初始条件、加载双椭球热源模型,定义热边界条件和换热边界条件,包括环境初始温度、换热单元以及装夹条件等,求解步骤四中温度场微分方程组,获得最终电弧增材制造过程的温度场结果(参见图7,其中黑色的部分是基板,基板上的长方体是增材区域,渐变色是温度场),如图5所述,选取基板2上5个红外测温点,在长度方向选取五个点进行测温,将测量的温度分布表示为Tt并记录测温时间,将测温时间作为冷却时间,进而确定模拟的冷却时间的依据,电弧增材制造冷却过程的模拟时间大于记录测温时间;本例子根据验证实验确定能接受的相对误差为30%,若结果不在误差范围之内则修正电弧形态模型以及熔滴形态模型,直至获得能接受的结果。
本发明针对电弧增材制造有限元建模前无法预先确定增材区域几何尺寸建模问题,采用电弧形态模拟技术以及熔滴形态模拟技术,可以准确获得增材区域生死单元宽度和增材区域生死单元高度的几何尺寸;本发明可以在电弧增材制造实验前,确定增材区域生死单元宽度与高度,减少电弧增材制造数值模拟过程不必要的实验量;本发明通过数值模拟与实验验证相结合的方法,验证模拟获得的温度场,不断修正模型,直至获得准确的有限元模型,提高计算机模拟的计算精度。
本发明的保护范围包括但不限于以上实施方式,本发明的保护范围以权利要求书为准,任何对本技术做出的本领域的技术人员容易想到的替换、变形、改进均落入本发明的保护范围。

Claims (3)

1.一种电弧增材制造有限元建模方法,其特征在于:包括以下步骤:
步骤1:计算电弧形态,确定增材区域生死单元宽度;所述步骤1包括以下步骤:
(1)定义电弧增材制造初始工艺参数,材料热物性参数,电势边界条件;
(2)根据电流连续方程及欧姆定律求解轴向和径向的电流密度:
其中σ是电弧的材料电导率,z是轴向坐标,r是径向坐标,是电势,Jz是轴向的电流密度,Jr是径向的电流密度;
根据安培环流定律求解电弧自感应磁场强度:
其中μ0是真空磁导率,B0是电弧自感应磁场强度;
根据轴向的电流密度Jz、径向的电流密度Jr以及公式(5)的能量方程,求解能量方程中的热量源项Q,热量源项Q为电阻热和电子流传递热之和:
其中KB是玻尔兹曼常数,e是电子电荷,T为电弧形态的温度场,Q是热量源项;
根据电弧自感应磁场强度B0、轴向的电流密度Jz、径向的电流密度Jr以及动量守恒方程,求解动量守恒方程中的动量源项:
FZ=Jr×B0+ρg (6)
Fr=-Jz×B0 (7)
其中ρ是电弧的材料密度,g是重力加速度,FZ是指电弧的轴向电磁力,Fr是指电弧的径向电磁力,FZ和Fr均属于电弧的动量源项;
根据动量方程和连续性方程,求解电弧的轴向速度场和电弧的径向速度场:
其中u是电弧的轴向速度场,v是电弧的径向速度场,t是时间,P是压强,μ是材料黏性系数,其中公式(8)为轴向动量方程,公式(9)为径向动量方程,公式(10)是质量连续性方程,公式(10)是对压强P的修正;
根据热量源项Q、电弧的轴向速度场u和电弧的径向速度场v并结合能量方程确定电弧形态的温度的分布T:
其中CP是电弧的材料等压比热容,k是电弧的材料热导率,T为电弧形态的温度场;
(3)根据电弧形态的温度场T获得加热斑点半径rH,则增材区域生死单元宽度为2rH
步骤2:计算熔滴形态,确定增材区域生死单元高度;所述步骤2包括以下步骤:
(1)确定熔滴的表面张力、重力、电磁力和等离子体流力:
Fγ=2πRγ (12)
其中R是焊丝半径,γ是熔滴表面张力系数,Fγ为熔滴的表面张力;
其中g是重力加速度,ρ是熔滴的材料密度,Fg为熔滴的重力;
其中表示电流密度,表示磁感应强度的矢量,Fm为熔滴的电磁力;
其中CD是等离子流系数,AP是等离子流力作用面积,ρf是等离子流密度,vf是等离子流的速度,FP为熔滴的等离子体流力;
(2)将熔滴的表面张力、重力、电磁力和等离子体流力在径向的分力表示为Fr1,熔滴的表面张力、重力、电磁力和等离子体流力在轴向的分力表示为FZ1,Fr1和FZ1为熔滴的动量源项;
(3)将Fr1代替步骤1中的Fr,FZ1代替步骤1中的FZ,带入步骤1中的公式(8)、公式(9)和公式(10)求解熔滴的轴向速度场u1和熔滴的径向速度场v1
(4)使用VOF模型,对熔滴运动过程自由界面行为进行追踪:
其中F是单元格流体体积分数;
(5)根据求解的熔滴形态即单元格流体体积分数F,获得熔滴半径rD,则增材区域生死单元高度2rD
步骤3:根据增材区域生死单元宽度和增材区域生死单元高度建立电弧增材制造有限元几何模型;
步骤4:通过验证电弧增材制造有限元几何模型的温度分布,进而修正电弧形态和熔滴形态,从而修正电弧增材制造有限元几何模型。
2.根据权利要求1所述的电弧增材制造有限元建模方法,其特征在于:所述步骤3包括以下步骤:
(1)根据实际基板几何尺寸,构建基板几何模型;
(2)根据增材区域生死单元宽度2rH、增材区域生死单元高度2rD以及增材区域长度L构建增材区域生死单元几何模型;
(3)对基板几何模型和增材区域生死单元几何模型进行网格划分。
3.根据权利要求2所述的电弧增材制造有限元建模方法,其特征在于:所述步骤4包括以下步骤:
(1)加载电弧增材制造的初始条件;
(2)加载双椭球热源模型,定义热边界条件和换热边界条件,热边界条件包括环境初始温度以及装夹条件,换热边界条件包括环境初始温度、换热单元以及装夹条件:
公式(17)中qf(x,y,z)为前半椭球的热流密度值,公式(17)中的a和b分别表示前半椭球的长度和深度,cf表示前半椭球的宽度,ff是前半椭球的热量分配系数,公式(18)中qr(x,y,z)为后半椭球的热流密度值,公式(18)中的a和b分别表示后半椭球的长度和深度,cr表示后半椭球的宽度,fr是后半椭球的热量分配系数,公式(17)和公式(18)中x,y,z是三维方向的坐标,Q2是增材制造热输入,e是自然数,fr和ff满足条件:
fr+ff=2 (19)
(3)根据导热偏微分方程,得到电弧增材制造过程的温度场Tc
其中CP是增材材料的材料等压比热容,ρ是增材材料的材料密度,k是增材材料的材料热导率,为内热源强度;
(4)在基板上选取5个红外测温点,沿增材区域的起弧处的垂直方向10mm处取第一红外测温点,沿增材区域的收弧处的垂直方向10mm处取第二红外测温点,再在第一红外测温点和第二红外测温点之间选取3个红外测温点,每个相邻的红外测温点之间的距离均为L/4,L为增材区域长度;
(5)将测量的温度分布表示为Tt并记录测温时间,将测温时间作为冷却时间,进而确定模拟的冷却时间的依据,电弧增材制造冷却过程的模拟时间大于记录测温时间;
(6)提取模拟温度场结果与实验测量结果,计算相对误差:
其中Tc是步骤(3)模拟的电弧增材制造过程的温度场结果、Tt是步骤(5)实验测温结果;
(7)预先设定相对误差范围,若步骤(6)获得的相对误差结果不在预先设定的相对误差范围之内,则修正增材区域生死单元宽度2rH和增材区域生死单元高度2rD进而修正前半椭球的长度、深度和宽度以及后半椭球的长度、深度和宽度,直至获得的模拟温度场结果与实验测量结果在预先设定的相对误差范围内,最终确定增材区域生死单元几何模型。
CN201710167587.9A 2017-03-21 2017-03-21 一种电弧增材制造有限元建模方法 Active CN107066700B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710167587.9A CN107066700B (zh) 2017-03-21 2017-03-21 一种电弧增材制造有限元建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710167587.9A CN107066700B (zh) 2017-03-21 2017-03-21 一种电弧增材制造有限元建模方法

Publications (2)

Publication Number Publication Date
CN107066700A CN107066700A (zh) 2017-08-18
CN107066700B true CN107066700B (zh) 2020-04-07

Family

ID=59620187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710167587.9A Active CN107066700B (zh) 2017-03-21 2017-03-21 一种电弧增材制造有限元建模方法

Country Status (1)

Country Link
CN (1) CN107066700B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108132075B (zh) * 2017-11-14 2019-07-12 西北工业大学 高能束增材制造有限元热力耦合模型的校验方法及其数据测量装置
CN107944154B (zh) * 2017-11-29 2021-03-30 西南交通大学 电弧填丝增材制造分层切片工艺变量自动计算方法
CN108334692B (zh) * 2018-01-30 2020-04-07 清华大学 预测增材制造零件变形的方法
CN109128435A (zh) * 2018-08-20 2019-01-04 南京理工大学 磁场控制式电弧机器人增材制造方法
CN111822828A (zh) * 2020-06-16 2020-10-27 南京航空航天大学 一种基于熔滴过渡的电弧增材成形预测建模方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106180986A (zh) * 2016-07-08 2016-12-07 湘潭大学 一种电弧增材制造成形质量的主动控制方法
CN106216862A (zh) * 2016-07-20 2016-12-14 华中科技大学 一种基于电弧增材和高能束流减材的复合制造方法及装置
CN106425490A (zh) * 2016-09-05 2017-02-22 华中科技大学 一种增减材复合加工设备及其应用

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8315474B2 (en) * 2008-01-18 2012-11-20 Sanyo Electric Co., Ltd. Image processing device and method, and image sensing apparatus

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106180986A (zh) * 2016-07-08 2016-12-07 湘潭大学 一种电弧增材制造成形质量的主动控制方法
CN106216862A (zh) * 2016-07-20 2016-12-14 华中科技大学 一种基于电弧增材和高能束流减材的复合制造方法及装置
CN106425490A (zh) * 2016-09-05 2017-02-22 华中科技大学 一种增减材复合加工设备及其应用

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
铝合金电弧增材制造焊道宽度尺寸预测;柏久阳等;《焊接学报》;20150930;第36卷(第9期);第87-91页 *

Also Published As

Publication number Publication date
CN107066700A (zh) 2017-08-18

Similar Documents

Publication Publication Date Title
US9950476B2 (en) Distortion prediction and minimisation in additive manufacturing
Lee et al. Modeling of heat transfer, fluid flow and solidification microstructure of nickel-base superalloy fabricated by laser powder bed fusion
Rayegani et al. Fused deposition modelling (FDM) process parameter prediction and optimization using group method for data handling (GMDH) and differential evolution (DE)
Zeng et al. A review of thermal analysis methods in laser sintering and selective laser melting
CN102962452B (zh) 基于红外测温图像的金属激光沉积制造扫描路径规划方法
Zeng et al. Comparison of 3DSIM thermal modelling of selective laser melting using new dynamic meshing method to ANSYS
Lopez-Botello et al. Two-dimensional simulation of grain structure growth within selective laser melted AA-2024
Liang et al. A modified method for estimating inherent strains from detailed process simulation for fast residual distortion prediction of single-walled structures fabricated by directed energy deposition
Chen et al. Influence of welding sequence on welding deformation and residual stress of a stiffened plate structure
Foteinopoulos et al. On thermal modeling of Additive Manufacturing processes
US9933255B2 (en) Process mapping of melt pool geometry
Chen et al. An inherent strain based multiscale modeling framework for simulating part-scale residual deformation for direct metal laser sintering
Chowdhury et al. Artificial neural network based geometric compensation for thermal deformation in additive manufacturing processes
Shayfull et al. Potential of conformal cooling channels in rapid heat cycle molding: a review
Zhang et al. A coupled finite element cellular automaton model to predict thermal history and grain morphology of Ti-6Al-4V during direct metal deposition (DMD)
Zhang et al. Inverse identification of interfacial heat transfer coefficient between the casting and metal mold using neural network
Park et al. Optimization of conformal cooling channels with array of baffles for plastic injection mold
Ali et al. Residual stress development in selective laser-melted Ti6Al4V: a parametric thermal modelling approach
Zhao et al. Process parameters optimization of injection molding using a fast strip analysis as a surrogate model
US6813533B1 (en) Method for simulation of laser material deposition
Huang et al. Prediction of laser welding induced deformation in thin sheets by efficient numerical modeling
Criales et al. Predictive modeling and optimization of multi-track processing for laser powder bed fusion of nickel alloy 625
US9283695B1 (en) Computer-implemented simulation method and non-transitory computer medium capable of predicting fiber orientation for use in a molding process
CN107368642A (zh) 金属增材制造多尺度多物理场耦合仿真方法
Salah et al. Numerical study of thermal aspects of electric discharge machining process

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