Meshing help

[gregthomas]

Hi,

I've got a very thin plate. It must be piezoelectric, which limits the number of elements I can use to mesh it. I'm having a lot of trouble meshing it and getting an analysis to run. The code I've pasted below gives warnings, and then calculates deflections of 0 everywhere after I apply some forces to the plate. Can anyone help? Thanks,

Greg



/prep7


mp,DENS,1,7700 ! Density, kg/m**3

tb,DPER,1 ! Relative permittivity at constant strain
tbdata,1,729,635,729

tb,ANEL,1 ! Anisotropic elastic stiffness, N/m^2
tbdata,1,13.9E10,7.43E10,7.78E10 ! c11,c13,c12
tbdata,7,11.5E10,7.43E10 ! c33,c13
tbdata,12,13.9E10 ! c11
tbdata,16,2.56E10 ! c44
tbdata,19,2.56E10 ! c44
tbdata,21,3.06E10 ! c66

tb,PIEZ,1 ! Piezoelectric stress coefficients, C/m^2
tbdata,2,-5.2 ! e31
tbdata,5,15.1 ! e33
tbdata,8,-5.2 ! e31
tbdata,10,12.7 ! e15
tbdata,15,12.7 ! e15


mp,dens,2,2330
mp,ex,2,1.65E+011
mp,nuxy,2,0.22


!
! Define a membrane
!
! Transducer size, m
block,0,500e-6,0,400e-6,0,0.5e-6 ! Define volume


et,1,solid5,0
MSHKEY,1
LESIZE,7,,,16
LESIZE,6,,,16
mat,1 $ type,1 ! Set element attributes


vmesh,1


nsel,s,loc,x,0 ! Impose displacement constraints
d,all,uz,0
d,all,ux,0
d,all,uy,0
nsel,s,loc,x,500e-6 ! Impose displacement constraints
d,all,uz,0
d,all,ux,0
d,all,uy,0



allsel

f,331,fz,-1
f,57,fz,1


fini

/SOLU
betad,5e-2
solve
FINISH