24小时热门版块排行榜    

查看: 3195  |  回复: 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的回帖

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的回帖
查看全部 9 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:38 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 05:50 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:26 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:01 by mibUvS8DDCwf
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:21 by mibUvS8DDCwf
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见