| 查看: 3133 | 回复: 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 |
回复此楼» 猜你喜欢
梦想
已经有5人回复
-
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
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
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
6楼2018-05-28 09:04:24
