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

Abstract

To solve the problem of parameter selection during the design of magnetically 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.

Sheng Lu 1 and Chenyang Zuo 1 and Changhao Piao 1, 2

Academic Editor:Salvatore Alfonzetti

1, Institute of Pattern Recognition and Applications, Chong Qing University of Posts and Telecommunications, Chongqing 400065, China
2, Department of Mechanical Engineering, Inha University, Incheon 402-751, Republic of Korea

Received 13 March 2015; Revised 26 May 2015; Accepted 28 June 2015; 25 October 2015

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. 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 [figure omitted; refer to PDF] 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-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 ( [figure omitted; refer to PDF] ) and mutation probability ( [figure omitted; refer to PDF] ) according to the variance of population fitness makes [figure omitted; refer to PDF] and [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] and [figure omitted; refer to PDF] 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.

2. 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.

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

[figure omitted; refer to PDF]

The analysis object of this paper is the transmitter series-the receiver series topology. The equivalent circuit is shown in Figure 2. [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , respectively, represent transmitting coil and receiving coil's equivalent inductances. [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , respectively, represent the equivalent resistances of transmitting coil and receiving coil. [figure omitted; refer to PDF] is the output electromotive force of source of the inverter circuit and [figure omitted; refer to PDF] is the output resistance of the inverter circuit. [figure omitted; refer to PDF] and [figure omitted; refer to PDF] indicate the resonant capacitor of transmitter and the receiver. [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are the current of the lunching circuit and receiving circuit. [figure omitted; refer to PDF] is load resistance. [figure omitted; refer to PDF] is the mutual inductance between transmitting coil and receiving coil. [figure omitted; refer to PDF] expresses the distance of transmitting coil and receiving coil. When we add [figure omitted; refer to PDF] 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.

Figure 2: The equivalent circuit of SS structure.

[figure omitted; refer to PDF]

According to the model of mutual inductance equivalent circuit, we can establish the following equations: [figure omitted; refer to PDF]

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 [figure omitted; refer to PDF] and [figure omitted; refer to PDF] is ignored. The self-inductance of the transmitting coil and receiving coil is [figure omitted; refer to PDF] , and impedance of the transmitting loop and receiving loop is, respectively, represented by [figure omitted; refer to PDF] and [figure omitted; refer to PDF] . The impedance reflecting from the receiving end to the transmitting end is [figure omitted; refer to PDF] . And the mutual inductance of the transmitting coil and receiving coil is [figure omitted; refer to PDF] [5]. Thus, there exists the following relationship: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] , respectively, represent the loop number, coil wire diameter, coil radius, and the distance of coils. The current of the receiver circuit is [figure omitted; refer to PDF] The current of the transmitter circuit is [figure omitted; refer to PDF] The power of the load is [figure omitted; refer to PDF] The power of the main input is [figure omitted; refer to PDF] The efficiency of the transmission system is [figure omitted; refer to PDF]

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 [figure omitted; refer to PDF] . When the coupling coefficient [figure omitted; refer to PDF] , transmission power reaches to maximum at resonant frequency point [figure omitted; refer to PDF] . But with the decrease of [figure omitted; refer to PDF] , the maximum power will be with much attenuation. If the coupling coefficient [figure omitted; refer to PDF] , the maximum transmitted power will appear on both sides of [figure omitted; refer to PDF] . By this time, with the decrease of [figure omitted; refer to PDF] , the maximum transmission power point will be closer to [figure omitted; refer to PDF] . A little deviation from [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] of SS system needs to satisfy the conditions as follows: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the quality factor of transmitter circuit and [figure omitted; refer to PDF] is the quality factor of receiver circuit. So, as is shown in (8), if the excitation frequency is controlled as [figure omitted; refer to PDF] , high transmission power and distance will be obtained. As mentioned above, a little deviation from [figure omitted; refer to PDF] will not influence the transmission power apparently.

Let [figure omitted; refer to PDF] and [figure omitted; refer to PDF] represent the rated voltage of [figure omitted; refer to PDF] and [figure omitted; refer to PDF] . The rated value of [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , in other words the minimum rated current of [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , respectively, is represented by [figure omitted; refer to PDF] and [figure omitted; refer to PDF] . The parameters of the system need to satisfy the following constraints: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , respectively, represent the voltage of [figure omitted; refer to PDF] and [figure omitted; refer to PDF] . According to the actual working situation, [figure omitted; refer to PDF] and [figure omitted; refer to PDF] have the maximum values [figure omitted; refer to PDF] and [figure omitted; refer to PDF] and the minimum values [figure omitted; refer to PDF] and [figure omitted; refer to PDF] . There are size restrictions in coil; that is to say, the turns and radius of coil have the maximum values [figure omitted; refer to PDF] and [figure omitted; refer to PDF] and the minimum values [figure omitted; refer to PDF] and [figure omitted; refer to PDF] . The switch tube of high-frequency inverter circuit has the biggest switching frequency which is the maximum frequency [figure omitted; refer to PDF] of the AC power excitation. Meanwhile, the minimum working frequency [figure omitted; refer to PDF] is set. As maximum distance value [figure omitted; refer to PDF] and minimum distance value [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF]

The maximum, minimum, and rated values of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and the above parameters are known. In order to reach the highest transmission efficiency and transmission distance of the system, the parameters [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are optimized. The optimized objective function is as follow: [figure omitted; refer to PDF]

Inserting (2)-(4) into (8)-(10), we get the constriction of the system: [figure omitted; refer to PDF]

3. 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.

3.1. 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 [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are chosen randomly to form the population. It can be shown as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are, respectively, the [figure omitted; refer to PDF] component values of individuals [figure omitted; refer to PDF] and [figure omitted; refer to PDF] and [figure omitted; refer to PDF] and [figure omitted; refer to PDF] represent the maximum and minimum of the [figure omitted; refer to PDF] component, respectively. By judging whether the [figure omitted; refer to PDF] is smaller than [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , the larger the value [figure omitted; refer to PDF] is, the smaller the similarity of the individuals [figure omitted; refer to PDF] and [figure omitted; refer to PDF] is. And the smaller the value [figure omitted; refer to PDF] is, the larger the similarity of the individuals [figure omitted; refer to PDF] and [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the objective function of the worst individual among the current feasible solutions and [figure omitted; refer to PDF] 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.

3.2. Improved Crossover and Mutation 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 [figure omitted; refer to PDF] is relatively high and mutation probability [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] is relatively small and mutation probability [figure omitted; refer to PDF] is relatively high. It can avoid precocious phenomenon. Combining with logistic mathematical equation, the methods of calculating [figure omitted; refer to PDF] and [figure omitted; refer to PDF] can be obtained by [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the variance of population fitness and crossover probability [figure omitted; refer to PDF] , whose available interval is [figure omitted; refer to PDF] , decreases according to the decrease of [figure omitted; refer to PDF] . On the contrary, the mutation probability [figure omitted; refer to PDF] , whose available interval is [figure omitted; refer to PDF] , increases according to the decrease of [figure omitted; refer to PDF] . With the deep evolution, the individuals are more integrated, which shows that [figure omitted; refer to PDF] is smaller and smaller and the variation tendency of [figure omitted; refer to PDF] and [figure omitted; refer to PDF] conforms to the rules of genetic algorithm.

The crossover operation uses heuristic crossover and selects [figure omitted; refer to PDF] and [figure omitted; refer to PDF] as the individuals needing crossover operation according to crossover probability [figure omitted; refer to PDF] [18]. The new individuals after crossover operation are [figure omitted; refer to PDF] and [figure omitted; refer to PDF] ; that is, [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] is the current algebra, and [figure omitted; refer to PDF] 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 floating-point encoding, to do the operation of individuals mutating [19]. The new individuals after mutation operation are [figure omitted; refer to PDF] : [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] is the lower bound of [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] is the upper bound of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] and [figure omitted; refer to PDF] 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.

3.3. 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 [figure omitted; refer to PDF] is generated from the domain [figure omitted; refer to PDF] of individual [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] . is calculated. If [figure omitted; refer to PDF] , the new individual [figure omitted; refer to PDF] replaces the old individual [figure omitted; refer to PDF] ; namely, [figure omitted; refer to PDF] . If [figure omitted; refer to PDF] , [figure omitted; refer to PDF] will be calculated. If [figure omitted; refer to PDF] is true, the new individual will be accepted.

The cooling operation in the simulated annealing algorithm is as follows: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is positive constant, slightly less than 1, and [figure omitted; refer to PDF] is cooling times. [figure omitted; refer to PDF] 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.

4. 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.

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

Parameters	Nominal value
[figure omitted; refer to PDF] /Ω	50
[figure omitted; refer to PDF] /Ω	5
[figure omitted; refer to PDF] /V	40
[figure omitted; refer to PDF] /mm	1.6
[figure omitted; refer to PDF] and [figure omitted; refer to PDF] /V	500
[figure omitted; refer to PDF] /A	4
[figure omitted; refer to PDF] /A	3
[figure omitted; refer to PDF] and [figure omitted; refer to PDF] /nF	[figure omitted; refer to PDF] 0.09, 3 [figure omitted; refer to PDF]
[figure omitted; refer to PDF]	[figure omitted; refer to PDF] 2, 10 [figure omitted; refer to PDF]
[figure omitted; refer to PDF] /mm	[figure omitted; refer to PDF] 50, 80 [figure omitted; refer to PDF]
[figure omitted; refer to PDF] /mm	[figure omitted; refer to PDF] 50, 100 [figure omitted; refer to PDF]
[figure omitted; refer to PDF] /KHz	[figure omitted; refer to PDF] 2000, 5000 [figure omitted; refer to PDF]

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

(a) The initial population of 100

[figure omitted; refer to PDF]

(b) The initial population of 150

[figure omitted; refer to PDF]

(c) The initial population of 200

[figure omitted; refer to PDF]

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 [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , 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 [figure omitted; refer to PDF] and the transmission efficiency is [figure omitted; refer to PDF] . And the corresponding optimal parameters are [figure omitted; refer to PDF] mm, [figure omitted; refer to PDF] rad/s ( [figure omitted; refer to PDF] kHz), [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] 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, [figure omitted; refer to PDF] μ H, [figure omitted; refer to PDF] pF, and [figure omitted; refer to PDF] μ 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 [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are 1.24 A, 0.92 A, 497.4 V, and 366.5 V, respectively. Obviously, the capacitor voltage and current satisfy the requirements.

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

[figure omitted; refer to PDF]

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 [figure omitted; refer to PDF] , [figure omitted; refer to PDF] mm, and [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] μ H and [figure omitted; refer to PDF] μ 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 5: The coil model with finite-element mesh generation.

[figure omitted; refer to PDF]

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.

Figure 6: The finite-element model of power transfer system.

[figure omitted; refer to PDF]

Figure 7: The magnetic field distributions of resonant point (a) and unresonant point (b), respectively.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

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.

Figure 8: The relationships of [figure omitted; refer to PDF] with [figure omitted; refer to PDF] under the different parameters.

[figure omitted; refer to PDF]

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.

5. 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.

Acknowledgments

This work is supported by the Foundation of Chongqing Science and Technology Commission (CSTC 2013jcyjjq60002, CSTC 2013yykfC60005, and CSTC 2013 jcsf-jcssX0022) and the National Natural Science Foundation of China (NSFC 11247325).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] Z. M. Zhao, Y. M. Zhang, K. N. Chen, "New progress of magnetically-coupled resonant wireless power transfer technology," Proceedings of the Chinese Society of Electrical Engineering , vol. 33, no. 3, pp. 1-13, 2013.

[2] Y. Xie, Z. G. Wang, Y. S. Zhang, "Development of battery system testing machine," International Journal of Future Engineering , vol. 2015, 2015.

[3] Y.-S. Seo, M. Q. Nguyen, Z. Hughes, S. Rao, J.-C. Chiao, "Wireless power transfer by inductive coupling for implantable batteryless stimulators," in Proceedings of the IEEE MTT-S International Microwave Symposium Digest, pp. 1-3, Montreal, Canada, June 2012.

[4] M. Ali, G. Yang, R. Dougal, "A miniature packaged rectenna for wireless power transmission and data telemetry," in Proceedings of the IEEE International Workshop on Antenna Technology Small Antennas and Novel Metamaterials (IWAT '06), pp. 225-228, IEEE, March 2006.

[5] A. Kurs, A. Karalis, R. Moffatt, J. D. Joannopoulos, P. Fisher, M. Soljacic, "Wireless power transfer via strongly coupled magnetic resonances," Science , vol. 317, no. 5834, pp. 83-86, 2007.

[6] S.-H. Lee, R. D. Lorenz, "Development and validation of model for 95%-efficiency 220-W wireless power transfer over a 30-cm air gap," IEEE Transactions on Industry Applications , vol. 47, no. 6, pp. 2495-2504, 2011.

[7] O. Jonah, S. V. Georgakopoulos, "Wireless power transfer in concrete via strongly coupled magnetic resonance," IEEE Transactions on Antennas and Propagation , vol. 61, no. 3, pp. 1378-1384, 2013.

[8] X. Zhang, Q. Yang, H. Chen, Y. Li, X. Zhang, L. Jin, "Research on characteristics of frequency splitting in electromagnetic coupling resonant power transmission systems," Proceedings of the Chinese Society of Electrical Engineering , vol. 32, no. 9, pp. 167-172, 2012.

[9] Y. Li, Q. X. Yang, Z. Yan, C. Zhang, H. Y. Chen, "Analysis on effective range of wireless power transfer and its impact factors," Transactions of China Electrotechnical Society , vol. 28, no. 1, pp. 106-112, 2013.

[10] C. Y. Xia, Y. Sun, N. Jia, Y. Zhuang, Y. Zhai, "Magnetic circuit parameter optimization for coupled magnetic resonances power transfer system," Transactions of China Electrotechnical Society , vol. 27, no. 11, pp. 139-145, 2012.

[11] L. H. Chen, S. Liu, Y. C. Zhou, T. J. Cui, "An optimizable circuit structure for high-efficiency wireless power transfer," IEEE Transactions on Industrial Electronics , vol. 60, no. 1, pp. 339-349, 2013.

[12] D. Ahn, S. Hong, "A transmitter or a receiver consisting of two strongly coupled resonators for enhanced resonant coupling in wireless power transfer," IEEE Transactions on Industrial Electronics , vol. 61, no. 3, pp. 1193-1203, 2014.

[13] J. B. Hao, Z. B. Li, L. P. Zhao, "Genetic algorithm with constraint model for sequencing workingsteps on machining center," Journal of Xi'an Jiaotong University , vol. 42, no. 7, pp. 860-864, 2008.

[14] C. H. Piao, C. Jiang, H. T. Qiao, C. Cho, S. Lu, "Modeling and implementation of fixed switching frequency sliding mode control to two-stage DC-DC converter," International Journal of Control, Automation and Systems , vol. 12, no. 6, pp. 1225-1233, 2014.

[15] C. H. Piao, C. Jiang, H. T. Qiao, C. D. Cho, S. Lu, "Study on the control algorithm of two-stage DC-DC converter for electric vehicles," Advances in Mechanical Engineering , vol. 6, 2014.

[16] X. K. Zhu Study on the equivalent circuit of wireless power transfer system based on Witricity technology [M.S. thesis] , Nanchang University, 2012.

[17] K. Deb, "An efficient constraint handling method for genetic algorithms," Computer Methods in Applied Mechanics and Engineering , vol. 186, no. 2-4, pp. 311-338, 2000.

[18] S. H. Tang, Q. J. Zhu, "Effects of the initial population, crossover and mutation rate to the results of genetic algorithms and a possible solution scheme," Bulletin of Science and Technology , vol. 17, no. 3, 2001.

[19] Z. Michalewicz, C. Z. Janikow, J. B. Krawczyk, "A modified genetic algorithm for optimal control problems," Computers & Mathematics with Applications , vol. 23, no. 12, pp. 83-94, 1992.

[20] K. Jung, Y. Kim, S. Lim, B. Ahn, G. Lee, "Effect of friction on critical damage value in compressive deformation," International Journal of Future Engineering , vol. 2015, 2015.