diff options
Diffstat (limited to '3835/CH3')
| -rw-r--r-- | 3835/CH3/EX3.1/Ex3_1.sce | 14 | ||||
| -rw-r--r-- | 3835/CH3/EX3.10/Ex3_10.sce | 17 | ||||
| -rw-r--r-- | 3835/CH3/EX3.11/Ex3_11.sce | 14 | ||||
| -rw-r--r-- | 3835/CH3/EX3.12/Ex3_12.sce | 35 | ||||
| -rw-r--r-- | 3835/CH3/EX3.13/Ex3_13.sce | 14 | ||||
| -rw-r--r-- | 3835/CH3/EX3.2/Ex3_2.sce | 35 | ||||
| -rw-r--r-- | 3835/CH3/EX3.3/Ex3_3.sce | 25 | ||||
| -rw-r--r-- | 3835/CH3/EX3.7/Ex3_7.sce | 35 | ||||
| -rw-r--r-- | 3835/CH3/EX3.8/Ex3_8.sce | 34 | ||||
| -rw-r--r-- | 3835/CH3/EX3.9/Ex3_9.sce | 29 |
10 files changed, 252 insertions, 0 deletions
diff --git a/3835/CH3/EX3.1/Ex3_1.sce b/3835/CH3/EX3.1/Ex3_1.sce new file mode 100644 index 000000000..a75cd142b --- /dev/null +++ b/3835/CH3/EX3.1/Ex3_1.sce @@ -0,0 +1,14 @@ +clear +// +//given +pi=3.14 +l=%pi*0.2 //l=mean length of the ring=%pi*mean diameter of the ring +A=400*10**-6 //A=cross sectional area of ring +u1=1000 //u1=relative permeability of steel +u2=4*%pi*10**-7 //relative permeability of air +R=l/(A*u1*u2) //reluctance of steel ring +printf("\n The reluctance of steel ring is= %0.0f AT/Wb",R) +//case b +flux=500*10**-6 +f=flux*R +printf("\n The magnetomotive force is= %0.0f AT",f) diff --git a/3835/CH3/EX3.10/Ex3_10.sce b/3835/CH3/EX3.10/Ex3_10.sce new file mode 100644 index 000000000..7af48aec2 --- /dev/null +++ b/3835/CH3/EX3.10/Ex3_10.sce @@ -0,0 +1,17 @@ +clear +// +//given +n=2000 //number of turns +flux=0.05*10**-3 //Wb +i=10 //A +lx=(n*flux)/i //self inductance in X +printf("\n lx= %0.5f H",lx) +//since coils are identical self inductance in Y=self inductance in x +fluxlinkingX=0.75*0.05*10**-3 //Wb flux linking due to current in coil X +fluxlinkingY=2000*0.05*0.75*10**-3 //Wb flux linkages in coil Y +m=fluxlinkingY/5 //mutual inductance +printf("\n m= %0.5f H",m) +//The rate of change in current di/dt=2000A/sec --> di/dt=(10-(-10))/0.01 +rate=2000 +ey=m*rate +printf("\n The induced emf in coil Y= %0.0f V",ey) diff --git a/3835/CH3/EX3.11/Ex3_11.sce b/3835/CH3/EX3.11/Ex3_11.sce new file mode 100644 index 000000000..6a4737c32 --- /dev/null +++ b/3835/CH3/EX3.11/Ex3_11.sce @@ -0,0 +1,14 @@ +clear +// +//given +//when currents are in same direction the total induction is: +//lt=l1+l2+2m +//when currents are in opposite direction the total emf is: +//lt=l1+l2-2m +//According to this problem +//l1+l2+2m=1.2 +//l1+l2-2m=0.2 +//Solving the above equations we get l1=0.4H M=0.25H +//on substituting we get l2=0.3H +//k=m/squareroot(l1*l2) +printf("\n k=0.72168") diff --git a/3835/CH3/EX3.12/Ex3_12.sce b/3835/CH3/EX3.12/Ex3_12.sce new file mode 100644 index 000000000..917043a7e --- /dev/null +++ b/3835/CH3/EX3.12/Ex3_12.sce @@ -0,0 +1,35 @@ +clear +// +//given +//case a +B=1 //Wb/m**2 +u1=4*3.14*10**-7 +A=10**-4 //cm**2 +per=800 //permeability +n=250 //number of turns +flux=B*A +printf("\n flux %0.5f Wb",flux) +r=781250 //AT/Wb calculated using formula for reluctance +mmf=flux*r //AT +i=mmf/n //exciting current required in A +printf("\n i %0.5f A",i) +l=(n*flux)/i //self inductance of the coil +printf("\n l= %0.5f H",l) +w=(l*i*i)/2 //energy stored +printf("\n w= %0.5f J",w) +//case b +airgap=1*10**-3 //air gap is assumed +rair=airgap/(u1*A) //reluctance of air gap in AT/Wb +mmfa=flux*rair //mmf of air in AT +printf("\n mmfa") +//rcore=((2.5*3.14)-0.1)/(32*3.14*10**-6) //reluctance of core +//mmfc=flux*rcore +mmfc=780 //AT +F=mmfc+mmfa +I=F/n //A +printf("\n exciting current= %0.2f A",I) +n=250 //number of turns +L=(n*flux)/I //self inductanc eof coil with air gap +printf("\n l= %0.5f H",L) +e=(L*I*I)/2 //energy stored in coil +printf("\n e= %0.5f J",e) diff --git a/3835/CH3/EX3.13/Ex3_13.sce b/3835/CH3/EX3.13/Ex3_13.sce new file mode 100644 index 000000000..14d732c4a --- /dev/null +++ b/3835/CH3/EX3.13/Ex3_13.sce @@ -0,0 +1,14 @@ +clear +// +//given +A=10**-1 //area +flux=0.1 //Wb +//case a +B=flux/A //flux density Wb/m**2 +u1=4*3.14*10**-7 +F=(B*B*A)/(2*u1) //force in N +printf("\n force= %0.5f N",F) +//case b +l=10**-2 //length of the air gap +w=(B*B*A*l*2)/(2*u1) //energy stored in two airgaps, 2=air gaps +printf("\n W= %0.5f J",w) diff --git a/3835/CH3/EX3.2/Ex3_2.sce b/3835/CH3/EX3.2/Ex3_2.sce new file mode 100644 index 000000000..5c1ee3de9 --- /dev/null +++ b/3835/CH3/EX3.2/Ex3_2.sce @@ -0,0 +1,35 @@ +clear +// +//given +l=0.5 +A=4*10**-4 +N=250 +I=1.5 +flux=0.25*10**-3 +fluxdensity=flux/A +f=N*I //magnetomotive force +H=(N*I)/l //magnetic field strength +pi=3.14 +u1=4*%pi*10**-7 +u2=fluxdensity/(u1*H) +printf("\n The flux density is= %0.3f Wb/m**2",fluxdensity) +printf("\n The magnetomotive force is= %0.0f AT",f) +printf("\n The magnetic field strength is= %0.0f AT/m",H) +printf("\n The relative permeability is= %0.0f ",u2) +//case b +//given +I=5 +flux=0.6*10**-3 +A=4*10**-4 +N=250 +l=0.5 +fluxdensity=flux/A +printf("\n The flux density is= %0.1f Wb/m**2",fluxdensity) +f=N*I //magnetomotive force +printf("\n The magnetomotive force is= %0.0f AT",f) +H=(N*I)/l //magnetic field stength +printf("\n Magnetic field strength= %0.0f AT/m",H) +pi=3.14 +u1=4*%pi*10**-7 +u2=fluxdensity/(u1*H) +printf("\n The relative permeability is= %0.1f ",u2) diff --git a/3835/CH3/EX3.3/Ex3_3.sce b/3835/CH3/EX3.3/Ex3_3.sce new file mode 100644 index 000000000..0c20977ac --- /dev/null +++ b/3835/CH3/EX3.3/Ex3_3.sce @@ -0,0 +1,25 @@ +clear +// +//given +pi=3.14 +N=250 +I=1.5 +ls1=0.627 //mean length of steel string +la=0.0001 //length of air gap +A=4.91*10**-4 //cross sectional area of magnetic circuit +f=N*I //magnetomotive force +printf("\n Magnetomotive force= %0.0f AT",f) +fa=1050 //fa=mmf of air gap=1050AT +fs=450 //fs=mmf of steel ring=450 +//case b +u1=4*%pi*10**-7 +ra=la/(u1*A) //reluctance of air gap +printf("\n The reluctance of air gap is= %0.3f AT/Wb",ra) +flux=fa/ra +printf("\n The flux is= %0.2f Wb",flux) +//case c +fluxdensity=flux/A +printf("\n The flux density is= %0.5f Wb/m**2",fluxdensity) +//case d +rs=fs/flux //reluctance of steel string +printf("\n The reluctance of steel string is= %0.6f AT/Wb",rs) diff --git a/3835/CH3/EX3.7/Ex3_7.sce b/3835/CH3/EX3.7/Ex3_7.sce new file mode 100644 index 000000000..5d2fa0d55 --- /dev/null +++ b/3835/CH3/EX3.7/Ex3_7.sce @@ -0,0 +1,35 @@ +clear +// +//given +fluxa=1.4*10**-3 +area=0.002 +B=fluxa/area //flux density in air gap +printf("\n B= %0.3f Wb/m**2",B) +//u1=4*3.14*10**-7 +//ha=B/u1 in AT/m //magnetic field in air gap +ha=55.7 +la=2 //length of air gap in m +mmf=ha*la //mmf of air gap +printf("\n mmf= %0.3f AT",mmf) +//since the flux density of central limb is 0.7 the corresponding field srength is h1=250AT/m +h1=250 +mmfl=112.45 //mmf for magnetic central limb-->mmf=250*(450-0.2)*10**-3 +totmmf=mmf+mmfl +printf("\n totmmf= %0.5f AT",totmmf) +//mean length of core CGHF=0.75m +ml=0.75 //as above +//since the central limb and magnetic core are in parallel they have same mmf that is 223.86AT +h2=totmmf/ml //magnetic intensity in CGHF +printf("\n h2= %0.5f AT",h2) +flux2=B*area +printf("\n flux2= %0.5f Wb",flux2) +totflux=fluxa+flux2 //Wb +Bfabc=totflux/area //flux density in magnetic core fabc in Wb/m**2 +H=3000 //AT/m +totmmffabc=H*ml //total mmf in fabc in AT +printf("\n total mmf in fabc= %0.5f Wb/m**2",totmmffabc) +totmmfm=totmmffabc+totmmf //total mmf in magnetic core in AT +printf("\n totmmfm= %0.5f AT",totmmfm) +N=500 +I=totmmfm/N //The required current to set up flux in air gap +printf("\n The total current required to set up flux in air gap is= %0.5f A",I) diff --git a/3835/CH3/EX3.8/Ex3_8.sce b/3835/CH3/EX3.8/Ex3_8.sce new file mode 100644 index 000000000..684507d7c --- /dev/null +++ b/3835/CH3/EX3.8/Ex3_8.sce @@ -0,0 +1,34 @@ +clear +// +//given +r1=3.98*10**6 //reluctance of air gap in AT/Wb and the value is same for r2 +r3=5.97*10**6 //reluctance of air gap in AT/Wb +//assume that current of 1A flows through 150 turns coil,for assumed directions of fluxes application of mesh current leads to matrix equations that can be simplified to: +//[flux1 flux2]=[2.36 1.41]*10**-5 Wb +//The self inductance and mutual inductance are obtained as follows: +n1=150 //number of turns +i1=1 //A +flux1=2.36*10**-5 //Wb +l1=(n1*flux1)/i1 //self inductance +printf("\n l1= %0.3f mH",l1) +n2=200 //number of turns +flux2=1.41*10**-5 +m12=(n2*flux2)/i1 //mutual inductance +printf("\n m12= %0.3f mH",m12) +//assume that 1A of current flows through 200 turns coil +//The self inductance of the coil is determined as above using the matrix and the result is as follows +//[flux1 flux2]=[1.89 3.14]*10**-5 Wb +//Hence self and mutual inductance are computed as follows +n2=200 //number of turns +flux2=3.14*10**-5 //Wb +i2=1 //A +l2=(n2*flux2)/i2 //self inductance +printf("\n l2= %0.3f mH",l2) +flux1=1.89*10**-5 +m21=(n1*flux1)/i2 //mutual inductance +printf("\n m21= %0.3f mH",m21) +//case b +//When the air gap l3 is closed the reluctance of the limb is zero since the permeability of the magnetic material is infinity.Thus,the limb behaves like short circuit and the entire flux passes through it.Thus,the flux linking 200 turns coil is zero and mutual inductance is zero +//case 3 +W=((3.5)/2)+((6.3)/2)+2.8 //work equation in joules +printf("\n Work done= %0.5f J",W) diff --git a/3835/CH3/EX3.9/Ex3_9.sce b/3835/CH3/EX3.9/Ex3_9.sce new file mode 100644 index 000000000..012cde006 --- /dev/null +++ b/3835/CH3/EX3.9/Ex3_9.sce @@ -0,0 +1,29 @@ +clear +// +//given +B=0.8 //Wb/m**2 +A=25*10**-4 //m**2 +flux=20*10**-4 //Wb +l=3.14*40*10**-2 //m +f=2000*3.14 //AT +n=800 //number of turns +//case a +i=f/n //A exciting current +printf("\n i= %0.3f A",i) +l=(n*flux)/i //self inductance in H +printf("\n l= %0.5f H",l) +//case b +fluxa=20*10**-4 //Wb +gap=1*10**-2 +u1=4*3.14*10**-7 +rair=gap/(u1*A) //reluctance of air in AT/Wb +printf("\n rair= %0.5f AT/Wb",rair) +fair=rair*flux //mmf for air gap in AT +printf("\n fair= %0.5f AT",fair) +fcore=6233.18 //AT--> 5000*((0.4*3.14)-0.01)=6233.18 +totmmf=fcore+fair +printf("\n total mmf= %0.5f AT",totmmf) +I=totmmf/n //A exciting current +//self inductance +L=(n*flux)/I +printf("\n L= %0.5f H",L) |