[交流] 我的高斯热源加上去就发散，半椭球的热源就能用，有大神指导一下是怎么回事吗？谢谢

我的高斯热源加上去就发散，半椭球的热源就能用，有大神指导一下是怎么回事吗？谢谢

回复此楼

» 猜你喜欢

【博士招生】太原理工大学2026化工博士已经有5人回复
什么是人一生最重要的？已经有5人回复
280求调剂已经有3人回复
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？已经有8人回复
面上可以超过30页吧？已经有11人回复
版面费该交吗已经有15人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有18人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有10人回复

» 本主题相关价值贴推荐，对您同样有帮助:

2种焊接移动高斯热源的程序，有什么区别吗已经有7人回复
【求助】关于高斯热源峰值功率随时间变化的问题已经有7人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2013-09-12 10:11:59

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

dingmanboy

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 377.5
帖子: 6
在线: 13.3小时
虫号: 1968001

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

你啥都没给出来，怎么给你看啊

赞一下

回复此楼

3楼2013-09-13 10:48:16

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

851377749

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 1857.4
帖子: 94
在线: 19.8小时
虫号: 2468260

引用回帖:

3楼: Originally posted by dingmanboy at 2013-09-13 10:48:16
你啥都没给出来，怎么给你看啊

#include "udf.h"
DEFINE_SOURCE(heat,c,t,dS,eqn)
{
real source,r0,Q,a,v,time,m,n,p,H,r;
real x[ND_ND];
r0=3e-4;
H=1.5e-3;
Q=1100;
a=0.56;
v=8e-3;
time=CURRENT_TIME;
C_CENTROID(x,c,t);
m=x[0];
n=x[1];
p=x[2];
r=sqrt(m*m+(p-v*time)*(p-v*time));
source=9*a*Q/(M_PI*r0*r0*H*(1-exp(-3)))*exp(9*r*r/(r0*r0*log((fabs(n)/H))));
dS[eqn]=0;
return source;
}
这是我的热源，等会发一下cas

赞一下

回复此楼

4楼2013-09-13 17:15:19

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