October 20, 2023 at 12:57 pm

Caterina Azais

Subscriber

thank you so much, that seems to be making much more sense now... I have changed the code as follows, however when i print the PRNSOL commands it says that the S, D, .. data are not available and hence its ignored. What is the problem? Thank you !

! Geometry

l = 5e-3 ! beam length, m

w = 5e-3 ! electrode width, m

t = 100e-6 ! thickness, m

! Material properties

Y = 140e9 ! Young modulus, Pa

mu=0.358 ! Poisson ratio (nearly incompressible rubber)

eps=8.17 ! electrical permittivity, relative

eps0=8.854e-12 ! free-space permittivity, F/m

/VUP,1,z

/VIEW,1,1,1,1

/PREP7

ET, 1, 5, 3

! == Material parameters

! -- Elastic coefficients, MPa

c11 = 209.71e9

c12 = 121.14e9

c13 = 105.36e9

c33 = 210.9e9

c44 = 42.5e9

c66 = 44.28e9

tb,ANEL,1

tbdata, 1, c11, c12, c13

tbdata, 7, c11, c13

tbdata, 12, c33

tbdata, 16, c66

tbdata, 19, c44

tbdata, 21, c44

! -- Piezoelectric coefficients, pC/um2

e13 = -0.61

e33 = 1.14

e15 = -0.59

tb,PIEZ,1

tbdata, 3, e13

tbdata, 6, e13

tbdata, 9, e33

tbdata,14, e15

tbdata,16, e15

! -- Dielectric constants

emunit,EPZRO,8.854e-6 ! pF/um

mp,PERx,1, 7.57

mp,PERy,1, 7.57

mp,PERz,1, 8.31

! -- model and mesh block

block,-l/2,l/2,-w/2,w/2,0,t

! Set the merging tolerance (adjust as needed)

KETOL = 1.0e-5

vmesh,1

! Structural BC

nsel,s,loc,z,0

d, all, volt, 0

d,all,ux,0

d,all,uy,0

d,all,uz,0

nsel,all

nsel,s,loc,z,t

SF,ALL,PRES,10000 ! APPLY PRESSURE LOAD

NSEL,ALL

FINISH

/SOLVE

OUTPR,,LAST

/OUT,SCRATCH

SOLVE

FINISH

/POST1

/OUT,

PRNSOL,S,COMP

PRNSOL,EPEL,COMP

PRNSOL,EF,COMP

PRNSOL,D,COMP