- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Advances in Mechanical Engineering

Volume 2013 (2013), Article ID 249121, 10 pages

http://dx.doi.org/10.1155/2013/249121

## Optimum Design of 3-3 Stewart Platform Considering Inertia Property

The State Key Laboratory of Tribology, Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China

Received 16 April 2013; Accepted 17 October 2013

Academic Editor: Jiansheng Li

Copyright © 2013 Zhu-Feng Shao 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

Optimum design is a pivotal approach to fulfill the potential advantages of the parallel manipulator in practical applications. This paper concerns the optimum design issue of the 3-3 Stewart platform considering the inertia property, in addition to the kinematic performance. On the basis of spherical usable workspace, global conditioning index (GCI) is analyzed. Atlases of the workspace and GCI are deduced with the established nondimensional design space. Further, after dynamic modeling, the global inertia index (GII) is deduced from the joint-space inertia matrix, and corresponding GII atlases are drawn. In particular, an example is presented to illustrate the process of obtaining the practical optimum results based on these non-dimensional atlases. Since both kinematic and dynamic properties are considered, the optimum result will possess comprehensive performance improvements.

#### 1. Introduction

The parallel manipulator usually consists of the base and the end effector, which are connected together with at least two identical kinematic chains. Its significant structure feature is the closed-loop kinematic chain. Compared with conventional serial manipulators, parallel manipulators possess some inherent potential advantages [1], such as high accuracy, improved stiffness, increased dynamics, compact structure with multiple degree of freedoms (DOF), and so on. As a result, parallel manipulators attract great academic attentions and are more and more widely adopted in industrial applications, namely, flight simulators [2], telescopes [3–5], and machine tools [6].

Optimum design of the parallel manipulator is recognized as an important theoretical issue, as well as a great practical approach to make full use of its potential advantages in practical applications [7, 8]. Optimization design of the parallel manipulator is always a challenging issue, which is mainly based on the kinematics now [9]. And several kinematic performance criteria could be involved, such as workspace [10], singularity [11], accuracy [12], stiffness [13], and conditioning index [14]. Among them, workspace and conditioning index are the two most important key criteria. On the other side, the dynamic performance of the parallel manipulator is analyzed and evaluated in a lot of works [15, 16]. However, little effort has been made on the dynamic optimum design.

Inertia property is an important part of the manipulator dynamics, which shows its impacts in many ways [17, 18], such as mechanical resonance frequency, acceleration capability, and dynamic response. Inertia property has been invested on both serial manipulators and parallel manipulators [19, 20]. The inertia index of parallel manipulator is proposed in [21], which gives accurate and concise description of the inertia property of the parallel manipulator. In this paper, in addition to kinematic criteria, the inertia index is adopted to carry out the optimum design of the parallel manipulator, which helps us to obtain the optimum results considering both dynamic and kinematic performances.

The study object of this paper is a 3-3 Stewart parallel manipulator, which is also called the octahedral platform [22]. Generally, the Stewart platform is a 6 DOF spatial parallel mechanism [23], which is composed of two bodies (the end effector and the base) and six extensible limbs. Each limb, driven by a set of servomotor and ball screw, is connected with the base by a universal joint at one end and the end effector by a spherical joint at the other end. The 3-3 Stewart manipulator contains three concentric spherical joints and three concentric universal joints. This special construction usually possesses better kinematic performance and larger workspace, which also makes closed-form direct kinematic solution feasible [24].

The remainder of this paper is organized as follows. In the next section, kinematics is analyzed, based on the established coordinate systems. Workspace and conditioning index atlases are obtained in Section 3. In order to obtain the inertia index, the dynamic model is deduced briefly in Section 4. The global inertia index is defined, and corresponding performance atlases are given in Section 5. In Section 6, atlases base optimum design of the parallel manipulator is carried out, considering kinematics as well as dynamics. Finally, conclusions are listed in Section 7.

#### 2. Kinematics and Jacobian Matrix

As shown in Figure 1, the base of the 3-3 Stewart platform is represented by three pairs of concentric universal joints, located at points , , and . And the end effector is defined by three pairs of concentric spherical joints, marked with points , , and . To facilitate analysis, a fixed Cartesian coordinate system is attached to the base as the base frame, whose origin is located at the geometric center, with the -axis pointing at the point and the -axis perpendicular to the base upward. The end-effector coordinate system is located at the geometric center of the end effector, with -axis perpendicular to the end effector upward and the -axis pointing at the midpoint of section . When the end effector is at the central configuration, where the end effector is not rotated relative to the base and limbs are of the same length, the end-effector frame is parallel to the base frame while the -axis and the -axis are collinear.

For convenience, important notations adopted in this paper are shown in Table 1. Generally, vectors and matrices are expressed in bold, while numeric variables are described in italic. Left superscripts of symbols indicate the coordinates systems where the physical quantities are described.

Position of the end effector can be described by the translation vector . The Euler representation is adopted to describe the end-effector posture, which is composed of three successive relative rotations, namely, first rotation of angle about the -axis, followed by second rotation of angle around the -axis, and finally followed by third rotation of angle around the -axis. With respect to the fixed global base coordinate system , the general rotation matrix can be expressed by a multiplication of three basic rotation matrices as follows: with its developed expression where stands for the sine operation and stands for the cosine operation.

According to kinematic schema of the limb shown in Figure 2, the closed-loop position vector equations can be written with rotation matrix and translation vector as where is the limb vector pointing from to , is the position vector of point in base frame, and is the position vector of point in the end-effector frame.

Taking the time derivation of (3), the velocity mapping function can be deduced in the base frame as where the skew-symmetric matrix operator is associated with the angular velocity of the end effector expressed in the fixed base frame . Without end-effector rotation, .

At the same time, can be expressed by the angular velocity and the sliding velocity of the th limb as where is the limb length and is the unit limb vector. The angular velocity could also be expressed as function of two orthogonal angular velocities in the university joint .

Considering (4) and (6) together gives the relation between the end-effector velocity and the limb velocity as Further, using the projections of six vector relations of (7) along the direction of six legs , after organization, we obtain the following matrix form: where and the Jacobian matrix can be written as

#### 3. Kinematic Performance Atlases

##### 3.1. Design Space

Considering the applicability in engineering, the generalized symmetric structure is adopted to carry out the optimum design of the 3-3 Stewart platform. In other words, both spherical joints and universal joints are circumferentially uniformly distributed. Accordingly, the architecture of the object 3-3 Stewart manipulator can be described by four geometric parameters. As shown in Figure 1(b), two parameters are distribution radii of rotational joints and . The initial height between the base and the end effector can be used as the third parameter, when the end effector is at the central configuration with half-stretched limbs. The fourth parameter is the limb length ratio , the proportion of the maximum length, and the minimum length, which can be obtained with practical engineering experiences. In this paper, . Then, the architecture of the object 3-3 Stewart platform can be totally defined by three length parameters, namely, , , and .

In order to embody all possible combinations of these three parameters in a finite area, parameters need to be normalized and physical dimensions must be eliminated. Define a dimension factor and Then, three dimensionless parameters () can be obtained as

To facilitate further analysis, parameters are illustrated in two-dimensional plane, as shown in Figure 3. -axis is parallel with the -axis, -axis and -axis are perpendicular to each other, and the angle between -axis and -axis is . Theoretically, these three dimensionless parameters can take any value from 0 to 3. However, for the Stewart platform, the distribution radius of the end effector should be no larger than the base distribution radius . Thus, the available design space is under the line of and marked with red lines in Figure 3.

##### 3.2. Workspace Atlas

In order to carry out the workspace analysis, let us adopt the Spherical Usable Workspace (SUW), which can be defined as the maximum continuous spherical workspace without singularity locus. Usually, the SUW of the Stewart manipulator is the maximum inscribed sphere of the reachable workspace, encircled with twelve envelope surfaces.

For the Stewart platform, there is only the second type singularity, which can be further classified into three categories, namely, architecture singularity, configuration singularity, and formulation singularity [25]. Without huge end-effector rotation (), only the configuration singularity when the base and the end effector coincide may appear for our object Stewart platform.

The radius is used to represent the size of the SUW and to draw the workspace atlas. Without considering the end-effector rotation, the SUW atlas is deduced and shown in Figure 4. We can find that(1)in the left region of the design space, radius of the SUW is generally proportional to . In this area, the SUW radius is mainly limited by the singularity condition. In other words, the SUW cannot intersect with the plane of ;(2)in the middle part, radius grows with the shrink of the approximately;(3)in the right area, radius curves are almost parallel to the -axis. When is specified, the SUW radius increases with the enlarging difference between and ;(4)the relative larger SUW appears at the lower parts of design space, where the value of is smaller.

##### 3.3. Atlases of Global Conditioning Index

The condition number of the Jacobian matrix proposes a measure of the control accuracy of the manipulator, which is also used to evaluate kinematic isotropy and dexterity of the manipulator. The condition number should be kept as small as possible to obtain good kinematic performance. Its reciprocal is usually called the local conditioning index (LCI), which is configuration dependent. Further, the global conditioning index (GCI) is defined to evaluate the global behavior in the entire workspace, which can be expressed as where is the SUW of the parallel manipulator. And larger GCI value promises better control accuracy and kinematic performance generally.

The dimensional homogeneity of the Jacobian matrix is obtained by adopting the distribution radius of the end effector as the characteristic length. Each element in lower three rows of the Jacobian matrix is divided by the characteristic length. Then, the GCI atlas for our object manipulator is obtained and shown in Figure 5, from which we can conclude that(1)the design space can be subdivided into three regions with lines of and ;(2)in region I, GCI curves are almost parallel to the -axis, and GCI value grows with the increase of (or decrease of );(3)regions II and III are approximately symmetrical about the line of ;(4)the maximum GCI value appears around the midpoint of the line (, ).

#### 4. Dynamics and Inertia Matrix

The joint-space inertia matrix is the fundamental to deduce the inertia index, which can be derived from the dynamic model. In this paper, the virtual work principle [26] is adopted to establish the dynamic model of the object Stewart manipulator. The concise expression of the joint-space inertia matrix is obtained, based on the limb Jacobian matrices.

##### 4.1. Limb Jacobian Matrices

Limb Jacobian matrices describe the relation between limb velocities and the end-effector velocity. As shown in Figure 2, limb local frame is attached to the rotation center of the th universal joint, under which the inertia of the lower part limb is defined. Parallel limb local frame is located at the rotation center of the th spherical joint, where the inertia of the upper part limb is defined. The pose relation between limb local frame and corresponding limb is identical. The rotation matrix between the base frame and the frame can be obtained.

Equation (4) can be expressed in the limb local frame as where .

Accordingly, the sliding and angular velocities of the limb can be deduced from (6) in the limb local frame as Further, the linear velocity of lower and upper part centroids in the th limb can be expressed as where and are centroid position vectors. Expressions of limb Jacobian matrices , , and are listed in Appendix.

##### 4.2. Acceleration

Acceleration parameters can be deduced by taking the time derivation of (14), (15), (16), and (17) as The acceleration of end-effector mass center can be expressed in the base frame as where and is the mass center vector of the end effector. Expressions of , , , , and are given in Appendix.

##### 4.3. Dynamic Model

In this part, inertia force and torque of the end effector are given in the base frame , while the inertia force and torque of each limb are deduced in the corresponding limb frame . The inertial force and moment of the end effector exerted at the mass center can be expressed as where , is the moment of inertia of the end effector with respect to the mass center in , and are the external force and torque, and is the gravity.

The force and inertia exerted to the th lower part limb about its mass center can be deduced in the limb local frame as where is the lower part mass and is its inertia matrix about the mass center in .

The inertial force and moment of the upper part limb about the mass center can be written in as where is the upper part mass and is its inertia matrix about the mass center in .

Given a virtual displacement of the end effector, the dynamic equation can be obtained, according to the virtual work principle: where is the limb driving force, and and are corresponding virtual placements of lower and upper parts of the th limb in . Further, the joint-space inertia matrix can be obtained as where is a 3 × 3 unit matrix and relation of , , and is given in Appendix.

#### 5. Global Inertia Index and Atlases

The inertia index is the mean of eigenvalues of the joint-space inertia matrix, which can show the load characteristics. The inertia index is important to evaluate the dynamic performance of the parallel manipulator. In particular, a small value of the inertia index ensures that the manipulator could possess high mechanical resonance frequency, good dynamic accuracy, and high acceleration capacity.

The inertia index is also dependent on configuration of the manipulator and is a local performance index, which can be called as local inertia index (LII). In order to study the global performance of the manipulator within the workspace, the global inertia index (GII) can be defined as Parameters and are the eigenvalues and main diagonal elements of the joint-space inertia matrix , respectively.

In order to facilitate practical applications, the GII atlases of different parts are illustrated separately, such as the end effector, upper part limbs, and lower part limbs. The GII atlas for the end effector is shown in Figure 6, and the inertia and mass parameters of the end effector are given as and , for example. We can see that(1)the minimum value of GII appears in the place where , , and equal 0.86, 1.28, and 0.86, respectively;(2)in the middle region of the design space, there is a quite large area where the GII value is relatively small;(3)in most area, the GII value grows with the decrease of and the increase of difference between and ;(4)with the parameter point approaching the lower left corner or the right corner, the GII value grows sharply.

Let us assume that limbs of the Stewart manipulator are thin-walled cylinders. Then, mass and inertia parameters of upper and lower parts of the limb can be deduced with the given line density, which can be obtained from engineering experience and size estimation. Here the unit line destiny is adopted. GII atlases for the lower and upper parts of the limb are shown in Figures 7 and 8, respectively. We can discover that(1)distributions of the GII for the upper and lower parts are similar;(2)the manipulator with the best limb GII is in the middle of the design space, while the largest limb GII is located at the lower left and right corners.

#### 6. Atlases Based Optimum Design

Based on performance atlases obtained above, we can clearly see the performance trends and realize the performance of 3-3 Stewart manipulator with any combination of , , and . Besides, the rationality of the optimum goal can be recognized easily. In this section, the optimum design is carried out. The first step is to determine the performance requirements and to draw corresponding performance atlases, which has already been accomplished above.

The second step is to obtain the candidate region for the optimum design. Let us consider the kinematic performance only, at first. For example, if the zero-rotation SUW radius is not less than 0.3 and GCI value is not less than 0.55, then the candidate region can be determined as shown in Figure 9. Parameter ranges in this region can be obtained as , , and . In the candidate region, the manipulator can meet the desired workspace and condition number requirements.

Then, the inertia index is considered on the basis of the above kinematic requirements. Usually, the Stewart manipulator is under heavy duty, and the inertia parameter of the end effector (include the load) is dominant. Under this assumption, only Figure 6 is considered here. If the GII value is required less than or equal 4, the deduced candidate region is shown in Figure 10. Corresponding parameter ranges are , , and . Compared with the above result only considering kinematic performance, the candidate region is reduced, and parameter ranges shrink accordingly. In this candidate region, the manipulator has the desired dynamic performance as well as the desired kinematic performance.

The next step is to determine a set of optimum parameters, in view of the most desired performance or engineering conditions. In order to demonstrate the process, we consider the manipulator with , , and to obtain the minimum GII value within the candidate region, for example.

At the last step, the dimension factor is determined, and the actual parameters of the manipulator can be obtained. The dimension factor can be determined by the workspace requirement or some specific local index. For example, if the objective workspace radius is , the dimension factor can be obtained as . Then, the practical dimensional parameters can be figured out as , , and . With the given limb length ratio , limb length can also be determined.

#### 7. Conclusions

In this paper, optimum design approach for the parallel manipulator, which considers both kinematic and dynamic performances, is proposed. The nondimensional design space is established. In kinematics, the spherical usable workspace and the global conditioning index are investigated, and corresponding atlases are obtained. In dynamics, the global inertia index and its atlases are given, based on the deduced dynamic model. The atlases accurately and clearly illustrate variation trends of the manipulator performance with respect to design parameters. The global inertia index can take the practical load into account and is especially suitable for the optimization of the manipulator under heavy duty.

The analysis of the GII atlases shows that the manipulator with the best GII for the end effector is located in the upper middle region, while the maximum GII value appears at lower left and right corners. The GII distributions for the upper and lower part limbs are similar. Generally, the GII variation is more significant in the design space than that of the SUW radius or GCI.

Finally, on the basis of deduced atlases, candidate region is determined in the design space, and a set of parameters are obtained. Through adopting inertia index, this optimum design method can ensure that the parallel manipulator has the desirable kinematic and dynamic performances.

#### Appendix

Consider

#### Acknowledgments

This research is sponsored by the National Natural Science Foundation of China (nos. 51205224, 11178012) and the National Outstanding Youth Science Foundation (no. 51225503).

#### References

- J. P. Merlet,
*Parallel Robots*, Springer, 2nd edition, 2006. - J. Wu, T. Li, and L. Guan, “Computed-torque control for a moving flight simulator platform,”
*Journal of Tsinghua University*, vol. 46, no. 8, pp. 1405–1413, 2006. View at Scopus - X. Tang and R. Yao, “Dimensional design on the six-cable driven parallel manipulator of FAST,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 133, no. 11, Article ID 111012, pp. 1–12, 2011. View at Publisher · View at Google Scholar · View at Scopus - R. Yao, X. Tang, J. Wang, and P. Huang, “Dimensional optimization design of the four-cable-driven parallel manipulator in fast,”
*IEEE/ASME Transactions on Mechatronics*, vol. 15, no. 6, pp. 932–941, 2010. View at Publisher · View at Google Scholar · View at Scopus - X. Tang and Z. Shao, “Trajectory planning and tracking control of a multi-level hybrid support manipulator in FAST,”
*Mechatronics*, vol. 23, no. 8, pp. 1113–1122, 2003. - X. Tang, J. Wang, and M. Gao, “Kinematic calibration of gantry hybrid machine tool based on estimation error and local measurement information,”
*International Journal of Advanced Manufacturing Technology*, vol. 26, no. 4, pp. 382–390, 2005. View at Publisher · View at Google Scholar · View at Scopus - X.-J. Liu, Z.-L. Jin, and F. Gao, “Optimum design of 3-DOF spherical parallel manipulators with respect to the conditioning and stiffness indices,”
*Mechanism and Machine Theory*, vol. 35, no. 9, pp. 1257–1267, 2000. View at Publisher · View at Google Scholar · View at Scopus - X. Tang, D. Sun, and Z. Shao, “Structure and dimensional design of reconfigurable PKM,”
*International Journal of Advanced Robotic Systems*, vol. 10, pp. 1–10, 2013. - X.-J. Liu, J. Wang, and H.-J. Zheng, “Optimum design of the 5R symmetrical parallel manipulator with a surrounded and good-condition workspace,”
*Robotics and Autonomous Systems*, vol. 54, no. 3, pp. 221–233, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Bai, “Optimum design of spherical parallel manipulators for a prescribed workspace,”
*Mechanism and Machine Theory*, vol. 45, no. 2, pp. 200–211, 2010. View at Publisher · View at Google Scholar · View at Scopus - H.-B. Choi and J. Ryu, “Singularity analysis of a four degree-of-freedom parallel manipulator based on an expanded 6x6 Jacobian matrix,”
*Mechanism and Machine Theory*, vol. 57, pp. 51–61, 2012. - J. Enferadi and A. A. Tootoonchi, “Accuracy and stiffness analysis of a 3-RRP spherical parallel manipulator,”
*Robotica*, vol. 29, no. 2, pp. 193–209, 2011. View at Publisher · View at Google Scholar · View at Scopus - C. Gosselin, “Stiffness mapping for parallel manipulators,”
*IEEE Transactions on Robotics and Automation*, vol. 6, no. 3, pp. 377–382, 1990. View at Publisher · View at Google Scholar · View at Scopus - C. Gosselin and J. Angeles, “Optimum kinematic design of a spherical three-degree-of-freedom parallel manipulator,”
*Journal of mechanisms, transmissions, and automation in design*, vol. 111, no. 2, pp. 202–207, 1989. View at Scopus - Y. Zhao, F. Gao, X. Dong, and X. Zhao, “Dynamics analysis and characteristics of the 8-PSS flexible redundant parallel manipulator,”
*Robotics and Computer-Integrated Manufacturing*, vol. 27, no. 5, pp. 918–928, 2011. View at Publisher · View at Google Scholar · View at Scopus - M. Li, T. Huang, J. Mei, X. Zhao, D. G. Chetwynd, and S. J. Hu, “Dynamic formulation and performance comparison of the 3-DOF modules of two reconfigurable PKM—the tricept and the trivariant,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 127, no. 6, pp. 1129–1136, 2005. View at Publisher · View at Google Scholar · View at Scopus - T. Huang, J. Mei, Z. Li, X. Zhao, and D. G. Chetwynd, “A method for estimating servomotor parameters of a parallel robot for rapid pick-and-place operations,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 127, no. 4, pp. 596–601, 2005. View at Publisher · View at Google Scholar · View at Scopus - O. Peter, J. Hongzhou, H. Jingfeng, Y. Chifu, and H. Junwei, “Analysis of coupling effects on hydraulic controlled 6 degrees of freedom parallel manipulator using joint space inverse mass matrix,” in
*Proceedings of the 2nd International Conference on Intelligent Computing Technology and Automation (ICICTA '09)*, pp. 845–848, Changsha, China, October 2009. View at Publisher · View at Google Scholar · View at Scopus - Y. Yao, S. Fu, and L. Han, “Block diagonal dominance analysis and judgement of stewart platform's joint-space inertia matrix,”
*Chinese Journal of Mechanical Engineering*, vol. 44, no. 6, pp. 101–106, 2008. View at Scopus - G. Zhang and J. Furusho, “Speed control of two-inertia system by PI/PID control,” in
*Proceedings of the 3rd IEEE International Conference on Power Electronics and Drive Systems (PEDS '99)*, pp. 567–572, July 1999. View at Scopus - Z.-F. Shao, X. Tang, X. Chen, and L.-P. Wang, “Research on the inertia matching of the Stewart parallel manipulator,”
*Robotics and Computer-Integrated Manufacturing*, vol. 28, no. 6, pp. 649–659, 2012. - N. Rojas, J. Borras, and F. Thomas, “The octahedral manipulator revised,” in
*IEEE International Conference on Robotics and Automation*, pp. 2293–2298, St. Paul, Minn, USA, 2012. - B. Dasgupta and T. S. Mruthyunjaya, “Stewart platform manipulator: a review,”
*Mechanism and Machine Theory*, vol. 35, no. 1, pp. 15–40, 2000. View at Publisher · View at Google Scholar · View at Scopus - Ş. Staicu, “Dynamic analysis of the 3-3 Stewart platform,”
*UPB Scientific Bulletin D*, vol. 71, no. 2, pp. 3–18, 2009. View at Scopus - C. Gosselin and J. Angeles, “Singularity analysis of closed-loop kinematic chains,”
*IEEE Transactions on Robotics and Automation*, vol. 6, no. 3, pp. 281–290, 1990. View at Publisher · View at Google Scholar · View at Scopus - Z. Shao, X. Tang, X. Chen, and L. Wang, “Inertia match of a 3-RRR reconfigurable planar parallel manipulator,”
*Chinese Journal of Mechanical Engineering*, vol. 22, no. 6, pp. 791–799, 2009. View at Publisher · View at Google Scholar · View at Scopus