thermal stress

Hi all,

I have a  20 node hexahedral "Thermomechanical" element in our inhouse software which i am trying to validate against Abaqus .

The problem is, the displacements, nodal forces match exactly ( Just compared the max min values in the legends! and checked a few random nodes).

The Abaqus input file is attached:

The problem is as follows:

One single element is chosen with a linear elastic material model (Some parameters are chosen in an unphysical manner but here the idea is

to just validate the implementation):