A Full Frequency-Dependent Cable Model for the Calculation of Fast Transients

The calculation of frequency-dependent cable parameters is essential for simulations of transient phenomena in electrical power systems. The simulation of transients is more complicated than the calculation of currents and voltages in the nominal frequency range. The model has to represent the frequency dependency and the wave propagation behavior of cable lines. The introduced model combines an improved subconductor method for the determination of the frequency-dependent parameters and a PI section wave propagation model. The subconductor method considers the skin and proximity effect in all conductors for frequency ranges up to few megahertz. The subconductor method method yields accurate results. The wave propagation part of the cable model is based on a cascaded PI section model. A modal transformation technique has been used for the calculation in the time domain. The frequency-dependent elements of the related modal transformation matrices have been fitted with rational functions. The frequency dependence of cable parameters has been reproduced using a vector fitting algorithm and has been implemented into an resistor-inductor-capacitor network (RLC network) for each PI section. The proposed full model has been validated with measured data.


Introduction
The safe and reliable operation of cable systems requires the analysis of transient behavior of such systems.Single core cables with a core and a shield are utilized prospectively.The first models for the calculations of wave propagation and their transient behavior were developed in the 1970s.The structure of the line models is divided into two main components, a module for calculation of frequency-dependent parameters and the wave-propagation model.
For calculation of the line parameters, Ametani's algorithm [1] (cable constants) uses the principal assumptions of Wedepohl [2].Ametani's algorithm uses a bundle of equations to determine the coupled cable impedances.The frequency-dependent earth return impedance is considered with an fictive earth conductor given by Carson [3] and Pollaczek [4].The skin effect within the cable core is implemented using Schelkunoff's equations [5].The proximity effect between the conductors and the skin effect in the shield are neglected.Morched et al. [6] also use principal assumptions similar to the ideas of Wedepohl [2] and Ametani [1].
The wave propagation models can be mainly classified into two categories: the models with distributed parameters and the models with lumped parameters.The models with distributed parameters are based on partial differential equations to describe the dependency of current and voltage on the time and distance.These differential equations are also known as telegraph equations [7].d'Alembert solved the telegraph equations [8].Using the Bergeron method [9], d'Alembert's solution is represented with forward and backward waves, where the line parameters are assumed to be frequency independent.In addition, the losses are neglected.
Dommel [10] developed the method of Bergeron by considering the line losses using lumped resistances.Snelson [11] was able to further develop Dommel's model by taking the frequency dependency of the line parameters into account.Based on Snelson's model, Marti J. [12] presented a model for overhead lines.In this model, a modal transformation has been used to simplify the modelling.The modal transformations are performed using modal transformation matrices, whose parameters are generally frequency dependent.However, Marti J. used constant (frequency independent) modal transformation matrices.The frequency dependency of the modal transformation matrices are considered later in the model of Marti L. [13].In order to avoid the difficulties by using frequency-dependent modal transformation matrices, Noda et al. [14] have developed the line model directly in the phase domain.Morched et al. [6] have further developed this model and introduced a universal line model.The model can be used for overhead lines and cables.In order to simplify the calculations, the universal line model uses a vector fitting algorithm [15,16], which approximates some of the line parameters (surge admittance and propagation constant) by rational functions.
Another approach (alternatively to the vector fitting algorithm) for the approximation of the line parameters is based on the algorithm developed in Noda [17].Noda has shown that this algorithm can be also applied to both overhead lines [18] and cables [19].The universal line model gives good results.However, it sometimes produces numerical instabilities.This happens by approximation of the propagation constant and especially in the case of short cables, since the associated time delays become closer, leading to poor fitting.This problem was overcome in [20] by introducing a new method to calculate the time delays.
A second type for modeling of the line is used by lumping of the line parameters.The lumped parameters can be combined to build a PI section [10].The behavior of the models with distributed parameters is approximated by connecting several PI sections in series, which results in a simplified representation of a line model with lumped parameters.The number of PI sections depends on the simulated cable length, the expected highest frequency in the investigated transient and the accepted calculation error.The principle of the model with lumped parameters or with PI sections is relatively simple compared to the models with distributed parameters.However, in the traditional PI section model the line parameters are calculated only at a specified frequency point.In other words, the frequency dependency of the line parameters is neglected.Therefore, the results of this model will only be valid for steady state simulations.For simulations of transients, the traditional PI section model has to be further developed [21].
In this paper, the main steps for the development of a full frequency-dependent three-phase PI section cable model (3PPI model) is introduced.The frequency-dependent impedances are determined with an improved sub-conductor method considering the skin and the proximity effect in the core, the shield and the earth.The frequency dependence of the cable parameters is implemented through an RLC (resistor-inductor-capacitor) network with a defined number of PI sections.To decouple the cable system, a modal transformation technique has been used.The calculation time necessary to work with a large number of PI sections is acceptable.For the simulation in the time domain, a recursive convolution technique has been used.The differential equations to describe the voltages and currents along the cable are of the first order, which can be easily solved as state space equations.

Fundamentals
For any inductive loop, the mutual inductance M of a conductor configuration characterizes the magnetic influencing of adjacent circuits on each other.The mutual inductance between two conductor loops M and N as shown in Figure 1 can be determined by Equation (1) (Kuepfmueller [22]).
Equation ( 1) is sufficiently accurate if the radii r i are negligibly small.Otherwise, the conductor cross-sectional area A must be taken into account.Maxwell [23] incorporated the cross-sectional area A through a mean geometric distance g.According to Bruederlink [24], for any two areas A 1 and A 2 , shown in Figure 2, their mean geometric distances g MN are given by Equation (2).The mean geometric distance of an area on itself g MM can be obtained by Equation (3).
For the characterization of conductor loops with common return conductor through earth, the analytical approaches of Carson [3] can be applied on cable systems over narrow frequency ranges only.To incorporate the return conductor into the model, it is reasonable to define a fictive hull cylinder with a radius R H . Figure 3 shows a simplified equivalent configuration introduced by Rees [25], illustrating the conductors m and n as well as the fictive hull cylinder which is the common return.The subconductor self impedance Z mm and the mutual impedances Z mn are given by Expression (4).
A resistive component R m is thereby included only in the self impedance.

Partial Subconductor Method
The principle of the partial subconductor method as published by Comellini [26] is based on the segmentation of the conductor cross-sectional area into sufficiently small subconductors so that the current density will be nearly homogeneous within each subconductor.The subconductors of the system are inductively coupled through the mutual inductance per unit length Z mn of all conductor loops.Therefore, it is possible to model skin and proximity effects in the conductors of the cable and in the earth.Figure 4 shows the principle of segmentation.The inductive coupling of current and voltage of all subconductors is characterized by the equation system of the subconductor-submatrices as specified in Equation (5) for three conductors A, B, and C.
It is possible to incorporate all conductors of the system into the algorithm, for instance any of the cable components such as the core conductor, shield, armor, earth, semiconducting layers and possible return conductors of the system.The conventional subconductor method realizes a "fixed" segmentation with arc segments.Dommel [27] discussed different shapes of segments, for instance circles, squares and so forth.The limitation to frequencies of f < 10 kHz is a significant disadvantage of the fixed segmentation.The reason for the limitation to frequencies is named by Lucas [28] described in detail in Schidt et al. [29,30].The limitation of the conventional sub-conductor method can be avoided by using frequency-dependent segmentation.

Frequency-Dependent Segmentation of Conductors
The segmentation of a cylindrical conductor is realized considering the skin and proximity effect.To determine the segmentation algorithm, the current density J( f , r) of the cylindrical conductor is analytically derived from Maxwell's equations as a function of the frequency f and the radius r.With increasing frequency f the skin effect grows stronger.At high frequencies, a significant current flow remains directly under the surface only.For instance, at a frequency f = 1000 Hz, 90% of the current I is concentrated in the outer 4.5 mm.Therefore, the cylindrical conductor is segmented across its radius r A as a function of the relative current density J/J max .It is assumed that a constant error of the current density distribution can be achieved by calculating the current density J( f , r) with the outer radius of the i-th skin r i directly dependent on the relative current density J/J max .The number of skins n can be set in accordance to the required accuracy.Nevertheless, it must be considered that the number of skins impacts the algorithm's computation time.The consideration of the proximity effect is covered by the skin subdivision into arc segments.
As an example, Figure 6 displays the allocation of the skins over the conductor's radius exemplary for five skins.The determination of the skin radii r 1 to r n follows a 1/n-increment of the relative current density J/J max .However, the subdivision of the conductors into skins as demonstrated in Figure 6 does not cover the proximity effect.This must be incorporated by another subdivision of the skins into arc segments (see Figure 7).Arnold [31] yields an analytical approach to calculate the proximity effect.Based on the error which occurs from the calculation of the subconductor impedances Z mm and Z mn , the segmentation of the skins can be defined.Nevertheless, calculating the mutual mean geometric distances g mn of two arc segments as indicated in Equation ( 2) and the self mean geometric distance g mm of one arc segment on itself as indicated in (3) causes a high numerical complexity.In order to simplify, Rees [25] and Brakelmann [32] introduced the centroid distance s of the corresponding arc segments into Expression (4), instead of using the mean geometric distance g.The replacement of the mean geometric distance g by the centroid distance s of the arc segments results in an error which cannot be ignored, but Rees [25] shows that at least for certain configurations it is essentially acceptable to use the centroid distance s.A numerical determination of the deviation between centroid distance s and mean geometric distance g is presented by Schmidt [29,30].Those calculations have shown that for the geometric factors k s = rs /d s smaller than 2, the relative deviation s mn /g mn is minimized.Figure 8 gives an example for the segmentation into three skins for a frequency of f = 50 Hz.

Earth Segmentation
The segmentation of the earth is based on the analytically calculated current distribution in the earth.The approach of Ruedenberg [33] can directly be applied to underground cable systems.The calculated current density J E ( f , r) is shown in Figure 9 and displays the relative current density J E /J max(E) as a function of the radius r at different frequencies between 50 Hz and 1 MHz.The specific earth resistance E is assumed to a typical value (in Central Europe) of 150 Ω m.The segmentation is again realized by analytical calculation of the conductor current density J C /J max(C) .For frequencies of f = 50 and 1000 Hz, Figure 10 shows how ten skins would be allocated.To provide a better overview, only every second skin is plotted in Figure 10.At the fictive hull cylinder with the radius R H the current density J must have declined to zero.As Dommel [34] was able to demonstrate, this condition can be fulfilled for a distance of ten times the penetration depth E .The segmentation of the skins is done according to the results of the errors computation discussed in Schmidt [30].

Shield Segmentation
The shields of single core cables are usually stranded rigid copper conductors illustrated in Figure 11.
cross section AD, single wire radius rD, single wire In this case, the frequency-dependent segmentation is not necessary since the skin effect is only effective in single wires.The self-impedance per unit length of the single wire with the cross section D can be calculated by Equation (6).

Calculation of Frequency-Dependent Impedances
Equation ( 4) describes the impedances Z of all subconductors in a certain area of all considered subconductors in a certain area and thus the inductive coupling between all subconductors.The reference conductor (or common return conductor) of all loops is the radius of the fictive hull cylinder R H .The complete impedance matrix of a system appears as represented in Equation ( 5) and contains the subconductor matrices of each individual phase.The rank of the impedance matrix Z is equal to the number of all subconductors in the enclosed area n all .A resistive component R , which characterizes the direct current resistance per unit length R m of the electrical subconductor m as a function of the specific resistance and the subconductor area A m , is merely included in the self impedances Z mm .The direct current resistances R m of all subconductors are comprised in a column vector R .The impedance matrix Z can be expressed according to Equation (7).
The indices of the centroid distance elements in Matrix (8) designate the respective subconductor.The subscripted index names the m-th subconductor of the current-carrying loop and the superscripted index names the n-th subconductor influenced by the current-carrying loop.
The subconductor matrices of the centroid distances s mn and the admittance matrix Y are symmetrical.Therefore, to calculate the centroid distances s, it is sufficient to determine the elements above the main diagonal and to reflect them at the main diagonal.For unknown currents i of the cable system and known voltages v , a system of equations (with reference to Equation ( 5)) can be obtained as given by Equation (9).
The admittance matrix Y is the inverse of the impedance matrix Z .The voltage across all subconductors in a conductor K is the same.It is: The sum of the subconductor admittance matrices Y mn and Y mm can be realized using basic matrix operations demonstrated in Relation (11).
The dimension (N × n all ) of the transformation matrix T is defined by the total number of conductors in the system N and the sum of all subconductors n all .Row K is associated with conductor K.All elements corresponding to the submatrix Y K are filled with 1, the remaining elements are filled with 0. For a seven-conductor system (3 × core, 3 × shield, 1 × earth) the transformation matrix T is of the form indicated by Equation (12).
The number of rows in the transformation matrix T according to Equation ( 12) is equal to the number of conductors N. The number of elements per row represents the number of subconductors n in conductor K.The reduced impedance matrix Z red can be obtained from the reduced admittance matrix Y red .

System Impedances
The reduced impedance matrix Z red (13) contains the submatrices of all conductors.It is useful to modify the matrix containing all submatrices to a (6 × 6) matrix.Expression (14) shows the structure of the reduced system with all submatrices for the core (C), shield (S) and earth (E).
The voltage vector v describes longitudinal voltages measured against the hull cylinder.The equation system ( 14) can be partitioned into impedances of the cable system and impedances of return conductors (15).
The reference potential is the potential of the fictive hull cylinder U H = 0.It can be expressed in Equation (16).(16) Equation ( 16) can be rearranged to obtain (17) for the current i  (17) Inserting Equation ( 17) into ( 16) results in Equation ( 18) for the voltage u p c .(18) The modified matrix Z mod has the rank of the sum of all cable conductors.The matrix a (6 × 6) matrix single core cable system (3 × core, 3 × shield).It is also possible to reduce the coefficient matrix of equation system (14) to a (3 × 3) matrix.In case of earthed shields on both cable ends, the return conductor submatrices of the partitioned matrix (15) contains the shield submatrices.The earth and shield influence is considered in all modified impedance matrices Z mod .In principle, this approach can also be realized for one-sided shield earthing (Schmidt [30]).

Impact of Semiconducting Layers
The influence of semiconducting layers is formulated in only very few publications.Ametani [35] determined the impact of semiconducting layers by modifying Schelkunoff's method [5].Semiconducting layers can also be considered with the subconductor method.A one-skin segmentation for the inner semiconducting layer is shown in Figure 12.
inner semiconducting layer core Segmentation can be realized also with more than one skin.The specific resistances of semiconducting layers can be sl = 0.01 . . . 10 Ω m to Ametani [35].Nevertheless, cable manufacturers specify particular values between sl = 10 and 100 Ω m.The calculated resistance per unit length R ac(C) for relevant specific semiconducting layer resistances sl is shown in Figure 13.

Main Insulation Capacitance
The single core cable capacitance shown in Figure 14 can be calculated by Equation (19) for the cross-linked polyethylene (XLPE) main insulation C CS : and for the sheath capacitance C sh (20).The capacitances depend on the cylinder geometry, the relative permittivity ε r of the insulation and the polarization.The frequency-dependent effects of polarization are investigated in several publications (Liu [36], Hadid and Schmidt [37]).The frequency range of polarization types are given in Table 1 (see also Liu [36]).The table shows that electronic, ionic, and orientational polarization will not contribute significantly to the capacitances of single core cables.A typical interfacial polarization occurs between the XLPE insulation and the semiconducting layers.This effect is of relevance in a frequency range f ≤ 1 Hz.The interfacial polarization can be ignored in frequency ranges f ≥ 10 Hz.Therefore, the frequency dependency of the capacitances can be generally neglected in an algorithm for frequency-dependent parameters.A typical interfacial polarization occurs between the XLPE insulation and the semiconducting layers.This effect is of relevance in a frequency range f ≤ 1 Hz.The interfacial polarization can be ignored in frequency ranges f ≥ 10 Hz.Therefore, the frequency dependency of the capacitances can be generally neglected in an algorithm for frequency-dependent parameters.The semiconducting layers have an influence on the capacitance C CS , but for transient calculations it can be neglected.Hadid and Schmidt [37] and Wagenaars [38] confirm this approach.

Earth Installation
Ametani [1], Wedepohl [2], Gustavsen [39], and Marti L. [13] neglected the earth resistance R E and the earth capacitance C E .In Schmidt [30], this approach is confirmed for specific earth resistances E < 500 Ω m and frequencies f < 100 kHz.The sheath capacitance per unit length C sh can be calculated with Equation (20).

Air Installation
For correct calculation of shield-shield capacitance C SS , the sheath capacitance C sh and the air capacitance C E have to be considered.The sheath capacitance C sh is significantly larger than the air capacitance C E .Therefore, the sheath capacitance C sh can be neglected (Schmidt [30]).This is also valid in case of small phase distances.The sheath capacitance C sh can be calculated with the charge simulation methods of Steinbigler [40], and Probst [41].

Approach
The wave propagation in a multi-conductor cable or transmission line can be represented in the frequency domain with second-order differential equations from Marti J. [42]: The longitudinal impedances in the phase domain Z are calculated using the partial subconductor method as in Equation (7).The shunt admittance in the phase domain Y L represents the capacitance equations given in Equations ( 19) and (20).The index "L" stands for lateral.In order to solve the differential Equations ( 21) and ( 22), a modal transformation technique has been proposed in Marti J. [42].As a result, the coupling between the cable conductors disappears and the cable system can be considered as decoupled single conductors.In the modal domain, the Equations ( 21) and ( 22) can be rewritten as: where λ is a diagonal matrix which represents the eigenvalues matrix of Z Y L and Y L Z .For the transformation into modal domain the modal transformation matrices T i and T v are used.The relation between the matrices is given by Marti J. [13]: Therefore, it is sufficient to calculate only one of them.In this paper, the results for T i will be calculated as in Chrysochos et al. [43].
In order to make the calculation in the modal domain, the quantities Z , Y L , V and I have to be transformed into this domain.From Equations ( 23) and ( 24), it can be derived: With simple matrix perturbation, the above equations are then rewritten as: By rearranging the terms and by using the Equation ( 25) in ( 28) and ( 29), the longitudinal impedance and the shunt admittance in the modal domain can be defined as: Both matrices in Equations ( 30) and ( 31) are diagonal.The frequency dependence of the elements of Y L,m cannot be neglected (as assumed for the elements of Y L ).This is related to the frequency dependent elements of T i .The elements of Z m are also strongly frequency-dependent.In order to take the frequency dependency of Z m and Y L,m into account, their diagonal elements need to be approximated by mathematical functions using the vector fitting algorithm (VF) of Gustavsen and Semlyen [15], and Gustavsen [16].For example, the first diagonal element in Z m can be approximated using VF as: The accuracy of Z m fit depends on the number of the real terms N and the number of the complex terms M. Z m fit can be reproduced with an electrical network.In Hoshmeh et al. [21], an RL network has been used.Using such a network allows us to employ only the real poles and residues in the approximation (the complex terms are discarded).Although the fitting is achieved with a low order of approximation, the accuracy is affected.In this paper, an RLC network shown in Figure 15 is introduced with regards to all the poles and residues in the approximation including the complex terms.
Representation of the frequency dependence of the cable longitudinal impedance.
As can be seen, the network has been divided into two sections.The elements of section I can be evaluated in a similar way as in Hoshmeh et al. [21].The elements of section II are obtained as follows: The first diagonal element of Y L,m is similarly approximated.Furthermore, R 0 and L 0 are assumed to be zeros in Y L,m fit .After calculating Z m fit and Y L,m fit , they will be combined to form a developed PI section.Figure 16 shows a representation for the developed PI section after regarding the cable length, where v m,s and v m,r are the voltages in the modal domain at the sending and receiving ends, respectively.One single conductor in the developed cable model can then be represented by cascading a number of PI sections.For the remaining diagonal elements of Z m and Y L,m , the same procedure is used to form the other conductors in the developed cable model.From the Equations ( 23)- (25), the voltages and the currents are transformed into the modal domain as follows: Equations ( 34) and (35) are still in the frequency domain.However, since the simulations in the power system are necessary in the time domain, these equations have to be transformed into this domain.For the voltages in the Equation (34), the convolution integral is applied as: Evaluation this integral point by point is time consuming.By approximating the elements of T tr i using rational functions, an efficient recursive convolution method is used from Marti J. [12] and the above integral can be rewritten as a sum of exponentials: where b and l are the approximations parameters of the rational functions.The wave travel time τ for the impulse from one end of the cable to the other one is related to the imaginary term of the propagation function γ: The voltage vector v m (t) is calculated recursively from v m (t − ∆ t) as: In this relation, a 1 , a 2 and a 3 are constants that mainly depend on b, l and ∆ t Marti J. [42].The same procedure is used to transform the currents in the Equation (35) into the time domain.In order to execute the time domain simulations, the differential equations for single conductors in the developed cable model have to be formulated.As shown in Hoshmeh et al. [21], the differential equations can be represented using state space techniques.A numerical backward Euler rule of integration can then be applied to find the solution for the related state space equations of Meyer and Dommel [44].To solve the cable equations with the rest of the network, which is always defined in phase quantities, the calculated voltages and currents in the modal domain must be transformed back to phase quantities.
Figure 17 shows an overview for the realization of the 3PPI model.Starting from defining Z and Y L , the cable model part is executed only once to provide the PI sections for the calculation part.In this part, the computation is evaluated every ∆ t.For the first one, the vector of the voltages v is calculated using the network model and the vector of the currents i.By applying the recursive convolution between v and the approximation of t * i ( t * i ), the vector of the voltages v m is calculated.Using the PI sections and v m , the state space equations are formulated and solved by means of numerical backward Euler rule of integration to yield the current vector i m .The recursive convolution between i m and the approximation of t i ( t i ) provides the vector of currents i for the next ∆ t.

Validation of the Full Frequency-Dependent Cable Model
For the validation of the 3PPI model, a DC-Test on a three-phase cable with a length of ≈ 2.5 km in trefoil formation has been used.The test is performed by charging the cable conductors up to 1 kV.After that, the charged conductors are connected to earth through the circuit breaker (s2) and the currents i core(A) , i core(B) , i shield(A) and i shield(B) at the sending end are measured.The test configuration is shown in Figure 18.The cable parameters are given in Table 2.The capacitances are calculated with Equations ( 19) and (20).The calculated frequency-dependent self resistance R self and self inductance L self are shown in Figure 19.The simulations with the 3PPI model have been executed with 75 PI sections.The order of approximation for Z m fit and for Y L,m fit is 8.The simulation time step is set to ∆ t = 10 −7 s.In this simulation, two facts should be replicated correctly: the wave travel time τ, which is related to the imaginary term of γ and the attenuation, which is related to the real term of γ.For the sake of clearness, Figure 20 shows only a comparison between the outcomes of 3PPI Model with the measured currents for the core and shield A. The comparison shows that the model results match very well with the measured data.This means that the wave travel time and attenuation are replicated correctly in the developed cable model which proves its validity.For example, in our case the simulated wave travel time is 14.8 µs and the measured value is 14.5 µs.The small deviation between the travel times (0.3 µs) can also be minimized by increasing the number of PI sections.However, this will increase the calculation time.
In addition, it should be noted that a certain variation between simulation and measurement is always to be expected in general.This is mostly related to the lack of accurate information about some cable data, for example the thickness of inner and outer insulation screens and also their conductivity and permittivity.The inaccuracy in the cable data will affect the accuracy of the calculated longitudinal impedance Z and shunt admittance Y L .This will be reflected in the simulation results of the cable model.
The required number of PI sections in the model depends on the expected frequencies in the transients, the cable length and the accepted calculation error.With the proposed model, calculation times of few seconds have been achieved.

Conclusions
In this paper, a full frequency-dependent cable model for calculations of transients was introduced.The frequency-dependent cable parameters are determined with an improved partial subconductor algorithm.The algorithm is able to handle a large number of conductors.Complete cable systems can be implemented including components like shields, armor, earth and other grounding conductors in parallel through the earth.If necessary, semiconducting layers can be realized in form of additional conductors.The skin and the proximity effect is considered in all conductors.
The improved partial subconductor method merely allows the calculation of the inductive coupling between conductor loops.The capacitive coupling is realized by other analytical methods.The frequency dependency of a XLPE insulation capacitance can be neglected for calculations in frequency ranges f < 5 MHz.
The wave propagation cable model based on lumped parameters.A modal transformation technique was employed in the cable model to transform the cable system into decoupled single conductors.Every single conductor in the modal domain is represented through cascaded developed PI sections, where the frequency dependence of cable parameters is implemented through an RLC network in every PI section.
For an efficient computation time, the simulations in the time domain are executed using a recursive convolution technique.The number of PI sections in the full frequency-dependent cable model depends on the simulated cable length, the expected highest frequency in the investigated transient, and the accepted calculation error.
A comparison of the outcomes of proposed 3PPI model with measured data has shown a good match, which validates the approach used in developing this model.
A next important step is to compare the results of the developed full Frequency-dependent PI section cable model with other existing cable models, which is the goal of our next publication.

Figure 1 .
Figure 1.Equivalent configuration for the determination of the mutual inductance M between linear conductor loops.

Figure 2 .
Figure 2. Principal configuration for the determination of mean geometric distance g.

Figure 3 .
Figure 3. Equivalent circuit for determination of mutual inductance M mn with common return (hull cylinder).

Figure 5 2 Figure 5 .
Figure5illustrates the current density J C /J max(C) at different frequencies for a solid, cylindrical copper conductor with a cross-section area A C of 1000 mm 2 .

Figure 8 .
Figure 8. Example for a segmentation of a conductor into three skins at frequency f = 50 Hz.

Figure 9 .
Figure 9. Current density J E /J max(E) in the earth at different frequencies f .

Figure 10 .
Figure 10.Segmentation of the earth into ten skins.

Figure 12 .
Figure 12.One-skin segmentation of the inner semiconducting layer.

Figure 13 .
Figure 13.Resistance per unit length R ac(C) considering the semiconducting layer.

Figure 14 .
Figure 14.Configuration of a single core cable.

Figure 17 .
Figure 17.Principle of the three-phase cable model.VF: vector fitting algorithm.

Figure 18 .
Figure 18.Test configuration for the current measurements.

Figure 19 .Figure 20 .
Figure 19.Simulated frequency-dependent self resistance R self and self inductance L self .

Table 1 .
Frequency ranges of polarization.