Download PDF
Research Article  |  Open Access  |  30 Nov 2021

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Views: 1897 |  Downloads: 858 |  Cited:  3
Intell Robot 2021;1(2):99-115.
10.20517/ir.2021.11 |  © The Author(s) 2021.
Author Information
Article Notes
Cite This Article

Abstract

This paper presents elastodynamic modeling and analysis for a five-axis lightweight robotic arm. Natural frequencies are derived and visualized within the dexterous workspace to show the overall performances and compare them to the frequencies when the robotics is with payload. The comparison shows that the payload has a relatively small influence to the first- and second-order frequencies. Sensitivity analysis is conducted, and the system's frequency is more sensitive to the second joint stiffness than the others. Moreover, observations from the displacement response analysis reveal that the robotics produces linear elastic displacements of the same level between the loaded and unloaded working modes but larger rotational deflections under the loaded working condition. The main contribution of this work lies in that a systematic approach of elastodynamic analysis for serial robotic manipulators is formulated, where the arm gravity and external load are taken into account to investigate the dynamic behaviors of the robotic arms, i.e., frequencies, sensitivity analysis, and displacement responses, under the loaded mode.

Keywords

Lightweight robotic arm, elastodynamics, natural frequency, displacement response

1. Introduction

Lightweight robotic arms and anthropomorphic assistive robots with high payload capacity are desired for applications of industry and welfare, among other fields, such as assisted daily living [1-3], pick-and-place operations [4], etc. Pick-and-place robots are well suited for a static environment where the task is repeated and precise tolerances are demanded[5]. As a mechanical system, the dynamic characteristics of the robotic arm is of importance to account for the requirements of application, such as high precision, speed, and payload. Henceforth, higher natural frequencies and low elastic displacements of a robotic manipulator will allow higher operational speeds and working cycles for efficient productivity[6]. Natural frequencies indicate the condition in which a mechanism tends to vibrate[7,8]. Differing from a structure or element, the dynamic behavior of a mechanism usually heavily depends on its architecture and configurations[9]; thus, it is not a trivial task to characterize the robot dynamics throughout the workspace, which calls for the kineto-elastodynamic analysis to provide the fundamentals of the modeling, design and control.

The elastodynamic modeling and analysis of a robotic manipulator have been reported previously[10,11], and they are roughly grouped into two categories: lumped modeling[12-15] and distributed-flexibilities modeling[9,16-19]. In general, with lumped modeling it is simpler to model the elastodynamic equation with acceptable computational accuracy, while the latter provides a more accurate model but with the high-dimensional generalized coordinate space and more complex procedure[20]. The commonly used method to study the elasticity of the robotic manipulators is the virtual joint method (VJM) as it can provide acceptable computation accuracy that is close to that of finite element analysis (FEA)[21]. Besides, VJM can be time efficient. VJM is based on pseudo-rigid body models with "virtual joints"[22-25]. Generally, the link flexibilities and linear/torsional springs take into account the bending contributions to the mechanism [26-29]. The stiffness formulated in the above approaches is limited to a subspace defined by the degrees of freedom (dofs) of the manipulator end-effector. Pashkevich et al.[30] overcame this issue by introducing a full-mobility lumped-parameter model by localizing 6-dof virtual springs to the links' ends and/or joints. In these models, the stiffness matrix is calculated in an unloaded equilibrium configuration of a robotic manipulator. On the other hand, the external loads directly influence the manipulator equilibrium configuration and, consequently, may modify the static properties. The lightweight design of the robotics accordingly decreases the link structural stiffness; thus, the robot geometry change due to external loads should be considered[31-33]. Consequently, elastodynamics of the robotic manipulators is an important concern in their design and applications. Based on the matrix structural analysis, Cammarata et al.[9,34] proposed an algorithm to assemble the stiffness matrix to investigate the manipulators with lower kinematic pairs. In this manner, the overall robotic manipulator in-parallel architecture can be split into substructures for modeling the elastodynamics[35]. Wu et al.[36] analyzed and compared the stiffness and natural frequencies of a 3-dof parallel manipulator with/without a redundant leg, where the joint deformations are ignored in the stiffness modeling. The small-amplitude deformations of the active joints can be considered as parameter uncertainties in terms of small variations to be integrated into the dynamic model[37]. Briot and Khalil[14] used the Newton-Euler recursive approach to develop a general symbolic elastodynamic calculation model for flexible parallel robots. Taghvaeipour et al.[15] derived the posture-dependent stiffness matrix in the elastodynamic modeling by resorting to the generalized spring concept. The previous models were established in the nominal configurations; hence, the geometry changes of the manipulator in this work are considered in the elastodynamic modeling and analysis.

In this paper, the elastodynamic characteristics of a lightweight robotic arm are investigated. The arm gravity and external load are taken into account to derive the stiffness matrix. Isocontours of natural frequencies over the dexterous workspace are formulated and sensitivity analysis is conducted. The frequencies and displacement responses of the robotics with payload are analyzed and compared with the dynamic behaviors of the unloaded mode. The main contribution of this work lies in that a systematic approach of elastodynamic analysis for serial robotic manipulators is formulated, where the arm gravity and external load are taken into account to investigate the dynamic behaviors of the robotic arms, i.e., frequencies, sensitivity analysis, and displacement responses, under the loaded mode.

2. Kinematics of the lightweight robotic arm

The lightweight robotic arm under study has five degrees of freedom (dof)[38], which adopts a modular design approach, as shown in Figure 1. The revolute joints are composed of CPU series gearboxes of Harmonic Drive and Maxon motor with gearhead to enhance the torque capabilities, except Joint 4 with geared motor. The actuators of joints are controlled by Maxon EPOS controllers. The Controller Area Network (CANopen) bus is adopted to build the communications between motors and controllers, and A CAN–USB interface is used to establish the communications between CANopen bus and the PC[38]. In accordance with the Denavit-Hartenberg (D-H) convention[39], the Cartesian coordinate systems are established accordingly.

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 1. The 5-dof lightweight robotic arm and its coordinate systems[38].

2.1. Kinematics of robotic arm

Throughout this work, i, j, and k stand for the unit vectors of the x-axis, y-axis, and z-axis, respectively. The transformation matrix in forward kinematics of the end-effector in reference frame is expressed as

$$ \begin{equation} \begin{aligned} { }^{0} \mathbf{A}_{5} = \left[\begin{array}{cc} \mathbf{R} & \mathbf{q} \\ \mathbf{0} & 1 \end{array}\right] = \prod_{i = 1}^{5}{ }^{i-1} \mathbf{A}_{i} ; \quad \quad { }^{i-1} \mathbf{A}_{i} = \left[\begin{array}{cc} { }^{i-1} \mathbf{R}_{i} & { }^{i-1} \mathbf{q}_{i} \\ \mathbf{0} & 1 \end{array}\right] \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} { }^{i-1} \mathbf{R}_{i} = \mathbf{R}\left(z_{i-1}, \theta_{i}\right) \mathbf{R}\left(x_{i}, \alpha_{i}\right) \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \begin{aligned} { }^{i-1} \mathbf{q}_{i} = \left[\begin{array}{lll} a_{i} \cos \alpha_{i} & a_{i} \sin \alpha_{i} & d_{i} \end{array}\right]^{T} \end{aligned} \end{aligned} \end{equation} $$

where D-H parameters are given in Table 1, and the inverse geometry problem for this robotics is well documented in the literature [8].

Table 1

D-H parameters of the 5-dof robotic arm

Joint iαiai[mm]di[mm]θi
1π/20250θ1
206000θ2
3π/200θ3
4–π/20600θ4
5π/20150θ5

2.2. Kinematic jacobian matrix

The velocities between the joints and end-effector are mapped with the Kinematic Jacobian matrix

$$ \begin{equation} \begin{aligned} \dot{\boldsymbol{\theta}} = \mathbf{J}^{-1} \mathbf{v}_{ef} \end{aligned} \end{equation} $$

where $$ \begin{array}{l} \dot{\theta}= [{\dot{\theta}}_1 & {\dot{\theta}}_2 & {\dot{\theta}}_3 & {\dot{\theta}}_4 & {\dot{\theta}}_5]^{T} \end{array}$$ denotes the joint angular velocities and $$ \begin{array}{l} \mathbf{v}_{ef} = [\omega ^{T} & \dot{\mathbf{q}}^{T}]^{T} \end{array}$$ is the velocity of the end-effector. Moreover, J is the kinematic Jacobian matrix of the robotic arm[40], namely,

$$ \begin{equation} \begin{aligned} \mathbf{J} = \left[\begin{array}{lllll}\mathbf{j}_{1} & \mathbf{j}_{2} & \mathbf{j}_{3} & \mathbf{j}_{4} & \mathbf{j}_{5}\end{array}\right] \quad \mathrm{where} \quad \mathbf{j}_{i} = \left[\begin{array}{c}\mathbf{z}_{i-1} \\ \mathbf{p}_{i-1} \times \mathbf{z}_{i-1}\end{array}\right] \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \mathbf{z}_{i-1}=\mathbf{R}_{i-1}\mathbf{k}; \quad \quad \mathbf{p}_{i-1}=\mathbf{q}_{i-1}-q \end{aligned} \end{equation} $$

where Ri–1 and qi–1 denote the rotation matrix and position vector of the transformation matrix from the reference coordinate system to the (i – 1)th coordinate system, respectively, which can be extracted from $$\Pi_{i=0}^{i-1} { }^{i-1} \mathbf{A}_{i}$$ in Equation (1).

2.3. Dexterous workspace

The reachable workspace of the robotic arm can be visualized by considering the limitation of the joint displacements and link dimensions. To effectively perform the kinematic performance, a dexterous workspace is defined, throughout which the inverse of the condition number of the Jacobian matrix is greater than 0.2, namely κ–1 (J) ≥ 0.2. Since the Jacobian matrix of Equation (4) is not homogeneous, a characteristic length[41] is introduced to normalize the Jacobian matrix as follows:

$$ \begin{equation} \begin{aligned} \begin{array}{l} \mathbf{j}_{i}^{\prime} = \left[\begin{array}{c} \mathbf{z}_{i-1} \\ \mathbf{p}_{i-1} \times \mathbf{z}_{i-1} / L \end{array}\right] ; \quad L^{2} = \frac{1}{5} \displaystyle{\sum_{i = 1}^{5}}\left\|\mathbf{p}_{i-1} \times \mathbf{z}_{i-1}\right\|^{2} \end{array} \end{aligned} \end{equation} $$

By constraining the condition number of the kinematic Jacobian matrix, a regular dexterous workspace is quarterly visualized in Figure 2.

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 2. The quarter of the reachable and dexterous workspace (red volume) for the robotic arm.

3. Elastodynamic model of robot

The elastodynamic modeling procedure pertains to the calculation of the stiffness and mass matrices of the manipulator, which is described in the following sections. Prior to the derivation of the elastodynamic model, the following assumptions are made:

  • The actuator stiffness is considered as an 1-dof torsional spring, while the link is considered as cantilever with a 6-dof spatial spring located at the end but treated as rigid.

  • The centers of mass of the regular components are coincident with their geometric centers.

  • The sum of moments of inertia of the actuators and Harmonic drivers are considered as lumped.

3.1. Stiffness matrix

To derive the elastodynamic equation of the robotic arm, the stiffness matrix is calculated with the virtual spring approach[42], based on the screw coordinates[43]. Hence, the component masses and external loads are taken into account to compute the Cartesian stiffness matrix. Figure 3 shows the VJM model of the robotic arm, where Gj, j = 1, 2, ..., 7, stands for the gravity and F for the external loads.

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 3. Virtual spring model of the 5-dof robotic arm with auxiliary loads, where Ac stands for the actuator and EE for end-effector.

Let θ and θ' be the original and the deformed displacements of the virtual springs, respectively, following the principle of virtual work, i.e. the work of the auxiliary forces is equal to the work of internal forces τθ, namely,

$$ \begin{equation} \begin{aligned} \sum\left(\mathbf{G}_{j}^{T} \delta \mathbf{t}_{j}\right)+\mathbf{F}^{T} \delta \mathbf{t} = \tau_{\theta}^{T}\left(\boldsymbol{\theta}^{\prime}-\boldsymbol{\theta}\right) \end{aligned} \end{equation} $$

where the virtual displacements δtj and δt can be computed from the linearized geometrical model derived from δtj =Jj(θ'θ) and δt =Jθ(θ'θ), respectively, Jj and Jθ being the Jacobians, namely,

$$ \begin{equation} \begin{aligned} \mathbf{J}_{\theta} = \left[\begin{array}{lllllll}\mathbf{j}_{1} & \mathbf{j}_{2} & \mathbf{J}_{u} & \mathbf{j}_{3} & \mathbf{J}_{l} & \mathbf{j}_{4} & \mathbf{j}_{5}\end{array}\right] \in \mathbb{R}^{6 \times 17} \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \mathbf{J}_{j} = \mathbf{J}_{\theta}(:, 1: k) \end{aligned} \end{equation} $$

where Jθ(:, 1 : k) stands for the first k columns in Jθ and k stands for the total degrees of freedom of the virtual springs from the base to Gj. Moreover,

$$ \begin{equation} \begin{aligned} \mathbf{J}_{u} = \left[\begin{array}{cccccc} \mathbf{x}_{1} & \mathbf{y}_{1} & \mathbf{z}_{1} & \mathbf{0} & \mathbf{0} & \mathbf{0} \\ \mathbf{q}_{2} \times \mathbf{x}_{1} & \mathbf{q}_{2} \times \mathbf{y}_{1} & \mathbf{q}_{2} \times \mathbf{z}_{1} & \mathbf{x}_{1} & \mathbf{y}_{1} & \mathbf{z}_{1} \end{array}\right] \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \mathbf{J}_{l} = \left[\begin{array}{cccccc} \mathbf{z}_{3} & \mathbf{x}_{3} & \mathbf{y}_{3} & \mathbf{0} & \mathbf{0} & \mathbf{0} \\ \mathbf{q}_{4} \times \mathbf{z}_{3} & \mathbf{q}_{4} \times \mathbf{x}_{3} & \mathbf{q}_{4} \times \mathbf{y}_{3} & \mathbf{z}_{3} & \mathbf{x}_{3} & \mathbf{y}_{3} \end{array}\right] \end{aligned} \end{equation} $$

Equation (7) is rewritten as

$$ \begin{equation} \begin{aligned} \sum\left(\mathbf{G}_{j}^{T} \mathbf{J}_{j}\left(\boldsymbol{\theta}^{\prime}-\boldsymbol{\theta}\right)\right)+\mathbf{F}^{T} \mathbf{J}_{\theta}\left(\boldsymbol{\theta}^{\prime}-\boldsymbol{\theta}\right)=\boldsymbol{\tau}_{\theta}^{T}\left(\boldsymbol{\theta}^{\prime}-\boldsymbol{\theta}\right) \end{aligned} \end{equation} $$

consequently, the force equilibrium equation is derived as

$$ \begin{equation} \begin{aligned} \tau_{\theta}=\sum\left(\mathbf{J}_{j}^{T} \mathbf{G}_{j}\right)+\mathbf{J}_{\theta}^{T} \mathbf{F}=\mathbf{J}_{g}^{T} \mathbf{G}+\mathbf{J}_{\theta}^{T} \mathbf{F} \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \mathbf{J}_{g}=\left[\begin{array}{llll} \mathbf{J}_{1}^{T} & \mathbf{J}_{2}^{T} & \ldots & \mathbf{J}_{7}^{T} \end{array}\right]^{T} ; \quad \mathbf{G}=\left[\begin{array}{llll} \mathbf{G}_{1}^{T} & \mathbf{G}_{2}^{T} & \ldots & \mathbf{G}_{7}^{T} \end{array}\right]^{T} \end{aligned} \end{equation} $$

Assuming that Kθ is the stiffness matrix in the joint space, with the linearized force–deflection relation, the equilibrium condition can be written as

$$ \begin{equation} \begin{aligned} \mathbf{J}_{g}^{T} \mathbf{G}+\mathbf{J}_{\theta}^{T} \mathbf{F}=\mathbf{K}_{\theta}\left(\boldsymbol{\theta}_{i}^{\prime}-\boldsymbol{\theta}_{i}\right) \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \mathbf{K}_{\theta}=\operatorname{diag}\left[\begin{array}{lllllll} K_{\text act,1 } & K_{\text act,2} & \mathbf{K}_{u} & K_{\text {act,3 }} & \mathbf{K}_{l} & K_{\text act,4} & K_{\text act,5} \end{array}\right] \end{aligned} \end{equation} $$

where Kact,i is the actuation stiffness and Ku and Kl are the upper and lower link stiffness matrices, respectively.

To calculate the stiffness matrix of the loaded mode, a neighborhood in the loaded configuration in which the external loads and the joint location are supposed to be incremented by small values δF and δθ, which can still satisfy the equilibrium conditions, is considered, leading to

$$ \begin{equation} \begin{aligned} \left(\mathbf{J}_{g}+\delta \mathbf{J}_{g}\right)^{T} \mathbf{G}+\left(\mathbf{J}_{\theta}+\delta \mathbf{J}_{\theta}\right)^{T}(\mathbf{F}+\delta \mathbf{F})=\mathbf{K}_{\theta}\left(\boldsymbol{\theta}^{\prime}-\boldsymbol{\theta}+\delta \boldsymbol{\theta}\right) \end{aligned} \end{equation} $$

and the linearized kinematic constraint

$$ \begin{equation} \begin{aligned} \delta \mathbf{t}=\mathbf{J}_{\theta} \delta \boldsymbol{\theta} \end{aligned} \end{equation} $$

Based on Equation (14), expanding Equation (16) yields

$$ \begin{equation} \begin{aligned} \mathbf{H}_{g}^{T} \otimes \mathbf{G} \delta \boldsymbol{\theta}+\mathbf{J}_{\theta}^{T} \delta \mathbf{F}+\mathbf{H}_{\theta}^{T} \otimes \mathbf{F} \delta \boldsymbol{\theta}=\mathbf{K}_{\theta} \delta \boldsymbol{\theta} \end{aligned} \end{equation} $$

where the symbol ⊗ represents the Kronecker product between matrices and Hg = ∂Jg/θ, Hθ = ∂Jθ/∂θ. Combining Equations (17) and (18), the stiffness model of the robotic manipulator is reduced to

$$ \begin{equation} \begin{aligned} {\left[\begin{array}{cc} \mathbf{0} & \mathbf{J}_{\theta} \\ \mathbf{J}_{\theta}^{T} & \mathbf{K}_{F}-\mathbf{K}_{\theta} \end{array}\right]\left[\begin{array}{c} \delta \mathbf{F} \\ \delta \boldsymbol{\theta} \end{array}\right]=\left[\begin{array}{c} \delta \mathbf{t} \\ \mathbf{0} \end{array}\right]} \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \mathbf{K}_{F}=\mathbf{H}_{g}^{T} \otimes \mathbf{G}+\mathbf{H}_{\theta}^{T} \otimes \mathbf{F} \end{aligned} \end{equation} $$

From δF = Kδt, the Cartesian stiffness matrix K of the robotic arm is calculated as

$$ \begin{equation} \begin{aligned} \mathbf{K}=\left(\mathbf{J}_{\theta}\left(\mathbf{K}_{\theta}-\mathbf{K}_{F}\right)^{-1} \mathbf{J}_{\theta}^{T}\right)^{-1} \end{aligned} \end{equation} $$

3.2. Mass matrix

The mass matrix can be derived from the expression of the systems kinetic energy, consisting of energies of the revolute joints, links, and end-effector. The energy of the five active joints are

$$ \begin{equation} \begin{aligned} E_{J}=\frac{1}{2}\left(\sum_{i=1}^{5} I_{\theta, i} \dot{\theta}^{2}+\sum_{i=3}^{5} m_{\theta, i} \mathbf{v}_{\theta, i}^{T} \mathbf{v}_{\theta, i}\right) \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \mathbf{v}_{\theta, 3}=\mathbf{E}_{3} \dot{\theta} ; \quad \mathbf{v}_{\theta, i}=\mathbf{E}_{45} \dot{\theta}, \quad i=4,5 \end{aligned} \end{equation} $$

and

$$ \begin{equation} \begin{aligned} \begin{array}{l} \mathbf{E}_{3}=\left[\begin{array}{lll} \mathbf{z}_{0} \times \mathbf{q}_{2} & \mathbf{z}_{1} \times\left(\mathbf{q}_{2}-\mathbf{q}_{1}\right) & \mathbf{0}_{3} \end{array}\right] \end{array} \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \begin{array}{l} \mathbf{E}_{45}=\left[\begin{array}{llll} \mathbf{z}_{0} \times \mathbf{q}_{4} & \mathbf{z}_{1} \times\left(\mathbf{q}_{4}-\mathbf{q}_{1}\right) & \mathbf{z}_{2} \times\left(\mathbf{q}_{4}-\mathbf{q}_{2}\right) & \mathbf{0}_{3 \times 2} \end{array}\right] \end{array} \end{aligned} \end{equation} $$

where Iθ,i is the moment of inertia of the ith joint, mθ,i is the mass, and vθ,i is the velocity in the Cartesian space. Let Iθ = diag[Iθ,1, Iθ,2, …, Iθ,5]; then, Equation (22) can be written in a compact form, namely,

$$ \begin{equation} \begin{aligned} E_{J}=\frac{1}{2} \dot{\theta}^{T} \mathbf{M}_{J} \dot{\theta} \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \mathbf{M}_{J}=\mathbf{I}_{\theta}+m_{\theta, 3} \mathbf{E}_{3}^{T} \mathbf{E}_{3}+\left(m_{\theta, 4}+m_{\theta, 5}\right) \mathbf{E}_{45}^{T} \mathbf{E}_{45} \end{aligned} \end{equation} $$

The kinetic energy of the upper/lower links and the wrist link can be expressed as

$$ \begin{equation} \begin{aligned} E_{L}=\frac{1}{2}\left(\mathbf{v}_{u}^{T} \mathbf{M}_{u} \mathbf{v}_{u}+\mathbf{v}_{l}^{T} \mathbf{M}_{l} \mathbf{v}_{l}+\mathbf{v}_{w}^{T} \mathbf{M}_{w} \mathbf{v}_{w}\right) \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \begin{aligned} \mathbf{M}_{u} =\left[\begin{array}{cc} \mathbf{R}_{1} \mathbf{I}_{u} \mathbf{R}_{1}^{T} & \mathbf{0} \\ \mathbf{0} & m_{u} \mathbf{1}_{3} \end{array}\right] \end{aligned} \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \begin{aligned} \mathbf{M}_{l} =\left[\begin{array}{cc} \mathbf{R}_{3} \mathbf{I}_{l} \mathbf{R}_{3}^{T} & \mathbf{0} \\ \mathbf{0} & m_{l} \mathbf{1}_{3} \end{array}\right] \end{aligned} \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \begin{aligned} \mathbf{M}_{w} =\left[\begin{array}{cc} \mathbf{R}_{4} \mathbf{I}_{w} \mathbf{R}_{4}^{T} & \mathbf{0} \\ \mathbf{0} & m_{w} \mathbf{1}_{3} \end{array}\right] \end{aligned} \end{aligned} \end{equation} $$

where the subscripted I, m, and v stand for the moment of inertia, mass, and velocities in the Cartesian space, respectively, and

$$ \begin{equation} \begin{aligned} \mathbf{v}_{u}=\mathbf{E}_{u} \dot{\theta} ; \quad \mathbf{v}_{l}=\mathbf{E}_{l} \dot{\boldsymbol{\theta}} ; \quad \mathbf{v}_{w}=\mathbf{E}_{w} \dot{\boldsymbol{\theta}} \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \begin{array}{l} \mathbf{E}_{u}=\left[\begin{array}{ccc} \mathbf{z}_{0} & \mathbf{z}_{1} & \mathbf{0}_{3} \\ \mathbf{z}_{0} \times \mathbf{q}_{u} & \mathbf{z}_{1} \times\left(\mathbf{q}_{u}-\mathbf{q}_{1}\right) & \mathbf{0}_{3} \end{array}\right] \end{array} \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \begin{array}{l} \mathbf{E}_{l}=\left[\begin{array}{cccc} \mathbf{z}_{0} & \mathbf{z}_{1} & \mathbf{z}_{2} & \mathbf{0}_{3 \times 2} \\ \mathbf{z}_{0} \times \mathbf{q}_{l} & \mathbf{z}_{1} \times\left(\mathbf{q}_{l}-\mathbf{q}_{1}\right) & \mathbf{z}_{2} \times\left(\mathbf{q}_{l}-\mathbf{q}_{2}\right) & \mathbf{0}_{3 \times 2} \end{array}\right] \end{array} \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \begin{array}{l} \mathbf{E}_{w}=\left[\begin{array}{ccccc} \mathbf{z}_{0} & \mathbf{z}_{1} & \mathbf{z}_{2} & \mathbf{z}_{3} & \mathbf{0}_{3 \times 1} \\ \mathbf{z}_{0} \times \mathbf{q}_{4} & \mathbf{z}_{1} \times\left(\mathbf{q}_{4}-\mathbf{q}_{1}\right) & \mathbf{z}_{2} \times\left(\mathbf{q}_{4}-\mathbf{q}_{2}\right) & \mathbf{z}_{3} \times\left(\mathbf{q}_{4}-\mathbf{q}_{2}\right) & \mathbf{0}_{3 \times 1} \end{array}\right] \end{array} \end{aligned} \end{equation} $$

where qu and ql are the position vector of the centers of the mass of the upper and lower links, respectively. Equation (27) can be cast in a matrix form as follows:

$$ \begin{equation} \begin{aligned} E_{L}=\frac{1}{2} \dot{\boldsymbol{\theta}}^{T} \mathbf{M}_{L} \dot{\boldsymbol{\theta}} \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \mathbf{M}_{L}=\mathbf{E}_{u}^{T} \mathbf{M}_{u} \mathbf{E}_{u}+\mathbf{E}_{l}^{T} \mathbf{M}_{l} \mathbf{E}_{l}+\mathbf{E}_{w}^{T} \mathbf{M}_{w} \mathbf{E}_{w} \end{aligned} \end{equation} $$

Similarly, the kinetic energy of the end-effector can be obtained as

$$ \begin{equation} \begin{aligned} E_{E}=\frac{1}{2} \mathbf{v}_{e f}^{T} \mathbf{M}_{E} \mathbf{v}_{e f} ; \quad \mathbf{M}_{E}=\left[\begin{array}{cc} \mathbf{R I}_{e} \mathbf{R}^{T} & \mathbf{0}_{3} \\ \mathbf{0}_{3} & m_{e} \mathbf{1}_{3} \end{array}\right] \end{aligned} \end{equation} $$

where Ie is the moment of inertia of the end-effector and me is the mass.

From the total kinetic energy of the robotic arm E = EJ + EL + EE, the mass matrix M for the robotic arm can be expressed as

$$ \begin{equation} \begin{aligned} \mathbf{M}=\mathbf{M}_{E}+\mathbf{J}^{-T}\left(\mathbf{M}_{J}+\mathbf{M}_{L}\right) \mathbf{J}^{-1} \end{aligned} \end{equation} $$

3.3. Dynamic equation and analysis

The dynamic equation of the robotic arm can be formulated as

$$ \begin{equation} \begin{aligned} \mathbf{M} \ddot{\mathbf{u}}+\mathbf{C} \dot{\mathbf{u}}+\mathbf{K} \mathbf{u}=\mathbf{f}-\mathbf{M} \dot{\mathbf{v}}_{e f}=\mathbf{F} \end{aligned} \end{equation} $$

where C is the damping matrix, F is the resultant force, and u and ü are the elastic displacement and acceleration, respectively. Since damping can only slightly influence the natural frequency and mode of free vibrations, the damp can be ignored to determine the natural frequencies. Simplification of Equation (35) results in the linearized elastodynamic equation below

$$ \begin{equation} \begin{aligned} \mathbf{M} \ddot{\mathbf{u}} + \mathbf{K} \mathbf{u}= 0 \end{aligned} \end{equation} $$

The rigidity of the system may be represented by the natural frequency The higher is the frequency the higher is the stiffness. From Equation (36), we get

$$ \begin{equation} \begin{aligned} \operatorname{det}\left(-\omega^{2} \mathbf{M}+\mathbf{K}\right)=0 \end{aligned} \end{equation} $$

where f = ω/2π denotes the natural frequency.

The displacement response analysis can be carried out from Equation (35) based on the initial conditions

$$ \begin{equation} \begin{aligned} \mathbf{u}_{0}=\mathbf{u}(0) ; \quad \dot{\mathbf{u}}_{0}=\dot{\mathbf{u}}(0) \end{aligned} \end{equation} $$

Here, the damping ratios are set to ς = 6% according to the manipulator structure. From Equation (37), the displacement vector u can be represented in terms of the modal contributions, namely,

$$ \begin{equation} \begin{aligned} \mathbf{u}=\mathbf{Q}\boldsymbol{\eta} \end{aligned} \end{equation} $$

where Q and η are the modal matrix and the vector of the displacements in each mode, respectively. Consequently, Equation (35) can be rewritten as

$$ \begin{equation} \begin{aligned} \ddot{\boldsymbol{\eta}}+\boldsymbol{\Phi} \dot{\boldsymbol{\eta}}+\boldsymbol{\Omega} \boldsymbol{\eta}=\mathbf{f}_{d} \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \boldsymbol{\Phi}=\mathbf{Q}^{T} \mathbf{C Q}=\operatorname{diag}\left[\begin{array}{llll} 2 \varsigma \omega_{1} & 2 \varsigma \omega_{2} & \ldots & 2 \varsigma \omega_{6} \end{array}\right] \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \boldsymbol{\Omega}=\mathbf{Q}^{T} \mathbf{K} \mathbf{Q}=\operatorname{diag}\left[\begin{array}{llll} \omega_{1}^{2} & \omega_{2}^{2} & \ldots & \omega_{6}^{2} \end{array}\right] \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \mathbf{f}_d=\mathbf{Q}^{T}\mathbf{F} \end{aligned} \end{equation} $$

Since the mass and stiffness matrices in Equation (35) are time-varying, the common way to solve such a problem is to divide the motion period into extremely short intervals, where the stiffness and mass matrices are considered as constant in each interval. Let T denote the complete motion period that is divided into N intervals, namely, Δt = T/N. In the nth time interval τ ∈ [tn–1, tn], the equation of motion in the ith mode is expressed as

$$ \begin{equation} \begin{aligned} \ddot{\eta}_{i}+2 \varsigma \omega_{i} \dot{\eta}_{i}+\omega_{i}^{2} \eta_{i}=f_{d i} \end{aligned} \end{equation} $$

Thus, the ith mode contributes to the displacement response[44] is

$$ \begin{equation} \begin{aligned} \eta_{i}\left(t_{n}\right)= & e^{\varsigma \omega_{i} \Delta t}\left(\cos \omega_{d i} \Delta t+\frac{\varsigma}{\sqrt{1-\varsigma^{2}}} \sin \omega_{d i} \Delta t\right) \eta_{i}\left(t_{n-1}\right) \\ & +\frac{1}{\omega_{d i}} \int_{t_{n-1}}^{t_{n}} f_{d i}(\tau) e^{-\varsigma \omega_{i}\left(t_{n}-\tau\right)} \sin \omega_{d i}\left(t_{n}-\tau\right) \mathrm{d} \tau \\ & +\left(\frac{1}{\omega_{d i}} e^{\varsigma \omega_{i} \Delta t} \sin \omega_{d i} \Delta t\right) \dot{\eta}_{i}\left(t_{n-1}\right) \end{aligned} \end{equation} $$

where

$$ \begin{equation} \begin{aligned} \omega_{d i}=\omega_{i} \sqrt{1-\varsigma^{2}} \end{aligned} \end{equation} $$

Differentiating Equation (43) with respect to time leads to

$$ \begin{equation} \begin{aligned} \dot{\eta}_{i}\left(t_{n}\right)=\dot{\eta}_{i, 1}\left(t_{n}\right)+\dot{\eta}_{i, 2}\left(t_{n}\right)+\dot{\eta}_{i, 3}\left(t_{n}\right) \end{aligned} \end{equation} $$

with

$$ \begin{equation} \begin{aligned} \dot{\eta}_{i, 1}\left(t_{n}\right)=e^{\varsigma \omega_{i} \Delta t}\left(\frac{2 \varsigma^{2}-1}{\sqrt{1-\varsigma^{2}}} \omega_{i} \sin \omega_{d i} \Delta t+\frac{2 \varsigma}{\sqrt{1-\varsigma^{2}}} \omega_{d i} \cos \omega_{d i} \Delta t\right) \eta_{i}\left(t_{n-1}\right) \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \dot{\eta}_{i, 2}\left(t_{n}\right)=\frac{1}{\omega_{d i}} e^{\varsigma \omega_{i} \Delta t}\left(\varsigma \omega_{i} \sin \omega_{d i} \Delta t+\omega_{d i} \cos \omega_{d i} \Delta t\right) \dot{\eta}_{i}\left(t_{n-1}\right) \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \dot{\eta}_{i, 3}\left(t_{n}\right)=\frac{1}{\omega_{d i}} \int_{t_{n-1}}^{t_{n}} f_{d i}(\tau) e^{-\varsigma \omega_{i}\left(t_{n}-\tau\right)}\left(\varsigma \omega_{i} \sin \omega_{d i}\left(t_{n}-\tau\right)-\omega_{d i} \cos \omega_{d i}\left(t_{n}-\tau\right)\right) \mathrm{d} \tau \end{aligned} \end{equation} $$

Hence, ηi(tn) and $$\dot{\eta}_{i}(t_{n})$$ can be solved as long as ηi(tn–1) and $$\dot{\eta}_{i}(t_{n-1})$$ are given, and

$$ \begin{equation} \begin{aligned} \eta_{i}(0)=\mathbf{e}_{i}^{T} \mathbf{M u}(0) ; \quad \dot{\eta}_{i}(0)=\mathbf{e}_{i}^{T} \mathbf{M} \dot{\mathbf{u}}(0) \end{aligned} \end{equation} $$

where ei is the ith column of the modal matrix. The total displacement response is calculated by the following addition

$$ \begin{equation} \begin{aligned} \mathbf{u}\left(t_{n}\right)=\sum_{i=1}^{6} \eta_{i}\left(t_{n}\right) \mathbf{e}_{i}\left(t_{n}\right) \end{aligned} \end{equation} $$

Consequently, the natural frequency and displacement response can be obtained with numerical calculations.

4. Numerical simulation

Elastodynamic characteristics of the robotic arm are investigated in this section. The properties of the robotics components are listed in Tables 2 and 3, respectively. Moreover, according to the output shaft of the gearbox, the actuation stiffnesses are calculated and set to Kact,i = 2 · 104 Nm/rad, i = 1, ..., 5, and the link stiffness matrices given in Appendix A are derived by means of FEA with ANSYS[45]. The numerical simulation was carried out with Matlab.

Table 2

Mass and moment of inertia of the active joints

Joint i12345
Iθ,i [kg · mm2]0.02100.00020.00010.00010.0002
mθ,i [kg]-2.22721.81962.24422.0053
Table 3

The properties of the links and end-effector

LinksMass [kg]Moment of inerita [kg cm2]
upper linkmu = 4.7995Iu = diag[1.1884, 25.0670, 24.4940]
lower linkmt= 1.7795It = diag[4.0802, 4.0861, 0.2345]
wrist link-Iw = diag[0.5556, 0.9154, 0.6119
end-effectorme = 1.2961Ie = diag [0.4563, 0.4382, 0.2347]

4.1. Natural frequency

To effectively measure the overall performance of the robotic arm, the distributions of natural frequencies over the dexterous workspace in Figure 2 are visualized, as displayed in Figures 4 and 5.

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 4. The natural frequency with constant-orientation [0, 0, 0] (in unit of rad): (a) first order; (b) second order. (The color bar stands for the numerical value of the term in the legend, which is applicable to Figs. 5 to 7.)

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 5. The natural frequency with constant-orientation [0, π/20]: (a) first order; and (b) second order.

Let the end-effector orientation follow the ZXZ Euler convention; the distributions of the first- and second-order natural frequencies over workspace are displayed in Figures 4 and 5 when the end-effector remains vertical and horizontal, respectively. It can be observed that the nonsymmetric distributions of the natural frequencies in Figure 5 are different from the symmetric ones in Figure 4. This is because the robot configurations are not axisymmetric about the vertical direction with the vertical end-effector, leading to different inverse kinematic solutions of such a 5-dof robotic arm, which are different from the axisymmetric robot configurations with horizontal end-effector. As the mass and stiffness matrices of the robot are configuration dependent, non-symmetric distributions of natural frequencies in Figure 5 occur. These two figures show that the first two orders of natural frequencies increase with the increasing z coordinates but with decreased x and y coordinates, namely both the first and second frequencies increases from the workspace boundaries to the origin of the global coordinate systems. As displayed in Figure 4, when the end-effector remains vertical, the natural frequencies have the same varying trend in any vertical cross-section of the workspace. By contrast, the first- and second-order frequencies become smaller counterclockwise within the workspace when the end-effector is in the horizontal configuration, as shown in Figure 5. Moreover, it is found that the differences among the frequencies of the manipulator in different configurations are not so large, which means that the robotic arm has close frequencies inside the overall workspace.

4.2. Sensitivity analysis

Sensitivity analysis can be used to evaluate the influence of the geometric parameters and design variables to the manipulator performances. Based on the elastodynamic equation, there exists

$$ \begin{equation} \begin{aligned} \left(-\omega_{i}^{2} \mathbf{M}+\mathbf{K}\right) \mathbf{e}_{i}=0 \end{aligned} \end{equation} $$

Upon differentiation of Equation (49), the derivative equation with respect to a variable δ is obtained as follows:

$$ \begin{equation} \begin{aligned} \left(-2 \omega_{i} \frac{\partial \omega_{i}}{\partial \delta} \mathbf{M}-\omega_{i}^{2} \frac{\partial \mathbf{M}}{\partial \delta}+\frac{\partial \mathbf{K}}{\partial \delta}\right) \mathbf{e}_{i}+\left(-\omega_{i}^{2} \mathbf{M}+\mathbf{K}\right) \frac{\partial \mathbf{e}_{i}}{\partial \delta}=0 \end{aligned} \end{equation} $$

Taking the dot-product on both sides of Equation (50) yields

$$ \begin{equation} \begin{aligned} \mathbf{e}_{i}^{T}\left(-2 \omega_{i} \frac{\partial \omega_{i}}{\partial \delta} \mathbf{M}-\omega_{i}^{2} \frac{\partial \mathbf{M}}{\partial \delta}+\frac{\partial \mathbf{K}}{\partial \delta}\right) \mathbf{e}_{i}+\mathbf{e}_{i}^{T}\left(-\omega_{i}^{2} \mathbf{M}+\mathbf{K}\right) \frac{\partial \mathbf{e}_{i}}{\partial \delta}=0 \end{aligned} \end{equation} $$

From

$$ \begin{equation} \begin{aligned} \mathbf{e}_{i}^{T} \mathbf{M e}_{i}=1 ; \quad \mathbf{e}_{i}^{T}\left(-\omega_{i}^{2} \mathbf{M}+\mathbf{K}\right)=\left(\left(-\omega_{i}^{2} \mathbf{M}+\mathbf{K}\right) \mathbf{e}_{i}\right)^{T}=0 \end{aligned} \end{equation} $$

we have

$$ \begin{equation} \begin{aligned} -2 \omega_{i} \frac{\partial \omega_{i}}{\partial \delta}-\omega_{i}^{2} \mathbf{e}_{i}^{T} \frac{\partial \mathbf{M}}{\partial \delta} \mathbf{e}_{i}+\mathbf{e}_{i}^{T} \frac{\partial \mathbf{K}}{\partial \delta} \mathbf{e}_{i}=0 \end{aligned} \end{equation} $$

or

$$ \begin{equation} \begin{aligned} \frac{\partial \omega_{i}}{\partial \delta}=-\frac{1}{2 \omega_{i}}\left(-\omega_{i}^{2} \mathbf{e}_{i}^{T} \frac{\partial \mathbf{M}}{\partial \delta} \mathbf{e}_{i}+\mathbf{e}_{i}^{T} \frac{\partial \mathbf{K}}{\partial \delta} \mathbf{e}_{i}\right) \end{aligned} \end{equation} $$

Figure 6 illustrates the sensitivity of the first-order natural frequency to the first two active joints with constant orientation [0, π/2, 0]. It is found that the first-order natural frequency is much more sensitive to the second joint, particularly in the upper and lower workspace regions, which implies that the robots dynamic performance can be improved by replacing the second joint with a stiffer actuator. It is noted that the distributions of sensitivity coefficients are not symmetric, which is because the robot configurations are not axisymmetric about the vertical direction when the robot end-effector moves with some constant orientations, since the robot under study is a 5-dof robotic arm. Moreover, if a payload with more mass were exerted to the robot, it could be predicted that the sensitivity coefficients will be increased with very tiny varying trends, compared to the present results.

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 6. Sensitivities of the first-order natural frequency to the joint stiffness: (a) Joint 1; and (b) Joint 2.

4.3. Dynamic analysis of loaded system

With the payload 5 kg applied to the end-effector of the robotic arm, they constitute a new dynamic system and the solved frequencies with constant-orientation [0, π/2, 0] are illustrated in Figure 7, from which it is observed that the frequencies of the loaded robotic system decrease about 20% compared to Figure 5. Table 4 lists the average frequencies[46] within the constant-orientation workspace defined by

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 7. The frequencies with payload at constant-orientation [0, π/2, 0]: (a) first order; and (b) second order.

Table 4

The mean frequency (Hz) within the dexterous workspace

$$\bar{f}_{1}$$$$\bar{f}_{2}$$$$\bar{f}_{3}$$$$\bar{f}_{4}$$$$\bar{f}_{5}$$$$\bar{f}_{6}$$
Natural frequency7.94048.957919.026384.9292136.1418300.6579
Frequency with payload5.95386.535815.277255.899886.6291276.5956

$$ \begin{equation} \begin{aligned} \bar{f}_{i}=\frac{\int f_{i} \mathrm{~d} \Omega}{\int \mathrm{d} \Omega} \end{aligned} \end{equation} $$

where Ω stands for the workspace volume. Different from the traditional industrial robots with low frequencies, the high order frequencies have large values to make the manipulator achieve high-speed motion. Compared to the average natural frequencies, the frequencies of the robotics with payload reduce 10%–40% for the six orders of frequencies. From the view of kineto-elastodynamic characteristics, the difference between the frequency of the loaded system and its natural frequency could be a consideration in the design of the mechanical system, where the smaller difference implies higher rigidity and higher payload capability.

Assuming that the motion of the robotic arm follows the trajectory (unit: mm) defined by

$$ \begin{equation} \begin{aligned} \begin{array}{l} x=750+750(\cos \pi \tau-1) \\ y=750(1-\cos \pi \tau) \\ z=600(1-\cos \pi \tau) \end{array} \end{aligned} \end{equation} $$

where the end-effector keeps constant-orientation [0, π, 0] and the motion period T = 0.5s is divided into 1024 intervals, Figure 8 shows the displacement responses of the end-effector, from which it is seen that the linear elastic displacement responses are close, whenever the robotic arm is under loaded and unloaded working modes. The angular displacements of the end-effector generate relatively large differences. The largest deformations appear around 0.3 s where the end-effector is located in the middle layer of the workspace, approximately z = 250 mm.

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 8. Displacement responses of the end-effector: (a) x direction; (b) y direction; (c) z direction; (d) ϕx direction; (e) ϕy direction; and (f) ϕz direction.

Figure 9 shows the comparison of the joint angular displacements between the numerical simulation and experimental measurements along previous trajectory, where the experimental data are read from the motor encoders. Due to the frictions and time-varying disturbance in the joints, the experimental curve profiles have more fluctuations and larger vibration amplitudes than the simulation ones. On the other hand, the comparison shows that the differences between these two curves are small, thus, the built analytical model can be acceptable for dynamic analysis of the robots.

On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm

Figure 9. Comparison of the joint angular displacements under loaded mode.

5. Conclusion

This paper presents the elastodynamic characteristics of a 5-dof lightweight robotic arm. The main contribution is that a systematic approach of elastodynamic analysis for serial robotic manipulators is formulated, where the arm gravity and external load are taken into account to investigate the dynamic behaviors of the robotic arms, i.e., frequencies, sensitivity analysis, and displacement responses, with auxiliary payloads exerted to the robot. The modeling in this work eases the evaluation of elastodynamics of the manipulator at a large number of postures as the elastodynamic aspect is usually time-consuming. As the mass and stiffness matrices are posture dependent, the proposed method can effectively provide a symbolic calculation and achieve the modal analysis along an operating trajectory. Moreover, such a model can compute the additional mass or evaluate the influence of an isolator to the system more precisely to eliminate/reduce vibration in the vibration control. The developed model can be used in either performance evaluation or design optimization.

The frequencies of the loaded robotics are visualized within the representative workspace regions to show the overall dynamic performance and compare them with the natural frequencies. The comparison reveals that the studied robot keeps relatively high rigidity with high payload ratio. It is found from sensitivity analysis that the natural frequency can effectively increase by improving the second joint stiffness. Based on the displacement responses analysis, the payload has a slight influence on the translational elastic displacements of this robotic system, although it leads to reduced frequencies, while the effect on the rotation deflections cannot be ignored. In the future, the developed model will be integrated into its control system and an optimum redesign of the robotics will be conducted.

Declarations

Authors' contributions

The author contributed solely to the article.

Availability of data and materials

Not applicable.

Financial support and sponsorship

This work was supported by Natural Science Foundation of Liaoning Province (Grant No. 20180520028).

Conflicts of interest

The author declared that there are no conflicts of interest.

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Copyright

© The Author(s) 2021.

Supplementary Materials

APPENDIX A

REFERENCES

1. Ivlev O, Martens C, Graeser A. Rehabilitation Robots FRIEND-I and FRIEND-II with the dexterous lightweight manipulator. Technology and Disabilily 2005;17:111-23.

2. Bien Z, Chung MJ, Chang PH, Kwon DS. Integration of a Rehabilitation Robotic System (KARES II) with Human-Friendly Man-Machine Interaction Units. Auton Robot 2004;16:165-91.

3. Mahoney RM. The Raptor wheelchair robot system. In: Mokhtari M, editor. Integration of Assistive Technology in the Information Age. IOS press; 2001. pp. 135-41.

4. Universal Robots. Available from https://www.universal-robots.com/GB/Cases.aspx.

5. Wu G, Shen H. In: Introduction. Singapore: Springer Singapore; 2021. pp. 1-15.

6. Hoevenaars AG, Krut S, Herder JL. Jacobian-based natural frequency analysis of parallel manipulators. Mech Mach Theory 2020;148:103775.

7. Briot S, Pashkevich A, Chablat D. On the optimal design of parallel robots taking into account their deformation and natural frequencies. In: ASME IDETC & CIE Conf. vol. DETC2009-86230. San Diego, California, USA: 2009. pp. 367-76.

8. Siciliano B, Khatib O. Springer Handbook of Robotics. Springer; 2016.

9. Cammarata A, Condorelli D, Sinatra R. An algorithm to study the elastodynamics of parallel kinematic machines with lower kinematic pairs. ASME J Mech Robot 2013;5:011004.

10. Dwivedy SK, Eberhard P. Dynamic analysis of flexible manipulators, a literature review. Mech Mach Theory 2006;41:749-77.

11. Briot S, Khalil W. In: Ceccarelli M, editor. Dynamics of Parallel Robots, vol. 35 of Mechanisms and Machine Science. Springer International Publishing AG Switzerland; 2015.

12. Khalil W, Gautier M. Modeling of mechanical systems with lumped elasticity. In: Proceedings 2000 ICRA. Millennium Conference. IEEE International Conference on Robotics and Automation. Symposia Proceedings (Cat. No. 00CH37065). vol. 4. San Francisco, CA, USA: IEEE; 2000. pp. 3964-69.

13. Wittbrodt E, Adamiec-Wójcik I, Wojciech S. Dynamics of Flexible Multibody Systems: Rigid Finite Element Method. Foundations of Engineering Mechanics. Springer Science & Business Media; 2007.

14. Briot S, Khalil W. Recursive and symbolic calculation of the elastodynamic model of flexible parallel robots. Int J Robot Res 2014;33:469-83.

15. Taghvaeipour A, Angeles J, Lessard L. Elastodynamics of a two-limb Schönflies motion generator. Proc Ins Mech Eng Part C J Mech Eng Sci 2015;229:751-64.

16. Boyer F, Coiffet P. Symbolic modeling of a flexible manipulator via assembling of its generalized Newton-Euler model. Mech Mach Theory 1996;31:45-56.

17. Rognant M, Courteille E, Maurine P. A systematic procedure for the elastodynamic modeling and identification of robot manipulators. IEEE Trans Robot 2010;26:1085-93.

18. Bauchau OA. Flexible Multibody Dynamics, vol. 176 of Solid Mechanics and Its Applications. Springer Science & Business Media; 2010.

19. de Jalon JG, Bayo E. Kinematic and Dynamic Simulation of Multibody Systems: the Real-time Challenge. Springer Science & Business Media; 2012.

20. Wu G, Shen H. In: Ding H, Sun R, editors. Parallel PnP Robots, vol. 7 of Research on Intelligent Manufacturing. Springer, Singapore; 2021.

21. Ku DM, Chen LW. Kineto-elastodynamic vibration analysis of robot manipulators by the finite element method. Comput Struct 1990;37:309-17.

22. Salisbury JK. Active stiffness control of a manipulator in cartesian coordinates. In: 1980 19th IEEE Conference on Decision and Control including the Symposium on Adaptive Processes. Albuquerque, NM, USA: 1980. pp. 95-100.

23. Gosselin C. Stiffness mapping for parallel manipulators. IEEE Trans Robot Autom 1990;6:377-82.

24. Wittbrodt E, Adamiec-Wójcik I, Wojciech S. Dynamics of Flexible Multibody Systems. Springer; 2006.

25. Quennouelle C, Gosselin C. Kinematostatic modeling of compliant parallel mechanisms. Meccanica 2011;46:155-69.

26. El-Khasawneh BS, Ferreira PM. Computation of stiffness and stiffness bounds for parallel link manipulators. Int J Mach Tool Manuf 1999;39:321-42.

27. Gosselin CM, Zhang D. Stiffness analysis of parallel mechanisms using a lumped model. Int J Robot Autom 2002;17:17-27.

28. Dai J, Ding X. Compliance analysis of a three-legged rigidly-connected platform device. ASME J Mech Des 2006;128:755-64.

29. Majou F, Gosselin C, Wenger P, Chablat D. Parametric stiffness analysis of the Orthoglide. Mech Mach Theory 2007;42:296-311.

30. Pashkevich A, Chablat D, Wenger P. Stiffness analysis of overconstrained parallel manipulators. Mech Mach Theory 2009;44:966-82.

31. Kövecses J, Angeles J. The stiffness matrix in elastically articulated rigid-body systems. Multi Syst Dyn 2007;18:169-84.

32. Quennouelle C, Gosselin CM. Stiffness Matrix of Compliant Parallel Mechanisms. In: Lenarčič J, Wenger P, editors. Advances in Robot Kinematics: Analysis and Design. Springer Netherlands; 2008. pp. 331-41.

33. Tyapin I, Hovland G. Kinematic and elastostatic design optimisation of the 3-DOF Gantry-Tau parallel kinematic manipulator. Modeling, Identification and Control 2009;30:39-56. DOI

34. Cammarata A, Caliò I, D'Urso D, et al. Dynamic stiffness model of spherical parallel robots. J Sound Vib 2016;384:312-24.

35. Wu L, Wang G, Liu H, Huang T. An approach for elastodynamic modeling of hybrid robots based on substructure synthesis technique. Mech Mach Theory 2018;123:124-36.

36. Wu J, Li T, Wang J, Wang L. Stiffness and natural frequency of a 3-DOF parallel manipulator with consideration of additional leg candidates. Robot Auton Syst 2013;61:868-75.

37. Lara-Molina FA, Koroishi EH, Costa TL. Elastodynamic Performance of a Planar Parallel Mechanism Under Uncertainties. In: International Symposiu on Multibody Systems and Mechatronics. Springer; 2017. pp. 183-92.

38. Wu G, Zhao W, Zhang X. Optimum time-energy-jerk trajectory planning for serial robotic manipulators by reparameterized quintic NURBS curves. Proc Ins Mech Eng Part C J Mech Eng Sci 2021;235:4382-93. DOI

39. Denavit J, Hartenberg RS. A kinematic notation for lower-pair mechanisms based on matrices. ASME J Appl Mech 1955;22:215-21.

40. Tsai LW. Robot Analysis: The Mechanics of Serial and Parallel Manipulators. John Wiley & Sons; 1999.

41. Ranjbaran F, Angeles J, González-Palacios MA, Patel RV. The mechanical design of a seven-axes manipulator with kinematic isotropy. J Intell Robot Syst 1995;14:21-41.

42. Pashkevich A, Klimchik A, Chablat D. Enhanced stiffness modeling of manipulators with passive joints. Mech Mach Theory 2011;46:662-79.

43. Wu G, Bai S, Kepler J. Mobile platform center shift in spherical parallel manipulators with flexible limbs. Mech Mach Theory 2014;75:12-26.

44. Rao SS. Mechanical Vibrations. 4th ed. Prentice Hall; 2003.

45. Dong C, Liu H, Huang T, Chetwynd DG. A screw theory-based semi-analytical approach for elastodynamics of the Tricept robot. ASME J Mech Robot 2019;11:031005.

46. Alessandro C, Rosario S. Elastodynamic optimization of a 3T1R parallel manipulator. Mech Mach Theory 2014;73:184-96.

Cite This Article

Research Article
Open Access
On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm
Guanglei WuGuanglei  Wu

How to Cite

Wu, G. On the elastodynamics of a five-axis lightweight anthropomorphic robotic arm. Intell. Robot. 2021, 1, 99-115. http://dx.doi.org/10.20517/ir.2021.11

Download Citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click on download.

Export Citation File:

Type of Import

Tips on Downloading Citation

This feature enables you to download the bibliographic information (also called citation data, header data, or metadata) for the articles on our site.

Citation Manager File Format

Use the radio buttons to choose how to format the bibliographic data you're harvesting. Several citation manager formats are available, including EndNote and BibTex.

Type of Import

If you have citation management software installed on your computer your Web browser should be able to import metadata directly into your reference database.

Direct Import: When the Direct Import option is selected (the default state), a dialogue box will give you the option to Save or Open the downloaded citation data. Choosing Open will either launch your citation manager or give you a choice of applications with which to use the metadata. The Save option saves the file locally for later use.

Indirect Import: When the Indirect Import option is selected, the metadata is displayed and may be copied and pasted as needed.

About This Article

Special Issue

© The Author(s) 2021. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Data & Comments

Data

Views
1897
Downloads
858
Citations
3
Comments
0
31

Comments

Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at support@oaepublish.com.

0
Download PDF
Share This Article
Scan the QR code for reading!
See Updates
Contents
Figures
Related
Intelligence & Robotics
ISSN 2770-3541 (Online)
Follow Us

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/