The operation of the RL circuit supplied by a half-wave controlled rectifier is illustrated in this program. Let vs(t) be defined to E * sin (wt). At wt = 0, the current through the circuit is zero. As wt becomes >0, vs becomes positive. If a diode instead of an SCR has been used, the diode would start conduction at wt = 0. With an SCR, the conduction does not start till the SCR is triggered. Let the SCR be triggered when wt = a. Then a is called the firing angle and the SCR continues to conduct.

When wt = p, the source becomes zero, but at this instant, the current through the circuit is not zero and there is some energy stored in the inductor. When vs becomes negative, the current through the circuit would not become zero suddenly because of the inductor. The inductor acts as a source and keeps the SCR forward-biased till the energy stored in the inductor becomes zero. Let the current through the circuit become zero at wt = b and the value of b > p. For b < wt < 2p, the current through the circuit is zero. Then

L (di/dt) + Ri = E * sin (wt ) for a < wt < b,

i(a) = i(b) = 0 , i(wt) = 0 for 0 < wt < a and

i(wt) = 0 for b < wt < 2p .

The complementary solution is;

i(wt) = A *exp [- (R(wt - a) / wL ) ], where A is a constant to be evaluated.

The particular solution is:

i(wt) = [E/Z] sin (wt - f ),
where f = tan ^{-1} (wL/R) , Z^{2}
= R^{2} + (wL)^{2} .

The total solution is then:

i(wt) = [E/Z] sin (wt - f ) + A *exp [- (R/wL *(wt - a) ].

Since i(a) = 0 , A = [E/Z * sin (f - a)]. The solution is illustrated below.

First specify the values of parameters as shown below.

E :=340 V R:= 10 W wL := 10 W f := atan(wL/R) f = 0.785 rad

Z := sqrt [(R^{2}
+ (wL)^{2)}]
Z = 141.42 W

Let firing angle a be: a := p/6 rad

Next compute the coefficient
A as shown below. Then define a range variable n to find solution at each
degree within a cycle and express the angle within a cycle in radians.
Find the load current. If the angle is less than the firing angle, assign
it zero value The expression for load current calculated thus would be
negative if wt

exceeds the extinction
angle, b.
The extinction angle is that angle at which load current becomes zero.
Create an array for current, called I_{n}. This array has the same
value as Cur_{n} if Cur_{n} >0.0. Otherwise the other
elements of array I_{n} are set to zero.

The voltage across
the SCR is zero when I_{n} > 0. Otherwise it is equal to the applied
voltage. Note that when wt < a
, the voltage across the SCR is positve. Then the SCR blocks forward conduction.
Next find the extinction angle. Assign a guess value first.

Next the average value of load current is found out. The RMS load current is evaluated after that.

A = 6.222 Amp d = 3.936 rad

extinction_angle = 225.489 deg

The extinction angle is computed to be 225.49^{o}.

Curavg = 8.339 Amp

CurRMS = 12.851 Amp

**The load current waveform**

**The voltage across the SCR**

**The Plot of Output Voltage**

Vout_{n} : = if (n > 0.001, E.Sin(q_{n}),
0.0)