| 查看: 359 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
fanerye2015新虫 (初入文坛)
|
[求助]
DEFINE_ADJUST编译总是报错
|
|
| 我的模型里面有气液两相，用DEFINE_ADJUST调整单元材料属性，将在一定范围内的单元赋予液体属性，代码如下 |
回复此楼» 猜你喜欢
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有5人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
» 本主题相关商家推荐: (我也要在这里推广)
fanerye2015
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160.5
- 散金: 7
- 红花: 1
- 帖子: 42
- 在线: 24小时
- 虫号: 3716536
- 注册: 2015-03-07
- 专业: 金属材料的凝固与结晶学
4楼2015-12-16 09:55:25
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