关于搅拌摩擦焊ANSYS模拟温度场

zhangshinan

小弟初学ANSYS，从论坛下了些命令流自己修改了下 进行搅拌摩擦焊温度场的模拟 但模拟出的温度分布云图全是同一个颜色：都是初始温度 请教各位大侠这究竟是哪出问题了 是不是 热加载有问题 谢谢 命令流如下：
finish
/clear
/FILNAME,fsw-summer
/TITLE,fsw-summer
save,fsw-summer
/PREP7
et,1,plane55
et,2,solid70
MPTEMP,1,20,50,100,200,300,400         
MPTEMP,7,500,600
MPDATA,dens,1,,2788,2781,2775,2750,2725,2700
MPDATA,dens,1,7,2675,2590  
MPDATA,c,1,,860,870,900,970,1020,1120      
MPDATA,c,1,7,1320,1400
MPDATA,kxx,1,,156,158.3,161,175,185,193   
MPDATA,kxx,1,7,197,207
!mpdata,hf,1,1,15,30,50,80,110,120,130
mpdata,hf,1,1,6,10.5,18,20,30,40  
mpdata,hf,1,7,50,147
mp,reft,1,20
/PREP7  
length=0.300  
width=0.095   
height=0.010  
r_shoulder=0.009   
r_tool=0.0045   
v_move=8.33e-4   
w_rotate=31.4  
u_fri=0.42
pu=2.6e9
ini_temp=20
w_width=r_shoulder+0.01   
effi=0.5
k,1,0,0,0
k,2,0,w_width,0
k,3,0,width,0
k,4,length,0,0
k,5,length,w_width,0
k,6,length,width,0
k,7,0,0,height
k,8,0,-w_width,0
k,9,0,-width,0
k,11,length,-width,0
k,10,length,-w_width,0
LSTR,1,7
a,1,2,5,4
a,2,3,6,5
a,1,8,10,4
a,8,9,11,10
lsize_h=0.0025  
lsize_w=0.0025
lsize_l=0.004  
nelem_h=height/lsize_h     
nelem_w=w_width/lsize_w
nelem_l=length/lsize_l        
lesize,1,,,nelem_h,,,,,1
lesize,3,,,nelem_l,,,,,1      
lesize,5,,,nelem_l,,,,,1
lesize,7,,,nelem_l,,,,,1
lesize,10,,,nelem_l,,,,,1
lesize,13,,,nelem_l,,,,,1
lesize,2,,,nelem_w,,,,,1   
lesize,4,,,nelem_w,,,,,1
lesize,9,,,nelem_w,,,,,1
lesize,11,,,nelem_w,,,,,1
lesize,6,,,10,4,,,,1
lesize,12,,,10,4,,,,1
lesize,8,,,10,0.25,,,,1
lesize,14,,,10,0.25,,,,1

lplot
amesh,1,2,1
amesh,3,4,1
type,2
mat,1     
extopt,esize,nelem_h,0      
extopt,aclear,1   
extopt,attr,1,0,0
allsel,all
VDRAG,all,,,,,,1
/VIEW,1,1,2,3   
eplot
/solu
pi=3.14159
hgen=effi*w_rotate*u_fri*pu
total=(length-0.02)/lsize_l+1  
n1_x=2*r_tool/lsize_l+1
n1_y=r_tool/lsize_w+1  
n1_z=nelem_h+1-1   
n2_x=2*r_shoulder/lsize_l+1
n2_y=r_shoulder/lsize_w+1   
!n2_z=nelem_h+1   
n2_z=1
*dim,hgen1,array,n1_x,2*n1_y,n1_z     
*do,i,1,n1_x,1  !4
*do,j,1,2*n1_y,1  !4
*do,k,1,n1_z,1  !5
x=-r_tool+(i-1)*lsize_l
y=(j-1)*lsize_w-r_tool
z=(k-1)*lsize_h
  dis=sqrt(x**2+y**2)   
  *if,dis,ge,(r_tool-0.0002),and,dis,le,(r_tool+0.0002),then  
     hgen1(i,j,k)=hgen*dis   
  *else
     hgen1(i,j,k)=0
  *endif
*enddo
*enddo
*enddo

*dim,hgen2,array,n2_x,2*n2_y,n2_z      
*do,i,1,n2_x,1  !11
*do,j,1,2*n2_y,1  !6
*do,k,1,n2_z,1  !1

  x=-r_shoulder+(i-1)*lsize_l
  y=(j-1)*lsize_w-r_shoulder
  z=(k-1)*lsize_h
  
  dis=sqrt(x**2+y**2)   
  *if,dis,le,r_shoulder,and,dis,ge,r_tool,then
     hgen2(i,j,k)=hgen*dis
  *else
     hgen2(i,j,k)=0
   *endif
*enddo
*enddo
*enddo

/solu
antype,4,new  
trnopt,full
lumpm,0
tunif,ini_temp,
thopt,full,0.05,64, , ,
nropt,full,,on  
tsres,erase
lnsrch,1   
time,1e-6
autots,1
deltim, , , ,1
kbc,1         
timint,1   
tintp,0.005, , ,1,0.5,0.2,
outres,all,
!save

allsel,all         
nsel,s,loc,x,0
sf,all,conv,-1,20
allsel,all           
nsel,s,loc,x,length
sf,all,conv,-1,20
allsel,all
nsel,s,loc,y,-width
sf,all,conv,-1,20
allsel,all
nsel,s,loc,y,width
sf,all,conv,-1,20
allsel,all
nsel,s,loc,z,0   
sf,all,conv,500,20
allsel,all
allsel,all
nsel,s,loc,z,hight
sf,all,conv,-1,20
allsel,all
     
dt=lsize_l/v_move
*do,steps,1,total,1  !total=71

  xcenter=lsize_l*(steps-1)+0.01
  
*do,i,1,n1_x,1
*do,j,1,2*n1_y,1
*do,k,1,n1_z,1
x=-r_tool+xcenter+(i-1)*lsize_l
y=(j-1)*lsize_w-r_tool
z=height-(k-1+1)*lsize_h
node_1=node(x,y,z)
bf,node_1,hgen,hgen1(i,j,k)
*enddo
*enddo
*enddo

*do,i,1,n2_x,1
*do,j,1,2*n2_y,1
*do,k,1,n2_z,1
  x=-r_shoulder+xcenter+(i-1)*lsize_l
  y=(j-1)*lsize_w-r_shoulder  
  z=height-(k-1)*lsize_h
   node_2=node(x,y,z)
   bf,node_2,hgen,hgen2(i,j,k)   
   
*enddo
*enddo
*enddo
allsel,all         
nsubst,1
time,steps*dt      
solve
bfdele,all,hgen
*enddo
save
allsel,all
bfdele,all,hgen
!sfdele,all,conv
time,3000
nsubst,5
allsel,all           
solve
save