summaryrefslogtreecommitdiff
path: root/3523/CH16
diff options
context:
space:
mode:
Diffstat (limited to '3523/CH16')
-rw-r--r--3523/CH16/EX16.7.1/Ex16_1.sce28
-rw-r--r--3523/CH16/EX16.7.10/Ex16_10.sce27
-rw-r--r--3523/CH16/EX16.7.11/Ex16_11.sce19
-rw-r--r--3523/CH16/EX16.7.2/Ex16_2.sce30
-rw-r--r--3523/CH16/EX16.7.3/Ex16_3.sce30
-rw-r--r--3523/CH16/EX16.7.4/Ex16_4.sce45
-rw-r--r--3523/CH16/EX16.7.5/Ex16_5.sce16
-rw-r--r--3523/CH16/EX16.7.6/Ex16_6.sce30
-rw-r--r--3523/CH16/EX16.7.7/Ex16_7.sce29
-rw-r--r--3523/CH16/EX16.7.8/Ex16_8.sce21
-rw-r--r--3523/CH16/EX16.7.9/Ex16_9.sce34
11 files changed, 309 insertions, 0 deletions
diff --git a/3523/CH16/EX16.7.1/Ex16_1.sce b/3523/CH16/EX16.7.1/Ex16_1.sce
new file mode 100644
index 000000000..f60e630df
--- /dev/null
+++ b/3523/CH16/EX16.7.1/Ex16_1.sce
@@ -0,0 +1,28 @@
+clear all
+clc
+close
+
+iload=5*1e-3;//Load current in A
+
+//Capacitances of Cockcroft-Waltobn type voltage doubler in F
+C1=0.01*1e-6;
+C2=0.05*1e-6;
+
+f=50;//frequency in Hz
+Vs=100*1e3//Supply voltage in V
+
+//Ripple voltage in volt
+dv=iload/(C2*f)
+printf('Ripple voltage in V %f',dv)
+
+//Voltage drop in Volt
+Vdrop=iload/f*(1/C1+1/(2*C2))
+printf('Voltage drop in V %f',Vdrop)
+
+//Average output voltage
+V_av=2*sqrt(2)*Vs-Vdrop//in V
+printf('Avarage voltage in V %f',V_av)
+
+//Ripple factor
+RF=Vdrop/(2*sqrt(2)*Vs)*100//in percentage
+printf('Ripple voltage in percentage %f',RF)
diff --git a/3523/CH16/EX16.7.10/Ex16_10.sce b/3523/CH16/EX16.7.10/Ex16_10.sce
new file mode 100644
index 000000000..6d675c739
--- /dev/null
+++ b/3523/CH16/EX16.7.10/Ex16_10.sce
@@ -0,0 +1,27 @@
+clear all
+clc
+close
+
+n=12;//no ofstage
+C1=0.125*1e-6;//Each stage capacitor in F
+C2=1000e-12;//Load capacitance in F
+R1=70;//Front resistance in ohm
+R2=400;//Tail resistance in ohm
+
+R1T=R1*n;
+R2T=R2*n;
+C1T=C1/n;
+
+theta=sqrt(C1T*C2*R1T*R2T);
+
+eta=1/(1+(1+R1T/R2T)*C2/C1T);
+
+alpha=R2T*C1T/(2*eta*theta);
+
+//Wavetail time in us
+T2=7*theta*1e6;
+printf('Wave tail time in us %f',T2)
+
+//Wave front time in us
+T1=T2/25;
+printf('Wave front time in us %f',T1)
diff --git a/3523/CH16/EX16.7.11/Ex16_11.sce b/3523/CH16/EX16.7.11/Ex16_11.sce
new file mode 100644
index 000000000..4617fb4a1
--- /dev/null
+++ b/3523/CH16/EX16.7.11/Ex16_11.sce
@@ -0,0 +1,19 @@
+clear all
+clc
+close
+
+C = 8*10^-6;//in Farad
+L = 8*10^-6;//in Henry
+V = 25*10^3;//in V
+T1 = 8;//in us time for the first peak
+ohmega = 0.02*10^-6;//in sec^-1
+R = sqrt((4*L/C)-(4*L^2*ohmega^2));
+printf("resistance is %f ohms \n",R)
+gama = R/(2*L);//in sec^-1
+printf("parameter gama is %f sec^-1 \n",gama)
+
+//Now eq for generated impulse pulse is I(t)= 156.25*10^3exp(-12.3*10^4t)sin(0.02*10^6t)A
+
+
+
+
diff --git a/3523/CH16/EX16.7.2/Ex16_2.sce b/3523/CH16/EX16.7.2/Ex16_2.sce
new file mode 100644
index 000000000..8970b928b
--- /dev/null
+++ b/3523/CH16/EX16.7.2/Ex16_2.sce
@@ -0,0 +1,30 @@
+
+clear all
+clc
+close
+
+iload=5*1e-3;//Load current in A
+
+//Capacitances of Cockcroft-Waltobn type voltage tripler in F
+C1=0.01*1e-6;
+C2=0.05*1e-6;
+C3=0.10*1e-6;
+
+f=50;//frequency in Hz
+Vs=100*1e3//Supply voltage in V
+
+//Ripple voltage in V
+dv=iload/f*(2/C1+1/C3)
+printf('Ripple voltage in V %f',dv)
+
+//Voltage drop in V
+Vdrop=iload/f*(1/C2+1/C1+1/(2*C3))
+printf('Voltage drop in V %f',Vdrop)
+
+//Average output voltage in V
+V_av=3*sqrt(2)*Vs-Vdrop
+printf('Avarage voltage in V %f',V_av)
+
+//Ripple factor in percentage
+RF=Vdrop/(3*Vs*sqrt(2))*100
+printf('Ripple voltage in percentage %f',RF)
diff --git a/3523/CH16/EX16.7.3/Ex16_3.sce b/3523/CH16/EX16.7.3/Ex16_3.sce
new file mode 100644
index 000000000..890ebaa64
--- /dev/null
+++ b/3523/CH16/EX16.7.3/Ex16_3.sce
@@ -0,0 +1,30 @@
+clear all
+clc
+close
+
+Vs=200*1e3//Supply voltage
+f=50//Frequency in Hz
+n=12//Number of stages
+
+C=0.15*1e-6//Each stage capacitance in F
+iload=5*1e-3//Load current in A
+
+//Ripple voltage in V
+dv=iload/(f*C*2)*n*(n+1)
+printf('Ripple voltage in V %f',dv)
+
+//Voltage drop in V
+Vdrop=iload/(f*C)*(2*n^3/3+n^2/2-n/6+n*(n+1)/4)
+printf('Voltage drop in V %f',Vdrop)
+
+//Average output voltage in V
+V_av=2*n*sqrt(2)*Vs-Vdrop
+printf('Avarage voltage in V %f',V_av)
+
+//Ripple factor in percentage
+RF=Vdrop/(2*n*Vs*sqrt(2))*100
+printf('Ripple voltage in percentage %f',RF)
+
+//Otimum number of stages
+nopt=sqrt(sqrt(2)*f*C*Vs/iload)
+printf('Optimum number of stgaes for minimum voltage drop %f',int(nopt))
diff --git a/3523/CH16/EX16.7.4/Ex16_4.sce b/3523/CH16/EX16.7.4/Ex16_4.sce
new file mode 100644
index 000000000..9f8d862a9
--- /dev/null
+++ b/3523/CH16/EX16.7.4/Ex16_4.sce
@@ -0,0 +1,45 @@
+clear all
+clc
+close
+
+f=50;//Power frequency
+xl=8/100;//leakage reactance
+r=3.5/100;//resistance
+Vc=500;//Charging voltage in kV
+Ic=4;//Charging current in A
+capc=100;//kVA rating of transformer
+vhigh=250;//Voltage rating of secondary of transformer in kV
+vlow=220;//Voltage rating of primary of transformer in V
+
+//Reactance of cable in kiloohm
+Xc=Vc/Ic
+
+//Leakage recatance of transformer in kiloohm
+XL=xl*(vhigh^2/capc)
+
+//Additional series inductance
+xh=Xc-XL;
+
+//Inductance of the required series inductor in Henry
+L=xh/(2*%pi*f)*1e3;
+printf('Inductance of the required series inductor in %f Henry \n',L)
+
+//Total circuit resistance in kiloohm
+R=r*(vhigh^2/capc)
+
+//The maxium current can be supplied by transformer in A
+I=capc/vhigh;
+
+Vsec = I*R;
+printf("exciting voltage on the transformer secondary %f kV \n",Vsec)
+
+//Exciting voltage of secondary of transformer in kV
+Vexsec=I*R;
+
+//Input voltage to primary of transformer in V
+Vin=Vexsec*1e3*vlow/(vhigh*1e3);
+printf('Input voltage to primary of transformer in %f V \n',Vin)
+
+//Input power to transformer in kW
+Pin=Vin*capc/vlow
+printf('Input power to primary of transformer in %f kW \n',Pin)
diff --git a/3523/CH16/EX16.7.5/Ex16_5.sce b/3523/CH16/EX16.7.5/Ex16_5.sce
new file mode 100644
index 000000000..ce5d76375
--- /dev/null
+++ b/3523/CH16/EX16.7.5/Ex16_5.sce
@@ -0,0 +1,16 @@
+clear all
+clc
+close
+
+u=10//speed of belt in m/s
+w=0.1//width of the belt in m
+rhos=0.5*1e-6//surface charge density on the belt in C/m^2
+Rleak=1e14//Resistanc ein ohm
+
+//Charging current in A
+I=rhos*u*w
+printf('Charging current in uA %f',I*1e6)
+
+//Potentail difference between the dome and the base in V
+V=I*Rleak
+printf('Potentail difference between the dome and the base in MV is %f',V/1e6)
diff --git a/3523/CH16/EX16.7.6/Ex16_6.sce b/3523/CH16/EX16.7.6/Ex16_6.sce
new file mode 100644
index 000000000..0bfdbc538
--- /dev/null
+++ b/3523/CH16/EX16.7.6/Ex16_6.sce
@@ -0,0 +1,30 @@
+clear all
+clc
+close
+
+C1=0.125*10^-6;//in Farad
+C2=1*10^-9;//in Farad
+R1=360;//in ohms
+R2=544;//in ohms
+theta = sqrt(C1*C2*R1*R2);//in usec
+n = 1/[1+(1+(R1/R2))*(C2/C1)];
+alpha = (R2*C1)/(2*theta*n);
+printf("theta parameter of wave eq %f us \n",theta*10^6)
+printf("n the parameter of circuit eq %f \n",n)
+printf("alpha parameter of circuit eq %f \n",alpha)
+T2 = 10.1*theta;//duration of lightning impulse pulse in us
+T1 = T2/45;//duration of lightning impulse pulse in us
+printf("duration of lightning impulse pulse %f us \n",T2*10^6)
+printf("duration of lightning impulse pulse %f us \n",T1*10^6)
+//answer in the book for T1 is wrong
+
+T = T1/T2;
+printf("generated lighting impulse is %f us \n",T)
+alpha1 = [alpha-sqrt((alpha^2)-1)]/theta;//in us^-1
+alpha2 = [alpha+sqrt((alpha^2)-1)]/theta;//in us^-1
+printf("aplha1 parameter of wave eq is %f us^-1 \n",alpha1*10^-6)
+printf("aplha1 parameter of wave eq is %f us^1 \n",alpha2*10^-6)
+
+//answer in the book is slightly different
+// Now eq of waveform of generated pulse is e(t)=99.75(e^-0.015t - e^-2.77t)
+
diff --git a/3523/CH16/EX16.7.7/Ex16_7.sce b/3523/CH16/EX16.7.7/Ex16_7.sce
new file mode 100644
index 000000000..69044cc5c
--- /dev/null
+++ b/3523/CH16/EX16.7.7/Ex16_7.sce
@@ -0,0 +1,29 @@
+clear all
+clc
+close
+
+C1=0.125*10^-6;//in Farad
+C2=1*10^-9;//in Farad
+R1=360;//in ohms
+R2=544;//in ohms
+theta = sqrt(C1*C2*R1*R2);//in usec
+n = 1/[1+(R1/R2)+(C2/C1)];
+alpha = (R2*C1)/(2*theta*n);
+printf("theta parameter of wave eq %f us \n",theta*10^6)
+printf("n the parameter of circuit eq %f \n",n)
+printf("alpha parameter of circuit eq %f \n",alpha)
+T2 = 16.25*theta;//duration of lightning impulse pulse in us
+T1 = T2/120;//duration of lightning impulse pulse in us
+printf("duration of lightning impulse pulse %f us \n",T2*10^6)
+printf("duration of lightning impulse pulse %f us \n",T1*10^6)
+//answer in the book for T1 is wrong
+
+T = T1/T2;
+printf("generated lighting impulse is %f us \n",T)
+alpha1 = [alpha-sqrt((alpha^2)-1)]/theta;//in us^-1
+alpha2 = [alpha+sqrt((alpha^2)-1)]/theta;//in us^-1
+printf("aplha1 parameter of wave eq is %f us^-1 \n",alpha1*10^-6)
+printf("aplha1 parameter of wave eq is %f us^1 \n",alpha2*10^-6)
+
+// Now eq of waveform of generated pulse is e(t)=60.2(e^-0.0088t - e^-4.62t)
+
diff --git a/3523/CH16/EX16.7.8/Ex16_8.sce b/3523/CH16/EX16.7.8/Ex16_8.sce
new file mode 100644
index 000000000..b9b5b921d
--- /dev/null
+++ b/3523/CH16/EX16.7.8/Ex16_8.sce
@@ -0,0 +1,21 @@
+clear all
+clc
+close
+
+//Elements of circuits
+C1=0.125*1e-6;//in F
+C2=1e-9;//in F
+
+T1=250*1e-6;
+T2=2500*1e-6;
+alpha=4;
+theta=T2/6;
+
+X=(1+C2/C1)*1/alpha^2;
+R1=alpha*theta/C2*(1-sqrt(1-X));//in ohm
+
+R2=alpha*theta/(C1+C2)*(1+sqrt(1-X));//in ohm
+
+//Circuit efficiency
+eta=1/(1+(1+R1/R2)*C2/C1)
+printf('Circuit efficiency %f',eta)
diff --git a/3523/CH16/EX16.7.9/Ex16_9.sce b/3523/CH16/EX16.7.9/Ex16_9.sce
new file mode 100644
index 000000000..9d656858b
--- /dev/null
+++ b/3523/CH16/EX16.7.9/Ex16_9.sce
@@ -0,0 +1,34 @@
+clear all
+clc
+close
+
+n=8;//no ofstage
+C1=0.16*1e-6;//Each stage capacitor in F
+C2=1e-9;//Load capacitance in F
+T2=50*1e-6;
+T1=1.2*1e-6;
+Vch=120;//Charging voltage in kV
+
+//Total capacitance in F
+CT=C1/n;
+
+alpha=6.4;
+theta=T2/9.5;
+
+X=(1+C2/C1)/alpha^2;
+R1=alpha*theta/C2*(1-sqrt(1-X));//in ohm
+
+R2=alpha*theta/(CT+C2)*(1+sqrt(1-X));//in ohm
+//Perstage shaping resistance in ohm
+printf('Perstage shaping resistance in %f ohm',R1/n)
+
+Vdc=n*Vch;
+eta=1/(1+(1+R1/R2)*C2/CT)
+
+//Maximum output voltage
+Vmax=eta*Vdc;
+printf('Maxium output voltage in %f kV',Vmax)
+
+//Energy rating in J
+E=0.5*CT*(Vdc*1e3)^2;
+printf('Energy rating in %f J',E)
generated by cgit (git 2.34.1) at 2024-12-19 12:39:10 +0000