Use of 3D-FEM Tools to Improve Loss Allocation in Three-Core Armored Cables

: Loss allocation through analytical expressions in three-core lead-sheathed armored cables is challenging due to the complex geometry of this type of cable, commonly employed in submarine energy transmission systems (involving twisted conductors, sheaths and armor). Most of the expressions of the IEC standard 60287 have not been properly adapted for three-core armored cables, leading to inaccurate values for the different losses, so important efforts are currently devoted to improving them. In this work, an improved ultra-shortened 3D ﬁnite element model (FEM) is employed for developing an in-depth analysis of the electromagnetic interactions that take place in 6 real cables, being especially focused on those aspects that are not considered in the IEC standard. As a result, important conclusions are derived regarding the losses in conductors and sheaths, which introduce different corrections for improving the accuracy of the IEC expressions. The new formulation is then employed to propose a simpliﬁed experimental armor loss allocation procedure. This is virtually applied through the FEM tool to more than 700 cable conﬁgurations, showing a remarkable improvement in the loss allocation over the IEC standard and previous experimental procedures.


Introduction
Loss allocation in power cables is a requirement for the thermal rating assessment that is performed at any stage of an underground/submarine high voltage transmission project (planning, design or operation). Due to the high number of offshore wind farms developed in recent years [1,2], the interest in this topic has grown significantly. In this sense, for three-core lead-sheathed armored cables the problem of loss allocation is even more challenging due to their three-heat source design (conductors, sheaths and armor, Figure 1). This was pointed out 10 years ago when [3,4] revealed that the IEC standard [5] greatly overestimates the power losses in this type of cable, hence obtaining oversized cables. The main reason is that [5] is partially based on studies developed in the 1930 s, when cables had a different design. Moreover, some erroneous simplifying assumptions were included to obtain feasible analytical expressions (e.g., assuming the armor as a solid cylinder), leading to inaccurate values for the armor losses [3,4].
Regarding the experimental methods [26,28], even though they use real measurements of the global losses, they resort to [5] to allocate them. However, the oversimplified assumptions involved in the formulae lead to inaccurate values of the different loss components [16]. Hence, due to the difficulties of effectively separating the losses by experimental methods, 3D-FEM modelling is, without any doubt, the most accurate loss allocation method. However, due to the challenges in the FEM modelling of this type of cable (quite In this sense, recent studies have suggested some modificati nor changes (e.g., increasing the additional length that every twis compared to the overall axial length of the cable [29]), adding correction factors to include the effect of the armor in the condu cited before, developing new formulae for determining the s [6,7,31]. However, the improvements obtained from most of properly reported and assessed in the literature. Therefore, since the use of analytical expressions as those f vanced tools, such as 3D-FEM, are mainly devoted to help in the d more accurate analytical expressions. To shed some light to the p in TCACs, this work exploits the power of ultra-shortened 3D-F those aspects where the formulae of [5] can be improved. This wa marizes the traditional approach from [5], highlighting its main an ultra-shortened but accurate 3D-FEM model is proposed, incl veloped to get an in-depth knowledge of the physical phenome are applied in Section 4 to analyze separately the conductor an cables. As a result, different corrections are proposed to enhance mulation of these losses in [5]. Section 5 proposes a simplified that leverages these corrections to obtain a new assessment of th ally, the conclusions are included in Section 6, also providing so research lines aimed at obtaining eventually universal and accur that could be incorporated in a future IEC standard version. In this sense, recent studies have suggested some modifications, as, for example, minor changes (e.g., increasing the additional length that every twisted phase conductor has compared to the overall axial length of the cable [29]), adding some empirically-based correction factors to include the effect of the armor in the conductor losses [29,30] or, as cited before, developing new formulae for determining the sheath and armor losses [6,7,31]. However, the improvements obtained from most of these proposals are not properly reported and assessed in the literature.

Traditional Approach to Losses Estimation in TCACs
Therefore, since the use of analytical expressions as those from [5] is preferred, advanced tools, such as 3D-FEM, are mainly devoted to help in the development of new and more accurate analytical expressions. To shed some light to the problem of loss allocation in TCACs, this work exploits the power of ultra-shortened 3D-FEM models to highlight those aspects where the formulae of [5] can be improved. This way, Section 2 briefly summarizes the traditional approach from [5], highlighting its main weaknesses. In Section 3 an ultra-shortened but accurate 3D-FEM model is proposed, including also new tools developed to get an in-depth knowledge of the physical phenomena involved. These tools are applied in Section 4 to analyze separately the conductor and sheath losses in 6 real cables. As a result, different corrections are proposed to enhance the accuracy of the formulation of these losses in [5]. Section 5 proposes a simplified experimental procedure that leverages these corrections to obtain a new assessment of the armor losses. Eventually, the conclusions are included in Section 6, also providing some remarks about open research lines aimed at obtaining eventually universal and accurate analytic expressions that could be incorporated in a future IEC standard version.

Traditional Approach to Losses Estimation in TCACs
In [5] the conductor losses (P c ) are derived as follows: R c (θ) = R DC c (θ) 1 + y s + y p , being R c (θ) and R DC c (θ) the AC and DC conductor resistances at temperature θ, respectively, I c the conductor current, y s the skin effect factor, y p the proximity effect factor, ρ c the Energies 2021, 14, 2434 3 of 23 conductor resistivity, S c the cross section of the conductor, α c the temperature constant and θ 0 the reference temperature. On the other hand, the sheath (P s ) and armor (P a ) losses are obtained by means of factors λ 1 (sheath loss factor) and λ 2 (armor loss factor), defined as Regarding sheath losses, [5] assumes that they are a consequence of eddy and circulating currents [32], so λ 1 can be defined as λ 1 = λ 1 + λ 1 , being λ 1 and λ 1 the loss factors of circulating and eddy-current losses, respectively. However, [5] only provides expressions for λ 1 and λ 1 in the case of three single-core cable arrangements, depending on whether the sheaths are in single-point bonding (SP) or solid bonding (SB) configuration. Hence, the same expressions have been applied for the case of TCACs, although only the corresponding to λ 1 has been modified to include the effect of the armor (by introducing an empirical factor 1.5): with being R s the sheath resistance in Ω/m, X s the sheath reactance in Ω/m, s the spacing between conductor axes, d the mean diameter of the sheath, t s the sheath thickness, µ 0 the vacuum permeability, ω the angular frequency and g s , λ 0 , ∆ 1 , ∆ 2 and β 1 parameters defined in [5]. It should be remarked that in the case of sheaths in SB, [5] also assumes that λ 1 λ 1 , so it is neglected. Nonetheless, there are some aspects of the design of this type of cable that are not adequately considered in this formulation. In this sense, previous works [3,4,10,14,25] showed that phase and armor twisting, as well as the magnetic properties of the armor wires, increase the positive-sequence impedance of these cables. Therefore, R c , R s and X s are affected by these aspects too. In this sense, [6][7][8] have suggested some modifications for improving the estimation of R c . On one hand, considering the additional length of the conductors due to their twisting, leading to an increase in the conductor resistance [29], is suggested. This is done by multiplying R DC c in Equation (3) by the lay factor LF, defined as where L c is the lay length of the conductors and R t denotes the distance cable centerconductor center. On the other hand, to include the influence of the armor on R c it has been suggested to multiply both skin and proximity factors in Equation (2) by 1.5 [30]. However, the source of this proposal is still to be published (technical brochure of CIGRE W.G. B1.56). Although both improvements are directly applied to R c , they also influence the sheath loss factor obtained by Equation (6), where R c plays an important role. In any case, the improvements derived from these proposals are not properly reported in the literature, so they must still be evaluated.
Another point of controversy is found in the bibliography regarding the source of the induced losses in the sheaths. Some studies [6,7,33,34] suggest that eddy-current losses may be proportional to circulating-current losses for large cables, so they should not be ignored in the computation of the SB-case sheath losses. However, there is no expression in [5] especially designed for evaluating eddy losses in this type of cable (Equation (7) is in fact only valid for SP single-core cables). Hence, the effects of the core and armor twisting, as well as their magnetic properties, are not appropriately considered.
It is to be remarked that all these aspects cannot be accurately evaluated through laboratory measurements due to the impossibility of individually measuring the three loss components [3,4,20,21,[23][24][25][26][27][28], there existing still a degree of uncertainty that assigns erroneously the sheath eddy losses to other metallic elements. Nonetheless, regarding the armor losses, the experimental results are conclusive enough [3,4], verifying that [5] highly overestimates λ 2 . In this sense, many efforts are now devoted to improving the accuracy of the expressions included in [5] (e.g., CIGRE W.G. B1.64 [31]).
As can be seen, although some improvements have been suggested in the formulation of the conductor losses, its accuracy still has to be properly evaluated. Moreover, the sheath and armor losses still require further improvements to a greater or lesser extent. All this is deeply analyzed in the following sections through the ultra-shortened 3D-FEM model presented next.

Ultra-Shortened 3D-FEM Model
In [14] it was proposed a new methodology to reduce the computational burden required for solving the 3D electromagnetic field in TCACs. For this task, rotated periodicity was introduced, based on the repetition of a certain electromagnetic field pattern shown by this type of cable in a length equal to the crossing pitch (CP), a parameter defined in [6,25] for cables with armor and phases twisted in opposite direction (contralay) as where L a and L c are, respectively, the armor wire and phase conductor lay length. Thanks to this approach it is possible to analyze complex three-core cables with numerous wires in the armor by solving the electromagnetic field in just a short portion of the cable length equal to CP. Nonetheless, the simulation time may still be high for the largest cables, and it may be necessary to make use of a coarser mesh to reduce the number of degrees of freedom (DoF), thus solving faster, but providing less accurate results. To overcome this issue, new advances are presented in next subsections. The software Comsol Multiphysics [35] has been employed for this task.

Ultra-Shortened 3D Electromagnetic FEM Model
An in-depth analysis of the results obtained in numerous TCACs showed that the full electromagnetic behavior of the cable can be found in just a small slice of the cable, so a much shorter geometry for the FEM model can be employed [36,37]. This is possible since all armor wires are twisted around the three sheath-conductor cores, going through all the possible positions relative to every core, hence each one collects the full electromagnetic pattern. Therefore, it is possible to find a shorter section of the cable (L) where the full electromagnetic pattern is distributed between all the armor wires included in that section. This length is where n is any divisor of the number of armor wires (N). An example is shown in Figure 2, where a CP-length section of a 115 kV, 240 mm 2 three-core cable is represented (for a clearer figure the loss color plot is only depicted for one armor wire). Considering that N = 98 in this cable, the possible values for n are 1, 7, 14, 49 and 98. In the case of Figure 2 (n = 7) it can be easily observed how the loss pattern shown by the selected armor wire is split into CP/7-length sections, each one having the same loss distribution as other armor wires located every N/7th position in Section 1 of the figure. This way, the full pattern of that where n is any divisor of the number of armor wires (N). An example is shown in Figure  2, where a CP-length section of a 115 kV, 240 mm 2 three-core cable is represented (for a clearer figure the loss color plot is only depicted for one armor wire). Considering that N = 98 in this cable, the possible values for n are 1, 7, 14, 49 and 98. In the case of Figure 2 (n = 7) it can be easily observed how the loss pattern shown by the selected armor wire is split into CP/7-length sections, each one having the same loss distribution as other armor wires located every N/7th position in Section 1 of the figure. This way, the full pattern of that armor wire can be fully reproduced by adequately aligning those CP/7-length armor wires of Section 1. All this can be done for all n possible values, as shown in Figure 3, where the armor wires required to reproduce the full loss pattern of the single CP-length armor wire of Figure 2 are depicted. Moreover, the length of the resulting ultra-shortened geometry (depending on n) is also indicated. To reproduce the full behavior of the armor, the same methodology must be applied to the remaining armor wires, except for n = N, since the resulting shortened geometry already contains the full information. Consequently, only a CP/n-length section of the cable is required to have all the phenomena involved in the three-core cable.
It is to be noticed that this pattern is also observed in vector magnitudes, such as the magnetic flux density inside the armor wires, so this approach is fully valid and contains all the physical interactions involved in the power cable. Taking all this into account, boundary conditions for applying rotated periodicity must be adequately set for a perfect matching between the source and destination boundaries (defined in Figure 4). This is achieved by assigning coordinate systems , , and , , to both boundaries, being the destination boundary system obtained by rotating the source coordinate system a certain angle that is a function of n, as All this can be done for all n possible values, as shown in Figure 3, where the armor wires required to reproduce the full loss pattern of the single CP-length armor wire of Figure 2 are depicted. Moreover, the length of the resulting ultra-shortened geometry (depending on n) is also indicated. To reproduce the full behavior of the armor, the same methodology must be applied to the remaining armor wires, except for n = N, since the resulting shortened geometry already contains the full information. Consequently, only a CP/n-length section of the cable is required to have all the phenomena involved in the three-core cable.
that armor wire can be fully reproduced by adequately aligning those CP/7-len wires of Section 1. All this can be done for all n possible values, as shown in Figure 3, where wires required to reproduce the full loss pattern of the single CP-length arm Figure 2 are depicted. Moreover, the length of the resulting ultra-shortened geo pending on n) is also indicated. To reproduce the full behavior of the armor methodology must be applied to the remaining armor wires, except for n = N resulting shortened geometry already contains the full information. Conseque CP/n-length section of the cable is required to have all the phenomena invo three-core cable.
It is to be noticed that this pattern is also observed in vector magnitudes, magnetic flux density inside the armor wires, so this approach is fully valid an all the physical interactions involved in the power cable. Taking all this into account, boundary conditions for applying rotated must be adequately set for a perfect matching between the source and destinat aries (defined in Figure 4). This is achieved by assigning coordinate systems and , , to both boundaries, being the destination boundary system o rotating the source coordinate system a certain angle that is a function of n, It is to be noticed that this pattern is also observed in vector magnitudes, such as the magnetic flux density inside the armor wires, so this approach is fully valid and contains all the physical interactions involved in the power cable.
Taking all this into account, boundary conditions for applying rotated periodicity must be adequately set for a perfect matching between the source and destination boundaries (defined in Figure 4). This is achieved by assigning coordinate systems where the minus sign is for phases twisted counterclockwise and the positive sign for phases twisted clockwise. where the minus sign is for phases twisted counterclockwise an phases twisted clockwise. being the magnetic vector potential, the conductivity, the and the external current density. Rotated periodicity is applied where ( , , ) and ( , , ) are the coordinates of matchin tion and source coordinate systems, respectively ( Figure 4). An example of the improvement achieved by the proposed next for the same cable of 115 kV. Its main dimensions and materi in Table 1, where dc is the conductor diameter, da is the armor wire mean diameter of the armor.  This way, the electromagnetic problem is analyzed by solving being → A the magnetic vector potential, σ the conductivity, µ the magnetic permeability and → J e the external current density. Rotated periodicity is applied by means of where x d , y d , z d and (x s , y s , z s ) are the coordinates of matching points in the destination and source coordinate systems, respectively ( Figure 4). An example of the improvement achieved by the proposed approach is presented next for the same cable of 115 kV. Its main dimensions and material properties are shown in Table 1, where d c is the conductor diameter, d a is the armor wire diameter and D a is the mean diameter of the armor. For this cable, Table 2 shows the most relevant simulation data and the main results obtained for different values of n, where R and X represent the cable positive-sequence resistance and the inductive reactance, I s is the sheath current, and P c , P s and P a are the conductor, sheath and armor losses, respectively. Additionally, relevant data regarding DoF, mesh size and quality, and the RAM memory employed for the simulations are included. Simulations have been performed in a laptop with an i7 processor and 64 GB of RAM memory. It must be remarked that for n = 1 a coarser mesh was necessary for keeping the RAM memory requirements below the available 64 GB. Instead, finer meshes were employed for the rest of the cases, increasing the mesh density with n. This improves the average and the minimum element quality ( Table 2, where 1 is the highest quality), and hence the overall quality of the mesh.  [5] 0.246 λ 2 (FEM) 0.0415 0.0413 0.0413 0.0414 0.0415 λ 2 [5] 0.221 As can be seen, differences in the results derived for each value of n are virtually negligible, and the reduction achieved in computation time for the case of n = N is to be remarked, with just 30 s for solving such a complex model and with no loss in accuracy. Consequently, a model length equal to CP/N is the best choice. Further, it is of interest to highlight how there is a good agreement for this cable between the λ 1 values from FEM and [5]. However, important differences are observed in λ 2 , as previously reported in [3,4].

Loss Allocation
For a better understanding of all the phenomena involved in loss generation, as well as analyzing its dependence with technical and material parameters, some additional improvements have been implemented in the 3D-FEM model.

Comparing Results
The FEM software employed in this study allows for making simple numerical operations with the solutions, of any spatial magnitude or field, obtained for the same geometry but considering two different sets of boundary conditions (C1 and C2). This feature is used here to define a ratio γ for evaluating, on each point of the cable geometry (x, y, z), Energies 2021, 14, 2434 8 of 23 the variation of a certain magnitude ξ when boundary conditions are changed from C1 to C2, as When plotted over the cable geometry this ratio allows for visualizing the regions of the cable where ξ increases (γ > 1) or decreases (γ < 1) in relation to the reference situation (C1). An example is shown in Figure 5 for a 275 kV, 2000 mm 2 three-core unarmored cable (Table 3), where it is represented the variation that takes place in the conductor losses (ξ = P c ) when switching from sheaths in SP (C1) to SB (C2).

1, 14, x FOR PEER REVIEW
C2, the average value of can be taken as the factor in affected by changes in the bonding configuration. In thi itatively and quantitatively how variations in the bound ductor resistance.   From Figure 5 it can be concluded that the current density in the conductors is affected by the bonding scheme, leading to an increase in losses of more than 20% in half of the conductor when changing to SB, while it decreases in less than 20% in the remaining regions. The extension of each region results in a global increase or decrease in ξ (P c in this example). In this case, since the current through the conductors is the same in C1 and C2, the average value of γ can be taken as the factor in which the conductor resistance is affected by changes in the bonding configuration. In this way it is easier to evaluate qualitatively and quantitatively how variations in the boundary conditions influence the conductor resistance.

Separation of Circulating and Eddy Losses
As mentioned before, there is a controversial point regarding loss allocation in the sheaths, since [5] assumes that eddy losses may be neglected when sheaths are in SB, although some studies suggest the opposite for large cables [6,7,33,34]. Therefore, to clarify this point it is required to find a way to evaluate these terms separately, especially in TCACs. In this case, the induced currents mainly follow a helical path due to the conductor twisting (Figure 6a for a 275 kV, 2000 mm 2 three-core unarmored cable). Therefore, it is expected that eddy losses differ from those derived from Equation (7), which is intended for three parallel and straight single-core cables. Moreover, the presence of a magnetic armor is not considered in this expression.  Having this in mind, it is here proposed a way to evaluate sheath circulating eddy losses separately by means of 3D-FEM simulations. To this aim, it is assumed the external flux density causing eddy-current losses is not substantially affected by bonding scheme, i.e., the same eddy-current distribution is expected in the sheaths. Th the 3D-FEM model of the cable is first solved for sheaths in SP and SB. Once done, FEM software allows subtracting the eddy-current density distribution obtained in SP-cable ( ) from that derived for the SB situation ( ), as shown in Figure 6. This res in a current density distribution ( ) that absolutely follows the helical path of sheaths ( Figure 6c). Furthermore, the resulting net current is also equal to that measu in the sheaths for the SB-cable, thus verifying our assumption. This current density re sents the cause of the circulating losses (P's) which can be evaluated integrating thro the three sheath volumes by where is the sheath resistivity. Having this in mind, it is here proposed a way to evaluate sheath circulating and eddy losses separately by means of 3D-FEM simulations. To this aim, it is assumed that the external flux density causing eddy-current losses is not substantially affected by the bonding scheme, i.e., the same eddy-current distribution is expected in the sheaths. Thus, the 3D-FEM model of the cable is first solved for sheaths in SP and SB. Once done, the FEM software allows subtracting the eddy-current density distribution obtained in the SP-cable ( → J SP ) from that derived for the SB situation ( → J SB ), as shown in Figure 6. This results in a current density distribution ( → J long ) that absolutely follows the helical path of the sheaths (Figure 6c). Furthermore, the resulting net current is also equal to that measured in the sheaths for the SB-cable, thus verifying our assumption. This current density represents the cause of the circulating losses (P s ) which can be evaluated integrating through the three sheath volumes by where ρ s is the sheath resistivity. On the other hand, since the net current flowing through the sheaths can be computed (I s ), it is also possible to obtain the AC resistance of the sheaths from P s as Finally, since the global sheath losses in SB can be also derived from → J SB by Equation (17), it is now possible to obtain eddy losses (P s ") as An in-depth analysis of the results derived from this method will be discussed in the following sections.

Analysis of Conductor and Sheath Losses
Inaccurate results are derived from [5] since R s , X s and R c are obtained by ignoring important aspects related to TCACs, such as the core and armor twisting, the bonding scheme, as well as the magnetic properties of the armor wires. The influence of these aspects on P c and P s is analyzed in this section, being these two losses individually evaluated by means of the improved 3D-FEM model ( Figure 7). Moreover, the proposals presented in previous works [6,10,29,31] are also evaluated. As a result, different approaches for improving both the conductor and sheath losses are proposed. To this aim, 6 test cables are employed (Table 3), although most of the analyses throughout the paper are illustrated using the largest cable (Cable 6) for clarity. The electrical resistivity of the conductive parts is taken from [5] for an ambient temperature of 20 °C. A complex relative permeability in the form of = − ′′ (ranging from 1 − 0 (no armor/stainless steel armor) to 600 − 350) is employed for the armor wires in order to consider hysteresis losses.

Conductor Losses
The effects of core twisting, sheath bonding and the armor on the conductor losses (not considered in [5]) are analyzed next.

Core Twisting
The influence of the core twisting is shown in Figure 8 for Cable 6 when neither armor nor sheaths are present, where the evolution of Pc (Figure 8a) and Rc (Figure 8b) with Lc are represented (Cable 6 Lc6 is marked with a vertical line). It is easily seen how the values of Rc and Pc, derived from FEM simulations, increase with the twisting of the conductors. This leads to important differences with those derived by the IEC standard (labelled as "Pc IEC"), especially when phases are more twisted (shorter lay length), where it is more evident that [5] underestimates their values. As a result, different approaches for improving both the conductor and sheath losses are proposed. To this aim, 6 test cables are employed (Table 3), although most of the analyses throughout the paper are illustrated using the largest cable (Cable 6) for clarity. The electrical resistivity of the conductive parts is taken from [5] for an ambient temperature of 20 • C. A complex relative permeability in the form of µ r = µ − jµ (ranging from 1 − j0 (no armor/stainless steel armor) to 600 − j350) is employed for the armor wires in order to consider hysteresis losses.

Conductor Losses
The effects of core twisting, sheath bonding and the armor on the conductor losses (not considered in [5]) are analyzed next.

Core Twisting
The influence of the core twisting is shown in Figure 8 for Cable 6 when neither armor nor sheaths are present, where the evolution of P c (Figure 8a) and R c (Figure 8b) with L c are represented (Cable 6 L c6 is marked with a vertical line). It is easily seen how the values of R c and P c , derived from FEM simulations, increase with the twisting of the conductors. This leads to important differences with those derived by the IEC standard (labelled as "P c IEC"), especially when phases are more twisted (shorter lay length), where it is more evident that [5] underestimates their values.
The influence of the core twisting is shown in Figure 8 for Cable 6 when neither armor nor sheaths are present, where the evolution of Pc (Figure 8a) and Rc (Figure 8b) with Lc are represented (Cable 6 Lc6 is marked with a vertical line). It is easily seen how the values of Rc and Pc, derived from FEM simulations, increase with the twisting of the conductors. This leads to important differences with those derived by the IEC standard (labelled as "Pc IEC"), especially when phases are more twisted (shorter lay length), where it is more evident that [5] underestimates their values.
(a) (b) As recommended in [29], the additional length of the conductor due to the twisting of the phases should be considered in the analytical computation of the DC resistance. This is done by multiplying Equation (3) by LF (previously defined in Equation (9)). This As recommended in [29], the additional length of the conductor due to the twisting of the phases should be considered in the analytical computation of the DC resistance. This is done by multiplying Equation (3) by LF (previously defined in Equation (9)). This correction in [5] leads to better results in Figure 8 (denoted by "P c IEC LF" and "R c IEC LF", respectively), since lower differences with FEM results are observed.
Nonetheless, some differences are still observed. This might be explained as an increase in the proximity effect due to the conductor twisting, which leads to a slightly higher AC resistance. In this sense, since y p is expressed as a function of d c /s, a first hypothesis is that the concepts of conductor-conductor distance and conductor diameter must be re-defined in case of twisting. To better visualize this, Figure 9 shows the resulting cross section obtained by a cut plane orthogonal to one of the phases (indicated by a green arrow). As can be seen, while the cross section of one conductor remains circular in that plane, the other two are in fact ellipses, where the major axis is higher that the conductor diameter in LF 3 , and the minor axis remains equal to d c . Moreover, in that plane the distance between the axes of circular and elliptical conductors is s· cos α/2, being 4, x FOR PEER REVIEW 11 of correction in [5] leads to better results in Figure 8 (denoted by "Pc IEC LF" and "Rc IE LF", respectively), since lower differences with FEM results are observed. Nonetheless, some differences are still observed. This might be explained as an i crease in the proximity effect due to the conductor twisting, which leads to a slight higher AC resistance. In this sense, since yp is expressed as a function of dc/s, a first h pothesis is that the concepts of conductor-conductor distance and conductor diamet must be re-defined in case of twisting. To better visualize this, Figure 9 shows the resulti cross section obtained by a cut plane orthogonal to one of the phases (indicated by a gre arrow). As can be seen, while the cross section of one conductor remains circular in th plane, the other two are in fact ellipses, where the major axis is higher that the conduct diameter in LF 3 , and the minor axis remains equal to dc. Moreover, in that plane the d tance between the axes of circular and elliptical conductors is • cos /2, being = atan − 1 .
(2 All this influences the magnetic flux concatenated by the elliptical conductors (Figure 9a), reinforcing the proximity effect. Consequently, we propose to multiply eve term dc/s in the expression of yp (given in [5]) by a coefficient ft as follows After applying this approach, differences with FEM results decrease even more, observed in Figure 8 (denoted by "Pc corr" and "Rc corr", respectively). All this influences the magnetic flux concatenated by the elliptical conductors (Figure 9a), reinforcing the proximity effect. Consequently, we propose to multiply every term d c /s in the expression of y p (given in [5]) by a coefficient f t as follows

Effect of Sheath Bonding
where x p , G x p and H x p are defined in [5]. After applying this approach, differences with FEM results decrease even more, as observed in Figure 8 (denoted by "P c corr" and "R c corr", respectively).

Effect of Sheath Bonding
The presence of sheaths in a three-core cable (either SP or SB) also influences the conductor losses. This is shown in Figure 10a, where 3D-FEM results for P c are represented for Cable 6 ( Table 3, but unarmored) as a function of L c . Additionally, the case where the three-core cable has no sheaths is included (denoted by "P c ns"). As can be observed, the value of Pc, and hence its resistance, is higher when no sheaths are installed, and it reduces in presence of sheaths, especially when connected in SB, being more noticeable when phases are more twisted. A possible answer to this behavior may be found in Figure 10b, where the ratio is depicted for the current density J from absence ( ) to presence ( ) of shields. Thus, it is observed how, in both SB and SP, in half of the conductor cross section the current density increases over the situation with no sheaths, while it decreases in the other half. However, these variations are more intense in the SB case, reaching up to ±15%. In this sense, an average evaluation of in the conductor volume results in a factor of 0.99 for the SP case and 0.98 for the SB situation, in accordance with the results observed in Figure 10a. Anyway, these values are quite low, so its influence can be neglected (this effect will not be negligible with armor, as analyzed next).

Combined Effect of Magnetic Armor and Sheath Bonding
As previously indicated, to analyze the influence of the magnetic armor on Rc and Pc, only sheaths in SB are considered, since it is the preferred bonding configuration in offshore applications [30]. This way, the results derived by FEM simulations for Cable 6 are represented in Figure 11, where different values of µr, La and Lc are considered (the actual values Lc6 and La6 of Cable 6 are also shown). Figure 11 also includes the results obtained As can be observed, the value of P c , and hence its resistance, is higher when no sheaths are installed, and it reduces in presence of sheaths, especially when connected in SB, being more noticeable when phases are more twisted. A possible answer to this behavior may be found in Figure 10b, where the ratio γ is depicted for the current density J from absence (J ns c ) to presence (J s c ) of shields. Thus, it is observed how, in both SB and SP, in half of the conductor cross section the current density increases over the situation with no sheaths, while it decreases in the other half. However, these variations are more intense in the SB case, reaching up to ±15%. In this sense, an average evaluation of γ in the conductor volume results in a factor of 0.99 for the SP case and 0.98 for the SB situation, in accordance with the results observed in Figure 10a. Anyway, these values are quite low, so its influence can be neglected (this effect will not be negligible with armor, as analyzed next).

Combined Effect of Magnetic Armor and Sheath Bonding
As previously indicated, to analyze the influence of the magnetic armor on R c and P c , only sheaths in SB are considered, since it is the preferred bonding configuration in offshore applications [30]. This way, the results derived by FEM simulations for Cable 6 are represented in Figure 11, where different values of µ r , L a and L c are considered (the actual values L c6 and L a6 of Cable 6 are also shown). Figure 11 also includes the results obtained in the unarmored case, together with those derived by [5] and the approaches proposed in previous subsections (Equations (9) and (21)) (labelled as "P c corr" in Figure 11). accordance with the results observed in Figure 10a. Anyway, these values are quite low, so its influence can be neglected (this effect will not be negligible with armor, as analyzed next).

Combined Effect of Magnetic Armor and Sheath Bonding
As previously indicated, to analyze the influence of the magnetic armor on Rc and Pc, only sheaths in SB are considered, since it is the preferred bonding configuration in offshore applications [30]. This way, the results derived by FEM simulations for Cable 6 are represented in Figure 11, where different values of µr, La and Lc are considered (the actual values Lc6 and La6 of Cable 6 are also shown). Figure 11 also includes the results obtained in the unarmored case, together with those derived by [5] and the approaches proposed in previous subsections (Equations (9) and (21)) (labelled as "Pc corr" in Figure 11). A first conclusion arises from the results regarding the stainless-steel armor (µ r = 1 − j0), since it has no effects on the conductor losses as the results overlap with those obtained for an unarmored cable. In both cases the results are quite close to those derived by [5], especially when the corrections previously suggested are applied ("P c corr"). However, in presence of a magnetic armor a different behavior is observed, since P c increases with µ r , especially when conductors are more twisted (Figure 11a), leading to important differences with regard to the values provided by [5]. Furthermore, the lay length of the armor also influences the conductor losses, but to a lesser extent (Figure 11b).
The reasons for these results come again from the use of γ. In this case, Figure 12a,b illustrates, respectively, the changes in the conductor losses and the current density distribution when switching from absence (P na c , J na c ) to presence (P a c , J a c ) of armor. The results show how both losses and current density increase in areas close to the center of the cable, whereas it strongly decreases in the most external regions of the conductors. This can be understood as an increase in the proximity effect inside the conductors since the current density moves towards the inner region of the cable. This conclusion is in line with recent recommendations [30] regarding the application of a factor of 1.5 to both y s and y p for including the effect of the armor on R c . However, from this analysis it seems reasonable to apply the correction factor only to y p , resulting in the curve denoted by "P c cor2" in Figure 11. By doing this, the differences relative to FEM results decrease, although conductor losses are generally overestimated.
recommendations [30] regarding the application of a factor of 1.5 to both ys and yp for including the effect of the armor on Rc. However, from this analysis it seems reasonable to apply the correction factor only to yp, resulting in the curve denoted by "Pc cor2" in Figure 11. By doing this, the differences relative to FEM results decrease, although conductor losses are generally overestimated. Hence, we propose the following corrections in Equations (2) and (3) [29]: where was defined in Equation (18). The errors in Rc (relative to FEM) obtained from this approach and [5] are represented in Figure 13 for the 6 cables and different values of µr. It is easily observed how the new proposal reduces the relative differences, showing a maximum below 3% for Cable 6, while the IEC standard underestimates the conductor resistance in up to −6% in Cable 5. Hence, we propose the following corrections in Equations (2) and (3) [29]: where y p was defined in Equation (18). The errors in R c (relative to FEM) obtained from this approach and [5] are represented in Figure 13 for the 6 cables and different values of µ r . It is easily observed how the new proposal reduces the relative differences, showing a maximum below 3% for Cable 6, while the IEC standard underestimates the conductor resistance in up to −6% in Cable 5.
ergies 2021, 14, x FOR PEER REVIEW Figure 13. Error in conductor resistance, relative to FEM solutions, of the results deri formulation provided by [5] and the proposed approaches.

Sheath Losses
Although important advances have been developed recently [6], some Figure 13. Error in conductor resistance, relative to FEM solutions, of the results derived from the formulation provided by [5] and the proposed approaches.

Sheath Losses
Although important advances have been developed recently [6], some aspects such as the core twisting or the presence of the armor are not considered in the available formulae [5,32]. The lack of accurate expressions for this type of cable, as well as the influence of the mentioned aspects on the sheath losses are described next (only sheaths in SB will be considered, as commented earlier).

Influence of Twisting and Armor µ r
The influence of the cores and armor twisting in the total sheath losses are depicted in Figure 14a for Cable 6, where CP collects in one factor different sets of values for these parameters (CP 6 of Cable 6 is marked with a vertical line). Additionally, the effect of µ r is represented, as well as the values for the sheath losses derived by [5] when twisting is taken into consideration by including LF in the computation of R s : where S s is the cross section and α s is the temperature constant of the sheath. Figure 13. Error in conductor resistance, relative to FEM solutions, of the results derived from the formulation provided by [5] and the proposed approaches.

Sheath Losses
Although important advances have been developed recently [6], some aspects such as the core twisting or the presence of the armor are not considered in the available formulae [5,32]. The lack of accurate expressions for this type of cable, as well as the influence of the mentioned aspects on the sheath losses are described next (only sheaths in SB will be considered, as commented earlier).

Influence of Twisting and Armor µr
The influence of the cores and armor twisting in the total sheath losses are depicted in Figure 14a for Cable 6, where CP collects in one factor different sets of values for these parameters (CP6 of Cable 6 is marked with a vertical line). Additionally, the effect of µr is represented, as well as the values for the sheath losses derived by [5] when twisting is taken into consideration by including LF in the computation of Rs: where is the cross section and is the temperature constant of the sheath.
(a) (b) Figure 14. Evolution of (a) the total sheath losses and (b) the sheath loss factors with CP and the armor permeability in Cable 6 (SB). From Figure 14a it can be concluded that both CP and µ r influence the sheath losses, especially for lower values of CP. Although results show a good agreement between FEM and the IEC expressions for the case of the unarmored cable, it is not the case for armored ones, resulting in values generally overestimated. This is a consequence of applying an empirical factor of 1.5 in Equation (6) for computing λ 1 .
To analyze how circulating and eddy losses are affected by the armor, the loss factors λ 1 and λ" 1 are derived by means of the methodology proposed in Section 3.2.2. The results are represented in Figure 14b for Cable 6 as a function of CP and the armor permeability. It is observed how both λ 1 and λ" 1 increase with µ r and CP −1 , especially λ 1 for lower values of CP. Similar results are derived for Cables 1 to 5 (not included for clarity).
On the other hand, the results obtained for λ 1 and λ" 1 by the analytical expressions included in [5], the new expressions proposed in [6], as well as those derived by FEM simulations, are compared in Figure 15 for the 6 cables employed in this work, where only two values for µ r are considered for clarity. mainly for µr = 100 − j50 (Figure 15a), although [5,6] overestimate the total losses in the larger cables (cables 5, 6). It is also to be remarked that only the expressions of [6] include the effect of µr in the sheath losses, that increase with µr as expected (Figure 15b) (results obtained by [5] remain unaffected). However, it provides greater values for λ'1 than FEM, leading to an overestimation of the total losses, especially for the largest cable. Moreover, FEM simulations provide slightly greater values for λ"1 than [5,6], although differences in λ'1 are more significant.
From these results it is also possible to analyze the relative importance of λ'1 and λ"1 in the sheath total losses. This is shown in Figure 16, where the ratios λ1/λ'1 ( Figure 16a) and λ1/λ"1 (Figure 16b), derived by FEM simulations, are represented for all 6 cables as a function of the real part of µr(µ'). Figure 15. Circulating, eddy and total losses for the 6 cables with (a) µ r = 100 − j50 and (b) µ r = 600 − j350 derived by 3D-FEM simulations and the analytical expressions presented in [5,6].
From the analysis of Figure 15 it is concluded that [5,6] and FEM simulations provide similar results for the smaller cables (cables 1, 3) and the aluminum conductor (cable 4), mainly for µ r = 100 − j50 (Figure 15a), although [5,6] overestimate the total losses in the larger cables (cables 5, 6). It is also to be remarked that only the expressions of [6] include the effect of µ r in the sheath losses, that increase with µ r as expected (Figure 15b) (results obtained by [5] remain unaffected). However, it provides greater values for λ 1 than FEM, leading to an overestimation of the total losses, especially for the largest cable. Moreover, FEM simulations provide slightly greater values for λ" 1 than [5,6], although differences in λ 1 are more significant.
From these results it is also possible to analyze the relative importance of λ 1 and λ" 1 in the sheath total losses. This is shown in Figure 16, where the ratios λ 1 /λ 1 (Figure 16a) and λ 1 /λ" 1 (Figure 16b), derived by FEM simulations, are represented for all 6 cables as a function of the real part of µ r (µ ). It can be seen how both ratios hardly vary for µ' > 100, ranging from 1.26 to 1.32 in the case of λ1/λ'1, while λ1/λ"1 does it from 4 to 4.7. Similar results are derived by [5,6], although a wider range is observed for both parameters, having λ'1 a greater relative weight in λ1 than that observed in FEM results. In any case, these values can be taken as references for the evaluation of λ1 in other armored cables. It can be seen how both ratios hardly vary for µ > 100, ranging from 1.26 to 1.32 in the case of λ 1 /λ 1 , while λ 1 /λ" 1 does it from 4 to 4.7. Similar results are derived by [5,6], although a wider range is observed for both parameters, having λ 1 a greater relative weight in λ 1 than that observed in FEM results. In any case, these values can be taken as references for the evaluation of λ 1 in other armored cables.

Sheath Current Evaluation
The discrepancies observed previously in the sheath losses may be caused by differences in the sheath current calculated by each method. In this sense, previous references like [5,38] only provide an expression for three single-core cables in SB, where neither the core twisting nor the presence of a magnetic armor are considered. In contrast, [6] proposes new expressions for computing X s that include all these aspects. In both cases the sheath current is derived as To show how µ r , as well as the twisting of phase/armor (both included in the CP), influence I s , some results derived by 3D-FEM simulations for Cable 6 are represented in Figure 17. The evolution of I s in Figure 17 suggests an important influence of both parameters, since it increases with µ r and CP −1 . This behavior is more remarkable for lower values of CP, where FEM results deviate far from the sheath current estimated by Equation (26). It is also noticed that Equation (26) overestimates the induced current in the cases with no armor and stainless-steel armor. This is a consequence of the non-uniform current density distribution in the conductors due to skin and proximity effects, something not considered in Equation (26), where I c is assumed to be applied in the center of the conductor [10,33]. Therefore, both Rs and Xs are also affected by all these aspects. lem is even more complicated since the induced currents are non in the sheaths, as shown in Figure 18a for Cable 6 with µr = 300 − j2 employed for evaluating the changes in the loss distribution when (Figure 18b), an important increase is observed in the sheath loss armor wires, where the electromagnetic coupling between both ele induced currents that contribute to increasing both eddy and circu Therefore, both R s and X s are also affected by all these aspects. Furthermore, the problem is even more complicated since the induced currents are non-uniformly distributed in the sheaths, as shown in Figure 18a for Cable 6 with µ r = 300 − j200. In this sense, if γ is employed for evaluating the changes in the loss distribution when the armor is considered (Figure 18b), an important increase is observed in the sheath losses in areas closer to the armor wires, where the electromagnetic coupling between both elements results in higher induced currents that contribute to increasing both eddy and circulating currents. lem is even more complicated since the induced currents are non-uniformly distr in the sheaths, as shown in Figure 18a for Cable 6 with µr = 300 − j200. In this sense, employed for evaluating the changes in the loss distribution when the armor is cons (Figure 18b), an important increase is observed in the sheath losses in areas closer armor wires, where the electromagnetic coupling between both elements results in induced currents that contribute to increasing both eddy and circulating currents. However, this does not necessarily lead to higher values of Rs, as depicted in 19a, where it is shown the influence of µr and CP on Rs obtained by Equation (1 resistance derived from Equation (24) is plotted as a reference). As can be seen, the of both µr and CP is weak, especially µr, resulting in = ≈ , although th ence of the armor seems to lead to slightly lower values for . This may be expla the armor alleviating the proximity effect in the sheaths. Similar results are obser Figure 19b, where the variation of Rs with µ' is represented for the 6 cables consid this study, confirming the weak influence of µr on . However, this does not necessarily lead to higher values of R s , as depicted in Figure 19a, where it is shown the influence of µ r and CP on R s obtained by Equation (18) (DC resistance derived from Equation (24) is plotted as a reference). As can be seen, the impact of both µ r and CP is weak, especially µ r , resulting in R s = R AC s ≈ R DC s , although the presence of the armor seems to lead to slightly lower values for R s . This may be explained as the armor alleviating the proximity effect in the sheaths. Similar results are observed in Figure 19b, where the variation of R s with µ is represented for the 6 cables considered in this study, confirming the weak influence of µ r on R s . Subsequently, since the armor has no important effects on Rs, the higher values observed for Is in Figure 17 suggest that Xs is then largely affected. This underscores the need of improving the formulation for computing Xs, not only for obtaining Is, but also for obtaining λ'1, where both CP and µr should play an important role. To this aim, the 3D ultrashortened FEM model developed in this work is paving the way to refine Equation (26) in the future next.

Summary
The main conclusions derived from all the analyses developed in this section are reflected in Figure 20, where the most relevant results, as well as the proposed solutions, are summarized for clarity. Their use for improving the experimental evaluation of Pa is developed next. Subsequently, since the armor has no important effects on R s , the higher values observed for I s in Figure 17 suggest that X s is then largely affected. This underscores the need of improving the formulation for computing X s , not only for obtaining I s , but also for obtaining λ 1 , where both CP and µ r should play an important role. To this aim, the 3D ultra-shortened FEM model developed in this work is paving the way to refine Equation (26) in the future next.

Summary
The main conclusions derived from all the analyses developed in this section are reflected in Figure 20, where the most relevant results, as well as the proposed solutions, are summarized for clarity. Their use for improving the experimental evaluation of P a is developed next. served for Is in Figure 17 suggest that Xs is then largely affected. This underscores of improving the formulation for computing Xs, not only for obtaining Is, but als taining λ'1, where both CP and µr should play an important role. To this aim, the shortened FEM model developed in this work is paving the way to refine Equatio the future next.

Summary
The main conclusions derived from all the analyses developed in this sectio flected in Figure 20, where the most relevant results, as well as the proposed solut summarized for clarity. Their use for improving the experimental evaluation of veloped next.

Estimation of Armor Losses
As indicated previously, it's challenging to determine accurately the differ ponents of the total losses in TCACs, especially through experimental measurem previous studies [26,28], different methods to reduce the measurement uncertai proposed, mainly based on measuring the total power dissipated in two three-cor one with armor and the other unarmored. During tests, different magnitudes a ured, such as the total power, the currents through the sheaths and conductors, a

Estimation of Armor Losses
As indicated previously, it's challenging to determine accurately the different components of the total losses in TCACs, especially through experimental measurements. In previous studies [26,28], different methods to reduce the measurement uncertainty were proposed, mainly based on measuring the total power dissipated in two three-core cables, one with armor and the other unarmored. During tests, different magnitudes are measured, such as the total power, the currents through the sheaths and conductors, as well as the DC resistance of metallic elements. Thus, P a is estimated by subtracting the power balance equations associated to both situations, where it is assumed that R c and R s are not influenced by the presence of the armor, hence having the same value in the two cables. The same applies for obtaining eddy losses in the sheaths by Equation (7). All this results in erroneous values for P a , since they collect all deviations derived from these assumptions.
To overcome these problems, new and more accurate expressions for evaluating all these parameters would be required, as commented earlier. Meanwhile, the uncertainties derived from these assumptions may be lessened by applying the approaches proposed in previous section, like Equations (21) and (23) for computing R c in TCACs, since the impact of the armor and the cores twisting on this parameter is considered. Additionally, the relative weight of circulating losses in the total sheath losses found previously (λ 1 ≈ 1.3λ 1 ) may help. All this, combined with the available experimental measurements (DC resistance and current in conductors and sheaths), allow us to propose the following expressions for estimating the sheath and conductor losses: where I c,meas and I s,meas represent the currents measured in conductors and sheaths, respectively. This way, the armor losses can be derived from the total power dissipated during the test (P meas total ) as P a = P meas total − P c − P s .
Even though the procedure includes experimental measurements, we have applied equivalent virtual tests in 750 simulations performed by means of the 3D-FEM model presented in this work, obtaining for each test the "virtual measurements" previously cited: R DC c , R DC s , I c,meas , I s,meas and P meas total . Regarding the 750 simulations, they consist in numerous situations obtained by modifying some of the main properties of a group of more than 50 TCACs ranging from 10 kV, 35 mm 2 to 275 kV, 2000 mm 2 . Changes in the number and thickness of the armor wires, the phase and armor lay length, as well as the sheath thickness and µ r = µ − jµ were applied to increase the number of cases. Further, copper and aluminum conductors were considered. Finally, several test cases include the increase in the electrical resistivity due to the temperature rise in the power cable. Only SB configuration has been considered in all the cases. Table 4 summarizes the range of the main numerical parameters employed in this study. Eventually, the values of R c and R s are computed by means of Equations (21)- (25). Having all these data, P c , P s and P a are obtained by means of Equations (28)- (30), and then compared to those obtained exclusively with FEM simulations, resulting in the error represented in Figure 21a. As can be easily observed, it is quite remarkable the good agreement shown between the estimated values and the FEM results, especially for P c and P s , since the maximum error barely exceeds 5% in just a few cases. In addition, a noticeable improvement is achieved in P a , with relative errors lower than 20% in most of the cases, improving the results previously observed in [16] where differences above 60% were usually obtained by the experimental procedures of [26,28].
Be aware that, in this case, a relative error below 20% is in practice negligible due to the small absolute value of P a , as can be observed in Figure 21b, where the result derived by FEM and [5] are represented, being the latter substantially overestimated. Precisely, the great differences between the absolute values of P a and P c , together with the highest value of µ r , are the cause of the highest errors observed in Figure 21a. In those situations, Equation (23) underestimates P c in about 1-2%. This results in 1-3 W/m of deviation that is later assigned to P a by Equation (30). However, this quantity is of the same order of magnitude of the actual value of P a obtained in those cases, thus causing the highest errors. Even so, less than 4% of the 750 cases analyzed in this work have a relative error greater than 20% in P a . value of µr, are the cause of the highest errors observed in Figure 21a. In those situations, Equation (23) underestimates Pc in about 1-2%. This results in 1-3 W/m of deviation that is later assigned to Pa by Equation (30). However, this quantity is of the same order of magnitude of the actual value of Pa obtained in those cases, thus causing the highest errors. Even so, less than 4% of the 750 cases analyzed in this work have a relative error greater than 20% in Pa. It is true that the application of the proposed approach in a laboratory setup will not be as accurate, mainly due to measurement and calibration errors. In any case, the approach presented in this work allows to avoid the erroneous assumptions considered in previous procedures [26,28], hence increasing the accuracy of the estimations. Furthermore, it also avoids the need of taking measurements in two cables, although this can be interesting for minimizing errors. It is true that the application of the proposed approach in a laboratory setup will not be as accurate, mainly due to measurement and calibration errors. In any case, the approach presented in this work allows to avoid the erroneous assumptions considered in previous procedures [26,28], hence increasing the accuracy of the estimations. Furthermore, it also avoids the need of taking measurements in two cables, although this can be interesting for minimizing errors.

Conclusions
This paper analyses the problem of properly quantifying and allocating the losses in the three metallic components (conductors, sheaths and magnetic armor) of TCACs commonly employed in submarine energy transmission. To this aim, an ultra-shortened and improved 3D-FEM model is proposed to perform fast and accurate analyses in 6 real cables. The influence of different aspects, such as the sheath bonding, conductor twisting, and the presence of a magnetic armor on the conductor and sheath losses is illustrated and quantified thanks to the new features of the proposed model.
From this analysis, it is concluded that both sheath and conductor losses increase with the cores and armor twisting, as well as with the magnetic permeability of the armor. Moreover, thanks to the new features implemented in the ultra-shortened 3D-FEM model, changes in the loss allocation can be visualized, concluding that the increase observed in the conductor losses may be identified as an increase in the proximity effect inside the conductor. On the other hand, the results also show how a magnetic armor leads to higher sheath currents, being this behavior a consequence of an increase in the sheath reactance.
As a result, different corrections are proposed for improving the estimation of the conductor and sheath losses, being applied in a new and simplified procedure suggested for the experimental evaluation of the armor losses, where the 3D-FEM model is employed as a "virtual laboratory". The procedure has been tested with more than 700 cable configurations showing excellent results, with differences below 5% in the sheath and conductor losses when compared to FEM results. Further, a good agreement is observed in the armor losses, with differences below 20%, substantially reducing the errors observed in previous experimental procedures.
Even though the procedure includes experimental measurements, the benefits of using FEM-based virtual tests are also remarked, such as reducing noticeably the experimental burden, and improving the knowledge of all the phenomena involved in this type of cable. In this sense, the proposed 3D-FEM model is currently employed in finding new and more accurate expressions, not only for evaluating the cable losses, but also for obtaining its electrical parameters, being all this in line with the remarkable efforts of on-going CIGRÉ working groups B1.64 and B1.56.