石英悬臂梁的静态仿真的问题

boyeexie

我做石英悬臂梁的静态仿真，如图1所示，梁的右端面固定约束，左端面施加Z方向的力

                               
登录/注册后可看大图

                               
登录/注册后可看大图

通过仿真得到的电势云图如图2所示，在梁的+-z面上存在的电势差，按照石英的压电方程，其在Z面上是不可能产生电荷的，为什么在+-z表面上有电势呢？
按照压电方程，由于剪应力XZ的作用，该梁应该在+-Y方向有电势差，为什么会得到这样的仿真结果呢？
请高手明示，是我的理解错了，还是我的仿真错了？
附仿真命令流在楼下

                               
登录/注册后可看大图

                               
登录/注册后可看大图

boyeexie

命令流


/PREP7
/TITLE, Quartz electric Field analysis

/COM,   Element1
ET,1,SOLID98
KEYOPT,1,1,3
KEYOPT,1,3,0
KEYOPT,1,5,0

/COM,  -- MATERIAL MATRICES : ANSYS INPUT
/COM,   - COMPLIANCE COEFFICIENTS, M2/N
S11=12.77E-12         
S12=-1.79E-12
S13=-1.22E-12
S14=  -4.5E-12            !有的资料为4.5E-12
S33=  9.6E-12
S44=20.04E-12
TB,ANEL,1,,,1                     ! ANISOTROPIC ELASTIC COMPLIANCE MATRIX
TBDATA,1,  s11,  s12,  s13,    0,  s14,     0  
TBDATA,7,        s11,  s13,    0, -s14,     0
TBDATA,12,             s33,    0,    0,     0
TBDATA,16,           2*(s11-s12),    0, 2*s14
TBDATA,19,                         s44,     0
TBDATA,21,                                s44

/COM,  - PIEZOELECTRIC STRAIN COEFFICIENTS, C/N
d11=2.31E-12
d14=0.73E-12             !有的资料为d11=-2.3E-12, d14=-0.67E-12
TB,PIEZ,1,,,1                     ! PIEZOELECTRIC STRAIN MATRIX
TBDATA,1,  d11,      0,   0
TBDATA,4, -d11,      0,   0
TBDATA,7,    0,      0,   0
TBDATA,10,   0, -2*d11,   0
TBDATA,13, d14,      0,   0
TBDATA,16,   0,    d14,   0

EMUNIT,MKS !Free-space permittivity is set to 8.85E-12 F/m.
!/COM,  - RELATIVE PERMITTIVITY AT CONSTANT STRESS
!EP11=4.52
!EP33=4.68         
!TB,DPER,1,,,1,       ! Permittivity matrix at constant stress
!TBDATA,1, EP11, EP11, EP33

/COM,  - RELATIVE PERMITTIVITY AT constant strain
EP11=4.5
EP33=4.6
MP,PERX,1,EP11                    ! Relative PERMITTIVITY (X DIRECTION)
MP,PERY,1,EP11                    ! Relative PERMITTIVITY (Y DIRECTION)
MP,PERZ,1,EP33                    ! Relative PERMITTIVITY (Z DIRECTION)

/COM,  - DENSITY, Kg/m3
RHO=2651
MP,DENS,1,RHO                     ! DENSITY

MAT,1
TYPE,1

/PNUM,KP,1     
KPLOT  
K,1,-0.0005,0.0001,0,  
K,2,0.0005,0.0001,0,   
K,3,0.0005,-0.0001,0,  
K,4,-0.0005,-0.0001,0,
K,5,-0.0005,0.0001,0.0002,  
K,6,0.0005,0.0001,0.0002,   
K,7,0.0005,-0.0001,0.0002,  
K,8,-0.0005,-0.0001,0.0002,
V,1,2,3,4,5,6,7,8
ALLSEL,ALL
SMRT,1      !Specifies meshing parameters for automatic (smart) element sizing
MSHAPE,1,3D !0 is quadrilateral elements when2D or hexahedral elements when3D;
            !1 is triangle elements 2D or tetrahedral elements when3D.
MSHKEY,0    !free mesh is 0; mapped mesh is 1; free and mapped mesh is 2
VMESH,ALL
EREF,ALL    !mesh more
FINISH

/SOL
ANTYPE,0
DA,3,ALL,0  
ASEL,S, , ,5
NSLA,S  
F,ALL,FZ,1e-006
ALLSEL,ALL
SOLVE

boyeexie

有做过石英梁压电分析的吗？帮忙看看！谢谢！