| 查看: 1559 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
我的高斯热源加上去就发散，半椭球的热源就能用，有大神指导一下是怎么回事吗？谢谢
|
|||
| 我的高斯热源加上去就发散，半椭球的热源就能用，有大神指导一下是怎么回事吗？谢谢 |
回复此楼» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有15人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 2种焊接移动高斯热源的程序，有什么区别吗
已经有7人回复
- 【求助】关于高斯热源峰值功率随时间变化的问题
已经有7人回复
|
#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; } 高斯热源 |
5楼2013-09-13 17:22:52
3楼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
6楼2020-05-25 19:48:04