Changing the value of density during thermal analysis

[Feddozz]

Hi everyone,
I need to perform a thermal analysis. A slab of solid material subjected to one side heat flux. My energy balance equation takes in account conduction and heat generation.
The conduction term is the usual term.
The heat generation term is function of the density, but the density is a function of temperature.
I think I should perform a transient analysis and manage to change the value of the density during the analysis after each substep, at least.

Does someone know how to do it with ABAQUS?

 

[xerf]

Look at HETVAL user subroutine.