| 查看: 1887 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
super雪儿新虫 (正式写手)
|
[求助]
尴尬做一次伸手党谁有高斯热源模型的UDF编程的程序啊，跪求
|
|
|
如题，一点C++基础都没有，我现在在学，但是希望能先有个高斯热源模型能加载上，哪位大神有的话可不可以让我用一下，谢谢了啊 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有11人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
疑惑？
已经有5人回复
super雪儿
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1736.4
- 散金: 928
- 红花: 3
- 帖子: 833
- 在线: 131小时
- 虫号: 3513169
- 注册: 2014-11-01
- 专业: 金属材料的微观结构
7楼2016-11-15 07:01:50
super雪儿
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1736.4
- 散金: 928
- 红花: 3
- 帖子: 833
- 在线: 131小时
- 虫号: 3513169
- 注册: 2014-11-01
- 专业: 金属材料的微观结构
2楼2016-11-14 18:21:47
|
1， # include"udf.h" DEFINE_PROFILE(heat_flux,t,i) { real xx[ND_ND]; real a, pi, r, ef,p, heat, x,y; face_t f; ef= 0.8; pi = 3.141596; r = 0.003; p=1400; heat = 3 *ef*p / (pi*r*r); begin_f_loop(f,t) { x = xx[0]; y = xx[1]; if (sqrt(x*x + y*y) <= 0.003) F_PROFILE(f, t, i) = heat*exp(-3 * (x*x + y*y) / (r*r)); else F_PROFILE(f,t,i) = 0; } end_f_loop(f,t) } 2， # include "udf.h" # define v 10 # define delt .05 # define delx .5 # define I0 2*10000000 # define n 39 DEFINE_PROFILE(g_profile,t,i) { int x[ND_ND]; int j,k; double xreg; double xl[39],I[39],y[39]; face_t f; /*xl= (v*delt)+(delx/2);*/ begin_f_loop(f,t) { F_CENTROID(x,f,t); y[0]=x[1]; for(j=1;j<=n;j++) { y[j]=(j-1)*delt; /*k=j;*/ xl[j]=(j-1)*v*delt; printf("%f\n",&xl[j]); /*xreg=(y[j]-xl[j]);*/ /*double fabs(double (y[j]-xl[j]));*/ { if ((y[j]-xl[j]) < 0) xreg = -1*(y[j]-xl[j]); else xreg = (y[j]-xl[j]); } if(xreg<delx) { I[j]=I0*exp((-2*((xreg)*(xreg)))/(delx*delx)); } else { I[j]=0; } printf( "%f\n", &I[j]); } } end_f_loop(f,t) } |
3楼2016-11-14 21:01:41
malasong9471
木虫 (正式写手)
- 应助: 58 (初中生)
- 金币: 7237
- 红花: 29
- 帖子: 639
- 在线: 117.1小时
- 虫号: 4100890
- 注册: 2015-09-25
- 性别: GG
- 专业: 工程热力学
4楼2016-11-14 22:51:54