- 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 2010 (2010), Article ID 143961, 12 pages

http://dx.doi.org/10.1155/2010/143961

## Modeling of Configuration-Dependent Flexible Joints for a Parallel Robot

^{1}Department of Mechanical and Materials Engineering, Queen's University, Kingston, ON, Canada K7L 3N6^{2}Department of Aerospace Engineering, Ryerson University, Toronto, ON, Canada M5B 2K3

Received 7 November 2008; Revised 13 April 2009; Accepted 12 July 2009

Academic Editor: Zhen Huang

Copyright © 2010 Zili Zhou 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

This paper provides a method to determine the variable flexible joint parameters which are dependent on configurations for a PRS Parallel Robot. Based on the continuous force approach, virtual springs were used between the joint components to simulate the joint flexibility. The stiffness matrix of the joint virtual springs was derived. The method uses system dynamic characteristics in different configurations to set the virtual spring stiffness for all the joints in the system. Modal testing was conducted on a set of selected robot configurations to determine the system natural frequencies and mode shapes along with their variation. To obtain the virtual spring stiffness, the system was condensed at the joint nodal coordinates. Then eigen-sensitivity analysis was conducted on the condensed system matrices with respect to the stiffness parameters of the joint virtual springs. Thus, the virtual spring parameters in the model can be set to match the variation of the system dynamic responses with the robot configuration changes. The virtual spring parameters between the selected robot configurations were obtained by interpolation. The research indicates that the method is effective and relatively easy to conduct, compared to other methods. The variable flexible joint model is applicable to flexible multibody systems with variable configurations.

#### 1. Introduction

Dynamic modeling of flexible multibody systems is a classical problem. The link flexibility and joint compliance (elasticity) of robots or mechanisms in general have been studied for a long time. This paper focuses on the difficult issue of joint flexibility which varies due to system configuration changes. Joint flexibility has been studied since the 1980s, mainly by relating it to contact/impact mechanics. Tzou et al. [1] proposed a stochastic approach to model the random feature of the dynamic contacts in the joints. Bowden and Dugundji [2] presented the linear and nonlinear analysis for the global dynamics of jointed structures. Moon and Li [3] conducted an experimental study on a pin-jointed truss. Tzou and Rong [4] provided a mathematical model for a three-dimensional spherical joint based on the contact force analysis. Kakizaki et al. [5] presented a dynamic modeling method for a SCARA manipulator with link flexibility and joint clearances. Seneviratne et al. [6] provided a combined massless-link and spring-damping model for modeling the joint clearance. Ravn et al. [7] presented the analysis of revolute joint clearances with and without lubricant. Schwab et al. [8] presented a study on the dynamic response of mechanisms affected by the joint clearance. Ting et al. [9] presented a novel and simple approach to identify the worst position and direction errors due to the joint clearance of a single loop linkage. Zhu and Ting [10] studied the uncertainty of planar and spatial robots with joint clearances. Wang et al. [11] presented the virtual spring method at joints to completely decouple the dynamic model of complex robotic systems with closed kinematic chains.

In the above investigations, two basic approaches were used in the dynamics of flexible joint systems. The first approach uses an impulse momentum model, where the pieced intervals are analyzed. The impulse momentum equations can be solved with the restitution condition for the jump discontinuity in the system velocities and the joint reaction forces. The second approach uses a continuous force model to represent the force of interaction between the impact surfaces. Stiffness and damping coefficients are set to account for the impact surface compliance and energy dissipation during the impact process. Obviously, the second approach is easier to handle than the first. No matter which method is used, a key step is to determine the parameters in the joint models. Though the above models are effective to different extents when used to represent the physics of an individual joint in a configuration-fixed system, they are difficult to use for multiple joints assembled together in a configuration-variable system (such as a robot). The difficulty is that the joint flexibility parameters may change when the system configuration changes. Also, the parameters may be different for the same type of joints at different locations in a system. It is the highlight for this paper to display how to set the flexibility parameters for multiple joints in a system with variable configuration.

Based on the continuous force model, this paper presents a modeling method on the joint whose flexibility results from the variable joint looseness due to the system configuration change. The flexible joint model was used in the simulation of a PRS parallel robot shown in Figure 1. Instead of studying the joint individually, the method uses system dynamic characteristics to set the flexibility parameters for all the joints in the system. Experimental modal testing was conducted at a set of selected robot configurations to determine the variation of the system natural frequencies and mode shapes. To easily obtain the joint flexibility parameters, the system was condensed at the joint coordinates. At each selected robot configuration, eigen-sensitivity analysis was conducted on the condensed system matrices with respect to the joint flexibility parameters. Thus, the joint flexibility parameters in the model can be updated at each selected robot configuration to be consistent with the variation of the system dynamic responses. The joint flexibility parameters between the selected configurations were obtained by interpolation.

#### 2. Flexible Joint Model

Figure 2 shows the system model of a PRS parallel robot, where is the global reference, is the platform body reference, and is the leg body reference, . As a continuous force model, all the joint constraints are released and replaced by virtual springs and virtual dampers correspondingly, except for the rotational constraint around in the revolute joints. Figures 3(a) and 3(b) show the flexible revolute joint model and the related coordinate system, where , , , and are the stiffness parameters of the virtual springs at the revolute joints in their original constrained directions: along , and around , respectively, . Similarly (not shown), let , , and be the stiffness parameters of the virtual springs at the spherical joints in the directions of the leg body reference : along , , and , respectively, . Altogether there are 21 stiffness parameters of the joint virtual springs, or simply called joint stiffness parameters, in the robot system model. In the revolute joint, the rotational constraint around (leg axis in undeformed state) is still kept, because releasing of this constraint is kinematically meaningless for the tripod motion (configuration determining), and the kinetic energy of this rotation is small enough to omit (compared to the rotation around ).

In this case it is recognized that the special revolute joint (not related to the leg body) is not a “general” revolute joint even for the flexible case, and so the symbol “R” in the phrase “PRS parallel robot” may not seem perfectly suitable. Rather than use a special symbol (such as ), we have maintained the use of R in order to reduce unnecessary confusion. For simplicity, this research considers the joint virtual dampers with modal damping in the entire system; only the joint virtual springs are modeled here.

The potential energy of the joint virtual springs in the th branch can be originally written as

where , , , and are the nodal coordinates of the leg at the revolute joint of the th branch along , , , and about directions respectively; , , and are the nodal coordinates of the leg at the spherical joint of the th branch along , , and directions respectively; , , and are the nodal coordinate projections in leg body reference from , , and , the nodal coordinates of the platform at the spherical joint of the th branch along , , and directions, respectively, and the projections can be calculated as
in which is the transformation matrix from platform body reference to global reference *OXYZ*, and from leg body reference to global reference *OXYZ*. Substituting (2) into (1) leads to the matrix form

where are the associated joint nodal coordinates in the th branch

and is the joint stiffness matrix in the th branch

where denotes , and denotesin which and are the revolute and spherical joint stiffness matrices, respectively, in the th branch:

The joint coordinates in the th branch can be expressed in terms of the system total nodal coordinates by joint coordinate connectivity matrix :

where the elements of matrix are 1 and 0. Thus, by substituting (7) into (3), the system joint stiffness matrix is obtained by summing the joint stiffness matrix from all three branches:

By writing the system potential energy, the system stiffness matrix can be obtained as the summation of system link stiffness matrix and system joint stiffness matrix , both associated with the system total nodal coordinates :

In (6), the joint stiffness parameters will be set variable to match the variation of the system dynamic responses at different configurations. Therefore, (5), (8), and (9) are variable matrices.

Using static condensation method given by Guyan [12], the condensed system stiffness matrix and mass matrix are obtained as follows: where is system mass matrix, and and are partitioning matrices with elements of 1 and 0 which divide system total nodal coordinates as master coordinates , mainly consisted of the joint nodal coordinates, and slave coordinates as follows:

The above and are associated with the master coordinates .

#### 3. Eigen-Sensitivity on the Condensed System

According to Fox and Kapoor [13], the natural frequency sensitivity and mode shape sensitivity with respect to joint stiffness parameter , *b* = 1 to 7, to 3, can be respectively, derived on the condensed system as

where *f _{i}* is the

*i*th natural frequency, is the

*i*th normalized mode shape vector, and . Using (9) and noticing that and are constant lead to Therefore, the terms and in (13) and (14) can be calculated using (10) and (11) as In light of (8), the term in (16) and (17) is where can be calculated using (5) as

where denotes and denotes in which, by using (6) and (7), the differentials of the revolute and spherical joint stiffness matrices with respect to the joint stiffness parameters in the th branch are

#### 4. Modal Testing and Simulation

MATLAB codes were written based on the above analysis. The link model was built by using finite elements and was used to investigate the joint flexibility. For coordinate condensation, the selected master coordinates are the linear joint nodal coordinates in (4) in three branches (removing angular coordinate due to the difficulty to measure), plus three linear nodal coordinates at the centre of the platform, where the cutting tool is supposed to be attached, along , , and directions. Thus there are 30 master coordinates in total. All other coordinates are set slave. Table 1 shows the physical parameters of the robot. Table 2 lists the two sets of robot configurations that were used in the modal testing and simulation. They feature both asymmetric and symmetric robot configurations.

Modal testing was conducted at each selected robot configuration. In order to understand the entire robot prototype, all measurement points were chosen to be evenly distributed on the platform and three legs. The excitation points were chosen close to the spherical joints on the moving platform. Accelerometers were used to detect the responses in the local body reference directions at the measurement points of each leg and the platform. The excitation and response signals were amplified and then recorded using LabVIEW where the FRF of each measurement is generated. The sampling frequency was 2000 Hz and sampling time was 1 second. The window function was set as force-exponential with 50 force window and 10 exponential window. The FRF measurements were imported to the postprocessing software ME’scope which extracts the system natural frequencies and operating deflection shapes (ODSs) that are theoretically close to the mode shapes.

With eigen-sensitivity (13) and (14) at hand, the natural frequencies and mode shapes in the model can be modified at the selected robot configurations according to the modal testing results by adjusting the joint stiffness parameters as follows:

where represents the modification values of the 21 joint stiffness parameters; represents the errors of natural frequencies between the measurement and the calculation; , to , represents the errors of components in the th mode shape between the measurement and the calculation; is the matrix whose elements are the values from (13); , to , is the matrix whose elements are the vector element values from (14); is the eigen-solution order number; the superscript + stands for the Moore-Penrose pseudoinverse. Based on (21), an iteration procedure is formed to update the joint stiffness parameters until the errors of natural frequencies and mode shapes are within specified limits.

Figures 4(a), 4(b), and 4(c) show the stiffness parameters of the joint virtual springs for the robot configurations from 661 to 667, and Figures 5(a), 5(b), and 5(c) from 261 to 267 The initial joint stiffness parameters are uniformly set as 10^{5} N/ (or N*/rad for , ) for all configurations. By using (21) for iteration, the joint stiffness parameters are obtained at the modal testing configurations 661 to 667 and 261 to 267, and linearly interpolated into the simulation configurations between these modal testing configurations. These figures indicate that the joint stiffness parameters are set configuration-dependent in order to match the variation of the system dynamic responses.

Figure 6 displays and compares the variations of the first 15 calculated and measured natural frequencies with the slider (prismatic joint) position changes starting at configuration 661 and ending at configuration 667 (referring to Table 1), and Figure 7 starting at configuration 261 and ending at configuration 267. The differences between the calculated and the measured frequencies are mostly less than 1, and the largest deviations are less than 7 at the 5th and 6th natural frequencies at configurations 661 and 267. The possible reasons for the deviations are that the solutions of the joint stiffness parameters are locally optimal and that these two frequencies may be related with some slave coordinates.

Figures 8(a), 8(b), and 8(c), respectively, display and compare the 5th, 7th, and 13th mode shapes at the node of the platform center (omitting mode shapes at other nodes) at configuration 661 to 667, and Figures 9(a), 9(b), and 9(c) at configuration 261 to 267 These mode shapes have the largest displacements at the platform center compared to other mode shapes: the 5th mode shape has the largest displacement in direction, the 7th in direction, and the 13th in direction. For easy comparison, the selected mode shape components are normalized. These figures show that the calculated mode shapes are close to the measured ones although they do not perfectly agree. The disagreement may result from the undamped joint model on which the joint stiffness parameters are calculated. For the simplicity of the method, however, it is worthy of the small loss of mode shape accuracy.

It must be indicated that, for the ideal rigid (tight) joint model of the same robot, natural frequencies and mode shapes will also be variable with the configuration changing (not shown here), but the natural frequencies are much higher than those of the flexible joint model. In simulation of the ideal joint model, the first natural frequency is about 200 Hz and the second is about 500 Hz. Therefore, Figures 6–9 display the flexible joint effect on the system due to configuration changes.

#### 5. Conclusions

This study presents a configuration-dependant flexible joint model for a parallel robot. The method is based on adaptation of virtual springs between the joint components to simulate the joint flexibility. The joint stiffness matrix of virtual springs was derived. The system was condensed at the joint nodal coordinates. Eigen-sensitivity analysis was conducted on the condensed system matrix with respect to the stiffness parameters of the joint virtual springs. Dynamic modification was conducted at a series of robot configurations for the virtual spring parameters to be set variable to match the variation of the system natural frequencies and mode shapes obtained from modal testing. The virtual spring parameters between the selected robot configurations were obtained by interpolation. The research indicates that the presented method is effective and relatively easy to conduct, compared to other methods. The variable flexible joint model is applicable to flexible multibody systems with variable configurations.

#### References

- H. S. Tzou, Y. Rong, and A. Nassirharand, “Design and stochastic simulation of elastically jointed mechanical systems,”
*Computer-Aided Design*, vol. 21, no. 7, pp. 435–440, 1989. - M. Bowden and J. Dugundji, “Joint damping and nonlinearity in dynamics of space structures,”
*American Institute of Aeronautics and Astronautics*, vol. 28, no. 4, pp. 740–749, 1990. - F. C. Moon and G. X. Li, “Experimental study of chaotic vibrations in a pin-jointed space truss structure,”
*American Institute of Aeronautics and Astronautics*, vol. 28, no. 5, pp. 915–921, 1990. - H. S. Tzou and Y. Rong, “Contact dynamics of a spherical joint and a jointed truss-cell system,”
*American Institute of Aeronautics and Astronautics*, vol. 29, no. 1, pp. 81–88, 1991. - T. Kakizaki, J. F. Deck, and S. Dubowsky, “Modeling the spatial dynamics of robotic manipulators with flexible links and joint clearances,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 115, no. 4, pp. 839–847, 1993. - L. D. Seneviratne, S. W. E. Earles, and D. N. Fenner, “Analysis of a four-bar mechanism with a radially compliant clearance joint,”
*Journal of Mechanical Engineering Science*, vol. 210, no. 3, pp. 215–223, 1996. - P. Ravn, S. Shivaswamy, B. J. Alshaer, and H. M. Lankarani, “Joint clearances with lubricated long bearings in multibody mechanical systems,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 122, no. 4, pp. 484–488, 2000. - A. L. Schwab, J. P. Meijaard, and P. Meijers, “A comparison of revolute joint clearance models in the dynamic analysis of rigid and elastic mechanical systems,”
*Mechanism and Machine Theory*, vol. 37, no. 9, pp. 895–913, 2002. View at Publisher · View at Google Scholar - K.-L. Ting, J. Zhu, and D. Watkins, “Effects of joint clearance on position and orientation deviation of linkages and manipulators,”
*Mechanism and Machine Theory*, vol. 35, no. 3, pp. 391–401, 2000. View at Publisher · View at Google Scholar - J. Zhu and K.-L. Ting, “Uncertainty analysis of planar and spatial robots with joint clearances,”
*Mechanism and Machine Theory*, vol. 35, no. 9, pp. 1239–1256, 2000. View at Publisher · View at Google Scholar - J. Wang, C. M. Gosselin, and L. Cheng, “Modeling and simulation of robotic systems with closed kinematic chains using the virtual spring approach,”
*Multibody System Dynamics*, vol. 7, no. 2, pp. 145–170, 2002. View at Publisher · View at Google Scholar - R. J. Guyan, “Reduction of stiffness and mass matrices,”
*American Institute of Aeronautics and Astronautics*, vol. 3, no. 2, p. 380, 1965. - R. L. Fox and M. P. Kapoor, “Rates of change of eigenvalues and eigenvectors,”
*American Institute of Aeronautics and Astronautics*, vol. 6, no. 12, pp. 2426–2429, 1968.