| 查看: 1029 | 回复: 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 这句话好像没有起到作用，求各位大神的帮忙！！！  |
回复此楼» 猜你喜欢
基金申报
已经有4人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有9人回复
-
溴的反应液脱色
已经有6人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
推荐一本书
已经有10人回复
-
参与限项
已经有5人回复
-
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
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