Abstract

This paper studies dynamics of a three-column centrifuge and tries to design a ball balancer to suppress its vibration and increase its operating speed. At first, two vibration models with and without the balancer are introduced. Then modal analysis of the centrifuge itself is carried out with two isolated natural frequency groups identified. Centrifugal forces acting on the system at the steady state are analyzed and a tip-over phenomenon is discovered and validated through simulations. A ball balancer is adopted in the centrifuge and several potential tip-over phenomena related to the balancer are studied. At last, a comparison between dynamics of the centrifuge before and after the balancer is installed is made with benefits and side effects of the balancer presented.

1. Introduction

A three-column centrifuge is named due to its three supporting legs. Up until now, it has been widely employed in industries of textile, food, pharmacy, chemistry, metallurgy, and so forth, for the purpose of dehydration, incrassation, separation, classification of solid particles, and so forth.

Researches on the design of the three-column centrifuge can be found as early as the 1960s. For satisfying organization assembly line production at chemical plants, Gerke and Noskov [1] discussed the shearing off and discharge of the centrifuged cake. Then they improved the centrifuge by employing a hydraulic motor [2]. In 1968, a continuous centrifuge in the metal working industry was introduced by Il’yashenko and Gerke [3] for continuous deoiling of metal chips contaminated by oil.

Some early studies on dynamics of the centrifuge can be found in the late 1960s. A linear vibration model was built by Likuchev [4] assuming that vibration amplitude of the centrifuge is quite small, and then frequencies of the centrifuge were studied. Later, the rigidity of several types of supports was discussed in [5]. The steady motion of the centrifuge was investigated by Vanichkin [6] with the reaction force of the suspension system calculated. However, most of these studies were based on linear models, and nonlinear dynamic characteristics of the centrifuge were not presented. Serious vibrations often happen once the centrifuge operates above a certain speed. Up until now, not much work has been done on the reason and a solution to this is still unclear.

The suspension structure of the centrifuge is quite similar to that of the vertical axis washing machine with a pendulant suspension system [711]. Thus, the balancer demonstrated in the washer may also be employed in the centrifuge. However, differences between the centrifuge and the washer have to be considered. The biggest difference lies in the stiffness of their suspension systems. The centrifuge employs rigid suspension rods while the washer uses soft springs. Thus, the suspension system of the centrifuge is much stiffer than that of the washer. Although this makes complete balance of the centrifugal forces much easier, side effects could be induced such as the tip-over phenomenon discovered in this paper.

The balancer has been studied for many years. When the rotational speed of a rotor is much higher than its critical frequency, eccentricity of the system can be partially or sometimes totally counteracted by the balancer. However, disastrous consequences can still happen even after the eccentricity is eliminated, especially when the eccentricity and the balancer are in different planes. The three-column centrifuge provides a very good example of such condition. The tip-over phenomenon can still occur even after the elimination of the eccentricity.

The structure of this paper is as follows. Vibration models of the centrifuge were introduced in Section 2. Then dynamics of the centrifuge itself are studied in Section 3.1. Influences of the balancer on dynamics of the centrifuge are discussed in Section 3.2. And a comparison between dynamics of the centrifuge before and after the balancer is installed is made in Section 3.3. At last, some conclusions are drawn in Section 4.

2. Vibration Models of the Three-Column Centrifuge

As can be seen from Figure 1(a), a three-column centrifuge mainly comprises a core assembly and a suspension system. The assembly is composed of a bowl, an outer shell, a supporting bed, and a motor. The motor and the outer shell are fixed on the supporting bed, while the bowl is installed on the bed through a revolving pivot. Figure 1(b) shows the suspension structure of the centrifuge; as can be seen, the supporting bed is connected to the supporting legs though three suspension rods with spherical joints at both ends. A compression spring is inserted between each leg and the supporting bed.

For describing movements of the core assembly, two reference frames are constructed in Figure 1: a global one on the ground and a dynamic one in the core assembly.

The structure of the three-column centrifuge is quite similar to that of the vertical axis washing machine discussed in [711]. And the vibration model of the centrifuge can be similarly constructed using Lagrange’s equation based on the kinetic and geopotential energies derived in [8, 11]; thus, only the final forms of the vibration models are summarized here. Two models with and without the balancer are built in the paper. The one without balancer can be expressed aswhere with denoting the translational displacement of (the origin of the dynamic reference frame ) and being the Cardan angles; is a vector relative to the rotational speed and rotational acceleration of the motor; is the generalized forces provided by the suspension system; and are the mass matrix and geopotential energy of the core assembly. A detailed expression of the above model can be found in [11].

For suppressing vibrations of the centrifuge, a balancer composed of a circular cavity and several balls is employed. The cavity is fixed on the bowl and the balls can run freely in the cavity. For exploring basic characteristics of the dynamic system, the balancer is assumed to be ideally manufactured and lubricated, and only viscous damping exists between the balls and the cavity; besides, a two-ball model is employed and the idea that no collision happens between the balls is assumed. Suppose that the rotational angle of the th ball relative to the cavity (or bowl) is ; the vibration model of the centrifuge with a ball balancer can be described aswhere ; with representing the coefficient of viscous damping between the balls and the cavity and denoting the rotational speed of the th ball relative to the bowl. The absolute rotational speed of the ball in can be expressed as where is the spin speed of the bowl in . A detailed illustration of (2) can be found in [8].

The generalized forces provided by the suspension system are a little different from those of the washing machine discussed in [711]. In the washing machine, the spring shown in Figure 1(b) is usually installed under the lower spherical joint. If the supporting bed moves lower along the rod, the spring will compress, which makes the suspension system very soft. However, in the centrifuge, the supporting bed is mainly suspended by the rod and the spring is only used for avoiding jumping of the supporting bed from the lower joint. Considering that the stiffness of the rod is much higher than that of the spring, the force acting on the supporting bed by the rod is estimated bywhere and denote the dynamic and initial distances between the upper joint and the bed, respectively; is the pretightening force acting on the bed by the spring; and are the stiffness coefficients of the rod and the spring, respectively; is the frictional force between the bed and the rod; is the translational velocity of the bed relative to the rod.

Considering that the moments acting on the spherical joints are influenced by axial forces, frictional damping is assumed in all joints. The descriptions of the moments in [11] caused by viscous damping are transformed into frictional models. Take the moment acting on the upper joint shown in Figure 1(b) for example; due to frictional damping [12], it can be described aswhere implies that the joint is detached; is the radius of the joint; denotes the frictional coefficient; for simplicity reason, the stick and slip frictional coefficients in all joints are assumed to be the same; is the rotational speed of the rod relative to the ground. The moments acting on other joints can also be transformed using the same way and no more derivations will be carried out here.

3. Discussions

This section explores dynamics of the three-column centrifuge with and without the ball balancer. The design parameters of the centrifuge including their notes, nominal values, and units are listed in Table 1.

3.1. Dynamics of the Three-Column Centrifuge
3.1.1. Modal Analysis

For a better understanding of the suspension characteristics, this section explores natural frequencies of the centrifuge itself.

The vibration model described in (1) is linearized at its static equilibrium position. Let the generalized velocities and accelerations in (1) be equal to zero; that is, and ; the static equilibrium condition can be expressed as The static equilibrium position can be obtained by Newton-Raphson methodwhere represents the static equilibrium position and is the Jacobian matrix. Natural frequencies can be obtained based on the following linearized system:where and are the mass and Jacobian matrices at , respectively.

Generalized eigenvalues and vectors can be solved from (8), and, during the process, nominal parameter values in Table 1 are employed. The natural frequencies are listed in Table 2 which are the square roots of the eigenvalues. As can be seen, the frequencies can be separated into two isolated groups: the first group comprises the first-, second-, and third-order natural frequencies, while the other one is composed of the fourth-, fifth-, and sixth-order ones.

The first-, second-, and third-order mode shapes are displayed in Figures 2(a), 2(b), and 2(c), respectively. As Figures 2(a) and 2(b) show, the three rods swing along - or -direction simultaneously and the centrifuge performs like a single pendulum. In fact, the first two natural frequencies can also be estimated by where is the acceleration of gravity and is the length of the rods at the equilibrium position which is very close to . Assuming that , the first two natural frequencies can be calculated as Hz. As can been seen, the result almost equals that listed in Table 2.

The third-order mode shape is shown in Figure 2(c) which behaves in a way that the core system rotates around -axis. As can be seen, no elongations or contractions of the rods occur in the first three mode shapes, which makes the corresponding natural frequencies so small.

The fourth-, fifth-, and sixth-order mode shapes are displayed in Figures 2(d), 2(e), and 2(f), respectively. It can be seen that the fourth- and fifth-order mode shapes are similar to each other: both of them rotate around an axis parallel to -axis. However, their rotational centers denoted by “” are different. Figure 2(e) displays the sixth-order mode shape; as the figure reveals, it rotates around -axis. As can be seen, elongations or contractions of the rods are involved in the three figures. Considering that the rods are very stiff, the corresponding natural frequencies have to be big.

3.1.2. Dynamics of the Three-Column Centrifuge

This section discusses dynamics of the three-column centrifuge itself. Considering that the rods are very stiff, vibration angles of the centrifuge around - and -axes are quite small; this makes dynamics of the system very similar to that of a planar rotor. Once above the third-order natural frequency, at the steady state, the centrifuge will rotate around the equivalent mass center of the whole system. We now analyze centrifugal forces acting on the system at the steady state. In Figure 3, , , and are the three hanging joints and denotes the pretightening force applied on the joints by the spring shown in Figure 1(b). Assuming that the three joints are separated evenly along circumferential direction, distance between any two of them is identical. The distance between the central axis of the assembly and joint can be derived as , and the distance between and is .

The equivalent mass center of the core assembly including the motor is denoted by and the mass center of the unbalance is represented by . For a general description, the asymmetry of the system caused by the motor is neglected and the equivalent mass center is assumed to be aligned with . After the unbalance is involved, the equivalent mass center of the whole system will be aligned with around which the whole system will rotate at the steady state. The distance between and is the eccentricity of the whole system which can be calculated aswhere and denote mass of the unbalance and that of the core assembly including the motor, respectively; represents the gyration radius of the unbalance. The gravity forces acting on and are and , respectively. And the centrifugal force acting on iswhere denotes the operating speed. The centrifugal force acting on can be expressed as

At the steady state, the centrifugal forces and will counteract each other and an equilibrium between them will occur. Considering that, in a real centrifuge, the height of relative to the suspension plane is usually higher than that of , moments generated by and cannot counteract each other completely, and the difference between them has to be compensated by other moments such as those generated by the gravity. Considering that both and are proportional to , if grows, the difference could be large enough to tip-over the core assembly. And that is the very reason why the centrifuge cannot run very high. We now discuss the critical speed above which the centrifuge will tip-over. The smallest tip-over speed appears in Figure 3 where all the moments relative to are in a critical equilibrium conditionwhere . Submitting (11), (12), , and into (13) leads to If , tip-over phenomenon of the assembly relative to will occur and will detach. This could cause serous vibrations of the system, especially when a small axial damping of the rods, such as the one shown in Figure 1(b), is employed. For a general description of the tip-over speed, the asymmetry of the system caused by the motor was not considered during the above derivations, and the result calculated by (14) may be a little higher than the real one. This can be seen from Figure 3 where, at this very moment, the gravity force acting on the motor makes the tip-over phenomenon more easily occur. Equation (14) can be adjusted by taking the gravity moment generated by the motor into consideration where with and denoting mass of the assembly excluding the motor and that of the motor itself, respectively; denotes the distance from the mass center of the motor to the central axis .

Figure 4 displays variations of the tip-over speed along with increases of , , and when height of the unbalance relative to the suspension plane . As can be seen from Figure 4(a), in the interval , the tip-over speed grows larger while increases. Different types of unbalance are considered; as the figure reveals, a bigger corresponds to a smaller . Figure 4(b) shows variations of along with the increase of when . Different pretightening forces are considered. As can be seen, a bigger can be obtained by a stronger .

We now validate the above results. Simulations are carried out using Runge-Kutta fourth-order method based on (1). During the process, nominal parameter values in Table 1 are used and the operating speed of the centrifuge is assumed to be where is the operating speed at the steady state and is a time constant. In this paper, the value is adopted.

Figure 5 shows displacement of the core assembly under different conditions. As Figure 4(a) reveals, the tip-over speed is 14.18 Hz when and . Time responses corresponding to and are displayed in Figures 5(a) and 5(b), respectively. As can be seen, a stable periodical solution is achieved in Figure 5(a) while an unstable one is obtained in Figure 5(b). Figure 5(c) displays elongations of the second suspension rod when . As the figure reveals, is sometimes smaller than in the interval which implies that the spherical joint is sometimes detached. Serious vibrations can be seen from Figure 5(b) in the interval . For alleviating the situation, a larger pretightening force can be adopted. According to Figure 4(b), a pretightening force equal to 200 N can guarantee a stable solution even after the operating speed increases to . The corresponding time response is shown in Figure 5(d); as can be seen, a stale solution is obtained.

We now discuss influences of the frictional damping in the spherical joints on dynamics of the system. Figure 6(a) shows bifurcation diagram of the system along with the variation of the frictional coefficient when , , , , , and . As can be seen, instability happens when decreases. The time response and power spectrum of displacement when are shown in Figures 6(b) and 6(c), respectively. Two isolated peaks can be identified in Figure 6(c). The one at corresponds to the operating speed of the centrifuge, and the other one at is very close to the first-order (or second-order) natural frequency of the system. This phenomenon is very similar to the instability that happens in a vertical axis washing machine studied in our previous paper [11] which is mainly caused by a Hopf bifurcation due to small damping in the spherical joints.

Simulations reveal that, for a three-column centrifuge, once the operating speed is higher than 20 Hz, in order to suppress the instability, both the frictional coefficient and the axial damping of the rods have to be big. Figure 7 displays displacement and the corresponding power spectrum when , , , , and . As can be seen, even though a large axial damping and a big are applied, the instability can still not be completely eliminated. According to (15), if equals , should be infinite; however, due to the instability born with the suspension structure, the maximum operating speed of the centrifuge itself might be constrained.

3.2. Dynamics of the Three-Column Centrifuge with a Ball Balancer

This section studies dynamics of the three-column centrifuge when a ball balancer is employed. A two-ball model is employed and no collision happening between the balls is assumed. Besides, for exploring fundamental dynamics of the system, only viscous damping is assumed in the balancer.

Instability due to the balancer usually happens in the neighborhood of the critical frequencies [1315] of a rotor. Simulations reveal that instability of the balancer does occur near the first-order (or second-order) natural frequency of the three-column centrifuge; however, the influence can be neglected because the operating speed when instability happens is only about 2 Hz; at this circumstance, vibrations caused by the instability are too small to be considered. Thus, instability characteristics of the balancer will not be presented in this paper.

The study in this section mainly focuses on the tip-over phenomenon of the three-column centrifuge after the balancer is employed. This phenomenon usually occurs at the steady state when the operating speed is much higher than the first-order (or second-order) natural frequency of the system. During the studies, no rolling friction and no collisions between the balls are assumed in the balancer and both of the two balls could reach their optimal positions.

Centrifugal forces acting on the system at the steady state are shown in Figure 8. As Figure 8(b) displays, the dashed circle represents the equivalent mass center of the two balls Ball 1 and Ball 2. Forces acting on are the centrifugal force and the gravity force where is the equivalent gyration radius of . If mass of the balls is small, at the steady state, the resultant centrifugal force is not enough to counteract , and, in this situation, equilibrium of the centrifugal forces is governed by the following equation:where , , and , from which the eccentricity of the system can be derived as

Figure 9 displays variations of , , , , and along with the increase of when and . As can be seen from Figure 9(a), when mass of the balls is small, the two-ball group together at the opposite side of the unbalance and equals . However, once grows bigger than 0.8 kg, gets smaller and smaller. Figure 9(b) reveals changes of the centrifugal force acting on when varies; as can be seen, at first, decreases as increases; however, once is beyond 0.8 kg, almost turns to zero; in this situation, the eccentricity of the system shown in Figure 9(c) almost disappears. Figure 9(d) shows variations of the resultant force . As can be seen, the resultant force is very small, which proves the rationality of (17).

According to Figure 9, if mass of the balls is big enough, vibrations of the system will vanish and the centrifugal force will be zero. In this situation, (17) can be rewritten as from which the critical value of can be derived after submitting into the above equationActually, corresponds to the pitch-fork bifurcation point discussed in a planar rotor [15]. As Figure 9(e) reveals, once , the product stays at a relatively stable value that equals ; this is caused by the idea that the two balls are separated from each other and the gyration radius decreases as increases in order to keep the equilibrium between and .

Based on the above discussions, it can be seen that when mass of the balls is smaller than , the two balls stick together and the equivalent gyration radius of the balls equals ; however, once , the equilibrium between and will appear which leads to . Thus, the equivalent gyration radius of the balls can be summarized as

Although an equilibrium between and can be obtained once , the tip-over phenomenon can still happen. As Figure 8(a) reveals, if the height of the balancer relative to the suspension system does not equal that of the unbalance , a resultant moment relative to will occur, which could lead to a tip-over phenomenon at a high operating speed. Actually, three different types of tip-over phenomena may be induced; the first one is shown in Figure 8(a) which is caused by and governed byThe second type happens when ; as shown in Figure 8(c), this situation is similar to that shown in Figure 3 and governed by A comparison can be made here between (23) and (13), it can be seen that two more items and are counteracting in (23) which is beneficial for obtaining a higher tip-over speed. It should be noted that, due to the pitch-fork bifurcation discussed in (19) and (20), is always smaller than or equal to ; if the height of the balancer is lower than , no matter how big the mass of the balls is, the moment generated by is always smaller than that by and thus a tip-over phenomenon will always happen, even if the eccentricity of the system could reach zero.

We now discuss the third type of tip-over phenomenon that might occur at the transient stage. For doing this, the most serious condition is considered in Figure 8(d). As can be seen, the unbalance and the balls are on the same side; at this circumstance, the equilibrium of the centrifugal moments can be estimated by

We now derive the tip-over speed. For a general description, asymmetry of the system caused by the motor is not considered in the above equations; this can be adjusted by taking the gravity moment generated by the motor into consideration. Different tip-over speeds can be derived based on the above equations and the one governed by (22) can be combined together with that governed by (23):where and the item is employed for adjusting asymmetry of the system due to the motor. The value of in (26) represents the status of the centrifugal moments. A positive implies that , while a negative one indicates that . A zero implies that equilibrium among the centrifugal moments has been achieved, and, at this circumstance, an infinite can be obtained.

The tip-over speed governed by (24) can be derived asConsidering that, at the transient stage, the equivalent gyration radius of the balls is hard to calculate, the most serious condition is assumed in (27).

Figure 10(a) displays variations of and along with the increase of when , , , and . As can be seen, decreases when grows bigger. This implies that, in order to avoid the tip-over phenomenon at the transient stage, a slower starting speed should be considered when a big balancer is employed. As the figure shows, in the interval , increases quickly with ; however, once , decreases rapidly in the interval and then changes slowly once is bigger than .

As Figure 10(a) reveals, when , an optimal mass of the balls exists which could guarantee a very big . Actually, this optimal mass can be obtained by letting (26) be equal to zero; namely, . However, calculating directly from is not easy, because the eccentricity is a function of . This problem can be solved by numerical algorithms such as the bisection method. Figure 10(b) displays variations of the optimal mass of the balls along with the increase of . Considering that , a bigger can be obtained by a closer distance between and , which is helpful for getting a smaller eccentricity of the system.

Figure 10(c) shows variations of when different installing heights of the balancer are considered: , , and . As can be seen, if , no matter how the balancer is designed, the tip-over speed is always limited; if , an optimal mass exists; and if , as long as , complete balance of the centrifugal moments can be obtained which could lead to an infinite .

We now validate the conclusions made in Figure 10. Simulations are carried out using Runge-Kutta fourth-order method. The viscous damping coefficient in the balancer is adopted and the initial positions of the balls are , .

Figures 11(a) and 11(b) display time responses when and . According to Figure 10(a), the corresponding tip-over speed should be . Figures 11(a) and 11(b) display simulation results when and , respectively. As can be seen, the displacement in Figure 11(a) almost turns to zero at the steady state, while, in Figure 11(b), serious vibrations caused by the tip-over of the centrifuge can be observed.

Figure 11(c) displays variations of the positions of the first ball relative to the bowl when , , and . As can be seen, when mass of the balls is much bigger than , although the tip-over phenomenon does not occur, oscillations of the balls appear. This could make the tip-over speed a little smaller than the theoretical one because may be sometimes bigger than during the oscillations of the balls. Thus, a too big balancer is not recommended for the centrifuge.

Figure 11(d) displays the result when , , and . As can be seen, the tip-over speed is much bigger than that in Figure 11(b) where is considered. However, due to the instability caused by a small shown in Figure 7, a big frictional coefficient in the spherical joints has to be applied.

Figure 11(e) shows the simulation result when , , and . For avoiding the tip-over phenomenon during the transient stage, the starting speed of the system is 7 Hz which is smaller than ; then the speed is increased slowly from 7 Hz to 100 Hz with a step of per second; the final speed is reached after 56 seconds; as can be seen, no tip-over phenomenon occurs during the process. This proves that if the installing height of the balancer equals , the operating speed of the three-column centrifuge can be dramatically improved.

3.3. Comparison between Dynamics of the Centrifuge before and after the Balancer Is Installed

For a better understanding of the effect of the balancer, this section makes a simple comparison between dynamics of the centrifuge itself and those after the balancer is installed.

Before the balancer is installed, the eccentricity of the centrifuge equals where denotes mass of the unbalance, is the gyration radius of the unbalance, and is the mass of the core assembly. After the balancer is employed, the eccentricity turns to where is mass of the balls and is the equivalent gyration radius of the balls. As can be seen, if mass of the balls increases, could become smaller and smaller. And once , an equilibrium between and will happen which makes . From this point of view, the balancer could reduce vibrations of the system directly.

We now compare the tip-over speeds. The tip-over phenomenon could still happen even after the eccentricity of the system is completely eliminated by the balancer. This is caused by the idea that the tip-over speed is not governed by equilibrium of the centrifugal forces; it is governed by the equilibrium of their corresponding moments relative to the suspension plane. In this situation, the installing height of the balancer plays a very important role. Taking the situation shown in Figure 10(c) for example, where , , , and are considered, different values of could lead to different results: (1) if is smaller than , the tip-over speed could be enhanced too, but the value of is always limited; (2) if is bigger than , there is an optimal mass of the balls that could make infinite; however, vibration of the system still exists because is always smaller than ; (3) if equals , as long as , that is, as long as the eccentricity of the system is completely eliminated by the balancer, an infinite could be achieved.

Side effects may appear after employing the balancer. The first one appears during the transient stage; as can be seen from Figure 8(d), if the positions of the balls and the unbalance are on the same side, the tip-over phenomenon could be induced more easily. The second one is caused by an inappropriate design of the balancer; as Figure 10(c) shows, if is higher than and a too big mass of the balls is employed, the tip-over speed could even be reduced.

Rolling friction acting on the balls may influence the performance of the balancer too. Although the balancer still works, its effect for counteracting the unbalance could be weakened. A detailed discussion of the imperfect location of the balls based on a planar rotor was made in [13] where countermeasures were also studied. To counteract the rolling friction, multiple balancers or a balancer containing multiple races might be adopted. For a three-column centrifuge, two balancers might be considered; the first one can be installed on the upper rim of the bowl and the second one may be fixed on the bottom rim. Theoretically, the effect of the two balancers is equivalent to that of the one fixed on the middle. However, it might have a better capability of coping with the rolling friction.

4. Conclusions

Modal analysis of a three-column centrifuge was carried out with two isolated natural frequency groups identified. The operating speed of the centrifuge should be higher than the highest frequency in the first group and lower than the smallest one in the second one.

Forces acting on the system at the steady state were studied. Based on the equilibrium of the centrifugal moments relative to the suspension system, the tip-over speed of the centrifuge without balancer was derived and then validated by simulations. The ideal operating speed of the centrifuge could be infinite if height of the unbalance equals that of the core assembly; however, due to instability caused by a small frictional coefficient in the spherical joints, the maximum operating speed of the centrifuge might be constrained.

Influences of the ball balancer on dynamics of the centrifuge were studied. If the installing height of the balancer is higher than that of the unbalance , there is an optimal mass of the balls that could keep the operating speed infinite. If is smaller than , no matter how the balancer is designed, the maximum operating speed of the centrifuge is always limited. The ideal condition is ; as long as the mass of the balls is bigger than or equal to , the operating speed can be infinite.

Rolling friction acting on the balls may have negative effect on the performance of the balancer. Some future works have to be carried out on how to cope with it. A possible way is employing two balancers: one on the top side of the bowl and the other one on the bottom side. It might have a better capability of coping with the rolling friction; however, its performance still needs to be proved by experiments.

Conflict of Interests

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

Acknowledgments

The authors would like to acknowledge the support by China Postdoctoral Science Foundation (2013M531271), National Natural Science Foundation of China (51205166), National Natural Science Foundation of China (51575236), and China Scholarship Council (201308320011).