| 查看: 1052 | 回复: 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 这句话好像没有起到作用，求各位大神的帮忙！！！  |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有4人回复
7楼2016-01-18 10:34:27
2楼2016-01-09 09:44:23
bennett_14
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 7204.8
- 红花: 4
- 沙发: 1
- 帖子: 932
- 在线: 421.5小时
- 虫号: 1073877
- 注册: 2010-08-13
- 性别: GG
- 专业: 计算机科学
6楼2016-01-09 17:40:20
8楼2016-03-06 16:29:31