Research article 01 Jul 2021
Research article  01 Jul 2021
Design optimization of vehicle asynchronous motors based on fractional harmonic response analysis
 State Key Laboratory of Automotive, Jinlin University, Jinlin, 130025, PR China
 State Key Laboratory of Automotive, Jinlin University, Jinlin, 130025, PR China
Correspondence: ChuanXue Song (schuanxue_jlu@sina.com)
Hide author detailsCorrespondence: ChuanXue Song (schuanxue_jlu@sina.com)
To make vehicles more reliable and efficient, many researchers have tried to improve the rotor performance. Although certain achievements have been made, the previous finite element model did not reflect the historical process of the motor rotor well, and the rigidity and mass in rotor optimization are less discussed together. This paper firstly introduces fractional order into a finite element model to conduct the harmonic response analysis. Then, we propose an optimal design framework of a rotor. In the framework, objective functions of rigidity and mass are defined, and the relationship between high rigidity and the firstorder frequency is discussed. In order to find the optimal values, an accelerated optimization method based on response surface (ARSO) is proposed to find the suitable design parameters of rigidity and mass. Because the higher rigidity can be transformed into the firstorder natural frequency by objective function, this paper analyzes the firstorder frequency and mass of a motor rotor in the experiment. The results proved that not only is the fractional model effective, but also the ARSO can optimize the rotor structure. The firstorder natural frequency of asynchronous motor rotor is increased by 11.2 %, and the mass is reduced by 13.8 %, which can realize high stiffness and light mass of asynchronous motor rotors.
Electric vehicles appear to be one of the viable choices in face of the world's increasing attention to environmental protection, energy shortage and other issues (Emadi et al., 2008; Arhun et al., 2018a; Migal et al., 2019). Therefore, it is necessary to design and optimize vehicle monitors. Among the various kinds of electric monitors, induction motors have been widely used because of better performance (Dvadnenko et al., 2018; Francis et al., 2019; Hnatov et al., 2019; Zarma et al., 2019). High reliability and high efficiency are its main characteristics of induction motor (Benbouzid, 1999). For induction motors, it is reported that the failure rate is from stator, bearing, rotor and other aspects (O'Donnell, 2007; Benbouzid and Kliman, 2003; Yildirim et al., 2014; Albrecht et al., 1986).
Reducing the vibration effect of motors can improve the reliability, durability and service life of a motor. In order to achieve this, it is usual to analyze the natural frequency of the motor. The natural frequency can be obtained by numerical methods or modal test experiments (Ewins, 1985; Arhun et al., 2018b). Much research has been proposed in recent years. Ma et al. (2015) proposed an analysis method to study the vibration characteristics of an asymmetric and anisotropic rotorbearing system. By establishing a linear differential equation with periodic coefficients, the calculation efficiency was improved, and the results were also verified by experimental studies. They suggested to investigate the frequency characteristics and stability of asymmetric anisotropic rotorbearing systems. Modifications are made to incorporate the effect of stator asymmetry into an existing threedimensional (3D) solid finite element procedure. In the research of Yin et al. (2020), a specially designed modal conversion horn with an oblique beam was applied. They transferred the longitudinal vibration to an additional bending vibration. Modal analysis was conducted to illustrate the elliptical vibration process by the finite element method. Then the output performance of the motor was evaluated via a series of experiments. Widdle et al. (2006) proposed a model for highfrequency torsional vibration analysis. The research of Li et al. (2008) studied the nonlinear vibration of a threephase AC motorlinkage mechanism system with links fabricated from threedimensional braided composite materials. Migal et al. (2021) described a method to select an asynchronous traction electric drive for an electric vehicle that enables the assessment of necessary technical, environmental and operational qualities.
In general, the prior studies have made progress on vibration modal analysis and design optimization of motor, but there is still room for improvement. Firstly, the built finite element model does not consider the historical process of the motor rotor well. Therefore, it cannot accurately reflect the actual dynamic characteristics of the motor. Moreover, most of the existing methods discussed the motor performance by considering frequency, but fewer involved the research of considering the mass and rigidity together. To improve the above situations, this paper combines the fractional harmonic response (Wang and Jiang, 2018; Yan et al., 2020) to the finite element model that can have the dynamic characteristic. Then, we set the high rigidity and the light mass as the optimization objectives. By analyzing the relationship between rigidity and firstorder frequency, we transform the optimization into the firstorder frequency and the mass. In order to find optimal values, an accelerated optimization method based on response surface (ARSO) is proposed in this paper, and three intelligent optimization algorithms – including traversal search algorithm (TS) (Fang and Xu, 2017), multiobjective genetic algorithm (MOGA) (Ponnambalam et al., 2000) and optimization algorithm based on response surface (RSO) (Munck et al., 2008) – are selected for comparison. Finally, experimental results are given to prove the effectiveness of our proposed approach.
The paper is organized as follows: Sect. 2 presents the fractional finite element model of rotor support casing. In Sect. 3, the design optimization of asynchronous motor rotor is described. Section 4 shows the experimental results. Section 5 is optimization analysis, and discussion is given in Sect. 6. Section 7 concludes this paper.
2.1 Definition of fractional derivative
As defined in the work by Scherer et al. (2011), the Grünwald–Letnikov definition shows an expression for qth derivative, where $q\in {R}^{+}$. In addition, it allows one to consider the socalled shortmemory principle as follows.
where $(\mathrm{1}{)}^{j}{(}_{j}^{q})$ represents binomial coefficients, ${c}_{j}^{q}$ for $j=\mathrm{0},\mathrm{1},\mathrm{2},\mathrm{\dots}$, which are given by
Based on Eq. (1), the solution of a fractional differential equation given by ${}_{a}{D}_{t}^{q}y\left(t\right)=f\left(y\right(t),t)$ is defined as $y\left({t}_{k}\right)=f\left(y\right({t}_{k}),{t}_{k}){h}^{q}{\sum}_{j=v}^{k}{c}_{j}^{q}y\left({t}_{kj}\right)$, where t_{k}=kh, h is the time step, and v follows the rule
with L_{m} being the memory length, which can be set according to the required accuracy.
2.2 Model description
Figure 1 describes the finite element rotor model. The rotor can be divided into common beam elements, including elastic shaft, distributed mass and stiff plate. N nodes and M plates are contained in the rotor. As seen in Fig. 1, the elastic modulus, area moment of inertia, shear modulus, Poisson ratio, shaft length, shaft density and shaft crosssectional area are represented by E, I, G, μ, L, ρ and A, respectively. For disc P_{i}, J_{ddi}, m_{rpi} and J_{pdi} are equatorial moment of inertia, disc mass and pole moment of inertia, respectively. F_{xi} and F_{yi} represent the force of the i node, and the moments of the i node are described by M_{xi} and M_{yi}. In the coordinate system, the position of the elastic center line is decided by the displacements of x(s,t) and y(s,t); (s,t) angles of x and y can be used to obtain the orientation of the crosssection.
Now the equation of rigid motion is introduced. m_{p}, J_{dd}, J_{pd} and ω are set to represent the plate mass, moment of inertia of the Equator, the moment of inertia of the pole and rotational angular velocity of plate, respectively. By Lagrangian theorem, the equation of rigid motion is given as follows.
where M_{td} and M_{rd} are mass matrix and mass inertia matrix. G_{d} is the gyro matrix, and ${q}_{\mathrm{d}}=[x,y,\mathit{\varphi},\mathit{\psi}]$ is the generalized displacement vector. Q_{d} is generalized external force vector. Then, M_{td}, M_{rd} and G_{d} are given:
Next, the equation of beam motion is described secondly. E, G, μ, d, D and L are set as the element modulus of elasticity, the shear modulus, the Poisson ratio, the inner diameter, the outer diameter and the length, and then I is calculated:
The crosssectional area A
The effective shear area A_{s}
or
Shear deformation coefficient ϕ_{s}
For a beam element, two nodes, 8 degrees of freedom, and 8 degrees of freedom are contained. Then the crosssectional displacement of the unit as a function of time can be seen as a function of the position along the axis of the unit. Therefore, the generalized displacement of the unit endpoint with time change can be described by
According to Lagrange theorem, the equation of motion of the beam is as follows:
where Q_{e} is the generalized external force vector, M_{te} and M_{re} are mass matrix and mass inertia matrix, G_{e} is gyro matrix, K_{be} is shearing stiffness matrix, and K_{ae} is the unit tensile stiffness matrix. The corresponding calculation process is given in Eqs. (19)–(24).
Finally, the equation of rotor motion is obtained by the equations of rigid motion and beam as follows:
where M_{s} is the system mass matrix, Q_{s} is the system generalized external force vector, G_{s} is the system gyro matrix, C_{s} is the system damping matrix and K_{s} is the system stiffness matrix.
2.3 Fractional harmonic response analysis
After the finite element rotor is built, now we conduct fractional harmonic response analysis. Here we assume the Q_{s} and C_{s}−ωG_{s} in Eq. (25) to be zero, and the free vibration can be expressed as the generalized eigenproblem,
where ω is the natural frequency and X the mode of vibration. By Eqs. (8)–(24), the global stiffness matrix K_{s} and mass matrix M_{s} are calculated. Then by Eq. (26), the natural frequency and mode of plate vibration can be obtained.
Next, the harmonic response analysis is used to solve the vibration problem. It is a method to determine the structural response of solid materials towards loads that change according to time and can predict the sustained dynamic performance of the system structure by the frequency response curve. Actually, the essence of harmonic response analysis is to solve the forced vibration equation of structure. Here, the equation of forced vibration is given as follows.
where F(t) is force vector matrix and F_{0} is the force load amplitude.
In the view of mechanical engineering, a rotor can be seen as an energy dissipation system, and the force relies on the operational process of the system. As known, the fractional derivative does not depend on the discrete points of fractional order finite element model, but it does depend on the historical process. Due to the special nature of fractional derivative, it has been successfully used to describe the viscoelastic characteristic of all kinds of systems (Yan et al., 2020). Thus, this paper also considers the fractional derivative (𝒟^{α}(*)) into the modeling of force. Since the fractional order usually acts on a damping item for a viscoelastic model, a fractional force model can be written as
No matter what finite element model is established, it is still essential to discuss the design optimization of the rotor because the dynamic stiffness and mass of a rotor are the important factors that affect its working characteristics and load efficiency. The firstorder natural frequency should be made as high as possible to effectively avoid the resonance region, and then the high rigidity of the asynchronous motor is required. While the rigidity of the asynchronous motor is increased, the mass can change. Therefore, it is necessary to consider rigidity and mass for the design of the asynchronous motor rotor.
3.1 Optimal design framework
The optimal design is based on the built fractional order finite element model. Figure 2 describes the optimal progress. Firstly, the geometric parameters of the motor rotor are used as design variables. Secondly, the maximized firstorder natural frequency and the lightest quality of the monitor are set as the optimization objective. Thirdly, the rotor structure of the motor is confined to a certain space. Finally, we apply intelligent algorithms to optimize the geometric parameters of the asynchronous motor rotor and perform simulation analysis on the optimized model. The simulation results will be compared with the previous rotor parameters to verify the optimal progress.
The establishment of the correct mathematical model is the key of the optimal design and the relevant parameters that are needed to discuss in the optimization process. The detailed process is given as follows:

determine the parameter type, initial value and variable range of the optimized design;

establish a parameterized model including design variables and solution parameters;

determine the objective function and its mathematical expression based on the optimization objective;

define the range of the rotor to ensure the movement range of the design point can be controlled within the feasible region;

check the rationality of the optimization design model to improve the efficiency and stability;

set the initial value, upper and lower bounds of the variables, optimize operating parameters, and complete computer programming.
3.1.1 Design variables
This paper uses ${x}_{i}(i=\mathrm{1},\mathrm{2},\mathrm{\dots},n)$ to present the deign variables, and they can be shown by a matrix as $x=[{x}_{\mathrm{1}},{x}_{\mathrm{2}},\mathrm{\dots},{x}_{n}{]}^{\mathrm{T}}$, where ${x}_{\mathrm{1}},{x}_{\mathrm{2}},\mathrm{\dots},{x}_{n}$ are n components of the vector x. And then R^{n} is an ndimensional European space, which can be generated by the design variables.
The function of the design variable is also called the state variable, which is used as the constraint of the design value. In this paper, the rotor mass (m), the firstorder frequency (f_{I}) and the maximum stress (σ) are selected as state variables, which can be presented as follows:
According to the definition, Fig. 3 shows the design variables including the shaft extension L_{1}, the shaft segment length L_{2} and the shaft diameter D. Then the optimization design problem of the motor rotor can be seen as an optimization problem of mathematical model.
3.1.2 Objective function
Objective function can be expressed as F(x). For the asynchronous motor, the rotor optimization design requires high rigidity and light mass, which can be recognized as multiobjective optimization problem:

High rigidity. If the high rigidity of the asynchronous motor rotor is the objective function, the dynamic stiffness of the rotor should be increased as much as possible. Due to the inevitable error during the assembling process of the asynchronous motor, it can cause an imbalance that is dynamically affected by the mass. Therefore, the dynamic stiffness of the rotor must be sufficiently large. We describe the stiffness by the vibration amplitude of the rotor under dynamic excitation, and the dynamic stiffness (K) is $K={F}_{\mathrm{e}}/A$, where F_{e} is dynamic excitation load and A is the corresponding displacement amplitude.
In reality, the dynamic excitation of the rotor is constantly changing. When the changing frequency is far from the first natural frequency of the rotor, the dynamic stiffness and static stiffness are basically the same; when the changing frequency of the dynamic excitation is close to the first natural frequency of the rotor, the corresponding displacement amplitude will increase sharply and resonance will appear. Therefore, in order to achieve the high stiffness of the rotor, the firstorder natural frequency of the rotor needed to be increased. So the relevant structural parameters will be changed, which make the firstorder natural frequency and the dynamic stiffness increase simultaneously. Therefore, the highstiffness problem can be transformed into the firstorder frequency solution.

Light mass. The mass of the rotor directly affects its dynamic unbalance state under highspeed operation. Therefore, it is required to reduce the rotor mass as much as possible. In this paper, the mass level is divided into higher, default and lower, which represents the priority and importance of the optimal solution. Therefore, the highrigidity target is determined as the higher level, and the light mass is positioned at original level.
3.1.3 Constraint condition
Besides design variables and objective functions, constraint condition is required, and it can be generally divided into equality constraint (Eq. 31) and inequality constraint (Eq. 32).
or
where g(x) and h(x) are the objective function with ndimensional vectors. According to the above definition, the constraint condition in the optimization design of asynchronous motor rotor can be defined that the maximum stress of the shaft under the torque does not exceed the allowable stress of the shaft material:
3.2 Optimization algorithms
In order to find the suitable parameters for the design of the asynchronous motor rotor, this section described a improved algorithm based on response surface optimization(RSO) and named it as accelerated response surface optimization (ARSO). Moreover, three algorithms – including traversal search algorithm (TS), multiobjective genetic algorithm (MOGA) and RSO – are chosen to evaluate the effectiveness of ARSO. The following sections mainly describe ARSO, but the other three algorithms (Fang and Xu, 2017; Ponnambalam et al., 2000; Munck et al., 2008) are also briefly described.
3.2.1 Accelerated optimization algorithm based on response surface (ARSO)
In order to improve poor generalization ability of traditional RSO (Ponnambalam et al., 2000), the least squares support vector machine (LSSVM) provides a new way of structural reliability analysis (Deng et al., 2003; Cao et al., 2014) on RSO. Firstly, the inequality constraint function of SVM is converted into equality constraint function:
And constraint condition
where ω is the weight vector, C is the penalty factor, e_{i} is the error, ϕ is the nonlinear function of lowdimensional to highdimensional mapping, and b is constant coefficient.
Then the squared error can be weighted to promote robustness:
And constraint condition
where v_{i} is the weight coefficient. Finally, the new regression prediction function can be obtained by Lagrangian polynomials:
where α_{i} is Lagrange multiplier, and K(x,x_{i}) is kernel function.
It is timeconsuming based on the above calculation. The reason is when (K(x,x_{i})) is used to construct shape functions, an n×n linear system should be computed for every computational point. If m monomials are added, an $(n+m)\times (n+m)$ linear system should be solved. Therefore, this paper proposes an accelerated optimization based on the previous works and name it ARSO.
For the LSSVM based on RSO, kernel function includes [K^{T}(x)p^{T}(x)] and K^{−1}. The time consumption of [K^{T}(x)p^{T}(x)] is less than that of K^{−1} because the computational complexities are o(n+m) and o((n+m)^{3}), respectively. Every computational point has its own [K^{T}(x)p^{T}(x)] and is different from the rest, but it may have the same K^{−1} as the other. Thus, we need a storage place, containing the public nodes and the corresponding points' information.
First of all, we construct an N_{com}×N_{cor} matrix S=[s_{ij}]. N_{com}×N_{cor} denotes the number of the computational nodes and the number of the corresponding points, respectively. The element S=[s_{ij}] is a δ function.
If the jth row elements are the same as the kth row elements, we confirm the jth and kth Gauss points have the same kernel functions K(x_{j}) K(x_{k}). It is a rule, which gives onetoone mappings between the two type sets by comparing the rows of S. During the calculation process, the improved method can preset the search area to reduce the number of repeated matrix inversions. To ensure the algorithm's effectiveness, the main parameter settings are shown in Table 1.
3.2.2 Three compared algorithms
We select three traditional algorithms to compare. These are traversal search algorithm (TS), multiobjective genetic algorithm (MOGA) and optimization algorithm based on response surface (RSO). TS is to randomly generate a certain number of sample points according to the design space of the independent variable and calculate the generated sample points one by one. MOGA is an algorithm with the characteristics of global random searching and implicit parallel searching. It simulates the problemsolving process as the genetic evolution process of biological populations to find the optimal solution. RSO is an optimization method based on limited design space. Firstly, a polynomial method is used to fit a response surface that is shown in Fig. 4; it is the response surface fitted by the quadratic polynomial. Secondly, all the points on the response surface form a subset of the design space, which reduces the design space of the entire optimized design and improves the computational efficiency.
3.3 Optimization processing
3.3.1 Optimal solution
The optimization design problem includes establishing the model and solving the model. Solving the model is to obtain the optimal solution. The next step is to analyze and judge the practicality of the optimal solution and finally determine the optimal design plan.
Generally, the n design variables can make the objective function reach the extreme value under the restriction of constraints:
When the constraint function is a nonconvex set or the objective function is a nonconvex function, there may be multiple optimal solutions. The optimal solution obtained is called a local optimal solution. Only when the set of constraints is a convex set or the objective function is a unimodal function in the domain can the obtained local optimal solution be judged as the global optimal solution.
3.3.2 Modify the original model parameters
Once the final optimization plan is determined, the original model can be modified according to the parameter values of the optimal solution, which could improve product performance and reduce manufacturing costs.
3.3.3 Simulation analysis of optimization scheme
During the initial establishment of the optimal design model, due to the complexity of the actual problem and the limitations of the mathematical model description, some design parameters and constraints were simplified. After the final optimization plan is determined, the simplified content should be added to the fractional finite element model again, and then the static and dynamic characteristics that are closer to the actual geometric model can be simulated.
In order to verify the correctness of the established fractional model, it is necessary to explore whether its dynamic response is consistent with the actual asynchronous motor. Therefore, a modal test of the asynchronous motor percussion method is designed to measure the natural frequency of the asynchronous motor in a free state. In the experiment, the asynchronous motor was suspended with a flexible rope to make it in a free state. The experimental instrument relied on the German Plufor vibration analyzer Vib XpertII. By pasting an acceleration sensor on one side of the GM7101 motor, it was hit with a hammer on the opposite side. The vibration response of the asynchronous motor can be measured by the hammering experiment module in Vib XpertII. The suspension mode and sensor arrangement of asynchronous motor are used in The State Key Laboratory of Automotive Simulation and Control (ASCL). The GM7101 test site is shown in Fig. 5. Then the modal test results are given in Table 2, and the comparison between simulation results and experimental results is given in Table 3.
In Table 3, it can be seen that the firstorder frequency obtained by the simulation calculation is slightly lower than the experimental results, and the error may be caused by the following:

When the threedimensional model and finite element model are established, many simplified principles are adopted, and many small chamfers, threads, ventilation ducts and other structures are ignored.

The influence of contact stiffness and damping factors of the contact surface is ignored.

Silicon steel entities are used to bringing inevitable errors.

The actual bearing has a clearance. The clearance of the bearing is ignored when the spring structure is used for equivalent replacement, which causes a deviation in the modal analysis.

The experimental instrument and its installation may bring a few errors.

The boundary condition of the simulation analysis is a completely free state, but the suspension method used in the experiment is not easy to achieve the completely free state.
This section introduces the optimized settings and analyzes the optimal results.
5.1 Optimization settings
5.1.1 Design variables
As shown in Table 4, three design variables are defined respectively: the shaft extension L_{1}, the shaft segment length L_{2} and the shaft diameter D.
5.1.2 State variables
There are three state variables for the rotor optimization design problem of asynchronous motors: the mass of the rotor m, the firstorder natural frequency f_{1} and the maximum stress σ.
5.1.3 Constraint conditions
The maximum bending stress does not exceed the allowable shear stress, which is σ≤50 MPa.
5.1.4 Objective function
The main goal of the rotor optimization design of the asynchronous motor is to improve the high rigidity of the rotor but to minimize the rotor mass. Therefore, there are two objective functions: the firstorder natural frequency of the rotor should be as large as possible when it exceeds the frequency corresponding to the maximum speed; the rotor mass is as small as possible under 10 % reduction of the mass. The relevant settings of state variables, constraint conditions and objective function are described in Table 5.
5.2 Optimization results
As seen in Table 6, it can be found that the optimization results obtained by the four algorithms are all acceptable. The firstorder natural frequencies obtained by TS and MOGA are both around 150.4 Hz, which is smaller than 152 Hz, which was obtained by the RSO and ARSO. Meanwhile, there are no large differences in the mass obtained by the four algorithms, but the calculation time of ARSO is the least compared to the other three algorithms. On the whole, the proposed ARSO performs better than other algorithms, and the firstorder natural frequency and the lightest mass are 152.17 Hz and 7.801 kg, respectively.
6.1 Model validation
According to Eq. (28), modal superposition method is used to discuss the harmonic response and obtain the modal solution and mode shape. The frequency range is set to be 0–200 Hz, and the number of operations is 500; that is, the interval of each solution is 0.5 Hz. And the simple harmonic exciting force is added to two vertical directions of the shaft section on the shaft flywheel rotor. The amplitude is 0.001 N, and the phase is 0 and 90^{∘}.
The flying wheel is selected as the response surface, and the maximum displacement response is set as the vertical coordinate and the excitation frequency is set as the horizontal coordinate. As seen in Fig. 6, the harmonic response simulation curve is obtained. By analyzing the displacement response frequency curve, the peak can be found at 123.5 Hz, which is also the distribution position of the first two natural frequencies of the system. At this time, the system resonates and the response increases sharply, which fulfill the dynamic characteristics of the system and also prove the accuracy of the natural frequency calculation. Meanwhile, it can also be seen that the system has a greater vibration response when the rotor is operating in the firstorder natural frequency range. After the harmonic response curve crosses the firstorder natural frequency, it shows a downward trend with the increase of the exciting force frequency; that is, the disturbance vibration caused by the unbalanced force attenuates with the increase of the frequency. Then, the system enters the normal working range, and its stability is gradually improved. Moreover, the peak value of the displacement response of the shafting rotor is about 10^{−7} and extremely small. It indicates that the shafting is running well and there will be no collision between the bearing rotor and the stator.
According to the test results of 129 Hz in Tables 2 and 3, it can known that the response curve can effectively represent the real frequency distribution, especially the first natural frequency.
6.2 Design variables discussion
In the optimal phase, it is necessary to choose suitable design variables. This section discusses the relationship among those variables and obtains the optimal values.
According to the relationship diagram shown in Fig. 7, the first three subfeatures show that L_{1} and L_{2} are not linear with frequency especially for L_{1}. On the contrary, it can be found that natural frequency is linear with diameter in the third subfigure, where it is consistent with the theoretical analysis. It is obvious that there is always a linear relationship between mass and length (L_{1} or L_{2}) from the last three pictures.
In Table 7, the optimal solution is rounded by being combined with the rationality requirements of the actual structure of the rotor. The L_{1} and L_{2} are close to each other when the shaft diameter becomes smaller. The numerical results show that the longer the rotor length is, the higher the natural frequency is. In terms of mass, a smaller diameter can not only reduce the quality because the mass is proportional to the square of diameter, but also effectively offset the influence of length.
According to the above discussion, Table 8 gives the comparison between the original values and the optimized results. It can be seen that the firstorder natural frequency is increased by 11.2 %, and the rotor mass is reduced by 13.8 %.
In this paper, the fractional model of the asynchronous motor rotor was firstly established with a peculiar memory characteristic, and the introduced harmonic response was able to fit the reality well. Then, we set high rigidity and less mass as optimization functions and transform them into the problem of the firstorder frequency and mass. In order to find the optimal parameters, an accelerated optimization method based on response surface is proposed. Finally, the GM7101 asynchronous motor was selected to evaluate our work. The experimental results show that the optimized firstorder natural frequency was increased by 11.2 %, and the mass was reduced by 13.8 %, where they verified the effectiveness and correctness of the fractional order and harmonic response analysis and realized the high rigidity and light mass of the asynchronous motor rotor.
However, an asynchronous motor is a complex system, and there are many factors that affect its vibration characteristics, such as frame type, ribs, ventilation ducts and so on. To study more complex asynchronous motors, more influencing factors will be considered in the future.
The data and code in this study can be obtained from the opensource platform Gitee (https://gitee.com/leiaojlu/designoptimization/tree/master, Lei, 2021).
AL was responsible for the conceptualization, methodology and writing of the paper draft. CXS and YLL were responsible for supervision. YF was responsible for editing.
The authors declare that they have no conflict of interest.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This paper was edited by Dario Richiedei and reviewed by Yana Rahun and one anonymous referee.
Albrecht, P. F., Appiarius, J. C., McCoy, R. M., Owen, E. L., and Sharma, D. K.: Assessment of the reliability of motors in utility applications – updated, IEEE T. Energy Conver., 1, 39–46, 1986. a
Arhun, S., Migal, V., Hnatov, A., Hnatova, H., and Ulyanets, O.: System approach to the evaluation of the traction electric motor quality, EAI Endorsed Transactions on Energy Web, 7, 162733, https://doi.org/10.4108/eai.1372018.162733, 2018a. a
Arhun, S., Migal, V., Hnatov, A., Ponikarovska, S., and Novichonok, S.: Determining the quality of electric motors by vibrodiagnostic characteristics, EAI Endorsed Transactions on Energy Web, 7, 164101, https://doi.org/10.4108/eai.1372018.164101, 2018b. a
Benbouzid, M. and Kliman, G.: What stator current processingbased technique to use for induction motor rotor faults diagnosis?, IEEE T. Energy Conver., 22, 62–62, 2003. a
Benbouzid, M.: Bibliography on induction motors faults detection and diagnosis, IEEE T. Energy Conver., 14, 1065–1074, 1999. a
Cao, J., Ding, W. Y., Zhao, D. S., Song, Z. G., and Yue, G.: Optimization design of supporting structure in foundation pit based on response surface method of LSSVM, Journal of Kunming University of Science and Technology (Natural Science Edition), 39, 43–47, 2014. a
Deng, J., Yue, Z. Q., Tham, L. G., and Zhu, H. H.: Pillar design by combining finite element methods, neural networks and reliability: a case study of the Feng Huangshan copper mine, China, Int. J. Rock Mech. Min., 40, 585–599, 2003. a
Dvadnenko, V., Arhun, S., Bogajevskiy, A., and Ponikarovska, S.: Improvement of economic and ecological characteristics of a car with a startstop system, International Journal of Electric and Hybrid Vehicles, 1, 209–222, 2018. a
Ewins, D. J.: Modal testing: theory and practice, Wiley, Chichester, England, 1985. a
Emadi, A., Lee, Y. J., and Rajashekara, K.: Power electronics and motor drives in electric, hybrid electric, and plugin hybrid electric vehicles, IEEE T. Ind. Electron., 55, 2237–2245, 2008. a
Francis, J., Biju, N. A., Johnson, A., Mathew, J., and Sankar, V.: Comparison of sixphase and threephase induction motors for electric vehicle propulsion as an improvement toward sustainable transportation: proceedings of GBSE 2018, Springer, Singapore, 2019. a
Fang, J. and Xu, Y.: Optimization of production line equilibrium based on traversal search algorithm and genetic algorithm, Computer Applications and Software, 34, 276–280, 2017. a, b
Hnatov, A., Arhun, S., Tarasov, K., Hnatova, H., and Patlins, A.: Researching the model of electric propulsion system for bus using matlab simulink, 2019 IEEE 60th International Scientific Conference on Power and Electrical Engineering of Riga Technical University (RTUCON), IEEE, Riga, Latvia, 2019. a
Lei, A.: Design optimization [data set and code], https://gitee.com/leiaojlu/designoptimization/tree/master, 2021. a
Li, Z., Cai, G., Huang, Q., and Liu, S.: Analysis of nonlinear vibration of a motorlinkage mechanism system with composite links, J. Sound Vib., 311, 924–940, 2008. a
Ma, W. M., Wang, J. J., and Wang, Z.: Frequency and stability analysis method of asymmetric anisotropic rotorbearing system based on threedimensional solid finite element method, J. Eng. Gas Turb. Power, 137, 102502, https://doi.org/10.1115/1.4029968, 2015. a
Migal, V., Arhun, S., Hnatov, A., Dvadnenko, V., and Ponikarovska, S.: Substantiating the criteria for assessing the quality of asynchronous traction electric motors in electric vehicles and hybrid cars, Journal of the Korean Society for Precision Engineering, 36, 989–999, 2019. a
Migal, V., Lebedev, A., Shuliak, M., Kalinin, E., and Korohodskyi, V.: Reducing the vibration of bearing units of electric vehicle asynchronous traction motors, J. Vib. Control, 27, 1123–1131, https://doi.org/10.1177/1077546320937634 2021. a
Munck, M. D., Moens, D., Desmet, W., and Vandepitte, D.: A response surface based optimisation algorithm for the calculation of fuzzy envelope FRFs of models with uncertain properties, Comput. Struct., 86, 1080–1092, 2008. a, b
O'Donnell, P.: Report of large motor reliability survey of industrial and commercial installations, Part II, IEEE T. Ind. Appl., 21, 865–872, 2007. a
Ponnambalam, S. G., Aravindan, P., and Naidu, G. M.: A multiobjective genetic algorithm for solving assembly line balancing problem, Int. J. Adv. Manuf. Tech., 16, 341–352, 2000. a, b, c
Scherer, R., Kalla, S. L., Tang, Y. F., and Huang, J. F.: The Grünwald–Letnikov method for fractional differential equations, Comput. Math. Appl., 62, 902–917, 2011. a
Wang, N. and Jiang, D.: Vibration response characteristics of a dualrotor with unbalancemisalignment coupling faults: Theoretical analysis and experimental study, Mech. Mach. Theory, 125, 207–219, 2018. a
Widdle, R. D., Krousgrill, C. M., and Sudhoff, S. D.: An induction motor model for highfrequency torsional vibration analysis, J. Sound Vib., 290, 865–881, 2006. a
Yan, D., Wang, W., and Chen, Q.: Fractionalorder modeling and nonlinear dynamic analyses of the rotorbearingseal system, Chaos Soliton. Fract., 133, 109640, https://doi.org/10.1016/j.chaos.2020.109640, 2020. a, b
Yildirim, M., Polat, M., and Kurum, H.: A survey on comparison of electric motor types and drives used for electric vehicles, 2014 16th International Power Electronics and Motion Control Conference (PEMC), IEEE, Antalya, Turkey, 2014. a
Yin, Z., Dai, C., Cao, Z., Li, W., and Li, C.: Modal analysis and moving performance of a singlemode linear ultrasonic motor, Ultrasonics, 108, 106216, https://doi.org/10.1016/j.ultras.2020.106216, 2020. a
Zarma, T. A., Galadima, A. A., and Aminu, M. A.: Review of motors for electrical vehicles, Journal of Scientific Research and Reports, 1–6, 2019. a