I am modeling a piezo electric slab with coupled field elements. I have activated the following dof's: ux,uy,uz,and volt. I believe that one of my constitutive equation then should be the following:
electric flux =(stress*compliance)+(piezo matrix*E field)
in the absence of electric flux...