#### Abstract

The structural scheme of mechanical elastic energy storage (MEES) system served by permanent magnet synchronous motor (PMSM) and bidirectional converters is designed. The aim of the research is to model and control the complex electromechanical system. The mechanical device of the complex system is considered as a node in generalized coordinate system, the terse nonlinear dynamic model of electromechanical coupling for the electromechanical system is constructed through Lagrange-Maxwell energy method, and the detailed deduction of the mathematical model is presented in the paper. The theory of direct feedback linearization (DFL) is applied to decouple the nonlinear dynamic model and convert the developed model from nonlinear to linear. The optimal control theory is utilized to accomplish speed tracking control for the linearized system. The simulation results in three different cases show that the proposed nonlinear dynamic model of MEES system is correct; the designed algorithm has a better control performance in contrast with the conventional PI control.

#### 1. Introduction

Power balance is a basis for stable operation of power system [1, 2], which is frequently disturbed by some factors such as faults and intermittent power generations. These threatening behaviors may lead to oscillation between the generating units and the interconnected systems, and even worse, which probably cause serious large-scale blackout in power system [3]. Currently, one of the principal methods solving the power imbalance is to reserve adequate spare capacity for the system, which not only increases the investment, but also reduces the utilization of electrical equipments. With the development of energy storage (ES) technologies, a new way of utilizing energy storage technology to solve the power imbalance has been paid more and more attentions [4, 5]. From state-of-the-art technologies of ES, a great variety of the technologies have been developed [6–11]. However, from a viewpoint of power system operation, the technologies of ES at least need to have the following features: large capacity, fast power response speed, high efficiency, long service life, low operation cost and environmental friendship, and so forth. In terms of these requests, most of ES technologies are still immature except pumped-storage. A new method utilizing spiral springs as ES medium is proposed in reference [12], which is called mechanical elastic energy storage (MEES); some of the basic characteristics for the technology are discussed in [13]. Due to the benefits of small size, low loss, and high efficiency [14], permanent magnet synchronous motor (PMSM) is chosen as the servomotor for the MEES system. Thus, the MEES system, which is integrated by the mechanical device (spiral springs and the correspondingly coating shell), the electrical installation (PMSM), and relevant control information elements, is a typically complex nonlinear system of electromechanical coupling. Parameter coupling, which is induced between electromagnetic parameters of PMSM and mechanical parameters of spiral springs, may affect the dynamic performance of the system. The preexisting studies have shown that revealing the law of electromechanical interaction from the point of view of energy is an effective solution to research dynamic characteristics of complex electromechanical system. Lagrange-Maxwell equation of flywheel energy storage is deduced in reference [15] by calculation of integral of magnetic density in the area of magnetic field.

The modeling and control of complex dynamic systems are the rapidly emerging researching topics [16], especially for the nonlinearities in the systems [17]. The paper concentrates on modeling and decoupling control of complex PMSM based MEES system in detail. The mechanical device is taken into account in the analysis of electromechanical coupling, which can make the research become closer to the reality. Especially for the mechanical spiral springs with a big quality and large inertia, neglecting their influence will make the results inaccurate, even unavailable. The overall aim of the paper is to model and control the complex electromechanical system based on the characteristics of mechanical spiral springs in order to improve the operating stability of the system; in particular, the main innovative points of the research are represented as follows:(a)previous researches mainly focus on the motor ontology for the electromechanical issue; the difference in the research is to consider the mechanical device as a node in generalized coordinate system; the terse nonlinear dynamic model for a complex electromechanical system is deduced and constructed through Lagrange-Maxwell energy method by the principle of electromechanical coupling. (b)the generalized coordinates of MEES system are determined. The functions of the kinetic energy, magnetic energy, the potential energy, and the energy dissipation of the whole system are obtained, and the differential equations with electromechanical parameters of MEES system are established by applying the extended Lagrange-Maxwell equation.(c)the theory of direct feedback linearization (DFL) is applied to decouple the nonlinear model and convert the developed model for MEES system from nonlinear to linear. An optimal linear controller is designed to accomplish the speed tracking control for the linearized model.

The paper is organized as the following: it starts with introductions to the ES technologies and the electromechanical characteristics of MEES system in Section 1. The constitution and operating principles of the MEES system are described in Section 2. Mathematical modeling of electromechanical coupling for the complex MEES system is deduced in Section 3. In Section 4, DFL based decoupling control of the dynamic model is proposed, and an optimal linear controller is designed to accomplish the speed tracking control for the linearized model. Section 5 gives the results of numerical simulation. Finally, the conclusions are summarized in Section 6.

#### 2. Preliminaries and Problem Description

##### 2.1. Constitution of MEES System

In order to give the readers a better understanding of the operating principle of the research object, the composition scheme of MEES system served by PMSM and bidirectional converters is described in Figure 1.

PMSM based MEES system consists of several elements: ES box (multiple spiral springs wrapped up in it), gear box, PMSM, bidirectional converters, breaker, and controller. Therefore, the electromechanical coupling is caused by the interaction between electromagnetic parameters of PMSM and mechanical parameters of spiral springs. A MEES power station can be built by getting tens or even more units together. Through a central controller, the units in the power station can be controlled orderly to realize the goal of large-scale energy storage and power generation. It is important to note that a linkage form structure for ES box is proposed in order to increase the energy storage capacity; for details see [13].

##### 2.2. Operating Principle of MEES System

The basic operation principle of the MEES system is concerned with two processes: energy storage and power generation. The load or drive source of MEES system in energy storage or power generation is spiral springs in ES box. In the process of energy storage, PMSM, which is driven by power grid, tightens the springs through gear box. Then, the electric energy is stored in the form of elastic deformation energy. Once the unit receives a signal of releasing energy, the tightened springs start to release energy and drive the motor to generate electricity.

#### 3. Modeling of Electromechanical Coupling for Complex MEES System

##### 3.1. Physical Analysis

The operating rotation speed of the spindle in elastic ES box in the process of energy storage or power generation is required to run basically stable owing to the large size and the big inertia of ES box. As mentioned above, the elastic ES box can be considered as an output load of PMSM in ES. Hence, supposing that the rotation speed of the shaft in ES box is a constant, the torque calculation formula of ES box in ES process can be given by where and denote the operating torque and initial torque, , , , and are the modulus of elasticity, width, thickness, and length of spring material, and denotes the time. It should be noted that must be ensured between the maximum torque and minimum torque. Additionally, it should be compatible with the electromagnetic torque of PMSM. In terms of formula (1), the characteristic curve for ES box in ES can be shown in Figure 2.

The servomotor of PMSM consists of two parts, stator and rotor. When three-phase alternating current (AC) with mutual phase difference of 120° is flowed into the stator windings of the motor, a rotating magnetic field with a uniform motion in space is induced. The rotation speed of the magnetic field is related to the frequency of sinusoidal wave in stator windings. The driven torque, which is produced by the interaction between the rotating magnetic field of stator and permanent magnetic field of rotor, makes the rotor rotate and achieves the conversion of electric energy to elastic energy.

In the process of ES, the torque of the springs raises with the operation of the system. The power angle between the axis of stator and rotor becomes larger, which augments the magnetic torque to balance the increasing torque of springs. With the rising of energy stored in ES box, the input power and the stator current increase accordingly. Simultaneously, the armature reaction causes the growth of flux linkage of air gap and back electromotive force (e.m.f) of stator and compels the stator voltage to rise. In other words, the interaction between the electromagnetic parameters of electrical system and the mechanical parameters of mechanical system composes the electromechanical coupling in complex MEES system.

##### 3.2. Mathematical Modeling

Before modeling the system of electromechanical coupling, six assumptions are described as follows:(a)the saturation of the core is neglected,(b)the losses of eddy and hysteresis are ignored,(c)the distribution of air gap is uniform,(d)the self-inductance and mutual inductance among the windings are independent of the position of the rotor,(e)the damping winding in the rotor and the damping effect of permanent magnet are both neglected,(f)back e.m.f. maintains the sinusoidal variation.

Based on the principle of electromechanical coupling dynamics, Lagrange-Maxwell energy method is applied to construct the dynamic mathematical model of MEES system as follows.

*(**1) Determining Generalized Coordinates of the System.* Owing to the large size, big quality, and high inertia of the mechanical spiral springs, neglecting their influence will make the results inaccurate, even unavailable. Hence, the mechanical device is taken into account and considered as a node in generalized coordinate system. Six generalized coordinates shown in Table 1 are constructed for the MEES system, four of which belong to electromagnetic system (), and the remaining two are mechanical system (). The parametric descriptions of six generalized coordinates shown in Table 1 are as follows: , , and are three-phase AC flowed in stator, , , and indicate three-phase voltage applied in stator, represents the magnetic flux of permanent magnet in rotor, and are the angular displacement of rotor and the spindle of ES box, respectively, and are the mechanical angular velocity of rotor and spindle of ES box, respectively, and and represent the torque of gear and load. is equal to the torque of the springs in the ES box in energy storage.

*(**2) Finding Out the Energy Relationship for the MEES System. *The kinetic energy of the MEES system can be given as
where is the moment of inertia of PMSM, and are the equivalent moment of inertia of the gear and ES box with respect to input shaft, and is the transformation ratio of gear.

For the MEES system, the magnetic energy , which is jointly produced by the stator current, permanent magnet of rotor, and interaction with each other, consists of three parts. The calculation of the magnetic energy produced by the current of stator is given as where and indicate the self-inductance and mutual inductance of three-phase windings in stator, for an exterior-gibbose PMSM, , .

For the distributed winding with a phase belt of 60°, due to the number of the slots for per pole per phase is often more than two, the magnetic energy of air gap can be considered not changed with the rotation of rotor, namely, the magnetic energy produced by permanent magnet of rotor is a constant, ( is a constant).

is the magnetic energy produced by the interaction between the stator current and stator flux generated by the rotor; the calculation formula of can be written as

The potential energy of the system is given by where is a torsion constant of the springs.

The Lagrange function can be written as

The dissipation function of magnetic system can be expressed as where , , and indicate three-phase resistance of the stator.

The dissipation function of mechanical system can be written as where , , and are the equivalent viscous damping coefficients of the motor, the gear, and the spring, respectively.

Then, the dissipation function of the whole system can be given as

Nonconservative generalized force of mechanical system () can be expressed as

External e.m.f. of electromagnetic system is given by

The Lagrange-Maxwell equation of MEES system can be expressed as

For electromagnetic system, firstly, the stator winding of A-phase () is considered as the voltage equation of stator winding of A-phase can be written as

Similarly, the voltage equation of stator winding of B-phase () can be written as the voltage equation of stator winding of C-phase () can be written as

For permanent magnet of rotor, the effective flux of permanent magnet is generally considered to be a constant, namely, = ( is a constant). Thus

For mechanical system, when has the motion equation of mechanical system can be expressed as Let and ; the motion equation of mechanical system can be rewritten as

For mechanical system, when has the motion equation of mechanical system can be written as In terms of , , that is, (22) can be rewritten as

The electromechanical dynamic model of MEES system in static ABC coordinates can be obtained through simultaneous equations of (14), (15), (16), (20), and (23).

Define as the transformation matrix from static coordinates of ABC to rotating frame of dq0 as follows: in terms of the transformation matrix , the relationship between current and voltage in ABC frame and dq0 frame can be written as Substituting (25) into equations (14), (15), and (16), respectively, the voltage equations of the electromagnetic system in dq0 frame can be obtained as where is the electric angular speed, is the number of pole-pairs, is the resistance of the stator, , are the equivalent inductance of direct axis and quadrature axis, and .

Similarly, the motion equations of mechanical system in rotating frame of dq0 are given as Due to and , two formulas in (27) express the same meaning essentially. Considering , namely, , (27) can be transformed into a new form: The electromechanical dynamic model of MEES system in rotating frame of dq0 can be obtained through simultaneous equations of (26) and (28).

#### 4. Decoupling Control of the Dynamic Model

##### 4.1. Fundamental Principle of Direct Feedback Linearization

For an affine nonlinear system in the standard form, where are state variables, are input variables, are output variables, , are smooth vectors fields defined on , and are smooth scalar functions defined on . The nonlinearity of MEES dynamic model in rotating frame of dq0 is caused by product terms of and , which will weaken the performance of conventional PI controller. Linearization and decoupling control of the system are effective solutions to improve the servo capability. The theory of direct feedback linearization (DFL) is known intimately up to now and is applied widespread to the varieties of nonlinear domains [18, 19]. The basic principle of DFL is to differentiate the output function repeatedly, until the th derivatives of output variable explicitly contain control variable , where is the relatively degree of the system.

##### 4.2. DFL for MEES System

The state variables , output variables , and control variables are chosen as , , ,* *, , and , , respectively. The output variables are differentiated with respect to time until at least an input variable explicitly appears by DFL theory.

Taking the derivative of , we can obtain (30) contains the input and the relatively degree . Introducing a new variable makes .

Taking the derivative of , we can obtain as (31) does not contain any input variable, second derivation calculus to can be obtained as follows: (32) contains the input and the relatively degree . Introducing another new variable makes .

The comprehensive order of the system equals , which is exactly equal to the degree of the original system. Hence, a nonlinear state feedback control law is devised as follows: where , , , Substituting (33) into the electromechanical dynamic model in rotating frame of dq0, the ultimate linearization equation of the MEES system can be obtained as follows: where , , and are three matrices.

##### 4.3. Design of Optimal Linear Control

Optimal linear control block diagram of speed tracking for linearized system (35) is designed in Figure 3.

For the linearized system (35), the theory of linear control can be applied to determine . In order to improve the control performance, in terms of optimal control theory of linear system, optimal control inputs are chosen as where , is the matrix of weight coefficient.

is the solution of Riccati equation as follows: where is a positive definite matrix.

#### 5. Simulation and Results

Based on a 0.16 kWh/0.8 kW MEES system being developed, the operation parameters of the system are given in Table 2. In terms of the description above, the control objectives in the paper are determined as the speed of rotor and the current of -axis, namely, the speed of rotor tracking the reference speed and the current of -axis being zero. The reference speed of the motor is set at 600 r/min. The work of modeling and control is implemented in Matlab environment.

In terms of the operation parameters in Table 2, the matrix parameters and in (35) are obtained as follows:

In the paper, the matrices of weight coefficient and positive definite are taken as 1 and a identity matrix, respectively. Then, the matrix of can be calculated as Considering (36), the inputs of the optimal control can be expressed as Substituting (41) into (33), the nonlinear state feedback controller is given as

In order to evaluate the effectiveness of the proposed algorithm properly and adequately, the simulation schemes are selected as three different cases, which are case I (always cutting off the energy storage device), case II (from the initial state of storing no energy to the final state of storing full of energy), and case III (from the initial state of storing half of full of energy to the final state of storing full of energy). The simulation results are shown in Figures 4, 5, and 6, respectively. In order to better display the control performances, the proposed algorithm is in contrast with the conventional PI control. The proportional and integral coefficients for -axis PI current controller are set as 2.6 and 50, respectively, and the proportional and integral coefficients for PI speed tracking controller are set as 0.24 and 10, respectively. Case I simulates PMSM operating in no load state, which is to test the no-load performance of the servomotor. The results in Figure 4 show that the proposed algorithm makes the controlled objective become stable more quickly and almost has no oscillation in stabilization process in comparison with the conventional PI control. Figure 5 shows the results of case II, which essentially represents the process from point to point in Figure 2. Figure 6 displays the results of case III, and it should be noted that the running time in the condition will be half of the whole time of energy storage, which is 5 s. The simulation results in case II and case III are similar to case I. Therefore, the simulation results reveal that in contrast with conventional PI control, the presented algorithm not only has a preferable ability of speed tracking in all of the three cases, but also makes the electromechanical system have a better dynamic performance in energy storage process.

**(a) Electric angular speed of rotor**

(b) Current of -axis |

**(a) Electric angular speed of rotor**

(b) Current of -axis |

**(a) Electric angular speed of rotor**

(b) Current of -axis |

#### 6. Conclusions

In terms of the characteristics of electromechanical coupling for PMSM based MEES system presented in the paper, dynamic modeling and control of the system are proposed from an electromechanical point of view. The conclusions were made as follows:(a)dynamic model for PMSM based MEES system is nonlinear with significant electromechanical coupling. The nonlinear model constructed by Lagrange-Maxwell energy method can reflect the characteristics of electromechanical coupling for the system;(b)DFL can perform the dynamic decoupling for the complex electromechanical system completely;(c)the speed tracking controller proposed in the paper in terms of optimal linear control theory for linearized system has a better performance in contrast with conventional PI control algorithm, and the simulation results in three different cases prove the validity and feasibility of the designed method.

#### Conflict of Interests

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

#### Acknowledgments

This research is supported by the Specialized Research Fund for the Doctoral Program of Higher Education of China under Grant no. 20120036130001, the Fundamental Research Funds for the Central Universities of China under Grant no. 11MG40, and the Independent Research Funds of State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources of China under Grant no. 201209.