24小时热门版块排行榜

>论坛更新日志 (3538)
>虫友互识 (681)
>导师招生 (131)
>论文投稿 (124)
>文献求助 (108)
>考博 (70)
>休闲灌水 (70)
>硕博家园 (65)
>博后之家 (57)
>公派出国 (47)
>招聘信息布告栏 (44)
>论文道贺祈福 (37)
>基金申请 (36)
>考研 (34)
>教师之家 (25)
>找工作 (23)

返回列表

wdwangcool

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 625.7
散金: 74
帖子: 181
在线: 100.1小时
虫号: 2649291
注册: 2013-09-12
性别: GG
专业: 金属材料的力学行为

[求助] 高斯出现对称热源，求助已有1人参与

兄弟姐妹们看一下这两幅图，我做的是高斯热源移动，但是为什么会在相对的两个面出现对称的热源呢？我只想要一个面有热源。请问该怎么做呢？
下面是命令流
finish
/clear
/filname,weldtem
/title,welding tempserature field simluation
/unit,si
/prep7
et,1,solid70
keyopt,1,2,1
et,2,surf152
r,2
keyopt,2,4,1
keyopt,2,5,0
keyopt,2,8,3
mptemp,1,0,400,600,800,1000,1500
mpdata,kxx,1,1,18.9,28.35,31.08,32.76,31.5,28.2
mpdata,c,1,1,580,756,840,882,756,797
mpdata,hf,1,1,9.66,59.43,84.32,109.2,180.6,359
mpdata,dens,1,1,7800,7550,7550,7550,7550,7550
pi=3.1415926
I=160
u=23
efficient=0.78
r=0.0045
Qm=3*i*u*efficient/(pi*r*2)
kk=0.85
Qmax=Qm*kk
ly=0.12
lx=0.4
lz=0.12
lsize=0.01
v=0.00185
tm=lx/v
tinc=lsize/v
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=1+ly/lsize
/view,1,1,1
block,0,lx,0,ly,0,lz
save
esize,lsize
mshkey,1
vatt,1,1,1
vmesh,all
type,2
asel,s,,,2
!asel,u,,,3
amesh,all
/pnum,line,off
/pnum,defa
eplot
finish

*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)*lsize
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*lsize
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0.03
distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=Qmax*exp(-3*kk*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
/solu
antype,trans
tunif,15
kbc,1
timint,on
autots,off
asel,s,area,,,2
nsla,s,all
sf,all,conv,-1
allsel
*do,t,1e-6,tm,tinc
time,t
autots,0
deltim,,,,1
kbc,0
outres,all,all
tsres,erase

sf,all,hflux,%flux2%
solve
*enddo
finish

1.png