cs.ROAug 28, 2026

Dual Park-Ravani Interpolation of Rigid Motions: Acceleration-Field Continuity and Holonomic Hermite Repair

Authors: Daniel Condurache

Organizations: “Gheorghe Asachi” Technical University of Ia¸i, Romania

Read PDF
Ask AI
Cite

Copy citation

DOINot available for this paper

Abstract

The Park-Ravani construction generates a twice continuously differentiable, frame-invariant spline on SO(3) by exponentiating cubic canonical-coordinate polynomials. We show that the construction transfers, without changing form, to the group of orthogonal dual tensors, a representation of rigid displacements. The transferred recurrence is stated compactly through the dual extension of the right Jacobian of the exponential map and its first Fréchet derivative. This yields interpolation of prescribed rigid poses and continuity of the body dual twist and its first derivative. Using the higher-order rigid-body kinematics of dual spatial twists, we then prove that the resulting curve has a continuous physical acceleration field, not merely a continuous quantity obtained by formally differentiating the dual part of a twist. We also distinguish algebraic dual transfer from temporal differential prolongation: their simultaneous first-order use takes place in a hyper-dual algebra, and interpolation of arbitrary prolonged nodal data need not be holonomic. A noncommuting three-pose example verifies the recurrence, all knot continuity statements, and dimensional covariance under a change from meters to millimeters. We define and analyze the first-order holonomy defect of a generic hyper-dual interpolant, exhibit an exact counterexample, and remove the defect by cubic or quintic Hermite interpolation in dual logarithmic coordinates.

Figures & tables

VerificationResidual
Compact versus expanded recurrence8.33×10−178.33\times 10^{-17}
Dual pose interpolation8.05×10−168.05\times 10^{-16}
Body dual twist continuity1.36×10−161.36\times 10^{-16}
Body dual-twist derivative continuity2.29×10−162.29\times 10^{-16}
Spatial dual twist continuity2.00×10−162.00\times 10^{-16}
Spatial dual-twist derivative continuity2.29×10−162.29\times 10^{-16}
Table 1: Residuals at the internal knot t1t_{1} .
Pointaxa_{x}aya_{y}aza_{z}
p1p_{1}1.336338-1.7924701.773506
p1+0.25e1p_{1}+0.25e_{1}1.284059-1.7544101.928334
p1+0.20e2p_{1}+0.20e_{2}1.256146-1.8409931.788991
p1+0.18e3p_{1}+0.18e_{3}1.267878-1.8440181.728332
Table 2: Physical point accelerations at t1t_{1} in m s−2\mathrm{m\,s^{-2}} .
ResidualInitial endpointFinal endpoint
Cubic pose008.05×10−168.05\times 10^{-16}
Cubic body twist7.14×10−137.14\times 10^{-13}4.04×10−124.04\times 10^{-12}
Quintic pose008.05×10−168.05\times 10^{-16}
Quintic body twist1.59×10−131.59\times 10^{-13}1.73×10−121.73\times 10^{-12}
Quintic body-twist derivative1.08×10−101.08\times 10^{-10}6.32×10−96.32\times 10^{-9}
Table 3: Independent curve-level verification of the Hermite interpolants.
More in this topic

Lie Group Methods

Lie group methods use transformation groups and their Lie algebras to model rotations, rigid motions, symmetries, and other structured dynamics.

See all 42 papers

Explore similar work

CardsList
Aug 27, 2026

Arbitrary-Order Hermite Interpolation of Rigid-Motion Jets via Hyper-Multidual Quaternions

We study bilateral interpolation of finite-order rigid-motion jets represented by unit dual quaternions. An order-nn multidual (MD) algebra is the truncated polynomial algebra R[ε]/(εn+1)\mathbb{R}[\varepsilon]/(\varepsilon^{n+1}); hyper-multidual (HMD) quaternions are dual quaternions with coefficients in this algebra. Temporal HMD transforms encode a pose and its derivatives, whereas a generic HMD curve need not be the temporal jet of its pose projection; we call this requirement holonomicity. We show that a temporal transform and its relative descriptor are unitary and derive recursive coefficient constraints, together with a local realizability converse in an admissible logarithm chart. We then extend screw linear interpolation (ScLERP) algebraically to unit HMD quaternions. Although it matches complete endpoint transforms, direct HMD--ScLERP is generically non-holonomic for arbitrary endpoint jets. We give a coefficient criterion and explicit endpoint and first-order interior contact defects. A holonomic alternative is obtained by mapping endpoint transforms to logarithmic dual-quaternion coordinates, applying the degree-(2n+1)(2n+1) Hermite polynomial that matches derivatives through order nn, and lifting by the exponential. HMD arithmetic also recovers higher-order rigid-motion acceleration fields without explicit differentiation of dexp\mathrm{dexp}. Rotation and full SE(3)\mathrm{SE}(3) tests through second order, with an additional third-order polynomial check, reproduce the stated defects and endpoint jets.

Daniel Condurache

Technical University of Iasi, D. Mangeron Street no. 59, 700050 Iasi, Romania · Center for Industrial Robots Simulation and Testing, Technical University of Cluj-Napoca, Memorandumului 14, 400114 Cluj-Napoca, Romania

Lie Group MethodsRepresentation Geometry
AskPDF
Oct 9, 2026

Cubic Spline Interpolation for Smooth Trajectory Generation in Autonomous Motion Systems

We presented a program that automates the construction and solution of the cubic spline interpolation problem for a sequence of position-time waypoints, producing trajectories that are continuous in position, velocity, and acceleration and therefore avoid the infinite-acceleration problem inherent to piecewise-linear waypoint following. By assembling the interpolation and continuity conditions into a single linear system and solving directly for the spline coefficients, the program removes the need to manually derive per-segment equations for each new waypoint set. We illustrated the approach on free-form curve drawing and signature digitization, two tasks that use the same underlying machinery required for smooth robot and drone trajectory generation. Future work includes extending the implementation to multi-axis trajectories with explicit velocity and acceleration bounds, adopting the more computationally efficient tridiagonal solution method for larger waypoint sets, and evaluating the approach quantitatively against alternative trajectory generation methods.

Aarush Mandayam, Sriram Aditya Gundla, Elliott Tammaro

110x Accel · Chestnut Hill College, Philadelphia, Pennsylvania

Trajectory GenerationRobot Motion Generation
AskPDF
Apr 11, 2026

A Coordinate-Invariant Local Representation of Motion and Force Trajectories for Identification and Generalization Across Coordinate Systems

Identifying the trajectories of rigid bodies and of interaction forces is essential for a wide range of tasks in robotics, biomechanics, and related domains. These tasks include trajectory segmentation, recognition, and prediction. For these tasks, a key challenge lies in achieving consistent results when the trajectory is expressed in different coordinate systems. A way to address this challenge is to utilize trajectory models that can generalize across coordinate systems. The focus of this paper is on such trajectory models obtained by transforming the trajectory into a coordinate-invariant representation. However, coordinate-invariant representations often suffer from sensitivity to measurement noise and the manifestation of singularities in the representation, where the representation is not uniquely defined. This paper aims to address this limitation by introducing the novel Dual-Upper-Triangular Invariant Representation (DUTIR), with improved robustness to singularities, along with its computational algorithm. The proposed representation is formulated at a level of abstraction that makes it applicable to both rigid-body trajectories and interaction-force trajectories, hence making it a versatile tool for robotics, biomechanics, and related domains.

Arno Verduyn, Erwin Aertbeliën, Maxim Vochten +1

Department of Mechanical Engineering and Flanders Make at KU Leuven, 3001 Leuven, Belgium. · Department of Mechanics, Royal Military Academy, 1000 Brussels, Belgium.

Motion Representation Learning
AskPDF