Open access peer-reviewed chapter

PID Controller Design Methods for Multi-Mass Resonance System

By Hidehiro Ikeda

Submitted: June 6th 2017Reviewed: January 23rd 2018Published: February 23rd 2018

DOI: 10.5772/intechopen.74298

Downloaded: 625


Motor drive systems are indispensable for applications in the industrial field. High-speed and high-accuracy control is required for motor drive systems. However, solutions to meet these requirements can cause mechanical resonance vibrations to occur in the system as a result of miniaturization and system weight reduction. It is therefore necessary to model these systems as multi-mass resonance systems with multiple masses and finite rigid shafts, gears, and loads. In addition, vibration suppression control should be applied to these systems. This chapter provides two off-line tuning methods for a digital proportional-integral-derivative (PID)-type controller for a two-mass resonance system to suppress its mechanical resonance vibrations. These methods include a coefficient diagram method and a fictitious reference iterative tuning method. The former method uses a nominal mathematical model of the object while the latter method uses only the initial experimental data without use of the mathematical model. In this chapter, the two methods are compared. A controller is proposed that consists of a modified integral-proportional derivative (I-PD) speed controller and a proportional-integral (PI) current controller, and requires no information about the load side state variables. Finally, the effectiveness of the proposed method is confirmed through computer simulations and experimental results.


  • two-mass resonance system
  • vibration suppression control
  • modified I-PD controller
  • coefficient diagram method
  • fictitious reference iterative tuning

1. Introduction

Motor drive systems are used in a wide range of applications, including industrial robots, home electrical appliances, automobiles, steel rolling mills, computers, and space work surfaces. In general, motor drive systems consist of electric motors, gears, belts, flexible shafts, and mechanical load equipment. Recently, the overall stiffness of these systems has been decreasing because of demands for high-speed and high-accuracy system responses, miniaturization, system weight reduction, and low system costs. Additionally, the system constructions have become more complex and the central processing unit (CPU) processing speeds that are required to perform the system calculations have increased exponentially. Consequently, torsional resonance vibrations occur between the motor and the load side. It is therefore necessary to model the system as a multi-mass resonance system, which is composed of several masses with finite rigid shafts, gears, and loads. In addition, a vibration suppression control method should be applied to the system.

The first-order approximation model of the multi-mass resonance system has the form of a two-mass resonance model. Several control methods are effective for control of a two-mass resonance system [15]. PID-type controllers are the most commonly used controllers for industrial applications because of their simplicity and their practicality for use with multi-mass systems. Various PID controller design methods have been proposed; examples of these methods include the limit sensitivity method, the Ziegler and Nichols tuning method, and methods that use the system polynomial. These controller design methods, which are called model-based design methods, may be able to produce the required results in cases where both the system equation and the real system’s parameters are unknown.

In the industrial fields in which many typical motor drive systems are used, experienced technicians often adjust the control system on site to suit the needs of the manufacturing equipment. However, engineer shortages in these fields are becoming a serious problem. It is therefore essential to develop a simple controller design method for industrial applications.

Under these circumstances, and to save both the time required and the cost of tuning the controllers for the motor drive systems, some direct controller tuning methods have been proposed based on the transient response data from closed-loop systems, without modeling of the plant. The fictitious reference iterative tuning (FRIT) method is one of the most promising candidate methods for practical direct parameter tuning [6, 7]. Using the FRIT method, the controller gains can be designed using only single-shot experimental input-output data without knowledge of the model parameters of the object to be controlled.

This chapter introduces two types of approaches to PID-type controller design for suppression of the two-mass resonance system. The first method is based on the assumption that the mathematical model of the object to be controlled is known. This design method is called the coefficient diagram method (CDM) [810]. The CDM is an algebraic approach designed to produce the characteristic polynomial directly in the parameter space. The design of the coefficient diagram of the control system is performed using a differential evolution (DE) procedure to obtain the optimal controller gains in a short time [1113]. The second design method is a FRIT method, which is a PID controller design approach based on one-shot experimental data only and does not use the mathematical model of the object to be controlled [1416]. The effectiveness of the two proposed design methods is confirmed using a combination of computer simulations and experiments.


2. Two-mass speed control system applying PID-type control

2.1. Description of two-mass resonance model

The model consists of two rigid masses and a torsional shaft. The multiple masses on the load side of the actual system are approximated as one inertial element. Similarly, several shafts and gears are approximated as a single torsional shaft. The typical two-mass model is depicted schematically in Figure 1 below [1417].

Figure 1.

Two-mass resonance model.

Here, Jand ωdenote the moment of inertia and the angular speed, respectively, and the suffixes Mand Lindicate the motor side and the load side, respectively. Tinis the input torque, Tdisis the torsional torque, TLis the load torque, and Ksis the shaft stiffness. The continuous state equations of this two-mass resonance model are shown as Eqs. (1) to (3). Additionally, a current loop is considered in this research for high-speed torque control. Eq. (4) is indicative of the voltage equation when using a permanent magnet dc servo motor as the driving motor. In the equations, Ktis the torque constant of the dc motor, Lais the armature current, Rais the total resistance, ucis the control input, Keis the back-electromotive force (back-EMF) constant, and the viscous friction and nonlinear friction sources such as the Coulomb torque are neglected. Therefore, the torque input is calculated using Tin = Ktia.


The research in this case deals with a normalized model to provide generality for the design of the proposed control system. Eqs. (5)(8) show the normalized state equations. The state equation parameters are normalized as shown in Eq. (9), where the suffix puindicates a normalized parameter, K0 [V/pu] is the converter gain, Ka [pu/A] is the current feedback coefficient, Kω[pu/(rad/s)] is the angular speed feedback gain, and τe[s] is an electrical time constant.


Therefore, the unit for all state variables is [pu]. Figure 2 shows a block diagram of the normalized two-mass resonance model.

Figure 2.

Block diagram of normalized two-mass resonance model.

Eq. (10) gives the resonance angular frequency, the anti-resonance angular frequency, and the inertia ratio, respectively.


The next equation gives the simplified transfer function for the two-mass mechanical element, in which the input and the output are iaand ωM, respectively.


The nominal parameters for the two-mass resonance model in this chapter are given in Table 1 below. In this chapter, the proposed CDM method is evaluated through computer simulations, while the proposed FRIT method is evaluated experimentally using the experimental setup shown below.

JM2.744 × 10−4 (kgm2)JL2.940 × 10−4 (kgm2)
Ks18.5 (Nm/rad)Ra2.884 (Ω)
La6.676 (mH)E25.0 (V)
Kt0.2778 (Nm/A)Ke0.2778 (V/(rad sec))

Table 1.

Nominal parameters of the two-mass resonance model.

Figure 3 shows a photograph of the experimental system that was constructed in this research. The two-mass resonance system is simulated using the dc servo motor and a dc generator with a finite rigid coupling. The controller is realized using a digital signal processor that calculates the pulse-width modulation (PWM) signal to send to a four-quadrant dc chopper [17].

Figure 3.

Overview of the experimental system.

The digital signal processor (DSP) board (PE-PRO/F28335 Starter Kit, Myway Plus Corp.), consists of the DSP (TMS320F28335PGFA), a digital input/output (I/O), ABZ counters for the encoder signals, analogue-to-digital (A/D) converters and digital-to-analogue (D/A) converters [18]. The motor and load angles and the angular speeds are detected using 5000 pulses-per-revolution encoders. The dc servo motor current is measured using a current sensor and the A/D converter.

The control period (Ts) and the detection period of the encoder are both 1 ms and the current detection period is 10 μs. While we considered the application of the system to specific apparatus, we then constructed a digital control system that contains a discrete controller. In addition, we used MATLAB/Simulink software to perform the proposed off-line tuning process based on simulations and constructed the PID-type control system as a continuous system [19]. A disturbance is added to the dc generator as a torque using the electric load device on constant current mode. Figure 4 shows the apparatus for the two-mass resonance model that was used in the experimental setup. Figure 5 shows the experimental system configuration [17].

Figure 4.

Photograph of the experimental two-mass resonance model.

Figure 5.

Configuration of the experimental system (for the two-mass resonance model).

2.2. Modified-IPD speed controller and PI current controller

In this chapter, classical PID speed and current controllers are used to suppress the resonance vibrations for the two-mass resonance model. In general, the PI controller, which consists of a proportional controller and an integral controller that are placed in parallel to determine the speed error, is used as the angular speed controller. However, because the resonance system has a complex structure, it is difficult to suppress the vibrations using the classic PI controller alone. Therefore, the I-PD controller is used in this chapter and a first-order lag element is also used to increase the degrees of freedom for the controller design. Additionally, a simple PI controller is used to realize the high-speed torque response for the current minor loop. The continuous control system proposed here is shown in Figure 6. In the figure, ωrefis the reference angular speed, Kp, Ki, Kd, and Tdrepresent the m-IPD speed controller gains, and Kapand Kaiare the PI current controller gains. This chapter proposes two design methods for these six controller gains (i.e., Kp, Ki, Kd, Td, Kap, and Kai). Then, during the simulations and experiments, a digital control system is used, as shown in Figure 7. In this case, the D-control element of the speed controller performs a z-transform in combination with the first lag element of the speed controller to construct a difference equation and avoid the need for a complete differentiation procedure.

Figure 6.

Proposed control system (continuous controller model).

Figure 7.

Proposed control system (discrete controller model).

3. Controller design using coefficient diagram method

3.1. Coefficient diagram method

First, this section explains the coefficient diagram method (CDM) for design of the proposed controller, which is required to suppress the resonance vibrations for the two-mass resonance model. The CDM is an algebraic approach that is used to design the characteristic polynomial directly in the parameter space. In the CDM, a coefficient diagram (CD) is used to perform the controller design. The CD provides the ability to analyze the time response, stability and robustness qualities of the controller using a diagram. In the CD, the vertical axis shows the coefficient of the characteristic polynomial (ai), the stability indices (γi), and the equivalent time constant (τ) logarithmically, while the horizontal axis shows the order ivalues that correspond to the coefficient. Here, the characteristic polynomial is as shown in Eq. (12).


In the CDM, the stability indices (γi), which are defined in Eq. (13) below, are indicators of the stability of the control system.


The equivalent time constant τ, which represents the transient response characteristic, is expressed using the following equation:


The coefficient aican then be calculated using τand the stability indices γias shown in Eq. (15).


In the CDM, use of the standard values of the stability indices is recommended, and these values are listed as follows:


This form is called “the standard form of the stability indices.”

3.2. Design method for the controller gains using the CDM

The CDM is a very simple and effective method for controller design. However, in higher order systems, it is difficult to complete the design by trial and error alone. In this work, the design of the CD for the control system is performed using the differential evolution (DE) method. The DE method is an optimization search method [12, 13].

In the proposed design method, the reference value of the equivalent time constant τrefis first specified. Then, the coefficients of the characteristic polynomial are calculated using the six controller gains with random initial settings. Each of the coefficients is determined using the following equations:


The stability indices are then computed using these calculated coefficients and the evaluation function Fin Eq. (25) is calculated using the terms from Eqs. (26)(30). In this case, the evaluation function Fconsists of an evaluation to match with the set reference value of the equivalent time constant, an evaluation to reduce the change in the next stability index, and an evaluation to match the standard forms of the stability indices γs,i. The weights w1 to w5 of the evaluation functions are set to have values of (w1, w2, w3, w4, w5) = (100, 2, 10, 1, 4), respectively. These steps are subsequently repeated to obtain the optimal controller gains by the DE method.


3.3. Simulation results

Figure 8 shows an example of the simulation results in the form of the angular speed step responses for the control input. The resonance vibrations can be seen in this figure.

Figure 8.

Angular speeds of the step responses of the two-mass model.

Figure 9 shows the frequency response characteristics of the two-mass resonance model from the control input ucto the motor angular speed ωM. The peak point of the mechanical resonance can be observed in this figure. It is therefore essential to construct the controller design method such that it reduces this resonance peak gain.

Figure 9.

Frequency responses of the two-mass resonance model fromuctoωM.

Table 2 shows the controller gains that were designed using the proposed method, where the reference time constant τrefis set to 0.05. Figure 10 shows the designed CD. In the figure, each coefficient is multiplied by 100n, where nis the order of the characteristic polynomial. The figure shows that the form of the diagram is very smooth and that the convex shape is appropriately upward. The designed stability indices are shown in Figure 11. The results in this figure confirm that the stability indices nearly fit the standard form of these indices, and the fluctuations in the numerical values of adjacent indices are also small.

Gain nameValueGain nameValueGain nameValue
Kp2.792 × 102Ki9.007 × 103Kd3.522

Table 2.

Controller gain results when designed using the proposed CDM method.

Figure 10.

Designed coefficient diagram.

Figure 11.

Designed and standard form stability indices.

Figure 12 shows the simulation results, where the speed reference command input ωrefis changed from 0 to 30 rad/s at t = 0 s, and the disturbance torque input is changed from 0 to 20% of the rated torque at t = 0.25 s. The figure indicates that the wave provides a good reference-following performance and illustrates the validity of the vibration suppression characteristics and the disturbance response simultaneously. Additionally, the gain characteristic that was derived using the proposed method over the range from the reference speed ωrefto the motor angular speed ωMis illustrated in Figure 13. The effectiveness of the proposed method is confirmed by the characteristic shown in this figure because the resonance peak is reduced considerably. Therefore, the results for the proposed method show that it is effective as a design method for the vibration suppression controller for the two-mass resonance system.

Figure 12.

Simulation results when using the proposed method.

Figure 13.

Gain characteristics when using the proposed method over the range fromωreftoωM.

4. Controller design using fictitious reference iterative tuning (FRIT)

In this section, an off-line tuning method for the vibration suppression-type speed and current controller gains for the two-mass system is proposed based on the FRIT method; this method uses only single-shot experimental input-output data and does not use either the model parameters or the state equation of the two-mass resonance model. While most FRIT designs only use one state variable, this method uses specific multiple state variables to design the controller gains when using the FRIT method [6, 7, 1416].

4.1. Frit

Figure 14 shows a typical control system, in which Gis the transfer function of the object to be controlled, ris the reference signal, ρrepresents the controller gains, C(ρ) represents the controller, uis the control input parameter, and yis the output parameter. In this case, the mathematical model of Gis not known in advance and is not required for this method.

Figure 14.

Typical closed-loop system.

Initially, as shown in Figure 15, a single-shot experiment is performed using the initial controller gains ρ0, and the control input u0 and output y0 are measured. Then, the reference model M(s), which matches the desired response, is determined. A fictitious reference signal is then generated using the controller, the control input u0, and the output y0, as shown in Eq. (31) below. This means that the initial data u0 and y0 can be obtained using any value of ρif r˜ρis input to the closed-loop system used to implement C(ρ).


Figure 15.

Measurement of the initial data.

The optimal controller gains that are required to achieve yM = y0 are then determined, as shown in Figure 16, using an optimization search method. Finally, these controller gains then represent the best available solutions that allow the desired control system response to be obtained. Therefore, the controller design process can be performed without any prior information about either the model parameters or the state equations.

Figure 16.

Reference model showing the input of the fictitious reference signal.

4.2. Vibration suppression controller design method by FRIT

Figure 17 shows a simplified form of the proposed vibration suppression control system, where Cω1, Cω2, Cω3, and Ciare the controllers. While the FRIT method generally uses one control input and one state variable for the initial experimental data, the proposed method uses one control input, uc, plus two state variables, ωMand ia. The controller gain vector ρis given as follows:


Figure 17.

Simplified proposed control system.

Therefore, the fictitious reference signal ω˜refρcan be calculated using the following equation without any need for the two-mass resonance model.


Here, uc0, ωM0 and ia0 represent the initial experimental data. The reference model M(s), as shown in Eq. (34), is then used depending on the purpose of the system, where the time constant τis a reference model parameter.


Here, τis calculated using the following equation with the 99% response time parameter T99.


The differential evolution method is then used to search for the optimal gains. The performance index function Fis then defined as shown in Eq. (36) below using ωM0 and yM=Msω˜refρ.


4.3. Experimental results

Figure 18 shows an example of the experimental results in the form of the angular speed step response of the control input. The resonance vibrations can again be observed in a similar manner to the case of the simulation results shown in Figure 8 above.

Figure 18.

Angular speeds of the step responses (DC voltage input).

Figure 19 shows the gain characteristics of the frequency responses from the experimental results shown in Figure 18, which relate the input voltage to the motor angular speed ωMand the load angular speed ωL. From these characteristics, the peak resonance vibrations at approximately 300 rad/s are also observed. These results were calculated based on the experimental input and output waves of the voltage step response using the method that was proposed in [20]. Additionally, both the resonance and anti-resonance points can be found in these figures.

Figure 19.

Calculated frequency responses to dc voltage input (left:uctoωM, right:uctoωL).

Figure 20 shows the experimental results obtained when using the general PI speed and current controller for comparison with the effects of the proposed control system. Figure 21 shows the gain characteristics for the frequency responses shown in Figure 20, which relate ωrefto ωMand ωL, where ωrefis 30 rad/s. These characteristics show that the peak gain of the resonance is not greatly attenuated. Figures 2224 show the initial experimental waves for ωM0, ia0 and uc0, respectively, that were obtained using the values of the initial controller gain ρ0, which are listed as follows:


Figure 20.

Experimental results obtained using conventional PI controller.

Figure 21.

Frequency responses of conventional PI speed control system (left:ωreftoωM; right:ωreftoωL).

Figure 22.

InitialωM0 data.

Figure 23.

Initialia0 data.

Figure 24.

Initialuc0 data.

Both the initial rise and the oscillation can be observed in these figures.

Table 3 shows the results for the controller gains determined using the proposed FRIT design method with searching by the DE method, where T99 was set at 0.2 s. Figure 25 shows a comparison of the experimental results obtained using the controller gains that were designed using the proposed method with the simulated results for the reference output yM.The results in the figure show that the proposed off-line tuning method works very well, despite the fact that the design was performed using the initial one-shot experimental data alone. Figure 26 shows the experimental results that were obtained for ωM, ωL, and iawhen the proposed FRIT method was used. Here, ωrefis stepped from 0 to 30 rad/s when tis 0 s. The figure shows the good response of the proposed vibration suppression speed controller. Figure 27 shows the experimental results, where ωrefis stepped from 30 to 50 rad/s when tis 0 s and the disturbance torque is increased from 0 to 10% of the rated torque when tis 0.5 s. As these figures show, good waves were observed in terms of their reference-following performance and disturbance response.

Gain nameValueGain nameValueGain nameValue
Kp1.24Ki17.0Kd6.53 × 10−3
T8.59 × 10−4Kap2.37Kai135

Table 3.

Controller gain results when designed using the proposed FRIT method.

Figure 25.

Initial experimental results forωM0 and simulation results foryM.

Figure 26.

Experimental results forωM,ωL, andiaobtained when the proposed FRIT method.

Figure 27.

Experimental results for speed step response and disturbance response.

Figure 28 shows the gain characteristics of the frequency responses of the proposed control system, where these characteristics are shown from the perspectives of ωrefrelative to ωMand ωL. The resonance vibration suppression effect can be observed in these figures. Therefore, the effectiveness of the proposed control system and the design method based on use of the FRIT method can be confirmed. Additionally, Figure 29 shows the experimental results (ωL) that were obtained for various values of the speed reference time parameter, where T99 = 0.15, 0.175, 0.2, 0.25, 0.3, and 0.35. As shown in this figure, the response times change satisfactorily and the proposed design method for the controller gains can thus also be used to design the response times arbitrarily.

Figure 28.

Frequency responses of the proposed control system (left:ωreftoωM; right:ωreftoωL).

Figure 29.

Experimental responses ofωLto various values ofT99.

5. Conclusion

This chapter has proposed two design methods for the controller gains required for vibration suppression control in a two-mass resonance system. The proposed controller consists of a modified-IPD speed controller and a PI current controller. The proposed controller design methods are based on application of the coefficient diagram method (CDM) and application of the fictitious reference iterative tuning (FRIT) method. Both methods use the motor side variables only, including the motor’s angular speed and the armature current. The CDM method uses the coefficient of the characteristic polynomial of the control system and can determine the control performance based on the shape of the coefficient diagram and the stability indices. In this chapter, the fitting performances for the standard form of the stability indices and the coefficient values were used to determine the controller gains, which were designed using the differential evolution method. The effectiveness of the proposed CDM was confirmed by the simulation results. The FRIT method can be used to design the controller without knowledge of the model state equations and their parameters. Furthermore, a fictitious signal that was calculated using the initial experimental data for multi-state variables was also proposed in this chapter. The effectiveness of the proposed FRIT method was confirmed using the experimental results. Consequently, the CDM and the FRIT method were shown to produce the same design performance. The CDM is useful for controller design when the mathematical model and the object parameters are known. The FRIT method is effective when the mathematical model is unknown but the initial experimental data can be observed.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Hidehiro Ikeda (February 23rd 2018). PID Controller Design Methods for Multi-Mass Resonance System, PID Control for Industrial Processes, Mohammad Shamsuzzoha, IntechOpen, DOI: 10.5772/intechopen.74298. Available from:

chapter statistics

625total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Data-Based Tuning of PID Controllers: A Combined Model- Reference and VRFT Method

By Jyh-Cheng Jeng

Related Book

First chapter

Wavelet PID and Wavenet PID: Theory and Applications

By José Alberto Cruz Tolentino, Alejandro Jarillo Silva, Luis Enrique Ramos Velasco and Omar Arturo Domínguez Ramírez

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us