!Verhalten einer Piezo-Scheibe

/Filename,Piezoblock_Sensor
/Title,Piezoblock_Sensor

/prep7

*afun,deg

/units,SI

/view,1,1,1,1

csys,0


et,1,solid226,1001

!Materialdaten

mp,dens,2,7800		!Dichte

tb,piez,2		!Material-Tabelle Piezoeigenschaften
tbdata,3,-7.15		!e13
tbdata,6,-7.15		!e23
tbdata,9,13.70		!e33
tbdata,14,11.90		!e52
tbdata,16,11.90		!e61


mp,perx,2,930		!eps11Sr
mp,pery,2,930
mp,perz,2,857		!eps33Sr

tb,anel,2		!Material-Tabelle anisotrop-elastische Eigenschaften: Steifigkeit c
tbdata,1,12.3e10,7.67e10,7.025e10
tbdata,7,12.3e10,7.025e10		!->habe ich nicht im Datenblatt!
tbdata,12,9.711e10
tbdata,16,2.315e10
tbdata,19,2.226e10
tbdata,21,2.226e10

!Geometrie

l=40/1000		!Laenge [m]
b=10/1000		!Breite [m]
h=0.5/1000		!Hoehe [m]


blc5,0,0,b,l,h

lsel,s,line,,6

lesize,all,,,14

allsel

lsel,s,line,,7

lesize,all,,,6

allsel


vsweep,all

/solu

antype,trans

trnopt,full

nsel,s,loc,z,0

d,all,uz,0
d,all,volt,0

allsel

tend=0.2

schritte=1000

dt=tend/schritte

*dim,kraftz,table,schritte,1

*do,s,1,schritte
	kraftz(s,0)=s*dt
	kraftz(s,1)=100*cos(s*dt)
*enddo




*do,s,1,schritte
	time,dt*s
	nsel,s,loc,z,h
	f,all,fz,-kraftz(dt*s)
	nsub,1
	solve
*enddo



/post26

nsol,2,320,u,z,uz
nsol,3,320,volt,z,Uz

plvar,2,3