| 查看: 1582 | 回复: 6 | |||
[交流]
我的高斯热源加上去就发散，半椭球的热源就能用，有大神指导一下是怎么回事吗？谢谢
|
|||
| 我的高斯热源加上去就发散，半椭球的热源就能用，有大神指导一下是怎么回事吗？谢谢 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有21人回复
-
越来越担心硕士毕业后没有竞争力了
已经有14人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有3人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
闲聊
已经有6人回复
-
申博求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 2种焊接移动高斯热源的程序，有什么区别吗
已经有7人回复
- 【求助】关于高斯热源峰值功率随时间变化的问题
已经有7人回复
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
|
#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
6楼2020-05-25 19:48:04
7楼2020-05-25 19:51:27
简单回复
天涯20122楼
2013-09-12 10:37
回复
