The Parameters Optimization of MCR-WPT System Based on the Improved Genetic Simulated Annealing Algorithm

To solve the problemof parameter selection during the design ofmagnetically coupled resonant wireless power transmission system (MCR-WPT), this paper proposed an improved genetic simulated annealing algorithm. Firstly, the equivalent circuit of the system is analysis in this study and a nonlinear programming mathematical model is built. Secondly, in place of the penalty function method in the genetic algorithm, the selection strategy based on the distance between individuals is adopted to select individual. In this way, it reduces the excess empirical parameters. Meanwhile, it can improve the convergence rate and the searching ability by calculating crossover probability and mutation probability according to the variance of population’s fitness. At last, the simulated annealing operator is added to increase local search ability of the method. The simulation shows that the improved method can break the limit of the local optimum solution and get the global optimum solution faster. The optimized system can achieve the practical requirements.


Introduction
Traditional electric power transmission device mainly transmits directly connected via the metal wire, and this kind of connection will cause many problems [1].With the electric vehicle booming becoming an important green industry, it has received more and more focus [2].The wireless charging technology could promote the popularity of electric vehicle and solve the problems met in wire charging process.So far, wireless power transmission is divided into electromagnetic induction, magnetic coupling resonant, and radio waves according to the power transmission principle [1].Electromagnetic induction wireless power through the loosely coupled transformer or discrete transformer is mainly based on the principle of electromagnetic induction [3].Radio waves wireless power transmission technology is mainly using the principle that the electromagnetic wave power can be transmitted and received through antenna [4].Magnetic coupling resonant is a new kind of wireless power transmission technology that is proposed by MIT's Professor Marin Soljacic in 2007 [5].To realize the high effective power transmission, Marin used a couple of transmitting and receiving coils, which have the same resonant frequency, and put the receiving coils in the near field of the transmitting coils.However, radio waves and electromagnetic induction wireless power transmission have the disadvantages of that the transmission power, transmission distance and transmission efficiency cannot reaches to maximum at the same time.By comparison, magnetic coupling resonant wireless power transmission is concerned with the far transmission distance, safety, reliability, high efficiency, and many other advantages.It is considered to be one of the most potential embodiments of wireless power transmission.Nowadays, more and more researchers contribute to this technology [6,7].
The mutual inductance of magnetic coupling resonant wireless power transmission system is very small.In order to improve transmission performance, it requires a high quality factor  and MHz level operating frequency.In addition, it needs to select the most appropriate circuit parameters, because unreasonable parameters selection will damage the device or cause frequency bifurcation [8].To make the system work stably, reliably, and efficiently, the parameters must satisfy some constraints.Papers [9][10][11] have studied the MCR-WPT parameters' design.These studies have not taken the constraints of the parameters into account comprehensively.Usually, low efficiency is generated by using stepwise design method.Magnetic coupling resonant wireless power transmission system is a kind of multiparameter and multiconstrained nonlinear system [12].The complexity of the system equation will increase with the number of the coils.When the traditional optimization methods are used to solve the system, the system model needs to be simplified.The traditional methods are based on the system function model derivative, which does not fully consider the constraints of the components' working conditions for the optimized parameters, so the parameters cannot be directly used in practice and need to be further amended.Nevertheless, only one parameter is optimized in every amendment, not all of them.As a result, it cannot guarantee getting global optimal solution when making a few parameters together in this way.
In order to achieve the maximums of transmission efficiency and distance, respectively, this paper fully considers the constraints of the system to establish the mathematical model and proposes an improved genetic simulated annealing algorithm to optimize parameters simultaneously.Genetic algorithm (GA), which does not require that the objective function is derivative, can be optimized in multiparameter and multiconstraint conditions.Therefore, it shows a good performance on robustness [13].GA performs better in global search than in local search.Although GA is able to search the globally optimal solution in a random way based on probability, in reality, there exist problems such as premature phenomenon and low ability in local optimization.To improve the searching abilities of GA, firstly, the selection strategy based on the distance between individuals was used.Thus, it can avoid the loss of population diversity and the fitness value can be calculated according to feasibility-based rules.Also, it simplifies the setting of too many parameters.Then, the method of calculating crossover probability (  ) and mutation probability (  ) according to the variance of population fitness makes   and   change with the increasing of iteration.It really makes sense of raising the rate of convergence.Simulated annealing algorithm (SA) simulates the process of metal annealing, which can reach a commendable result by setting a high initial temperature and a slow process of annealing.SA has a good ability of local search as well as convergence of globally optimal solution.Meanwhile, SA is able to break away from the local optimal solution and restrain the premature phenomenon.Improved genetic simulated annealing algorithm (IGSA) combines the advantages of GA and SA.At the same time, it adds an improved selection strategy and   and   strategy which can increase the ability of searching and the rate of convergence.IGSA is used for parameters optimization design for the series-series system structure.The comparison shows that the optimization ability of IGSA is better than the general algorithm and ANSYS finite-element software is used for simulation verification.

The Mathematical Model of the MCR-WPT System
As shown in Figure 1, the transmitter of magnetically coupled resonant wireless power transfer system is composed of AC power, rectifier filter circuit, high-frequency inverter circuit, transmitting coil, and resonant capacitor.Rectifier filter circuit converts alternating current (AC) to direct current (DC) for high-frequency inverter circuit, and then it will output high-frequency current to transmitting coil.The receiver consists of receiving coil, resonant capacitor, rectifier circuit, and load.Rectifier filter circuit converts alternating current (AC) occurring in receiving coil to direct current (DC) for load.On special occasion, DC-DC converter may be introduced into the rectifier filter circuit to convert the high voltage and small current input into large current and low voltage output [14,15].The power between transmitting terminal and receiving end transfers through the magnetic coupling resonance.The analysis object of this paper is the transmitter seriesthe receiver series topology.The equivalent circuit is shown in Figure 2.   and   , respectively, represent transmitting coil and receiving coil's equivalent inductances. 1 and  2 , respectively, represent the equivalent resistances of transmitting coil and receiving coil.  is the output electromotive force of source of the inverter circuit and   is the output resistance of the inverter circuit.  and   indicate the resonant capacitor of transmitter and the receiver.  and   are the current of the lunching circuit and receiving circuit.  is load resistance. is the mutual inductance between transmitting coil and receiving coil. expresses the distance of transmitting coil and receiving coil.When we add   to the ends of the transmitting coil, it can produce the same frequency in the coil magnetic field.Receiving coil is in the near field region of transmitting coil and these two coils have the same frequency.As a result, it can create magnetic coupling resonance.According to the model of mutual inductance equivalent circuit, we can establish the following equations: In the actual system, for ease of design and analysis, the transmitting coil and receiving coil are all designed as the parallel and coaxial structure.Meanwhile, the influence of coil resistances  1 and  2 is ignored.The self-inductance of the transmitting coil and receiving coil is , and impedance of the transmitting loop and receiving loop is, respectively, represented by   and   .The impedance reflecting from the receiving end to the transmitting end is   .And the mutual inductance of the transmitting coil and receiving coil is  [5].Thus, there exists the following relationship: where , , , and , respectively, represent the loop number, coil wire diameter, coil radius, and the distance of coils.The current of the receiver circuit is The current of the transmitter circuit is The power of the load is The power of the main input is The efficiency of the transmission system is According to the frequency and quality factor of system, paper [16] has studied the influence of transmission power.The critical coupling coefficient is proposed as When the coupling coefficient  <   , transmission power reaches to maximum at resonant frequency point  0 .But with the decrease of , the maximum power will be with much attenuation.If the coupling coefficient  <   , the maximum transmitted power will appear on both sides of  0 .By this time, with the decrease of , the maximum transmission power point will be closer to  0 .A little deviation from  0 will not make transmission power decline.In order to control the working frequency of the system easily at a resonance point and ensure the transmission distance and high transmission efficiency, the coupling coefficient  of SS system needs to satisfy the conditions as follows: where   =  0   /  is the quality factor of transmitter circuit and   =  0   /  is the quality factor of receiver circuit.So, as is shown in (8), if the excitation frequency is controlled as  0 , high transmission power and distance will be obtained.As mentioned above, a little deviation from  0 will not influence the transmission power apparently.
Let    -rate and    -rate represent the rated voltage of   and   .The rated value of   and   , in other words the minimum rated current of   and   , respectively, is represented by  -rate and  -rate .The parameters of the system need to satisfy the following constraints: where The maximum, minimum, and rated values of   ,   ,   , and the above parameters are known.In order to reach the highest transmission efficiency and transmission distance of the system, the parameters , , , and  0 = 2 0 are optimized.The optimized objective function is as follow: Inserting ( 2)-( 4) into ( 8)- (10), we get the constriction of the system:

MCR-WPT System Parameters Optimization Based on Improved Genetic Simulated Algorithm
Genetic algorithm achieves a good performance in global search, but it is easy to trap into local optimal solution.Simulated annealing algorithm performs well in local search, but its convergence speed is slow.To solve the parameters optimization problems of systems with multiple constraints and parameters, the advantages of the two algorithms are combined.The idea of using genetic simulated algorithm to optimize parameters of MCR-WPT system is described as follows.Firstly, the feasibility-based rules are applied to selecting operation to prevent the loss of population diversity.Secondly, to improve the search capability and convergence speed of the algorithm, population fitness variance is used to calculate crossover probability and mutation probability.
The heuristic cross and nonuniform mutations are used for individual crossover and mutation operation.Finally, the simulated annealing operation is carried out to prevent the algorithm from early convergence.

Improved Selection Strategy.
Deb proposed a method to handle GA feasibility selection operation [17].The key point is that the method can solve feasible solutions.The main requirements of the method are as follows: (1) the feasible solution is superior to the nonfeasible solution; (2) comparing the two feasible solutions, the one with optimal objective value is desired; (3) comparing the two nonfeasible solutions, the one with low degree of default function is desired.Compared with penalty function method, there are no extra empirical parameters used in the method, and a championship method is used to select the superior individuals.But the method falls into local solution easily.Because of the neglect of nonfeasible solution, the evolutional population is mainly composed of feasible solutions.Thus, the similarity of population individuals is quite high and the individuals are relatively single.Through every selecting operation, the individuals of next generation are mainly feasible ones, which reduce the population diversity.To keep population diversity, the selected individuals should avoid being all feasible ones or feasible and nonfeasible ones.In this paper, the selection strategy is improved.The selecting operation of feasibility for individuals is decided by the distance-based method.Firstly, to calculate Euclidean distance and complete normalization processing, two individuals  and  are chosen randomly to form the population.It can be shown as where    and    are, respectively, the  th component values of individuals  and  and  max  and  min  represent the maximum and minimum of the  th component, respectively.By judging whether the  , is smaller than ℎ, ℎ ∈ (0, 1), the larger the value ℎ is, the smaller the similarity of the individuals  and  is.And the smaller the value ℎ is, the larger the similarity of the individuals  and  is.If the condition does not satisfy the judgment, the individual is to be selected again.If it does, the individual fitness will be calculated according to where  max is the objective function of the worst individual among the current feasible solutions and  is the number of constraint conditions.Lastly, according to Deb's feasible rule, these two individuals are compared and we can determine which individual will enter the crossover and mutation operations of next generation.Because similar individuals are compared every time, population diversity is kept.Thus, it will be conducive to individuals distributing in a wide range space and contribute to global search.

Improved Crossover and Mutation
where  is the variance of population fitness and crossover probability   , whose available interval is [0.4,0.88], decreases according to the decrease of .On the contrary, the mutation probability   , whose available interval is [0.16, 0.3], increases according to the decrease of .With the deep evolution, the individuals are more integrated, which shows that  is smaller and smaller and the variation tendency of   and   conforms to the rules of genetic algorithm.
The crossover operation uses heuristic crossover and selects   and   as the individuals needing crossover operation according to crossover probability   [18].The new individuals after crossover operation are    and    ; that is, where  =   / max ,   is the current algebra, and  max is the maximum algebra.
The crossover operation changes with the deep evolution, and the crossover effect gradually decreases.This crossover operation is conducive to global search in the early evolution and local optimization solution search in the late mutation operation.Meanwhile, it can prevent the excellent individuals from being damaged.
Mutation operation determines the local search ability of algorithm.In this paper, we use nonuniformity mutation operator, which is proposed by Michalewicz for floatingpoint encoding, to do the operation of individuals mutating [19].The new individuals after mutation operation are    : where   ∈ [  ,   ],   is the lower bound of , and   is the upper bound of ,  1 and  2 is random number between 0 and 1.The mutation step in the operator decreases adaptively with the deep evolution.In the early stage of iteration, the big mutation step improves search efficiency.In the late stage, the small mutation step search keeps the optimal solution.

Simulated Annealing Operation.
The individuals after crossover and mutation operation use simulated annealing operation.The Metropolis criterion is used to accept those new individuals.The specific operations are as follows.The new individual   is generated from the domain (  ) of individual   , and Δ = (  )−(  ). is calculated.If Δ ≤ 0, the new individual   replaces the old individual   ; namely,   :=   .If Δ > 0,  −Δ/  > random(0, 1) will be calculated.If  −Δ/  > random(0, 1) is true, the new individual will be accepted.
The cooling operation in the simulated annealing algorithm is as follows: where  is positive constant, slightly less than 1, and  is cooling times.(0) indicates the initial temperature and usually is a very large value.The simulated annealing operation does not always negate nonfeasible solution but would accept it in a certain probability.This is conducive to solve the problem of search stagnation in the late stage of the algorithm and prevent the optimization results from locally optimal solution.

MCR-WPT System of Simulations and Finite-Element Verification
In this section, three methods including simple genetic algorithm (SGA), simple genetic simulated annealing algorithm (SGSA), and IGSA have been applied to optimize system's parameters.With the initial population of 100, 150, and 200, respectively, the optimal optimized results are shown in Figure 3.The specifications of the magnetic coupling resonant wireless power transmission system based on SS topology are provided in Table 1.

Effective value
Effective value  From Figure 3, it is known that the rate of convergence of the system with SGA is small.It plunges into local optimum easily and the optimization is useless.Also, we can see that the local search ability of the system is enhanced by SGSA.However, the rate of convergence is small and the optimization effect is worse than using IGSA.Thanks to improved selection strategy of individuals and calculation method of   and   , not only the local search ability of the system improves observably, but also the rate of convergence changes faster and the optimization is better than using the other two methods.
With the initial population of 200, the optimized maximum is  obj = 1.2744 and the transmission efficiency is  = 84.4%.And the corresponding optimal parameters are  = 80 mm,  0 = 27.73 × 10 6 rad/s ( 0 = 4414.2kHz),  = 6, and  = 71.5 mm.By substituting these parameters into (10), it is known that they meet the constraint conditions.
According to optimal parameters, the other parameters are obtained; that is,   =   =  = 14.5 H,   =   =  = 89.7 pF, and  = 1.33 H.Further, based on these values, the model and simulations are achieved by Simulink.Figure 4 shows the voltage and current waveforms of the transmitting and receiving coils.It is clear that   ,   ,    , and    are 1.24 A, 0.92 A, 497.4 V, and 366.5 V, respectively.Obviously, the capacitor voltage and current satisfy the requirements.
The finite-element method is usually applied in force and electromagnetism analysis [20].In order to verify the consistency between coil inductance and mutual inductance when the coil has a feature with  = 6,  = 80 mm, and  = 1.6 mm, we use ANSYS to simulate the system.In this paper, PLANE53 element is selected as the element attribute of coil and air, and an axisymmetric model was employed for system modeling.The finite-element mesh generation for the established model is shown in Figure 5.After exerting boundary conditions and excitation, we get  = 12.9 H and  = 1.3 H by LMATIRX.Obviously, the differences between the simulated values and the theoretical values are very small.The reason for small differences is to simplify the coil model.The cross section of the coil model is rectangle, which is larger than the actual coil.
Figure 6 shows the system's finite-element model with PLANE53 and CIRCUIT124 elements.Based on the optimal parameters and adopting ANSYS electromagnetic analysis, the simulation result is shown in Figure 7.As the figure shows, the coupling degree of the resonance point is weaker than that of the unresonant points.But the overall intensity of the magnetic field in the resonance point is more than 2 times larger than that in the unresonant point.
Through ANSYS postprocessing function, the values of input power and load power in each frequency point are extracted.So the curve of the relationship between the transmission efficiency and frequency can be obtained.As  is shown in Figure 8, the transmission efficiency reached the maximum at the resonant frequency.Due to the small simulated inductance value, the resonant frequency would be larger than the theoretical value.
In order to validate that the transmission efficiency can reach the maximum based on the optimal parameters, without considering the constraint conditions, one of the parameters is changed while the others are kept constant.In this way, three transmission efficiency curves can be obtained as shown in Figure 8. Obviously, system transmission efficiency reaches the maximum by using the optimal parameters.In the meantime, the efficiency curve crest is relatively flat.So the attenuation changes inconspicuously with the deviation of the resonance point at ±100 Khz.But, for the other three parameters, the transmissions efficiencies never reached the maximum and the efficiency curve crests are all very sharp.So the attenuation is influenced apparently if a little deviation of resonant frequency point appears.According to the analysis above, the system can achieve a good transmission performance by using the optimal parameters.Meanwhile, the system can control the frequency easily.It fully meets the system design requirements.

Conclusion
In this paper, at first, by analyzing the SS topology of the MCR-WPT system, the parameters expressions are obtained and the corresponding boundary conditions and objective functions are established.Then, IGSA is used to optimize the system's parameters.The proposed algorithm can decide which individuals would be chosen for the next genetic manipulation by calculating the distances among individuals.The method increases the population individual multiplicity and prevents local convergence.Besides, it can adjust the crossover and mutation probabilities based on the fitness variances of the individuals in every generation.Thus, the convergence speed increases.Finally, simulated annealing operation for progeny is implemented to improve the search ability of the algorithm.The ANSYS finite-element model is designed to verify the optimized results.It verifies that the system with optimal parameters can achieve optimal performances and it can control the frequency easily.

Figure 1 :
Figure 1: The schematic diagram of MCR-WPT system.

Figure 2 :
Figure 2: The equivalent circuit of SS structure.

Figure 3 :
Figure 3: The solved objective function curves with different initial population.

Figure 4 :
Figure 4: The voltage and current waveforms of the transmitting and receiving coils with optimal parameters.

Figure 5 :
Figure 5: The coil model with finite-element mesh generation.
and    , respectively, represent the voltage of   and   .According to the actual working situation,   and   have the maximum values  -max and  -max and the minimum values  -min and  -min .There are size restrictions in coil; that is to say, the turns and radius of coil have the maximum values  max and  max and the minimum values  min and  min .The switch tube of high-frequency inverter circuit has the biggest switching frequency which is the maximum frequency  0-max of the AC power excitation.Meanwhile, the minimum working frequency  0-min is set.As maximum distance value  max and minimum distance value  min exist between transmission coil and receiving coil according to the working condition of lord and transmitting terminal, the constricts of each parameter can be expressed in  -min <   <  -max ,  -min <   <  -max ,  min <  <  max , Operation.At the beginning, the algorithm needs strong global search ability and weak local search ability, which are helpful for searching in global range.Thus, at the early stage of the algorithm, the crossover probability   is relatively high and mutation probability   is relatively small.At the late stage of the algorithm, in order to converge to optimization solution, the global search ability needs to reduce and local search to be strengthened.So, at the late stage of the algorithm, the crossover probability   is relatively small and mutation probability   is relatively high.It can avoid precocious phenomenon.Combining with logistic mathematical equation, the methods of calculating   and   can be obtained by

Table 1 :
The parameter values of WCR-WPT system.