求助：关于温度场二次开发，如何提取节点的温度梯度？

sunlll1980

首先，程序如下，在计算完模型的温度场以后运行：
ttl=1150                     !共晶温度
tts=1090                     !共晶结束温度
fsc=0.75                    !临界固相率
*get,nnode,node,,num,max    !提取节点数
*dim,nodtemp,array,nnode,2000
*dim,nnumb,array,node
m=1
*do,st,10,14400,10
SET, , ,1, ,st, ,               
n=1
  *do,i,1,nnode,1              ！节点循环
  *get,ntemp,node,i,temp        ！提取节点的温度
  nodtemp(n,m)=ntemp           ！存储节点温度
  n=n+1
  *enddo
m=m+1
*enddo
!
m=2
*do,n,1,nnode,1
ntemp=nodtemp(n,m)
*if,ntemp,gt,tts,and,ntemp,lt,ttl,then
dmtemp=(ttl-ntemp)/(ttl-tts)
fs=sqrt(dmtemp)
  *if,fs,gt,fsc,then
  *get,mm,node,n,tg,sum        ！提取节点的温度梯度
  r1=(ntemp-nodtemp(n,m-1))/10
  r2=(nodtemp(n,m+1)-ntemp)/10
  r=(r1+r2)/2
  xsh=mm/sqrt(r)
    *if,xsh,lt,0.7,then      !判断是否小于临界数，小于则出现缺陷
    nnumb(u)=n                ！记录了节点号
    NSEL,S, , ,n
    ESLN,S
    /COLOR,ELEM,BLAC,       n
    *endif
  *endif
*endif
m=m+1
*enddo
eplot
finish

问题:系统提示数据库中没有温度梯度存储,是不时我的提取节点温度梯度语句有问题,或者在ansys温度计算过程中要加什么选项来存储温度梯度,请各位帮忙解答!谢谢先!!

sunlll1980

等待帮助

scqxaie

*get,par,Node,n,TG,X,
*get,par,Node,n,TG,Y,
*get,par,Node,n,TG,Z,

lk_dragon

有限元中梯度应该值单元内的把，如果是节点的梯度值，是否需要进行每个单元求解出来后进行后处理？

hai00011

期待...