Published: 31 March 2014

Discrete dynamics analysis for nonlinear collocated multivariable mass-damper-spring intelligent mechanical vibration systems

Cheng Zeng1
Shan Liang2
Yongsheng Sun3
Yingying Su4
1, 2, 3, 4College of Automation, Chongqing University, Chongqing, 400044, China
Corresponding Author:
Shan Liang
Views 154
Reads 68
Downloads 1185

Abstract

A new time-discretization method for the development of nonlinear collocated multivariable mass-damper-spring (MDS) intelligent mechanical vibration systems is proposed. It is based on the Runge-Kutta series expansion method and zero-order hold assumption. In this paper, we show that the mathematical structure of the new discretization scheme is explored and characterized in order to represent the discrete dynamics properties for nonlinear collocated multivariable MDS intelligent mechanical vibration systems. In particular, the decent effects of the time-discretization method on key properties of nonlinear multivariable MDS mechanical vibration systems, such as discrete zero dynamics and asymptotic stability, are examined. The resulting time-discretization provides discrete dynamics behavior for nonlinear MDS mechanical vibration systems, which enabling the application of existing controller design techniques. The ideas presented here generalize well-known results from the linear case to nonlinear plants.

1. Introduction

It is generally known that the classical mass-damper-spring (MDS) intelligent mechanical vibration systems have been studied by some researchers [1-3]. The MDS system is a simple mechanical vibration model, and is widely applied in many fields [4-9]. However, to analyze or utilize such systems in practice invariably requires discretization. Thus, the discrete-time plant or sampled-data model play an important role in the intelligent mechanical vibration systems. For example, it is well known that modern controllers and signal processing devices invariably operate in discrete time, in particular, employing digital technology.

In the linear case, some famous scholars, such as Ishitobi [1], Lin [2, 3] and Zhu [5, 6], have represented the discrete dynamics behavior of the MDS vibration systems. In this context, topics such as discrete-time plant (normal form), relative degree and degree-of-freedom (DOF), chaos and bifurcations, and zero-pole have been definitely shown. And these ideas presented generalize well-known notion from the continuous-time case to discrete-time and sampled systems, although the theory for the discrete-time case is less well developed than for the continuous-time case, both linear and nonlinear systems.

One would reasonably expect similar results to hold for nonlinear MDS systems. However, the situation for the nonlinear case is more complex than for linear case. Indeed, to the best of our knowledge, the discrete dynamics properties for nonlinear MDS systems are still open problems. Moreover, we feel deeply that this research is well very important, especially for discrete dynamics analysis of intelligent mechanical vibration systems and their applications.

Zero dynamics of the nonlinear MDS plant, corresponding to the zero of linear case, are fundamental characteristics of nonlinear intelligent mechanical vibration systems [10]. The occurrence of nonlinear zero dynamics is relevant to the discrete dynamics behavior of nonlinear intelligent mechanical systems. However, the previous results cannot be applied to the nonlinear MDS mechanical vibration systems. Hence, it is an important research topic to find criteria which guarantee that the discrete dynamics properties of nonlinear MDS systems are obtained.

The objective of this paper is to analyze the discrete dynamics behavior of nonlinear MDS intelligent mechanical vibration systems with collocated sensors and actuators, irrespectively of whether these systems are single-input single-output (SISO) plant or multivariable case. Our aim is only to analyze the multivariable MDS vibration systems because the SISO plant is a particular example of multivariable case. More importantly, we show how a particular strategy can be used to present the discrete dynamics properties, which is accurate to some order in the sampling period. An insightful interpretation of the obtained discrete dynamics behavior can be made in terms of additional discrete zero dynamics. Moreover, the current paper extends the well-known idea from the linear case to nonlinear systems, which including extra discrete zeros or zero dynamics due to the sampling process. Finally, numerical experimental example is given to illustrate the theories developed in this paper.

2. Discrete-time models of nonlinear collocated multivariable vibration MDS systems

The dynamics of such intelligent mechanical vibration system can be described by Newton’s Second Law, and further consider an n-mode, m-input, m-output nonlinear collocated multivariable MDS intelligent mechanical vibration system described by [11]:

1
my¨+Ff+Fsp=F,

where Ff is the frictional damping force, Fsp is the spring dynamics force and F is the force acting on the mass. In order to reduce the computation complexity of such mechanical system, and in particular, without loss of generality, the two-DOF (2-DOF) MDS mechanical vibration system is primarily considered, and also depicted in Figure 1. More importantly, 2-DOF MDS mechanical vibration plant has been mainly employed as a research subject since it is used most commonly in practice. In particular, a 2-DOF mechanical vibration MDS model, which constitutes a slider, spring, damping components and a pendulum is consider, as shown in Fig. 1. Therefore, the whole idea is given for the case of a two-input-two-output (2-DOF) mechanical vibration system for simplicity of description.

Fig. 1Multivariable 2-DOF MDS intelligent mechanical vibration model

Multivariable 2-DOF MDS intelligent mechanical vibration model

In this case, we can represent the discrete dynamics behavior of nonlinear 2-DOF MDS mechanical vibration model, which discretized zero dynamics have been also discussed owing to their significance. Furthermore, it is straightforward to extend the particular case (2-DOF) to the general case, such as multi-DOF.

The discrete dynamics analysis of nonlinear 2-DOF MDS mechanical vibration system is derived below. First, a 2-input and 2-output nonlinear multivariable MDS mechanical vibration system can be expressed in its so-called normal form [11]:

2
z˙11=z21,z˙21=-k1m1z11-z12-k-1m1z11-z123-c1m1z21-z22+u1m1,z˙12=z,z˙22=k1m2z11-z12+k-1m2z11-z123+c1m2z21-z22 -k2m2z11-z12-k-2m2z11-z123-c2m2z21-z22+u2m2,y1=z11.

Remark 1: It is obvious that the relative degrees r1, r2 of the system (2) are two. In fact, the relative degrees of many nonlinear mechanical vibration systems in the practical field are two.

We are interested in the discrete-time plant or sampled-data model for the nonlinear MDS mechanical vibration system (2) when the input is a piecewise constant signal generated by a zero-order hold (ZOH); i.e.:

3
ut=ukT, kTt<k+1T, k=0,1,,

where T is a sampling period.

For small sampling periods, a more accurate sampled-data model of nonlinear 2-DOF multivariable MDS mechanical vibration system is definitely obtained when applying a higher-order Runge-Kutta (RK) expansion such as:

4
zl+1,k+1i=yi,k+1lyi,kl+Tyi,kl+1+T22!yi,kl+2+Tri-lri-l!yi,kri+Tri-l+1ri-l+1!yi,kri+1,
i=1,,m, l=0,,ri-1.

Remark 2: Consider here a higher-order RK expansion (4) instead of the uniformly valid asymptotic expansion with the continuous-time system relative degrees ri. The reason is that the closed-loop system becomes unstable in the case of the latter when a discrete-time controller design method based on the assumption of the stability of the zero dynamics is applied [12]. The sampled-data model cannot be used for discrete-time controller design because of the existence of the unstable sampling zero dynamics. More precisely, when the input is imposed through a ZOH on the original continuous-time system, the output does not converge to the origin.

Before proceeding, the following assumption is needed here for the preservation of an affine property in the process of sampling.

Assumption 1: The unique equilibrium point lies on the origin.

Next, note here that:

5
u˙it=0, u¨it=0, tkT,k+1T, i=1, 2,

and that:

6
y˙i=z˙1i=z2i,y¨i=z˙2i=bi+aiui,yi(3)=b˙i+a˙iui=b-i+a-iui,

where i= 1, 2,

b1=-k1m1z11-z12-k-1m1z11-z123-c1m1z21-z22,
b2=k1m2z11-z12+k-1m2z11-z123+c1m2z21-z22-k2m2z11-z12-k-2m2z11-z123-c2m2z21-z22,
a1=1m1, a2=1m2,

and

7
b-i=biz1iz2i+biz2ibi+aiz2ibi+l=3nbizlicl, a-i=aiz1iz2i+biz2iai+aiz2iai+l=3naizlicl.

Thus, a more accurate discrete-time model of nonlinear 2-DOF multivariable MDS mechanical vibration system is precisely derived by substituting (6) into the RK expansion forms of yi((k+1)T) and y˙i((k+1)T) and neglecting higher order terms as follows:

8
z1,k+11=z1,k1+Tz˙1,k1+T22!b1,k+a1,ku1,k+T33!b-1,k+a-1,ku1,k,z2,k+11=z˙1,k1+Tb1,k+a1,ku1,k+T22!b-1,k+a-1,ku1,k,z1,k+12=z1,k2+Tz˙1,k2+T22!b2,k+a2,ku2,k+T33!b-2,k+a-2,ku2,k,z2,k+12=z˙1,k2+Tb2,k+a2,ku2,k+T22!b-2,k+a-2,ku2,k,y1,k=z1,k1,y2,k=z1,k2.

We next analyze the local and global truncation errors between the true system output and the ith output of the obtained sampled-data model (8), assuming that, at t=KT, the state z is equal to the true system state z(kT). We compare the true system output y(kT+T) with the first (third) state of the approximate sampled-data model in (8) at the end of the sampling interval. This yields the following local truncation output error:

9
yikT+T-qz1,ki==z1,ikT-z1,ki+Tz2,ikT-z2,ki +T22!bi,k+ai,kui,kt=ξ1-bi,k+ai,kui,kt=kT+T33!b-i,k+a-i,kui,kt=ξ1-b-i,k+a-i,kui,kt=kT.

From the assumptions used in the local truncation error we have [13]:

10
z1,i0-z1,ki0OT4, Tz2,i0-z2,ki0OT4.

For the last two terms in (9), since the control input is generated by a ZOH, we know that:

11
ei=T22!bi,kς,η+ai,kς,ηui,kt=ξ1-bi,kς,η+ai,kς,ηui,kt=kT+T33!b-i,kς,η+a-i,kς,ηui,kt=ξ1-b-i,kς,η+a-i,kς,ηui,kt=kTT22!L1zξ1-zkT+T33!L2zξ1-zkT,

where L1 and L2 are Lipschitz constants. Furthermore, according to [13], the Lipschitz condition guarantees that the variation of the state trajectory z(t) can be bounded as:

12
zξ1-zkTCeL2ξ1-kT-1L2<CeL2T-1L2OT.

Therefore, the local truncation error between the ith output of the resulting sampled-data model and the corresponding true continuous-time system output is of order T4 by following the equation (11).

In the following section, the focus is on the global truncation error of the approximate discrete-time model for nonlinear 2-DOF multivariable MDS mechanical vibration system. The proof for the global truncation error is completely analogous to the proof for the local truncation error. In particular, the key difference between the local and global truncation errors arise from the fact that the global truncation error is associated with the number of steps N, and has practical relevance to applications in which the asymptotic discrete dynamics behaviors of the proposed model is of interest. Consider:

13
θk+1=Aθk+ϕE,

where:

A=1T01, ϕ=T22!+T33!T+T22!, θk=yikT+T-qz1,kiy˙ikT+T-qz2,ki,
E=(bi,k(ς,η)+ai,k(ς,η)ui,k)t=ξ1-bi,kς,η+ai,kς,ηui,kt=kT.

Then, after N steps, we have:

14
θk+N=ANθk+i=0N-1AiϕE.

As shown in the proof of the local truncation error, EO(T), and by simple straightforward calculation, every element in the last column of the matrix i=0N-1Ai is of order O(1). In addition, given that the last row of ϕ is T+T22! such that there will always be a term of order O(T2) for every state, and note that all the elements of AN are of order O(1). So the global truncation error is obvious that:

15
i=0N-1AiϕEOT3θk+NOT3.

On the other hand, the known Euler integration leads to a model having local truncation error T2 and global truncation error T. The local and global truncation errors between the ith output of the resulting sampled-data model (8) and the corresponding true continuous-time output are of order T4 and T3 respectively. This fact means that the approximate sampled-data model (8) is closer to the true system than that by Euler.

Remark 3: A key result from nonlinear systems theory is that a nonlinear system of degree and uniform relative degree can be described by a model in normal form. The approximate sampled-data model of interest here is obtained by using a truncated RK series expansion of the states in normal form. Moreover, it has been shown that the model possesses previously unknown properties, regarding the propagation of truncation errors.

Remark 4: The Euler model is well-known as an approximate model, it is difficult to obtain good discrete dynamics performance because of the poor approximation. An insightful observation is that the proposed model (8) is more accurate than the Euler model in terms of the discretization of nonlinear 2-DOF multivariable MDS mechanical vibration system. Therefore, a controller design by the Euler model is easier but better discrete dynamics performance could be obtained if the discrete-time model (8) is used for controller design.

3. Discrete dynamics analysis of nonlinear collocated multivariable 2-DOF MDS intelligent mechanical vibration systems

In this section, we obtain the discrete dynamics analysis of nonlinear collocated multivariable 2-DOF MDS intelligent mechanical vibration system. First, we present the result that shows the discrete zero dynamics, which play an important role in the area of discrete dynamics analysis. Substituting z1,k+1i=z1,ki=0 into (8) yields:

16
0=Tz˙1,k1+T22!b1,k0+a1,k0u1,k+T33!b-1,k0+a-1,k0u1,k,z2,k+11=z˙1,k1+Tb1,k0+a1,k0u1,k+T22!b-1,k0+a-1,k0u1,k,0=Tz˙1,k2+T22!b2,k0+a2,k0u2,k+T33!b-2,k0+a-2,k0u2,k,z2,k+12=z˙1,k2+Tb2,k0+a2,k0u2,k+T22!b-2,k0+a-2,k0u2,k.

Here, ai,k0 and bi,k0 denote the values of ai,k and bi,k, respectively. Therefore, noting that constant terms bi,k0 is definite higher order with respect to T than the corresponding one of variable terms zi,k. In addition, when applying z-transform to (16), then it leads to:

17
Φ(z)v=0,
18
Φz=z-10Ta1,k0+T22a-1,k000z-10Ta2,k0+T22a-2,k0-T0T22a1,k0+T36a-1,k000-T0T22a2,k0+T36a-2,k0, v=z2,k1z2,k2u1,ku2,k.

As a result, the sampled counterpart of the discrete zero dynamics (sampling zero dynamics) are obtained from |Φ(z)|=0.

Next, the matrix S is defined as:

19
S=10-2T0010-2T00100001.

Multiplying (18) by the matrix S from the left hand side leads to:

20
SΦz=z+10T26a-1,k000z+10T26a-2,k0-T0T22a1,k0+T36a-1,k000-T0T22a2,k0+T36a-2,k0.

Using the computational formula for determinant of block matrix:

21
EFGH=HE-FH-1G.

We can obtain:

22
SΦz=T22a1,k0+T36a-1,k000T22a2,k0+T36a-2,k0
×z+100z+1+T36a-1,k000a-2,k0T22a1,k0+T36a-1,k000T22a2,k0+T36a-2,k0-1.

Further, it is obvious that the discrete zero dynamics of nonlinear collocated multivariable 2-DOF MDS intelligent mechanical vibration system are given by |SΦ(z)|=0. In particular, these discrete zero dynamics are mainly determined in the second half of (22) due to the existence of variable z. Thus:

23
SΦz=z+100z+1+T36a-1,k000a-2,k0T22a1,k0+T36a-1,k000T22a2,k0+T36a-2,k0-1=z+1+a-1,k3a1,k+Ta-1,kT00z+1+a-2,k3a2,k+Ta-2,kT=0.

Remark 5: The discrete zero dynamics of nonlinear sampled-data collocated multivariable 2-DOF MDS intelligent mechanical vibration model are classified into two categories. The former ones have counterpart in the underlying continuous-time system and they are called intrinsic zero dynamics. The latter ones are generated in the sampling process and they are called sampling zero dynamics. More importantly, we are focused primarily on the latter ones in these two kinds of the discrete zero dynamics in terms of the discrete dynamics analysis of nonlinear collocated multivariable 2-DOF MDS intelligent mechanical vibration system.

Finally, we consider the stability of nonlinear discrete zero dynamics because of the importance of discrete zero dynamics in the nonlinear collocated multivariable 2-DOF MDS intelligent mechanical vibration system. The determinant (23) can be expanded as follows when supposing β1=z+1+a-1,k3a1,k+Ta-1,k, β2=Tz+1+a-2,k3a2,k+Ta-2,kT:

24
(z+1)2+T(β1+β2)(z+1)+β1β2T2=0.

When we perform the bilinear transformation z=λ+1λ-1 on the above formula, the equation is written as:

25
[4+2T(β1+β2)+T2β1β2]λ2-[2T(β1+β2)+2T2β1β2]λ+T2β1β2=0.

It is clear that the two roots of (25) lie in the open left half of λ-plane if:

26
-2Tβ1+β2+2T2β1β24+2Tβ1+β2+T2β1β2-Tβ1+β22>0,
27
T2β1β24+2T(β1+β2)+T2β1β2β1β24>0.

Furthermore, for sufficiently small sampling periods, the discrete zero dynamics of the nonlinear collocated multivariable 2-DOF MDS intelligent mechanical vibration system are stable if:

28
β1β2(β1+β2)<0.

Remark 6: Discrete dynamics behaviors of nonlinear collocated multivariable MDS intelligent mechanical vibration systems, like in the 2-DOF analysis above, can be readily extended to more complex MDS models such as three-DOF (3-DOF) (see also Fig. 2.), and so on. On the other hand, 2-DOF case also reduces a single DOF immediately.

Fig. 2Multivariable 3-DOF MDS intelligent mechanical vibration model

Multivariable 3-DOF MDS intelligent mechanical vibration model

4. Numerical experimental results of nonlinear discrete dynamics analysis

In the above discussion, we have definitely focused on the discrete dynamics analysis of the nonlinear collocated multivariable 2-DOF MDS intelligent mechanical vibration systems, and the performance of the proposed time-discretization method for nonlinear input-driven systems by using RK approach is evaluated by applying it to a nonlinear collocated multivariable MDS mechanical vibration system. In particular, discrete-time open-loop responses to a MDS mechanical vibration system are considered. In this paper, the discrete dynamics exact or approximate values obtained at every time step using the proposed time-discretization method are represented to the values obtained using the MATLAB solver and RK method at the corresponding time steps.

Fig. 3State response of the 2-DOF mechanical MDS vibration system

State response of the 2-DOF mechanical MDS vibration system
State response of the 2-DOF mechanical MDS vibration system
State response of the 2-DOF mechanical MDS vibration system
State response of the 2-DOF mechanical MDS vibration system

This numerical experiment is simulated under various combinations of sampling periods and the local truncation error by using MATLAB. Acceptably accurate results are obtained when the sampling period is T= 0.01 with a local truncation error. If the sampling period is further reduced, then more accurate results could be obtained. Results obtained from the MATLAB solver and the RK series discretization method are shown in Fig. 3 when the sampling period is T= 0.001. Thus, Figure 3 shows the state response of the nonlinear collocated multivariable 2-DOF MDS mechanical vibration system by using this method. Additionally, this special system clearly presents that the proposed RK series expansion method discretizes a nonlinear MDS intelligent mechanical vibration system quite accurately.

Next, we design a discrete-time model following controller on the basis of the nonlinear 2-DOF MDS mechanical multivariable discrete-time vibration model with ZOH, and apply it to the original continuous-time MDS system to further analyze the discrete dynamics behavior of the nonlinear collocated multivariable 2-DOF MDS mechanical vibration system. It is obvious that the desired tracking objective of the output error results can be readily approached by generated from the plant control input sequences in simulation (see also Figure 4 and 5).

In addition, the good performance of the discrete dynamics behavior for the nonlinear collocated multivariable 2-DOF MDS mechanical vibration controlled system can be achieved because the stability of the discrete zero dynamics of nonlinear multivariable 2-DOF MDS mechanical vibration system are significantly improved in the Figure 6.

Fig. 4Inputs of model following control for the 2-DOF MDS mechanical vibration system

Inputs of model following control for the 2-DOF MDS mechanical vibration system
Inputs of model following control for the 2-DOF MDS mechanical vibration system

Fig. 5Error between the desired outputs and the true plant outputs for the 2-DOF MDS mechanical vibration system

Error between the desired outputs and the true plant outputs for the 2-DOF MDS mechanical vibration system
Error between the desired outputs and the true plant outputs for the 2-DOF MDS mechanical vibration system

Finally, as an example of such application, these above results about nonlinear collocated multivariable 2-DOF MDS mechanical vibration controlled system can be also applied to 2-DOF half-vehicle suspension model with nonlinear springs and dampers shown in Figure 7.

Fig. 6The magnitudes of discrete zero dynamics for nonlinear MDS mechanical vibration system

The magnitudes of discrete zero dynamics for nonlinear MDS mechanical vibration system

Based on the above analysis, we can now utilize the MTATLAB to build up the discrete dynamics behavior of the nonlinear 2-DOF half-vehicle suspension model. The corresponding response curves of acceleration signal in time domain with different amplitude of input are shown in Figure 8.

As can be seen from the Figure 8, as soon as pulse signals begin, the acceleration signals change acutely and further decay rapidly under the action of damping, though low frequency oscillation will still maintain a period of time. With the input signals amplitude increasing, the amplitude of the output acceleration signals also increase accordingly, and vice versa.

Fig. 7Nonlinear 2-DOF half-vehicle suspension model

Nonlinear 2-DOF half-vehicle suspension model

Fig. 8The corresponding response curve of acceleration signal with different amplitude of input

The corresponding response curve of acceleration signal with different amplitude of input

(a) The response curve with the input of +0.1

The corresponding response curve of acceleration signal with different amplitude of input

(b) The response curve with the input of –0.1

5. Conclusions

This paper has developed a new approach to propose the discrete dynamics behavior of nonlinear collocated multivariable MDS intelligent mechanical vibration systems. It is based on the RK series expansion method which uses a more sophisticated derivative approximation than the simple Euler approach. Moreover, an insightful interpretation is given in terms of an explicit characterization of the nonlinear sampling zero dynamics of the obtained discrete-time models, which can be used effectively to analyze and process the discrete dynamics properties. An application is given here for the nonlinear 2-DOF half-vehicle suspension model with springs and dampers. This extends the well-known results for discrete dynamics behavior of linear systems to the nonlinear case. These results are believed to give important insights which are relevant to many aspects of nonlinear intelligent mechanical vibration systems theory.

References

  • M. Ishitobi, S. Liang On zeros of discrete-time models for collocated mass-damper-spring systems. Journal of Dynamics Systems, Measurement and Control, Vol. 3, 2004, p. 205-209.
  • J. L. Lin On transmission zeros of mass-dashpot-spring systems. ASME Journal of Dynamic Systems, Measurement and Control, Vol. 121, 1999, p. 179-183.
  • G. S. Calafiore, S. Carabelli, B. Bona Structural interpretation of transmission zeros for matrix second-order systems. Automatic, Vol. 33, 1997, p. 745-748.
  • J. L. Lin, J. N. Juang Sufficient conditions for minimum-phase second-order linear systems. Journal of Vibration and Control, Vol. 1, 1995, p. 183-199.
  • Q. Zhu, M. Ishitobi Chaos and bifurcations in a nonlinear vehicle model. Journal of Sound and Vibration, Vol. 275, 2004, p. 1136-1146.
  • Q. Zhu, M. Ishitobi Chaotic vibration of a nonlinear full-vehicle model. International Journal of Solids and Structures, Vol. 43, 2006, p. 747-759.
  • S. Liang, C. Li, Q. Zhu, Q. Xiong The influence of parameters of consecutive speed control humps on the chaotic vibration of a 2-DOF nonlinear vehicle model. Journal of Vibroengineering, Vol. 13, 2011, p. 406-413.
  • A. Sezgin, Y. Arslan Analysis of the vertical vibration effects on ride comfort of vehicle driver. Journal of Vibroengineering, Vol. 14, 2012, p. 559-571.
  • S. Liang, Y. Sun, Q. Zhu, Z. Yang, C. Zeng Ride comfort analysis of a nonlinear vehicle excited by the consecutive speed-control humps. Journal of Vibroengineering, Vol. 15, 2013, p. 1668-1676.
  • J. Yuz Sampled-data models for linear and nonlinear systems. School of Electric Engineering and Computer Science, University of Newcastle, Ph. D. dissertation, NSW, Australia, 2005.
  • H. Khalil Nonlinear systems. 3rd ed., Upper Saddle River, NJ, Pretice-Hall, 2002.
  • M. Ishitobi, M. Nishi, S. Kunimatsu Stability of zero dynamics of sampled-data nonlinear systems. Preprints of the 17th IFAC World Congress, 2008, p. 5969-5973.
  • J. Butcher The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods. Wiley-Interscience New York, NJ, USA, 1987.

About this article

Received
17 September 2013
Accepted
21 December 2013
Published
31 March 2014
Keywords
discrete dynamics
nonlinear MDS vibration system
discrete-time model
sampling zero dynamics
Acknowledgements

The authors would like to thank anonymous referees for their helpful comments and suggestions, and this work is partially supported by the National Basic Research Program of China (2013CB328903), the National Natural Science Foundation of China (60574003), the Natural Science Foundation Project of CQ CSTC (cstc2012jjA40026) and Project No. CDJXS12170006 Supported by the Fundamental Research Funds for the Central University.