石英悬臂梁的静态仿真的问题
我做石英悬臂梁的静态仿真，如图1所示，梁的右端面固定约束，左端面施加Z方向的力http://forum.simwe.com/templates/simwe/images/attachimg.gif http://forum.simwe.com/attachments/day_090415/0904152316315832ac78bb8b93.jpg.thumb.jpg
通过仿真得到的电势云图如图2所示，在梁的+-z面上存在的电势差，按照石英的压电方程，其在Z面上是不可能产生电荷的，为什么在+-z表面上有电势呢？
按照压电方程，由于剪应力XZ的作用，该梁应该在+-Y方向有电势差，为什么会得到这样的仿真结果呢？
请高手明示，是我的理解错了，还是我的仿真错了？
附仿真命令流在楼下
http://forum.simwe.com/templates/simwe/images/attachimg.gif http://forum.simwe.com/attachments/day_090415/0904152316a11106efad95a473.jpg.thumb.jpg
命令流
命令流/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 有做过石英梁压电分析的吗？帮忙看看！谢谢！
页:
[1]