| 查看: 1201 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
jiangge2585铜虫 (小有名气)
|
[求助]
udf 热源 初始化
|
|
|
求大神帮我看看我的程序哪里出问题了?这是一个激光高斯面热源程序，我的意图是：激光脉宽0.008s，频率是10Hz，就是1秒激光打10次，每次持续0.008s。但是 模拟的结果却是不能初始化，初始化时就卡主了，一直停在那。求大神救我！ 这是程序： DEFINE_PROFILE(heat_flux,t,i) { real xx[ND_ND]; real y,x,r,a,time_step; int n; face_t f; begin_f_loop(f,t) { F_CENTROID(xx,f,t); y=xx[1]; x=xx[0]; r=0.001; a=0; n=0; if(a<=1) { if(a>=0.108*n&&a<=0.108*n+0.08) { if(sqrt(x*x+y*y)<=0.001) { F_PROFILE(f,t,i)=1160*0.7/(3.1415926*r*r)*exp(-(x*x+y*y)/(r*r)); } else { F_PROFILE(f,t,i)=0; } RP_Get_Real("flow-time"  ;a=RP_Get_Real("flow-time" ;n=(int)a/0.108; } else { F_PROFILE(f,t,i)=0; } } else { Error("error reading file" ;} } end_f_loop(f,t) } |
回复此楼» 猜你喜欢
求调剂推荐
已经有7人回复
-
289 分105500药学专硕求调剂(找B区学校)
已经有4人回复
-
0854求调剂
已经有13人回复
-
初试324 中药学 一志愿天中医 求调剂
已经有3人回复
-
药学求调剂
已经有14人回复
-
327求调剂
已经有27人回复
-
急需调剂
已经有5人回复
-
300求调剂
已经有7人回复
-
271求调剂
已经有33人回复
-
273求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- FLUENT中编写UDF时指针变量是否要初始化
已经有4人回复
- FLUENT UDF 编译通过，初始化出错的原因？？？求助
已经有4人回复
- UDF 热源随温度变化的处理
已经有4人回复
- 请教：模拟带内热源的管道流动
已经有4人回复
- 使用UDF后，初始化错误
已经有11人回复
- UDF应用之后的初始化错误
已经有6人回复
- 【求助】fluent计算中，udf编译通过，初始化出错
已经有11人回复
perfectgem
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1138.9
- 散金: 6
- 帖子: 92
- 在线: 11.7小时
- 虫号: 2101103
- 注册: 2012-10-29
- 性别: MM
- 专业: 零件加工制造
5楼2014-01-06 15:49:31
jiangge2585
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 75.5
- 帖子: 70
- 在线: 55.4小时
- 虫号: 1508657
- 注册: 2011-11-24
- 性别: GG
- 专业: 传动机械学
|
这些天想了想，好像是begin_f_loop(f,t) 和循环程序冲突了，所以才会卡主，begin_f_loop(f,t) 扫描出一个面网格中心P,然后把符合(sqrt(x*x+y*y)<=0.001) 的P丢进循环语句执行，执行过一次之后，a和n的值就已经超出了循环条件所规定的范围，不能再进行循环了。 我就改了一下，但下面这个程序却没有发热，不知道哪有问题了？求大神救我！ #include "udf.h" DEFINE_PROFILE(heat_flux,t,i) { real xx[ND_ND]; real y,x,r,a; int n; face_t f; begin_f_loop(f,t) { F_CENTROID(xx,f,t); y=xx[1]; x=xx[0]; r=0.001; if(sqrt(x*x+y*y)<=0.001) for(a=n=0;a>=0.1*n&&a<=0.1*n+0.008&&a<=1;n=(int)a/0.1) { F_PROFILE(f,t,i)=1160*0.7/(3.1415926*r*r)*exp(-(x*x+y*y)/(r*r)); a=RP_Get_Real("flow-time" ;} } end_f_loop(f,t) } |
2楼2013-06-25 15:14:35
★
1592203609: 金币+1, 谢谢回帖，建议给出更为详细的答案 2013-06-25 17:58:02
1592203609: 金币+1, 谢谢回帖，建议给出更为详细的答案 2013-06-25 17:58:02
|
本帖内容被屏蔽 |
4楼2013-06-25 15:27:00
5