| 查看: 408 | 回复: 3 | ||
fanerye2015新虫 (初入文坛)
|
[求助]
DEFINE_ADJUST编译总是报错
|
| 我的模型里面有气液两相，用DEFINE_ADJUST调整单元材料属性，将在一定范围内的单元赋予液体属性，代码如下 |
回复此楼» 猜你喜欢
081700学硕，323分，一志愿中国海洋大学求调剂学校
已经有8人回复
-
283求调剂
已经有4人回复
-
材料工程310专硕调剂
已经有3人回复
-
一志愿北京化工085600 310分求调剂
已经有10人回复
-
081700，311，求调剂
已经有6人回复
-
282求调剂
已经有5人回复
-
一志愿9材料学硕297已过六级求调剂推荐
已经有9人回复
-
278求调剂
已经有10人回复
-
301求调剂
已经有4人回复
-
085600，专业课化工原理，321分求调剂
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
fanerye2015
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160.5
- 散金: 7
- 红花: 1
- 帖子: 42
- 在线: 24小时
- 虫号: 3716536
- 注册: 2015-03-07
- 专业: 金属材料的凝固与结晶学
|
#include "udf.h" #include "sg_mphase.h" #include "math.h" #define v2 0.033333 #define r2 0.0004 DEFINE_ADJUST(powder,domain) { Thread *t; cell_t c; Thread *pri_th, *sec_th; real xx[ND_ND]; real x2,y2,z2,r; real Time2=RP_Get_Real("flow_time"  thread_loop_c(t,domain) { pri_th=THREAD_SUB_THREAD(t,0); sec_th=THREAD_SUB_THREAD(t,1); begin_c_loop_all(c,t) { C_CENTROID(xx,c,t); x2=xx[0]-0.003-v2*Time2; y2=xx[1]; z2=xx[2]-0.0005; r=sqrt(x2*x2+y2*y2+z2*z2); if(r<=r2) { C_VOF(c,sec_th)=1; C_T(c,t)=1928; } } end_c_loop_all(c,t); } } |
2楼2015-12-16 09:54:18
fanerye2015
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160.5
- 散金: 7
- 红花: 1
- 帖子: 42
- 在线: 24小时
- 虫号: 3716536
- 注册: 2015-03-07
- 专业: 金属材料的凝固与结晶学
|
#include "udf.h" /*包括常规宏*/ #include "sg_mphase.h" /*包括体积分数宏*/ #include "math.h" #define v2 0.033333 #define r2 0.0004 DEFINE_ADJUST(powder,domain) { Thread *t; cell_t c; Thread *pri_th, *sec_th; real xx[ND_ND]; real x2,y2,z2,r; real Time2=RP_Get_Real("flow_time" thread_loop_c(t,domain) { pri_th=THREAD_SUB_THREAD(t,0); sec_th=THREAD_SUB_THREAD(t,1); begin_c_loop_all(c,t) { C_CENTROID(xx,c,t); x2=xx[0]-0.003-v2*Time2; y2=xx[1]; z2=xx[2]-0.0005; r=sqrt(x2*x2+y2*y2+z2*z2); if(r<=r2) { C_VOF(c,sec_th)=1; C_T(c,t)=1928; } } end_c_loop_all(c,t); } } |
3楼2015-12-16 09:54:54
fanerye2015
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160.5
- 散金: 7
- 红花: 1
- 帖子: 42
- 在线: 24小时
- 虫号: 3716536
- 注册: 2015-03-07
- 专业: 金属材料的凝固与结晶学
4楼2015-12-16 09:55:25
