关于搅拌摩擦焊ANSYS模拟温度场
小弟初学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
页:
[1]