压电双晶体驱动悬臂梁变形的abaqus数值模拟

问题描述：聚乙烯悬臂梁，左端固定，右端自由，悬臂梁的上、下两个表面固定两个PZT-4压电体，在这两个压电体上输入电载荷驱动悬臂梁变形。

文献Haojiang Ding, Jian  Liang：  The  fundamental  solutions  for  transversely  isotropic  piezoelectricity  and boundary element method给出了材料参数。极化方向为3方向，具体材料参数如下。

这些压电材料的材料参数该如何输入呢？我们知道压电方程为

用Abaqus中的参数表示为

通过对比可以得到它们之间的关系，比如

其他就不在列举。

最终计算结果如下