1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
clc;
N=1500; // number of turns in coil
i=3; // current carried by coil
uo=4*%pi*10^-7; // free space permeability
l=0.04; // side of plunger
A=%pi*(l/2)^2; // cross section area of plunger
disp('case a');
disp('for gap length G=2 cm');
g=0.02; // air gap length in cm
x=0; // displacement of plunger
G1=g-x; // gap length
B=(uo*i*N)/G1; // air gap flux density
printf('Air gap flux density is %f Wb/m^2\n',B);
L1=(N^2*uo*A)/G1;
printf('Coil inductance is %f H',L1);
disp('for gap length G=1.5 cm');
g=0.02; // air gap length in cm
x=0.005; // displacement of plunger
G2=g-x; // gap length
B=(uo*i*N)/G2; // air gap flux density
printf('Air gap flux density is %f Wb/m^2\n',B);
L2=(N^2*uo*A)/G2;
printf('Coil inductance is %f H',L2);
disp('for gap length G=1 cm');
g=0.02; // air gap length in cm
x=0.01; // displacement of plunger
G3=g-x; // gap length
B=(uo*i*N)/G3; // air gap flux density
printf('Air gap flux density is %f Wb/m^2\n',B);
L3=(N^2*uo*A)/G3;
printf('Coil inductance is %f H',L3);
disp('for gap length G=0.5 cm');
g=0.02; // air gap length in cm
x=0.015; // displacement of plunger
G4=g-x; // gap length
B=(uo*i*N)/G4; // air gap flux density
printf('Air gap flux density is %f Wb/m^2\n',B);
L4=(N^2*uo*A)/G4;
printf('Coil inductance is %f H',L4);
disp('case b');
disp('for gap length G=2 cm');
W=(i^2*L1)/2;
printf('Energy stored is %f watt-sec\n',W);
disp('for gap length G=1.5 cm');
W=(i^2*L2)/2;
printf('Energy stored is %f watt-sec\n',W);
disp('for gap length G=1 cm');
W=(i^2*L3)/2;
printf('Energy stored is %f watt-sec\n',W);
disp('for gap length G=0.5 cm');
W=(i^2*L4)/2;
printf('Energy stored is %f watt-sec\n',W);
disp('case c');
disp('for gap length G=2 cm');
fe=round((i^2*g*L1)/(2*G1^2));
printf('Electromagnetic force is %f N\n',fe);
disp('for gap length G=1.5 cm');
fe=(i^2*g*L1)/(2*G2^2);
printf('Electromagnetic force is %f N\n',fe);
disp('for gap length G=1 cm');
fe=round((i^2*g*L1)/(2*G3^2));
printf('Electromagnetic force is %f N\n',fe);
disp('for gap length G=0.5 cm');
fe=round((i^2*g*L1)/(2*G4^2));
printf('Electromagnetic force is %f N\n',fe);
disp('case 4');
// for g=2 cm and g=0.5cm, displacement is given by
xi=0;
xf=0.015;
Wm=integrate('(i^2*g*L1)/(2*(g-x)^2)','x',xi,xf);
printf('Mechanical work done is %f watt-sec\n',Wm);
disp('case e');
We=integrate('(i^2*g*L1)/(g-x)^2','x',xi,xf);
printf('Electrical energy supplied by source is %f watt-sec\n',We);
|
