| 查看: 409 | 回复: 3 | |||
fanerye2015新虫 (初入文坛)
|
[求助]
DEFINE_ADJUST编译总是报错
|
| 我的模型里面有气液两相，用DEFINE_ADJUST调整单元材料属性，将在一定范围内的单元赋予液体属性，代码如下 |
回复此楼» 猜你喜欢
302求调剂一志愿华中师范大学
已经有8人回复
-
338求调剂
已经有6人回复
-
化学调剂
已经有12人回复
-
342求调剂
已经有5人回复
-
0835学硕299求调剂 08大类可接受
已经有4人回复
-
085602调剂 初试总分335
已经有3人回复
-
288环境专硕,求调材料方向
已经有5人回复
-
申博
已经有5人回复
-
081200-11408-276学硕求调剂
已经有5人回复
-
085600调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
