{"id":312511,"date":"2023-10-20T12:13:14","date_gmt":"2023-10-20T12:13:14","guid":{"rendered":"\/forum\/forums\/topic\/ansys-piezoelectric-simulation-help\/"},"modified":"2023-10-20T12:13:14","modified_gmt":"2023-10-20T12:13:14","slug":"ansys-piezoelectric-simulation-help","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/ansys-piezoelectric-simulation-help\/","title":{"rendered":"ansys piezoelectric simulation help"},"content":{"rendered":"

Hi everyone! <\/p>\n

I am new to ansys and have been really struggling to get my code working…. I want to simulate the voltge output of a dielectic material under a uniform pressure acting on the top lectrode. Can someone please please please help me? Thanks !<\/p>\n

! Geometry<\/div>\n
l = 5e-3 ! beam length, m<\/div>\n
w = 5e-3 ! electrode width, m<\/div>\n
t = 100e-6 ! thickness, m<\/div>\n
 <\/div>\n
!Loading<\/div>\n
PRES= 10000 !pressure loading, Pa<\/div>\n
 <\/div>\n
! Material properties<\/div>\n
Y = 140e9 ! Young modulus, Pa<\/div>\n
mu=0.358 ! Poisson ratio (nearly incompressible rubber)<\/div>\n
eps=8.17 ! electrical permittivity, relative<\/div>\n
eps0=8.854e-12 ! free-space permittivity, F\/m<\/div>\n
 <\/div>\n
\/VUP,1,z<\/div>\n
\/VIEW,1,1,1,1<\/div>\n
 <\/div>\n
 <\/div>\n
\/PREP7<\/div>\n
ET, 1, 5, 3<\/div>\n
 <\/div>\n
! == Material parameters<\/div>\n
! — Elastic coefficients, MPa<\/div>\n
c11 = 209.71e9<\/div>\n
c12 = 121.14e9<\/div>\n
c13 = 105.36e9<\/div>\n
c33 = 210.9e9<\/div>\n
c44 = 42.5e9<\/div>\n
c66 = 44.28e9<\/div>\n
 <\/div>\n
tb,ANEL,1<\/div>\n
tbdata, 1, c11, c12, c13<\/div>\n
tbdata, 7, c11, c13<\/div>\n
tbdata, 12, c33<\/div>\n
tbdata, 16, c66<\/div>\n
tbdata, 19, c44<\/div>\n
tbdata, 21, c44<\/div>\n
 <\/div>\n
! — Piezoelectric coefficients, pC\/um2<\/div>\n
e13 = -0.61<\/div>\n
e33 = 1.14<\/div>\n
e15 = -0.59<\/div>\n
 <\/div>\n
tb,PIEZ,1<\/div>\n
tbdata, 3, e13<\/div>\n
tbdata, 6, e13<\/div>\n
tbdata, 9, e33<\/div>\n
tbdata,14, e15<\/div>\n
tbdata,16, e15<\/div>\n
 <\/div>\n
! — Dielectric constants<\/div>\n
emunit,EPZRO,8.854e-6 ! pF\/um<\/div>\n
mp,PERx,1, 7.57<\/div>\n
mp,PERy,1, 7.57<\/div>\n
mp,PERz,1, 8.31<\/div>\n
 <\/div>\n
! — model and mesh block<\/div>\n
block,-l\/2,l\/2,-w\/2,w\/2,0,t<\/div>\n
 <\/div>\n
! Set the merging tolerance (adjust as needed)<\/div>\n
KETOL = 1.0e-5<\/div>\n
vmesh,1<\/div>\n
 <\/div>\n
 <\/div>\n
! Structural BC – fix the displacement of bottom face<\/div>\n
nsel,s,loc,x,-l\/2<\/div>\n
d,all,ux,0<\/div>\n
nsel,r,loc,y,-w\/2<\/div>\n
d,all,uy,0<\/div>\n
nsel,r,loc,z,0<\/div>\n
d,all,uz,0<\/div>\n
nsel,all<\/div>\n
 <\/div>\n
! Electrical BC<\/div>\n
nsel,s,loc,z,0<\/div>\n
nsel,r,loc,x,-l\/2,l\/2<\/div>\n
nsel,r,loc,y,-w\/2,w\/2<\/div>\n
cp,1,volt,all<\/div>\n
ng=ndnext(0) ! ground node<\/div>\n
 <\/div>\n
nsel,all<\/div>\n
nsel,s,loc,z,t<\/div>\n
cp,2,volt,all<\/div>\n
nl=ndnext(0) ! load node<\/div>\n
nsel,all<\/div>\n
 <\/div>\n
\/SOLU<\/div>\n
antype,static<\/div>\n
cnvtol,f,1,1.e-6<\/div>\n
cnvtol,amps,1,1.e-6<\/div>\n
cnvtol,volts,1,1.e-6<\/div>\n
d,ng,volt,0 !put ground node to zero volts<\/div>\n
nsel,s,loc,z,t<\/div>\n
nsel,r,loc,x,-l\/2,l\/2<\/div>\n
nsel,r,loc,y,-w\/2,w\/2<\/div>\n
sf,all,PRES,0  !apply pressure<\/div>\n
 <\/div>\n
nsel, all<\/div>\n
solve<\/div>\n
fini<\/div>\n
 <\/div>\n
\/POST1<\/div>\n
nsel,all<\/div>\n
pldisp,1 ! display deformed\/undeformed shape<\/div>\n
andscl ! animate deformed\/undeformed shape<\/div>\n
nsel,s,loc,x,l\/2<\/div>\n
nd=ndnext(0) ! pick node for display<\/div>\n
nsel,r,node,,nd<\/div>\n
prnsol,epel ! print strain<\/div>\n
prnsol,volt ! print voltage<\/div>\n
nsel,all<\/div>\n
fini<\/div>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n","protected":false},"template":"","class_list":["post-312511","topic","type-topic","status-closed","hentry"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t