| 查看: 2005 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
super雪儿新虫 (正式写手)
|
[求助]
尴尬做一次伸手党谁有高斯热源模型的UDF编程的程序啊，跪求
|
|
|
如题，一点C++基础都没有，我现在在学，但是希望能先有个高斯热源模型能加载上，哪位大神有的话可不可以让我用一下，谢谢了啊 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有5人回复
-
重要来源：本周末出结果
已经有10人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有11人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
-
Filecode 又变了，巨变
已经有4人回复
-
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
FileCode能看出啥？
已经有26人回复
-
关于Filecode分析方法
已经有11人回复
malasong9471
木虫 (正式写手)
- 应助: 58 (初中生)
- 金币: 7303
- 红花: 29
- 帖子: 639
- 在线: 117.1小时
- 虫号: 4100890
- 注册: 2015-09-25
- 性别: GG
- 专业: 工程热力学
4楼2016-11-14 22:51:54
super雪儿
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1736.4
- 散金: 928
- 红花: 3
- 帖子: 833
- 在线: 131.1小时
- 虫号: 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
super雪儿
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1736.4
- 散金: 928
- 红花: 3
- 帖子: 833
- 在线: 131.1小时
- 虫号: 3513169
- 注册: 2014-11-01
- 专业: 金属功能材料
5楼2016-11-15 06:58:53
