I suggess you use an artifical temperature field, that you can define easily in ABAQUS or ANSYS, then you can prescribe elastic or another properties temperature dependant. I have used this approach to modell bamboo, which is a FGM material.

Yu can also use a user-defined field variable and even update that using ufield/usdfld. Make the material properties depend on that variable (keyword depvar)

i dnt hve an idea about abacus , but in ansys i already did a macro which be able to change the element material properties in elements of a cylinder according to the power low distribution, the macro is given as

*create, mat_fgm_K

/nopr

R_coord = arg1

grad_common = ((R_coord-cylin_radin)/(cylin_radout-cylin_radin))**B

E_fgm = (E_c-E_m)*grad_common + E_m

nu_fgm = (nu_c-nu_m)*grad_common + nu_m

den_fgm = (dens_c-dens_m)*grad_common + dens_m

/gopr

*end !mat_fgm_K

to change the element material i used this macro

CSYS,1

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!changing element material !!!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

counter = 0

*do, j, R1,R2,diff

R_coord = j

counter = counter+1

*use, mat_fgm_K,R_coord

mp, ex, counter, E_fgm

mp, prxy, counter, nu_fgm

mp, dens, counter, den_fgm

nSEL,S,LOC,X,j

ESLN,S

emodif, All, mat, counter

esel, all

*enddo

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!end-changing element material !!!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

hope this help

regards

I think that, Bendine Kouider Macro will be very helpful and sufficient to calculate the modulus of elasticity, Poisson's ratio and density required to make material properties graduation. Moreover, it helps in the required modification to the concerned elements. 

Regards