| 查看: 1074 | 回复: 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 这句话好像没有起到作用，求各位大神的帮忙！！！  |
回复此楼» 猜你喜欢
805600专硕材料与化工348分求调剂
已经有6人回复
-
重庆大学材料与化工085600，初试370+，求求调剂建议
已经有5人回复
-
309求调剂
已经有3人回复
-
一志愿北交大材料工程，总分358
已经有4人回复
-
302求调剂一志愿北航070300，本科郑大化学
已经有11人回复
-
化学工程专硕324分，一志愿中国矿业大学求调剂
已经有7人回复
-
材料求调剂
已经有11人回复
-
【求调剂】新能源材料本科，一志愿211，初试321
已经有4人回复
-
336求调剂
已经有5人回复
-
材料求调剂
已经有8人回复
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
