Serviços Personalizados
Journal
Artigo
Indicadores
 Citado por SciELO
 Acessos
Links relacionados
Compartilhar
Sba: Controle & Automação Sociedade Brasileira de Automatica
versão impressa ISSN 01031759
Sba Controle & Automação v.18 n.4 Natal out./dez. 2007
https://doi.org/10.1590/S010317592007000400009
SISTEMAS DE POTÊNCIA
Distribution transformers modeling with angular displacement  actual values and per unit analysis
Dario E. Rodas R.^{I}; Antonio PadilhaFeltrin^{II}; Luis F. Ochoa^{III}
^{I}Universidad Tecnológica de Pereira (UTP)  Pereira, Colombia  Phone: +57 63212044  drodas@utp.edu.co
^{II}UNESP  Campus de Ilha Solteira  Av. Brasil 56, Ilha Solteira  SP, Brazil  Phone: +55 18 3743 1169  padilha@dee.feis.unesp.br
^{III}University of Edinburgh  Edinburgh EH9 3JL, UK  Phone: +44 131 6505575  luis_ochoa@ieee.org
ABSTRACT
In this work a detailed modeling of threephase distribution transformers aimed at complementing wellknown approaches is presented. Thus, incidence of angular displacement and tapping is taken into account in the proposed models, considering both actual values and per unit. The analysis is based on minimal data requirement: solely shortcircuit admittance is needed since threephase transformers are treated as nonmagneticallycoupled singlephase transformers. In order to support the proposed methodology, results obtained through laboratory tests are presented.
Keywords: Transformer modeling, angular displacement, distribution systems.
1 NOMENCLATURE
V_{x}: Voltage corresponding to winding x, x Î [1, 6].
V_{a}, V_{b}, V_{c}: Voltages corresponding to primaryside phases a, b and c.
V_{A}, V_{B}, V_{C}: Voltages corresponding to secondaryside phases A, B and C.
Z_{1}, Z_{2}: Impedances corresponding to windings 1 and 2, respectively.
Z_{m}: Mutual impedance between windings 1 and 2.
I_{1}, I_{2}: Currents corresponding to windings 1 and 2, respectively.
n_{1}, n_{2}: Turns corresponding to windings 1 and 2, respectively.
Yt: Shortcircuit admittance referred to the primary side.
a: Transformation ratio.
[Yp]: Primitive admittance matrix.
S_{n}: Base power.
V_{n}_{1}, V_{n}_{2}: Base voltages corresponding to windings 1 and 2, respectively.
Z_{base}: Base impedance.
R, X: Resistance and reactance, respectively.
a, b: Tap variations at the primary and secondary sides respect to their nominal values, respectively.
[Ybus]: Admittance matrix.
[Y_{pp}], [Y_{ps}], [Y_{sp}], [Y_{ss}]: Submatrices of matrix [Ybus].
[N]: Connection matrix.
[N_{pp}], [N_{ss}]: Submatrices of matrix [N].
g: Grounded.
2 INTRODUCTION
It is well known that transformers are fundamental elements of power systems, requiring proper models so accurate analyses of the network can be performed. In transmission system, transformer modeling is singlephase based since its operation is considered to be balanced. On the other hand, given the inherent unbalance nature of distribution systems, a detailed modeling of distribution transformers requires a threephase approach taking into account a  b  c components, including neutral and ground effects. Thus, threephase transformer modeling has been the focus of several researches (Chen et alli, 1991; Chen and Chang, 1992; Chen et alli, 1996; Chen and Chang, 1996; Baran and Staton, 1997; Kersting et alli, 1999; Kersting, 2002; Dugan, 2003; Wang et alli, 2004; Xiao et alli, 2006). Nonetheless, models found in the literature still lack of details regarding some connections, angular displacement and tapping.
Depending on the power flow algoithm adopted for the distribution network analysis, transformer models may lead to convergence problems. Consequently, some proposals are given in the literature in order to minimize this effect (Chen and Chang, 1996; Baran and Staton, 1997; Kersting et alli, 1999; Wang et alli, 2004; Xiao et alli, 2006). Evaluation of the impacts that certain transformer models present on the power flow analysis is out of the scope of this paper.
Although transformer models presented in 1991 by Chen et alli have been widely used by researchers and commercial software, their applicability outside the US is limited since the angular displacement is commonly restricted to two clock hours: 1 and 11, depending on whether it is a stepup or stepdown arrangement.
The objective of this paper is to improve and to make as general as possible widelyaccepted transformer models, considering angular displacement as well as tapping in both primary and secondary windings. It is also demonstrated that distribution transformers can be modeled as power transformers when analyzing variable taps. Proposed threephase transformer models allow including angular displacement, which is a basic tool for generalizing transformer models to other different clock hours or even vector groups. A discussion about circuit and matrix modeling considering the variables in actual values and in per unit is also included. Models in actual values (Siemens) are fully described since they allow obtaining directly all major parameters involved in the transformer modeling, leading to a better understanding of the per unit analysis. Moreover, in order to validate the proposed methodology, results from laboratory tests are presented.
Generally speaking, there two ways to model distribution transformer. One requires a large quantity of data (Gorman and Grainger, 1992), which is often difficult to obtain. The other way considers the use of minimal information (only the shortcircuit admittance is needed) whereas the core is modeled as a load (Dugan, 2003; Chen et alli, 1991). The modeling presented in this work is classified in the latter group. Moreover, this proposal allows the modeling of all kinds of winding connections given by the American National Standards Institute (ANSI) and the International Electrotechnical Commission (IEC).
The structure of this paper is as follows: section 3 presents the fundamental concepts for singlephase transformers, whereas in section 4 the incidence of taps on singlephase transformers is addressed. Distribution transformer analysis is presented in section 5. Finally, in section 6 transformer bank modeling is addressed taking into account actual values and per unit, as well as angular displacement and tapping incidence.
3 SINGLEPHASE TRANSFORMER FUNDAMENTAL CONCEPTS
A singlephase transformer model based on matrices and, initially, considering its nominal values (winding 1 with n_{1} turns and winding 2 with n_{2} turns), as presented in figure 1, is given by:
These values are winding based, i.e. primitive impedance is used.
Equations (1) and (2) could be either based on actual values or p.u. Regarding the case of actual values in equation (2), Yt is the shortcircuit admittance referred to the primary. Then, the primitive admittance matrix is:
where a = .
Matrix [Yp] in p.u. is obtained by dividing each term of equation (3), whose units are Siemens, by its corresponding base values:
where p: primary and s: secondary.
Then, substituting a = , the primitive admittance [Yp] in p.u. is obtained:
4 GENERAL MATHEMATICAL MODEL FOR A SINGLEPHASE TRANSFORMER WITH VARIABLE TAPS
When taps vary, fundamentally two parameters experience alterations:
a) The voltage at the secondary side of the transformer due to the new transformation ratio, and;
b) The transformer shortcircuit impedance since the inductive reactance is a function of the squared number of turns. When modifying the number of turns, the resistance is also affected.
When nominal values are considered (taps at their nominal positions in both primary and secondary sides), shortcircuit tests are performed in order to obtain the transformer's nominal admittance: Yt_{nom}.
4.1 Offnominal Taps incidence on Transformer's Impedance
In power transformers, where the resistance can be neglected (Z » X), tap variations modify the impedances according to the new number of turns and quadratically since the inductance is a function of the squared number of turns.
Let a and b be the tap variations at the primary and secondary sides respect to their nominal values, respectively. Thus, when taps are altered, the new transformer ratio is:
Consequently, Z_{1  new} = a^{2}Z_{1  nom}, Z_{2  new} = b^{2}Z_{2  nom} and Z_{m  new} = abZ_{m  nom}.
This is valid for both actual values and p.u. as long as the nominal base values of Z_{base} are kept constant. The transformer's new admittance will be:
Then, using equation (2) and considering nominal values:
4.2 Offnominal Taps incidence on Transformer's Primitive Admittance
In order to avoid mistakes in the analysis, the process performed in section 3 is repeated taking into account a new number of turns for both primary and secondary sides. Thus,
When the shortcircuit test is carried out, the voltage in the secondary side is null ( = 0). Then, according to equation (5), the transformer's shortcircuit admittance is:
In addition, the following expressions are also obtained:
Finally, it is obtained:
Then, substituting Yt' according to equation (5):
Equation (13) is in Siemens and Yt_{nom} stands for the nominal shortcircuit admittance referred to the primary side. For analysis in p.u., the following expression should be used:
5 DISTRIBUTION TRANSFORMERS
In this section a reallife distribution transformer (20 KVA, 13200 V/240 V) is analyzed by means of laboratory tests. Taps at the primary side are:
Tap1 = 13530 V (2.5%);
Tap2 = 13200 V (nominal);
Tap3 = 12870 V (2.5%);
Tap4 = 12540 V (5.0%);
Tap5 = 12210 V (7.5%).
The following analysis will consider solely variations at the primary side, which is the case of distribution network transformers. Consequently, b = 1. These shortcircuit tests were carried out considering the transformer's nominal current.
5.1 Laboratory Tests
Tests carried out varying the available taps of the transformer led to the results shown in table 1.
Considering the impedance (Z_{t} = R_{t} + jX_{t}) of the most distant tap from the nominal position (Tap 5, a = 0.925), one can compute the corresponding shortcircuit admittance referred to the primary side. Thus,
= 90.34 + j231.742W, then
= ()^{1} = 1.46 · 10^{3}  j3.746 · 10^{3} Siemens.
Taking into account that in this case a = 13200 / 240 = 55 and using equation (12), the transformer's primitive admittance matrix is obtained:
In order to verify the performance of the power transformer model applied on distribution transformers, equation (13) is used to compute the primitive admittance matrix for Tap 5 (table 1). The nominal shortcircuit admittance is calculated as follows:
Zt_{nom} = 101.61 + j273.287W, then
Yt_{nom} = (Zt_{nom})^{1} = 3.4298 · 10^{3 }Ð  69.6^{°} Siemens.
Hence, considering also a = 55 and a = 0.925, the following matrix is calculated:
The comparison of both results prove the good performance of the method based on power transformers since the errors in magnitude do not exceed 0.3%, whereas for the angle this value is 1.3%. Consequently, given the accuracy obtained with a transformer of relatively low capacity (20 KVA), it is expected that most distribution transformers (commonly with capacities higher than 20KVA) present better results. Offnominal load cenarios are also expected to exhibit similar performances.
5.2 Modeling with Taps at the Primary Side
When taps are located at the primary side, the winding at the secondary side (winding 2) is not modified. Therefore, the resistance referred to the secondary side (R_{2}) is not modified either. Nevertheless, when taps vary, R_{2} referred to the primary side should take into account the transformation ratio aa. In addition, considering that the resistance of the primary side varies according to the number of turns, in other words, the resistance varies linearly with a and » R_{1} , the new resistance is obtained as follows:
On the other hand, the reactance varies quadratically:
Thus, the distribution transformer model for any a is:
Another way to approach the transformer modeling with taps at the primary side is that when distribution transformers are considered as power transformers. In this case, selfquadratic variation of the reactance (a^{2}) is applied to the whole impedance, i.e. it is applied on the resistance as well.
In order to compare both approaches presented above, and given that the resistance is the parameter being considered differently, the multiplying factors of R_{t  nom} are presented in table 2:
Since results are similar, approximations can be used for the resistance. Moreover, generally speaking, the whole impedance can be approximated. Thus, Y_{new} = Y_{nom}/a^{2}.
Results obtained through laboratory measurements (table 1) and through equations (15), (16), and the approximation above presented, are shown in table 3.
Strictly speaking, the reactive part of the model adjusts properly the estimated values compared to those measured in laboratory. It is also observed that, regarding the resistive part, the first approach, equation (15), follows more accurately (almost perfect adjustment) the measured values. Nevertheless, due to systematic problems when forming matrix [Ybus], it is recommendable to consider the distribution transformer as power transformer (a^{2}) since embedded errors are relatively small.
6 TRANSFORMER BANKS MODELING
This section presents in a general manner the mathematical models of threephase transformers obtained from using singlephase transformers as a threephase bank. Incidence of the transformer's angular displacement on its corresponding model is also addressed.
Magnetization impedance of a transformer is quite large and is neglected for voltage regulation studies. This approximation does not affect power flow results. Iron losses obtained through the opencircuit test can be modeled by adding them to the actual load.
Modeling a threephase transformer considering all magnetic couplings may require data from several shortcircuit tests. Therefore, it is a common procedure to consider threephase distribution transformers as a bank of three singlephase transformers, magnetically isolated from each other and connected properly to achieve the connection modeling under analysis (Kersting, 2002). Thus, required information for modeling the threephase transformer is reduced exclusively to the shortcircuit test for solely one transformer, information that is always available.
In general, it is recommended to obtain matrix [Yp] directly from the shortcircuit tests since, otherwise, in order to avoid numerical instability problems it would be necessary to obtain extremely accurate measurements of the Z parameters.
The following subsections 6.1 and 6.2 are aimed at presenting the procedures needed to obtain matrix [Ybus] for both approaches, Siemens and per unit. Two application examples will be presented thereafter. Subsections 6.4 and 6.5 present the connections summary for various types of connections and the transformer modeling considering the ideal ordering for the primitive admittance and connection matrices, respectively. The latter is important since conclusions will be drawn in order to model the incidence of the angular displacement, presented in subsection 6.6. Finally, the incidence of taps in the threephase transformer modeling is addressed in subsection 6.7.
6.1 Approach in Siemens
To obtain matrix [Ybus] in actual values the following steps should be performed:
1. A 6 by 6 primitive admittance matrix [Yp] is formed in Siemens, composed by three decoupled singlephase transformers. First, third and fifth rows and columns represent the primary sides of transformers, whereas the remaining rows and columns represent the corresponding secondary sides.
This step is general for all possible connections in Siemens and, meanwhile, considering nominal taps. At this stage the threephase transformer is not yet connected.
2. For a particular connection under study the corresponding connection matrix [N], which relates windings voltages to nodal voltages and considers the angular displacement, should be selected.
3. The following formula is applied:
6.2 Approach in Per Unit
Matrix [Ybus] can be also computed in per unit following a similar procedure to that presented for actual values:
1. Same as step 1 of previous subsection, but considering matrix [Yp] in p.u., using three submatrices as in equation (4).
2. Same as step 2 of previous subsection.
3. Same as step 3 of previous subsection.
4. Matrix [Ybus] can be presented in terms of its submatrices:
where [Y_{pp}] stands for the primary side, [Y_{ss}] stands for the secondary side and [Y_{ps}] and [Y_{sp}] for the primarysecondary and secondaryprimary, respectively. These submatrices should be corrected due to changes in base voltages according to the corresponding connection, as follows:
Wyewye: [Ybus] remains without alteration.
Wyedelta: Base voltage at the secondary changes, consequently [Ybus] in p.u. should be adapted in terms of its submatrices. Thus, [Y_{pp}] remains the same, [Y_{ps}] and [Y_{sp}] are divided by , and [Y_{ss}] is divided by 3.
Deltadelta: Both base voltages at the primary and secondary sides change, therefore the whole matrix [Ybus] should be divided by 3.
Deltawye: Base voltage at the primary side changes. Thus, [Y_{pp}] is divided by 3, [Y_{ps}] and [Y_{sp}] are divided by , and [Y_{ss}] remains the same.
6.3 Examples
In this subsection two distribution transformers, with grounded and isolated neutral, will be analyzed aimed at illustrating the calculation details of matrix [Ybus] considering both approaches, actual values and per unit.
6.3.1 Groundedneutral Yd1 Transformer
The electrical scheme of the wyedelta transformer with angular displacement of 30º (1 o'clock) is presented in figure 2.
Firstly, winding voltages are presented as a function of node voltages. Numbers are referred to winding voltages (1, 3 and 5 for the primary side and 2, 4 and 6 for the secondary side) and letters (a, b, c, A, B, C) are referred to nodal voltages respect to the reference (ground), as follows:
The matrix which is multiplying the nodal voltages is called connection matrix [N]. Then, solving equation (19) for this particular connection, the admittance matrix for those six nodes is expressed by:
At this point, the problem has been solved in Siemens, considering a equal to the nominal winding ratio and Y_{t} as the shortcircuit admittance (in Siemens) referred to the primary side.
If values in p.u. are required, firstly, equation (19) should consider a = 1 and Y_{t} in p.u. Secondly, due to the wyedelta connection, submatrices [Y_{ps}] and [Y_{sp}] are divided by , and [Y_{ss}] is divided by 3. Thus, it is calculated:
6.3.2 Isolatedneutral Yd1 Transformer
In this case there exist 7 nodes: a, b and c at the primary side; A, B and C at the secondary side and the floating point n (common point of wye connection).
Then, considering V_{1} = V_{a}  V_{n} , V_{3} = V_{b}  V_{n} and V_{5} = V_{c}  V_{n} , a 7 by 7 admittance matrix in p.u. can be computed by using equation (19). Finally, using Kron's reduction for eliminating the floating point n, the following matrix is obtained:
6.4 Connections Summary
6.4.1 In Siemens
The connections summary based on the four submatrices of [Ybus], i.e. [Y_{pp}], [Y_{ps}], [Y_{sp}] and [Y_{ss}], is shown in table 4. This is actually an extension of a table presented in Chen et alli (1991) for per unit analysis and including the corresponding angular displacement.
Matrices presented in table 4, considering Y_{t} in Siemens, are given by:
and
Notice that in general [Y_{sp}] = [Y_{ps}]^{t}. Given that [Y_{I}] = [Y_{I}]^{t} and [Y_{II}] = [Y_{II}]^{t} in some cases [Y_{sp}] = [Y_{ps}].
6.4.2 In Per Unit
When an analysis in per unit is required, table 4 can be reduced, obtaining table 5 as presented in Chen et alli (1991).
Matrices [Y_{I}], [Y_{II}] and [Y_{III}] are those of the previous approach (Siemens) while considering Y_{t} in p.u.
Table 5 can be extended by including the angular displacement, as shown in table 4. It is important to remark that for those transformers fabricated under the ANSI 30 degrees, angular displacement for stepup Deltawye or Wyedelta connections is 330º, whereas for stepdown Deltawye or Wyedelta connection is 30º. Consequently, stepdown Deltawye and stepup Wyedelta transformers that do not match the base angular displacement should be corrected as described thereafter in subsection 6.7.
6.5 Ideal ordering for the primitive admittance and connection matrices
For threephase transformers the primitive admittance matrix [Yp], as presented in equation (18), is ordered by grouping firstly the primary side parameters and then the secondary side. Thus, the new order of the rows and columns will be 135246:
For any connection, it can be noticed that the four submatrices in equation (25) are diagonal:
[Y_{pp}]_{prim} = Yt[I],
[Y_{ps}]_{prim} = aYt[I],
[Y_{sp}]_{prim} = aYt[I]
and
[Y_{ss}]_{prim} = a^{2}Yt[I].
In addition, given the diagonal nature of [Yp], it is verified that [Y_{ps}]_{prim} = [Y_{sp}]_{prim} .
Thus,
Analysis in per unit requires a = 1 and Yt in p.u., otherwise a is the windings ratio and Yt the admittance referred to the primary side (in Siemens).
This new reordering makes the connection matrix [N] to present other particularities:
Submatrices [0] are due to the absence of connections between primary and secondary sides in threephase transformers. However, this would not be the case for autotransformers.
In order to compute [Ybus], equation (19) is solved by using equations (25) and (26):
Equation (27) is valid in Siemens. Per unit analysis should consider a = 1, Yt in p.u., and some submatrices should be modified according to subsection 6.2.
Given that, from matrices theory, (A · B)^{t} = B^{t} · A^{t} then:
([N_{pp}]^{t} [N_{ss}])^{t} = [N_{ss}]^{t} ([N_{pp}]^{t})^{t} = [N_{ss}]^{t} [N_{pp}].
Consequently, it is observed that for [Ybus] always [Y_{sp}] = [Y_{sp}]^{t}.
6.6 Incidence of the Angular Displacement
There are two ways for including the angular displacement in the transformer modeling:
a. A connection matrix [N] is constructed taking into account the particular angular displacement under study. Then (19) is solved for either p.u. or siemens.
b. First, a connection summary is used (table 4 or table 5). Then, by considering the corresponding basic clock hour (angular displacement) and the appropriate methodology, any other clock hour can be achieved.
When equation (19) is used in order to compute [Ybus] for each of the 26 IEC connections (excluding 10 zigzag connections) it is observed that submatrices [Y_{pp}] and [Y_{ss}] do not suffer any modification due to changes in angular displacement. Consequently, given that [Y_{sp}] = [Y_{ps}]^{t} it is necessary to solely focus on analyzing [Y_{ps}].
Starting from the basic clock hour, in order to "move" a given clock hour to the next one within the same vector group (see table 6), rows of the submatrix [Y_{ps}] should be "rotated" in such a way that the first row becomes the second, the second becomes the third and the third becomes the first. If the desired clock hour was not achieved yet, the process is repeated considering the new matrix (obtained in the previous step) as the basic one.
If a vector group change is required, firstly the [Y_{ps}] should be obtained considering a connection 6 clock hours above the basic clock hour, which is the negative of it. In this way:

For Yy or Dd connection, whose basic clock hour is 0 (table 4), the corresponding 6clockhoursabove submatrix is: [Y_{ps}]_{6 o'clock }= [Y_{ps}]_{0 o'clock}.

For Yd connection, whose basic clock hour is 1: [Y_{ps}]_{7 o'clock} = [Y_{ps}]_{1 o'clock}.

For Dy connection, whose basic clock hour is 11: [Y_{ps}]_{5 o'clock} = [Y_{ps}]_{11 o'clock}.
If after the vector group change another clock hour (within the same vector group) is required, the procedure explained previously should be performed.
6.7 Incidence of taps in the threephase transformer modeling
When considering offnominal taps, matrix [Ybus] from equation (20) is modified into:
This can be demonstrated by extending the procedure presented in subsection 4.2, for singlephase transformer, to the threephase transformer analysis. Thus, the corresponding primitive admittance matrix is:
Reordering rows and columns based on primary and secondary sides, the new primitive admittance matrix is:
Given that the four primitive submatrices [Y_{pp}]_{prim}, [Y_{ps}]_{prim}, [Y_{sp}]_{prim} and [Y_{ss}]_{prim} are diagonal:
By using equation (26) and solving equation (19), it is obtained:
Thus, submatrix ëY_{pp}û is divided by a^{2}, submatrices ëY_{ps}û and ëY_{sp}û are divided by ab and submatrix [Y_{ss}] is divided by b^{2}. This is valid for all clock hours, both in per unit or Siemens.
7 CONCLUSIONS
A methodology that allows, with minimal data (one shortcircuit test), modeling threephase distribution and power transformers based on a singlephase transformers bank was presented. Proposed models provide detailed information required for implementing most distribution transformer connections. In addition, a methodology to include angular displacement (basic tool for generalizing transformermodels to other different clock hours or even vector groups) and tapping incidence, commonly not addressed in the specialized literature (Wang et alli, 2004; Xiao et alli, 2006), was presented. Both approaches, actual values and per unit, were also addressed aimed at specifying the corresponding differences.
REFERENCES
Baran, M. E. and E. A. Staton. (1997). Distribution transformer models for branch current based feeder analysis. IEEE Transactions on Power Systems, Vol. 12, No. 2, pp. 698703. [ Links ]
Chen, T.H.; Chen, M.S.; Inoue, T.; Kotas, P. and Chebli, E. A. (1991). Three–phase cogenerator and transformer models for distribution system analysis. IEEE Transactions on Power Delivery, Vol. 6, No. 4, pp. 16711681. [ Links ]
Chen, T.H. and Chang, J.D. (1992). Open wyeopen delta and open deltaopen delta transformer modelsfor rigorous distribution system analysis. IEE Proceedings–Generation, Transmission and Distribution, Vol. 139, No. 3, pp. 227234. [ Links ]
Chen, T.H.; Chang, J.D. and Chang, Y.L. (1996). Models of grounded midtap openwye and opendelta connectedtransformers for rigorous analysis of a distribution system. IEE Proceedings –Generation, Transmission and Distribution, Vol. 143, No. 1, pp. 13502360, Jan. 1996. [ Links ]
Chen, T.H. and Chang, Y.L. (1996). Integrated models of distribution transformers and their loads forthreephase power flow analyses. IEEE Transactions on Power Delivery, Vol. 11, No. 1, pp. 507513. [ Links ]
Dugan, R. C. (2003). A perspective on transformer modeling for distribution system analysis. IEEE Power Engineering Society General Meeting, pp. 114119. [ Links ]
Gorman, M. J. and Grainger, J. J. (1992). Transformer Modelling for Distribution System Studies, Part I: Linear Modelling Basics. Part II: Addition ofModels to Ybus and Zbus. IEEE Transactions on Power Delivery, Vol. 7, No. 2, pp. 575580. [ Links ]
Kersting, W. H.; Phillips, W. H. and Carr, W. (1999). A new approach to modeling threephase transformer connections. IEEE Transactions on Industry Applications, Vol. 35, No. 1, pp. 169175. [ Links ]
Kersting, W. H. (2002). Distribution System Modeling and Analysis. CRC Press LLC, Florida. [ Links ]
Wang, Z.; Chen, F.; and Li, J. (2004). Implementing transformer nodal admittance matrices into backward/forward sweepbased power flow analysis for unbalanced radial distribution systems. IEEE Transactions on Power Systems, Vol. 19, No. 4, pp. 18311836. [ Links ]
Xiao, P.; Yu, D. C.; and Yan, W. (2006). A unified threephase transformer model for distribution load flow calculations. IEEE Transactions on Power Systems, Vol. 21, No. 1, pp. 153159. [ Links ]
Artigo submetido em 04/05/2006
1a. Revisão em 22/08/2006
2a. Revisão em 15/05/2007
Aceito sob recomendação do Editor Associado Prof. Denis Vinicius Coury