我做的三维焊接温度场模拟例子，可以运行

Strosailor

感觉温度场分布不对，似乎不怎么传热，欢迎大家交流。

/prep7
/units,si
toffst,273
*set,rb,0.001
length=0.1
width=0.05
height=0.005
et,1,solid70
et,2,solid90
et,3,solid87
mptemp,1,20,250,500,750,1000,1500,1700,2500
mpdata,kxx,1,1,50,47,40,27,30,35,40,300
mpdata,c,1,1,460,480,530,675,670,660,780,820
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,hf,1,1,24,46,73,140,125,224,400,450,500,500
block,0,0.1,0,0.05,-0.003,0
block,0,0.1,0,0.004,-0.003,0
block,0,0.1,0,0.01,-0.003,0  
vovlap, all                     
numcmp,node
numcmp,area
numcmp,volume

Lesize,4,rb  ！定义焊缝区单元尺寸
Lesize,5,rb
Lesize,9,rb
Lesize,10,rb
Lesize,14,rb
Lesize,19,rb
Lesize,23,rb
Lesize,24,rb
Lesize,37,rb
Lesize,38,rb
Lesize,39,rb
Lesize,40,rb

Lesize,45,rb  ！定义过渡区单元尺寸
Lesize,46,rb
Lesize,47,rb
Lesize,48,rb

Lesize,26,3*rb   ！定义远缝区单元尺寸
Lesize,35,3*rb
Lesize,31,3*rb
Lesize,36,3*rb
Lesize,41,3*rb  
Lesize,42,3*rb
Lesize,43,3*rb
Lesize,44,3*rb
Lesize,2,3*rb
Lesize,7,3*rb
Lesize,11,3*rb
Lesize,12,3*rb
Save

type,1  ！划分网格
mat,1
mshkey,1
mshape,0
vmesh,1
vmesh,2

type,2
mat,1
mshape,1,3D
mopt,pyra,1
mshkey,0
mshape,1
vmesh,3
mopt,split,2
vimp,3,1,3
save

/solu

v=0.0045
tm=length/v
tinc=rb/v
U=25
I=100
ef=0.75  
p=25*100*ef
r=0.003
qmax=3*p/(3.14159*r**2)
max_time=length/rb+1  
max_x=length/rb+1  
max_y=width/rb+1
  
*dim,flux2,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux2(i,0,k)=(i-1)*rb
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*rb
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
a=(i-1)*rb-xcenter
b=(j-1)*rb-ycenter
distance=sqrt(a**2+b**2)
*if,distance,le,r,then
flux2(i,j,k)=qmax/exp(3*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo

/solu

antype,trans
ic,all,temp,20
*do,t,1e-6,tm,tinc
time,t
autos,0
deltim,,,,1
kbc,1
nsel,s,loc,x,0,0.1   ！选择要加载热源的节点
nsel,r,loc,y,0,0.05
nsel,r,loc,z,0,0
sf,all,hflux,%flux2%

FLST,5,3,5,ORDE,3  !选择要加载对流载荷的表面
FITEM,5,1   
FITEM,5,3   
FITEM,5,-4  
ASEL,U, , ,P51X
ALLSEL,BELOW,AREA   
NSLA,S,1
sf,all,conv,-1,20 !在节点上施加对流载荷(与空气的对流换热系数）
allsel,below,volu
timint,on                                        ！打开时间积分
tintp,,,,1                                        ！将THETA设定为1
outres,all,all                                        ！输入每个子步的结果
solve
*enddo
finish

shen2mzy

正好需要，十分感谢

eli-lu

我也遇到了这样的问题，中心温度场过高

fangdian

谢谢楼主。下了

acailaogu

ansys的非线性分析能力够呛，热应力够呛，转向msc。marc了

tangxiushan

参考下，谢谢分享

crruose

问题解决没  原因在呢

shezhang

命令流有什么用啊

chutianyao

是什么原因？

wxpand2004

我做完也感觉不对，不知道楼主解决问题了没！

wyksuda

学习中，谢谢楼主。不过可惜，暂时还帮不上忙

守护蓝天

有没有做激光切削的啊？激光切削和焊接有什么区别呀？

niulian614

学习中，谢谢Lz

kakaroteqiqi

来看看哈！！

181809216

楼主解决没啊

200802172

可能是热源施加的问题

cumt410

谢谢啦，正好需要瞬态传热的东西

守护蓝天

17# cumt410 你好，能告诉下不同温度下材料的密度如何查吗？

mmmm123

谢谢楼主。

永恒

感谢楼主，我也在做这方面的。