| 查看: 1013 | 回复: 8 | ||
[求助]
求助关于热流密度子程序
|
|
SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS, 1 JLTYP,TEMP,PRESS,SNAME) C INCLUDE 'ABA_PARAM.INC' C DIMENSION FLUX(2), TIME(2), COORDS(3) CHARACTER*80 SNAME x=COORDS(1) y=COORDS(2) z=COORDS(3) t=TIME(1) c q功率 v车速 R半径 t0刹车时间 R=460.0 r1=315.0 r2=180.0 PI=3.1415 r=(x**2+y**2)**0.5 if (0<=t<=36.9) then v=291.5525-202.6642*exp(0.0036484*t) q=5700000000*v*(0.0036484*v-1.0637) else v=100.9948-1.11914*t q=5700000000*1.1191*v endif if (r>=r2 .and. r<=r1 ) then heat=0.8*q/PI/(r1**2-r2**2)/2 shape=r/R flux(1)=heat*shape else flux(1)=0 endif RETURN END 这是一个关于热流密度的子程序，我想要让它在180＜r＜320的半径上输出热流密度flux，热流密度与时间还有圆环上各位置的半径相关，但是输出的结果好像有问题，半径相关的：shape=r/R 这句话好像没有起到作用，求各位大神的帮忙！！！  |
回复此楼» 猜你喜欢
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
职称评审没过，求安慰
已经有22人回复
-
垃圾破二本职称评审标准
已经有15人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
2楼2016-01-09 09:44:23
3楼2016-01-09 09:51:38
4楼2016-01-09 11:10:09
|
|
5楼2016-01-09 17:26:34
bennett_14
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 7204.8
- 红花: 4
- 沙发: 1
- 帖子: 932
- 在线: 421.5小时
- 虫号: 1073877
- 注册: 2010-08-13
- 性别: GG
- 专业: 计算机科学
6楼2016-01-09 17:40:20
7楼2016-01-18 10:34:27
8楼2016-03-06 16:29:31
9楼2016-09-23 19:22:44