summaryrefslogtreecommitdiff
path: root/1319/CH1/EX1.41/1_41.sce
blob: 95aa107a583b861b00829e238f34523806b01519 (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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
// To determine the parameters of a toroid

clc;
clear;

//Arc length of different materials
lni=0.3; // nickel iron alloy
lss=0.2;//Silicon Steel
lcs=0.1;//Cast steel

A=1*(10^-3); // Area of cross section of Toroid

N=100; // Hundred Turns

phi=6*(10^-4); // Flux
B=0.6; // Flux densities
m0=4*%pi*(10^-7);

//Field intensities
Hni=10; // Ni alloy
Hss=77;// Si Steel
Hcs=270; // CAst steel

mmf=(Hni*lni)+(Hss*lss)+(Hcs*lcs); // Total mmf

B=phi/A;

I=mmf/N;

//P = Permeability, RP = Relative permeablity

//For Nickel Alloy
mni=B/Hni; // P
mrni=mni/m0; // RP

//For Si Steel
mss=B/Hss; //P
mrss=mss/m0; // RP

//For Cast Steel
mcs=B/Hcs; //P
mrcs=mcs/m0; // RP

deff('x=rel(y,z)','x=(y*z)/phi'); // Fucntion to find out Reluctance

//Reluctances
Rni=rel(Hni,lni);
Rss=rel(Hss,lss);
Rcs=rel(Hcs,lcs);

printf('i) The mmf required to establish a magnetic flux of 6 *(10^-4) Wb = %g AT\n',mmf)
printf('ii) The Current through the coil = %g A\n',I)
printf('iii) Relative permeability and reluctance of each ferro magnetic material\n')
printf('    Nickel Iron alloy is %g and %g ohms \n',mrni,Rni)
printf('    Silicon Steel is %g and %g ohms \n',mrss,Rss)
printf('    Cast Steel is %g and %g ohms \n',mrcs,Rcs)