Search

Exploring the design and control strategies of soft robots through simulation is highly attractive due to its cost-effectiveness. Although many existing models (e.g., finite element analysis) are effective for simulating soft robotic dynamics, there remains a need for a general and efficient numerical simulation approach in the soft robotics community. In this paper, we develop a discrete differential geometry-based numerical framework to achieve the model-based inverse design of a novel snap-actuated jumping robot. It is found that the dynamic process of a snapping beam can be either symmetric or asymmetric, such that the trajectory of the jumping robot can be tunable (e.g., horizontal or vertical). By employing this novel mechanism of the bistable beam as the robotic actuator, we next propose a physics-data hybrid inverse design strategy for the snap-jump robot with a broad spectrum of jumping capabilities. We first use the physical engine to study the influences of the robot's design parameters on the jumping capabilities, then generate extensive simulation data to formulate a data-driven inverse design solution. The inverse design solution can rapidly explore the combination of design parameters for achieving a target jump, which provides valuable guidance for the fabrication and control of the jumping robot. The proposed methodology paves the way for exploring the design and control insights of soft robots with the help of simulations., Comment: 8 pages, 6 figures

Serpentine structures, composed of straight and circular strips, have garnered significant attention as potential designs for flexible electronics due to their remarkable stretchability. When subjected to stretching, these serpentine strips buckle out of plane, and previous studies have identified two distinct buckling modes whose order of appearance may interchange in serpentine structures with a single cell. In this study, we employ anisotropic rod theory to model serpentine strips as a multi-segment boundary value problem (BVP), with continuity conditions enforced at the interface between the straight and curved strips. We solve the BVP using methods of continuation, and our results reveal that: 1) the exchange of the two buckling modes in a single-cell serpentine strip is induced by a double-eigenvalue and associated secondary bifurcations, which also alter the stability of the two buckling modes; 2) a variety of stable states with reversible symmetry can be manually obtained in tabletop models and are found to be disconnected from the planar branch in numerical continuation. Furthermore, we demonstrate that modulating the strip thickness across different cells leads to the initiation of buckling in the thinnest section, thereby allowing for the tuning of buckling modes in serpentine strips. In structures with two cells, the sequence of the two buckling modes can also be controlled by designing serpentine strips with nonuniform height. This work could enhance the mechanical design of serpentine-interconnect-based flexible structures and could have applications in multistable actuators and mechanical memory devices.

Slender structures, such as rods, often exhibit large nonlinear geometrical deformations even under moderate external forces (e.g., gravity). This characteristic results in a rich variety of morphological changes, making them appealing for engineering design and applications, such as soft robots, submarine cables, decorative knots, and more. Prior studies have demonstrated that the natural shape of a rod significantly influences its deformed geometry. Consequently, the natural shape of the rod should be considered when manufacturing and designing rod-like structures. Here, we focus on an inverse problem: can we determine the natural shape of a suspended 2D planar rod so that it deforms into a desired target shape? We begin by formulating a theoretical framework based on the statics of planar rod equilibrium that can compute the natural shape of a planar rod given its target shape. Furthermore, we analyze the impact of uncertainties (e.g., noise in the data) on the accuracy of the theoretical framework. The results reveal the shortcomings of the theoretical framework in handling uncertainties in the inverse problem, a fact often overlooked in previous works. To mitigate the influence of the uncertainties, we combine the statics of the planar rod with the adjoint method for parameter sensitivity analysis, constructing a learning framework that can efficiently explore the natural shape of the designed rod with enhanced robustness. This framework is validated numerically for its accuracy and robustness, offering valuable insights into the inverse design of soft structures for various applications, including soft robotics and animation of morphing structures., Comment: 21 pages, 9 figures

Soft robots utilizing inflatable dielectric membranes can realize intricate functionalities through the application of non-mechanical fields. However, given the current limitations in simulations, including low computational efficiency and difficulty in dealing with complex external interactions, the design and control of such soft robots often require trial and error. Thus, a novel one-dimensional (1D) discrete differential geometry (DDG)-based numerical model is developed for analyzing the highly nonlinear mechanics in axisymmetric inflatable dielectric membranes. The model captures the intricate dynamics of these membranes under both inflationary pressure and electrical stimulation. Comprehensive validations using hyperelastic benchmarks demonstrate the model's accuracy and reliability. Additionally, the focus on the electro-mechanical coupling elucidates critical insights into the membrane's behavior under varying internal pressures and electrical loads. The research further translates these findings into innovative soft robotic applications, including a spherical soft actuator, a soft circular fluid pump, and a soft toroidal gripper, where the snap-through of electroelastic membrane plays a crucial role. Our analyses reveal that the functional ranges of soft robots are amplified by the snap-through of an electroelastic membrane upon electrical stimuli. This study underscores the potential of DDG-based simulations to advance the understanding of the nonlinear mechanics of electroelastic membranes and guide the design of electroelastic actuators in soft robotics applications., Comment: 27 pages, 8 figures

In this paper, we propose a novel one-dimensional (1D) discrete differential geometry (DDG)-based numerical method for geometrically nonlinear mechanics analysis (e.g., buckling and snapping) of axisymmetric shell structures. Our numerical model leverages differential geometry principles to accurately capture the complex nonlinear deformation patterns exhibited by axisymmetric shells. By discretizing the axisymmetric shell into interconnected 1D elements along the meridional direction, the in-plane stretching and out-of-bending potentials are formulated based on the geometric principles of 1D nodes and edges under the Kirchhoff-Love hypothesis, and elastic force vector and associated Hession matrix required by equations of motion are later derived based on symbolic calculation. Through extensive validation with available theoretical solutions and finite element method (FEM) simulations in literature, our model demonstrates high accuracy in predicting the nonlinear behavior of axisymmetric shells. Importantly, compared to the classical theoretical model and three-dimensional (3D) FEM simulation, our model is highly computationally efficient, making it suitable for large-scale real-time simulations of nonlinear problems of shell structures such as instability and snap-through phenomena. Moreover, our framework can easily incorporate complex loading conditions, e.g., boundary nonlinear contact and multi-physics actuation, which play an essential role in the use of engineering applications, such as soft robots and flexible devices. This study demonstrates that the simplicity and effectiveness of the 1D discrete differential geometry-based approach render it a powerful tool for engineers and researchers interested in nonlinear mechanics analysis of axisymmetric shells, with potential applications in various engineering fields., Comment: 36 pages, 11 figures

In this paper, we systematically study the dynamic snap-through behavior of a pre-deformed elastic ribbon by combining theoretical analysis, discrete numerical simulations, and experiments. By rotating one of its clamped ends with controlled angular speed, we observe two snap-through transition paths among the multiple stable configurations of a ribbon in three-dimensional (3D) space, which is different from the classical snap-through of a two-dimensional (2D) bistable beam. Our theoretical model for the static bifurcation analysis is derived based on the Kirchhoff equations, and dynamical numerical simulations are conducted using the Discrete Elastic Rods (DER) algorithm. The planar beam model is also employed for the asymptotic analysis of dynamic snap-through behaviors. The results show that, since the snap-through processes of both planar beams and 3D ribbons are governed by the saddle-node bifurcation, the same scaling law for the delay applies. We further demonstrate that, in elastic ribbons, by controlling the rotating velocity at the end, distinct snap-through pathways can be realized by selectively skipping specific modes, moreover, particular final modes can be strategically achieved. Through a parametric study using numerical simulations, we construct general phase diagrams for both mode skipping and selection of snapping ribbons. The work serves as a benchmark for future investigations on dynamic snap-through of thin elastic structures and provides guidelines for the novel design of intelligent mechanical systems., Comment: 22 pages, 16 figures

A snap-through bifurcation occurs when a bistable structure loses one of its stable states and moves rapidly to the remaining state. For example, a buckled arch with symmetrically clamped ends can snap between an inverted and a natural state as the ends are released. A standard linear stability analysis suggests that the arch becomes unstable to asymmetric perturbations. Surprisingly, our experiments show that this is not always the case: symmetric transitions are also observed. Using experiments, numerics, and a toy model, we show that the symmetry of the transition depends on the rate at which the ends are released with sufficiently fast loading leading to symmetric snap-through. Our toy model reveals that this behavior is caused by a region of the system's state space in which any initial asymmetry is amplified. The system may not enter this region when loaded fast (hence remaining symmetric) but will traverse it for some interval of time when loaded slowly, causing a transient amplification of asymmetry. Our toy model suggests that this behavior is not unique to snapping arches, but rather can be observed in dynamical systems where both a saddle-node and a pitchfork bifurcation occur in close proximity., Comment: 6 pages main text + 19 pages of supplementary material

Deformable linear objects (DLOs), such as rods, cables, and ropes, play important roles in daily life. However, manipulation of DLOs is challenging as large geometrically nonlinear deformations may occur during the manipulation process. This problem is made even more difficult as the different deformation modes (e.g., stretching, bending, and twisting) may result in elastic instabilities during manipulation. In this paper, we formulate a physics-guided data-driven method to solve a challenging manipulation task -- accurately deploying a DLO (an elastic rod) onto a rigid substrate along various prescribed patterns. Our framework combines machine learning, scaling analysis, and physical simulations to develop a physics-based neural controller for deployment. We explore the complex interplay between the gravitational and elastic energies of the manipulated DLO and obtain a control method for DLO deployment that is robust against friction and material properties. Out of the numerous geometrical and material properties of the rod and substrate, we show that only three non-dimensional parameters are needed to describe the deployment process with physical analysis. Therefore, the essence of the controlling law for the manipulation task can be constructed with a low-dimensional model, drastically increasing the computation speed. The effectiveness of our optimal control scheme is shown through a comprehensive robotic case study comparing against a heuristic control method for deploying rods for a wide variety of patterns. In addition to this, we also showcase the practicality of our control scheme by having a robot accomplish challenging high-level tasks such as mimicking human handwriting, cable placement, and tying knots., Comment: YouTube video: https://youtu.be/OSD6dhOgyMA?feature=shared

The inflation of hyperelastic thin shells is an important and highly nonlinear problem that arises in multiple engineering applications involving severe kinematic and constitutive nonlinearities in addition to various instabilities. We present an isogeometric approach to compute the inflation of hyperelastic thin shells, following the Kirchhoff-Love hypothesis and associated large deformation. Both the geometry and the deformation field are discretized using Catmull-Clark subdivision bases which provide the C1-continuous finite element framework required for the Kirchhoff-Love shell formulation. To follow the complex nonlinear response of hyperelastic thin shells, the inflation is simulated incrementally, and each incremental step is solved via the Newton-Raphson method enriched with arc-length control. Eigenvalue analysis of the linear system after each incremental step allows for inducing bifurcation to a lower energy mode in case stability of the equilibrium is lost. The proposed method is first validated using benchmarks, and then applied to engineering applications, where we demonstrate the ability to simulate large deformation and associated complex instabilities., Comment: 25 pages, 12 figures

Mechanical interactions between rigid rings and flexible cables find broad application in both daily life (hanging clothes) and engineering systems (closing a tether-net). A reduced-order method for the dynamic analysis of sliding rings on a deformable one-dimensional (1D) rod-like object is proposed. In contrast to the conventional approach of discretizing joint rings into multiple nodes and edges for contact detection and numerical simulation, a single point is used to reduce the order of the model. To ensure that the sliding ring and flexible rod do not deviate from their desired positions, a new barrier function is formulated using the incremental potential theory. Subsequently, the interaction between tangent frictional forces is obtained through a delayed dissipative approach. The proposed barrier functional and the associated frictional functional are C2 continuous, hence the nonlinear elastodynamic system can be solved variationally by an implicit time-stepping scheme. The numerical framework is initially applied to simple examples where the analytical solutions are available for validation. Then, multiple complex practical engineering examples are considered to showcase the effectiveness of the proposed method. The simplified ring-to-rod interaction model has the capacity to enhance the realism of visual effects in image animations, while simultaneously facilitating the optimization of designs for space debris removal systems., Comment: We have reorganized the structures of this manuscript

In this paper, a flexible tether-net system is applied to capture the space debris and a numerical framework is established to explore its nonlinear dynamic behaviors, which comprises four principal phases: folding, spreading, contacting, and closing. Based on the discretization of the whole structure into multiple nodes and connected edges, elastic force vectors and associated Jacobian matrix are derived analytically to solve a series of equations of motion. With a fully implicit method applied to analyze the nonlinear dynamics of a slender rod network, the involved mechanical responses are investigated numerically accounting for the interactions. Contact between the deformable net and a rigid body is handled implicitly through a cost-effective modified mass algorithm while the catenary theory is utilized to guide the folding process (from planar configuration to origami-like pattern). The dragging and spreading actions for the folded hexagon net could be realized by shooting six corner mass toward a specific direction; next, the six corners would be controlled to move along a prescribed path producing a closing gesture, when touch between the flying net and the target body is detected, so that for the space debris could be captured and removed successfully. We think the established discrete model could provide a novel insight in the design of active debris removal (ADR) techniques, and promote further development of the model-based control of tether tugging systems., Comment: 19 pages, 18 figures

In this paper, we analyze the inverse dynamics and control of a bacteria-inspired uniflagellar robot in a fluid medium at low Reynolds number. Inspired by the mechanism behind the locomotion of flagellated bacteria, we consider a robot comprised of a flagellum -- a flexible helical filament -- attached to a spherical head. The flagellum rotates about the head at a controlled angular velocity and generates a propulsive force that moves the robot forward. When the angular velocity exceeds a threshold value, the hydrodynamic force exerted by the fluid can cause the soft flagellum to buckle, characterized by a dramatic change in shape. In this computational study, a fluid-structure interaction model that combines Discrete Elastic Rods (DER) algorithm with Lighthill's Slender Body Theory (LSBT) is employed to simulate the locomotion and deformation of the robot. We demonstrate that the robot can follow a prescribed path in three dimensional space by exploiting buckling of the flagellum. The control scheme involves only a single (binary) scalar input -- the angular velocity of the flagellum. By triggering the buckling instability at the right moment, the robot can follow an arbitrary path in three dimensional space. We also show that the complexity of the dynamics of the helical filament can be captured using a deep neural network, from which we identify the input-output functional relationship between the control inputs and the trajectory of the robot. Furthermore, our study underscores the potential role of buckling in the locomotion of natural bacteria.

Solving the software dependency issue under the HPC environment has always been a difficult task for both computing system administrators and application scientists. This work would like to tackle the issue by introducing the modern container technology, the Docker, to be specific. By integrating the auto-scaling feature of service discovery with the light-weight virtualization tool, the Docker, the construction of a virtual cluster on top of physical cluster hardware is attempted. Thus, through the isolation of computing environment, a remedy of software dependency of HPC environment is possible., Comment: PRAGMA-ICDS 15