Deformation gradient

When we only have mechanical loading, I understand that abaqus calculates the total deformation gradient using F = I + partial_U/partial_X and U = N_i*U_i where U_i is the displacement at the nodes  and N is the shape function .

If I use a thermal expansion in addition to mechanical loading , Does abaqus give the total deformation gradient incorporating both mechanical loading and  thermal expansion ?
F = I + partial_Utot/partial_X ; Utot = U_mechanical + U_thermal

Hello,

I have a 2D mesh that has 4 node solid elements. I am aware that to determine the components of deformation gradient tensor F, I will require displacement functions, u(x,y) and v(x,y) in a current configuration. Taking a partial derivative with respect to x and y respectively yields components of def. gradient tensor F.