beam189 单元施加完单元内外温差之后，怎样查看温度场产生的轴力弯矩剪力？

好不容易把温差加载上去了，单独对温度场求解之后发现梁单元内力无结果，求大神赐教！！ 附上求解部分代码

alls

lsel,s,loc,y,-tol,tol

lsel,u,loc,x,x2-tol,x3+tol

esll,s

eplot

!yPreApp,-pyz,-pyz   !压重层

tempApp,td,10,0  ! temApp,截面高度，外表面温度，内表面温度 （适用于BAAM189）

alls

!acel,,9.8

/SOLU

!antype,0

SOLVE

save,WDL,DB,,all

! 温度荷载施加函数

h=arg1

To=arg2

Ti=arg3

dTh= (Ti-To)/h !温度梯度

Tk= (To+Ti)/2

TL= To-(1-0.5*h)*dTh

bfe,all,temp,1,Tk,Tk,TL !施加温差