如何添加流场分析得出来的温度值？
浏览：44395 回答：6
这个问题曾经发了很多次，但至今仍未得到答案，很着急。再次恳请明白人给予指点！
做管道热应力分析，先用FLOTRAN计算管内对流换热系数从而求出管壁内的温度分布，然后将温度场施加到结构分析单元上进行热应力分析。但是温度好像只存在于rth文件中，只见过ldread,temp,,,,,,rth，而流场计算结果只有rfl文件，虽然流场计算也能得到温度分布，但温度应当如何提取出来呢？我试过ldread,temp,,,,,,rfl，是错的，真不知道这种问题应该怎么处理。
这个问题困扰了我有一个多月，也没有人可以讨论的，非常沮丧。还望有高人给我指一条明路，万分感谢！
做管道热应力分析，先用FLOTRAN计算管内对流换热系数从而求出管壁内的温度分布，然后将温度场施加到结构分析单元上进行热应力分析。但是温度好像只存在于rth文件中，只见过ldread,temp,,,,,,rth，而流场计算结果只有rfl文件，虽然流场计算也能得到温度分布，但温度应当如何提取出来呢？我试过ldread,temp,,,,,,rfl，是错的，真不知道这种问题应该怎么处理。
这个问题困扰了我有一个多月，也没有人可以讨论的，非常沮丧。还望有高人给我指一条明路，万分感谢！
Finish ！（从此处到“——”处命令无变化）
/clear
/filname,example2,1
/prep7
et,1,plane55,,,1
mp,kxx,1,1.2
rectng,0.2,0.8,0,0.05
type,1
lsel,s,line,,1,3,2
lesize,all,,,10
lsel,s,line,,2,4,2
lesize,all,,,2
amesh,1
finish
/solu
antype,static
lsel,s,line,,4
nsll,s,1
d,all,temp,80
lsel,s,line,,2
nsll,s,1
d,all,temp,10
allsel
outpr,basic,all
solve
finish
/post1
plnsol,temp
—————————————————
allsel
nsel,all
*get,n1,node,,num,min ！最小节点值
*get,n2,node,,num,max ！最大节点值
n3=n2-n1+1 ！总节点数
*dim,abc,array,n3 ！定义一个记录所有节点温度的数组
*do,i,1,n3
*get,abc(i),node,n1+i-1,temp ！给数组赋温度值
*enddo
——————————————————
/prep7 ！（从此处到“——”处命令与example1中的相同）
etchg,tts
keyopt,1,3,1
keyopt,1,6,1
mp,ex,1,120e9
mp,alpx,,1.3e-6
mp,nuxy,1,.3
lsel,s,line,,3
nsll,s,1
cp,1,uy,all
allsel
cp,2,ux,1,14,24
cp,3,ux,2,12,13
finish
/solu
antype,static
d,1,uy,,,11
——————————————————
*do,i,1,n3
bf,n1+i-1,temp,abc(i) ！给结构单元上的所有节点赋温度值
*enddo
——————————————————
Solve ！（从此处到结束的命令与example1中的相同）
finish
/post1
plnsol,s,x
nsel,s,loc,x,0.2
nsel,r,loc,y,0.05
*get,node1,node,,num,max
nsel,s,loc,x,0.8
nsel,r,loc,y,0.05
*get,node2,node,,num,max
allsel
path,spath1,2,,30
ppath,1,node1
ppath,2,node2
pdef,sx,s,x
plpath,sx
大家有兴趣的话可以分别运行一下这两段程序，会发现两种解法得出来的应力分布图是一样的。在这里我只张贴出沿矩形最上面那条边的温度分布曲线：
1、(说我图片文件无效，有兴趣的话大家自己运行下吧)
2、
我们发现它们的数值是完全一样的。这说明我使用的这种方法可行。
这仅仅是个简单的例子，用来说明我的思路，有不足或错误的地方还望大家不吝赐教!