summaryrefslogtreecommitdiff
path: root/1319/CH1/EX1.25/1_25.sce
blob: 6ff58b51a6c5c6f9035b0dd14de5699b6278aedf (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
25
26
27
28
29
30
31
32
//To find current in the 600 turn exciting coil

clc;
clear;

N=600;
mr=800;
m0=4*%pi*(10^-7);

phi=100*(10^-6);// Flux in air gap

l1=10*(10^-2);
l2=18*(10^-2);
lg=2*(10^-3); // Air gap length
Ac=(6.25)*(10^-4);// Central limb area
As=3*(10^-4);// Side limb area

Ra=lg/(m0*Ac);
Ri=l1/(mr*m0*Ac);

R=l2/(m0*mr*As);

Rt=Ra+Ri; // Total reluctacne of the central limb

AT1=Rt*phi; // MMF or Ampere turns for the central limb
AT2=R*phi/2; // MMF, Two identical limbs hence flux becomes half and only one limb is considered

AT=AT1+AT2;//Total MMF

I=AT/N; // Current in the 600 turns

printf('The current flowing in the 600 turns exciting coil = %g A\n',I)