求助:关于温度场二次开发,如何提取节点的温度梯度?
程序如下: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(tt)=n !记录了节点号
tt=tt+1
NSEL,S, , ,n
ESLN,S
/COLOR,ELEM,BLAC, n
*endif
*endif
*endif
m=m+1
*enddo
eplot
finish
我是在计算完模型的温度场后运行以上程序,而且显示了第一载荷图形。
程序中的提取节点温度梯度语句,在运行的过程中提示,数据库中没有可以提取的数值,请各位帮忙,温度梯度这时候为什么不能提取,有没有好的办法?谢谢各位先!!
Re:求助:关于温度场二次开发,如何提取节点的温度梯度?
等待
页:
[1]