Elasticity and stability of shape-shifting structures

As we enter the age of designer matter — where objects can morph and change shape on command — what tools do we need to create shape-shifting structures? At the heart of an elastic deformation is the combination of dilation and distortion or stretching and bending. The competition between the latter can cause elastic instabilities, and over the last fifteen years, these instabilities have provided a multitude of ways to prescribe and control shape change. Buckling, wrinkling, folding, creasing, and snapping have become mechanisms that when harmoniously combined enable mechanical metamaterials, self-folding origami, ultralight and ultrathin kirigami, and structures that appear to grow from one shape to another. In this review, I aim to connect the fundamentals of elastic instabilities to the advanced functionality currently found within mechanical metamaterials.     

The flip side of buckling

Buckling of slender structures under compressive loading is a failure of infinitesimal stability due to a confluence of two factors: the energy density non-convexity and the smallness of Korn’s constant. The problem has been well understood only for bodies with simple geometries when the slenderness parameter is well defined. In this paper, we present the first rigorous analysis of buckling for bodies with complex geometry. By limiting our analysis to the “near-flip” instability, we address the universal features of the buckling phenomenon that depend on neither the shape of the domain nor the degree of constitutive nonlinearity of the elastic material.      

Buckling and post-buckling of long pressurized elastic thin-walled tubes under in-plane bending

The present paper focuses on the structural stability of long uniformly pressurized thin elastic tubular shells subjected to in-plane bending. Using a special-purpose non-linear finite element technique, bifurcation on the pre-buckling ovalization equilibrium path is detected, and the post-buckling path is traced. Furthermore, the influence of pressure (internal and/or external) as well as the effects of radius-to-thickness ratio, initial curvature and initial ovality on the bifurcation moment, curvature and the corresponding wavelength, are examined. The local character of buckling in the circumferential direction is also demonstrated, especially for thin-walled tubes. This observation motivates the development of a simplified analytical formulation for tube bifurcation, which considers the presence of pressure, initial curvature and ovality, and results in closed-form expressions of very good accuracy, for tubes with relatively small initial curvature. Finally, aspects of tube bifurcation are illustrated using a simple mechanical model, which considers the ovalized pre-buckling state and the effects of pressure.     

Experimental investigation of a buckled beam under high-frequency excitation

In this paper, we investigate theoretically and experimentally dynamics of a buckled beam under high-frequency excitation. It is theoretically predicted from linear analysis that the high-frequency excitation shifts the pitchfork bifurcation point and increases the buckling force. The shifting amount increases as the excitation amplitude or frequency increases. Namely, under the compressive force exceeding the buckling one, high-frequency excitation can stabilize the beam to the straight position. Some experiments are performed to investigate effects of the high-frequency excitation on the buckled beam. The dependency of the buckling force on the amounts of excitation amplitude and frequency is compared with theoretical results. The transient state is observed in which the beam is recovered from the buckled position to the straight position due to the excitation. Furthermore, the bifurcation diagrams are measured in the cases with and without high-frequency excitation. It is experimentally clarified that the high-frequency excitation changes the nonlinear property of the bifurcation from supercritical pitchfork bifurcation to subcritical pitchfork bifurcation and then the stable steady state of the beam exhibits hysteresis as the compressive force is reversed. This work was partially supported by the Japanese Ministry of Education, Culture, Sports, Science, and Technology, under Grants-in-Aid for Scientific Research 16560377.     

Flexural-torsional bifurcations of a cantilever beam under potential and circulatory forces I: Non-linear model and stability analysis

The stability of a cantilever elastic beam with rectangular cross-section under the action of a follower tangential force and a bending conservative couple at the free end is analyzed. The beam is herein modeled as a non-linear Cosserat rod model. Non-linear, partial integro-differential equations of motion are derived expanded up to cubic terms in the transversal displacement and torsional angle of the beam. The linear stability of the trivial equilibrium is studied, revealing the existence of buckling, flutter and double-zero critical points. Interaction between conservative and non-conservative loads with respect to the stability problem is discussed. The critical spectral properties are derived and the corresponding critical eigenspace is evaluated.     

Local buckling of wide-flange thin-walled anisotropic composite beams

The present contribution deals with the onset of local buckling of compressively loaded thin-walled beams with open I, C, Z, T and L-cross-sections made of laminated composite materials. The method employs a discrete plate analysis approach in the course of which each structural subelement of interest—which presently is the flange—of the thin-walled cross-section is considered as a separate composite plate with elastic rotational restraints at those edges where an adjacent substructural element is located. While in many investigations the lamination schemes of webs and flanges are considered to be purely orthotropic, in the present paper the laminate layups are allowed to be of an arbitrary non-orthotropic nature, which also allows for the analysis of laminates with inherent bending–torsion coupling. The analysis of the buckling loads of the flanges of thin-walled composite beams is performed using the Ritz-method for which some especially adjusted displacement shape functions are employed. For the case of pure orthotropy, a novel closed-form solution is described. The accuracy of the employed approaches is established by comparison with accompanying finite element simulations of thin-walled composite beams. It is revealed that the presented methodology is highly efficient in terms of computational effort and yet performs with satisfying accuracy, which makes it very attractive for actual practical applications whenever the local stability behaviour of wide-flange thin-walled composite beams is to be considered.     

Effect of Coulomb Damping on Buckling of a Two-Rod System

This paper describes a significant influence of a slight Coulomb damping on buckling, using a simple two rods system. Coulomb damping produces equilibrium regions around the well-known stable and unstable steady states under the pitchfork bifurcation which occurs in the case without Coulomb damping. Also, the stability of the states in the equilibrium regions is examined by using the phase portrait. As a consequence, due to the slight Coulomb damping, it is theoretically clarified that the states in the equilibrium regions are locally stable, even in the neighborhood of the unstable steady states under the pitchfork bifurcation in the case without Coulomb damping, i.e., even in the neighborhood of the unstable trivial steady states in the postbuckling and the unstable nontrivial steady states under the subcritical pitchfork bifurcation. Furthermore, the experimental results are in qualitative agreement with the theoretically predicted phenomena.     

Elastic stability of all edges simply supported,stepped and stiffened rectangular plate under Biaxial loading

In this paper using finite difference method the lower bound buckling load for simply supported (a) stepped and stiffened rectangular thin plate (b) linear and non-linear variation of thickness (c) uniformly distributed compressive forces in both directions (d) uniformly distributed compressive force in y direction and non-uniform distribution of compressive force in x-direction is discussed. The thin plate is divided into 900 rectangular meshes. The partial derivatives are approximated using central difference formula. Eight hundred and forty one equations are formed and using the program developed and the least eigenvalue is obtained. The buckling coefficients are calculated for different types of stepped and non prismatic plates and the results are presented in tables and graphs for ready use by designers. Buckling factors for some cases are presented in the form of three separate tables and compared with the values obtained by Xiang, Wei and Wang. The results are in close agreement.     

Supercritical and subcritical buckling bifurcations for a compressible hyperelastic slab subjected to compression

We study the buckling bifurcation of a compressible hyperelastic slab under compression with sliding–sliding end conditions. The combined series-asymptotic expansions method is used to derive the simplified model equations. Linear bifurcation analysis yields the critical stress value of buckling, which gives a non-linear correction to the classical Euler buckling formula. The correction is due to the geometrical non-linearities coupled with the material non-linearities. Then through non-linear bifurcation analysis, the approximate analytical solutions for the post-buckling deformations are obtained. The amplitude of buckling is expressed explicitly in terms of the aspect ratio, the incremental dimensionless engineering stress, the mode of buckling and the material constants. Most importantly, we find that both supercritical and subcritical buckling could occur for compressible materials. The bifurcation type depends on the material constants, the geometry of the slab and the mode numbers.     

Finite Element Analysis of Dynamically Loaded Homogeneous Viscous Beams

The paper is concerned with impulsively loaded beams in which the material is treated as homogeneous viscous as an approximation of a rigid-viscoplastic constitutive relation. As opposed to the standard displacement method finite element formulation, where interpolation functions describing the velocity field across elements is given, a mixed formulation is used in which nodal velocities and nodal moments are carried as parameters. At each instant the accelerations (by the Tamuzh principle) and the rates of change of moment (by a virtual velocities formulation) are found, and velocities and moments are integrated forward independently. The properties of the mode solution are also introduced, and the forward integration is carried through only for the difference between the mode solution and the actual solution. This leads to a very efficient scheme for the numerical solution of a cantilever beam problem shown as an illustration.