3 Engineering Example 1

3.1 Detecting a train on a track

Introduction

One means of detecting trains is the ‘track circuit’ which uses current fed along the rails to detect the presence of a train. A voltage is applied to the rails at one end of a section of track and a relay is attached across the other end, so that the relay is energised if no train is present, whereas the wheels of a train will short circuit the relay, causing it to de-energise. Any failure in the power supply or a breakage in a wire will also cause the relay to de-energise, for the system is Òfail safeÓ. Unfortunately, there is always leakage between the rails, so this arrangement is slightly complicated to analyse.

Problem in words

A 1000 m track circuit is modelled as ten sections each 100 m long. The resistance of 100 m of one rail may be taken to be 0.017 ohms, and the leakage resistance across a 100 m section taken to be 30 ohms. The detecting relay and the wires to it have a resistance of 10 ohms, and the wires from the supply to the rail connection have a resistance of 5 ohms for the pair. The voltage applied at the supply is 4 V . See diagram below. What is the current in the relay?

Figure 1

No alt text was set. Please request alt text from the person who provided you with this resource.

Mathematical statement of problem

There are many ways to apply Kirchhoff’s laws to solve this, but one which gives a simple set of equations in a suitable form to solve is shown below. i 1 is the current in the first section of rail (i.e. the one close to the supply), i 2 , i 3 , i 10 , the current in the successive sections of rail and i 11 the current in the wires to the relay. The leakage current between the first and second sections of rail is i 1 i 2 so that the voltage across the rails there is 30 ( i 1 i 2 ) volts. The first equation below uses this and the voltage drop in the feed wires, the next nine equations compare the voltage drop across successive sections of track with the drop in the (two) rails, and the last equation compares the voltage drop across the last section with that in the relay wires.

30 ( i 1 i 2 ) + ( 5.034 ) i 1 = 4 30 ( i 1 i 2 ) = 0.034 i 2 + 30 ( i 2 i 3 ) 30 ( i 2 i 3 ) = 0.034 i 2 + 30 ( i 3 i 4 ) 30 ( i 9 i 10 ) = 0.034 i 10 + 30 ( i 10 i 11 ) 30 ( i 10 i 11 ) = 10 i 11

These can be reformulated in matrix form as A i ̲ = v ̲ , where v ̲ is the 11 × 1 column vector with first entry 4 and the other entries zero, i ̲ is the column vector with entries i 1 , i 2 , , i 11 and A is the matrix

A = 35.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 60.034 30 0 0 0 0 0 0 0 0 0 30 40

Find the current i 1 in the relay when the input is 4 V , by Gaussian elimination or by performing an L-U decomposition of A .

Mathematical analysis

We solve A i ̲ = v ̲ as above, although actually we only want to know i 11 . Letting M be the matrix A with the column v ̲ added at the right, as in Section 30.2, then performing Gaussian elimination on M , working to four decimal places gives

M = 35.0340 30.0000 0 0 0 0 0 0 0 0 0 4.0000 0 34.3447 30.0000 0 0 0 0 0 0 0 0 3.4252 0 0 33.8291 30.0000 0 0 0 0 0 0 0 2.9919 0 0 0 33.4297 30.0000 0 0 0 0 0 0 2.6532 0 0 0 0 33.1118 30.0000 0 0 0 0 0 2.3810 0 0 0 0 0 32.8534 30.0000 0 0 0 0 2.1572 0 0 0 0 0 0 32.6396 30.0000 0 0 0 1.9698 0 0 0 0 0 0 0 32.4601 30.0000 0 0 1.8105 0 0 0 0 0 0 0 0 32.3077 30.0000 0 1.6733 0 0 0 0 0 0 0 0 0 32.1769 30.0000 1.5538 0 0 0 0 0 0 0 0 0 0 12.0296 1.4487

from which we can calculate that the solution i ̲ is

i ̲ = 0.5356 0.4921 0.4492 0.4068 0.3649 0.3234 0.2822 0.2414 0.2008 0.1605 0.1204

so the current in the relay is 0.1204 amps, or 0.12 A to two decimal places.

You can alternatively solve this problem by an L-U decomposition by finding matrices L and U such that M = L U . Here we have

L = 1.0000 0 0 0 0 0 0 0 0 0 0 0.8563 1.0000 0 0 0 0 0 0 0 0 0 0 0.8735 1.0000 0 0 0 0 0 0 0 0 0 0 0.8868 1.0000 0 0 0 0 0 0 0 0 0 0 0.8974 1.0000 0 0 0 0 0 0 0 0 0 0 0.9060 1.0000 0 0 0 0 0 0 0 0 0 0 0.9131 1.0000 0 0 0 0 0 0 0 0 0 0 0.9191 1.0000 0 0 0 0 0 0 0 0 0 0 0.9242 1.0000 0 0 0 0 0 0 0 0 0 0 0.9286 1.0000 0 0 0 0 0 0 0 0 0 0 0.9323 1.0000

and

U = 35.0340 30.0000 0 0 0 0 0 0 0 0 0 0 34.3447 30.0000 0 0 0 0 0 0 0 0 0 0 33.8291 30.0000 0 0 0 0 0 0 0 0 0 0 33.4297 30.0000 0 0 0 0 0 0 0 0 0 0 33.1118 30.0000 0 0 0 0 0 0 0 0 0 0 32.8534 30.0000 0 0 0 0 0 0 0 0 0 0 32.6395 30.0000 0 0 0 0 0 0 0 0 0 0 32.4601 30.0000 0 0 0 0 0 0 0 0 0 0 32.3076 30.0000 0 0 0 0 0 0 0 0 0 0 32.1768 30.0000 0 0 0 0 0 0 0 0 0 0 12.0295

Therefore U i ̲ = 4.0000 3.4240 2.9892 2.6514 2.3783 2.1547 1.9673 1.8079 1.6705 1.5519 1.4464  and hence i ̲ = 0.5352 0.4917 0.4487 0.4064 0.3644 0.3230 0.2819 0.2411 0.2006 0.1603 0.1202 and again the current is found to be 0.12 amps.

Mathematical comment

You can try to solve the equation A i ̲ = v ̲ by Jacobi or Gauss-Seidel iteration but in both cases it will take very many iterations (over 200 to get four decimal places). Convergence is very slow because the norms of the relevant matrices in the iteration are only just less than 1. Convergence is nevertheless assured because the matrix A is diagonally dominant.