Shock and Vibration
Volume 2017 (2017), Article ID 1582547, 13 pages
https://doi.org/10.1155/2017/1582547

Research Article

Dynamic Model and Vibration Characteristics of Planar 3-RRR Parallel Manipulator with Flexible Intermediate Links considering Exact Boundary Conditions

Lianchao Sheng, Wei Li, Yuqiao Wang, Mengbao Fan, and Xuefeng Yang

School of Mechatronic Engineering, China University of Mining and Technology, Xuzhou 221116, China

Correspondence should be addressed to Wei Li; cmeecumt512@yahoo.com

Received 11 January 2017; Revised 28 March 2017; Accepted 24 May 2017; Published 12 July 2017

Academic Editor: Tai Thai

Copyright © 2017 Lianchao Sheng et al. 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.

Abstract.

Due to the complexity of the dynamic model of a planar 3-RRR flexible parallel manipulator (FPM), it is often difficult to achieve active vibration control algorithm based on the system dynamic model. To establish a simple and efficient dynamic model of the planar 3-RRR FPM to study its dynamic characteristics and build a controller conveniently, firstly, considering the effect of rigid-flexible coupling and the moment of inertia at the end of the flexible intermediate link, the modal function is determined with the pinned-free boundary condition. Then, considering the main vibration modes of the system, a high-efficiency coupling dynamic model is established on the basis of guaranteeing the model control accuracy. According to the model, the modal characteristics of the flexible intermediate link are analyzed and compared with the modal test results. The results show that the model can effectively reflect the main vibration modes of the planar 3-RRR FPM; in addition the model can be used to analyze the effects of inertial and coupling forces on the dynamics model and the drive torque of the drive motor. Because this model is of the less dynamic parameters, it is convenient to carry out the control program.

1. Introduction

Compared with the serial manipulator, the parallel manipulator has the advantages of high speed, high acceleration, high load capacity, low energy consumption, no-accumulation error, and high precision, which is widely used in aerospace, precision optical instruments, high-precision real-time measuring instruments, fast precision machine tools, and other fields [1]. Simultaneously, FPMs have attracted great social concerns and have been applied extensively. However, owning to the inertia force, dynamic load carrying capacity, and other factors, flexible intermediate links are prone to elastic deformation, which results in the decrease of the overall dynamic performance of the manipulator and influences the normal working of the FPM under the high speed with load condition. Elastic vibration has a greater impact on the accuracy and quality of the operation in the system [2, 3], for example, the elastic vibration of pick-and-place parallel manipulator in food packaging, directly affecting the quality of food [4]; the vibration of operating mechanism at the end of picking manipulator, leading to instability in grasping and damage to strawberry, tomatoes, and other soft appearance fruits [5, 6]. Therefore, the study of vibration characteristics of the FPM has important significance.

In order to further study the vibration problem of the planar 3-RRR FPM with flexible intermediate links, the dynamic model needs to be established which can accurately reflect the rigid-flexible coupling characteristics of the system. Furthermore, the vibration characteristics and the influences should be studied between the rigid body and the flexible body based on the model. In the existing literatures, the dynamic modeling method of the FPM is mainly from the energy point of view to obtain the partial differential equation. Then an approximation finite-dimensional ordinary differential equation is obtained by discretizing the complex continuous system model through the finite element method (FEM) [7, 8] and the AMM [9, 10]. Yu and Hong [11] described dynamic modeling, modal selection and modal synthesis, and dynamic stiffening of flexible multibody systems. Yu et al. [12] presented a dynamic model of the 3-RRR FPM and verified the correctness of the model through experiments. Zhang et al. [13, 14] established the dynamic model of the 3-RRR FPM using the finite element method and Lagrange equation and analyzed the influence of temperature and other factors on the elastic deformation of each flexible link. Zhang et al. [15] established the dynamic model of the 3-RRR FPM based on the energy method and illustrated the influence of the joint gap on the model. Fattah et al. [16] studied the dynamics of a 3-DOF parallel manipulator. Based on the finite element method and the Euler-Lagrange equation, the dynamic equations of the system were established. In order to emphasize the influence of system flexibility, the dynamics simulation experiment of rigid-flexible system was carried out. In spite of the high precision of the above models, the dynamic model described in the past literature has too many dynamic parameters, which makes the model difficult to solve and is disadvantageous to the design of subsequent controllers. In order to meet the requirements of the model premise, as far as possible to make the model simplification, low-level, easy to solve, and controller design, Mostafavi Yazdi and Irani [17] described the basic principle and modeling process of the AMM in detail and investigated the dynamic equations of the beam model with two teeth which were established by the AMM. Compared with the general beam model, the elastic vibration of the tooth beam is smaller. Chen [18] established a linear dynamic model of a multiflexible linkage manipulator by Lagrange equation and AMM. During the dynamic modeling of planar FPM with flexible intermediate links, Zhang et al. [19, 20] elaborated the dynamic model of FPM with flexible intermediate links which was established using the boundary conditions of pinned-pinned at both ends of the flexible intermediate links based on AMM and the Lagrange equation. However, the experiment in [21, 22] proved that, in the motion control, because of the elastic deformation of the flexible link, the dynamic platform produced violent vibration. It can not explain this phenomenon well using pinned-pinned boundary conditions vibration mode to discretize the flexible deformation.

Different from current literatures, in this paper, based on the extended Hamilton principle and AMM, with the pinned-free boundary condition for the flexible intermediate link, the dynamic model of the 3-RRR FPM with flexible intermediate links is established with considering the effect of rotor and the concentrated moment of inertia at both ends of link and the rigid-flexible coupling. The dynamic model established can not only satisfy the model precision but also solve the dynamic model easily. The dynamic response and the natural frequency of the flexible intermediate link are obtained by solving the dynamic equation. Based on the model, the influence of the coupling force and inertial force on the system and the driving torque of the motor are analyzed. Finally, the correctness of the model is verified by modal experiment. This provides a theoretical basis for the optimal design of the system structure and shock absorbers and the selection of motors.

The structure of the paper is as follows. The second chapter introduces the system overview and dynamic modeling of the planar 3-RRR FPM. The numerical simulations and results discussions are shown in third chapter. The forth chapter is the experimental verification section. Finally, the paper is concluded with a brief summary.

2. The System Overview and Dynamic Modeling of the Planar 3-RRR FPM

The system structure and the coordinate system of the planar 3-RRR FPM are shown in Figure 1. is the global coordinate system. is the moving coordinate system which follows the moving platform. The system is composed of the moving platform and the fixed base as well as three identical branches , , and . As Figure 1 shows, , , and are the three vertices of the regular triangle whose side length is 500 mm, where the active revolute joints are set. The points of , , , , , and are the passive revolute joints, where  mm. The radius of the circle of joints in moving platform is  mm. The 3-DOF planar motion can be achieved by using the three driving motors in the three active joints.

Figure 1: The structure diagram of the planar 3-RRR FPM.

In the actual structure, the cross-sectional area of the drive rod is larger to ensure its rigidity, while the cross-sectional area of the intermediate link is smaller to reduce its mass. Therefore, in this paper, only the flexible intermediate links of the flexible planar 3-RRR parallel robots are considered as the flexible part; others are the rigid part. The vibration of the mechanism is mainly caused by the motor drive, the inertial force, and the payload of the moving platform.

2.1. Dynamic Modeling of Flexible Intermediate Link

The schematic diagram of the th branch and moving platform are shown in Figure 2. The transverse deformation is smaller compared with the length of the flexible intermediate link, so the equation can be modeled linearly according to the linear elastic deformation. The flexible intermediate link is considered as a pinned-free Euler-Bernoulli beam. The radial and lateral stiffness is ignored, and the uncertainty of modeling and other parameters is regarded as a kind of parameter perturbation. Since the flexible intermediate link moves in the horizontal plane, the effect of gravity is not considered. The kinetic energy mainly includes the rotational kinetic energy of the motor rotor and the rotational kinetic energy of the joint , the translational kinetic energy of driving rod, the concentrated translational kinetic energy of the joint and its own rotational kinetic energy, and the translational kinetic energy of the intermediate link, where the rotational kinetic energy of joints and is mainly produced by its internal bearing. It can be expressed as

where is the rotational kinetic energy of the th drive joint and the motor rotor; is the kinetic energy of the th driving rod; is the kinetic energy of the joint ; is the kinetic energy of the th flexible intermediate link; is the transverse elastic displacement of points on flexible intermediate link; is the unit length mass of the th flexible intermediate link; is the mass of the joint ; is the rotational inertia of the driving joint and the motor rotor; is the rotational inertia of joint ; is the vector coordinates of the th drive rod; + is the vector coordinates of the th flexible intermediate link; and are the differential for the time and displacement, respectively.

Figure 2: Schematic diagram of the th branch and moving platform.

The potential energy of the th branch of the system is expressed as

where is the elastic modulus and is the moment of inertia of section of flexible intermediate link.

The virtual work of the th branch of the system can be expressed as

where is the constraint force between the joint and the end of the flexible intermediate link. is the drive torque of the th servomotor. is the virtual displacement vector at the end of the th flexible intermediate link. is the virtual angle displacement of the th servomotor.

Based on the extend Hamilton principle [23], the boundary conditions and the kinetic differential equation can be derived as follows:

where is the generalized variable of the th branch.

According to the extended Hamilton principle, the total energy of the th branch is varied to the generalized variable. Considering the interchangeability between variational and differential, the similar items are combined after integrating the variational components by fractional integration over time to . Since can be assigned any value, three differential equations for generalized variables and four boundary conditions for flexible intermediate links are obtained. Since the flexible intermediate link is pinned to joint , the deformation of the pinned joint is zero according to [24, 25]. This conclusion is used to simplify the differential equation.

The differential equation about is as follows:

The differential equation about is as follows:

The differential equation about is as follows:

The four nonhomogeneous boundary conditions of the flexible intermediate link are as follows:

In order to obtain the homogeneous boundary condition for flexible deformation discrete based on AMM, the Dirac function deformation is introduced according to [22].

Substituting (9) into (7) and (8),

As can be seen from (11), the deflection and bending moment at the joint are zero, and the bending moments and shear forces at the joint are zero. Therefore, the pinned-free boundary conditions are satisfied for the flexible intermediate link.

2.2. Elastic Deformation Discretization of Flexible Intermediate Link

After the above analysis, we can find that the established equation is a highly nonlinear differential equation with rigid motion and elastic motion coupling, and the analytical solution is difficult. Therefore, it is considered to use the AMM to discretize the elastic deformation of the flexible link and then to solve the numerical solution for dynamic analysis. According to AMM, the deformation of the flexible intermediate link can be expressed as follows:

where denotes the unknown generalized elastic variable of the th flexible intermediate link. is the mode function corresponding to the known boundary condition.

From the analysis of the previous section, the mode shape function with pinned-free boundary conditions of the flexible intermediate link is used to discrete elastic deformation. As can be seen from [24], the modal functions under the pinned-free boundary condition are shown as follows:

Equations (12) and (13) are substituted into (5), (6), and (7), respectively.

2.3. Dynamic Model of Moving Platform

The moving platform is considered to be a rigid body in the modeling process. The dynamic equation of the moving platform is as follows:

where and is the quality of the moving platform. is the moment of inertia of the moving platform. is the coordinate of the joint in the local coordinate system.

2.4. Rigid-Flexible Coupling Dynamics Model

Combining the dynamic equation of the moving platform with the kinetic equation of the branch, the rigid-flexible coupling ordinary differential equation of the 3-RRR FPM can be obtained as follows:

where is the generalized coordinate variable. , , is the rigid body motion coordinate. is the elastic coordinate of three flexible intermediate links. represents the driving torque of the three motors. is the binding matrix. The detailed form of the matrix is shown in Appendix.

is the constraints of internal force of the joint . is a positive definite symmetric mass matrix. is the centrifugal force and Coriolis matrix. is the stiffness matrix.

3. Numerical Simulation

The specific material and dimension parameters of the planar 3-RRR FPM are as follows: the material is aluminum alloy 7075. The elastic modulus of the material is . The density is . The length, width and thick of the drive rod are 210 mm, 25 mm, and 10 mm, respectively. The length, width, and thick of the flexible intermediate link are 210 mm, 25 mm, and 2 mm, respectively. The moment of inertia of the motor rotor and drive joint is the . The concentrated moment of inertia at the end of flexible intermediate link is . The moment of inertia of section of flexible intermediate link is .

In order to simulate the end of the operating conditions, a typical circular equation is given for its trajectory as follows:

The numerical simulation results are as follows.

Figure 3 shows that the amplitude of the first-order end vibration response of the flexible intermediate link is larger. This vibration is sufficient to affect the stability and dynamic performance of the system. Therefore, with the flexible intermediate link 1 as an example, the second-order vibration response and frequency characteristics are further studied. The amplitude-frequency characteristics can be obtained through fast Fourier transform. The results are shown in Figures 4 and 5.

Figure 3: Vibration response curves of the end of the flexible intermediate links.

Figure 4: The first-order response and amplitude-frequency characteristics of the first flexible intermediate link.

Figure 5: The second-order response and amplitude-frequency characteristics of the first flexible intermediate link.

As shown in Figure 4, the fast Fourier transform of the first-order response curve shows that the frequency components are mainly 24.88 Hz, 49.75 Hz, and 104.5 Hz. Through the same treatment, the second-order frequency components are mainly 24.88 Hz, 49.75 Hz, and 288.6 Hz as shown in Figure 5. In order to further study the composition of frequency components, the inertial force and coupling force analysis are shown in Figure 6.

Figure 6: The response and amplitude-frequency characteristics of inertial force and coupling force.

The mainly frequency of 24.88 Hz and 49.75 Hz caused by inertia force and coupling force can be seen from Figure 6, which also represents that the frequency 104.5 Hz in Figure 4 is natural frequency and the frequency 104.5 Hz in Figure 5 is second-order frequency of the flexible intermediate link. The driving torque of the motor can also be obtained through the dynamic equation, which provides the basis for the motor selection.

It can be seen from Figure 7 that the maximum driving torque of the motor does not exceed 20 N·m.

Figure 7: Three motor drive torque.

4. Modal Test Verification

In order to verify the correctness of the dynamic model established in this paper, the modal test bench is built as shown in Figure 8. The test bench is mainly composed of computer, DH5927N dynamic tester, hammer, piezoelectric acceleration sensor, charge adapter, and flexible 3-RRR FPM. A total of 11 test points were set up, among which the measuring points 1 and 11 were at the bearing, so no measurement was made. During the course of the experiment, the three drive motors are locked in the initial position and the end moving platform is in the central position. The multipoint excitation and the pick-up one point test method is implemented. In order to achieve a more accurate result, the test is performed 5 times repeatedly, and then the results are averaged. The test data is collected and recorded by the DH5927N dynamic tester. The dynamic response is analyzed by PolyLscf method, and the results are shown in Figure 9.

Figure 8: Modal test bench.

Figure 9: Modal test results.

The results show that the first-order frequency is mainly around 103.1 Hz and the damping ratio is 0.056. In addition the second-order vibration frequency is around 286.3 Hz and the damping ratio is 0.032. Compared with the frequency characteristics obtained through modal experiment and theoretical calculation, the first-order frequency 103.1 Hz and the second-order frequency 286.3 Hz obtained through modal test are basically consistent with the theoretical results 104.5 Hz and 288.6 Hz, which proves the correctness of the dynamic model established using the extended Hamilton principle and the AMM with the pinned-free boundary condition. Compared with the results obtained by other researchers [20], the model established in this paper has higher accuracy. The results’ comparison is shown in Table 1.

Table 1: The results comparison.

Project
The results of this article
The results of others

First-order simulation frequency ()
103.1 Hz
70.5 Hz

First-order experiment frequency ()
104.5 Hz
76.6 Hz

Second-order simulation frequency ()
286.3 Hz
280.8 Hz

Second-order experiment frequency ()
288.6 Hz
231.2 Hz

First-order relative error %
1.36%
8.65%

Second-order relative error %
0.8%
17.66%

The source of error between present work and others in the results is mainly in the selection of boundary conditions. This paper chooses the pinned-free boundary condition for the flexible intermediate link in the establishment of the mathematical model, while the boundary condition in the literature [20] is the pinned-pinned, and the joint gap caused by the machining and installation accuracy in the actual test bench is also part of the source of the error.

5. Conclusions

In this paper, based on the extended Hamilton principle and the AMM, the coupled dynamic model of planar 3-RRR FPM with high precision and less dynamic parameters is established with pinned-free boundary condition. This model takes into account the concentrated moment of inertia and rigid-flexible coupling effect on system dynamic. The first-order frequency 103.1 Hz and the second-order frequency 286.3 Hz obtained through modal experiment are very close to the theoretical results 104.5 Hz and 288.6 Hz. The research indicates that the established dynamic model has a high accuracy. Based on the established dynamic model, it is found that the inertial force and the coupling force can cause the forced vibration, the vibration frequencies are 24.88 Hz and 49.75 Hz, and the driving torque of the motor does not exceed 20 N·m which can provide a basis for motor selection. The rigid-flexible coupling dynamics modeling method proposed in this paper is applicable to all kinds of FPM with flexible link. The model has high precision and high efficiency, which solves the problem of complexity of dynamics model and is of great significance to the controller design based on the dynamic model.

Appendix

Nomenclature

: | Active revolute joint
, : | Passive revolute joint
: | Rod length
: | Radius of the circle of joints in moving platform
: | Angle between the driving link and the -axis
: | Angle between the flexible intermediate link and the -axis
, : | Distance between a point and joint on driving link in the , direction
, : | Distance between a point and joint on flexible intermediate link in the , direction
: | Total kinetic energy of the th branch
: | Rotational kinetic energy of the th drive joint and the motor rotor
: | Kinetic energy of the th driving rod
: | Kinetic energy of the joint
: | Kinetic energy of the th flexible intermediate link
: | Transverse elastic displacement of points on flexible intermediate link
: | Assumed mode
: | Material density
: | Unit length mass of the th flexible intermediate link
: | Mass of the joint
: | Rotational inertia of the driving joint and the motor rotor
: | Rotational inertia of joint
: | Vector coordinates of the th drive rod
: | Vector coordinates of the th flexible intermediate link
(): | Differential for the time and displacement, respectively
: | Potential energy of the th branch
: | Elastic modulus
: | Moment of inertia of section of flexible intermediate link
: | Virtual work of the th branch
: | Constraint force between the joint and the end of the flexible intermediate link
: | Drive torque of the th servomotor
: | Virtual angle displacement
: | Generalized variable of the th branch
: | Unknown generalized elastic variable of the th flexible intermediate link
: | Mode function corresponding to the known boundary condition
: | Quality of the moving platform
: | Moment of inertia of the moving platform
: | Coordinate of the joint in the local coordinate system
: | Positive definite symmetric mass matrix
: | Centrifugal force and Coriolis matrix
: | Stiffness matrix
: | Motion coordinate of moving platform
: | Generalized coordinate variable
: | Rigid body motion coordinate
: | Elastic coordinate of three flexible intermediate links
: | Driving torque of the three motors
: | Binding matrix
: | Constraints of internal force of the joint
: | The branch, respectively.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research work was supported by the National Natural Science Foundation of China (no. 51305444, no. 51307172, and no. U1610111), the Scientific and Technological Project of Jiangsu Province (BY201402806), the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), and Six Talent Peaks Project in Jiangsu Province (ZBZZ-041).

References

B. Li, “Dynamic analysis of a 3-DOF parallel manipulator,” Journal of Northeastern University (Natural Science), vol. 31, no. 11, pp. 1608–1610, 2010.

S. K. Dwivedy and P. Eberhard, “Dynamic analysis of flexible manipulators, a literature review,” Mechanism and Machine Theory, vol. 41, no. 7, pp. 749–777, 2006.

J. Chen, P. C. Hagan, and S. Saydam, “Parametric study on the axial performance of a fully grouted cable bolt with a new pull-out test,” International Journal of Mining Science and Technology, vol. 26, no. 1, pp. 53–58, 2016.

K. M. Zheng and Q. J. Zhang, “Elastic dynamics and analysis of vibration characteristics of delta robot with joint clearance,” Transactions of the Chinese Society of Agricultural Engineering, vol. 31, no. 14, pp. 39–48, 2015.

J. Song, T. Zhang, L. Xu, and X. Tang, “Research actuality and prospect of picking robot for fruits and vegetables,” Transactions of the Chinese Society of Agricultural Machinery, vol. 37, no. 5, pp. 158–162, 2006.

Q. C. Li, T. Hu, and C. Y. Wu, “Review of end-effectors in fruit and vegetable harvesting robot,” Transactions of the Chinese Society for Agricultural Machinery, vol. 39, no. 3, pp. 175–179, 2008.

S. Arabzadeh Jarkani, H. Khoshdast, E. Shariat, and A. Sam, “Modeling the effects of mechanical parameters on the hydrodynamic behavior of vertical current classifiers,” International Journal of Mining Science and Technology, vol. 24, no. 1, pp. 123–127, 2014.

X. Wang and J. K. Mills, “Dynamic modeling of a flexible-link planar parallel platform using a substructuring approach,” Mechanism and Machine Theory, vol. 41, no. 6, pp. 671–687, 2006.

M. Loudini, D. Boukhetala, and M. Tadjine, “Comprehensive mathematical modelling of a lightweight flexible link robot manipulator,” International Journal of Modelling, Identification and Control, vol. 2, no. 4, pp. 313–321, 2007.

L. Meirovitch and R. Parker, “Fundamentals of vibrations,” Applied Mechanics Reviews, vol. 54, no. 6, article 100, 2001.

Q. Yu and J. Z. Hong, “Some topics on flexible multi-body system dynamics,” Advances in Mechanics, vol. 29, no. 2, pp. 145–154, 2000.

Y. Yu, Z. Du, J. Yang, and Y. Li, “An experimental study on the dynamics of a 3-RRR flexible parallel robot,” IEEE Transactions on Robotics, vol. 27, no. 5, pp. 992–997, 2011.

Q. H. Zhang and X. M. Zhang, “Dynamic analysis of planar 3-RRR flexible parallel robots under uniform temperature change,” Journal of Vibration and Control, vol. 21, no. 1, pp. 81–104, 2015.

Q. H. Zhang, X. R. Fang, and X. M. Zhang, “Dyanmic analysis of planar 3-RRR flexible parallel robots with dynamic stiffening,” Shock and Vibration, vol. 2014, Article ID 370145, 13 pages, 2014.

X. C. Zhang, X. M. Zhang, and Z. Chen, “Dynamic analysis of a 3- R RR parallel mechanism with multiple clearance joints,” Mechanism & Machine Theory, vol. 78, no. 78, pp. 105–115, 2014.

A. Fattah, J. Angeles, and A. K. Misra, “Dynamics of a 3-DOF spatial parallel manipulator with flexible links,” in Proceedings of the 1995 IEEE International Conference on Robotics and Automation, vol. 1, pp. 627–632, May 1995.

S. J. A. Mostafavi Yazdi and S. Irani, “Transverse vibration of double cracked beam using assumed mode method,” in Proceedings of the 4th International Conference on Recent Advances in Space Technologies 2009, RAST '09, pp. 156–160, June 2009.

W. Chen, “Dynamic modeling of multi-link flexible robotic manipulators,” Computers & Structures, vol. 79, no. 2, pp. 183–195, 2001.

X. Zhang, J. K. Mills, and W. L. Cleghorn, “Dynamic modeling and experimental validation of a 3-PRR parallel manipulator with flexible intermediate links,” Journal of Intelligent and Robotic Systems, vol. 50, no. 4, pp. 323–340, 2007.

X. Y. Wang, “Experimental modal analysis of flexible linkages in a smart parallel platform,” in Proceeding of the 7th Cansmart Meeting - International Workshop on Smart Materials and Structures, vol. 1, pp. 37–46, Montreal, Canada, 2004.

B. Kang and J. K. Mills, “Dynamic modeling of structurally-flexible planar parallel manipulator,” Robotica, vol. 20, no. 3, pp. 329–339, 2002.

X. Zhang, J. K. Mills, and W. L. Cleghorn, “Multi-mode vibration control and position error analysis of parallel manipulator with multiple flexible links,” Transactions of the Canadian Society for Mechanical Engineering, vol. 34, no. 2, pp. 197–213, 2010.

L. Meirovitch, Fundamentals of Vibrations, Waveland Press, 2010.

H. Y. Hu, Mechanical Vibration Foundation, Beijing University of Aeronautics and Astronautics Press, Beijing, China, 2005.

S. Rao Singiresu, Mechanical Vibration, Tsinghua University Press, Beijing, China, 4th edition, Compiled by X. Y. Li and M. L. Zhang, 2009.