CN101824514A - 高精度预置式激光熔覆涂层温度场仿真方法 - Google Patents
高精度预置式激光熔覆涂层温度场仿真方法 Download PDFInfo
- Publication number
- CN101824514A CN101824514A CN 201010177126 CN201010177126A CN101824514A CN 101824514 A CN101824514 A CN 101824514A CN 201010177126 CN201010177126 CN 201010177126 CN 201010177126 A CN201010177126 A CN 201010177126A CN 101824514 A CN101824514 A CN 101824514A
- Authority
- CN
- China
- Prior art keywords
- material properties
- thermal
- physical parameter
- base material
- initialization layer
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000004372 laser cladding Methods 0.000 title claims abstract description 24
- 239000011248 coating agent Substances 0.000 title claims abstract description 15
- 238000000576 coating method Methods 0.000 title claims abstract description 15
- 238000004088 simulation Methods 0.000 title abstract description 13
- 238000005272 metallurgy Methods 0.000 claims abstract description 13
- 239000000463 material Substances 0.000 claims description 145
- 239000011148 porous material Substances 0.000 claims description 28
- 238000005253 cladding Methods 0.000 claims description 24
- 238000002844 melting Methods 0.000 claims description 12
- 230000008018 melting Effects 0.000 claims description 12
- 239000000843 powder Substances 0.000 claims description 8
- 238000003723 Smelting Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 239000011800 void material Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000002349 favourable effect Effects 0.000 abstract 1
- 102100023471 E-selectin Human genes 0.000 description 2
- 241000283074 Equus asinus Species 0.000 description 2
- 101000622123 Homo sapiens E-selectin Proteins 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010187 selection method Methods 0.000 description 2
- 238000002076 thermal analysis method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明针对目前预置式激光熔覆过程温度场数值模拟方法所存在的不足,公开了一种高精度预置式激光熔覆涂层温度场仿真方法,它通过将几何模型分为气孔模型和冶金化模型及动态地改变接触热阻进行仿真,充分利用了ANSYS软件强大的模拟功能,所述的结果精度高,真实性好,为激光熔覆研究提供了切实可行的手段。
Description
技术领域
本发明涉及一种热加工温度场的仿真方法,属于计算机技术在工程技术中的具体应用研究,尤其涉及一种预置式激光熔覆过程中的温度场仿真方法,具体地说是一种高精度预置式激光熔覆涂层温度场仿真方法。
背景技术
众所周知,激光熔覆是一种快速加热、快速冷却的工艺方法,利用红外测温仪和热电偶都很难全面测量出激光熔覆过程的温度变化情况,因此利用有限元求解温度场的方法已成为求解激光熔覆过程温度场的一种有效手段。由于压片存在着气孔以及压片与基材之间存在着接触热阻,当压片中某区域发生熔化后,就会把原有的气孔剔除出去,气孔的剔除影响着压片的物性参数和厚度,导致熔化前后压片的材料属性发生变化和厚度的减小;当基材表面产生熔化后,压片与基材之间接触形式发生改变,即间隙接触转变成冶金化结合,导致了接触热阻发生变化,因此气孔的剔除和压片与基材之间接触形式的改变必然影响着温度场(压片预置式激光熔覆过程如图2)。在压片预置式温度场仿真中,还没有人在仿真过程中动态地更改单元的材料属性、以及实现压片厚度变化和接触热阻变化对温度的影响。
在《基于ANSYS的压片预置式激光熔覆温度场仿真》(作者:马浩,王珉等,2009,38(4)《机械制造与自动化》-92-94)中,采用面积等效厚度法把整个压片的厚度折算出熔化后的厚度,这与压片中未熔化部分的厚度不变相矛盾,还忽略了熔覆前后的材料属性变化对其温度场的影响,压片的气孔率较大(文中为37.5%),熔覆前后的材料属性间差别明显较大,忽略材料属性的变化将极大地影响仿真结果。因此,根据激光熔覆的特点发明一种高精度的仿真方法是提高激光熔覆质量,改善复合材料特性的关键。
发明内容
本发明的目的是针对现有的预置式激光熔覆温度场仿真精度不高,难以在工业实际应用的问题,根据预置式激光熔覆中存在气孔和接触热阻变化较大的特点,发明一种高精度预置式激光熔覆涂层温度场仿真方法。
本发明的技术方案是:
一种高精度预置式激光熔覆涂层温度场仿真方法,它包括以下步骤:
(1)在ANSYS软件中建立预置层和基材的几何模型;
(2)根据熔覆前后熔覆材料的物性参数和基材物性参数,定义材料属性;
(3)根据定义的材料属性进行网格划分;
(4)设置接触热阻;
(5)热源加载及温度场求解;
其特征是:
(1)在建立预置层的几何模型时分别建立预置层的冶金化几何模型和气孔几何模型,它们的几何形状都是长方体,并且冶金化几何模型在气孔几何模型的上面,它们的高度关系满足以下条件:
H冶+H气=H
式中:H治为冶金化几何模型的高度,其大小等于熔覆层的厚度;
H气为气孔模型的高度;
H为预置层的厚度;
(2)在预置层和基材的接触表面上设置“间隙接触”和“冶金化结合接触”两个接触对,两个接触对分别对应着两个接触热导率,接触热导率等于接触热阻的倒数,利用生死单元技术使较大的接触热阻起作用,而较小的接触热阻不起作用,
间隙接触和冶金化结合接触热阻计算公式如下:
式中:Rt为接触热阻(m2·℃/W)
hf为接触表面的轮廓最大高度(m)
λ1、λ2和λ3分别为粉末片、基材和空气的热导率(W/m·℃)
k1、k2和k3是常数,k1和k2在0~0.5之间,分别反映粉末片和基材的热导率,热导率越大,其值就越小;对于间隙接触k3=0.696,对于冶金结合k3=0。
所述的根据熔覆前后熔覆材料的物性参数和基材物性参数定义的材料属性为:
对于预置层,定义三种材料属性,分别为材料属性1、2和3,每种材料属性包括比热容、密度、热导率和表面换热系数;材料属性1对应着熔化前的预置层材料属性;材料属性2对应着熔化后的预置层材料性;材料属性3为物性参数,并定义密度和比热容为0,表面换热系数为0,热导率为正交类型,即只进行垂直方面的传热,没有横向的传热;即TXX=TYY=0,TZZ=1×109;
其中材料属性2为已知值,根据材料属性2的物性参数,材料属性1的物性参数数值通过以下公式折算得到:
式中:X1为熔化前预置层材料的物性参数,
X2为预置层熔化后的物性参数,
φ是压片的气孔率;
对于基材,只定义一种材料属性,材料属性编号为4,基材熔化前后对应同一种物性参数。
具体进行网格化划分时应根据选取的材料属性进行,并在选取材料属性1时利用8节点的三维热分析单元SOLID70单元类型,用VMESH命令对冶金化模型和气孔模型进行网格划分;选取定义的材料属性4对基材模型进行网格划分。
所述的热源加载及温度场求解过程为:
一次加载热源,设定一个热源加载的时间,进行一次温度场求解;
一次求解结束后,进行材料属性改变和接触热阻改变的判断,材料属性改变过程:判断冶金化仿真模型中单元的节点是否全部达到预置层的熔点,若达到,则单元材料属性发生改变,否则不发生改变;若与气孔仿真模型单元相连接的冶金化仿真模型单元已改变材料属性,则气孔仿真模型单元也改变材料属性;接触热阻改变过程:若预置层或基材的接触表面上的原死单元中全部节点温度达到预置层或基材的熔点,原死单元改变成活单元,以同样的条件把原活单元变成死单;。
再次载热源,设定一个热源加载时间,进行温度场求解;
以上热源加载,温度场求解,材料属性改变和接触热阻改变的判断,再次加载热源和进行温度场求解,这是一个多次循环过程,通过循环语句实现以上循环功能;
将每个时刻求解出的温度场将存入结果文件中,可以根据时间来查看某时刻的温度场。
本发明的有益效果:
本发明提供了一种压片预置式激光熔覆温度场求解的数值模拟方法。在熔覆材料几何模型建立时首次提出了单独建立冶金化模型和气孔模型,在冶金化仿真模型中体现熔化前后的材料属性变化的影响,在气孔仿真模型中体现了厚度变化的影响。
在以往的压片预置式激光熔覆的温度场求解中,通常把压片建立成一个几何模型,从而忽略了熔化前后的压片厚度变化,本发明提出建立冶金化几何模型和气孔几何模型的思想,通过更改有限元仿真过程的单元材料属性解决了熔化前后熔覆材料厚度的变化对温度场影响的问题,因此本发明得到的温度场更加真实准确地反映粉末预置式激光熔覆的温度变化情况,进而有助于分析激光熔覆过程的温度对组织和性能的影响,并且为讨论激光熔覆的工艺参数对温度甚至组织和性能的影响提供了理论依据。
附图说明:
图1是本发明的热源加载和求解循环过程
图2是本发明的粉末压片激光熔覆过程示意图
图3是本发明的熔覆样品中心部分的实物界面
图4是本发明的有限元仿真模型
图5是激光中心位于X=15mm处的温度场
图6是在X=15截面上6个节点的温度曲线。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
一种高精度预置式激光熔覆涂层温度场仿真方法,它包括以下步骤:
(1)根据待熔覆压片和基材的厚度,建立压片和基材的几何模型
在建立压片几何模型时,分别建立熔覆材料的冶金化模型和气孔模型,并且冶金化模型在气孔模型的上面,它们的高度关系如下:
H冶+H气=H
式中:H冶为冶金化模型的高度,是通过测量熔覆层厚度得到
H为预置压片的厚度
H气为气孔模型的高度
由于基材在融化前后没有体积变化,只需建立一个基材几何模型,具体模型如图3。
(2)根据熔覆前后的熔覆材料物性参数和基材物性参数,定义材料属性
利用MAT命令定义三种压片的材料属性,分别为材料属性1,2和3。材料属性1对应着熔化前的压片材料属性;材料属性2对应着熔化后的压片材料性;材料属性3的物性参数:密度和比热容为0,表面换热系数为0,热导率为正交类型,即TXX=TYY=0,TZZ=1×109。
根据材料属性2的物性参数,材料属性1的物性参数通过折算公式折算得到,折算公式如下:
式中:X1为熔化前熔覆材料的物性参数
X2为熔化后熔覆材料的物性参数
φ是压片的气孔率
对于基材,只定义一种材料属性,为材料属性4。基材熔化前后对应同一种物性参数,即材料属性4。
(3)选取(2)定义的材料属性对模型进行网格划分
选取(2)定义的材料属性1,利用8节点的三维热分析单元SOLID70单元类型,用VMESH命令对冶金化模型和气孔模型进行网格划分;选取(2)定义的材料属性4对基材模型进行网格划分。
(4)根据不同的接触状态,设置不同的接触热阻
间隙接触和冶金化结合的接触热阻计算公式如下:
式中:Rt为接触热阻(m2·℃/W)
hf为接触表面的轮廓最大高度(m)
λ1、λ2和λ3分别为粉末片、基材和空气的热导率(W/m·℃)
k1、k2和k3是常数,k1和k2在0~0.5之间,分别反映粉末片和基材的热导率,热导率越大,其值就越小;对于间隙接触k3=0.696,对于冶金结合k3=0。
具体实施时k1、k2和k3的最佳取值为:
对于间隙接触:k1=0.289,k2=0.113,k3=0.696;
对于冶金化结合:k1=0.289,k2=0.113,k3=0。
利用REAl命令定义两个实常数,实常数的编号为1和2,分别对应间隙接触和冶金化结合时的接触热导率,接触热导率的大小等于接触热阻的倒数。
生成间隙接触所对应的接触对:选取编号为1的实常数和(2)中定义的材料属性1,利用ESURF命令在压片的接触面上生成CONTA174单元;选取(2)所定义的材料属性4,在基材接触面上生成TARGE170单元,CONTA174是8节点的三维面面接触单元,TARGE170是三维目标单元。
生成冶金化结合所对应的接触对:选取编号为2的实常数,按照上面生成间隙接触的接触对的方法来生成接触单元和目标单元。
利用ESEL命令选取编号为2的实常数所对应的接触对,使用EKILL命令杀死该接触对所包含的全部接触单元和目标单元。
(5)热源加载及温度场求解
使用SF命令加载热源,使用DELTIM设定热源加载的时间,用SOLVE命令进行温度场求解。
求解结束后,进行材料属性改变和接触热阻改变的判断,材料属性改变过程:利用VSEL命令选取冶金化模型,从已选取模型中用NSLV命令选取所有的节点,从NSLV选取的节点中再用NSEL命令选取温度达到或超过压片熔点的节点,用ESLN选取由NSEL选取的节点所组成的单元(单元中全部节点必须在NSEL选取的节点内),用EMODIF命令改变ESLN选取的单元的材料属性为(2)所定义的材料属性2;以相同的单元选取方法改变气孔模型单元的材料属性为(2)定义的材料属性3,只要ESLN选取的单元中任何一个节点在NSEL选取的节点内,气孔模型单元的材料属性就发生改变。
接触热阻改变过程:用ESEL命令选取实常数编号为1的接触对,用NSLE命令选取接触对中所有节点,从NSLE选取的节点中用NSEL选取温度达到或超过压片或基材的熔点的节点,从接触对的单元中用ESLN选取接触单元和目标单元(单元中全部节点必须在NSEL选取的节点内),用EKILL命令杀死ESLN选取的单元;以相同的单元选取方法用EALIVE命令使ESLN选取的单元变活。
使用DO循环语句实现热源加载、温度场求解、材料属性改变和接触热阻改变的循环过程。
以下是利用本发明的方法所进行的一个具体实例:
如图1-6所示。
选取的压片材料为NiCoCrAlY,其尺寸为30mm×6mm×0.4mm;基材为GH4033,其尺寸为30mm×6mm×4mm;选用CO2激光器,激光熔覆的工艺参数:功率为1.25KW,扫描速度为400mm/min,光斑直径2mm。具体步骤如下:
(1)几何模型的建立
选取熔覆样品中心部分进行截面截取,截面如附图4,测量出熔覆层的厚度为0.36mm,则气孔模型的高度为0.04mm,冶金化模型在气孔模型的上面。采用对称建模的方法,具体的仿真模型如图4。
(2)定义材料属性
定义三种压片材料属性,分别为材料属性1、2和3
材料属性1为熔化前的压片材料属性,材料属性2为熔化后的压片材料属性,对应的具体物性参数数值如表1。
表1熔化前后的压片的物性参数
材料属性3的物性参数数值:密度、比热容和表面换热系数都为0;热导率为正交类型,即TXX=TYY=0,TZZ=1×109。
定义基材的材料属性,为材料属性4。
(3)选取(2)中定义的材料属性对模型进行网格划分
选取(2)中定义的材料属性1对冶金化模型和气孔模型进行网格划分;选取(2)中定义的材料属性4对基材模型进行网格划分。
(4)根据不同的接触状态,设置不同的接触热阻
利用前述公式求解两个接触热阻,进而得到接触热导率,分别是间隙接触对应的热导率0.88×106W/m2·℃和冶金化结合对应的接触热导率1.3×106W/m2·℃。
(5)进入热源加载和温度场求解的循环过程,过程如图1所示。
求解的温度场如图5,6个节点的温度曲线如图6(6个节点位于X=15mm截面上)。
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。
Claims (3)
1.一种高精度预置式激光熔覆涂层温度场仿真方法,它包括以下步骤:
(1)在ANSYS软件中建立预置层和基材的几何模型;
(2)根据熔覆前后熔覆材料的物性参数和基材物性参数,定义材料属性;
(3)根据定义的材料属性进行网格划分;
(4)设置接触热阻;
(5)热源加载及温度场求解;
其特征是:
(1)在建立预置层的几何模型时分别建立预置层的冶金化几何模型和气孔几何模型,它们的几何形状都是长方体,并且冶金化几何模型在气孔几何模型的上面,它们的高度关系满足以下条件:
H冶+H气=H
式中:H冶为冶金化几何模型的高度,其大小等于熔覆层的厚度;
H气为气孔模型的高度;
H为预置层的厚度;
(2)在预置层和基材的接触表面上设置“间隙接触”和“冶金化结合接触”两个接触对,两个接触对分别对应着两个接触热导率,接触热导率等于接触热阻的倒数,利用生死单元技术使较大的接触热阻起作用,而较小的接触热阻不起作用,
间隙接触和冶金化结合接触热阻计算公式如下:
式中:Rt为接触热阻(m2·℃/W)
hf为接触表面的轮廓最大高度(m)
λ1、λ2和λ3分别为粉末片、基材和空气的热导率(W/m·℃)
k1、k2和k3是常数,k1和k2在0~0.5之间,分别反映粉末片和基材的热导率,热导率越大,k1和k2就越大;对于间隙接触k3=0.696,对于冶金结合k3=0。
2.根据权利要求1所述的高精度预置式激光熔覆涂层温度场仿真方法,其特征是所述的根据熔覆前后熔覆材料的物性参数和基材物性参数定义的材料属性为:
对于预置层,定义三种材料属性,分别为材料属性1、2和3,每种材料属性包括比热容、密度、热导率和表面换热系数;材料属性1对应着熔化前的预置层材料属性;材料属性2对应着熔化后的预置层材料性;材料属性3为物性参数,并定义密度和比热容为0,表面换热系数为0,热导率为正交类型,即只进行垂直方面的传热,没有横向的传热;
其中材料属性2为已知值,根据材料属性2的物性参数,材料属性1的物性参数数值通过以下公式折算得到:
式中:X1为熔化前预置层材料的物性参数,
X2为预置层熔化后的物性参数,
φ是压片的气孔率;
对于基材,只定义一种材料属性,材料属性编号为4,基材熔化前后对应同一种物性参数。
3.根据权利要求1所述的高精度预置式激光熔覆涂层温度场仿真方法,其特征是所述的热源加载及温度场求解过程为:
一次加载热源,设定一个热源加载的时间,进行一次温度场求解;
一次求解结束后,进行材料属性改变和接触热阻改变的判断,材料属性改变过程:判断冶金化仿真模型中单元的节点是否全部达到预置层的熔点,若达到,则单元材料属性发生改变,否则不发生改变;若与气孔仿真模型单元相连接的冶金化仿真模型单元已改变材料属性,则气孔仿真模型单元也改变材料属性;接触热阻改变过程:若预置层或基材的接触表面上的原死单元中全部节点温度达到预置层或基材的熔点,原死单元改变成活单元,以同样的条件把原活单元变成死单;。
再次载热源,设定一个热源加载时间,进行温度场求解;
以上热源加载,温度场求解,材料属性改变和接触热阻改变的判断,再次加载热源和进行温度场求解,这是一个多次循环过程,通过循环语句实现以上循环功能;
将每个时刻求解出的温度场将存入结果文件中,可以根据时间来查看某时刻的温度场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101771268A CN101824514B (zh) | 2010-05-19 | 2010-05-19 | 高精度预置式激光熔覆涂层温度场仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101771268A CN101824514B (zh) | 2010-05-19 | 2010-05-19 | 高精度预置式激光熔覆涂层温度场仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101824514A true CN101824514A (zh) | 2010-09-08 |
CN101824514B CN101824514B (zh) | 2012-07-04 |
Family
ID=42688675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101771268A Expired - Fee Related CN101824514B (zh) | 2010-05-19 | 2010-05-19 | 高精度预置式激光熔覆涂层温度场仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101824514B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105389435A (zh) * | 2015-11-12 | 2016-03-09 | 南京航空航天大学 | 激光束与粉末颗粒光固耦合过程的介观模拟方法 |
CN105868461A (zh) * | 2016-03-28 | 2016-08-17 | 中国科学院力学研究所 | 一种激光熔覆多组分合金涂层的成分计算方法 |
CN108251602A (zh) * | 2017-12-29 | 2018-07-06 | 中船海洋动力部件有限公司 | 激光加热金属部件表面得到指定温度场的方法 |
CN108251603A (zh) * | 2017-12-29 | 2018-07-06 | 中船海洋动力部件有限公司 | 激光加热活塞表面得到指定温度场的方法 |
CN109284524A (zh) * | 2018-07-19 | 2019-01-29 | 西北工业大学 | 一种创建高精度增材制造有限元模型的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006185173A (ja) * | 2004-12-27 | 2006-07-13 | Brother Ind Ltd | ヒータのシミュレーション方法、プログラム、記録媒体、及び画像形成装置 |
CN101288888A (zh) * | 2008-06-25 | 2008-10-22 | 中国科学院力学研究所 | 一种基于预应力的激光和喷丸复合成形方法 |
-
2010
- 2010-05-19 CN CN2010101771268A patent/CN101824514B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006185173A (ja) * | 2004-12-27 | 2006-07-13 | Brother Ind Ltd | ヒータのシミュレーション方法、プログラム、記録媒体、及び画像形成装置 |
CN101288888A (zh) * | 2008-06-25 | 2008-10-22 | 中国科学院力学研究所 | 一种基于预应力的激光和喷丸复合成形方法 |
Non-Patent Citations (2)
Title |
---|
《上海航天》 20020831 沈军 等 一种接触热阻的数值计算方法 第33-36页 1-3 , 第4期 * |
《信息技术》 20090831 马浩 等 基于ANSYS的压片预置式激光熔覆温度场仿真 第92-94页 1-3 第33卷, 第4期 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105389435A (zh) * | 2015-11-12 | 2016-03-09 | 南京航空航天大学 | 激光束与粉末颗粒光固耦合过程的介观模拟方法 |
CN105868461A (zh) * | 2016-03-28 | 2016-08-17 | 中国科学院力学研究所 | 一种激光熔覆多组分合金涂层的成分计算方法 |
CN105868461B (zh) * | 2016-03-28 | 2019-01-08 | 中国科学院力学研究所 | 一种激光熔覆多组分合金涂层的成分计算方法 |
CN108251602A (zh) * | 2017-12-29 | 2018-07-06 | 中船海洋动力部件有限公司 | 激光加热金属部件表面得到指定温度场的方法 |
CN108251603A (zh) * | 2017-12-29 | 2018-07-06 | 中船海洋动力部件有限公司 | 激光加热活塞表面得到指定温度场的方法 |
CN108251602B (zh) * | 2017-12-29 | 2019-05-17 | 中船动力研究院有限公司 | 激光加热金属部件表面得到指定温度场的方法 |
CN108251603B (zh) * | 2017-12-29 | 2019-05-17 | 中船动力研究院有限公司 | 激光加热活塞表面得到指定温度场的方法 |
CN109284524A (zh) * | 2018-07-19 | 2019-01-29 | 西北工业大学 | 一种创建高精度增材制造有限元模型的方法 |
CN109284524B (zh) * | 2018-07-19 | 2022-04-08 | 西北工业大学 | 一种创建高精度增材制造有限元模型的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101824514B (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101824514B (zh) | 高精度预置式激光熔覆涂层温度场仿真方法 | |
Hou et al. | Finite element analysis for the mechanical features of resistance spot welding process | |
Gao et al. | Numerical model of a thermoelectric generator with compact plate-fin heat exchanger for high temperature PEM fuel cell exhaust heat recovery | |
CN104182568B (zh) | 一种基于ansys有限元热分析的芯片温度预测方法 | |
CN104899351B (zh) | 基于离散元模拟的沥青混合料试件温度场分析方法 | |
Sun et al. | Numerical simulation of thermal insulation and longevity performance in new lightweight ladle | |
CN101153837A (zh) | 润湿性能动态测试方法与测试系统 | |
CN112182921A (zh) | 一种高性能合金钢选区激光熔化热力耦合行为预测方法 | |
CN103440355A (zh) | Invar钢板多层多道焊接的变形模拟方法 | |
CN101221428A (zh) | 一种步进式板坯加热炉在线三维温度场数模系统 | |
CN104985298A (zh) | 一种预测旋转电弧低合金结构钢小角度焊接温度场的方法 | |
CN106383928B (zh) | 一种粉末烧结材料高温塑性流变致密化本构模型构建方法 | |
CN111745276A (zh) | 一种复杂钣金构件电阻点焊残余应力和变形的快速模拟方法 | |
Quan et al. | A study on formation process of secondary upsetting defect in electric upsetting and optimization of processing parameters based on multi-field coupling FEM | |
CN106874591A (zh) | 一种方坯加热过程温度分布的计算方法 | |
CN101840447B (zh) | 预测旋锻过程锻压力的有限元建模方法 | |
Huang et al. | Effects of surface morphology on thermal contact resistance | |
Zhang et al. | Analytical solution of the melting process of phase-change materials in thermal energy storage system | |
Liu et al. | Thermal characteristics of Li‐ion battery based on phase change material‐aluminum plate‐fin composite heat dissipation | |
Tan et al. | Effect of Flow Channel Cross Section Shape on the Performance of Heat Storage Structure | |
Yu et al. | Dynamic evolution simulation of sintering pool in selective laser sintering walnut shell/Co-PES composite | |
Si et al. | Development of simulation program for temperature field of mass concrete structures | |
Shi et al. | Energy savings through thermally efficient crucible technology: fundamentals, process modeling, and applications | |
Li et al. | Prediction and control of weld width of deep drawing steel laser tailored blanks | |
CN103065001A (zh) | 一种基于非稳态cfd模拟的烘房工艺设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120704 Termination date: 20180519 |