blob: 9e5a7ccefd909150b410842d60c5f33c9e3610cf (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
clc();
clear;
//Given :
E1 = 8.55*10^10; //Modulus of elasticity in N/m^2
E2 = 21*10^10; // Modulus of elasticity in N/m^2
rho1 = 2650; // density of Quartz in kg/m^3
rho2 = 8800;// density of Nickel in kg/m^3
t = 2; // thickness of crystal in mm
l = 50; // rod length in mm
//Piezoelectric generator
printf("Piezoelectric generator \n\n");
for n = 1:3
// 1 mm = 1.0*10^-3 m
nu1 = (n/(2*t*10^-3))*sqrt(E1/rho1);// frequency in Hz
printf("For n = %d , Frequency = %.2f MHz\n",n,nu1*10^-6);
end
//Magnetostriction generator
printf("Magnetostriction generator\n\n");
for n1 = 1:3
// 1 mm = 1.0*10^-3 m
nu2 = (n1/(2*l*10^-3))*sqrt(E2/rho2);// frequency in Hz
printf("For n = %d , Frequency = %.1f kHz\n",n1,nu2*10^-3);
end
//Results differ from those in textbook, because in the formulae (n/(2*t))*sqrt(E/rho) and (n/(2*l))*sqrt(E/rho) , 2 is not multiplied with either t or l.
|
