Why am I getting the error message below about negative permittivity for this piezoelectric property?
*** ERROR ***
The calculated permittivity at constant strain matrix [epS] for
material 1 has negative diagonal terms.
MATERIAL INPUT ASSOCIATED WITH TB,DPER COMMAND FOR MATERIAL NUMBER 1
Dielectric (relative) permittivity at constant strain
computed from piezoelectric strain and stress matrices.
EP11 EP22 EP33 EP12 EP23 EP13
69.025 52.629 0.0000 0.0000 0.0000 0.0000
Looking at the manufacturer's data, the set of compliance coefficients is incomplete and is inconsistent:  s12 and s13 are missing and are assigned small values  s11 is measured at constant electric field (E), while s33 and s55 are measured at constant electric flux (D). SE and SD are used in different sets of piezoelectric constitutive equations. You cannot intermix them in one input. The relationship between them is the following: [SD]=[SE] [d]T[g], where [d] is the matrix of piezoelectric strain coefficients and [g] is another set of piezoelectric coefficients associated with another set of constitutive equations. There are four sets of piezoelectric matrices in total: [e], [d], [h] and [g]. The latter piezoelectric matrix can be derived from [d] and the inverse of [epT] as follows: [g]=[epT]1[d] where [epT] is the permittivity matrix at constant stress. These input inconsistencies are probably the source of nonphysical values of permittivity [epS]. Below is APDL to perform the calculation as it is done internally. It confirms the values. epzro=8.854e12 c11=0.66686e11 c13=0.18833e9 c12=0.11117e10 c33=0.11111e12 d31= 0.18e9 !0.15e9 d33=0.4e9 e31=c11*d31+c13*d33+c12*d31 e33=c13*d31+c33*d33+c13*d31 /com e31=%e31% /com e33=%e33% ep22t=2400*epzro ep22s=ep22t2*e31*d31e33*d33 ep22s=ep22s/epzro /com ep22s=%ep22s% 

