压电耦合分析

可以啊
以下是压电矩阵的输入
/com -- Material matrices for PZT4 (polar axis along Y-axis): ANSYS input
/com
/com [c11 c13 c12 0 0 0 ] [ 0 e31 0 ]
/com [c13 c33 c13 0 0 0 ] [ 0 e33 0 ]
/com [c12 c13 c11 0 0 0 ] [ 0 e31 0 ]
/com [ 0 0 0 c44 0 0 ] [e15 0 0 ]
/com [ 0 0 0 0 c44 0 ] [ 0 0 e15]
/com [ 0 0 0 0 0 c66] [ 0 0 0 ]
c11=13.2e10
c12=7.1e10
c13=7.3e10
c33=11.5e10
c44=2.6e10
c66=3e10
tb,anel,1 !Define structural table
tbdata,1,c11,c13,c12 !Input elastic stiffness matrix [c]
tbdata,7,c33,c13
tbdata,12,c11
tbdata,16,c44
tbdata,19,c44
tbdata,21,c66
e13=-4.1
e33=14.1
e15=10.5
tb,piez,1 !Define Piez. table
tbdata,2,e13 !Input Piezoelectric stress matrix [e]
tbdata,5,e33
tbdata,8,e13
tbdata,10,e15
tbdata,15,e15

ze怎么输入压电矩阵