24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

波文99

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 8
在线: 2.2小时
虫号: 3796063
注册: 2015-04-08

[求助] 求助关于热流密度子程序

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一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有5人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有7人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有6人回复
课题组招2027级博士上海工程技术大学激光智能制造方向已经有6人回复
现代”学阀”该如何界定已经有14人回复
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？已经有4人回复
有机合成以后会不会被AI改变？做科研的虫友怎么看已经有8人回复
我的奶奶已经有3人回复

1楼 2016-01-09 09:42:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiaoxu9206

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 146.1
帖子: 26
在线: 7小时
虫号: 4541138
注册: 2016-03-25
专业: 机械工程

引用回帖:

7楼: Originally posted by 波文99 at 2016-01-18 10:34:27
是啊，C为注释，第一行为ABAQUS的接口，问题找到了，原来只是fortran不区分大小写...

你好，在你的问题的基础上能向你请教几个问题吗？谢谢
1.你这个定义热源的时候是体热源还是面热源啊？
2.因为我不知道你的模型是什么样子的，所以第一个if里面的语句不能理解，能不能给解释一下啊？
3.heat=0.8*q/PI/(r1**2-r2**2)/2这个是你的热源模型吗？是自己推倒的还是查文献得到的啊？
希望您看到了能给小弟解答一下，万分感谢