24小时热门版块排行榜    

查看: 709  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wdwangcool

金虫 (小有名气)

[求助] 高斯出现对称热源,求助 已有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


高斯出现对称热源,求助-1
2.png
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
试试看:在施加热源的部分,除了由 xcenter 和 ycenter 计算 distance 的判断外,增加一个对于单元的 z 向中心位置的判断。
3楼2015-07-30 08:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

wdwangcool

金虫 (小有名气)

同学们帮帮忙给些提示也好啊
2楼2015-07-28 13:34:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wdwangcool

金虫 (小有名气)

引用回帖:
3楼: Originally posted by shgao20 at 2015-07-30 08:15:15
试试看:在施加热源的部分,除了由 xcenter 和 ycenter 计算 distance 的判断外,增加一个对于单元的 z 向中心位置的判断。

问题找到了 是因为多了一个壳单元,但是多了壳单元就会出现这种情况吗?
4楼2015-08-16 12:20:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +8 天下熯 2026-02-28 8/400 2026-03-02 00:15 by 暮雨星晴
[考研] 材料学硕318求调剂 +11 February_Feb 2026-03-01 13/650 2026-03-01 23:53 by ccp273206157
[考研] 一志愿郑大材料学硕298分,求调剂 +5 wsl111 2026-03-01 5/250 2026-03-01 23:45 by 暮雨星晴
[考研] 0856材料与化工,270求调剂 +6 YXCT 2026-03-01 6/300 2026-03-01 23:21 by 向上的胖东
[基金申请] 成果系统访问量大,请一小时后再尝试。---NSFC啥时候好哦,已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 材料类求调剂 +10 wana_kiko 2026-02-28 12/600 2026-03-01 22:10 by 海嵙Y
[考研] 一志愿中南大学理学化学 +4 15779376950 2026-03-01 5/250 2026-03-01 19:00 by Fff-1
[考研] 295求调剂 +7 19171856320 2026-02-28 7/350 2026-03-01 18:54 by 18137688336
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 290求调剂 +9 材料专硕调剂; 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
信息提示
请填处理意见