skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Solving Cosserat Rod Models via Collocation and the Magnus Expansion
Choosing a kinematic model for a continuum robot typically involves making a tradeoff between accuracy and computational complexity. One common modeling approach is to use the Cosserat rod equations, which have been shown to be accurate for many types of continuum robots. This approach, however, still presents significant computational cost, particularly when many Cosserat rods are coupled via kinematic constraints. In this work, we propose a numerical method that combines orthogonal collocation on the local rod curvature and forward integration of the Cosserat rod kinematic equations via the Magnus expansion, allowing the equilibrium shape to be written as a product of matrix exponentials. We provide a bound on the maximum step size to guarantee convergence of the Magnus expansion for the case of Cosserat rods, compare in simulation against other approaches, and demonstrate the tradeoffs between speed and accuracy for the fourth and sixth order Magnus expansions as well as for different numbers of collocation points. Our results show that the proposed method can find accurate solutions to the Cosserat rod equations and can potentially be competitive in computation speed.  more » « less
Award ID(s):
1734461
PAR ID:
10281088
Author(s) / Creator(s):
;
Date Published:
Journal Name:
2020 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS)
Page Range / eLocation ID:
8653 to 8660
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Controlling soft continuum robotic arms is challenging due to their hyper-redundancy and dexterity. In this paper we experimentally demonstrate, for the first time, closed-loop control of the configuration space variables of a soft robotic arm, composed of independently controllable segments, using a Cosserat rod model of the robot and the distributed sensing and actuation capabilities of the segments. Our controller solves the inverse dynamic problem by simulating the Cosserat rod model in MATLAB using a computationally efficient numerical solution scheme, and it applies the computed control output to the actual robot in real time. The position and orientation of the tip of each segment are measured in real time, while the remaining unknown variables that are needed to solve the inverse dynamics are estimated simultaneously in the simulation. We implement the controller on a multi-segment silicone robotic arm with pneumatic actuation, using a motion capture system to measure the segments' positions and orientations. The controller is used to reshape the arm into configurations that are achieved through combinations of bending and extension deformations in 3D space. Although the possible deformations are limited for this robot platform, our study demonstrates the potential for implementing the control approach on a wide range of continuum robots in practice. The resulting tracking performance indicates the effectiveness of the controller and the accuracy of the simulated Cosserat rod model. 
    more » « less
  2. For many problems in science and engineering, it is necessary to describe the collective behavior of a very large number of grains. Complexity inherent in granular materials, whether due the variability of grain interactions or grain-scale morphological factors, requires modeling approaches that are both representative and tractable. In these cases, continuum modeling remains the most feasible approach; however, for such models to be representative, they must properly account for the granular nature of the material. The granular micromechanics approach has been shown to offer a way forward for linking the grain-scale behavior to the collective behavior of millions and billions of grains while keeping within the continuum framework. In this paper, an extended granular micromechanics approach is developed that leads to a micromorphic theory of degree n. This extended form aims at capturing the detailed grain-scale kinematics in disordered (mechanically or morphologically) granular media. To this end, additional continuum kinematic measures are introduced and related to the grain-pair relative motions. The need for enriched descriptions is justified through experimental measurements as well as results from simulations using discrete models. Stresses conjugate to the kinematic measures are then defined and related, through equivalence of deformation energy density, to forces conjugate to the measures of grain-pair relative motions. The kinetic energy density description for a continuum material point is also correspondingly enriched, and a variational approach is used to derive the governing equations of motion. By specifying a particular choice for degree n, abridged models of degrees 2 and 1 are derived, which are shown to further simplify to micro-polar or Cosserat-type and second-gradient models of granular materials. 
    more » « less
  3. Abstract. We consider a nonlinear, moving boundary, fluid-structure interaction problem between a time dependent incompressible, viscous fluid flow, and an elastic structure composed of a cylindrical shell supported by a mesh of elastic rods. The fluid flow is modeled by the time-dependent Navier- Stokes equations in a three-dimensional cylindrical domain, while the lateral wall of the cylinder is modeled by the two-dimensional linearly elastic Koiter shell equations coupled to a one-dimensional system of conservation laws defined on a graph domain, describing a mesh of curved rods. The mesh supported shell allows displacements in all three spatial directions. Two-way coupling based on kinematic and dynamic coupling conditions is assumed between the fluid and composite structure, and between the mesh of curved rods and Koiter shell. Problems of this type arise in many ap- plications, including blood flow through arteries treated with vascular prostheses called stents. We prove the existence of a weak solution to this nonlinear, moving boundary problem by using the time discretization via Lie operator splitting method combined with an Arbitrary Lagrangian-Eulerian approach, and a non-trivial extension of the Aubin-Lions-Simon compactness result to problems on moving domains. 
    more » « less
  4. Numerous soft and continuum robotic manipulators have demonstrated their potential for compliant operation in highly unstructured environments or near people. Despite their recent popularity, modeling of their smooth bending deformation remains a challenge. For soft continuum manipulators, the widespread, constant curvature approach to modeling is inadequate for modeling some deformations that occur in practice, such as combined bending and twisting deformations. In this paper, we extend the classical Cosserat rod approach to model a variable-length, pneumatic soft continuum arm. We model the deformation of a pneumatically driven soft continuum manipulator, and the model is then compared against experimental data collected from a three degree of freedom, pneumatically actuated, soft continuum manipulator. The model shows good agreement in capturing the overall behavior of the bending deformation, with mean Euclidean error at the tip of the robot of 2.48 cm for a 22 cm long robot. In addition, the model shows good numerical stability for simulating long duration computations. 
    more » « less
  5. Abstract The method of kinematic synthesis requires finding the solution set of a system of polynomials. Parameter homotopy continuation is used to solve these systems and requires repeatedly solving systems of linear equations. For kinematic synthesis, the associated linear systems become ill-conditioned, resulting in a marked decrease in the number of solutions found due to path tracking failures. This unavoidable ill-conditioning places a premium on accurate function and matrix evaluations. Traditionally, variables are eliminated to reduce the dimension of the problem. However, this greatly increases the computational cost of evaluating the resulting functions and matrices and introduces numerical instability. We propose avoiding the elimination of variables to reduce required computations, increasing the dimension of the linear systems, but resulting in matrices that are quite sparse. We then solve these systems with sparse solvers to save memory and increase speed. We found that this combination resulted in a speedup of up to 250 × over traditional methods while maintaining the same accuracy. 
    more » « less