24小时热门版块排行榜    

查看: 935  |  回复: 5

wang12349

新虫 (小有名气)

[求助] 求大神看看内热源的程序有什么问题,算出来的温度场偏小50%左右?

我是做电热耦合的,物料通电产生内热源,内热源强度q=电场强度的平方/电阻率/2000,后面的2000是物料的比热容,不知道需不需要加?还有我感觉if语句没有用上呢?谢谢您了、我用ANSYS电热耦合模块做对比,发现温度场偏小50%左右。
#include "stdio.h"
#include "udf.h"

DEFINE_SOURCE(v_source, c, t, dS, eqn)
{
        real row_cal;
        real source;
        real r;
        real tt;
        tt=C_T(c,t);
        row_cal=C_R(c,t);
       
    if(tt<1073.0)
    {
        r=0.03;
        }
        if(tt>=1073.0&&tt<1473.0)
    {
        r=0.03-(tt-1073.0)/400.0*(0.03-0.024129);
        }
        if(tt>=1473.0&&tt<1773.0)
    {
        r=0.024129-(tt-1473.0)/300.0*(0.024129-0.019726);
        }
        if(tt>=1773.0&&tt<2073.0)
    {
        r=0.019726-(tt-1773.0)/300.0*(0.019726-0.015322);
        }
        if(tt>=2073.0&&tt<2473.0)
    {
        r=0.015322-(tt-2073.0)/400.0*(0.015322-0.009451);
        }
        if(tt>=2473.0)
    {
        r=0.009451;
        }
        source=row_cal*NV_MAG2(C_UDSI_G(c,t,0))/r/2000;
        dS[eqn]=0.0;
        return source;
}@wuming524
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingfuww

专家顾问 (正式写手)

后面的if 可以用else if 来。
2楼2016-08-08 09:29:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wang12349

新虫 (小有名气)

引用回帖:
2楼: Originally posted by xingfuww at 2016-08-08 09:29:05
后面的if 可以用else if 来。

谢谢,我再试算一下,这个内热源编程没问题吧。现在温度场分布的形状是对的,就是大小不对。
3楼2016-08-08 19:48:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wang12349

新虫 (小有名气)

引用回帖:
2楼: Originally posted by xingfuww at 2016-08-08 09:29:05
后面的if 可以用else if 来。

我昨天又试算了一下,发现不是if语句错误,还是内热源的大小错了
4楼2016-08-09 09:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingfuww

专家顾问 (正式写手)

引用回帖:
4楼: Originally posted by wang12349 at 2016-08-09 09:56:52
我昨天又试算了一下,发现不是if语句错误,还是内热源的大小错了...

................................
5楼2016-08-09 13:34:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wang12349

新虫 (小有名气)

内热源错误了,改成就对了 source=NV_MAG2(C_UDSI_G(c,t,0))/r。另外有用fluent做电热耦合的文献吗,求推荐
6楼2016-08-11 11:09:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wang12349 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +4 南迦720 2026-02-28 4/200 2026-02-28 16:31 by etapple
[考研] 0856材料求调剂 +8 hyf hyf hyf 2026-02-28 9/450 2026-02-28 16:27 by etapple
[考研] 285求调剂 +3 满头大汗的学生 2026-02-28 3/150 2026-02-28 16:22 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 312求调剂 +4 吃宵夜1 2026-02-28 5/250 2026-02-28 16:11 by gjm133
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 材料类求调剂 +3 wana_kiko 2026-02-28 3/150 2026-02-28 15:03 by lature00
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 290求调剂 +4 材料专硕调剂; 2026-02-28 5/250 2026-02-28 13:32 by houyaoxu
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[考研] 280求调剂 +6 Qq206./ 2026-02-21 6/300 2026-02-28 11:20 by 我!要一战成硕
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 10/500 2026-02-27 13:22 by ichall
[基金申请] 什么是人一生最重要的? +10 瞬息宇宙 2026-02-21 10/500 2026-02-27 08:46 by tfang
[教师之家] 版面费该交吗 +14 苹果在哪里 2026-02-22 17/850 2026-02-26 11:55 by AGI智创机器人
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见