Abstract

Asynchronous motor system has the characteristics of high order, strong coupling, and nonlinearity. From the dynamical model, it is the underactuated mechanical system, which means that the dimension of its input space is fewer than the degree of freedom. Following this perspective, the energy based nonlinear control technology-CL (controlled Lagrangians) method is used to solve the control problem in this paper. Based on the expected controlled energy and its derivative with respect to time, controlled Lagrangians and generalized force are constructed, and they produce the controlled equations. In order to ensure the complete matching between the controlled equation and the original equation, the gyroscopic forces containing the first-order term of velocity are innovatively introduced into the generalized force, and the matching conditions are obtained. By solving the matching conditions composed of some partial differential equations, the nonlinear smooth feedback control law can realize the global asymptotic stabilization of not only velocity but also position. Finally, the controlled energy is selected as the Lyapunov function, and the stability is proved according to the LaSalle invariant theorem. The effectiveness of the designed control law is demonstrated in the results of the simulation.

1. Introduction

With the advantages of low price, simple structure, convenient maintenance, and reliable operation, the asynchronous motor has always been in a leading position in today’s social industrial production. Under the concept of advocating production environmental protection and low-carbon economy, the research on the control performance of the asynchronous motor has important theoretical significance and practical value [1].

The asynchronous motor is a nonlinear system. However, the traditional linear control method cannot reveal its nonlinear nature. Therefore, the research on the control method of nonlinear theory is of great significance to improve the dynamic and static performance of the AC asynchronous motor. At present, the nonlinear control methods applied to the asynchronous motor mainly include feedback linearization control [2], backstepping control [3], sliding mode control [4], active disturbance rejection control [5, 6], and passive control theory [7, 8]. The control performance of the system has been significantly improved for application of the above method.

Sun developed chopping control and energy-saving controller of a three-phase AC asynchronous motor [2]. Yu et al. designed the nonlinear adaptive controller of the asynchronous motor system by using the subsystem separation method and backstepping technology to ensure the stability of the system [3]. Lekhchine et al. designed a renewable energy storage electrical system for asynchronous motors. In this system, the motor is driven by sliding mode control, which can overcome the chattering phenomenon through the sliding surface based on fuzzy logic [4]. Li et al. proposed a second-order ADRC and AC excitation control system based on stator flux oriented control to control the active and reactive power of the variable-speed pumped storage unit [5, 6]. Wu et al. discussed the problem of asynchronous passive control of Markov jump systems and obtained three equivalent sufficient conditions to ensure the random passivity of hidden Markov jump systems by using matrix inequality technology. Based on the established conditions, an asynchronous controller is designed [79]. Yu et al. studied the tracking control of the underactuated dynamic system and proposed a six-step motion strategy of the pendulum driven vehicle rod system [10, 11]. By implementing feedback and other measures, the asymptotic stability of the control Hamiltonian system can be realized. The literature [1215] reports some new results about the control of underactuated dynamic systems.

In terms of mathematical equivalence, some studies analyze the controlled Lagrange (CL) method [1619]. Compared with PBC, the CL method has a simpler mathematical form and clearer physical meaning, which is easy to understand. Usoro et al. described a Lagrangian method for solving nonlinear constrained optimization problems in set theory control problems. By introducing the matrix Lagrange multiplier, the problem is simplified to solve a set of nonlinear simultaneous matrix equations [16]. Műller et al. proved the possibility of maximizing the torque without exceeding the limit value of magnetic flux and stator current, which is independent of the number of revolutions of an asynchronous motor [17]. Lindgren et al. gave the exact slope distribution and other characteristic distributions of symmetric and asymmetric Lagrangian spatiotemporal waves at level crossings [18, 19]. In addition, some studies are extended to the general PBC method from the perspective of robust control and optimal control of general PCH systems [2026]. These results have been proved to be expressed in a Lagrangian form.

This paper applys CL method to analyze the asynchronous motor system from the perspective of an underdrive mechanical system. We will use the electromagnetic energy generated by the stator and rotor windings and the mechanical energy generated by the rotor to construct a controlled energy controller [12]. The controlled system maintains the Lagrangian mechanical structure in form, obtains the smooth nonlinear feedback control law, has a large convergence range, and helps to realize the CL method robust control and optimal control [14]. Compared with the port controlled dissipative Hamiltonian system, the nonlinear smooth feedback control law obtained in this paper can realize the global asymptotic stabilization of position and velocity at the same time.

2. Mathematical Model

For the convenience of writing, we will indicate the independent variables of the functions and matrices that appear, which will be omitted when they appear below. , and , means a collection consisting of the first n quantity of natural numbers. Let represents the function of the vector , corresponds to the element at the kth row and the jth column of the function matrix vector , denotes the ith element of function vector : , where , means the five-order identity matrix. Besides, some notations as given below:

The generalized coordinate variables of the AC asynchronous motor is , where are the components of stator inductance charge and rotor inductance charge on axis, and is the angular position of motor rotor. is the original control input, where are the elements of the stator voltage and rotor voltage on the axis. In addition, , where the input coupling matrix and .

In view of the mathematic model of a three-phase asynchronous motor, the following assumptions are expressed as follows [9]:(1)The spatial harmonic and the spatial difference between three-phase windings are ignored. Meanwhile, it is assumed that the generated magnetomotive force is distributed sinusoidally along the circumference of the air gap.(2)Magnetic circuit saturation and core loss are ignored. At the same time, it is assumed that the inductance parameters of every phase winding, not only self-inductance but also mutual one, are constant.(3)It is not considered of the influence of frequency and temperature changes on the variation of winding resistance.

According to the above assumptions for the AC asynchronous motor, its mathematical model in the coordinate system can be obtained by

In (2), , , L11, and L33 denote the equivalent self-inductance of the stator and rotor phase windings, respectively. And L13 is the equivalent mutual inductance of the stator and rotor phase windings. The load torque is denoted by TL, and , where includes no-load torque and the external one, and denotes the torsional torque generated when the motor and the mechanical load are connected with a relatively long shaft, where is the deformation coefficient. The equation (2) is abbreviated aswhere is the equilibrium point, and its input must satisfy , where and , so .

3. Design of Asynchronous Motor Controller Based on CL Method

3.1. Construction of Controlled Energy and Generalized Force

According to (3), the controlled kinetic energy of the controlled system is . The controlled kinetic energy matrix satisfies and .

Take the controlled potential energy as , and generalized force , then controlled Lagrangian function and controlled energy of the system are as follows:

Sometimes, the controlled energy has physical meaning, such as controlled kinetic energy or controlled potential energy, and maybe it has only a mathematical meaning which is sufficient and necessary. According to and , we obtain the controlled equations of the system as follows:

Using (4) and (5), we get

The generalized force of the system consists of two parts, namely gyroscopic forces and dissipation force , where gyroscopic forces matrix is .

Furthermore, represents the element at the th row and th column of function matrix , and for , is the th component function of the element of the gyroscopic forces matrix.

Since the gyroscopic forces matrix is an anti-symmetric one, there is . Similarly, the constant elements of the gyroscopic forces matrix also satisfy .

Remark 1. Since there is term in the original system, constant term is introduced into the matrix to construct the gyroscopic forces consisting of the first term of the velocity, which matches of the original system. Maybe these forces do not exist in real world, and only mathematical meaning is necessary for them. As we know, this introduction is for the first time.
When generalized force , we obtain from (5) thatwhich indicates that the energy of the closed system is decreasing. Multiplying −N(q) at the two sides of (4) at the same time, we obtainAccording to (2) and (7), the original control input and the control input of the controlled system are obtained. The relationship between them can be given asAccording to (9)–(14) in [14] and (8) in this article, we obtainThe deduction of (9) is tedious, so we borrow the similar process in literature [12] for abbreviation.

3.2. Determination of the Matching Conditions

If the controlled equation (5) matches the original (3), the control input determined by (9) is true; that is, is true for any point .

In the same form as gyroscopic forces matrix , matrix is given as follows:where gives out functions and similar to and . Substituting and (11) into (10), we get

Let , multiplying the line vector from the left to (12) and taking the obtained left side zero constantly, we acquire the matching conditions as follows (13)–(16):where and each j represents an equation, and

In equation (14), each pair of (j, f) corresponds to an equation. Otherwise, and j > f, When f + j ≠ 5, h = 0; then h = (−1)f.where and each j corresponds to an equation.

Multiplying (13) by , and multiplying (14) by , , , , , , , , , and , then the sum to obtain a equation has nothing to do with gyroscopic forces:

Remark 2. Equation (17) is obtained without gyroscopic forces terms due to anti-symmetric property of gyroscopic forces matrix, which indicates that the quadratic form of the anti-symmetric matrix is equal to zero.
Let , then (17) can be concisely expressed asFor the controlled kinetic energy matrix, its regular condition is . According to literature [12], is known, so can be obtained from the definition of matrixes , and K. If are zero, then is available. Therefore, cannot be zero at the same time. In order to facilitate the subsequent calculations, suppose . Multiplying (15) by and summing them, we getIn summary, the combination of (16), (18), and (19) and is the condition under which the controlled system matches the original one.

Remark 3. Equations (16) and (18) are partial differential equations, and equations (13)–(15) related to gyroscopic forces terms are linear algebraic equations. They are cascaded. PDEs are resolved at first to decrease difficulty, and then there are only linear algebraic equations which could be solved explicitly with introducing the previous solution.
Except algebraic equations, there are only two PDEs contained in the matching condition for the underactuation degree one system. They could be solved with involving enough independent variables according to some examples applied in CL methods.

3.3. Determination of the Matching Controller

For convenience, some notations can be expressed as follows:

Multiplying from the left to (12), the matching control of the system is

4. Matching Conditions Solution

For the asynchronous motors model, there is no control input at the fifth degree of freedom except the ath degree of freedom, where , at the same time, the following notations are defined as

Introduce the function vector from the definition of matrixes , , and , we get , where , and thus the element in the fifth row of matrix is . To ensure , we choose , where and are constants. The system matching conditions expressed by , , and are

Assume , and find a special solution to (23) as follows:

In the above (23), is a constant.

The matrix of the system is

It can be seen from the matrix that its determinant is , so a sufficient condition for is

From (24), the positive definite solution to the controlled potential energy and the satisfied conditions arewhere are constants. After some calculations, the Hessian matrix of is

It is clear that the Hessian matrix is positive definite. And from , we know that is the minimum point of the controlled potential energy.

Remark 4. The controlled energy consists of controlled kinetic energy and potential energy. On the one hand, the controlled kinetic energy is in quadratic form of velocity and achieves positive definiteness with supplied condition. On the other hand, the controlled potential energy is constructed in the form to be positive definite conveniently.
According to (30), the dissipation matrix is chosen in diagonal form as follows:where can be any value greater than zero.

5. Matching Control Law and Stability Analysis

Based on the above analysis, is known, so the matrix of the system is

In equation (32),  =  and  = , because , where , so we get

From (13)–(15) and from (30)–(32), we have the gyroscopic forces component function as follows:

In equation (33), one part of the value of gyroscopic forces component function needs to satisfy and be arbitrary values, the other part needs to be taken according to the following (35), and the rest can be arbitrary values.

To keep the control law simple, take the value of , to be one and the rest to be zero.

Remark 5. A large part of variables, and in equation (34), take the value zero for convenience to get the simpler control law. The different evaluation for them could still assure stability of the system and maybe affect the convergence rate of the system. Further, there is a room to search control law for better performance in some respects.
Substituting (20), (27)–(35) into (21), the obtained matching control law of the motor system is described asIn summary, the conclusion is listed as follows:.

Proposition 1. For AC asynchronous motor systems, if the parameters of controller meet the following conditions:then the smooth feedback control law expressed by equations (35)–(38) can stabilize the motor globally asymptotically at . It is the desired equilibrium point at which the controlled potential energy achieves the minimum.

Proof. Let Lyapunov candidate function . If the system controller parameters are selected according to (39), then the function is positive definite. From equation (6), there is . Therefore, the control law given by equations (35)–(38) enables the induction motor to achieve global asymptotic stabilization at .
For the asymptotic stability, it can be proved that there is no trajectory of isolated points other than equilibrium points in the set of .
Assuming that there is such a trajectory in the set, it can be obtained asThere is a certain point on this trajectory, and at the same time, (39) also holds in a certain area of this point . Differentiating and integrating (40) along the trajectory, we getwhere are constant.
Substituting (35)–(38), (40), and (41) into the first four equations of (1), we getIt can be seen from observation that (43) is true if and only if , and vice versa. Therefore, the assumption that there are isolated points belonging to the point group of does not hold.
From (4) and (41), we obtain . This shows that the trajectory in the set can only be the equilibrium point, and the Hessian matrix of is positive definite, so point is the sole extreme point of . According to the LaSalle principle, the system can achieve global asymptotic stabilization by the proposed control law.

6. Simulation Result Analysis

The simulation parameters of the system are selected as L11 = 0.45 H, L13 = 0.42 H, L33 = 0.45 H, J = 0.2 kg m2H, rad/s, T1 = 8Nm, H = 8Nm, np = 8, and R1 = R2 = 0.97 Ω. The parameters of the controller that satisfy the positive definite of the controlled energy are k5 = 0.1, k6 = 4, k7 = 6, k8 = 0.5, and k9 = 0.4.

When d1 = 7, d2 = 7, d3 = 5, and d4 = 5, the simulation results are as follows:

The desired equilibrium point of the system is , where . It can be seen from Figures 13 that when the AC asynchronous motor is disturbed by some uncertain factors, the control target of the system deviates from the expected balance point. Under the control input of Figure 4, the control target of the system can return to the desired equilibrium point as soon as possible. In this process, the change of electromagnetic torque, original, and controlled energy are showed in Figures 5 and 6.

Remark 6. Parameters k5, k6, k7, k8, and k9 in controller are connected with the controlled kinetic energy, and their evaluation should ensure controlled kinetic energy positive definite. Furthermore, they will affect primarily the amplitude of inputs.

Remark 7. Parameters d1, d2, d3, and d4 are related only to the dissipated forces, and thus the changing of them could be able to improve the convergence rate.
Compared with other control methods [8, 9], the nonlinear smooth feedback control law obtained in this article has a larger convergence range for its global asymptotic stabilization.

7. Conclusion

The control of an asynchronous motor is studied in this paper. By applying the controlled Lagrange function method to high-order, strongly coupled, and time-varying nonlinear systems, the controlled equation matching the original equation of the system is derived by using the expected controlled energy and its time derivative. Because the primary term of velocity exists in the original equation, the gyroscopic forces of the generalized primary term of velocity is innovatively introduced into the controlled equation, and the condition of complete matching between the original equation and the controlled equation can be obtained. By solving the matching condition composed of some partial differential equations, the specific matching control law of the system is obtained, and the global asymptotic stabilization of not only velocity but also position can be realized at the desired equilibrium point at the same time. Finally, the controlled energy of the system is chosen as the Lyapunov function for its property, which facilitates the proof of stability.

In the research process, we can see that the CL method analyzes the asynchronous motor system from the perspective of underactuated mechanical system, so that the controlled system maintains the Lagrangian mechanical structure in form, and the nonlinear smooth feedback control law can be obtained, which has a large convergence range and is helpful to realize robust control and optimal control. On the basis of the work in this paper, the nonlinear CL method will be improved by the introducing observer and be intended to solve tracking problem in following work.

Data Availability

The data supporting the findings of this study are available within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding this research.

Acknowledgments

The authors gratefully acknowledge the financial support provided by the National Natural Science Foundation of China (no. 61164010), and the Natural Science Foundation of Gansu Province (no. 20JR10RA210), and Natural Science Foundation of Gansu Province (no. 20JR10RA273).