应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 新手求教UDF问题~
171/2返回列表12下一页
查看: 1652  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiaoguai3

新虫 (初入文坛)

[求助] 新手求教UDF问题~

能解释，但一添加到边界或初始化就报错。
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()

我的UDF如下，二维模型，简单的固体传热问题，想在模拟过程的同时累加整个固体区域的所有得热量，达到一定值后改变热流边界值。由于刚接触UDF不久，不知道问题出在何处。请指教！d= Get_Domain(1);这个语句有问题吗？
#include "udf.h"
real Q_tot=0;
DEFINE_PROFILE (unsteady_heatflux, thread, position)
{
Domain *d;
real cp=880;
real density=2180;
real tempn, templ, volume, Q;
Thread *t;
cell_t c;
d= Get_Domain(1); /*Get the domain using Fluent utility *//*Loop over all cell threads in the domain*/
thread_loop_c(t,d) /*Compute Q */ /*Loop over all cells */
begin_c_loop(c,t)
{
real time= RP_Get_Real("flow-time";
real b=(int)(time/3600)+1;
int i=(int)((b-1)/24);
real Heat[3]={278.5, 280.7, 278.5};
volume=C_VOLUME(c,t); /* get cell volume */
templ=C_T_M1(c,t); /*Get cell tempertuare of previous step*/
tempn=C_T(c,t); /*Get cell tempertuare*/
Q=cp*density*volume*(tempn-templ);
Q_tot +=Q;
if (Q_tot<= Heat )
F_PROFILE(c,t,position)=60;
else F_PROFILE(c,t,position)=0;
}
end_c_loop(c,t)
}
回复此楼

» 本帖@通知

@wuming524

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-01-04 20:57:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuming524

专家顾问 (著名写手)

石化蛀虫

【答案】应助回帖

★
感谢参与，应助指数 +1
xiegangmai(金币+1): 谢谢参与！ 2012-01-04 22:45:44
楼主太抬举俺了，出现这种问题你需要检查你的语句，大多出现调用无效语句出现的错误。本人也不是太懂。待高人解答吧~好运。
赞一下(1人) 回复此楼
2楼2012-01-04 22:06:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
xiegangmai(金币+1): 鼓励讨论交流！ 2012-01-04 22:45:54
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？
赞一下(1人) 回复此楼
3楼2012-01-04 22:06:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoguai3

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by kylafree at 2012-01-04 22:06:24:
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？

不好意思打错了，是if (Q_tot <= Heat )
Q_tot是得热之和，由下求得
Q=cp*density*volume*(tempn-templ);
Q_tot +=Q;
赞一下 回复此楼
4楼2012-01-04 22:15:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoguai3

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by kylafree at 2012-01-04 22:06:24:
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？

为什么帖子中我数组后头的方括号就是打不出来？？？
赞一下 回复此楼
5楼2012-01-04 22:16:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主
★
臭水沟(金币+1): 谢谢交流~~ 2012-01-09 11:21:29
引用回帖:
5楼: Originally posted by xiaoguai3 at 2012-01-04 22:16:16:
为什么帖子中我数组后头的方括号就是打不出来？？？

请使用插入代码的功能
赞一下(1人) 回复此楼
明德厚学、求是创新
6楼2012-01-04 22:46:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoguai3

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by kylafree at 2012-01-04 22:06:24:
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？

不好意思，是这样的
CODE:
#include "udf.h"
real Q_tot=0;
DEFINE_PROFILE (unsteady_heatflux, thread, position)
{
Domain *d;
real cp=880;
real density=2180;
real tempn, templ, volume, Q;
Thread *t;
cell_t c;
d= Get_Domain(1); /*Get the domain using Fluent utility *//*Loop over all cell threads in the domain*/
thread_loop_c(t,d) /*Compute Q */ /*Loop over all cells */
begin_c_loop(c,t)
{
real time= RP_Get_Real("flow-time";
real b=(int)(time/3600)+1;
int i=(int)((b-1)/24);
real Heat[3]={278.5, 280.7, 278.5};
volume=C_VOLUME(c,t); /* get cell volume */
templ=C_T_M1(c,t); /*Get cell tempertuare of previous step*/
tempn=C_T(c,t); /*Get cell tempertuare*/
Q=cp*density*volume*(tempn-templ);
Q_tot +=Q;
if (Q_tot<= Heat[i] )
F_PROFILE(c,t,position)=60;
else F_PROFILE(c,t,position)=0;
}
end_c_loop(c,t)
}
赞一下 回复此楼
7楼2012-01-04 22:52:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoguai3

新虫 (初入文坛)
引用回帖:
楼: Originally posted by xiegangmai at 2012-01-04 22:46:29:
请使用插入代码的功能

谢谢版主
回复此楼
8楼2012-01-04 22:53:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoguai3

新虫 (初入文坛)
引用回帖:
楼: Originally posted by wuming524 at 2012-01-04 22:06:19:
楼主太抬举俺了，出现这种问题你需要检查你的语句，大多出现调用无效语句出现的错误。本人也不是太懂。待高人解答吧~好运。

还是谢谢~
回复此楼
9楼2012-01-04 22:53:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoguai3

新虫 (初入文坛)
引用回帖:
楼: Originally posted by kylafree at 2012-01-04 22:06:24:
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？

你好，我改了下CODE，还是出错，能帮我看下吗？谢谢了~
http://muchong.com/bbs/viewthread.php?tid=4026752&pid=2&page=1#pid2
赞一下 回复此楼
10楼2012-01-08 15:30:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoguai3 的主题更新
171/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕318求调剂 +9 February_Feb 2026-03-01 11/550 2026-03-01 19:47 by 无懈可击111
[考研] 265分求调剂不调专业和学校有行学上就 +5 礼堂丁真258 2026-02-28 7/350 2026-03-01 19:12 by Js512888
[考研] 一志愿中南大学理学化学 +4 15779376950 2026-03-01 5/250 2026-03-01 19:00 by Fff-1
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 291分工科求调剂 +9 science饿饿 2026-03-01 10/500 2026-03-01 18:55 by 18137688336
[考研] 0857调剂 +3 一ll半 2026-02-28 3/150 2026-03-01 18:32 by 热情沙漠
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 7/350 2026-03-01 17:20 by 2026计算机_诚心
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[基金申请] 刚录用，没有期刊号，但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[考研] 材料工程274求调剂 +3 Lilithan 2026-03-01 3/150 2026-03-01 14:58 by ms629
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 284求调剂 +6 天下熯 2026-02-28 6/300 2026-03-01 14:19 by Ducount.Y
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 272求调剂 +4 田智友 2026-02-28 4/200 2026-03-01 06:43 by 刘兵
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
信息提示
关闭
请填处理意见
关闭