Research Article  Open Access
Dynamics Analysis on Piezoelectric Laminated Vibrator and Optimization of PZT Position
Abstract
Piezoelectric laminated structure is widely used as actuator’s drive part. The different position of PZT on a piezoelectric vibrator causes different incentive effects. Therefore, seeking an optimal PZT position is of great significance to improve actuator’s drive forces and electromechanical conversion efficiency. In this research, the optimization of PZT position was studied using the approximate solution of piezoelectric vibrator mode shape with mutation sections. The vibration mode function was expressed as a linear superposition of the admissible function according to Rayleigh Ritz method. Then solving of functional variation was converted into the solving of the coefficient matrix of the admissible function by Hamilton’s principle. Through analyzing the forms of admissible functions, the admissible functions that satisfied the boundary conditions of displacement were chosen. For a given vibrator, approximate function for natural frequency and vibration mode was calculated in given admissible functions. Calculated values and experimental results were compared. Results showed that the more items an admissible function sequence had the closer the calculated results were to the experimental values. The errors of calculations were analyzed based on the selection of admissible functions and electromechanical coupling effect. Optimization of PZT position was achieved by analyzing the mode forces of the piezoelectric laminated vibrator.
1. Introduction
Piezoelectric laminated structure is widely used as actuator’s drive part as it is well characterized with simple structure, convenient control, quick response, and no electromagnetic pollution [1]. Researchers from all over the world have been involved in the study on theoretical models of such structures extensively [2–5]. Research on how to paste PZT on traveling wave ultrasonic motor indicated that the electromechanical conversion efficiency of piezoelectric ceramics was closely related to the force coefficient which depended on the size and position of the ceramic transducer as well as the mode shape of the stator [6]. Thus, to improve the electromechanical conversion efficiency of the vibrator theoretically, the first thing that must be done is to obtain the vibrator’s mode shape and then to determine optimal size and location of a piezoelectric ceramics transducer accordingly. Koplow et al. conducted dynamic analysis on the free vibration of discontinuous beam with onestep change in cross section and predicted the dynamic response of the beam [7]. Sadri et al. established a theoretical model for a thin plate covered with piezoelectric sheets and derived control equations for the plate based on RayleighRitz method; based on genetic algorithms the optimal placement of piezoelectric actuators was determined, but followup experiment verification has not been reported [8]. Bashash et al. divided a variable cross section stepped homogeneous material beam with arbitrary boundary into several segments and built models for them. As a result, characteristic matrix was formed by taking boundary conditions and continuity conditions into consideration [9]. For nonhomogeneous piezoelectric beam, however, this approach created a highorder algebraic equation, which is difficult to solve. Crawley and de Luis studied the strain transmitting mechanism from piezoelectric actuator to substructure and conducted researches on PZT optimization position of constant cross section piezoelectric beam and ignored the influence of PZT on its vibration mode [10]. From the perspective of vibration suppression, Ren and Jiang presented a systematic approach for the free vibration analysis and forced response of the beam bonded with PZT employing the travelling wave method [11]. They found that the PZT bonded position near the fixed end in beam had the powerful actuated capability. However, experiments demonstrated that this conclusion was only applicable to vibration suppression. Cho et al. established an equivalent electric circuit for piezoelectric bimorph beam by using piezoelectric PZT as equivalent to an impedance element [12], and the researches provided the new idea for mathematical model of piezoelectric laminated structure. However, the PZT position optimization was not mentioned. Xu and Zhou presented a twodimensional analysis for piezoelectric beam with variable thickness under the boundary condition of the simple support [13]. Mao proposed a solution for the problem of finding the shape of piezoelectric mode sensors for nonuniform EulerBernoulli beams with rectangular cross sections [14]. Cupiał discussed a perturbation solution of the natural frequencies and mode shapes of a piezoelectric rectangular plate [15].
Daraji and Hale applied piezoelectric elements in active vibration reduction and proposed fitness and objective functions to seek optimal distribution of segmented sensors on a flexible plate [16, 17]. These researchers also developed a model for complex isotropic plate using FEM and Hamilton’s principle [18].
Cazzulani et al. used the independent mode space control logics to achieve the optimal actuator and sensor placement [19]. This most popular model control logic used independently controlled modes’ damping and stiffness. In this method, Ambrosio et al. proposed a negative derivative feedback strategy, which feeds velocity back to the mode and also generates a control input proportional to the compensator velocity. The feedback was able to filter out both the higher and lower uncontrolled components. In active vibration control of smart structures, the actuator and sensor placement was a key point of the control system design. Ambrosio et al. [20] also proposed an H2 norm approach to solve actuator and sensor placement. This approach was examined using a finite element model of a square plate. Piezoelectric patch actuators and acceleration sensors were mounted on the three sides of the square plate. To seek for optimal actuator and sensor placement, researchers have realized that more research effort should be made to driving forces and deformation of a piezoelectric vibrator in piezoelectric motor, micro pump, and scanner.
PZT’s excitation to the vibrator is achieved through the substrate deformation caused by the piezoelectric effect. Thus, the matching relations between PZT’s position and mode shape of the vibrator have great influences upon driving forces and electromechanical conversion efficiency. From the point of mode shape, this paper describes the position optimization of a PZT on a piezoelectric laminate vibrator. First, the vibrator is divided into several segments according to the cross section mutation, and then the dynamic analysis for the piezoelectric laminate structure was carried out. Then, approximate solutions of vibration modes were obtained using RayleighRitz method. These approximate solutions were compared to experimental results at their resonant frequencies. Furthermore, force coefficients at different PZT positions were discussed. These results will provide a theoretical foundation for optimizing the design of the piezoelectric laminated beam. In this paper the research objectives were to (1) conduct dynamic analysis of the piezoelectric laminated structure, (2) derive approximate solutions of vibration modes using RayleighRitz method, and (3) propose an optimization method for the design of a piezoelectric laminated beam.
2. Dynamic Analysis on the Shape of Vibration Modes
Assume that coordinates of cross section change positions of the stepped beam of piezoelectric laminate structure in direction are , , and , as shown in Figure 1. The piezoelectric laminate structure was divided into several segments with equal cross section and differential equations were given for them, respectively. A characteristic matrix for the whole structure was assembled based on the continuity relationship between these segments [9].
Differential equation of the segment is where , , and are the mode shape, bending rigidity, and mass of a unit long beam, respectively. Define
Then (1) can be expressed as
The general solution for (3) iswhere , , , and are undetermined coefficients.
For a cantilever beam, deflection and rotation angle on the left end are zero. That is to say,
In this paper, this variable cross section beam is divided into three segments. At the right end of the beam, the boundary condition is free, and both bending moment and shear force are zero. Consider
At the mutation sections of the beam, the deflection, rotation angle, bending moment, and shear force are continuous. That is,
The characteristic equation of the system is obtained by substituting (5)~(10) into (4). Obviously, is a function of the whole beam’s natural frequency . can be expressed aswhere and .
The characteristic equation can also be expressed in the form of a matrix: where is a characteristic matrix and is a feature vector which means a coefficient vector of mode shape. Consider consists of three parts, as shown in (14). The first part and the third part represent the boundary conditions at the positions and , respectively. Part represents the continuous relationship between these segments at the mutation sections. Considerwhere and are order matrices and is an order matrix.
The determinant of is zero since the feature vector has a nonzero solution. Consider
Analytical expression for is obtained by solving (15), and then the coefficients of mode function in (4) can be solved. It is difficult to get a nonzero solution since the equation is a higherorder transcendental equation after being unfolded. In practice, the equation is often approximately solved by making some reasonable simplifications.
3. Approximate Solution of the Mode Shape
The essence for solving mode shape is to obtain the real vibrational track of the beam under a certain order mode shape. According to Hamilton’s principle, the real motion’s variation for Hamilton’s processing object equals zero:where is a variational symbol, is Hamilton’s magnitude, and is a Lagrange function defined as , where is system kinetic energy, is potential energy, and is generalized coordinates. Hamilton’s principle gives a way to find a true movement from all possible movements and it is generally applicable in mechanics. It has become a basic principle in mechanics field.
RayleighRitz method is a direct solution for this kind of variation problems, which extends Hamilton’s principle into elastic dynamics. Its basic concept is to assume the approximate solution of a vibration mode function as follows: where is defined as admissible function, is undermined coefficient, and is the item number of the admissible function sequence.
Equation (17) shows that approximate solution of a vibration mode function can be expressed by a linear combination of function sequences . The variational function can be converted into a function of coefficients . Coefficient should meet the condition that the function gets a stationary value. ConsiderSupposing that the beam’s deflection is , then the kinetic energy of the system is Potential energy of the system can be written as The fonctionelle of Hamilton’s processing object for the beam’s bending vibration can be expressed as where is the mass of a beam in unit length and is flexural rigidity of the beam.
When the piezoelectric beam is doing a certainorder mode vibration, Substitute formula (22) into formula (21):
Equation (23) is the fonctionelle of the function depending on the independent variable and . In one cycle, which is in a range of , (23) can be simplified as a fonctionelle of if, integrating on ,
Substitute (17) into (24); then the fonctionelle can be changed into a function of . Consider
Equation (25) can also be written in a quadratic format about . Consider
Defining , it is called generalized mass; defining , it is called generalized stiffness.
Equation (26) can be written in the form of a quadratic matrix: where
Calculate the derivatives on in formula (27), respectively, and let them be equal to zero. It can be expressed by a matrix as follows:
In order to get the nonzero solutions of , its determinant of feature matrix should be equal to zero. Consider
Unfold the elements of matrix and into a segmented form: where and are mass of the unit length beam and flexural rigidity of the segment of the uniform beam (). And here
Then generalized mass matrix and generalized stiffness matrix can be solved using (31).
The choice for an admissible function , which is shown in (17), has significant influence on calculating accuracy and convergence rate. Ideal admissible function should not only meet the displacement boundary conditions of structural vibration but also meet the force boundary conditions. However, it is very difficult to obtain an admissible function which meets both boundary conditions. In this research, an admissible function can be acceptable as long as it meets the displacement boundary conditions. The larger the number of items is the more accurate the approximate solution becomes. When , the approximate solution is its exact solution.
An admissible function sequence that meets the displacement boundary conditions can be assumed as
When , , and , approximate vibration mode functions are as follows:
4. Comparison of Theoretical and Experimental Results
In this section, the approximation method for the mode used in earlier sections was verified through vibration tests. As the first step, for the given vibrator structure, approximate function for frequency and mode shape was calculated with different numbers of items, . The surface of the vibrator was scanned with a laser Doppler vibrometer to obtain measured value of resonant frequency. The comparison between calculated values and experimental values was executed, and causes of error were analyzed.
The structure of the piezoelectric laminated vibrator in the experiment is shown in Figure 1. The PZT material was P81, and substrate material was phosphor bronze. The properties of the materials and the structural parameters are shown in Tables 1 and 2. The width of the vibrator was 16 mm.


These dimensions indicate that the midpoint of the ceramic sheet was exactly coincident with the midpoint of the substrate. When the vibrator works in higher orders of vibration modes (over the third order), the amplitude was very small. Thus, only the first two modes were considered in the paper. By taking data of Tables 1 and 2 into formulas (31)(32), elements of matrix and were obtained as shown in (30). Solving (30) value was then obtained while brings the information of resonant frequency (). Then, taking the value into (29), coefficients of mode function can be solved as , , and . The calculation results showed that only resonance frequency of the firstorder bending vibration is obtained as Hz when took a combination of admissible function with the first term in (34). When took a combination of admissible function with the first two terms in (35), resonance frequencies of the firstorder and secondorder bending vibrations were obtained as Hz and Hz. When was a combination of admissible function with the first three items in (36), resonance frequencies of the firstorder and secondorder bending vibrations were Hz and Hz as shown in Table 3.

To verify the results from the method described above, a laser based noncontact instrument, Doppler vibrometer, was used to perform the scan test on the vibrator. This laser Doppler vibrometer system’s components and working principle are shown in Figure 2. Operating parameters (e.g., scan range and excitation voltage) are specified on the terminal device. In this experiment, scanning frequency ranges from 10 to 3000 Hz. Excitation voltage was 60 V, and the scanning signal was assigned by the master controller, which motivates the vibrator to work. Scanning laser head (OFV56) sends out frequencystabilized laser to scan the surface of the vibrator and collects the reflected light. Doppler’s frequency shift signal, which is proportional to the target’s velocity, is generated through interference phenomenon. After the signal is processed by a decoded controller, velocity and displacement analog values for the measured object are received. In addition, the resonant frequencies and vibration modes of the vibrator can also be obtained. All of these data will go through a highspeed A/D converter to be processed and displayed on a PC.
The calculated values and experimental results are shown in Table 3.
It can be seen from Table 3 that with just taking the first term of the admissible function into consideration , the error of the theoretical value relative to the experimental value is 46.7% and that with just taking the first two terms of the admissible function sequence to form a basic function (), the error of the theoretical value relative to the experimental value is 28.3% for the firstorder bending vibration. Similarly, it is 39.5% for the secondorder bending vibration. With taking the first three terms to form a function (), the relative error between them is 23.6% for the firstorder bending vibration and 2.1% for the secondorder bending vibration. These results indicate that with more terms the approximate solution becomes more accurate. The calculated values are always greater than experimental values. These errors are mainly due to following factors.
(1) Approximation of the Admissible Function. The admissible functions only meet the displacement boundary conditions, not the force boundary conditions. The beam’s ends are free in its physical model, but there are bending moment and shear force on the right end of the beam when deriving the admissible function, which is equivalent to the imposed additional constraints.
(2) Influence Caused by the External Electric Field. With the action of an electric field, the stressstrain relationship of the PZT beam can be approximately expressed using following equation:where is the stress of PZT in direction; is the elastic modulus of PZT in direction; is strain of PZT in direction; is the piezoelectric constant of PZT; is the voltage of external electric field; and is the thickness of PZT.
Equation (37) can be simplified as follows:
Equations (37) and (38) indicate that the stress of the PZT generated from a certain strain not only is associated with the strain itself but also is related to the piezoelectric effect induced from external electric field. Under the same strain, the stress caused by PZT with external electric field became smaller than that with no external electric field. That is to say that the effective elastic modulus is smaller than that when the external electric field is open. The higher the applied voltage in external electric field is the faster the stress decreases.
The theoretical calculation in this paper was carried out without consideration of electromechanical coupling effect from the PZT. While measuring these values at the driving voltage of 60 V, the effective elastic modulus reduces due to the action of external electric field. Thus, measured mode frequencies were less than the theoretical values.
When , for example, the firstorder and secondorder approximate vibration mode functions can be written as
For the piezoelectric bimorph vibrator described in Figure 1, using the secondorder mode can get higher electromechanical coupling effect [21]. Thus, optimized analysis of later part is aimed at the secondorder mode in this study. After the length of piezoelectric beam is normalized, the curve of the secondorder vibration mode shape drawn using (39) is shown in Figure 3. When comparing this curve with the vibration curve of a uniform beam, it indicates that (1) when the midpoint of the PZT beam is exactly coincident with that of the substrate the secondorder vibration curve’s extreme point of multiplestep beam is almost coincident with that of the uniform beam and that (2) after the vibration amplitudes on the two ends are normalized the deformation of multiplestep beam is smaller than that of uniform beam due to the influence caused by PZT.
5. Optimization of PZT Position
The secondorder vibration mode of the vibrator discussed above was the case when the length of PZT was nearly equal to the length of the substrate and the position of PZT was located at the center of the substrate. Next scenario to be discussed is when the length of is smaller than the half wavelength of the vibration mode curve. In addition, optimal pasting position for a PZT with ideal paste condition will be discussed.
Under the ideal paste condition of the PZT, the thickness of adhesive layer between the PZT and the substrate tends to zero, and then the force exerted on the substrate by the PZT focuses on two endpoints of the ceramic sheet [10]; the relationship between the force and the parameters of PZT is expressed as
At the left end of the ceramic sheet, the force is expressed as
At the right end of the ceramic sheet, the force is expressed as
In (40), (41), and (42), and are concentrated force exerted on the substrate by PZT; is the elastic modulus of substrate; is the width of PZT and is the thickness of the substrate; is the stiffness ratio between the substrate and PZT; is the surface strain of the substrate; is a constant that reflects the strain distribution of the beam, and, for a linear EulerBernoulli beam, ; , where is the piezoelectric constant, is the voltage of external electric field, and is the thickness of PZT.
The action force and moment of PZT on the piezoelectric vibrator are shown in Figure 4, in which is the distance from fixed end to the midpoint of PZT, is the moment exerted from PZT to the vibrator. Generalized force of the piezoelectric laminated vibrator is expressed as follows: where is vibration mode of the vibrator and is the thickness of the substrate.
When the system is vibrating, the strain on top and bottom surfaces of the substrate is
Substituting (44) into (41) and (42), respectively, and then substituting (41) and (42) into (43), respectively, thus generalized vibrating force can be achieved as follows:For the first half of (45), except for , it can be expressed as Then (45) can be simplified aswhere
When the piezoelectric vibrator works near the resonant frequency, it can be treated as a single degree of freedom vibration, and the vibration equation is
Substituting (45) into (49), then (49) became
Equations (50) and (48) indicate that the mode force not only is associated with piezoelectric material properties , voltage of external electric field, and structural parameter but also related to the position of PZT relative to the mode shapes of the vibrator. In this study, the part in brackets of (48) is defined as the positional force coefficient of PZT. Consider
Then (51) can be expressed as
From (51) and (53) it can be seen that, with the material and geometry sizes fixed, the vibration shape curve of the piezoelectric vibrator and the positional force coefficient of PZT will change along with the change of PZT position. Therefore, the position of PZT on the piezoelectric vibrator will have optimal solution ever if the influence of PZT on the piezoelectric vibrator is considered.
According to the conclusion above, the first three admissible functions in sequence of the function shown in (33) were selected to study the variation of the positional force coefficient of PZT.
The structure of the vibrator is shown in Figure 1 and the parameters are given in Table 4. The width of the vibrator is 16 mm.

Through calculation, it indicates that the positional force coefficient of the secondorder bending vibration changes along with the changing of PZT position on the vibrator, as shown in Figure 5. From the changing curve, it can be seen that when mm the value of reaches its minimum and PZT’s incentive effect is lowest and when mm the value of reaches its maximum, and then the PZT’s incentive effect is the best.
6. Conclusions
In this paper the researches on the influence laws of PZT position and mode shape on performance of piezoelectric vibrator are conducted. And a new optimization method for the design of a piezoelectric laminated beam is proposed. The outcomes of this study are summarized as follows.(1)After dividing and assembling the piezoelectric laminated vibrator with mutation sections, the dynamics equation is deduced based on the continuity boundary conditions of the vibrator.(2)The form of vibration mode function was assumed based on the Hamilton principle and RayleighRitz method, and a sequence of admissible functions was chosen which met the displacement boundary conditions.(3)Approximate functions of frequency and mode shape were calculated by giving different values to the item number of the admissible function sequence, and the calculated results were compared with experimental results. It was concluded that with more items the approximate solution became more accurate. The cause of these errors was analyzed based on the selection of admissible function and the external electric field.(4)Mode force of the piezoelectric laminated vibrator was analyzed with the action of PZT’s inverse piezoelectric effect. The result indicates that mode force not only was associated with piezoelectric material properties , voltage of external electric field, and structural parameter but also related to the position of PZT relative to the mode shape curve of the vibrator.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This project was supported by National Natural Science Foundation of China (nos. 51475221 and 51375227), Zhejiang Provincial Natural Science Foundation of China (no. LY13E050015), and Application Type University Construction Project of Qingdao Agricultural University.
References
 C. S. Zhao, Ultrasonic Motor Technology and Application, Chinese Science Press, Beijing, China, 2007.
 R. R. Ramos, J. A. O. Hernandez, J. B. Castillero, and F. J. Sabina, “Electromechanical properties of laminated piezoelectric composites,” Mechanics of Composite Materials, vol. 32, no. 3, pp. 286–291, 1996. View at: Publisher Site  Google Scholar
 P. Dash and B. N. Singh, “Nonlinear free vibration of piezoelectric laminated composite plate,” Finite Elements in Analysis and Design, vol. 45, no. 10, pp. 686–694, 2009. View at: Publisher Site  Google Scholar
 J. A. Otero, R. RodríguezRamos, G. Monsivais, and R. PérezAlvarez, “Dynamical behavior of a layered piezocomposite using the asymptotic homogenization method,” Mechanics of Materials, vol. 37, no. 1, pp. 33–44, 2005. View at: Publisher Site  Google Scholar
 Y. Q. Huang, X. Y. Zhao, and S. X. Chen, “RayleighRitz method in elastic dynamics,” Journal of Hebei Institute of Technology, vol. 18, no. 2, pp. 62–66, 1996. View at: Google Scholar
 H. Zhu, Research on Key Techniques of Traveling Wave Type and Rod Shape Ultrasonic Motor, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2006.
 M. A. Koplow, A. Bhattacharyya, and B. P. Mann, “Closed form solutions for the dynamic response of EulerBernoulli beams with step changes in cross section,” Journal of Sound and Vibration, vol. 295, no. 12, pp. 214–225, 2006. View at: Publisher Site  Google Scholar
 A. M. Sadri, J. R. Wright, and R. J. Wynne, “Modelling and optimal placement of piezoelectric actuators in isotropic plates using genetic algorithms,” Smart Materials and Structures, vol. 8, no. 4, pp. 490–498, 1999. View at: Publisher Site  Google Scholar
 S. Bashash, A. SalehiKhojin, and N. Jalili, “Forced vibration analysis of flexible EulerBernoulli beams with geometrical discontinuities,” in Proceedings of the American Control Conference (ACC '08), pp. 4029–4034, IEEE, Seattle, Wash, USA, June 2008. View at: Publisher Site  Google Scholar
 E. F. Crawley and J. de Luis, “Use of piezoelectric actuators as elements of intelligent structures,” AIAA Journal, vol. 25, no. 10, pp. 1373–1385, 1987. View at: Publisher Site  Google Scholar
 J. T. Ren and J. S. Jiang, “Wave mode characteristics on piezoelectric stepped beam,” Journal of Mechanics, vol. 36, no. 5, pp. 540–548, 2004. View at: Google Scholar
 Y. S. Cho, Y. E. Pak, C. S. Han, and S. K. Ha, “Fiveport equivalent electric circuit of piezoelectric bimorph beam,” Sensors and Actuators A: Physical, vol. 84, no. 1, pp. 140–148, 2000. View at: Publisher Site  Google Scholar
 Y. Xu and D. Zhou, “Twodimensional analysis of simply supported piezoelectric beams with variable thickness,” Applied Mathematical Modelling, vol. 35, no. 9, pp. 4458–4472, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Q. Mao, “Design of piezoelectric mode sensor for nonuniform Euler–Bernoulli beams with rectangular crosssection by using differential transformation method,” Mechanical Systems and Signal Processing, vol. 33, pp. 142–154, 2012. View at: Publisher Site  Google Scholar
 P. Cupiał, “Threedimensional perturbation solution of the natural vibrations of piezoelectric rectangular plates,” Journal of Sound and Vibration, vol. 35, no. 10, pp. 143–160, 2015. View at: Google Scholar
 A. H. Daraji and J. M. Hale, “Reduction of structural weight, costs and complexity of a control system in the active vibration reduction of flexible structures,” Smart Materials and Structures, vol. 23, no. 9, Article ID 095013, 2014. View at: Publisher Site  Google Scholar
 A. H. Daraji and J. M. Hale, “Active vibration reduction of a flexible structure bonded with optimised piezoelectric pairs using half and quarter chromosomes in genetic algorithms,” Journal of Physics: Conference Series, vol. 382, no. 1, Article ID 012039, 2012. View at: Publisher Site  Google Scholar
 A. H. Daraji and J. M. Hale, “Active vibration reduction by optimally placed sensors and actuators with application to stiffened plates by beams,” Smart Materials and Structures, vol. 23, no. 11, Article ID 115018, 2014. View at: Publisher Site  Google Scholar
 G. Cazzulani, F. Resta, F. Ripamonti, and R. Zanzi, “Negative derivative feedback for vibration control of flexible structures,” Smart Materials and Structures, vol. 21, no. 7, Article ID 075024, 2012. View at: Publisher Site  Google Scholar
 P. Ambrosio, F. Resta, and F. Ripamonti, “An ${H}_{2}$ norm approach for the actuator and sensor placement in vibration control of a smart structure,” Smart Materials and Structures, vol. 21, no. 12, Article ID 125016, 2012. View at: Publisher Site  Google Scholar
 X. Q. Hu, J. H. Zhang, Y. Huang et al., “Bionic valveless piezoelectric pump,” Journal of Vibration, Measurement & Diagnosis, supplement 1, pp. 132–135, 2012. View at: Google Scholar
Copyright
Copyright © 2016 Caiqi Hu 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.