Snap buckling of bistable beams under combined mechanical and magnetic loading.

We investigate the mechanics of bistable, hard-magnetic, elastic beams, combining experiments, finite-element modelling (FEM) and a reduced-order theory. The beam is made of a hard magneto-rheological elastomer, comprising two segments with antiparallel magnetization along the centreline, and is set into a bistable curved configuration by imposing an end-to-end shortening. Reversible snapping is possible between these two stable states. First, we experimentally characterize the critical field strength for the onset of snapping, at different levels of end-to-end shortening. Second, we perform three-dimensional FEM simulations using the Riks method to analyse high-order deformation modes during snapping. Third, we develop a reduced-order centreline-based beam theory to rationalize the observed magneto-elastic response. The theory and simulations are validated against experiments, with an excellent quantitative agreement. Finally, we consider the case of combined magnetic loading and poking force, examining how the applied field affects the bistability and quantifying the maximum load-bearing capacity. Our work provides a set of predictive tools for the rational design of one-dimensional, bistable, magneto-elastic structural elements. This article is part of the theme issue 'Probing and dynamics of shock sensitive shells'.

Erratum: Twist-Induced Snapping in a Bent Elastic Rod and Ribbon [Phys. Rev. Lett. 122, 114301 (2019)].

This corrects the article DOI: 10.1103/PhysRevLett.122.114301.

Wetting dynamics of viscoelastic solid films.

We study the wetting phenomena of a soft viscoelastic solid film on a smooth and flat substrate. A poly-dimethylsiloxane (PDMS) rubber film is suspended from a stage at both ends, and the wetting behavior of the film against a glass substrate is observed while lowering the stage at a constant velocity. We find that the dynamics of the rubber-glass-air contact lines vary with the lowering velocity of the stage. When the stage velocity is sufficiently low, the film wets the substrate smoothly and the contact lines are straight throughout. Consequently, the contact line velocity is proportional to the lowering velocity. As the stage velocity is increased, the contact line velocity reaches a maximum at the critical stage velocity and then subsequently decreases. The contact lines are wavy and sensitive to the defects above the critical velocity, resulting in the trapping of air bubbles at the interface. We reproduce the wetting behavior using a simple numerical model, assuming an upper limit for the contact line velocity. The wetting behavior observed in our experiments is attributed to the transition in the in-plane stress state from tensile to compressive along the film, leading to buckling of the film above the critical stage velocity. Our results suggest the existence and importance of the maximum wetting velocity for viscoelastic solids.

Mechanics of two filaments in tight orthogonal contact.

Networks of flexible filaments often involve regions of tight contact. Predictively understanding the equilibrium configurations of these systems is challenging due to intricate couplings between topology, geometry, large nonlinear deformations, and friction. Here, we perform an in-depth study of a simple, yet canonical, problem that captures the essence of contact between filaments. In the orthogonal clasp, two filaments are brought into contact, with each centerline lying in one of a pair of orthogonal planes. Our data from X-ray tomography (µCT) and mechanical testing experiments are in excellent agreement with finite element method (FEM) simulations. Despite the apparent simplicity of the physical system, the data exhibit strikingly unintuitive behavior, even when the contact is frictionless. Specifically, we observe a curvilinear diamond-shaped ridge in the contact-pressure field between the two filaments, sometimes with an inner gap. When a relative displacement is imposed between the filaments, friction is activated, and a highly asymmetric pressure field develops. These findings contrast to the classic capstan analysis of a single filament wrapped around a rigid body. Both the µCT and FEM data indicate that the cross-sections of the filaments can deform significantly. Nonetheless, an idealized geometrical theory assuming undeformable tube cross-sections and neglecting elasticity rationalizes our observations qualitatively and highlights the central role of the small, but nonzero, tube radius of the filaments. We believe that our orthogonal clasp analysis provides a building block for future modeling efforts in frictional contact mechanics of more complex filamentary structures.

Loopy Lévy flights enhance tracer diffusion in active suspensions.

Brownian motion is widely used as a model of diffusion in equilibrium media throughout the physical, chemical and biological sciences. However, many real-world systems are intrinsically out of equilibrium owing to energy-dissipating active processes underlying their mechanical and dynamical features1. The diffusion process followed by a passive tracer in prototypical active media, such as suspensions of active colloids or swimming microorganisms2, differs considerably from Brownian motion, as revealed by a greatly enhanced diffusion coefficient3-10 and non-Gaussian statistics of the tracer displacements6,9,10. Although these characteristic features have been extensively observed experimentally, there is so far no comprehensive theory explaining how they emerge from the microscopic dynamics of the system. Here we develop a theoretical framework to model the hydrodynamic interactions between the tracer and the active swimmers, which shows that the tracer follows a non-Markovian coloured Poisson process that accounts for all empirical observations. The theory predicts a long-lived Lévy flight regime11 of the loopy tracer motion with a non-monotonic crossover between two different power-law exponents. The duration of this regime can be tuned by the swimmer density, suggesting that the optimal foraging strategy of swimming microorganisms might depend crucially on their density in order to exploit the Lévy flights of nutrients12. Our framework can be applied to address important theoretical questions, such as the thermodynamics of active systems13, and practical ones, such as the interaction of swimming microorganisms with nutrients and other small particles14 (for example, degraded plastic) and the design of artificial nanoscale machines15.

A mathematical model explores the contributions of bending and stretching forces to shoot gravitropism in Arabidopsis.

Plant shoot gravitropism is a complex phenomenon resulting from gravity sensing, curvature sensing (proprioception), the ability to uphold self-weight and growth. Although recent data analysis and modelling have revealed the detailed morphology of shoot bending, the relative contribution of bending force (derived from the gravi-proprioceptive response) and stretching force (derived from shoot axial growth) behind gravitropism remains poorly understood. To address this gap, we combined morphological data with a theoretical model to analyze shoot bending in wild-type and lazy1-like 1 mutant Arabidopsis thaliana. Using data from actual bending events, we searched for model parameters that minimized discrepancies between the data and mathematical model. The resulting model suggests that both the bending force and the stretching force differ significantly between the wild type and mutant. We discuss the implications of the mechanical forces associated with differential cell growth and present a plausible mechanical explanation of shoot gravitropism.

Twist-Induced Snapping in a Bent Elastic Rod and Ribbon.

Snapping of a slender structure is utilized in a wide range of natural and manmade systems, mostly to achieve rapid movement without relying on musclelike elements. Although several mechanisms for elastic energy storage and rapid release have been studied in detail, a general understanding of the approach to design such a kinetic system is a key challenge in mechanics. Here we study a twist-driven buckling and fast flip dynamics of a geometrically constrained ribbon by combining experiments, numerical simulations, and an analytical theory. We identify two distinct types of shape transitions: A narrow ribbon snaps, and a wide ribbon forms a pair of localized helices. We construct a phase diagram and explain the origin of the boundary, which is determined largely by the geometry. We quantify the effects of gravity and clarify the timescale dictating the rapid flipping. Our study reveals the unique role of geometric twist-bend coupling in the fast dynamics of a thin constrained structure, which has implications for a wide range of biophysical and applied physical problems.

Snap-buckling in asymmetrically constrained elastic strips.

When a flat elastic strip is compressed along its axis, it is bent in one of two possible directions via spontaneous symmetry breaking, forming a cylindrical arc. This is a phenomenon well known as Euler buckling. When this cylindrical section is pushed in the other direction, the bending direction can suddenly reverse. This instability is called "snap-through buckling" and is one of the elementary shape transitions in a prestressed thin structure. Combining experiments and theory, we study snap-buckling of an elastic strip with one end hinged and the other end clamped. These asymmetric boundary constraints break the intrinsic symmetry of the strip, generating mechanical behaviors, including largely hysteretic but reproducible force responses and switchlike discontinuous shape changes. We establish the set of exact analytical solutions to fully explain all our major experimental and numerical findings. Asymmetric boundary conditions arise naturally in diverse situations when a thin object is in contact with a solid surface at one end. The introduction of asymmetry through boundary conditions yields new insight into complex and programmable functionalities in material and industrial design.

Slip Morphology of Elastic Strips on Frictional Rigid Substrates.

The morphology of an elastic strip subject to vertical compressive stress on a frictional rigid substrate is investigated by a combination of theory and experiment. We find a rich variety of morphologies, which-when the bending elasticity dominates over the effect of gravity-are classified into three distinct types of states: pinned, partially slipped, and completely slipped, depending on the magnitude of the vertical strain and the coefficient of static friction. We develop a theory of elastica under mixed clamped-hinged boundary conditions combined with the Coulomb-Amontons friction law and find excellent quantitative agreement with simulations and controlled physical experiments. We also discuss the effect of gravity in order to bridge the difference in the qualitative behaviors of stiff strips and flexible strings or ropes. Our study thus complements recent work on elastic rope coiling and takes a significant step towards establishing a unified understanding of how a thin elastic object interacts vertically with a solid surface.

Granular rotor as a probe for a nonequilibrium bath.

This study numerically and analytically investigates the dynamics of a rotor under viscous or dry friction as a nonequilibrium probe of a granular gas. In order to demonstrate the role of the rotor as a probe for a nonequilibrium bath, the molecular dynamics (MD) simulation of the rotor is performed under viscous or dry friction surrounded by a steady granular gas under gravity. A one-to-one map between the velocity distribution function (VDF) of the granular gas and the angular distribution function for the rotor is theoretically derived. The MD simulation demonstrates that the one-to-one map accurately infers the local VDF of the granular gas from the angular VDF of the rotor, and vice versa.

Minimal model of stochastic athermal systems: origin of non-Gaussian noise.

For a wide class of stochastic athermal systems, we derive Langevin-like equations driven by non-Gaussian noise, starting from master equations and developing a new asymptotic expansion. We found an explicit condition whereby the non-Gaussian properties of the athermal noise become dominant for tracer particles associated with both thermal and athermal environments. Furthermore, we derive an inverse formula to infer microscopic properties of the athermal bath from the statistics of the tracer particle. We apply our formulation to a granular motor under viscous friction and analytically obtain the angular velocity distribution function. Our theory demonstrates that the non-Gaussian Langevin equation is the minimal model of athermal systems.

Roles of dry friction in the fluctuating motion of an adiabatic piston.

The motion of an adiabatic piston under dry friction is investigated to clarify the roles of dry friction in nonequilibrium steady states. We clarify that dry friction can reverse the direction of the piston motion and causes a discontinuity or a cusplike singularity for velocity distribution functions of the piston. We also show that the heat fluctuation relation is modified under dry friction.

Simulation of granular jets: is granular flow really a perfect fluid?

We perform three-dimensional simulations of the impact of a granular jet for both frictional and frictionless grains. Small shear stress observed in the experiment [X. Cheng et al., Phys. Rev. Lett. 99, 188001 (2007)] is reproduced through our simulation. However, the fluid state after the impact is far from a perfect fluid, and thus the similarity between granular jets and quark gluon plasma is superficial because the observed viscosity is finite and its value is consistent with the prediction of the kinetic theory.