CN114154359A - 基于ansys模拟热传导过程的双激光冲击锻打加工优化方法 - Google Patents
基于ansys模拟热传导过程的双激光冲击锻打加工优化方法 Download PDFInfo
- Publication number
- CN114154359A CN114154359A CN202111108893.8A CN202111108893A CN114154359A CN 114154359 A CN114154359 A CN 114154359A CN 202111108893 A CN202111108893 A CN 202111108893A CN 114154359 A CN114154359 A CN 114154359A
- Authority
- CN
- China
- Prior art keywords
- laser
- forging
- cladding layer
- cladding
- temperature
- 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.)
- Pending
Links
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/18—Manufacturability analysis or optimisation for manufacturability
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
本发明公开了基于ANSYS数值模拟激光熔覆层热传导过程优化双激光冲击锻打加工参数的方法,该方法包括简化模拟条件;选择并定义有限元单元;建立熔覆基体和熔覆层的三维几何模型;采用生死单元法去除熔覆层模型中的有限元单元并逐步激活;计算得到有限元单元激活时各节点输入的热量载荷;得到熔覆层模型已激活单元节点在激光熔覆过程中的瞬态温度;得到熔覆层的合理锻打时间,对锻打激光和熔覆激光的行进速度进行优化;根据合理锻打区域的大小随时间的变化以及合理锻打区域与最新激活的有限元单元的相对位置,对实际双激光冲击锻打加工的参数以及锻打激光和熔覆激光的相对距离进行优化。由此实现对激光熔覆层热传导过程的可视化监测。
Description
技术领域
本发明涉及激光冲击锻打领域,尤其涉及一种基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法。
背景技术
激光熔覆技术作为代表性的增材制造技术,能将合适的材料涂覆在失效金属的表面上形成熔覆层改善金属的耐磨性,硬度,强度和表面质量等性能,但传统激光熔覆层在冷却过程中与基体之间有较大的温度梯度且两者的收缩速率不同,导致熔覆层内部产生较大的残余拉应力,另外金属冷却过程中相组织的转变和生长也引起金属内部的组织应力,这些应力的叠加造成熔覆层内部气孔和微裂纹的产生。
为了减少熔覆层内部的缺陷,双激光冲击锻打技术在传统激光熔覆技术的基础上,在熔覆层处于高温状态时引入激光冲击的锻打力作用,使熔覆层的组织结构由铸态转变为锻态,熔覆层在快速冷却过程中产生的微裂纹在压应力的作用下焊合,并在压应力的作用下使熔覆层的组织晶粒细化,提高熔覆层的使用寿命,然而熔覆层冷却十分迅速,内部温度场变化、热量传导过程在实际加工过程中无法直接测量和实现可视化,难以合理地控制锻打激光的能量、脉宽、光斑大小、以及与熔覆激光的相对位置等工艺参数,无法很好匹配激光熔覆效果的同时保证激光冲击锻打效果,或者保证激光冲击锻打效果的同时确保熔覆质量,使得熔覆层在其锻造温度区间内得到充分锻造,由此,两者的协同控制以及加工参数匹配仍然是限制加工质量的最大问题。
发明内容
本发明的目的在于提出一种基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,以解决上述背景技术中的一个或多个技术问题。
为达此目的,本发明采用以下技术方案:
基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,该方法包括以下步骤:
步骤S100:简化双激光冲击锻打加工中熔覆层热传导数值模拟条件;
步骤S200:选择有限元单元类型,并定义有限元单元属性;
步骤S300:建立熔覆基体和熔覆层的三维几何模型形成熔覆基体模型和熔覆层模型,并建立空间直角坐标系;
步骤S400:采用生死单元法去除熔覆层模型中的有限元单元并逐步激活;
步骤S500:建立熔池热源模型,计算得到有限元单元激活时各节点输入的的热量载荷;
步骤S600:在最新激活的有限元单元各节点上施加热量载荷,得到熔覆层模型节点在激光熔覆过程中的瞬态温度;
步骤S700:选取熔覆层模型中的瞬态温度数据采集点绘制温度-时间变化曲线,得到熔覆层的合理锻打时间,对实际双激光冲击锻打加工过程中锻打激光和熔覆激光的行进速度进行优化;
步骤S800:将温度处于合理锻打温度区间的有限元单元作为合理锻打区域,根据合理锻打区域的大小随时间的变化以及合理锻打区域与最新激活的有限元单元的相对位置,对实际双激光冲击锻打加工的参数以及锻打激光和熔覆激光的相对距离进行优化。
优选的,步骤S400包括以下步骤:
步骤S401:将熔覆层有限元单元的密度、比热容、导热系数、弹性模量、线膨胀系数和泊松比的值设置为0;
步骤S402:根据实际加工过程中熔覆激光的行进路线、熔覆层的截面尺寸和生长速率计算得到熔覆层模型有限元单元激活的先后顺序和时间间隔;
步骤S403:按步骤S402计算得到的有限元单元激活的先后顺序和时间间隔逐个将熔覆层的有限元单元的质量、比热容、导热系数、弹性模量、线膨胀系数和泊松比恢复为原始值。
优选的,步骤S600包括以下步骤:
步骤S601:通过温度场分析的热平衡定理和非线性瞬态热传导方程并根据步骤S100简化后的热传导数值模拟条件,设定热传导计算初始条件与边界条件;
步骤S602:在最新激活的有限元单元各节点上施加热量载荷;
步骤S603:代入熔覆层模型进行求解,得到熔覆层模型节点在激光熔覆过程中的瞬态温度。
优选的,步骤S601中,设定热传导计算初始条件为:以有限元单元激活时刻为热传导开始时间t=0,基体初始温度为环境温度:
T|t=0=293K
设定热传导计算边界条件为:基体与工作台无热量交换,为绝热边界:
优选的,所述温度场分析的热平衡定理为:
在任意时刻,从熔覆层单位面积向基体内部传递的热量等于熔覆层单位面积对流换热向外界散发的热量以及基体单位面积通过表面辐射向外界散发的热量之和减去熔覆层单位面积吸收激光束的能量,即:
式中,qf为对流换热,qf=hf(Tw-Te),hf为对流换热系数;qr为辐射散热,ε为基体表面辐射率,σ为史蒂芬森玻尔兹曼常数, 5.67×10-8W/(m2 K4),Tw为基体表面温度;Te为环境温度;η为基体与熔覆层材料的激光吸收率;q(x,y,z)为激光束的能量密度,即单位时间内单位截面积上流过的热量;
温度场分析的非线性瞬态热传导方程为:
式中,λx、λy、λz为材料沿x、y、z三个方向的热导率;ρ(T)为材料密度; c(T)为材料比热容;H(T)为材料相变潜热;T为温度分布函数T(x,y,z,t);t为热传导时间。
优选的,步骤S700包括以下步骤:
步骤S710:根据实际加工过程中单道加工时熔覆层的宽度和厚度,将熔覆层模型等距划分成m层n道;
步骤S720:选取合适的单元节点作为瞬态温度数据采集点,得到各道、各层区域在激光熔覆过程中的温度-时间变化曲线;
步骤S730:根据熔覆层材料的合理锻打温度区间在各道、各层区域的温度- 时间变化曲线上进行截取,选取温度-时间变化曲线沿时间正轴方向第一段斜率为负的截取曲线所对应的时间差值作为熔覆层各道、各层区域的合理锻打时间;
步骤S740:通过步骤S730得到的熔覆层各道、各层区域的合理锻打时间对实际双激光冲击锻打加工过程中锻打激光和熔覆激光的行进速度进行优化。
优选的,步骤S720包括以下步骤:
步骤S721:在熔覆层模型的每道上等距截取多段,熔覆层模型每段的中心节点作为数据采集点;将熔覆层模型的各段采集温度值的平均值作为单道区域的温度值;
步骤S722:将熔覆层模型每道的中心节点作为数据采集点,将熔覆层模型的各道采集温度值的平均值作为单层区域的温度值;
步骤S723:将各时间点的温度值光滑连接后得到该区域的温度-时间变化曲线。
优选的,步骤S730中所述合理的锻打温度区间设定为:始锻温度低于该材料的熔点,终锻温度高于该材料的再结晶温度。
优选的,锻打激光加工参数的优化方法为:
若某道或某层区域的合理锻打时间比该有限元单元激活的时间间隔更长,则降低在该道或该层区域的锻打激光和熔覆激光的行进速率;
若合理锻打区域在坐标轴x0y平面投影面积大于锻打激光的光斑直径,则增大锻打激光的光斑;
若合理锻打区域在z轴上的厚度大于熔覆层模型单层的厚度,则增加锻打激光的激光功率;
将所有时间点的合理锻打区域的中心单元节点连接起来,得到锻打激光优化的行进路线;
取所有合理锻打区域中心单元节点与最新激活的有限元单元中心单元节点的距离平均值作为锻打激光和熔覆激光的相对距离。
与现有技术相比,本发明具有以下有益效果:
(1)本发明采用ANSYS数值模拟软件创新性地实现对激光熔覆过程的模拟,对难以测量的激光熔覆层热传导过程进行精确高效的模拟,可实现对激光熔覆层热传导过程的可视化监测;
(2)本发明可以通过非试验的方式对双激光冲击锻打参数进行优化,节省了大量的时间和成本,具有极大的经济效益和社会效益;
(3)采用本发明方法得出的数值模拟结果可以直接指导实际加工,对提高产品,尤其是焊缝的疲劳性能和力学性能,具有十分高效、有益的作用。
附图说明
附图对本发明做进一步说明,但附图中的内容不构成对本发明的任何限制。
图1是本发明其中一个实施例的总体方法流程图;
图2是本发明其中一个实施例双激光冲击锻打示意图;
图3是本发明其中一个实施例的数值模拟示意图。
上述附图中的标号说明:
1-激光锻打后的熔覆层,2-锻打激光光束,3-激光锻打前的熔覆层,4- 熔覆激光光束,5-金属粉末,6-熔池,7-熔覆基体,8-熔覆基体有限元单元, 9-未激活的熔覆层有限元单元,10-激活后的熔覆层有限元单元,11-合理锻打区域,12-最新激活的有限元单元,13-节点,14-笛卡尔直角坐标系,15- 熔覆层模型。
具体实施方式
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。
本实施例的一种基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,参考附图1-3,该包括以下步骤:
步骤S100:为了排除干扰因素,提升数值模拟的运算速度,相对于实际加工的复杂情况简化双激光冲击锻打加工中熔覆层热传导数值模拟条件;
步骤S200:根据数值模拟所需计算速度与精度确定熔覆基体的有限元单元和熔覆层的有限元单元的节点数量与各节点的自由度,确定对应的熔覆基体有限元单元8的单元类型以及熔覆层有限元单元的单元类型;
通过红外温度测量仪测出在实际加工过程中熔覆层的温度变化范围,在熔覆基体和熔覆层材料的物理参数表中查找处于熔覆层温度变化范围内不同温度下熔覆基体与熔覆层材料的密度、比热容、导热系数、弹性模量、线膨胀系数和泊松比的数值,通过插值法获得上述物理参数随温度变化的连续函数,根据物理参数随温度变化的连续函数对熔覆基体有限元单元的属性以及熔覆层有限元单元的属性进行定义;
步骤S300:根据实际双激光冲击锻打加工过程中熔覆基体和熔覆层的形状尺寸在ANSYS软件中建立相应规格尺寸的熔覆基体和熔覆层的三维几何模型形成熔覆基体模型和熔覆层模型15,然后根据数值模拟所需计算速度与精度确定网格尺寸,对熔覆基体模型和熔覆层模型进行网格划分并建立笛卡尔空间直角坐标系14;
步骤S400:采用生死单元法在熔覆层模型中去除所有有限元单元,然后根据实际加工过程中熔覆激光的行进路线、熔覆层的截面尺寸和生长速率计算得到熔覆层模型有限元单元激活的先后顺序和时间间隔并逐步激活未激活的有限元单元9,以模拟在熔覆激光作用下金属粉末熔化后逐步填充形成熔覆层的过程,激活后的熔覆层有限元单元10紧接着参与之后的热传导数值模拟计算;熔覆激光作用下金属粉末熔化后逐步填充形成熔覆层的过程参考附图2,熔覆激光光束4和金属粉末5在熔覆基体7上形成熔池6,熔池6固化形成激光锻打前的熔覆层3,激光锻打前的熔覆层3被锻打激光光束2锻打后最终形成激光锻打后的熔覆层1。
步骤S500:根据熔覆激光光斑的形状大小、能量分布,建立合适的熔池热源模型,再将激光功率、熔覆基体和熔覆层材料的激光吸收率代入熔池热源模型模型计算得到激活的熔覆层模型有限元单元各节点13的热量载荷。
步骤S600:每次激活新的有限元单元后,在最新激活的有限元单元12各节点上施加热量载荷,代入熔覆层有限元模型进行求解,得到熔覆层模型节点在激光熔覆过程中的瞬态温度;
步骤S700:选取熔覆层模型中的瞬态温度数据采集点绘制温度-时间变化曲线,得到熔覆层的合理锻打时间,对实际双激光冲击锻打加工过程中锻打激光和熔覆激光的行进速度进行优化;
步骤S800:在熔覆层热传导数值模拟过程中,将温度处于合理锻打温度区间的有限元单元作为合理锻打区域11并在ANSYS软件中用绿色高亮显示,正在激活的有限元单元12则用红色高亮显示,以将合理锻打区域与正在激活的有限元单元区域在ANSYS软件中作区分;根据绿色高亮区域范围大小随时间的变化以及绿色高亮区域与红色高亮区域的相对位置,对分别对实际双激光冲击锻打加工过程中锻打激光的激光功率、光斑大小、行进路线以及锻打激光和熔覆激光的相对距离进行优化;
步骤S900:判断激光熔覆与锻打效果,若优化结果符合加工要求则结束数值模拟,若不符合加工要求,则返回步骤S2重新进行模拟和优化分析。
优选的,步骤S1中对熔覆层热传导数值模拟简化的条件有:
(a)将基体材料与金属粉末材料设置为具有各向同性的同种材料;
(b)忽略基体表面和金属粉末熔化后形成的熔池内部流体的流动;
(c)忽略金属粉末熔覆过程中的气化作用;
(d)忽略熔池内部因元素扩散而导致的物理性质和化学性质的改变;
(e)忽略基体与加工平台间因温度差造成的热传导;
(f)忽略激光能量损耗,即假定激光能量被金属粉末和基体材料完全吸收并用于升温。
优选的,步骤S4具体包括以下步骤:
步骤S401:将熔覆层模型有限元单元的密度、比热容、导热系数、弹性模量、线膨胀系数和泊松比乘上一个极小的因子(如1×10-6),使它们的值为0或近似为0,此时熔覆层有限元单元可视为不存在,即完成去除所有有限元单元的处理;
步骤S402:然后根据实际加工过程中熔覆激光的行进路线、熔覆层的截面尺寸和生长速率计算得到熔覆层模型有限元单元激活的先后顺序和时间间隔;
步骤S403:按步骤S402计算得到的先后顺序和时间间隔将熔覆层有限元单元的质量、比热容、导热系数、弹性模量、线膨胀系数和泊松比恢复为原值,逐步激活熔池所在处的有限元单元,由此来模拟熔覆激光作用下金属粉末熔化后逐步填充形成熔覆层的过程,被激活的熔覆层有限元单元紧接着参与之后的热传导数值模拟计算。
优选的,步骤S600包括以下步骤:
步骤S601:通过温度场分析的热平衡定理和非线性瞬态热传导方程并根据步骤S100简化后的热传导数值模拟条件,设定热传导计算初始条件与边界条件;
步骤S602:在正在激活的有限元单元各节点上施加热量载荷;
步骤S603:代入熔覆层模型进行求解,得到熔覆层模型节点在激光熔覆过程中的瞬态温度。
进一步的,所述步骤S6中设定热传导计算初始条件,以有限元单元12“激活”时刻为热传导开始时间t=0,基体7初始温度为环境温度:
T|t=0=293K
设定热传导计算边界条件,基体7与工作台无热量交换,为绝热边界:
进一步的,所述步骤S6中参考温度场分析的热平衡定理,即在任意时刻,从熔覆层单位面积向熔覆基体内部传递的热量等于熔覆层单位面积对流换热向外界散发的热量以及熔覆基体单位面积通过表面辐射向外界散发的热量之和减去熔覆层单位面积吸收激光束的能量,即:
式中,qf为对流换热,qf=hf(Tw-Te),hf为对流换热系数;qr为辐射散热,ε为基体表面辐射率,σ为史蒂芬森玻尔兹曼常数, 5.67×10-8W/(m2 K4),Tw为熔覆基体表面温度;Te为环境温度;η为熔覆基体与熔覆层材料的激光吸收率;q(x,y,z)为熔覆激光束的能量密度,即单位时间内单位截面积上流过的热量。
进一步的,温度场分析的非线性瞬态热传导方程为:
式中,λx、λy、λz为材料沿x、y、z三个方向的热导率;ρ(T)为材料密度;c(T)为材料比热容;H(T)为材料相变潜热;T为温度分布函数T(x,y,z,t);t为热传导时间。
优选的,如图3所示,所述步骤S7具体包括以下步骤:
步骤S710:根据实际加工过程中单道加工时熔覆层的宽度和厚度,将熔覆层模型15等距划分成m层n道;在熔覆层模型有限元单元堆叠而成的网格阵列中,沿z轴厚度方向将熔覆层模型再划分为3层,每一层再沿x轴方向划分为3 列,每1列为1道;
步骤S720:选取合适的单元节点作为瞬态温度数据采集点,得到各道、各层区域在激光熔覆过程中的温度-时间变化曲线;
步骤S730:根据熔覆层材料的合理锻打温度区间在各道、各层区域的温度- 时间变化曲线上进行截取,选取温度-时间变化曲线沿时间正轴方向第一段斜率为负的截取曲线所对应的时间差值作为熔覆层各道、各层区域的合理锻打时间;
步骤S740:通过步骤S730得到的熔覆层各道、各层区域的合理锻打时间对实际双激光冲击锻打加工过程中锻打激光和熔覆激光的行进速度进行优化。
进一步的,步骤720具体还包括以下步骤:
步骤S721:在熔覆层模型的每道上等距截取4段,熔覆层模型每段的中心节点作为数据采集点;将熔覆层模型的各段采集温度值的平均值作为单道区域的温度值;
步骤S722:将熔覆层模型每道的中心节点作为数据采集点,将熔覆层模型的各道采集温度值的平均值作为单层区域的温度值;
步骤S723:将各时间点的温度值光滑连接后得到该区域的温度-时间变化曲线。
进一步的,所述步骤S730中合理锻打温度区间的确定,是结合熔覆层材料的合金相图、塑性图、变形抗力图和再结晶图,从塑性、变形抗力和锻件组织性能三个方面进行综合分析获得,当温度过高时,金属材料过热产生奥氏体晶粒极速长大,氧原子与其形成网状氧化物及易熔氧化物共晶,导致晶粒之间的结合力降低,此时锻打易使金属材料产生开裂;当温度过低时,金属材料不具备足够的塑性产生变形,没有锻打后的强化效果,应保证材料在温度区间内具有较高的塑性和较小的刚度,由此,合理的锻打温度区间为始锻温度低于该材料的熔点,终锻温度必须高于该材料再结晶温度。
优选的,所述步骤S8中双激光锻打加工参数的优化方法为:
(a)若某道或某层区域的合理锻打时间比该有限元单元激活的时间间隔更长,则在该道或该层区域的实际加工过程中进行锻打激光和熔覆激光的行进速率的调整降低2~10%,以将锻打激光和熔覆激光的行进速率调整至与模拟结果接近;
(b)如果合理锻打区域(即绿色高亮区域)在坐标轴x0y平面投影面积大于锻打激光的光斑直径,则在实际加工过程中将锻打激光的光斑大小增大 5~10%,以使锻打激光的光斑能覆盖合理锻打区域。
(c)如果合理锻打区域在z轴上的厚度大于熔覆层模型15单层的厚度,则在实际加工过程中对锻打激光的激光功率增加2~4%,以保证锻打效果;
(d)最后将所有时间点的合理锻打区域的中心单元节点连接起来,得到锻打激光优化的行进路线;
(e)取所有时间点合理锻打区域中心单元节点与最新激活的有限元单元中心单元节点的距离平均值作为实际加工过程中锻打激光和熔覆激光的相对距离。
以上结合具体实施例描述了本发明的技术原理。这些描述只是为了解释本发明的原理,而不能以任何方式解释为对本发明保护范围的限制。基于此处的解释,本领域的技术人员不需要付出创造性的劳动即可联想到本发明的其它具体实施方式,这些方式都将落入本发明的保护范围之内。
Claims (9)
1.基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,其特征在于,该方法包括以下步骤:
步骤S100:简化双激光冲击锻打加工中熔覆层热传导数值模拟条件;
步骤S200:选择有限元单元类型,并定义有限元单元属性;
步骤S300:建立熔覆基体和熔覆层的三维几何模型形成熔覆基体模型和熔覆层模型,并建立空间直角坐标系;
步骤S400:采用生死单元法去除熔覆层模型中的有限元单元并逐步激活;
步骤S500:建立熔池热源模型,计算得到有限元单元激活时各节点输入的的热量载荷;
步骤S600:在最新激活的有限元单元各节点上施加热量载荷,得到熔覆层模型节点在激光熔覆过程中的瞬态温度;
步骤S700:选取熔覆层模型中的瞬态温度数据采集点绘制温度-时间变化曲线,得到熔覆层的合理锻打时间,对实际双激光冲击锻打加工过程中锻打激光和熔覆激光的行进速度进行优化;
步骤S800:将温度处于合理锻打温度区间的有限元单元作为合理锻打区域,根据合理锻打区域的大小随时间的变化以及合理锻打区域与最新激活的有限元单元的相对位置,对实际双激光冲击锻打加工的参数以及锻打激光和熔覆激光的相对距离进行优化。
2.根据权利要求1所述的基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,其特征在于,步骤S400包括以下步骤:
步骤S401:将熔覆层有限元单元的密度、比热容、导热系数、弹性模量、线膨胀系数和泊松比的值设置为0;
步骤S402:根据实际加工过程中熔覆激光的行进路线、熔覆层的截面尺寸和生长速率计算得到熔覆层模型有限元单元激活的先后顺序和时间间隔;
步骤S403:按步骤S402计算得到的有限元单元激活的先后顺序和时间间隔逐个将熔覆层的有限元单元的质量、比热容、导热系数、弹性模量、线膨胀系数和泊松比恢复为原始值。
3.根据权利要求1所述的基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,其特征在于,步骤S600包括以下步骤:
步骤S601:通过温度场分析的热平衡定理和非线性瞬态热传导方程并根据步骤S100简化后的热传导数值模拟条件,设定热传导计算初始条件与边界条件;
步骤S602:在最新激活的有限元单元各节点上施加热量载荷;
步骤S603:代入熔覆层模型进行求解,得到熔覆层模型节点在激光熔覆过程中的瞬态温度。
5.根据权利要求3所述的基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,其特征在于,所述温度场分析的热平衡定理为:
在任意时刻,从熔覆层单位面积向基体内部传递的热量等于熔覆层单位面积对流换热向外界散发的热量以及基体单位面积通过表面辐射向外界散发的热量之和减去熔覆层单位面积吸收激光束的能量,即:
式中,qf为对流换热,qf=hf(Tw-Te),hf为对流换热系数;qr为辐射散热,ε为基体表面辐射率,σ为史蒂芬森玻尔兹曼常数,5.67×10-8W/(m2 K4),Tw为基体表面温度;Te为环境温度;η为基体与熔覆层材料的激光吸收率;q(x,y,z)为激光束的能量密度,即单位时间内单位截面积上流过的热量;
温度场分析的非线性瞬态热传导方程为:
式中,λx、λy、λz为材料沿x、y、z三个方向的热导率;ρ(T)为材料密度;c(T)为材料比热容;H(T)为材料相变潜热;T为温度分布函数T(x,y,z,t);t为热传导时间。
6.根据权利要求1所述的基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,其特征在于,步骤S700包括以下步骤:
步骤S710:根据实际加工过程中单道加工时熔覆层的宽度和厚度,将熔覆层模型等距划分成m层n道;
步骤S720:选取合适的单元节点作为瞬态温度数据采集点,得到各道、各层区域在激光熔覆过程中的温度-时间变化曲线;
步骤S730:根据熔覆层材料的合理锻打温度区间在各道、各层区域的温度-时间变化曲线上进行截取,选取温度-时间变化曲线沿时间正轴方向第一段斜率为负的截取曲线所对应的时间差值作为熔覆层各道、各层区域的合理锻打时间;
步骤S740:通过步骤S730得到的熔覆层各道、各层区域的合理锻打时间对实际双激光冲击锻打加工过程中锻打激光和熔覆激光的行进速度进行优化。
7.根据权利要求6所述的基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,其特征在于,步骤S720包括以下步骤:
步骤S721:在熔覆层模型的每道上等距截取多段,熔覆层模型每段的中心节点作为数据采集点;将熔覆层模型的各段采集温度值的平均值作为单道区域的温度值;
步骤S722:将熔覆层模型每道的中心节点作为数据采集点,将熔覆层模型的各道采集温度值的平均值作为单层区域的温度值;
步骤S723:将各时间点的温度值光滑连接后得到该区域的温度-时间变化曲线。
8.根据权利要求6所述的基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,其特征在于,步骤S730中所述合理的锻打温度区间设定为:始锻温度低于该材料的熔点,终锻温度高于该材料的再结晶温度。
9.根据权利要求1所述的基于ANSYS模拟热传导过程的双激光冲击锻打加工优化方法,其特征在于,锻打激光加工参数的优化方法为:
若某道或某层区域的合理锻打时间比该有限元单元激活的时间间隔更长,则降低在该道或该层区域的锻打激光和熔覆激光的行进速率;
若合理锻打区域在坐标轴x0y平面投影面积大于锻打激光的光斑直径,则增大锻打激光的光斑;
若合理锻打区域在z轴上的厚度大于熔覆层模型单层的厚度,则增加锻打激光的激光功率;
将所有时间点的合理锻打区域的中心单元节点连接起来,得到锻打激光优化的行进路线;
取所有合理锻打区域中心单元节点与最新激活的有限元单元中心单元节点的距离平均值作为锻打激光和熔覆激光的相对距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111108893.8A CN114154359A (zh) | 2021-09-22 | 2021-09-22 | 基于ansys模拟热传导过程的双激光冲击锻打加工优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111108893.8A CN114154359A (zh) | 2021-09-22 | 2021-09-22 | 基于ansys模拟热传导过程的双激光冲击锻打加工优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114154359A true CN114154359A (zh) | 2022-03-08 |
Family
ID=80462708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111108893.8A Pending CN114154359A (zh) | 2021-09-22 | 2021-09-22 | 基于ansys模拟热传导过程的双激光冲击锻打加工优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114154359A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117862531A (zh) * | 2024-03-12 | 2024-04-12 | 西安空天机电智能制造有限公司 | 用于锻打印的双束激光协同控制方法、装置、设备及介质 |
-
2021
- 2021-09-22 CN CN202111108893.8A patent/CN114154359A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117862531A (zh) * | 2024-03-12 | 2024-04-12 | 西安空天机电智能制造有限公司 | 用于锻打印的双束激光协同控制方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Somashekara et al. | Investigations into effect of weld-deposition pattern on residual stress evolution for metallic additive manufacturing | |
Mukherjee et al. | Mitigation of thermal distortion during additive manufacturing | |
Nie et al. | Experimental study and modeling of H13 steel deposition using laser hot-wire additive manufacturing | |
CN106909714B (zh) | 一种薄壁构件电弧填丝增材制造温度场预测方法 | |
CN105436505B (zh) | 一种用于提高零件表面质量的热等静压成形方法 | |
CN105718690A (zh) | 基于时间和空间主动跟踪的激光3d打印熔池凝固行为数值模拟方法 | |
CN114154359A (zh) | 基于ansys模拟热传导过程的双激光冲击锻打加工优化方法 | |
WO2019049981A1 (ja) | 積層造形物の解析方法及び積層造形物の解析装置、並びに積層造形物の製造方法及び積層造形物の製造装置 | |
CN108038266A (zh) | 一种选择性激光修复过程的数值模拟方法 | |
CN112464533B (zh) | 基于预制坯和锻造全过程数值仿真的涡轮盘渐变组织调控方法 | |
CN107066700A (zh) | 一种电弧增材制造有限元建模方法 | |
Živčák et al. | FEA simulation of thermal processes during the direct metal laser sintering of Ti64 titanium powder | |
CN114818437B (zh) | 钛合金整体叶盘等温锻造工艺的优化方法 | |
CN111062121A (zh) | 基于高度函数-格子玻尔兹曼法的粉末熔化数值模拟方法 | |
Mamalis et al. | Simulation of the precision forging of bevel gears using implicit and explicit FE techniques | |
CN115618503A (zh) | 一种舵翼类结构增材工艺仿真及工艺优化方法 | |
CN103559334A (zh) | 一种层流冷却温度场的建模方法及系统 | |
CN113111549B (zh) | 高炉炉缸浇注修复后的侵蚀模型建模方法及建模系统 | |
CN110110366A (zh) | 基于非线性规划模型的高温作业专用服装设计研究方法 | |
CN116756870A (zh) | 基于数值模拟仿真的电弧堆焊工艺参数优化方法 | |
CN115026136B (zh) | 一种无头轧制楔形过渡过程板形预测方法 | |
Li et al. | Deformation behavior of powder metallurgy connecting rod preform during Hot forging based on Hot compression and finite element method simulation | |
Zhao et al. | Evolution analysis of process-induced residual stress during the manufacture of diesel engine block | |
CN107245680B (zh) | 一种提高复合钎焊铝箔抗下垂性能的热处理工艺 | |
CN111695216A (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 |