想问一下压电换能器中的一个小问题模拟，小女先谢谢各位大侠了

zxjyjr

各位大侠，最近遇上一个压电换能器的问题，小女对ansys不是很懂，想请教一下
问题很简单，想做一个立方体压电快，再在上下各加一块导体片，然后用导线连起来上下的导体片，想看看挤压导体片的时候导线中的电流

我的想法是这样的，把上下表面的电压自由度分别couple，模拟上下导体片
再将上下表面用resistor连起来，然后下表面固定，上表面加位移载荷求解

但是它提示单元中有负电荷，就不知道怎么回事了
以下是建模，不完全，还需要手动把电阻加上，然后就可以求解了，但是求解会有问题，想请教一下各位大侠

/PREP7
/GRAPHICS,FULL   
  
!############################## parameter #######################
*AFUN,DEG
  
VOL=0
VOH=20
  
!###########################  Material property #################
  
!****************** 3---PZT **************************
et,3,solid5  ! UX,UY,UZ,TEMP,VOLT,MAG
!mp,nuxy,3,0.27
mp,dens,3,7.64e3  ! kg/m3
mp,perx,3,8.93e-9  ! F/m
mp,pery,3,8.93e-9  ! F/m
mp,perz,3,6.92e-9  ! F/m
mp,kxx,3,1
mp,murx,3,1
  
tb,anel,3   ! Anisotropic elastic matrix  72A??? 99.3.4
tbdata,1,11.565e10,6.489e10,6.229e10,0,0,0
tbdata,7,11.565e10,6.229e10,0,0,0
tbdata,12,9.298e10,0,0,0
tbdata,16,1.786e10,0,0
tbdata,19,1.786e10,0
tbdata,21,1.786e10
  
tb,piez,3   ! Piezoelectric matrix  72A??? 99.3.4
tbdata,1,0,0,-12.3082
tbdata,4,0,0,-12.3082
tbdata,7,0,0,20.7574
tbdata,10,0,0,0
tbdata,13,0,17.0357,0
tbdata,16,17.0357,0,0
  
!------------------------ model---------------------------
NUMSTR,KP,1
NUMSTR,LINE,1
NUMSTR,AREA,1

block,0,10,0,5,0,2
  
type,3
mat,3

esize,1
  
vmesh,all
  
!------------boundary condition-----------------------
  
nsel,s,loc,z,0  
d,all,uz
  
nsel,s,loc,z,2
d,all,uz,-0.3

!couple
nsel,s,loc,z,0
cp,1,volt,all

nsel,s,loc,z,2
cp,2,volt,all