| 查看: 3102 | 回复: 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 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
时间戳又变了8-15
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有57人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
关于Filecode分析方法
已经有13人回复
MoonLee777
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2274.6
- 红花: 17
- 帖子: 839
- 在线: 78小时
- 虫号: 5947368
- 注册: 2017-03-12
- 性别: GG
- 专业: 流体力学
3楼2018-05-26 00:36:31
2楼2018-05-23 17:17:59
4楼2018-05-27 10:27:14
MoonLee777
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2274.6
- 红花: 17
- 帖子: 839
- 在线: 78小时
- 虫号: 5947368
- 注册: 2017-03-12
- 性别: GG
- 专业: 流体力学
5楼2018-05-27 11:37:15
