<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.3 20210610//EN" "JATS-journalpublishing1-3-mathml3.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" dtd-version="1.3" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Smart Grids</journal-id>
<journal-title-group>
<journal-title>Frontiers in Smart Grids</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Smart Grids</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">2813-4311</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/frsgr.2026.1623347</article-id>
<article-version article-version-type="Version of Record" vocab="NISO-RP-8-2008"/>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Original Research</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Analysis of theoretical line loss in metro distribution network</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Zhao</surname> <given-names>Ying</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Chen</surname> <given-names>Zhe</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhang</surname> <given-names>Zhuo</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Lu</surname> <given-names>Yuxin</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
<uri xlink:href="https://loop.frontiersin.org/people/2388800"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhao</surname> <given-names>Yun</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Liu</surname> <given-names>Xipeng</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label><institution>Guangzhou Power Supply Bureau,  Guangdong Power Grid Co.,  Ltd.</institution>, <city>Guangzhou</city>, <country country="cn">China</country></aff>
<aff id="aff2"><label>2</label><institution>CSG Electric Power Research Institute</institution>, <city>Guangzhou</city>, <country country="cn">China</country></aff>
<aff id="aff3"><label>3</label><institution>Guangdong Provincial Key Laboratory of Intelligent Measurement and Advanced Metering of Power Grid</institution>, <city>Guangzhou</city>, <country country="cn">China</country></aff>
<author-notes>
<corresp id="c001"><label>&#x0002A;</label>Correspondence: Yuxin Lu, <email xlink:href="mailto:luyx2@csg.cn">luyx2@csg.cn</email></corresp>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="2026-03-12">
<day>12</day>
<month>03</month>
<year>2026</year>
</pub-date>
<pub-date publication-format="electronic" date-type="collection">
<year>2026</year>
</pub-date>
<volume>5</volume>
<elocation-id>1623347</elocation-id>
<history>
<date date-type="received">
<day>05</day>
<month>05</month>
<year>2025</year>
</date>
<date date-type="rev-recd">
<day>01</day>
<month>12</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>29</day>
<month>01</month>
<year>2026</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2026 Zhao, Chen, Zhang, Lu, Zhao and Liu.</copyright-statement>
<copyright-year>2026</copyright-year>
<copyright-holder>Zhao, Chen, Zhang, Lu, Zhao and Liu</copyright-holder>
<license>
<ali:license_ref start_date="2026-03-12">https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License (CC BY)</ext-link>. The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</license-p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>Novel cable traction power supply can eliminate phase splits and improve regenerative energy utilization in metro systems. However, the cable-based traction network has a relatively complex structure, making internal energy consumption and train-network line loss evaluation insufficiently studied.</p></sec>
<sec>
<title>Methods</title>
<p>To quantify train-network line losses, an equivalent solid-circuit model of the traction power supply system is established using an external cascade strategy. Combined with actual train operation load modeling, the traction transformer output power and network losses are computed and compared under both traction and regenerative braking operating conditions.</p></sec>
<sec>
<title>Results</title>
<p>The proposed modeling-and-calculation framework enables accurate evaluation of power output and line losses across operating modes, and effectively captures the electrical characteristics of the cable traction energy supply system.</p></sec>
<sec>
<title>Discussion</title>
<p>The developed theoretical line-loss calculation approach provides a quantitative basis for energy-consumption assessment and operational analysis of cable traction networks, supporting performance evaluation and planning of metro traction power supply systems.</p></sec></abstract>
<kwd-group>
<kwd>cable traction network</kwd>
<kwd>distribution network</kwd>
<kwd>lineloss calculation</kwd>
<kwd>subway systems</kwd>
<kwd>transformer output power</kwd>
</kwd-group>
<funding-group>
<funding-statement>The author(s) declared that financial support was received for this work and/or its publication. This work was financially supported by the key Science and Technology Project of China Southern Power Grid Co., Ltd. (Research and Application of Key Technologies of Theoretical Line Loss Calculation and Analysis Software for Distribution Network Considering Multi-node Distributed Energy Accesses, 030100KC22120004/030100KC22120006/ 030100KC22120008/030100KC22120010). China Southern Power Grid Co., Ltd. was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.</funding-statement>
</funding-group>
<counts>
<fig-count count="9"/>
<table-count count="2"/>
<equation-count count="21"/>
<ref-count count="18"/>
<page-count count="9"/>
<word-count count="4450"/>
</counts>
<custom-meta-group>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Grid Efficiency</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<label>1</label>
<title>Introduction</title>
<p>Conventional DC power supply systems in urban rail transit, although capable of supporting frequent start-stop operations, face persistent challenges such as stray currents causing infrastructure corrosion, low efficiency in regenerative braking energy utilization, negative sequence currents affecting power quality, and phase splits that limit train speed (<xref ref-type="bibr" rid="B4">Dolara et al., 2012</xref>). To overcome these issues, a novel AC traction power supply scheme has been proposed. This scheme integrates new contact supply devices, AC-DC-AC locomotives, and co-phase compensation devices, thereby eliminating phase splits and stray currents, reducing construction costs, and improving the utilization of regenerative braking energy (<xref ref-type="bibr" rid="B16">Steczek et al., 2023</xref>; <xref ref-type="bibr" rid="B3">Chen et al., 2016</xref>).</p>
<p>For the effective operation of this system, accurate theoretical line loss calculation is essential (<xref ref-type="bibr" rid="B2">Cai et al., 1996</xref>; <xref ref-type="bibr" rid="B6">Goodman and Kulworawanichpong, 2002</xref>). It provides the basis for optimizing operational strategies, planning network upgrades, and guiding equipment maintenance. However, in metro systems the stochastic and dynamic nature of train loads has a significant impact on line losses. As trains accelerate, decelerate, and run at constant speed across different substations, the traction load continuously shifts along the network, complicating the loss evaluation (<xref ref-type="bibr" rid="B18">Zhang Y. et al., 2025</xref>).</p>
<p>Existing research has addressed traction power supply in various forms, including direct supply, BT supply, and AT supply (<xref ref-type="bibr" rid="B11">Kulworawanichpong, 2015</xref>), with equivalent circuit models often employed for engineering analysis (<xref ref-type="bibr" rid="B1">Al-Khalidi and Kalam, 2008</xref>). Line loss calculation methods in distribution networks can be broadly categorized into: (1) traditional approaches based on substation output data, which are simple but neglect time-domain characteristics (<xref ref-type="bibr" rid="B12">Li et al., 2020</xref>); (2) methods using FTU/TTU measurements, which are convenient but cannot directly capture the effect of load variations (Zhang W. et al., <xref ref-type="bibr" rid="B17">2025</xref>; <xref ref-type="bibr" rid="B14">Ni et al., 2019</xref>); (3) real-time data-driven methods, which improve timeliness but remain limited in analyzing dynamic load impacts (<xref ref-type="bibr" rid="B7">Gulakhmadov et al., 2022</xref>).</p>
<p>This paper aims to extend these studies to urban rail traction networks by explicitly considering the dynamic characteristics of AC-DC-AC locomotives. By combining traction characteristic curves with train operation diagram is established. The main contributions of this paper are summarized as follows: (1) <bold>A traction-oriented cascade equivalent model</bold> is established, covering the substation&#x02013;feeder&#x02013;section&#x02013;train structure of metro AC traction networks. (2) <bold>An improved impedance aggregation method</bold> for multi-core cables and traction transformers is proposed, enabling position-dependent equivalent impedance evaluation. (3) <bold>A time-sequential loss calculation framework</bold> is developed, coupling dynamic train load, nodal-voltage iteration, and component-level loss decomposition along the line.</p>
<p>The rest of this paper is organized as follows. System structure and network circuit equivalence process are introduced in Section 2. Section 3 presents a power distribution calculation model. Under the condition that load is in traction and regenerative braking while single train and multi-train are working, power distribution in system is explored in Section 4. At last, corresponding conclusions are induced in Section 5.</p></sec>
<sec id="s2">
<label>2</label>
<title>Traction power supply system</title>
<sec>
<label>2.1</label>
<title>System structure</title>
<p>As illustrated in <xref ref-type="fig" rid="F1">Figure 1</xref>, the traction power supply system consists of a central substation (SS), traction cables, single-phase traction transformers, an overhead contact system (OCS), and a current return line (CRL). The central substation connects to the 110 kV three-phase grid and delivers 35 kV single-phase power through traction cables. Each traction substation is equipped with a single-phase transformer (<italic>T</italic><sub><italic>il</italic>1</sub>, <italic>T</italic><sub><italic>il</italic>2</sub>, <italic>T</italic><sub><italic>il</italic>3</sub>), stepping down the voltage to the catenary level for train operation (<xref ref-type="bibr" rid="B13">Li, 2016</xref>).</p>
<fig position="float" id="F1">
<label>Figure 1</label>
<caption><p>Schematic of traction power supply system.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0001.tif">
<alt-text content-type="machine-generated">Schematic diagram showing a power grid distributing electricity through a substation labeled SS to multiple transformers Ti1, Ti2, and Ti3, which supply power to three parallel railway trains labeled PSS via overhead contact system OCS and current return line CRL.</alt-text>
</graphic>
</fig>
<p>The central substation employs combined co-phase power supply technology, consisting of a traction transformer and a back-to-back negative sequence compensation device (NCD). The NCD transfers active power from the grid side to the traction side while compensating harmonics and reactive power (<xref ref-type="bibr" rid="B15">Shu et al., 2013</xref>). In practice, 35 kV two-core cables are used for traction power, while three-core cables are reserved for auxiliary supply (<xref ref-type="bibr" rid="B9">Huang et al., 2023</xref>). The traction network supplies locomotives with same-phase power along the entire line, improving supply reliability and reducing the installed capacity of traction transformers.</p>
</sec>
<sec>
<label>2.2</label>
<title>Network equivalent circuit</title>
<sec>
<label>2.2.1</label>
<title>Traction cables</title>
<p>Urban rail AC traction systems typically employ 35 kV cables, either two single-core or one double-core type. A single-core armored cable can be modeled as a linear passive element at the operating frequency, and a uniform line segment is equivalent to a &#x003A0;-type circuit, as shown in <xref ref-type="fig" rid="F2">Figure 2</xref>.</p>
<fig position="float" id="F2">
<label>Figure 2</label>
<caption><p>Equivalent &#x003A0;-type circuit.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0002.tif">
<alt-text content-type="machine-generated">Block diagram showing a central impedance matrix Z_L connected between two nodes, each node connected to ground via identical admittance matrices Y_L divided by two, representing a symmetrical multi-port network with matrix notations.</alt-text>
</graphic>
</fig>
<p>The impedance and admittance matrices of the cable are expressed as</p>
<disp-formula id="EQ1"><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>Z</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>11</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mo>&#x022EF;</mml:mo></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>&#x022EE;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x022F1;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x022EE;</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mo>&#x022EF;</mml:mo></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>Y</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Y</mml:mi></mml:mrow><mml:mrow><mml:mn>11</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mo>&#x022EF;</mml:mo></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mi>Y</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>&#x022EE;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x022F1;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x022EE;</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Y</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mo>&#x022EF;</mml:mo></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mi>Y</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(1)</label></disp-formula>
<p>Here, <italic>Z</italic><sub><italic>ii</italic></sub> and <italic>Y</italic><sub><italic>ii</italic></sub> denote the self-impedance and self-admittance of phase <italic>i</italic>, while <italic>Z</italic><sub><italic>ij</italic></sub> represents mutual impedance. Due to shielding, <italic>Y</italic><sub><italic>ij</italic></sub> &#x0003D; 0. The capacitances <italic>C</italic><sub><italic>ics</italic></sub>, <italic>C</italic><sub><italic>isa</italic></sub>, <italic>C</italic><sub><italic>iae</italic></sub> describe coupling between core&#x02013;shield, shield&#x02013;armor, and armor&#x02013;ground (<xref ref-type="bibr" rid="B5">Ghassemi, 2019</xref>). These parameters allow calculation of impedance for arbitrary line length and configuration, enabling accurate current and voltage distribution analysis in traction networks. This equivalent model also provides the basis for subsequent power flow and line loss calculations.</p></sec>
<sec>
<label>2.2.2</label>
<title>Central substation</title>
<p>As shown in <xref ref-type="fig" rid="F3">Figure 3</xref>, the traction port consists of the traction transformer and the co-phase compensation device (CPD), which share the load and ideally eliminate the negative-sequence current. The cable-side current satisfies</p>
<disp-formula id="EQ2"><mml:math id="M2"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(2)</label></disp-formula>
<p>where <italic>K</italic><sub><italic>N</italic></sub>&#x02208;[0, 1] is the compensation factor.</p>
<fig position="float" id="F3">
<label>Figure 3</label>
<caption><p>Power supply structure of the traction port.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0003.tif">
<alt-text content-type="machine-generated">Diagram showing two parallel electrical lines entering a block labeled &#x0201C;CPD.&#x0201D; Voltages u_r and u_c are marked on the left and right, with currents i_r and i_c flowing toward the block, and corresponding negative currents flowing away. Current i_L branches downward at two points along the lines, connecting to terminals labeled PS and PR below.</alt-text>
</graphic>
</fig>
<p>Applying Thevenin&#x00027;s theorem, the traction port has open-circuit voltage equal to the rated secondary voltage, and equivalent resistance</p>
<disp-formula id="EQ3"><mml:math id="M3"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>S</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(3)</label></disp-formula>
<p>which depends only on transformer and system short-circuit capacity, and is independent of cable and load.</p>
</sec>
</sec>
<sec>
<label>2.3</label>
<title>Traction supply system equivalent model</title>
<p>The AC traction power network of urban rail transit is essentially a 35 kV cable-fed system. The central substation connects to the contact line through single-phase traction transformers in parallel, supplying locomotives. Since the network involves multiple voltage levels, all circuits are converted to a unified 35 kV voltage level for modeling purposes.</p>
<p>As shown in <xref ref-type="fig" rid="F4">Figure 4</xref>, the traction network can be represented by a multi-loop equivalent circuit. The central substation is modeled by an ideal voltage source &#x01EBC;<sub>0</sub> with internal impedance <italic>Z</italic><sub><italic>s</italic></sub>. The line is divided into two parts: - a <italic>local loop</italic>, which contains the locomotive and the nearby traction substations, and - a <italic>distal loop</italic>, which represents the remaining sections equivalently.</p>
<fig position="float" id="F4">
<label>Figure 4</label>
<caption><p>Equivalent traction supply circuit.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0004.tif">
<alt-text content-type="machine-generated">Single-line schematic diagram illustrating a traction power supply system for railway electrification, including a central substation, local and distal loops, electrical impedances, current paths, catenary and rail conductors, and a train represented with relevant current and distance parameters.</alt-text>
</graphic>
</fig>
<p>In the local loop, each traction substation is represented by its transformer impedance (<italic>Z</italic><sub><italic>T</italic>1</sub>, <italic>Z</italic><sub><italic>T</italic>3</sub>, <italic>Z</italic><sub><italic>T</italic>4</sub>, &#x02026;). The catenary and return conductors are represented by impedances <italic>Z</italic><sub><italic>OCS</italic></sub> and <italic>Z</italic><sub><italic>CS</italic></sub>. The locomotive is modeled as a current source <italic>i</italic><sub><italic>n</italic></sub> located at position <italic>x</italic>, with the distance between adjacent substations denoted as <italic>D</italic> and the total line length as <italic>L</italic>. The distal loop is reduced to an equivalent impedance <inline-formula><mml:math id="M4"><mml:msubsup><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> obtained by recursive parallel conversion of the remote sections.</p>
<p>This equivalent model explicitly describes the current distribution across substations (<italic>I</italic><sub>1</sub>, <italic>I</italic><sub>2</sub>, &#x02026;) and forms the basis for subsequent derivation of line loss equations.</p>
</sec>
<sec>
<label>2.4</label>
<title>Line loss calculation method</title>
<p>Consider <italic>n</italic> traction substations uniformly distributed along the line. Let <italic>D</italic> denote the spacing between two adjacent substations along the track. Let <italic>l</italic> be the longitudinal position of the locomotive measured from the line origin (0 &#x02264; <italic>l</italic> &#x02264; <italic>nD</italic>). Thus the locomotive lies within the <italic>m</italic>-th supply section, where <italic>m</italic> &#x0003D; &#x02308;<italic>l</italic>/<italic>D</italic>&#x02309;, and the relative position inside this section is <italic>x</italic> &#x0003D; <italic>l</italic>&#x02212;(<italic>m</italic>&#x02212;1)<italic>D</italic>, satisfying 0 &#x02264; <italic>x</italic> &#x02264; <italic>D</italic>.</p>
<p>Let [<italic>I</italic><sub>1</sub>, <italic>I</italic><sub>2</sub>, &#x02026;, <italic>I</italic><sub><italic>n</italic></sub>] be the vector of substation currents. The current distribution coefficient is defined as</p>
<disp-formula id="EQ4"><mml:math id="M5"><mml:mrow><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
<p>where <italic>I</italic><sub><italic>L</italic></sub> is the total traction load current. The distal loop can be equivalently reduced to <inline-formula><mml:math id="M6"><mml:msubsup><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, obtained by recursive parallel conversion of transformer and contact-line impedances (<xref ref-type="bibr" rid="B10">Kim et al., 2010</xref>).</p>
<p>The nodal equations of the system can be expressed compactly as</p>
<disp-formula id="EQ5"><mml:math id="M7"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>A</mml:mi><mml:mo>&#x000B7;</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(4)</label></disp-formula>
<p>where matrix <italic>A</italic> is determined by system parameters and locomotive position <italic>l</italic>. Normalizing by <italic>I</italic><sub><italic>L</italic></sub> yields the current distribution vector</p>
<disp-formula id="EQ6"><mml:math id="M8"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>A</mml:mi><mml:mo>&#x000B7;</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(5)</label></disp-formula>
<p>Given locomotive position <italic>l</italic>, the local circuit index <italic>m</italic> and distance <italic>x</italic> within that circuit are obtained from</p>
<disp-formula id="EQ7"><mml:math id="M9"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>m</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>&#x02308;</mml:mo><mml:mrow><mml:mi>l</mml:mi><mml:mo>/</mml:mo><mml:mi>D</mml:mi></mml:mrow><mml:mo>&#x02309;</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:mi>x</mml:mi><mml:mo>=</mml:mo><mml:mi>l</mml:mi><mml:mo>-</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>m</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi>D</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(6)</label></disp-formula>
<p>with <italic>D</italic> the spacing between substations. Thus the equivalent impedance at the locomotive is</p>
<disp-formula id="EQ8"><mml:math id="M10"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>Z</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>C</mml:mi><mml:mi>S</mml:mi></mml:mrow></mml:msub><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>h</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mi>h</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mfrac><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>D</mml:mi></mml:mrow></mml:mfrac><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>O</mml:mi><mml:mi>C</mml:mi><mml:mi>S</mml:mi></mml:mrow></mml:msub><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:munderover></mml:mstyle><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(7)</label></disp-formula>
<p>and the corresponding load voltage</p>
<disp-formula id="EQ9"><mml:math id="M11"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>E</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mi>Z</mml:mi><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(8)</label></disp-formula>
<p>Finally, the line loss is obtained as the sum of resistive losses of cables and transformers and reactive losses due to power factor. The procedure is summarized as:</p>
<list list-type="order">
<list-item><p>Obtain phasors of three-phase voltage and current from transformer measurements and convert to the HV side;</p></list-item>
<list-item><p>Compute equivalent impedance <italic>Z</italic> from network parameters and locomotive position <italic>l</italic>;</p></list-item>
<list-item><p>Determine current distribution {<italic>k</italic><sub><italic>i</italic></sub>} and calculate losses of each load;</p></list-item>
<list-item><p>Sum all load losses to obtain the total network loss and the contribution of each substation.</p></list-item>
</list></sec></sec>
<sec id="s3">
<label>3</label>
<title>Dynamic traction load modeling based on train operation diagrams</title>
<sec>
<label>3.1</label>
<title>Methodology for establishing a traction load model</title>
<p>This section establishes a dynamic traction load model by describing the forces acting on the locomotive and converting them into power demand on the contact wire side.</p>
<sec>
<label>3.1.1</label>
<title>Traction and braking forces</title>
<p>Locomotives are subject to traction, braking, and resistance forces (<xref ref-type="fig" rid="F5">Figure 5</xref>), where traction and braking do not occur simultaneously. For the CRH2 locomotive (<xref ref-type="bibr" rid="B8">Huang et al., 2016</xref>), the traction force <italic>TE</italic> (kN) and regenerative braking force <italic>F</italic><sub><italic>b</italic></sub> are</p>
<disp-formula id="EQ10"><mml:math id="M12"><mml:mtable class="eqnarray" columnalign="right"><mml:mtr><mml:mtd><mml:mi>T</mml:mi><mml:mi>E</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mn>176</mml:mn><mml:mo>-</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>36</mml:mn><mml:mi>v</mml:mi><mml:mo>,</mml:mo></mml:mtd><mml:mtd><mml:mi>v</mml:mi><mml:mo>&#x02264;</mml:mo><mml:mn>125</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>16250</mml:mn><mml:mo>/</mml:mo><mml:mi>v</mml:mi><mml:mo>,</mml:mo></mml:mtd><mml:mtd><mml:mi>v</mml:mi><mml:mo>&#x0003E;</mml:mo><mml:mn>125</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mn>172</mml:mn><mml:mo>,</mml:mo></mml:mtd><mml:mtd><mml:mi>v</mml:mi><mml:mo>&#x02264;</mml:mo><mml:mn>70</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>172</mml:mn><mml:mo>-</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>v</mml:mi><mml:mo>-</mml:mo><mml:mn>70</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>/</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo></mml:mtd><mml:mtd><mml:mi>v</mml:mi><mml:mo>&#x0003E;</mml:mo><mml:mn>70</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(9)</label></disp-formula>
<fig position="float" id="F5">
<label>Figure 5</label>
<caption><p>Forces acting on a locomotive.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0005.tif">
<alt-text content-type="machine-generated">Diagram showing forces on a vehicle moving uphill and downhill on inclined planes. Labeled vectors indicate gravitational forces, velocity, angle theta, normal force, and driving or resisting forces for both uphill and downhill directions.</alt-text>
</graphic>
</fig>
</sec>
<sec>
<label>3.1.2</label>
<title>Resistance</title>
<p>The running resistance consists of a basic term and additional terms due to gradients, curves, and tunnels:</p>
<disp-formula id="EQ11"><mml:math id="M14"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>w</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>A</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mi>B</mml:mi><mml:mi>v</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mi>C</mml:mi><mml:msup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>w</mml:mi></mml:mrow><mml:mrow><mml:mi>g</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>w</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>w</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(10)</label></disp-formula>
<p>with <italic>w</italic><sub><italic>g</italic></sub> &#x0003D; 1000tan(&#x003B8;), <italic>w</italic><sub><italic>c</italic></sub> &#x0003D; <italic>K</italic><sub><italic>c</italic></sub>/<italic>R</italic>, and <italic>w</italic><sub><italic>t</italic></sub> &#x0003D; 0.00013<italic>L</italic><sub><italic>t</italic></sub>. Thus the net force on the locomotive is</p>
<disp-formula id="EQ12"><mml:math id="M15"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>F</mml:mi><mml:mo>=</mml:mo><mml:mi>T</mml:mi><mml:mi>E</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi>w</mml:mi><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(11)</label></disp-formula>
</sec>
</sec>
<sec>
<label>3.2</label>
<title>Power calculation for locomotive operation</title>
<p>The wheelset output power is</p>
<disp-formula id="EQ13"><mml:math id="M16"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>P</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>F</mml:mi><mml:mo>&#x000B7;</mml:mo><mml:mi>v</mml:mi><mml:mo>&#x000B7;</mml:mo><mml:mn>1000</mml:mn></mml:mrow><mml:mrow><mml:mn>3</mml:mn><mml:mo>.</mml:mo><mml:mn>6</mml:mn></mml:mrow></mml:mfrac><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(12)</label></disp-formula>
<p>and the equivalent power on the contact line side, considering transmission efficiency &#x003B7;, auxiliary load <italic>P</italic><sub><italic>aux</italic></sub>, and power factor angle &#x003C6;<sub><italic>T</italic></sub>, is</p>
<disp-formula id="EQ14"><mml:math id="M17"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mtext class="textrm" mathvariant="normal">Traction:</mml:mtext></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x003B7;</mml:mi></mml:mrow></mml:mfrac><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi><mml:mi>u</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:mi>Q</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>/</mml:mo><mml:mo class="qopname">tan</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext class="textrm" mathvariant="normal">Coasting:</mml:mtext></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi><mml:mi>u</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:mi>Q</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>/</mml:mo><mml:mo class="qopname">tan</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext class="textrm" mathvariant="normal">Braking:</mml:mtext></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x003B7;</mml:mi></mml:mrow></mml:mfrac><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi><mml:mi>u</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:mi>Q</mml:mi><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mi>e</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>/</mml:mo><mml:mo class="qopname">tan</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(13)</label></disp-formula>
</sec>
<sec>
<label>3.3</label>
<title>Dynamic load generation procedure</title>
<p>Using the train operation diagram, the dynamic load profile is obtained by:</p>
<list list-type="order">
<list-item><p>Integrating operational and locomotive data (traction/braking curves, power factor, <italic>P</italic><sub><italic>aux</italic></sub>, &#x003B7;, running time and speed).</p></list-item>
<list-item><p>Determining position <italic>X</italic><sub><italic>n</italic></sub> and speed <italic>v</italic><sub><italic>n</italic></sub> at time <italic>T</italic><sub><italic>n</italic></sub>, then computing force <italic>F</italic>.</p></list-item>
<list-item><p>Calculating active power <italic>P</italic><sub><italic>n</italic></sub> and apparent power <italic>S</italic><sub><italic>n</italic></sub> using <xref ref-type="disp-formula" rid="EQ14">Equation 13</xref>, forming data rows [<italic>T</italic><sub><italic>n</italic></sub>, <italic>X</italic><sub><italic>n</italic></sub>, <italic>v</italic><sub><italic>n</italic></sub>, <italic>P</italic><sub><italic>n</italic></sub>, <italic>S</italic><sub><italic>n</italic></sub>].</p></list-item>
</list>
<p>This process yields the time-varying traction load along the line.</p></sec></sec>
<sec id="s4">
<label>4</label>
<title>Example calculation of metro distribution network line losses</title>
<sec>
<label>4.1</label>
<title>Calculation procedure of line losses</title>
<p>At each sampling instant <italic>T</italic><sub><italic>n</italic></sub>, the algorithm proceeds as follows:</p>
<list list-type="order">
<list-item><p><bold>Input and preprocessing</bold>. The equivalent impedance parameters are imported, including the main substation impedance, cable unit-length impedance, transformer equivalent impedance, and contact-line impedance. Based on Section 2, the branch-coefficient matrices [<italic>X, K</italic>] and the equivalent impedance matrix [<italic>X, Z</italic>] are generated.</p></list-item>
<list-item><p><bold>Time-step initialization</bold>. The simulation clock is initialized at <italic>T</italic><sub>0</sub>. For the current locomotive position <italic>X</italic><sub><italic>n</italic></sub>, the corresponding complex traction power <italic>S</italic><sub><italic>n</italic></sub> is obtained from the load model.</p></list-item>
<list-item><p><bold>Equivalent-impedance extraction</bold>. Using the position <italic>X</italic><sub><italic>n</italic></sub> and the matrix [<italic>X, Z</italic>], the equivalent impedance <italic>Z</italic><sub><italic>n</italic></sub> seen from the locomotive is computed.</p></list-item>
<list-item><p><bold>Initialization of nodal voltage and current</bold>. The voltage and current are initialized as</p>
<p><disp-formula id="EQ15"><mml:math id="M18"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mn>35000</mml:mn><mml:mtext class="textrm" mathvariant="normal">V</mml:mtext><mml:mo>,</mml:mo><mml:mtext>&#x02003;&#x000A0;</mml:mtext><mml:msubsup><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p></list-item>
<list-item><p><bold>Iterative voltage update</bold>. For iteration index <italic>k</italic>:</p>
<p><disp-formula id="EQ16"><mml:math id="M19"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mn>35000</mml:mn><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:msubsup><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>If <italic>k</italic> &#x0003D; <italic>k</italic><sub>max</sub>, the point is marked as not convergent and the algorithm proceeds to the next time step. Otherwise, if</p>
<p><disp-formula id="EQ17"><mml:math id="M20"><mml:mrow><mml:mo class="qopname">max</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mi>&#x0211C;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup></mml:mstyle><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>|</mml:mo><mml:mo>,</mml:mo><mml:mo>|</mml:mo><mml:mi>&#x02111;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0003C;</mml:mo><mml:mi>&#x003B5;</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>convergence is declared.</p></list-item>
<list-item><p><bold>Current update</bold>. Once <inline-formula><mml:math id="M21"><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup></mml:math></inline-formula> converges, the current is updated as</p>
<p><disp-formula id="EQ18"><mml:math id="M22"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msubsup><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p></list-item>
<list-item><p><bold>Acquisition of additional nodal information</bold>. Using the converged voltage and current, the voltages, currents, and losses of all other target nodes along the line are computed. These are used for the loss-calculation steps described in Section 4.2.</p></list-item>
<list-item><p><bold>Time advancement</bold>. The simulation time is updated by</p>
<p><disp-formula id="EQ19"><mml:math id="M23"><mml:mrow><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mtext>&#x00394;</mml:mtext><mml:mi>T</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p></list-item></list>
<p>If <italic>T</italic><sub><italic>n</italic>&#x0002B;1</sub>&#x0003E;<italic>T</italic>, the simulation terminates; otherwise, the next time step is computed.</p>
</sec>
<sec>
<label>4.2</label>
<title>Raw data</title>
<p>This section selects a line from that in Section 2 for validation. Based on the model of the traction power supply system, <xref ref-type="fig" rid="F6">Figure 6</xref> shows the traction load distribution at a specific moment on a particular substation along the line. In this system, the primary side of the traction substation is connected to a 220kV grid, with the system short-circuit capacity at the connection point being 4,000 MVA. The secondary side voltage level is 35kV. The traction transformer adopts a Vv connection, with the phase sequence connected as ABC. The transformer capacity is 80 MVA, and the percentage of short-circuit impedance is 10.5%. Each individual traction transformer has a capacity of 10 MVA, with a percentage impedance of 8.5%. In this simulation, the contact wire is selected as the power supply device, and the voltage level is set to 3.5kV. To ensure that the traction power supply system meets operational requirements, the local circuit length is set to 5 km. Considering various on-site factors, the simulation parameters are set as shown in <xref ref-type="table" rid="T1">Table 1</xref>, with traction transformers placed at intervals of 0 km, 5 km, 10 km, 15 km, and 20 km.</p>
<fig position="float" id="F6">
<label>Figure 6</label>
<caption><p>Traction load speed-time diagram.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0006.tif">
<alt-text content-type="machine-generated">Line graph showing train speed in kilometers per hour versus time in seconds. Speed fluctuates repeatedly between near zero and approximately eighty-five kilometers per hour in a repeated pattern throughout the timeline.</alt-text>
</graphic>
</fig>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Line parameters.</p></caption>
<table frame="box" rules="all">
<thead>
<tr>
<th valign="top" align="left"><bold>Item</bold></th>
<th valign="top" align="center"><bold>Parameter</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Equivalent external power voltage</td>
<td valign="top" align="center">36 kV</td>
</tr>
<tr>
<td valign="top" align="left">Main substation equivalent impedance</td>
<td valign="top" align="center">0.05&#x0002B;<italic>j</italic>0.31 &#x003A9;</td>
</tr>
<tr>
<td valign="top" align="left">Cable unit length impedance</td>
<td valign="top" align="center">0.1137&#x0002B;<italic>j</italic>0.8117 &#x003A9;/km</td>
</tr>
<tr>
<td valign="top" align="left">Traction transformer equivalent impedance</td>
<td valign="top" align="center">0.48&#x0002B;<italic>j</italic>7.75 &#x003A9;</td>
</tr>
<tr>
<td valign="top" align="left">Contact supply device equivalent impedance</td>
<td valign="top" align="center">0.1635&#x0002B;<italic>j</italic>0.7443 &#x003A9;/km</td>
</tr></tbody>
</table>
</table-wrap>
<p>At present, there is a lack of actual data on the AC traction power supply systems of urban rail transit. Therefore, this paper selects the CJ6 locomotive as the operating locomotive, adopts the operation curve of the locomotive on this line calculating line losses. The speed-time curve of a single load operation on this line is shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. Validation is performed for both traction and braking conditions.</p>
</sec>
<sec>
<label>4.3</label>
<title>Calculation of losses generated by the load</title>
<p>When the resistance and power supply model are determined, the branch coefficients of each traction substation are obtained based on the method from Section 2. The variation of branch coefficients with train position is shown in <xref ref-type="fig" rid="F7">Figure 7</xref>. By combining line information, the equivalent impedance <italic>Z</italic> between each position and the 35 kV external power supply is calculated according to (7). After obtaining the distribution of <italic>Z</italic>, the power flow is solved using the equivalent load model from Section 3, considering the train&#x00027;s operation along the entire line. Based on the locomotive voltage obtained from the power flow analysis, the power characteristics of the main substation and traction substations, as well as the system losses, are evaluated.</p>
<fig position="float" id="F7">
<label>Figure 7</label>
<caption><p>Power distribution among traction substations along the line.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0007.tif">
<alt-text content-type="machine-generated">Line graph showing the output percentages of five traction transformers, each peaking sequentially and symmetrically along the locomotive position axis from zero to approximately seventeen kilometers, with overlapping triangular patterns for each transformer.</alt-text>
</graphic>
</fig>
<p>As shown in <xref ref-type="fig" rid="F7">Figure 7</xref> and <xref ref-type="table" rid="T2">Table 2</xref>, traction substations mainly supply power to their adjacent local loops during train operation. The power distribution exhibits a Gaussian-like pattern, being higher in the middle and lower at both ends, with maximum output when the train passes near a substation. Since substations 1 and 5 only feed one adjacent section, while substations 2&#x02013;4 supply two neighboring sections, the average output power of substations 2&#x02013;4 is nearly twice that of substations 1 and 5. Furthermore, as illustrated in <xref ref-type="fig" rid="F8">Figure 8</xref>, because the reactance of the traction network equivalent impedance is much larger than its resistance, reactive power losses dominate over active power losses.</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Statistical power characteristics of each traction substation.</p></caption>
<table frame="box" rules="all">
<thead>
<tr>
<th valign="top" align="left"><bold>Metric</bold></th>
<th valign="top" align="center"><bold>Substation 1</bold></th>
<th valign="top" align="center"><bold>Substation 2</bold></th>
<th valign="top" align="center"><bold>Substation 3</bold></th>
<th valign="top" align="center"><bold>Substation 4</bold></th>
<th valign="top" align="center"><bold>Substation 5</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Average Apparent Power (kVA)</td>
<td valign="top" align="center">180.5</td>
<td valign="top" align="center">455.2</td>
<td valign="top" align="center">318.6</td>
<td valign="top" align="center">248.7</td>
<td valign="top" align="center">132.1</td>
</tr>
<tr>
<td valign="top" align="left">Maximum Apparent Power (kVA)</td>
<td valign="top" align="center">4,050.7</td>
<td valign="top" align="center">4,690.3</td>
<td valign="top" align="center">3,790.4</td>
<td valign="top" align="center">3,980.5</td>
<td valign="top" align="center">3,885.7</td>
</tr>
<tr>
<td valign="top" align="left">Minimum Apparent Power (VA)</td>
<td valign="top" align="center">46.2</td>
<td valign="top" align="center">0.521</td>
<td valign="top" align="center">21.8</td>
<td valign="top" align="center">0.310</td>
<td valign="top" align="center">0.005</td>
</tr>
<tr>
<td valign="top" align="left">Average Active Power (kW)</td>
<td valign="top" align="center">75.8</td>
<td valign="top" align="center">172.4</td>
<td valign="top" align="center">160.9</td>
<td valign="top" align="center">112.7</td>
<td valign="top" align="center">40.2</td>
</tr>
<tr>
<td valign="top" align="left">Maximum Active Power (kW)</td>
<td valign="top" align="center">3,988.6</td>
<td valign="top" align="center">4,620.5</td>
<td valign="top" align="center">3,720.1</td>
<td valign="top" align="center">3,815.8</td>
<td valign="top" align="center">3,779.4</td>
</tr>
<tr>
<td valign="top" align="left">Minimum Active Power (kW)</td>
<td valign="top" align="center">&#x02013;2,950.2</td>
<td valign="top" align="center">&#x02013;4,420.6</td>
<td valign="top" align="center">&#x02013;3,215.8</td>
<td valign="top" align="center">&#x02013;3,720.4</td>
<td valign="top" align="center">&#x02013;3,490.2</td>
</tr>
<tr>
<td valign="top" align="left">Average Reactive Power (kVar)</td>
<td valign="top" align="center">20.3</td>
<td valign="top" align="center">38.7</td>
<td valign="top" align="center">33.5</td>
<td valign="top" align="center">24.1</td>
<td valign="top" align="center">9.3</td>
</tr>
<tr>
<td valign="top" align="left">Maximum Reactive Power (kVar)</td>
<td valign="top" align="center">850.2</td>
<td valign="top" align="center">927.5</td>
<td valign="top" align="center">772.6</td>
<td valign="top" align="center">761.8</td>
<td valign="top" align="center">95.2</td>
</tr>
<tr>
<td valign="top" align="left">Minimum Reactive Power (kVar)</td>
<td valign="top" align="center">&#x02013;570.8</td>
<td valign="top" align="center">&#x02013;830.4</td>
<td valign="top" align="center">&#x02013;590.7</td>
<td valign="top" align="center">&#x02013;720.2</td>
<td valign="top" align="center">&#x02013;670.1</td>
</tr></tbody>
</table>
</table-wrap>
<fig position="float" id="F8">
<label>Figure 8</label>
<caption><p>Distribution of active power losses with varying load positions.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0008.tif">
<alt-text content-type="machine-generated">Line graph illustrating system losses in kilowatts versus metro location in kilometers, with multiple sharp peaks and valleys. System losses range from 0 to 80 kilowatts across the 0 to 18 kilometer span.</alt-text>
</graphic>
</fig>
<p>The loss calculation follows the definitions in <xref ref-type="disp-formula" rid="EQ20">Equations 14</xref>, <xref ref-type="disp-formula" rid="EQ21">15</xref>. At each simulation step <italic>i</italic>, the system loss is obtained as</p>
<disp-formula id="EQ20"><mml:math id="M24"><mml:mtable class="eqnarray" columnalign="center"><mml:mtr><mml:mtd><mml:mtext>&#x00394;</mml:mtext><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>35000</mml:mn><mml:mo>&#x000B7;</mml:mo><mml:msup><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mtext>&#x00394;</mml:mtext><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mi>j</mml:mi><mml:mtext>&#x00394;</mml:mtext><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(14)</label></disp-formula>
<p>where <italic>S</italic><sub><italic>i</italic></sub> is the complex power delivered to the locomotive load, and <italic>I</italic><sup>&#x0002A;</sup> is the conjugate of the source current from the equivalent circuit. The relative loss level is then expressed as</p>
<disp-formula id="EQ21"><mml:math id="M25"><mml:mtable class="eqnarray" columnalign="center"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mo>|</mml:mo><mml:mtext>&#x00394;</mml:mtext><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow><mml:mrow><mml:mo>|</mml:mo><mml:mn>35000</mml:mn><mml:mo>&#x000B7;</mml:mo><mml:msup><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup><mml:mo>|</mml:mo></mml:mrow></mml:mfrac><mml:mo>&#x000D7;</mml:mo><mml:mn>100</mml:mn><mml:mi>%</mml:mi><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(15)</label></disp-formula>
<p>These equations are applied to compute the active loss &#x00394;<italic>P</italic><sub><italic>i</italic></sub>, reactive loss &#x00394;<italic>Q</italic><sub><italic>i</italic></sub>, and total loss ratio &#x003B4;<sub><italic>i</italic></sub> during train operation. The results under different train positions are presented in <xref ref-type="fig" rid="F8">Figures 8</xref>, <xref ref-type="fig" rid="F9">9</xref>.</p>
<fig position="float" id="F9">
<label>Figure 9</label>
<caption><p>Daily active power losses at different voltage levels vs. the number of traction substations.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="frsgr-05-1623347-g0009.tif">
<alt-text content-type="machine-generated">Line graph showing the relationship between the number of traction substations and daily active power loss in kilowatt-hours for six voltage levels ranging from three kilovolts to eleven kilovolts. Power loss decreases as the number of substations increases, with higher voltage levels exhibiting lower power losses throughout.</alt-text>
</graphic>
</fig>
<p>The supply capability of the traction system is mainly constrained by the permissible voltage range of locomotives. By defining the ratios of maximum and minimum line voltages to the nominal level, the acceptable range is set as [0.76, 1.16]. Simulation results show that as voltage level increases, the minimum number of traction substations required decreases. For instance, when the voltage is no less than 5 kV, two substations are sufficient to cover the entire line. Conversely, for a fixed number of substations, the minimum feasible voltage level decreases as more substations are added.</p>
<p>In terms of efficiency, the daily active power losses between the main substation and traction loads decrease with both higher voltage levels and more substations, as illustrated in <xref ref-type="fig" rid="F9">Figure 9</xref>. Note that the calculation in this section is purely mechanism-based. The accuracy of the results depends on the correctness of cable and transformer parameters. In real systems, measurement-based or data-driven models can be used to calibrate uncertain parameters, and thus can serve as a complementary tool to enhance the accuracy of the proposed mechanism-based method.</p>
<p><bold>In summary</bold>, increasing the supply voltage and adding more traction substations are both effective strategies for reducing system losses. Higher voltage levels reduce the required number of substations, lowering construction costs, while fewer substations may increase overall losses. This highlights the essential role of <italic>theoretical line loss calculations</italic>, which directly reflect system mechanisms and provide a reliable foundation for planning and optimization of metro traction networks.</p></sec>
</sec>
<sec sec-type="conclusions" id="s5">
<label>5</label>
<title>Conclusion</title>
<p>This paper investigates a theoretical line loss calculation method for metro traction distribution networks, using actual metro load data as the case study. The research highlights how theoretical line loss models, rather than statistical methods, directly capture the mechanisms of losses and thus provide a solid foundation for system evaluation and planning. The main findings are as follows:</p>
<list list-type="order">
<list-item><p>By analyzing the operational characteristics of the traction network, the line loss calculation formula was derived using the loop current method. This clarified the composition of losses and explained their relationship with dynamic metro loads.</p></list-item>
<list-item><p>Considering the network&#x00027;s structural features, a method was established to allocate the cross terms of line loss to individual loads, enabling more precise loss attribution across substations and circuits.</p></list-item>
</list>
<p>Future work can be expanded in the following directions:</p>
<list list-type="order">
<list-item><p>Develop optimization schemes for targeted loss reduction indicators, thereby providing actionable strategies for system operation and planning.</p></list-item>
<list-item><p>Explore the long-term relationship between line losses and load curve characteristics, building line loss models expressed through curve parameters and comparing them with traditional statistical approaches.</p></list-item>
</list></sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>YiZ: Writing &#x02013; original draft, Writing &#x02013; review &#x00026; editing. ZC: Conceptualization, Writing &#x02013; review &#x00026; editing. ZZ: Writing &#x02013; review &#x00026; editing. YL: Writing &#x02013; review &#x00026; editing. YuZ: Data curation, Writing &#x02013; review &#x00026; editing, Conceptualization. XL: Writing &#x02013; original draft.</p>
</sec>
<ack><title>Acknowledgments</title><p>The authors would like to thank the colleagues from Guangzhou Power Supply Bureau and CSG Electric Power Research Institute for their valuable technical support and insightful discussions during the course of this research.</p>
</ack>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>YiZ, ZC, and ZZ were employed by Guangdong Power Grid Co., Ltd.</p>
<p>The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="ai-statement" id="s9">
<title>Generative AI statement</title>
<p>The author(s) declared that generative AI was not used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p></sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Al-Khalidi</surname> <given-names>H.</given-names></name> <name><surname>Kalam</surname> <given-names>A.</given-names></name></person-group> (<year>2008</year>). <article-title>&#x0201C;Reduction in sub-transmission loss using underground power cables,&#x0201D;</article-title> in <source>2008 Australasian Universities Power Engineering Conference</source> (<publisher-loc>IEEE</publisher-loc>), <fpage>1</fpage>&#x02013;<lpage>6</lpage>.</mixed-citation>
</ref>
<ref id="B2">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cai</surname> <given-names>Y.</given-names></name> <name><surname>Irving</surname> <given-names>M.</given-names></name> <name><surname>Case</surname> <given-names>S.</given-names></name></person-group> (<year>1996</year>). <article-title>Compound matrix partitioning and modification for the solution of branched autotransformer traction feeds</article-title>. <source>IEE Proc. Electr. Power Applic</source>. <volume>143</volume>, <fpage>251</fpage>&#x02013;<lpage>257</lpage>. doi: <pub-id pub-id-type="doi">10.1049/ip-epa:19960375</pub-id></mixed-citation>
</ref>
<ref id="B3">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>M.</given-names></name> <name><surname>Li</surname> <given-names>Q.</given-names></name> <name><surname>Roberts</surname> <given-names>C.</given-names></name> <name><surname>Hillmansen</surname> <given-names>S.</given-names></name> <name><surname>Tricoli</surname> <given-names>P.</given-names></name> <name><surname>Zhao</surname> <given-names>N.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Modelling and performance analysis of advanced combined co-phase traction power supply system in electrified railway</article-title>. <source>IET Gener. Trans. Distr</source>. <volume>10</volume>, <fpage>906</fpage>&#x02013;<lpage>916</lpage>. doi: <pub-id pub-id-type="doi">10.1049/iet-gtd.2015.0513</pub-id></mixed-citation>
</ref>
<ref id="B4">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dolara</surname> <given-names>A.</given-names></name> <name><surname>Gualdoni</surname> <given-names>M.</given-names></name> <name><surname>Leva</surname> <given-names>S.</given-names></name></person-group> (<year>2012</year>). <article-title>Impact of high-voltage primary supply lines in the 2 &#x000D7; 25 kv&#x02013;50 hz railway system on the equivalent impedance at pantograph terminals</article-title>. <source>IEEE Trans. Power Deliv</source>. <volume>27</volume>, <fpage>164</fpage>&#x02013;<lpage>175</lpage>. doi: <pub-id pub-id-type="doi">10.1109/TPWRD.2011.2167635</pub-id></mixed-citation>
</ref>
<ref id="B5">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghassemi</surname> <given-names>M.</given-names></name></person-group> (<year>2019</year>). <article-title>High surge impedance loading (hsil) lines: a review identifying opportunities, challenges, and future research needs</article-title>. <source>IEEE Trans. Power Deliv</source>. <volume>34</volume>, <fpage>1909</fpage>&#x02013;<lpage>1924</lpage>. doi: <pub-id pub-id-type="doi">10.1109/TPWRD.2019.2910210</pub-id></mixed-citation>
</ref>
<ref id="B6">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goodman</surname> <given-names>C.</given-names></name> <name><surname>Kulworawanichpong</surname> <given-names>T.</given-names></name></person-group> (<year>2002</year>). <article-title>Sequential linear power flow solution for ac electric railway power supply systems</article-title>. <source>WIT Trans. Built Environ</source>. <volume>61</volume>:<fpage>10</fpage>.</mixed-citation>
</ref>
<ref id="B7">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gulakhmadov</surname> <given-names>A.</given-names></name> <name><surname>Asanova</surname> <given-names>S.</given-names></name> <name><surname>Asanova</surname> <given-names>D.</given-names></name> <name><surname>Safaraliev</surname> <given-names>M.</given-names></name> <name><surname>Tavlintsev</surname> <given-names>A.</given-names></name> <name><surname>Lyukhanov</surname> <given-names>E.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Power flows and losses calculation in radial networks by representing the network topology in the hierarchical structure form</article-title>. <source>Energies</source> <volume>15</volume>:<fpage>765</fpage>. doi: <pub-id pub-id-type="doi">10.3390/en15030765</pub-id></mixed-citation>
</ref>
<ref id="B8">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>S.</given-names></name> <name><surname>Hemida</surname> <given-names>H.</given-names></name> <name><surname>Yang</surname> <given-names>M.</given-names></name></person-group> (<year>2016</year>). <article-title>&#x0201C;Numerical calculation of the slipstream generated by a crh2 high-speed train,&#x0201D;</article-title> in <source>Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit</source>, 103&#x02013;116. doi: <pub-id pub-id-type="doi">10.1177/0954409714528891</pub-id></mixed-citation>
</ref>
<ref id="B9">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>X.</given-names></name> <name><surname>Wang</surname> <given-names>H.</given-names></name> <name><surname>Li</surname> <given-names>Q.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>Electrical characteristics of new three-phase traction power supply system for rail transit</article-title>. <source>Rail. Eng. Sci.</source> <volume>31</volume>, <fpage>75</fpage>&#x02013;<lpage>88</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s40534-022-00290-1</pub-id></mixed-citation>
</ref>
<ref id="B10">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>M.-K.</given-names></name> <name><surname>Kim</surname> <given-names>M.-S.</given-names></name> <name><surname>Kim</surname> <given-names>D.-H.</given-names></name> <name><surname>Lee</surname> <given-names>J.-W.</given-names></name></person-group> (<year>2010</year>). <article-title>A study on the impedance calculation by using equivalent model in catenary system</article-title>. <source>Int. J. Railway</source> <volume>3</volume>, <fpage>46</fpage>&#x02013;<lpage>53</lpage>.</mixed-citation>
</ref>
<ref id="B11">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kulworawanichpong</surname> <given-names>T.</given-names></name></person-group> (<year>2015</year>). <article-title>Multi-train modeling and simulation integrated with traction power supply solver using simplified newton-raphson method</article-title>. <source>J. Modern Transport</source>. <volume>23</volume>, <fpage>241</fpage>&#x02013;<lpage>251</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s40534-015-0086-y</pub-id></mixed-citation>
</ref>
<ref id="B12">
<mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>F.</given-names></name> <name><surname>Song</surname> <given-names>L.</given-names></name> <name><surname>Cong</surname> <given-names>B.</given-names></name></person-group> (<year>2020</year>). <article-title>&#x0201C;Power line loss analysis and loss reduction measures,&#x0201D;</article-title> in <source>Proceedings of the International Conference on Advanced Intelligent Systems and Informatics 2019</source> (<publisher-loc>Cham</publisher-loc>: <publisher-name>Springer International Publishing</publisher-name>), <fpage>972</fpage>&#x02013;<lpage>980</lpage>. doi: <pub-id pub-id-type="doi">10.1007/978-3-030-31129-2_90</pub-id></mixed-citation>
</ref>
<ref id="B13">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>Q.-z.</given-names></name></person-group> (<year>2016</year>). <article-title>Industrial frequency single-phase ac traction power supply system for urban rail transit and its key technologies</article-title>. <source>Journal of Modern Transp.</source> <volume>24</volume>, <fpage>103</fpage>&#x02013;<lpage>113</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s40534-016-0097-3</pub-id></mixed-citation>
</ref>
<ref id="B14">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ni</surname> <given-names>L.</given-names></name> <name><surname>Yao</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Zhang</surname> <given-names>J.</given-names></name> <name><surname>Yuan</surname> <given-names>J.</given-names></name> <name><surname>Zhou</surname> <given-names>Y.</given-names></name></person-group> (<year>2019</year>). <article-title>&#x0201C;A review of line loss analysis of the low-voltage distribution system,&#x0201D;</article-title> in <source>2019 IEEE 3rd International Conference on Circuits, Systems and Devices (ICCSD)</source>, 111&#x02013;114. doi: <pub-id pub-id-type="doi">10.1109/ICCSD.2019.8843146</pub-id></mixed-citation>
</ref>
<ref id="B15">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shu</surname> <given-names>Z.</given-names></name> <name><surname>Xie</surname> <given-names>S.</given-names></name> <name><surname>Lu</surname> <given-names>K.</given-names></name> <name><surname>Zhao</surname> <given-names>Y.</given-names></name> <name><surname>Nan</surname> <given-names>X.</given-names></name> <name><surname>Qiu</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>Digital detection, control, and distribution system for co-phase traction power supply application</article-title>. <source>IEEE Trans. Ind. Electr</source>. <volume>60</volume>, <fpage>1831</fpage>&#x02013;<lpage>1839</lpage>. doi: <pub-id pub-id-type="doi">10.1109/TIE.2012.2190959</pub-id></mixed-citation>
</ref>
<ref id="B16">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Steczek</surname> <given-names>M.</given-names></name> <name><surname>Jefimowski</surname> <given-names>W.</given-names></name> <name><surname>Maciok</surname> <given-names>T.</given-names></name> <name><surname>Bara&#x00144;ski</surname> <given-names>S.</given-names></name></person-group> (<year>2023</year>). <article-title>Analysis of energy losses in the novel distributed power supply system for trams</article-title>. <source>IEEE Trans. Transport. Electrific</source>. <volume>9</volume>, <fpage>1498</fpage>&#x02013;<lpage>1508</lpage>. doi: <pub-id pub-id-type="doi">10.1109/TTE.2022.3165702</pub-id></mixed-citation>
</ref>
<ref id="B17">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>W.</given-names></name> <name><surname>Xian</surname> <given-names>R.</given-names></name> <name><surname>Liu</surname> <given-names>H.</given-names></name> <name><surname>Xian</surname> <given-names>R.</given-names></name> <name><surname>Zheng</surname> <given-names>X.</given-names></name> <name><surname>Liu</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2025</year>). <article-title>&#x0201C;Impact of outlet short-circuit fault on the electrical-thermal characteristics of distribution transformer windings,&#x0201D;</article-title> in <source>Southern Power System Technology</source>, <fpage>1</fpage>&#x02013;<lpage>9</lpage>.</mixed-citation>
</ref>
<ref id="B18">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Qin</surname> <given-names>S.</given-names></name> <name><surname>Sun</surname> <given-names>J.</given-names></name> <name><surname>Shi</surname> <given-names>X.</given-names></name> <name><surname>Liu</surname> <given-names>T.</given-names></name> <name><surname>Lei</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2025</year>). <article-title>&#x0201C;Research review on transparent technology of low-voltage substation area,&#x0201D;</article-title> in <source>Southern Power System Technology</source>, Southern Power System Technology (&#x05357;&#x065B9;&#x07535;&#x07F51;&#x06280;&#x0672F;) (Guangzhou: Editorial Department of Southern Power System Technology), <fpage>1</fpage>&#x02013;<lpage>12</lpage>.</mixed-citation>
</ref>
</ref-list>
<fn-group>
<fn fn-type="custom" custom-type="edited-by" id="fn0001">
<p>Edited by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1986895/overview">Chengxi Liu</ext-link>, Wuhan University, China</p>
</fn>
<fn fn-type="custom" custom-type="reviewed-by" id="fn0002">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3066439/overview">Lijun Yang</ext-link>, Yanshan University, China</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3128093/overview">Thai Nguyen</ext-link>, Ho Chi Minh City University of Transport, Vietnam</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3130641/overview">Dong Doan Van</ext-link>, Ho Chi Minh City University of Transport, Vietnam</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3211610/overview">Jian Zhang</ext-link>, The University of Tennessee, Knoxville, United States</p>
</fn>
</fn-group>
</back>
</article>