Abstract
Impedance mismatching between different phases of a multiphase transformer is generally observed e.g., in a three-phase to seven-phase transformer, due to an unequal number of turns in different coils. This mismatching introduces error in the study of per phase equivalent circuit diagrams as well as induces an imbalance in output voltages and currents. Therefore, it is a challenging task to develop a per-phase equivalent circuit for the secondary and primary sides (In some cases) too. This paper proposes an artificial intelligence optimization technique like PSO based modeling of the per-phase equivalent circuit of the secondary side. This paper deals with the modeling and simulation of a three-phase to seven-phase power transformer using Artificial Intelligence technique like particle swarm optimization (PSO) and Genetic Algorithm (GA). The proposed model is optimized through PSO and GA algorithms and tested for minimum voltage error in each phase. The proposed model is designed and the objective function is optimized by PSO & GA in MATLAB environment. It is found that the optimized model can be effectively implemented as a per-phase equivalent circuit for the secondary side.
Introduction
The literature of the Multiphase system was found since 1960 s and flourished after 1960 [12]. Multiphase systems are known for reliability, fault tolerance, and high power carrying capabilities. Literatures are found for its application in transmission and utilization [3, 19]. The various benefits offered by the multiphase system have pulled researchers’ interests are higher machine torque density, lower torque ripples, lower harmonics, better transient and steady-state performances and robust control of multiphase electrical drives [3, 13]. Multiphase AC-DC rectifiers produce lower ripples and higher ripple frequency. Therefore Multiphase AC-DC rectifiers produce DC of higher quality.
Multiphase electric power from a three-phase power is generated by three techniques. One is by converting three-phase AC to a DC power and then converting it back to multiphase AC (i.e. AC-DC-AC); the second is through multiphase power converting transformers; the third one is by employing matrix converter (AC-AC). The second method utilizes transformer winding in a special connection and turns ratios while the first and third method needs rectifiers and inverters [1–4, 22–24].
The literature of 6- and 12-phase transformers based multi-pulse rectifier system is found in [14, 16]. Recently, 24-phase and 36-phase transformer systems have been designed for supplying a multi-pulse rectifier system [2, 26]. The reason of choice for a 6-, 12-, or 24-phase system is that these numbers are multiples of three, and designing this type of system is simple and straight forward. The explanation of the decision for a 6-, 12-, or 24-phase transformer is that these numbers are products of three and the design is straightforward [19].
As far as Multiphase transformer of odd numbers of phases is concerned, limited pieces of literature are found for the odd number of phases viz. 5-, 7- and 11-phase [1, 19–21]. Prototypes were designed and tested for a three to five-phase and three to seven-phase power converting transformer and reported in the literature. One of which is Scott based three to five-phase converting power transformer design is reported in [1]. The approach of [1] is extended and reported in [24] for a three-phase to seven-phase transformer. The work reported in [24] requires further research. This paper deals with the mathematical modeling of the transformer reported in [24].
Per-phase equivalent circuit is vital for the analysis of a multiphase transformer. Although three-phase to seven-phase conversion techniques has been already reported [18, 24] but the research that addresses the secondary winding mismatching, per-phase equivalent circuit modeling is not available yet. This work is the first of its kind that addresses the mismatching issue and developed a per-phase equivalent circuit of the secondary side using intelligent techniques like GA and PSO.
Section 2 deals with the transformer design. Section 3 is problem formulation while section 4 is optimization results. The conclusion is being put in section 5.
Materials and Methods
General Theory of two-phase to n-phase transformation
A more simplified and graphical two-phase to n phase theory of transformation is given in [20]. It is shown in [20] that a three-phase AC can be converted to any number of phases by using a special transformer design given and turn ratios of primary and secondary can be calculated from the generalized equations.
Refer Fig. 1 taken from [20]; let’s assume voltage V α and V β be the output of three-phase to two-phase conversion out of primary Scott connection. Vα and Vβ are voltage induced in winding “RO” and “YB” and have equal RMS value as that of the supply phase but they are 90 degrees phase apart as both constitute a two-phase AC power.

Phasor diagram for 2 phase to n phase conversion.
Let’s assume a phase conversion of two-phase to n balanced phase (refer Fig. 3). Let V
i
be the i
th
phasor of a balanced n phases balanced system at an angle of iδ and have a magnitude of V where
Similarly,
Where i may be any phase of balanced n-phase system i.e., i = 0,1,2 ... ..(n-1),
V i α ,-Projection of ith phasor along x-axis/
V i β -Projection of ith phasor along y-axis/
a α - Turns ratio for ith phase along x-axis or with respect to winding wound on Core 1.
a β -Turns ratio for ith phase along y-axis or with respect to winding wound on Core 2 which produces a flux 90 degrees out of phase w.r.t core 1.
Hence an n phase transformation from a two-phase balanced supply could be achieved by employing two core transformer with a pulsating flux of 90 degrees apart in phase. For n-number of phases, we can get turn ratios of individual windings by varying i from 0 to n-1 and putting
In matrix form, the voltage relationship can be generalized as:
Whereas Turn ratios can be generalized as below
A generalized theory of three-phase to multiphase transformation is developed using the phasor diagram approach in [24]. Using that approach a three-phase to seven-phase transformer can be designed for 1:1 turns rations if windings are connected according to the schematic diagram shown in Fig. 2. Required Turns ratios are shown in Table 1. Proper connections of the multiple windings of the requisite number of turns are the key to generate a seven-phase current from a three-current. The input and output currents are shown in Figs. 3, 4.

Winding arrangement of the transformer.
Turns ratios of each windings

Three-phase input current.

Seven-phase output current.
Per phase equivalent circuit model or any polyphase machine is required for transient as well as steady-state analysis of the machine. It is very clear from Table 1 taken from [24] that each phase is composed of two winding except phase 1. It is clear from Fig. 1 that the windings of two limbs are connected in series with proper polarity and turns ratios so resistance and inductance of each phase shall include both the windings. Resistance and inductance of a coil is calculated by the following expression.
N = Number of turns of the coil
μ=Permeability of the core
A = Area of the coil
l = Average length of the coil
Here A is fixed to the core cross-sectional are for all the windings but l is proportional to the number of turns. So it can be said that
If Resistances and self-inductance of the reference phase is assumed to be R2a and X2a then resistances and inductances of the rest of the phases can be expressed in terms of the reference phase as below.
Equivalent circuit and parameters determination
Per phase equivalent circuit of secondary winding poses some issues as each phase has a different number of windings. An intelligent system is applied to estimate the approximate value of secondary winding parameters for an acceptable error.
The per-phase equivalent circuit model of the secondary side of the multiphase transformer is shown in Fig. 5. Let’s assume that EMF induced in all seven secondary windings is E2, Load current of I L flowing through each circuit. The terminal voltage of each secondary windings will be different as each winding has a different value of impedance.

Exact equivalent circuit of secondary windings.
The aim is to formulate a per-phase equivalent circuit of primary side and secondary side as it has been done for three-phase transformers. The process would be simple and straightforward if all the seven-phases has the same parameters. Unfortunately due to the very nature of transformer design, all the seven-phase secondary as well primary three-phase too has different parameters. This paper is indeed an extension of [24] and aimed to develop a per-phase equivalent circuit modeling.
It is well known that in Scott connection the three-phase primary windings are connected in such a way to produce a two-phase flux separated by 90 degrees in the time domain. That is windings are in the ration of 1:0.57:0.57 and connected as shown in Fig. 1 primary section.. That too may impose imbalance but can be approximated as a two out of three-phase has equal parameters.
Secondary side per phase circuit
Let’s assume that Per phase equivalent circuit is shown in Fig. 5 while exact seven-phase circuit is shown in Fig. 6, where k is the optimum value of the factor which will decide the parameters of per phase equivalent circuit. It will deduced by suitable optimization technique.

Per-phase equivalent circuit of secondary windings.
Output voltages of all seven-phases under steady-state condition can be deduced from the circuit diagram shown in Fig. 5; it is as shown by Equation 9.
The output voltage of per-phase equivalent secondary circuit is
Let us formulate the error voltage (difference of secondary voltages) as objective function,
Similarly
It is clear from Equation 10 that f1 (k), f2 (k), f3 (k), and f4 (k) are the only objective functions to be optimized for this multi-objective optimization problem. This problem is solved by making a single objective function and then using particle swarm optimization (PSO) algorithm to find the best value of variable, k. The whole process is discussed in the section 3.3.
The aim of the optimization problem is to get the best value of k for which the functions (given in Equation 10) simultaneously have the least values. For this, a common objective function (or fitness function) is defined and is written as:
The terms for functions f5(k), f6(k) & f7(k) will have the same effect as that due to functionsf4(k), f3(k) & f1(k) respectively, so these terms have not been used in the objective function.
The optimization problem is a minimization problem such that ofun has the minimum value for a particular value of k.
There are various artificial intelligence (AI) techniques like Artificial Neural Network (ANN), adaptive neuro-fuzzy inference system (ANFIS), Genetic Algorithm (GA), Particle swarm optimization (PSO), etc [6]. Particle swarm optimization (PSO) is one of the best evolutionary techniques that is being used rigorously used by researchers throughout the world [8, 21]. Unlike GA, PSO does not use genetic operators like crossover and mutation. Particles get updated using internal velocity; and have a memory vital to the algorithm. In addition, PSO is a single direction information sharing method, as here only the ‘best’ particle gives out the information to others, and thus the process looks for the best solution. Also as compared to GA, PSO is easy to implement with few parameters to be adjusted.
In order to get the best value of k, PSO algorithm is used in this work. Particle swarm optimization (PSO) is an artificial intelligence (AI) technique modeled as per the performance of bird flocks [14, 27]. It is useful to find the keys to exceptionally tough optimization problems. In PSO algorithm, a swarm of particles is maintained, and each particle’s position represents a solution [10]. The particles change their position by getting influenced by the neighboring best particle (P
best
) as well as the overall best particle (G
best
). The position (
Where, ϕ i , is the step size or the velocity of the particle. Its value is obtained by
Where, ω represents inertia weight, acceleration coefficients are c1 and c2 ; r1 and r2 are the random numbers between 0 & 1; ωmax, ωmin are the maximum and minimum values of inertia weight, and, maxite is the predefined maximum number of iterations. The typical movement of particles during the optimization process is shown in Fig. 3.
In the problem discussed in this paper, the value of k, itself represents the position of a particle and thus Equation 4 can be rewritten as
In the case of GA, which was formerly presented by Holland [28] as an EC method, the optimization process depends on three important operators (namely the selection, crossover, and mutation) on every generation. Here, a number of chromosomes (the candidate solutions) are assessed firstly. Then a fitness function is used to evaluate each chromosome. Higher fitness valued chromosomes are most likely selected and included in the next population of generation and then they go for recombination. Two parents are combined through crossover in order to yield ith offspring using Equations 16 and 17.
where r is a random number, r ∈ [0,1].
Mutation is used to get faster convergence. It preserves genetic variety from one age group to the next and efforts to realize some stochastic unevenness of GA [29, 30].
The application of PSO for finding the best value of k using the defined fitness function (ofun) has been represented in the flowchart in Fig. 8. Similarly, GA was applied to find the best value of k using the defined fitness function (ofun).

Movement of particles during the optimization process.

Flowchart of PSO applied for finding the best value of k.
The values of parameters selected for PSO and GA are given in Tables 2 and 3 respectively.
Values of PSO parameters
Simulation parameters of PSO
Figure 9 shows the 3D representation of the results obtained for the PSO convergence characteristic using MATLAB. It can be observed that the value of fitness function converges to its minimum value in before the 50th iteration. Figure 10 shows the GA results.

PSO convergence characteristic.

GA convergence characteristic and result.
Figure 10 shows the variation of the value of objective function with respect to the value of k during the PSO process. Similarly, Figs. 11–14, show the variation of the values of functions (f1, f2, f3 and f4 respectively) with respect to the value of k during the PSO process.
From Fig. 15, it can be observed that the value of the fitness function is minimum near k = 1.23. The best-optimized value of k obtained using PSO is found to be 1.2364.

Variation of function f1 w.r.t. k.

Variation of function f1 w.r.t. k.

Variation of function f3 w.r.t. k.

Variation of function f4 w.r.t. k.

Variation of objective function w.r.t k during the optimization process.
This result can be verified by comparing it with the result obtained by plotting the ofun w.r.t k and observing its minimum value, the value of k corresponding to the minimum value of ofun is found to be 1.236 & that with GA is 1.2363 (Table 4). This value is very close to that obtained by PSO, but since PSO gives a more précis value so it can be said that PSO is a much better technique to obtain the value of k.
Comparison of results
The main dissimilarity between the PSO tactics as compared with EC like GA is that PSO does not use genetic operators like crossover and mutation. Particles get updated using internal velocity; and have a memory vital to the algorithm. Also, PSO is a single direction information sharing method, as here only the ‘best’ particle gives out the information to others, and thus the process looks for the best solution. Also as compared to GA, PSO is easy to implement with few parameters to be adjusted.
There per phase equivalent circuit of the secondary side of the transformer can be expresses as shown in Fig. 12.
The developed model can be verified by estimating output voltage under-compensation as well as without compensation. Let’s assume a load current of 5A(300Lag) is flowing through the secondary, an EMF of 150 V is induced in the secondary winding. Let’s assume Reference winding resistance equals to 1 ohms and reactance as per norms(X/R = 7). Results are tabulated in Table 4.
The R & L parameters are changed and the data is generated similar to that presented in Table 5. The standard deviation of compensated and non-compensated is calculated and plotted in Fig 17. It can be seen that compensated data is changed less as compared to non-compensated data in response to parameter variation.

Per phase equivalent circuit of the secondary side of the transformer.

Standard deviation of the seven-phase output voltage with parameter variation.
Output voltage under compensation and no compensation
Artificial intelligence-based algorithms like Particle swarm optimization techniques like (PSO) and GA can be used effectively to model the multiphase transformer circuits that pose unbalancing as a challenge. The obtained result can be verified by comparing it with the result obtained by plotting the ofun w.r.t k and observing its minimum value, the value of k corresponding to the minimum value of ofun is found to be 1.236. This value is very close to that obtained by PSO, but since PSO gives a more precise value so it can be said that PSO is a much better technique to obtain the value of k. It can be observed that the value of the fitness function converges to its minimum value before the 50th iteration.
Footnotes
Acknowledgment
World Bank Technical Education Quality Improvement Program (TEQIP)-III Project of Rajkiya Engineering College, Ambedkar Nagar (U. P., India), financially supported this work.-224122.
