马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我把代码贴上去,加移动热流率只能加在初始的位置,但是我加的是移动的,并且我还没加对流,加了对流就更不对,红字部分就是我要加的热源,但是就是不移动,望大虾看一下
/prep7
et,1,solid70
keyopt,1,1,0
et,2,mesh200
keyopt,2,1,6
mp,dens,1,7800
mp,c,1,502.0
mptemp,1,20,200,400,600
mpdata,ALPX,1,1,1.23e-5,1.27e-5,1.37e-5,1.45e-5
mpdata,ex,1,1,2.07e11,2.00e11,1.90e11,1.05e11
mpdata,nuxy,1,1,0.3,0.3,0.3,0.3
mpdata,kxx,1,1,130,135,140,145
rectng,0.005,0.050,0,0.030
rectng,0,0.003,0,0.110
rectng,0.003,0.005,0,0.110
aglue,1,2,3
allsel
lsel,s,loc,y,0.110
lesize,all,,,1
allsel
lsel,s,loc,y,0
lsel,r,loc,x,0,0.005
lesize,all,,,1
allsel
lsel,s,loc,x,0
lsel,a,loc,x,0.003
lesize,all,,,11
allsel
lsel,s,loc,x,0.005
lsel,u,loc,y,0,0.030
lesize,all,,,8
allsel
lsel,s,loc,x,0.005
lsel,r,loc,y,0,0.040
lesize,all,,,3
lsel,s,loc,y,0.030
lesize,all,,,3
allsel
lsel,s,loc,y,0
lsel,r,loc,x,0.005,0.050
lesize,all,,,3
allsel
lsel,s,loc,x,0.050
lesize,all,,,3
amesh,all
allsel
esize,,10
type,2
vext,all,,,,,0.150
allsel
asel,s,loc,x,0
aclear,all
k=0
*do,i,262,271,1
esel,s,elem,,i+50
*do,j,1,2
esel,a,elem,,i+60+(j-1)*10
*enddo
k=k+1
cm,element%k%,elem
allsel
*enddo
cmgrp,egrope,element1,element2,element3,element4,element5,element6,element7,element8
cmedit,egrope,add,element9,element10
/solu
antype,trans
nropt,full,,on
TROPT,FULL
esel,s,live
ALLSEL
outres,all,all
nsel,s,loc,y,0
! d,all,ux,0
! d,all,uy,0
! d,all,uz,0
d,all,temp,21
allsel
ic,all,temp,21
toffst,273
t=0.1
*do,i,0.1,1,0.1
k=i*10
cmedit,egrope,dele,element%k%
nsel,s,loc,x,0.003,0.005
nsel,r,loc,y,0.080,0.110
nsel,r,loc,z,(k-1)*0.015,k*0.015
bf,all,hgen,3765600000/12
! f,all,fx,200/12
allsel
time,i
kbc,1
allsel
deltim,i/15
solve
ekill,element%k%
allsel
fdele,all,all
TIME,i+0.001
SOLVE
*enddo
finish
/post1
esel,s,live |