Source-linked AI summary
Enhanced stiffness modeling of manipulators with passive joints
Anatoly Pashkevich, Alexandr Klimchik, Damien Chablat
TL;DR
The paper addresses stiffness analysis for serial and parallel manipulators with passive joints under external and internal loading. It proposes a nonlinear VJM-based method that computes loaded equilibria, stability, and potentially rank-deficient Cartesian stiffness, with Ortholide examples demonstrating loading-dependent nonlinear behavior and buckling-related changes.
Problem
Conventional stiffness analysis often uses unloaded configurations, omitting loading-induced changes in manipulator geometry, Jacobians, and Hessians.
Method
The method models links and actuators as pseudo-rigid bodies connected by multidimensional virtual springs and perfect passive joints, then linearizes loaded force-deflection relations near equilibrium.
Results
The examples demonstrate force- and torque-dependent stiffness, sudden stiffness changes beyond critical loading, and configurations associated with buckling.
Takeaways & Limitations
Loaded-mode stiffness effects and uncommon elastic behaviors should be included in manipulator design and accuracy analysis.
Takeaways & Limitations
The method assumes linear link elasticity and concentrated loading at the end-effector, limiting its scope for large link deflections or joint and distributed loads.
Abstract
from arXiv · showhide
The paper presents a methodology to enhance the stiffness analysis of serial and parallel manipulators with passive joints. It directly takes into account the loading influence on the manipulator configuration and, consequently, on its Jacobians and Hessians. The main contributions of this paper are the introduction of a non-linear stiffness model for the manipulators with passive joints, a relevant numerical technique for its linearization and computing of the Cartesian stiffness matrix which allows rank-deficiency. Within the developed technique, the manipulator elements are presented as pseudo-rigid bodies separated by multidimensional virtual springs and perfect passive joints. Simulation examples are presented that deal with parallel manipulators of the Ortholide family and demonstrate the ability of the developed methodology to describe non-linear behavior of the manipulator structure such as a sudden change of the elastic instability properties (buckling).
1 Introduction
Manipulator stiffness is central to positioning accuracy, safe dynamic operation, and application-specific compliance. Existing stiffness analysis commonly neglects loading effects, motivating a nonlinear loaded-mode approach that can capture buckling.
- Manipulator stiffness governs deformation and positioning errors under external, inertial, and gravitational loading across industrial and medical applications.
- Service robots instead require low stiffness and compliant architectures to reduce collision risks during human interaction.
- Stiffness analysis commonly uses a Cartesian stiffness matrix relating small translational or rotational displacements to applied forces or torques.
- FEA, SMA, and VJM provide complementary stiffness-modeling approaches, trading physical detail, computational cost, and parametric transparency.
- Conventional analysis focuses on unloaded equilibrium, so it cannot generally capture loading-induced configuration changes, buckling, or other nonlinear elastic phenomena.
2 Related works
Prior work developed VJM-based stiffness models for various manipulator architectures, including passive joints, but loaded-mode analysis remains insufficiently general. The unresolved challenges include nonlinear equilibrium, stability, geometrical stiffness, and rank-deficient Cartesian stiffness.
- VJM extends rigid manipulator models with localized virtual springs representing distributed and lumped flexibility, enabling preliminary-stage Cartesian stiffness analysis.
- Existing passive-joint methods can produce rank-deficient Cartesian stiffness matrices and avoid manual elimination of redundant springs through dedicated matrix inversion.
- Prior parallel-manipulator treatments mainly address pure architectures whose stiffness can be assembled from independently analyzed serial chains.
- Loaded-mode stiffness requires accounting for geometrical stiffness caused by configuration changes, whereas most prior work assumes quasi-static unloaded configurations.
- The literature identifies a need for robotic numerical techniques that compute potentially non-unique loaded equilibria, assess stability, and locally linearize force-deflection relations.
3 Problem of stiffness modeling
The paper models serial and parallel manipulators with flexible links, actuated joints, and passive joints using virtual springs and explicit passive-joint coordinates. It formulates loaded stiffness as a nonlinear force-deflection problem followed by local linearization and stability analysis.
- 3.1 Basic Assumptions: The studied manipulators combine flexible actuated joints and links with passive joints, forming serial chains that can serve as components of parallel architectures.
- 3.1 Basic Assumptions: The VJM represents link and actuator elasticity with localized virtual springs while retaining perfect passive joints as separate kinematic variables.
- 3.1 Basic Assumptions: Preloaded passive joints are modeled with auxiliary springs whose torque or force depends linearly on deviation from a non-loaded joint coordinate.
- 3.2 Problem statement: The stiffness model maps external force or torque to end-effector deflection and analyzes singularities that may induce nonlinear phenomena such as buckling.
- 3.2 Problem statement: Loaded-mode analysis evaluates force increments near a load-shifted equilibrium, then linearizes the relation to compute a Cartesian stiffness matrix that may be singular.
- 3.2 Problem statement: A central task is determining critical forces that can cause buckling or sudden changes in the loaded manipulator configuration.
4 Static equilibrium for the loaded mode
The loaded-mode analysis computes nonlinear static equilibria before linearizing the load–deflection relation, allowing stability and buckling behavior to be assessed under external and internal loading.
- Static equilibrium determines the loaded manipulator configuration required to linearize the load–deflection relation.
- The equilibrium formulation treats external loading and configuration as unknowns for a prescribed end-effector location while incorporating internal preloading.
- Virtual-work conditions combine passive-joint equilibrium with spring reactions and the Jacobians mapping generalized coordinates to end-effector motion.
- The numerical algorithm linearizes kinematics and equilibrium iteratively using Jacobians, the geometry function, and the current configuration.
- 3–5 iterations are usually sufficient for typical deformations, but buckling and multiple equilibria can make convergence highly dependent on the initial guess.
- The dual problem is meaningful only when the stiffness matrix is nonsingular; it is always singular for a separate serial chain with passive joints but usually nonsingular for parallel manipulators.
- Stability is assessed through the loaded stiffness relation, with stable configurations requiring the relevant matrix to be positive-definite.
5 Stiffness matrix for the loaded mode
The loaded-mode stiffness matrix is obtained by linearizing neighboring loaded equilibria, so configuration-dependent Jacobians and Hessians capture external and internal loading effects.
- The method computes a nonlinear elastic force–deflection relation and derives the engineering stiffness matrix by linearizing it near a loaded configuration.
- The linearized kinematic relation expresses end-effector displacement through Jacobians with respect to the flexible and passive coordinates.
- Expanding Jacobian differentials with Hessians yields a coupled linear system whose direct inverse provides the Cartesian stiffness matrix.
- Internal preloading affects the equilibrium coordinates, Jacobians, and Hessians even though it is not included explicitly in the supplementary linear relations.
- The approach also computes supplementary stiffness matrices that measure the sensitivity of joint coordinates to external loading.
- For parallel manipulators, the resulting Cartesian stiffness can be assembled by summing the contributions of separate kinematic chains.
- The technique generalizes previous serial-chain and parallel-manipulator results and computes Cartesian stiffness under external and internal loading.
6 Illustrative examples
The examples show that loaded stiffness analysis captures nonlinear force–deflection behavior, configuration changes, bifurcations, and buckling in serial chains and the Orthoglide manipulator.
- Serial-chain examples: The serial-chain examples compare passive-joint stiffness models across S-, Π-, and Z-configurations under external loading.The chain includes passive universal and spherical joints plus actuated rotating joints, and the study evaluates multiple virtual-spring assumptions.
- Model A: For the initial S-configuration, the stable straight posture persists below F ≤ K_θ/L, becomes unstable for K_θ/L < F ≤ 3K_θ/L, and may coexist with unstable straight and zig-zag equilibria above 3K_θ/L.The behavior is analogous to axial compression of a straight column, and the numerical algorithm reproduces the stable-equilibrium curve.
- Serial-chain examples: Π- and Z-configurations are mostly nonlinear but generally retain one stable and one unstable equilibrium; Π can exhibit stable-equilibrium bifurcation, while Z can buckle into a nonsymmetrical posture.For the Z-configuration, four stable and two unstable equilibria may differ in potential energy.
- Model B: Model B exhibits high initial stiffness under axial compression, followed by stiffness loss after buckling; Z-configuration buckles near F_0 = 1.07K_θ/L.Its qualitative behavior resembles Model A, but the stiffness coefficient is slightly lower because the elastic elements are arranged differently.
- Model C: Model C predicts lower critical forces and buckling in the z-direction for all examined postures, including a 25-fold stiffness reduction for Π above 0.20K_θ/L.Critical forces are 0.16K_θ/L for S and 0.17K_θ/L for Z; the Π result represents buckling absent from Models A and B.
- Conclusions from examples: The developed technique detects and numerically evaluates loading-induced nonlinear effects, while conventional unloaded analysis remains reasonable before buckling but cannot identify the applicable small-deflection loading range.For Orthoglide, geometric buckling occurs at a loading 2.5 times below the Euler local-link buckling force.
- Parallel-manipulator example: For Orthoglide points Q0–Q4, force–deflection relations are nonlinear and contain radical stiffness reductions at buckling, with critical forces from 2.2 to 7.8 kN and deflections from 0.5 to 4.9 mm.The isotropic point Q0 has a critical force of 4.6 kN and critical deflection of 1.4 mm; another buckling mode arises from the spatial arrangement of the parallel chains.
7 Discussion
The approach incorporates loading-dependent Jacobians and Hessians into a nonlinear stiffness model, improving accuracy and enabling buckling detection. Its VJM formulation also offers efficient load-deflection tracing, while retaining linear link-elasticity assumptions and scope restrictions.
- Loading-dependent Jacobians and Hessians improve stiffness-model accuracy and enable detection of buckling phenomena.
- The method accounts for nonlinear manipulator geometry but models each link's elasticity with linear expressions, a fundamental VJM assumption.
- The loading is assumed concentrated at the end-effector center, and architectures are restricted to strictly parallel serial chains without cross-linking.The paper leaves joint or distributed loading and more general architectures for further work.
- The technique efficiently traces complete load-deflection paths and supports iterative assessment of critical force and deflection.Its computational efficiency comes from reducing matrix size through the VJM representation.
- The technique captures uncommon nonlinear behavior relevant to robot accuracy and therefore applicable to manipulator design and analysis.
8 Conclusions
The proposed method models stiffness in loaded manipulators with passive joints, computes equilibrium, stability, and load-deflection behavior, and linearizes Cartesian stiffness. Examples reveal force-dependent stiffness changes and buckling-related risks, while the current scope is pure parallel mechanisms with similar chains and actuator placement.
- 8 Conclusions: The method represents links as pseudo-rigid bodies with multidimensional virtual springs and computes loaded equilibrium, stability, and full load-deflection paths.Its Cartesian stiffness linearization uses inversion of a larger dedicated matrix.
- 8 Conclusions: Examples on Ortholide parallel manipulators and associated serial chains demonstrate stiffness dependence on applied force or torque and sudden changes beyond critical wrench values.The examples also detect serial-chain configurations potentially dangerous with respect to buckling.
- 8 Conclusions: The method is applied to pure parallel mechanisms with similar kinematic chains and actuators located between the base and foot.It can be extended to other actuator locations and dissimilar chain geometries.
- 8 Conclusions: Future work targets cross-linked and more sophisticated parallel architectures, as well as heavy manipulators with gravity loading distributed within their links.These directions extend the current application scope.