clear all, clc

%Parameters of PZT
E=1/1.81e-11 ;% elastic modulus of the piezoelectric stack(Yangův modul)
d33= 6.5e-10;% piezoelectric constant
e33_e0=1750;% relativní permitivita 
Q=30;% quality factor

% Catalog PI product P-016.80
% http://www.piezo.ws/pdf/Piezo.pdf page18
C=1.3e-6;% capacitance of the piezoelectric stack
d=16e-3;% diameter of piezoelectric stack 
ome=69115;% natural frequency of the piezoelectric stack [rad/s]
k=56e+6;% stiffness of the piezoelectric stack

% Other parameters
h= 1e-4;% thickness of each piezoelectric layers
m=180.4e-3;% Stack mass
m_eff=(1/3)*m ;% effective mass of the piezoelectric stack 
A= (pi*d^2)/4;% cross-sectional area of each piezoelectric layer
n=floor(C*h/(8.854188e-12*e33_e0*A));% number of layers
%Damping
b_p=1/(2*Q); % relative dumping (poměrný útlum)
b=2*b_p*m*ome;% Damping of stack

 
 