Macroscopic/Mesoscopic Computational Materials Science Modeling and Engineering
View this Special IssueResearch Article  Open Access
Multiobjective Design of Turbo Injection Mode for Axial Flux Motor in Plastic Injection Molding Machine by Particle Swarm Optimization
Abstract
This paper proposes a turbo injection mode (TIM) for an axial flux motor to apply onto injection molding machine. Since the injection molding machine requires different speed and force parameters setting when finishing a complete injection process. The interleaved winding structure in the motor provides two different injection levels to provide enough injection forces. Two wyewye windings are designed to switch two control modes conveniently. Wyewye configuration is used to switch two force levels for the motor. When only one set of wyewinding is energized, field weakening function is achieved. Both of the torque and speed increase under field weakening operation. To achieve two control objectives for torque and speed of the motor, fuzzy based multiple performance characteristics index (MPCI) with particle swarm optimization (PSO) is used to find out the multiobjective optimal design solution. Both of the torque and speed are expected to be maximal at the same time. Three control factors are selected as studied factors: winding diameter, winding type, and airgap. Experimental results show that both of the torque and speed increase under the optimal condition. This will provide enough large torque and speed to perform the turbo injection mode in injection process for the injection molding machine.
1. Introduction
Axialflux motor is widely applied in the technical field of directdrive mechanism. Different from the radialflux motor, the axialflux type motor has larger torque, shorter axial length, and more compact structure. Therefore, the axialflux type motor is applied in many electric vehicle and industrial automation applications [1].
The fractional concentrated winding and fractional distributed winding have advantages of shorter wiring length, lower copper loss and lower m.m.f. harmonics [2]. Therefore, the axialflux motor has lower inductance values [3, 4]. If the motor requires operating under wider speed range, large demagnetizing current is required to generate the weak field effect [4].
The permanent magnet is embedded into rotor structure to provide more compact force [5]. Fractional slot design is selected to provide appropriate design for the stator structure [6]. Consequently, axialflux motor is usually used in high torque application [7].
This paper proposes a compact design for the injection molding process. Axialflux motor is used to provide the required injection force. Twowyewye winding structure is proposed to provide two different levels of speed range.
That means wider speed range is designed. By setting two energizing current commands and switch two modes for two wyewye windings, wider speed and higher torque are achieved. Turbo injection mode is designed by energizing winding current in one specific wye winding.
This will be very helpful to provide larger torque for the injection molding machine. The conventional injection molding machine requires a huge hydraulic pressure structure to provide enough larger force. However, the hydraulic pressure structure occupies very large space. The overall machine cannot be designed as a compact one. By using directdrive type axialflux motor, the injection molding machine can be designed more compactly.
2. Design of AxialFlux Motor in Plastic Molding Machine
The axialflux motor in the plastic injection molding machine is shown in Figures 1 and 2. The magnetic flux direction is in the axial direction. Two flat plates are placed facetoface. There are threephase windings on the stator plate. There are permanent magnets on the rotor plate.
24 slots and 16 poles are designed to provide the required magnetic propulsion force. The 3D model can be built up as shown from Figures 1 to 3. After the 3D model is built up, the geometric model is imported into 3D software to analyze the magnetic field and assess the associated magnetic torque.
(a) Stator
(b) Rotor
3. Motor Operation and Turbo Injection Mode
There are two techniques to design the turbo injection mode (TIM) characteristics of the motor: one is to regulate the field weakening and the other one is to regulate the torque constant [8].
Figure 4 shows the wyedelta winding configuration and control scheme for six power switches [9, 10]. The wye configuration can provide the required initial speedup and required climbing torque demand. When the motor goes to the steady state, the motor changes its connection state from the wye configuration into delta configuration [11, 12].
The electric impedance of wye configuration is 3 times of the one of delta configuration. The back E.M.F. of wye configuration is 1.732 times of the one of delta configuration. Therefore, the torque of wye configuration is 1.732 times of the one of delta configuration. The speed of delta configuration is 1.732 times of the one of wye configuration.
In order not to generate too higher instant acceleration, dualduty control method is used. It is easy for the wyewye configuration to prevent from generating too large transient current when switching from the wyewye configuration to only one wye winding configuration. The current control loop is also used to reduce the transient current for safety and protection [13, 14]. Figure 5 shows the block diagram of the overall control system [15].
4. Response Surface Method with Fuzzy MPCI Method
From practical applications, there are many control factors to increase the speed of the axialflux motor as shown in Table 1. Three factors are considered and selected in this paper. The first one is wiring dimension. The second one is winding type. The third one is airgap. The three factors are used to define the three control factors for the response surface method. The optimal solution is obtained by using response surface method. Figure 6 shows the membership function definition for the described design problem.



(a) Membership functions for speed
(b) Membership functions for torque
(c) Membership functions for MPCI
(d) Flow chart for MPCI method
The multiobjective design problem for the axialflux motor is to maximize the torque and to maximize speed at the same time. As shown in Table 1, there are three control factors selected for the associated motor design: the diameter of the wiring, the type of the wiring, and the width of airgap. For the response surface method, two noise factors are required for the associated ANOVA. They are insulating coefficient and the stator temperature. The required experimental data are provided in Tables 1~5.
5. Analysis of Response Surface Method
The average of the four centralpoint experimental values in Table 5 can be calculated as
The average of the response data for the sixteen experimental runs shown 6 in Table 4 can be calculated as


Therefore, the curvature sum of squares can be further calculated as follows:
The error sum of squares is computed as
statistical value is calculated as
The ANOVA analysis for the firstorder response surface model is evaluated by softwarebased regressive analysis. Due to limited pages, the ANOVA tables are neglected.
From the above equation, the curvature is very small. Therefore, the firstorder model can be used to evaluate the optimal parameters. The firstorder regressive model can be represented as follows. Assume that noise factors are zero, , and . The mean response surface can be obtained into
The mean response surface can be plotted as 3D surface plot in Figure 7(a). Also, the corresponding contour plot can be plotted in Figure 7(b).
(a) 3D surface plot
(b) Contour plot
The variance response surface can be formulated as follows:
Assume that and , the variance response surface can be rewritten into
The 3D surface plot of variance response surface can be plotted as shown in Figure 8(a). Also, the corresponding contour plot can be shown in Figure 8(b). The constraint of the variance response surface can be described intowhere independent variable ranges are ; ; .
(a) 3D surface plot
(b) Contour plot
After the optimization problem is formulated, the optimization process by particle swarm method is used to search for the optimal solution.
6. Particle Swarm Optimization Method
With the above system model derived by response surface method, the optimization process is performed to derive the optimal solution set for the described problem. However, the optimal solution may be located anywhere in the range between −1.0 and +1.0. The local optimal solution may not be the optimal in the entire global region. Therefore, the PSO process is used to derive the optimal solution in an efficient way. The local searching and global searching go simultaneously. Therefore, global optimal solution can be found.
By adding the random seeds into formulation, the PSO can jump out of the local optimal solution if the global solution is more optimal than the local solution. In the following, the OPSO (orthogonal particle swarm optimization) formulation is performed. Since in the response surface method, the nonlinear problem is formulated as firstorder model problem. However, the curvature for this model affects the precision when finding the optimal solution. That means the nonlinearity property exists in this problem. This influences the searching process when finding the optimal solution.
The particle swarm optimization originates from the emulation of the group dynamic behavior of animal. For each particle in a group, it is not only affected by its respective particle, but also affected by the overall group. There are position and velocity vectors for each particle. The searching method combines the experience of the individual particle with the experience of the group. For a particle or a point in a searching space with dimension [16, 17], the th particle associated with the problem can be defined aswhere and , PS is the population size. The respective particle objective function and group objective function associated with each are defined as [18, 19]
The refreshing speed vector can be defined as
The refreshing position and velocity vectors can be expressed as [20, 21]where .
When the searching begins, the initial solution is set. In the iteration process, the particle is updated by the value coming from the group particles and individual particle. The convergence condition is dependent on the minimum of the average square error of the particle. Both the experience of the individual particle and the experience of the group are mixed into the searching process [22, 23].
In the optimization problem, there might be a local minimum problem. The optimal solution might jump into a local trap and cannot jump out of the trap [24, 25]. Actually, a local minimum point does not represent a global solution in a wide range. In the group experience, random function is used to jump out of the local interval. An inertia weighting factor is considered in this algorithm to increase the convergence rate. An inertia weighting factor is added in the following expression [26]. The modified formula can be expressed aswhere the and are both constants. is the initial weighting value. is the final weighting value. gen is the number of current generation. is the number of final generation. However, the above mentioned is actually a kind of linear modification. To make the algorithm suitable for nonlinear searching problem, there is many nonlinear modification methods proposed to refresh the velocity vector. The modified term is defined as the key factor. By setting the learning factors and which are larger than 4.0, the modification for the speed vector is expressed into
A modified PSO method called orthogonal PSO (OPSO) is proposed to solve the update problem effectively. A simple orthogonal array in Taguchi method is used in this algorithm to help with the update [26].
7. Orthogonal Array Algorithm in OPSO Method
To run the Taguchi method, two functions are defined first. The particle swarms are composed of individual particle swarm and group particle swarm
These two functions are specified as two control factors in Taguchi method. Two levels are defined for the control factors. Therefore, the orthogonal array has two factors and two levels. The objective function calculating from the problem is used as the measured value in orthogonal array. By comparing the objective function value in the orthogonal array, the optimal solution can be obtained. Assume that the optimal solution is expressed as . The is adopted to refresh the particle position and velocity vectors as shown in the following expression. The particle refreshing process in OPSO optimization is illustrated in Figure 9
(a) Refreshing process from first generation to third generation
(b) Error percentage from first generation to fifth generation
(c) Illustrated results for optimal solution
The optimization problem is formulated and the optimal solution can be found. The optimal speed and optimal torque can be found at the same time by using fuzzy MPCI method. The derived optimal solution can provide optimal condition for the associated motor design.
By using the response surface method combined with PSO method, the mathematical model for this problem is provided and verified. This will be very helpful to assess the associated motor design application. The optimal solution can be found and located at the ends of the entire searching range. Global solution is found instead of local solution. Results show that the proposed mathematical method has the capability of providing appropriate searching process.
The local optimal solution is avoided and the proposed PSO algorithm can find the optimal solution at the endpoints globally. In PSO, global and local optimal ranges are searched at the same time. The convergence rate is faster than conventional searching method. The related confirmation experiments show that the proposed methodology can provide good prediction with the practical experimental case.
It is convinced that the proposed optimal parameter solution solved by PSO algorithm can maximize both torque and speed for the motor design problem. In the response surface method, the 3D view of mean response surface and 3D view of variance response surface are plotted in Figures 7 and 8. The maximal values of mean and variance can be observed in the figures. Also, the PSO process can be illustrated in Figure 9.
The final optimal solution is derived as . That is to say, the wiring dimension is set to 4 × 1.2 mm, stator wiring type is set to wyewye, and the airgap is set to 3.0 mm. The optimal solution can have better torque and speed performance. The optimal solution set can also be shown up in the combinational table.
The speed is 1120 rpm and the torque is 1.862 N·m. The optimal MPCI index value is obtained as 0.931. By averaging the 8th and 16th experimental data in combinational table, the optimal torque and speed can be observed in combination table.
8. Discussion
The optimal solution set is located at the (1.0, 1.0, 1.0). The optimal solution set means that the control factors with A_{3}, B_{3}, and C_{3} are optimal for the described problem. The resultant data for three studied cases are compared in Table 6. Fuzzy MPCI method is used to achieve the multiobjectives design. Both the torque and the speed keep at maximal values at the same time. Based on the derived multiobjective solution set, the maximal speed is 1119.5 rpm and the maximal force is 1.852 N·m. Also, in the combination table, the experimental speed 1120 rpm and experimental force 1.862 N·m can be found. Therefore, the optimal solution has good coincidence with the experimental data in combination table.

9. Conclusion
Maximal optimal values are obtained for both of the two multiobjective functions when the optimal control factor solution set is selected. The torque and speed characteristics keep at maximal values at the same time. Fuzzy MPCI method is used to perform the multiobjective design. Response surface method is used to derive the associated mathematical model for the studied system. It is very useful to squeeze the plastic material into molding area and provide enough squeezing pressure for the injection molding machine. The plastic material can be further condensed. The design purpose of turbo injection mode (TIM) for the axialflux motor is achieved. The improved TIM molding process can provide better product quality for the associated molding manufacturing process.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 A. Mebarki, K. Wejrzanowski, M. Shanel, and N. L. Brown, “A highpower, totally enclosed, permanent magnet, axial flux machine for engine integration,” in Proceedings of the 5th IET International Conference on Power Electronics, Machines and Drives, (PEMD '10), pp. 1–6, April 2010. View at: Publisher Site  Google Scholar
 A. M. ELRefaie, “Fractionalslot concentratedwindings synchronous permanent magnet machines: opportunities and challenges,” IEEE Transactions on Industrial Electronics, vol. 44, no. 1, pp. 107–121, 2010. View at: Publisher Site  Google Scholar
 W. L. Soong and T. J. E. Miller, “Fieldweakening performance of brushless synchronous AC motor drives,” IEE Proceedings: Electric Power Applications, vol. 141, no. 6, pp. 331–340, 1994. View at: Publisher Site  Google Scholar
 M. Aydin, S. Huang, and T. A. Lipo, “A new axial flux surface mounted permanent magnet machine capable of field control,” in Proceedings of the 37th IAS Annual Meeting of the Industry Applications Conference, vol. 2, pp. 1250–1257, Pittsburgh, Pa, USA, October 2002. View at: Publisher Site  Google Scholar
 A. Cavagnino, M. Lazzari, F. Profumo, and A. Tenconi, “Axial flux interior PM synchronous motor: parameters identification and steadystate performance measurements,” IEEE Transactions on Industry Applications, vol. 36, no. 6, pp. 1581–1588, 2000. View at: Publisher Site  Google Scholar
 T.S. Kwon, S.K. Sul, L. Alberti, and N. Bianchi, “Design and control of an axial flux machine for a wide fluxweakening operation region,” in Proceedings of the 42nd IAS Annual Meeting on Industry Applications Conference, pp. 2175–2182, New Orleans, La, USA, September 2007. View at: Publisher Site  Google Scholar
 A. M. ELRefaie, T. M. Jahns, P. J. McCleer, and J. W. McKeever, “Experimental verification of optimal flux weakening in surface PM machines using concentrated windings,” IEEE Transactions on Industry Applications, vol. 42, no. 2, pp. 443–453, 2006. View at: Publisher Site  Google Scholar
 J. L. Colton, Design of an integrated starteralternator for a series hybrid electric vehicle: a case study in axial flux permanent magnet machine design [A Dissertation Presented to the Faculty of the Graduate College at the University of Nebraska in Partial Fulfillment of Requirements for the Degree of Doctor of Philosophy], Electrical Engineering Theses and Dissertations, paper 12, pp. 27–31, 2010.
 J. Colton, D. Patterson, J. Hudgins, and K. Vacha, “Generator design for existent windmills: from water pumping to electricity generation,” in Proceedings of the IEEE Power Electronics and Machines in Wind Applications (PEMWA '09), June 2009. View at: Publisher Site  Google Scholar
 B. G. Wang, S. H. Wang, J. Q. Huang, M. Y. Cheng, and M. Q. Cai, “Light electric hub motor vehicle design with electronic variable speed control technology,” Electric Monthly, vol. 20, no. 8, pp. 140–148, 2010. View at: Google Scholar
 Z. Y. Lin, With reliability and highperformance electric vehicle drive system design and implementation [M.S. thesis], Department of Mechanical Engineering, National Cheng Kung University, 2009.
 B. Sun and Z. Wu, “Optimization of injection molding process parameters based on response surface methodology and genetic algorithm,” in Proceedings of the IEEE 2nd International Conference on Computer Engineering and Technology, pp. 397–400, 2010. View at: Google Scholar
 M. S. Yeoman, “Multiphase, dual polymer injection molding and cooling of an open cavity to form both distinct and graduated material properties within a complex threedimensional body,” in Proceedings of the COMSOL Conference, Bostan, Iran, 2009. View at: Google Scholar
 P. B. S. Reddy, K. Nishina, and A. S. Babu, “Taguchi's methodology for multiresponse optimization: a case study in the Indian plastics industry,” International Journal of Quality and Reliability Management, vol. 15, no. 6, pp. 646–668, 1998. View at: Publisher Site  Google Scholar
 C. W. Hung and C. Y. Huang, Optimization for plastic injection molding operations by dual response surface method and nonlinear programming [M.S. thesis], Department of Industrial Engineering, Tunghai University, 1998.
 Z.L. Gaing, “A particle swarm optimization approach for optimum design of PID controller in AVR system,” IEEE Transactions on Energy Conversion, vol. 19, no. 2, pp. 384–391, 2004. View at: Publisher Site  Google Scholar
 C. C. Hang, K. J. Astrom, and W. K. Ho, “Refinements of the ZieglerNichols tuning formula,” IEE Proceedings D: Control Theory and Applications, vol. 138, no. 2, pp. 111–118, 1991. View at: Publisher Site  Google Scholar
 Z. Malik and K. Rashid, “Comparison of optimization by response surface methodology with neurofuzzy methods,” IEEE Transactions on Magnetics, vol. 36, no. 1, pp. 241–257, 2000. View at: Publisher Site  Google Scholar
 H. J. Zimmermann, Fuzzy Set Theory and Its Applications, Kluwer Academic Publishers, London, UK, 1985.
 Z. L. Gaing and J. A. Chiang, “Robust design of inwheel PM motor by fuzzybased Taguchi method,” in Proceedings of the IEEE Power and Energy Society General Meeting, pp. 1–7, Minneapolis, Minn, USA, July 2010. View at: Publisher Site  Google Scholar
 Y. Shi and R. Eberhart, “A modified particle swarm optimization,” in Proceedings of the IEEE International Conference on Evolutionary Computation, pp. 69–72, May 1998. View at: Google Scholar
 R. C. Eberhart and Y. Shi, “Comparing inertia weights and constriction factors in particle swarm optimization,” in Proceedings of the Congress on Evolutionary Computation (CEC '00), pp. 84–88, July 2000. View at: Google Scholar
 S.Y. Ho, H.S. Lin, W.H. Liauh, and S.J. Ho, “OPSO: orthogonal particle swarm optimization and its application to task assignment problems,” IEEE Transactions on Systems, Man, and Cybernetics Part A:Systems and Humans, vol. 38, no. 2, pp. 288–298, 2008. View at: Publisher Site  Google Scholar
 J. Wang, S. Jin, S. Qin, and H. Jiang, “Swarm intelligencebased hybrid models for shortterm power load prediction,” Mathematical Problems in Engineering, vol. 2014, Article ID 712417, 17 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 J. Li, X. Wang, J. Jia, P. Wang, Y. Zhou, and Z. Zhao, “Location predictionbased data dissemination using swarm intelligence in opportunistic cognitive networks,” Mathematical Problems in Engineering, vol. 2014, Article ID 453564, 15 pages, 2014. View at: Publisher Site  Google Scholar
 L. Li, H. Huang, J. Lian et al., “Research of ant colony optimized adaptive control strategy for hybrid electric vehicle,” Mathematical Problems in Engineering, vol. 2014, Article ID 239130, 10 pages, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 JianLong Kuo and MengTi Chang. 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.