| 查看: 458 | 回复: 3 | |||
fanerye2015新虫 (初入文坛)
|
[求助]
DEFINE_ADJUST编译总是报错
|
| 我的模型里面有气液两相，用DEFINE_ADJUST调整单元材料属性，将在一定范围内的单元赋予液体属性，代码如下 |
回复此楼» 猜你喜欢
让我中一个面上吧！
已经有5人回复
-
科研孤儿太难了
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
