| 查看: 3105 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
用UDF定义电池的热源，以下是电池发热表达式以及程序，结果发散了，问题在哪？
|
||
|
#include "udf.h" DEFINE_SOURCE(heat_gen,c,t,dS,eqn) { real i=15.; real c1=15.; real F=96485.; real source; real t1=RP_Get_Real("flow-time" ); real T=C_T(c,t); real soc=1-i*t1/c1; if(0<=soc<=0.77) { if(293<=T<303) source=2.258e-6*pow(soc,-0.3952)*i*i-i*T*(99.88*soc-76.67)/F; else if(303<=T<313) source=1.857e-6*pow(soc,-0.2787)*i*i-i*T*(99.88*soc-76.67)/F; else source=1.659e-6*pow(soc,-0.1692)*i*i-i*T*(99.88*soc-76.67)/F; } else if(0.77<soc<=0.87) { if(293<=T<303) source=2.258e-6*pow(soc,-0.3952)*i*i+i*T*30/F; else if(303<=T<=313) source=1.857e-6*pow(soc,-0.2787)*i*i+i*T*30/F; else source=1.659e-6*pow(soc,-0.1682)*i*i+i*T*30/F; } else { if(293<=T<303) source=2.258e-6*pow(soc,-0.3952)*i*i+i*T*20/F; else if(303<=T<313) source=1.857e-6*pow(soc,-0.2787)*i*i+i*T*20/F; else source=1.659e-6*pow(soc,-0.1692)*i*i+i*T*20/F; } dS[eqn]=0; return source; }  ~]M@IE8{WRB79@99SPKIW_T.png  5O]L2TN5DYKXJFH733W}LPD.png  )}6FJ5$8DTIZSDUCX9GP({R.png  {8([`H`A{52%JV5I~K]0Q3G.png |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有63人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
4楼2018-05-27 10:27:14
2楼2018-05-23 17:17:59
MoonLee777
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2274.6
- 红花: 17
- 帖子: 839
- 在线: 78小时
- 虫号: 5947368
- 注册: 2017-03-12
- 性别: GG
- 专业: 流体力学
3楼2018-05-26 00:36:31
MoonLee777
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2274.6
- 红花: 17
- 帖子: 839
- 在线: 78小时
- 虫号: 5947368
- 注册: 2017-03-12
- 性别: GG
- 专业: 流体力学
5楼2018-05-27 11:37:15
