JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 603063 10.1155/2013/603063 603063 Research Article Dynamic Modeling and Control of Electromechanical Coupling for Mechanical Elastic Energy Storage System http://orcid.org/0000-0002-3856-0535 Yu Yang Mi Zengqiang Chen Michael State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources North China Electric Power University Baoding 071003 China ncepu.edu.cn 2013 22 10 2013 2013 12 07 2013 04 09 2013 2013 Copyright © 2013 Yang Yu and Zengqiang Mi. 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.

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 . 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 . 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 , which is called mechanical elastic energy storage (MEES); some of the basic characteristics for the technology are discussed in . Due to the benefits of small size, low loss, and high efficiency , 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  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 , especially for the nonlinearities in the systems . 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:

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.

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.

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.

The composition scheme of MEES system served by PMSM and bidirectional converters.

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 .

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 (1)T=EBH3π6Lωt+T0, where T and T0 denote the operating torque and initial torque, E, B, H, and L are the modulus of elasticity, width, thickness, and length of spring material, and t denotes the time. It should be noted that T 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.

Characteristics curve for ES box in ES.

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:

the saturation of the core is neglected,

the losses of eddy and hysteresis are ignored,

the distribution of air gap is uniform,

the self-inductance and mutual inductance among the windings are independent of the position of the rotor,

the damping winding in the rotor and the damping effect of permanent magnet are both neglected,

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 (k=1,2,3,4), and the remaining two are mechanical system (j=1,2). The parametric descriptions of six generalized coordinates shown in Table 1 are as follows: ia, ib, and ic  are three-phase AC flowed in stator, ua, ub, and uc  indicate three-phase voltage applied in stator, ψf represents the magnetic flux of permanent magnet in rotor, ϑ1 and ϑ2 are the angular displacement of rotor and the spindle of ES box, respectively, Ω1 and Ω2 are the mechanical angular velocity of rotor and spindle of ES box, respectively, and Tg and TL represent the torque of gear and load. TL is equal to the torque T of the springs in the ES box in energy storage.

Six generalized coordinates for MEES system.

Generalized coordinates Electromagnetic system Mechanical system
Stator Rotor Gear ES box
k = 1 k = 2 k = 3 k = 4 j = 1 j = 2
e k ( q j ) ψ f θ 1 θ 2
e ˙ k ( q ˙ j ) i a i b i c Ω 1 Ω 2
U k ( Q j ) u a u b u c - T g - T L

( 2) Finding Out the Energy Relationship for the MEES System. The kinetic energy of the MEES system Wk can be given as (2)Wk=12J1Ω12+12JeΩ12+12JspΩ22=12J1Ω12+12JeΩ12+12Jsp(Ω1i)2, where J1 is the moment of inertia of PMSM, Je and Jsp are the equivalent moment of inertia of the gear and ES box with respect to input shaft, and i is the transformation ratio of gear.

For the MEES system, the magnetic energy Wm, 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 Wm1 produced by the current of stator is given as (3)Wm1=12Laia2+12Lbib2+12Lcic2+Miaib+Mibic+Micia, where L and M indicate the self-inductance and mutual inductance of three-phase windings in stator, for an exterior-gibbose PMSM, Maa=Mbb=Mcc=Mab=Mbc=Mca=M, La=Lb=Lc.

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

W m 12 is the magnetic energy produced by the interaction between the stator current and stator flux generated by the rotor; the calculation formula of Wm12 can be written as (4)Wm12=iaψfcosθ1+ibψfcos(θ1-23π)+icψfcos(θ1+23π).

The potential energy of the system V is given by (5)V=12αspθ22, where αsp is a torsion constant of the springs.

The Lagrange function can be written as (6)L=Wk+Wm-V=12J1Ω12+12JeΩ12+12JspΩ22+12Laia2+12Lbib2+12Lcic2+Miaib+Mibic+Micia+C1+iaψfcosθ1+ibψfcos(θ1-23π)+ibψfcos(θ1+23π)-12αspθ22.

The dissipation function of magnetic system Fe can be expressed as (7)Fe=12Raia2+12Rbib2+12Rcic2, where Ra, Rb, and Rc indicate three-phase resistance of the stator.

The dissipation function of mechanical system Fm can be written as (8)Fm=12(B1+Be)Ω12+12Bsp(Ω1i)2, where B1, Be, and Bsp are the equivalent viscous damping coefficients of the motor, the gear, and the spring, respectively.

Then, the dissipation function of the whole system F can be given as (9)F=Fe+Fm=12Raia2+12Rbib2+12Rcic2+12(B1+Be)Ω12+12Bsp(Ω1i)2.

Nonconservative generalized force of mechanical system (j=1,2) can be expressed as (10)Q1=-Tg,j=1,Q2=-TL,j=2.

External e.m.f. of electromagnetic system is given by (11)Uk=-uk.

The Lagrange-Maxwell equation of MEES system can be expressed as (12)ddt(Le˙k)-Lek+Fe˙k=Uk,ddt(Lq˙j)-Lqj+Fq˙j=Qj.

For electromagnetic system, firstly, the stator winding of A-phase (k=1) is considered as (13)Le1=0,Le˙1=Lia=Lia+Mib+Mic+ψfcosθ1,ddt(Le˙1)ddt(Lia)=Ldiadt+Mdibdt+Mdicdt+d(ψfcosθ1)dt,Fe˙1Fia=Raia,U1=ddt(Le˙1)-Le1+Fe˙1,ua=U1; the voltage equation of stator winding of A-phase can be written as (14)ua=U1=ddt(Le˙1)-Le1+Fe˙1=Ldiadt+Mdibdt+Mdicdt+d(ψfcosθ1)dt+Raia.

Similarly, the voltage equation of stator winding of B-phase (k=2) can be written as (15)ub=U2=ddt(Le˙2)-Le2+Fe˙2=Ldibdt+Mdicdt+Mdiadt+d(ψfcos(θ1-(2/3)π))dt+Rbib; the voltage equation of stator winding of C-phase (k=3) can be written as (16)uc=U3=ddt(Le˙3)-Le3+Fe˙3=Ldicdt+Mdiadt+Mdibdt+d(ψfcos(θ1+(2/3)π))dt+Rcic.

For permanent magnet of rotor, the effective flux ψf of permanent magnet is generally considered to be a constant, namely, ψf = C2(C2 is a constant). Thus (17)Lψf=0,Lψ˙f=0.

For mechanical system, when j=1 has (18)Lq1Lθ1=-iaψfsinθ1-ibψfsin(θ1-23π)-icψfsin(θ1+23π)-αspθ1i2,Lq˙1LΩ1=(J1+Je)Ω1+JspΩ1i2,ddt(Lq˙1)ddt(LΩ1)=(J1+Je+Jspi2)dΩ1dt,Fq˙1FΩ1=(B1+Be+Bspi2)Ω1,Q1=ddt(Lq˙1)-Lq1+Fq˙1,Q1=-Tg, the motion equation of mechanical system j=1 can be expressed as (19)(J1+Je+Jspi2)dΩ1dt=-iaψfsinθ1-ibψfsin(θ1-23π)-icψfsin(θ1+23π)-αspθ1i2-(B1+Be+Bspi2)Ω1-Tg. Let J=J1+Je+(Jsp/i2) and B=B1+Be+(Bsp/i2); the motion equation of mechanical system j=1 can be rewritten as (20)JdΩ1dt=-iaψfsinθ1-ibψfsin(θ1-23π)-icψfsin(θ1+23π)-αspθ1i2-BΩ1-Tg.

For mechanical system, when j=2 has (21)Lq2Lθ2=-iiaψfsinθ2-iibψfsin(θ2-23π)-iicψfsin(θ1+23π)-αspθ2,Lq˙2LΩ2=i2(J1+Je)Ω2+JspΩ2,ddt(Lq˙2)ddt(LΩ2)=[i2(J1+Je)+Jsp]dΩ2dt,Fq˙2FΩ2=[i2(B1+Be)+Bsp]Ω2,Q2=ddt(Lq˙2)-Lq2+Fq˙2,Q2=-TL, the motion equation of mechanical system j=2 can be written as (22)i2(J1+Je+Jspi2)dΩ2dt=-iiaψfsinθ2-iibψfsin(θ2-23π)-iicψfsin(θ2+23π)-αspθ2-i2(B1+Be+Bspi2)Ω2-TL. In terms of J=J1+Je+(Jsp/i2), B=B1+Be+(Bsp/i2), that is, (22) can be rewritten as (23)i2JdΩ2dt=-iiaψfsinθ2-iibψfsin(θ2-23π)-iicψfsin(θ2+23π)-αspθ2-i2BΩ2-TL.

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 Cabc/dq0 as the transformation matrix from static coordinates of ABC to rotating frame of dq0 as follows: (24)Cabc/dq0=23[cosθcos(θ-23π)cos(θ+23π)-sinθ-sin(θ-23π)-sin(θ+23π)121212]; in terms of the transformation matrix Cabc/dq0, the relationship between current and voltage in ABC frame and dq0 frame can be written as (25)[idiqi0]=Cabc/dq0[iaibic],[uduqu0]=Cabc/dq0[uaubuc]. Substituting (25) into equations (14), (15), and (16), respectively, the voltage equations of the electromagnetic system in dq0 frame can be obtained as (26)[uduq]=[Ld00Lq][diddtdiqdt]+[Rsid-LqωeiqRsiq+Ldωeid+ωeψf], where ωe=pnΩ1=pnθ˙1 is the electric angular speed, pn is the number of pole-pairs, Rs is the resistance of the stator, Rs=Ra=Rb=Rc, Ld,Lq are the equivalent inductance of direct axis and quadrature axis, and Ld=Lq=Ls=3(L-M)/2.

Similarly, the motion equations of mechanical system in rotating frame of dq0 are given as (27)[J00i2J][θ¨1θ¨2]+[B00i2B][θ˙1θ˙2]+[αspi200αsp][θ1θ2]+[-32[iqψf+(Ld-Lq)idiq]-32i[iqψf+(Ld-Lq)idiq]]=[-Tg-TL]. Due to θ2=θ1/i and Tg=TL/i, two formulas in (27) express the same meaning essentially. Considering ωe=pnΩ1=pnθ˙1, namely, θ¨1=ω˙e/pn, (27) can be transformed into a new form: (28)Jpnω˙e+Bpnωe-32[iqψf+(Ld-Lq)idiq]+αspi2θ1=-Tg, 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, (29)x=f(x)+i=1pgi(x)ui,yi=hi(x),222222222222222222i=1,2,,p, where xRn are state variables, uRp are input variables, yRp are output variables, f,g are smooth vectors fields defined on Rn, and h are smooth scalar functions defined on Rn. The nonlinearity of MEES dynamic model in rotating frame of dq0 is caused by product terms of ωeid and ωeiq, 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 rth derivatives y(r) of output variable explicitly contain control variable u, where r is the relatively degree of the system.

4.2. DFL for MEES System

The state variables xR3, output variables yR2, and control variables uR2 are chosen as id, iq, ωe,   id, ωe, and ud, uq, 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 id, we can obtain (30)i˙d=-RsLdid+ωeiq+1Ldud; (30) contains the input ud and the relatively degree r1=1. Introducing a new variable v1 makes i˙d=v1.

Taking the derivative of ωe, we can obtain (31)ω˙e=-BJωe+3pn2Jiqψf-pnαspJi2θ1+pnJTg; as (31) does not contain any input variable, second derivation calculus to ωe can be obtained as follows: (32)ω¨e=BJω˙e+3pn2Ji˙qψf-pnαspJi2θ˙1+pnJT˙g=BJ(-BJωe+3pn2Jiqψf-pnαspJi2θ1+pnJTg)+3pnψf2J[-1Lq(Rsiq+Ldωeid+ωeψf)+1Lquq]-pnαspJi2ωepn+pnJT˙g; (32) contains the input uq and the relatively degree r2=2. Introducing another new variable v2 makes ω˙e=v2.

The comprehensive order of the system equals r1+r2=3, which is exactly equal to the degree of the original system. Hence, a nonlinear state feedback control law is devised as follows: (33)u=E(x)+D(x)v, where u=[ud,uq]T, v=[v1,v2]T, D(x)=[Ld002JLq/3pnψf], (34)E(x)=[Rsid-Ldωeiq[2LqB3pnψf(-BJωe+3pn2Jiqψf-pnαspJi2θ1+pnJTg)+(Rsiq+Ldωeid+ωeψf)+αspJi2ωe-pnJT˙g]]. 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: (35)[i˙di˙qω˙e]=[0000BpnJ-(2B23pnψfJ+2αsp3ψfi2)03pnψf2J-BJ][idiqωe]+[1002J3pnψf00][v1v2]+[02B3ψfJTg+23ψfT˙g-2Bpnαsp3ψfJi2θ1pnJTg-pnαspJi2θ1]=Ax+Ev+C, where A, E, and C 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.

Control block diagram of speed tracking for linearized system.

For the linearized system (35), the theory of linear control can be applied to determine v. In order to improve the control performance, in terms of optimal control theory of linear system, optimal control inputs are chosen as (36)v=-K*·x, where K*=R-1ETP*, R=RT is the matrix of weight coefficient.

P * is the solution of Riccati equation as follows: (37)ATP+PA-PER-1ETP+Q=0, where Q 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 d-axis, namely, the speed of rotor tracking the reference speed and the current of d-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.

Operation parameters of a 0.16 kWh/0.8 kW MEES system.

 Rated work voltage Un 380 V Rated work torque 2.1 N·m Rated rotating speed of ES box 15 r/min Number of pole pairs for rotor pn 4 Rated rotating speed of rotor 600 r/min Transformation ratio of the gear i 40 : 1 Moment of inertia for rotor J1 0.8 × 10−3 kg m2 Equivalent moment of inertia for the gear Je 0.0013 kg m2 Equivalent moment of inertia for the ES box Jsp 4.3264 kg m2 Torsion constant of spring αsp 5.0178 Effective flux of permanent magnet ψf 0.1442 Wb Resistance of stator Rs 1.95 Ω Equivalent inductance of direct axis Ld 0.2541 H Equivalent inductance of quadrature axis Lq 0.2541 H Viscous damping coefficient of the motor B1 0.015 N/rad/s Viscous damping coefficient of the gear Be 0.013 N/rad/s Viscous damping coefficient of the spring Bsp 0.005 N/rad/s

In terms of the operation parameters in Table 2, the matrix parameters A and E in (35) are obtained as follows: (38)A=[0000BpnJ(2B23pnψfJ+2αsp3ψfi2)03pnψf2J-BJ]=[00000.0019-0.232100.19990.0076](39)E=[1002J3pnψf00]=[1005.002900].

In the paper, the matrices of weight coefficient R and positive definite Q are taken as 1 and a 3×3 identity matrix, respectively. Then, the matrix of P* can be calculated as (40)P*=[1000-0.000228-0.0173180-4.985561-0.000265]. Considering (36), the inputs of the optimal control can be expressed as (41)v=-K*·x=-R-1ETP*·x=[-id0.001141iq+0.086640ωe]. Substituting (41) into (33), the nonlinear state feedback controller is given as (42)u=[1.6969id-0.2541iqωe1.953387iq+0.265858ωe+0.2541idωe]+[00.008956Tg-0.924108T˙g-0.000449θ1].

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 d-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 J to point E 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.

The simulation results of case I (always cutting off the energy storage device).

Electric angular speed of rotor

Current of d-axis

The simulation results of case II (from the initial state of storing no energy to the final state of storing full of energy).

Electric angular speed of rotor

Current of d-axis

The simulation results of case III (from the initial state of storing half of full of energy to the final state of storing full of energy).

Electric angular speed of rotor

Current of d-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:

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;

DFL can perform the dynamic decoupling for the complex electromechanical system completely;

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.

Tan K. T. So P. L. Chu Y. C. Chen M. Z. Q. Coordinated control and energy management of distributed generation inverters in a microgrid IEEE Transactions on Power Delivery 2013 28 2 704 713 Tan K. T. Peng X. Y. So P. L. Chu Y. C. Chen M. Z. Q. Centralized control for parallel operation of distributed generation inverters in microgrids IEEE Transaction on Smart Grid 2012 3 4 1977 1987 Bernstein A. Bienstock D. Hay D. Meric U. Gil Z. Sensitivity analysis of the power grid vulnerability to large-scale cascading failures ACM Sigmtrics Performance Evaluation Review 2012 40 3 33 37 Black M. Strbac G. Value of bulk energy storage for managing wind power fluctuations IEEE Transactions on Energy Conversion 2007 22 1 197 205 2-s2.0-33947685476 10.1109/TEC.2006.889619 Divya K. C. Østergaard J. Battery energy storage technology for power systems—an overview Electric Power Systems Research 2009 79 4 511 520 2-s2.0-59349091439 10.1016/j.epsr.2008.09.017 Cohen A. I. Wan S. H. An algorithm for scheduling a large pumped storage plant IEEE Transactions on Power Apparatus System 1985 104 8 2099 2104 2-s2.0-0022103614 Pimm A. J. Garvey S. D. Drew R. J. Shape and cost analysis of pressurized fabric structures for subsea compressed air energy storage Proceedings of the Institution of Mechanical Engineers C 2011 225 5 1027 1043 2-s2.0-79956020981 10.1177/0954406211399506 Hebner R. Beno J. Walls A. Flywheel batteries come around again IEEE Spectrum Magazine 2002 39 4 46 51 2-s2.0-0036556938 10.1109/6.993788 Ali M. H. Wu B. Dougal R. A. An overview of SMES applications in power and energy systems IEEE Transactions on Sustainable Energy 2010 1 1 38 47 2-s2.0-77955699143 10.1109/TSTE.2010.2044901 Rufer A. Hotellier D. Barrade P. A supercapacitor-based energy storage substation for voltage compensation in weak transportation networks IEEE Transactions on Power Delivery 2004 19 2 629 636 2-s2.0-1942472070 10.1109/TPWRD.2004.824408 Senthil K. Mitra S. Roy A. Sharma A. Chakravarthy D. P. Compact inductive energy storage pulse power system Review of Scientific Instruments 2012 83 5 4 054703 Yu Y. Mi Z. Q. Nonlinear dynamic model and chaotic characteristics of mechanical elastic energy storage unit in energy storage process Acta Physica Sinica 2013 62 3 8 038403 Mi Z. Q. Yu Y. Wang Z. Q. Tang J. Q. Preliminary exploration on permanent magnet motor based mechanical elastic energy storage unit and key technical issues Automation of Electric Power Systems 2013 37 1 26 30 Mohamed Y. A. I. Lee T. K. Adaptive self-tuning MTPA vector controller for IPMSM drive system IEEE Transactions on Energy Conversion 2006 21 3 636 644 2-s2.0-33947354917 10.1109/TEC.2006.878243 Jiang S. Ju L. Study on electromechanical coupling nonlinear vibration of flywheel energy storage system Science in China, Series E 2006 49 1 61 77 2-s2.0-33750179257 10.1007/s11431-004-5295-2 Gao Z. W. Kong D. X. Gao C. H. Modeling and control of complex dynamic systems: applied mathematical aspects Journal of Applied Mathematics 2012 2012 5 869792 10.1155/2012/869792 Mahmoodabadi M. J. Bagheri A. Nariman-Zadeh N. Jamali A. Maafi R. A. Pareto design of decoupled sliding-mode controllers for nonlinear systems based on a multiobjective genetic algorithm Journal of Applied Mathematics 2012 2012 22 639014 10.1155/2012/639014 MR2956519 ZBL1254.90214 Kenné G. Goma R. Kwawo H. Lamnabhi-Lagarrigue F. Arzandé A. Vannier J. C. An improved direct feedback linearization technique for transient stability enhancement and voltage regulation of power generators International Journal of Electrical Power and Energy Systems 2010 32 7 809 816 2-s2.0-77955227052 10.1016/j.ijepes.2010.01.018 Abootorabi Z. H. Markadeh G. R. A. Soltani J. Direct torque and flux regulation of synchronous reluctance motor drives based on input-output feedback linearization Energy Conversion and Management 2010 51 1 71 80 2-s2.0-74249111899 10.1016/j.enconman.2009.08.031