24小时热门版块排行榜    

查看: 2833  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +4 吃吃吃才有意义 2026-03-19 4/200 2026-03-19 13:48 by houyaoxu
[考研] 085601材料工程专硕求调剂 +8 慕寒mio 2026-03-16 8/400 2026-03-19 12:30 by peike
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
请填处理意见