Dear Prof. Taylor,
1) The equation I used is ρc∂T =−divq−Tref βTε ̇+Q, which is same in mannual86:section7.2.5. I add command "FOURier REFErence 1'' setting reference temperature to be 1. But the issue is still there. Does it work for version 8.5?
So far, I am not sure the problem is with my input or there is no suitable element I could use in FEAP8.5
2) Could FEAP 8.5 run a fully coupled thermo-mechanical 2D problem with governing equation:
ρc∂T =−divq−Tref βTε ̇+Q
ρa =div(sigma)+ρb
Best,
Eleven