UMAT piezoelectric UEL ABAQUS

Hello,

 I am using abaqus to do research on polarization procedure of piezoelectric materials. I want to develop a user subroutine that can allow me to input variable piezoeletric material properties and do coupled problem. I want to know if UMAT can solve this problem or I need to use other subroutine. If there is that subroutine, what is it.  Any suggestion is appreciated.

 Thanks in advance,

Linhui Zhang