Hi, I am trying to create a cantilever piezoelectric bimorph. I have successfully created the structure and found the master nodes of the electrodes. I am trying to connect the bimorph with a parallel circuit. My load is just a single resistor. One end has to connect to the middle electrode and the other end has to connect with both the top and bottom electrodes. I'm not really sure how to connect the resistor. Should I couple the voltage for the top and bottom electrodes?

My code is here:
finish
/clear
/COM,Bimorph Model

/FILNAME, Bimorph,0
/TITLE, Bimorph Power Analysis


*SET,Lm,17.3E-3 !Length of mass
*SET,Hm,7.7E-3 !Height of mass
*SET,Wm,3.6E-3 !Width of mass
*SET,Lb,10.7E-3 !Length of beam
*SET,Wb,3.2E-3   !Width of beam
*SET,Tp,0.278E-3 !Pizeoelectric Thickness
*SET,Tsh,0.102E-3 !Shim Thickness
*SET,Rload,151E3  !Resistance
/UNITS,MKS ! MKS system (m, kg, s, °C#

/PREP7
local,11
local,12,,,,,,180,,

ET,1,SOLID45  !Set Element Types
ET,2,SOLID5,3
ET,3,CIRCU94,0

CSYS,11
BLOCK,0,Lb+Lm,0,Wb,0,Tsh,  !Draw Shim
BLOCK,0,Lb,0,Wb,Tsh,Tsh+Tp, !Drawing Top Pizeoelectric Layer
BLOCK,0,Lb,0,Wb,0,-Tp,  !Drawing Bottom Pizeoelectric Layer
BLOCK,Lb,Lb+Lm,#(Wb/2#-(Wm/2##,##Wb/2#+#Wm/2##,Tsh,Tsh+Hm !Drawing Mass


EMUNIT,    MKS,  ! Define Piezoelectric #PZT PSI-5H4E) Material Model
MPTEMP, 1, 25, , , , , ! PZT Density
MPDATA,DENS,1,,7500
TB,PIEZ,1,,,0   ! Piezoelectric matrix
TBDATA,    3,    -12
TBDATA,    6,    -12
TBDATA,    9,    18.2
TBDATA,    14,    21.9
TBDATA,    16,    21.9

MPTEMP, 1, 25, , , , , ! Clamped Dielectric Constant
MPDATA,PERX,1,,1550  
MPDATA,PERY,1,,1550  
MPDATA,PERZ,1,,1390  

TB,    ANEL,    1,1,0 ! Anisotropic elastic matrix
TBTEMP,25
TBDATA,    1,    16.9E+10    ,    11.8E+10    ,10.9E+10
TBDATA,    7,    16.9E+10    ,    10.9E+10        
TBDATA,    12,    12.3E+10                
TBDATA,    16,    2.5E+10                
TBDATA,    19,    2.7E+10                
TBDATA,    21,    2.7E+10                

!*
  
MPTEMP, 1, 25, , , , , ! Define Shim #Brass) Material Model
MPDATA,DENS,2,,8530 ! Density of Shim
MPDATA,EX,2,,1.11E11  ! Elastic Modulus of Shim
MPDATA,PRXY,2,,0.375 ! Poissons Ratio of Shim
  
  
MPTEMP, 1, 25, , , , , ! Define Mass #Tungsten Alloy) Material Model
MPDATA,DENS,3,,17000 ! Density of Mass
MPDATA,EX,3,,4.11E11 ! Elastic Modulus of Mass
MPDATA,PRXY,3,,0.28  ! Poissons Ratio of Mass

VSEL,All
VGLUE,ALL    ! Combine all Volumes

VSEL,S,LOC,Z,0,Tsh, , ! Defining Element Type & Material Properties of layers
VSEL,R,LOC,X,0,Lb+Lm, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT,      2, ,   1,      11  !Shim Layer
VSEL,S,LOC,Z,Tsh,Tsh+Tp, ,
VSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT,    1, ,   2,       11  !Top Pizeoelectric Layer
VSEL,S,LOC,Z,0,-Tp, ,
VSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT,     1, ,   2,      11 !Bottom Pizeoelectric Layer
VSEL,S,LOC,Z,Tsh,Tsh+Hm, ,
NSEL,R,LOC,X,Lb,Lb+Lm, ,
VATT,      3, ,  1,       11 ! Mass
 
ESIZE,4E-4,0, ! Mesh Size

VSEL,ALL
VSWEEP,All !Sweep Volumes
 
NSEL,S,LOC,Z,-Tp,-Tp, ,
NSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
CP,1,VOLT,ALL    ! Couple nodes of bottom
*GET,ELEC_BOT,node,0,num,min

NSEL,S,LOC,Z,Tsh+Tp,Tsh+Tp, ,
NSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
CP,2,VOLT,ALL  ! Couple nodes of top
*GET,ELEC_TOP,node,0,num,min

NSEL,S,LOC,Z,0,0
NSEL,R,LOC,X,0,Lb, ,  ! Creating components
NSEL,R,LOC,Y,0,Wb, ,
ESLN,S,1
CM,cm_1,NODE
NSEL,S,LOC,Z,Tsh,Tsh
NSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
ESLN,S,1
CM,cm_2,NODE
CMSEL,S,CM_1
CMSEL,A,CM_2
CP,3,VOLT,ALL  ! Couple nodes of middle
*GET,ELEC_MID,node,0,num,min

ALLSEL,ALL  
/REPLOT
 
R,1,Rload, , , , ,   ! Define Resistor
RMORE, , , , , ,
RMORE, , , ,
!*  
TYPE,   3    
MAT,      0
REAL,       1   
ESYS,       11   
SECNUM,
TSHAP,LINE  
!*  
EN,99999,ELEC_TOP,ELEC_MID,, ! Creating Resistor
!*
ALLSEL,ALL  
/REPLOT

ASEL,S,LOC,Z,0,Tsh,,  ! Fixing End of Shim
ASEL,R,LOC,X,0,0,,  
DA,ALL,UX,0,,  
ASEL,S,LOC,Z,0,Tsh,,  
ASEL,R,LOC,X,0,0,,  
DA,ALL,UY,0,,
ASEL,S,LOC,Z,0,Tsh,,  
ASEL,R,LOC,X,0,0,,
DA,ALL,UZ,0,,

ALLSEL,ALL
/REPLOT
FINISH  

/SOL
ANTYPE,HARMIC ,NEW         ! Harmonic Analysis   
HROPT,FULL
EQSLV,SPARSE ! Sparse Solver
KBC,1
D,ELEC_MID,VOLT,0,  ! Setting Voltage on Middle to zero
ACEL,0,0,2.25, ! Acceleration
HARFRQ,104,106, ! Frequency Range
NSUBST,3,  !Substeps
   
!*  
ALPHAD,0,   
BETAD,0,
DMPRAT,0.02, ! Constant Damping Ratio
!*
ALLSEL,ALL  
/REPLOT  
/STATUS,SOLU
SOLVE

          
/POST26 ! Plot Power Graph
FILE,'Bimorph','rst','.'   
/UI,COLL,1  
NUMVAR,200  
SOLU,191,NCMIT  
STORE,MERGE
PLCPLX,0
PRCPLX,1
FILLDATA,191,,,,1,1
REALVAR,191,191
!*  
!*  
!*  
FORCE,TOTAL
SHELL,TOP   
LAYERP26,0  
ESOL,2,99999, ,NMIS,1,PWR
!*  
PLVAR,PWR, , , , , , , , , ,