## Contents

clc;
clear;
format compact;


## Example Cable Calculation d.c. - version

For a.c. cables example please refer to Calculating Cable Performance (a.c.).

% <https://mycableengineering.com/knowledge-base/estimating-cable-life
% Estimating Cable Life>.

%
% This calculation provides a simple look at how to calculate various
% cable parameters for a d.c. cable.
% The inputs are assumed to be the cable sustained
% current capacity (at maximum operating temperature), the maximum
% operating temperature, ambient temperature, system voltage, cable
% current, length, cable resistance and reactance.
%
% The calculations are based on a directly buried 16 mm2, Leoni
% BETAflam Solar cable (copper conductor), and using data from the
% manufacturer's % catalogue. Cable is installed in air,


## Cable Inputs - Manufacturer's Information

Uo = 720;           % V. System voltage (a.c.).
Iz = 152;           % A. Sustained current capacity (surface, w/o contact).
Tc = 90;            % °C. Maximum operating temperature of the cable.
R20 = 1.24e-3;      % ohm/m. Resistance of the cable at 20 °C.
X = 0;              % ohm/m. d.c. cable have no reactance.


## Cable Inputs - Operating Conditions

Ib = 132;   % A. Cable design current.
L = 56;     % m. Cable length.
Ta = 30;    % °C. Ambient temperature where the cable is installed.


## Current Capacity

In this instance we read the current capacity from the manufacturer's tables. This may also be calculated using standards, or thermal methods, which are beyond the scope of this example.

## Resistance at Operating Temperature

The resistance of the cable depends on the operating temperature. And, the operating temperature of the cable changes depending on the magnitude of current flowing within the cable. The temperature adjusted value of resistance is used in all calculations. To calculate the conductor operating temperature we use:

t = (Ib/Iz) * (Tc - Ta) + Ta

The cable resistance at 20 °C is adjusted for operating temperature. The amount of adjustment is dependent on the temperature coefficient of resistance (a20), with typical values being 3.93e-3 for copper conductors and 4.03e-3 for aluminium conductors. The adjusted resistance is given by:

R = R20 * [1+a20*(t-20)]

a20 = 3.93e-3;  % Temperature coefficent aluminium cable.
t = (Ib/Iz)*(Tc - Ta)+Ta;
R = R20*(1+a20*(t-20));
disp(['Conductor operating temperature = ', num2str(t, '%0.1f' ), ' °C']);
disp(['Temperature adjusted resistance = ', num2str(R, '%0.3e' ), ' ohm/m']);

Conductor operating temperature = 82.1 °C
Temperature adjusted resistance = 1.543e-03 ohm/m


## Voltage Drop

For voltage drop we use the formulae give in CENELEC technical report CLC/TR 50480. This give the voltage drop, in percent as:

dU = 100* b * R * Ib / Uo

The factor 'b' is taken as 2 for d.c. systems and 1 for three phase systems.

b = 2;  % d.c. system.
dU = 100*b*R*Ib*L/Uo;
disp(['Voltage drop = ', num2str(dU,'%0.2f'), ' %']);

Voltage drop = 3.17 %


## Power Loss

Calculation of the power loss within the cable is given by:

P = Ib^2 * R * n (i.e. the current squared, times the resistance)

Again with our resistance and reactance being per unit length, we need to multiple by the cable length to get the final power loss. Factor 'n' relates to the type of circuit/number of conductors. For a d.c. circuit we have two conductors ( n=2 )

n = 2;
P = Ib^2*R*n*L;
disp(['Power loss = ', num2str(P,'%0.0f'), ' W']);

Power loss = 3010 W


## Life Expectancy

Cable life expectancy can be estimated using the Arrhenius Equation. This is statistical estimation and care needs to be taken in interpreting the results. The estimated life, k in hours is given by:

k = Ae^(E/Rt)

In this equation, A is a factors, E the activation energy, and R the Boltzmann constant. For more details see on this equations, see Estimating Cable Life. Alternatively the Arrhenius Equation can be presented more practically as a logarithmic straight line (with slope and intercept coefficients of m and b):

ln(k) = m*t + b

For our cable, with XLPE insulation m and b are found to be -0.0693 and 18.5678 respectively. Taking the inverse ln of k, yields the life in hours.

m = -0.0693;    % Arrhenius slope factor (XLPE)
b = 18.5678;    % Arrhenius intercept factor (XLPE)
k = exp(1)^(m*t+b); ky = k/(24*365);
disp(['Expeted Life = ', num2str(k,'%0.0f'), ' hours (', num2str(ky,'%0.0f'), ' years)']);

Expeted Life = 391564 hours (45 years)