24小时热门版块排行榜    

查看: 3099  |  回复: 8

hiqun

禁虫 (正式写手)

[求助] DEFINE_SOURCE(energy,c,t,ds,eqn)条件赋值的相关问题

最近做固体发热,采用DEFINE_SOURCE(energy,c,t,ds,eqn)写的udf函数,程序l列举如下:
#include "udf.h"
DEFINE_SOURCE(energy,c,t,ds,eqn)
{real source;
real temp;
temp=C_T(c,t);

if(temp>2000)
{source=(temp-10)*1000;
dS[eqn]=0;
}
else
{source=temp*1000;
dS[eqn]=0;
}
return source;
}

鉴于水平有限,写这个的目的,就是希望,温度T>2000K的时候,加热功率自动减小;低于2000K时,维持此加热功率;目前调试下来,加热功率是一直降低的,无法实现我要的功能,恳请帮忙指导一下,在此谢过。

[ Last edited by hiqun on 2012-11-22 at 15:53 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)



xiegangmai: 金币+1, 谢谢参与 2012-11-26 19:11:50
从程序描述的函数来看,若温度低于2000k,散热又太快,功率就是越来越小的。

[ 发自手机版 http://muchong.com/3g ]
2楼2012-11-22 23:16:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与 2012-11-24 20:49:29
hiqun: 金币+25, ★★★很有帮助, 谢谢。 2012-11-27 21:30:29
大概看了一下你的程序,你无法通过source=temp*1000保证加热功率维持,试想,如果你散热很强,温度一直降低,那么你的source值是不是就一直降低,而不是你想要的维持。
   所以,设定的source可以和T无关。
可以改为:


#DEFINE K 90
if(temp>2000)
{source=CURRENT_TIME*K*1000;
dS[eqn]=0;
}
else
{source=1000;
dS[eqn]=0;
}

其中,k是一个因子,具体可调。
   当然,这个程序是抛砖引玉,只是一种思路,将功率和时间等联系起来。低于2000一下固定值等。
3楼2012-11-23 21:03:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

上面的DEFINE应该小写。。。。。
4楼2012-11-23 21:04:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哥本哈根月光

铁杆木虫 (正式写手)

引用回帖:
3楼: Originally posted by fangdian at 2012-11-23 21:03:33
大概看了一下你的程序,你无法通过source=temp*1000保证加热功率维持,试想,如果你散热很强,温度一直降低,那么你的source值是不是就一直降低,而不是你想要的维持。
   所以,设定的source可以和T无关。
可以 ...

我有一个疑问, dS[eqn]=0.0是此UDF必须要写进去的吗?
你眉目如当年,流转我心间
5楼2012-11-25 11:47:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanxueyi

金虫 (小有名气)

【答案】应助回帖


xiegangmai: 金币+1, 谢谢参与 2012-11-26 19:11:58
dS[]=0.0必须写,但是你最后在return前写一个就够了。
苦闷的博士
6楼2012-11-25 20:23:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
6楼: Originally posted by fanxueyi at 2012-11-25 20:23:07
dS[]=0.0必须写,但是你最后在return前写一个就够了。

dS[]=0在这段程序中起啥作用呢?不写好像也能执行这段UDF啊

[ 发自手机版 http://muchong.com/3g ]
7楼2012-11-25 23:09:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

引用回帖:
7楼: Originally posted by 哥本哈根月光 at 2012-11-25 23:09:21
dS[]=0在这段程序中起啥作用呢?不写好像也能执行这段UDF啊
...

dS[]这个是和收敛速度相关的,一般不影响结果,但是影响收敛速度
8楼2012-11-26 09:32:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哥本哈根月光

铁杆木虫 (正式写手)

引用回帖:
8楼: Originally posted by fangdian at 2012-11-26 09:32:42
dS[]这个是和收敛速度相关的,一般不影响结果,但是影响收敛速度...

谢谢,学到了
你眉目如当年,流转我心间
9楼2012-11-26 13:40:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hiqun 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么网上很多人说本周 12号出结果 +4 瞬息宇宙 2026-08-10 4/200 2026-08-11 02:02 by 旁观2020
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见