Showing 1–43 of 43 results for author: Bigoni, D

Effects of different loading on the bifurcation of annular elastic rods: theory vs. experiments

Authors: Matteo Gaibotti, Davide Bigoni, Arsenio Cutolo, Massimiliano Fraldi, Andrea Piccolroaz

Abstract: The bifurcation problem of a circular Euler-Bernoulli rod subject to a uniform radial force distribution is investigated under three distinct loading conditions: (i.) hydrostatic pressure, (ii.) centrally-directed, and (iii.) dead load. Previous studies on this apparently 'familiar' structural problem have yielded controversial results, necessitating a comprehensive clarification. This study shows… ▽ More The bifurcation problem of a circular Euler-Bernoulli rod subject to a uniform radial force distribution is investigated under three distinct loading conditions: (i.) hydrostatic pressure, (ii.) centrally-directed, and (iii.) dead load. Previous studies on this apparently 'familiar' structural problem have yielded controversial results, necessitating a comprehensive clarification. This study shows that results previously labelled as 'correct' or 'wrong' simply refer to different external constraints, whose presence becomes necessary only for the two latter loads, (ii.) and (iii.). Moreover, the paper presents the first experimental realization of a circular rod subjected to centrally-directed loads. The experimental findings align with the theoretical predictions and show the exploitation of a new type of load acting on a continuous structural element. The feasibility of this load is demonstrated through the use of inextensible cables and opens the way to applications in flexible robotics when cables are used for actuation. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 25 pages, 10 figures

Journal ref: Int. J. NonLinear Mech. 165 (2024) 104820

Effects of prestress in the coating of an elastic disk

Authors: Matteo Gaibotti, Sofia G. Mogilevskaya, Andrea Piccolroaz, Davide Bigoni

Abstract: An elastic disk is coated with an elastic rod, uniformly prestressed with a tensile or compressive axial force. The prestress state is assumed to be induced by three different models of external radial load or by 'shrink-fit' forcing the coating onto the disk. The prestressed coating/disk system, when loaded with an additional and arbitrary incremental external load, experiences incremental displa… ▽ More An elastic disk is coated with an elastic rod, uniformly prestressed with a tensile or compressive axial force. The prestress state is assumed to be induced by three different models of external radial load or by 'shrink-fit' forcing the coating onto the disk. The prestressed coating/disk system, when loaded with an additional and arbitrary incremental external load, experiences incremental displacement, strain, and stress, which are solved via complex potentials. The analysis incorporates models for both perfect and imperfect bonding at the coating/disk interface. The derived solution highlights the significant influence not only of the prestress but also of the method employed to generate it. These two factors lead, in different ways, to a loss or an increase in incremental stiffness for compressive or tensile prestress. The first bifurcation load of the structure (which differs for different prestress generations) is determined in a perturbative way. The results emphasize the importance of modelling the load and may find applications in flexible electronics and robot arms subject to pressure or uniformly-distributed radial forces. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 17 pages, 6 figures

Bifurcations of an elastic disc coated with an elastic inextensible rod

Authors: Matteo Gaibotti, Sonia G. Mogilevskaya, Andrea Piccolroaz, Davide Bigoni

Abstract: An analytical solution is derived for the bifurcations of an elastic disc that is constrained on the boundary with an isoperimetric Cosserat coating. The latter is treated as an elastic circular rod, either perfectly or partially bonded (with a slip interface in the latter case) and is subjected to three different types of uniformly distributed radial loads (including hydrostatic pressure). The pr… ▽ More An analytical solution is derived for the bifurcations of an elastic disc that is constrained on the boundary with an isoperimetric Cosserat coating. The latter is treated as an elastic circular rod, either perfectly or partially bonded (with a slip interface in the latter case) and is subjected to three different types of uniformly distributed radial loads (including hydrostatic pressure). The proposed solution technique employs complex potentials to treat the disc's interior and incremental Lagrangian equations to describe the prestressed elastic rod modelling the coating. The bifurcations of the disc occur with modes characterized by different circumferential wavenumbers, ranging between ovalization and high-order waviness, as a function of the ratio between the elastic stiffness of the disc and the bending stiffness of its coating. The presented results find applications in various fields, such as coated fibres, mechanical rollers, and the growth and morphogenesis of plants and fruits. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 21 pages, 9 figures

Journal ref: Proc. R. Soc. A., 480, 20230491 (2024)

Homogenization of elastic grids containing rigid elements

Authors: Luca Viviani, Davide Bigoni, Andrea Piccolroaz

Abstract: The inclusion of rigid elements into elastic composites may lead to superior mechanical properties for the equivalent elastic continuum, such as, for instance, extreme auxeticity. To allow full exploitation of these properties, a tool for the homogenization of two-dimensional elastic grids containing rigid elements is developed and tested on elaborate geometries, such as, for instance, Chinese lat… ▽ More The inclusion of rigid elements into elastic composites may lead to superior mechanical properties for the equivalent elastic continuum, such as, for instance, extreme auxeticity. To allow full exploitation of these properties, a tool for the homogenization of two-dimensional elastic grids containing rigid elements is developed and tested on elaborate geometries, such as, for instance, Chinese lattices. The rigid elements are assumed to be either jointed with full continuity of displacement or hinged to the elastic rods. It is shown that the two different constraints induce strongly different mechanical characteristics of the equivalent elastic solid. The presented results open the way to the design of architected materials or metamaterials containing both elastic and rigid parts. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Comments: 24 pages, 18 figures

Bimodal buckling governs human fingers luxation

Authors: Massimiliano Fraldi, Stefania Palumbo, Arsenio Cutolo, Angelo Rosario Carotenuto, Davide Bigoni

Abstract: Equilibrium bifurcation in natural systems can sometimes be explained as a route to stress shielding for preventing failure. Although compressive buckling has been known for a long time, its less-intuitive tensile counterpart was only recently discovered and yet never identified in living structures or organisms. Through the analysis of an unprecedented all-in-one paradigm of elastic instability,… ▽ More Equilibrium bifurcation in natural systems can sometimes be explained as a route to stress shielding for preventing failure. Although compressive buckling has been known for a long time, its less-intuitive tensile counterpart was only recently discovered and yet never identified in living structures or organisms. Through the analysis of an unprecedented all-in-one paradigm of elastic instability, it is theoretically and experimentally shown that coexistence of two curvatures in human finger joints is the result of an optimal design by nature that exploits both compressive and tensile buckling for inducing luxation in case of traumas, so realizing a unique mechanism for protecting tissues and preventing more severe damage under extreme loads. Our findings might pave the way to conceive complex architectured and bio-inspired materials, as well as next generation artificial joint prostheses and robotic arms for bio-engineering and healthcare applications. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Journal ref: Proceedings of the National Academy of Sciences 120.44 (2023)

Flutter instability in solids and structures, with a view on biomechanics and metamaterials

Authors: Davide Bigoni, Francesco Dal Corso, Oleg N. Kirillov, Diego Misseroni, Giovanni Noselli, Andrea Piccolroaz

Abstract: The phenomenon of oscillatory instability called \lq flutter' was observed in aeroelasticity and rotor dynamics about a century ago. Driven by a series of applications involving nonconservative elasticity theory at different physical scales, ranging from nanomechanics to the mechanics of large space structures and including biomechanical problems of motility and growth, research on flutter is expe… ▽ More The phenomenon of oscillatory instability called \lq flutter' was observed in aeroelasticity and rotor dynamics about a century ago. Driven by a series of applications involving nonconservative elasticity theory at different physical scales, ranging from nanomechanics to the mechanics of large space structures and including biomechanical problems of motility and growth, research on flutter is experiencing a new renaissance. A review is presented of the most notable applications and recent advances in fundamentals, both theoretical and experimental aspects, of flutter instability and Hopf bifurcation. Open problems, research gaps, and new perspectives for investigations are indicated. △ Less

Submitted 13 January, 2024; originally announced January 2024.

Comments: 29 pages, 9 figures

Journal ref: Proc. R. Soc. A. 479: 20230523 (2023)

Dynamics of elastic lattices with sliding constraints

Authors: L. Cabras, D. Bigoni, A. Piccolroaz

Abstract: This study investigates the impact of sliders -- constraints acting on elastic rods allowing for a transverse displacement jump while maintaining axial and rotational displacement continuity -- on the dynamics of a periodic elastic grid, including the effects of axial preload. The grid is linearly elastic and subject to in-plane incremental deformation, involving normal and shear forces and bendin… ▽ More This study investigates the impact of sliders -- constraints acting on elastic rods allowing for a transverse displacement jump while maintaining axial and rotational displacement continuity -- on the dynamics of a periodic elastic grid, including the effects of axial preload. The grid is linearly elastic and subject to in-plane incremental deformation, involving normal and shear forces and bending moment. The periodicity of the infinite grid permits a Floquet-Bloch wave analysis and a rigorous dynamic homogenization, leading to an equivalent prestressed elastic solid. The investigation is complemented by {\it ad hoc} developed F.E. simulations and perturbations with a pulsating Green's function. Results show that the sliders create band gaps, flat bands, and Dirac cones in the dispersion diagrams and generate macro-instability even for tensile prestress. The latter corresponds to the loss of ellipticity at the parabolic boundary in the equivalent elastic solid and provides a rare example of an almost unexplored form of material instability. Therefore, our results offer design strategies for metamaterials and architected materials showing reversible material instabilities and filtering properties for mechanical signals. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: 21 pages, 20 figures

Stabilization against gravity and self-tuning of an elastic variable-length rod through an oscillating sliding sleeve

Authors: Panagiotis Koutsogiannakis, Diego Misseroni, Davide Bigoni, Francesco Dal Corso

Abstract: An elastic rod, straight in its undeformed state, has a mass attached at one end and a variable length, due to a constraint at the other end by a frictionless sliding sleeve. The constraint is arranged with the sliding direction parallel to a gravity field, in a way that the rod can freely slip inside of the sleeve, when the latter is not moving. In this case, the free fall of the mass continues u… ▽ More An elastic rod, straight in its undeformed state, has a mass attached at one end and a variable length, due to a constraint at the other end by a frictionless sliding sleeve. The constraint is arranged with the sliding direction parallel to a gravity field, in a way that the rod can freely slip inside of the sleeve, when the latter is not moving. In this case, the free fall of the mass continues until the rod is completely injected into the constraint. However, when the sliding sleeve is subject to a harmonic transverse vibration, it is shown that the fall of the mass and the rod injection are hindered by the presence of a configurational force developing at the sliding sleeve and acting oppositely to gravity. During the dynamic motion, such a configurational force is varying in time because it is associated with the variable bending moment at the sleeve entrance. It is (experimentally, analytically, and numerically) demonstrated that, in addition to the states of complete injection or ejection of the elastic rod (for which the mass falls down or is thrown out), a stable sustained oscillation around a finite height can be realized. This 'suspended motion' is the signature of a new attractor, that arises by the constraint oscillation. This behaviour shares similarities with parametric oscillators, as for instance the Kapitza inverted pendulum. However, differently from the classical parametric oscillators, the 'suspended' configuration of the rod violates equilibrium and the stabilization occurs through a transverse mechanical input, instead of a longitudinal one. By varying the sliding sleeve oscillation amplitude and frequency within specific sets of values, the system spontaneously adjusts the sustained motion through a self-tuning of the rod's external length. This self-tuning property opens the way to the design of vibration-based devices with extended frequency range. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Journal ref: Journal of the Mechanics and Physics of Solids, Volume 181, 2023, 105452,ISSN 0022-5096

Phase-field modelling of failure in ceramics with multiscale porosity

Authors: Riccardo Cavuoto, Pietro Lenarda, Anna Tampieri, Davide Bigoni, Marco Paggi

Abstract: Many stiff biological materials exhibiting outstanding compressive strength/weight ratio are characterized by high porosity, spanning different size-scales, typical examples being bone and wood. A successful bio-mimicking of these materials is provided by a recently-obtained apatite, directly produced through a biomorphic transformation of natural wood and thus inheriting its highly hierarchical s… ▽ More Many stiff biological materials exhibiting outstanding compressive strength/weight ratio are characterized by high porosity, spanning different size-scales, typical examples being bone and wood. A successful bio-mimicking of these materials is provided by a recently-obtained apatite, directly produced through a biomorphic transformation of natural wood and thus inheriting its highly hierarchical structure. This unique apatite (but also wood and bone) is characterized by two major distinct populations of differently-sized cylindrical voids, a porosity shown in the present paper to influence failure, both in terms of damage growth and fracture nucleation and propagation. This statement follows from failure analysis, developed through in-silico generation of artificial samples (reproducing the two-scale porosity of the material) and subsequent finite element modelling of damage, implemented with phase-field treatment for fracture growth. It is found that small voids promote damage nucleation and enhance bridging of macro-pores by micro-crack formation, while macro-pores influence the overall material response and drive the propagation of large fractures. Our results explain the important role of multiscale porosity characterizing stiff biological materials and lead to a new design paradigm, by introducing an in-silico tool to implement bio-mimicking in new artificial materials with brittle behaviour, such as carbide or ceramic foams. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 12 pages, 10 figures

Fusion of two stable elastic structures resulting in an unstable system

Authors: Marco Rossi, Andrea Piccolroaz, Davide Bigoni

Abstract: It is shown that a compound elastic structure, which displays a dynamic instability, may be designed as the union (or 'fusion') of two structures which are stable when separately analyzed. The compound elastic structure has two degrees of freedom and is made up of a rigid rod connected with two springs to a smooth support, which evidences a jump in the curvature at the equilibrium configuration. I… ▽ More It is shown that a compound elastic structure, which displays a dynamic instability, may be designed as the union (or 'fusion') of two structures which are stable when separately analyzed. The compound elastic structure has two degrees of freedom and is made up of a rigid rod connected with two springs to a smooth support, which evidences a jump in the curvature at the equilibrium configuration. Instability is proven in a linearized context and is related to the application of a non-conservative load of the follower type, so that the instability disappears under dead loads. In the fully nonlinear range, the instability is also confirmed through numerical simulations. The obtained results may be useful in the design of new mechanical sensors, or devices for energy harvesting, or architected materials. In addition, our findings have conceptual implications on piecewise-linear theories of mechanics such as for instance plasticity or frictional contact. △ Less

Submitted 13 January, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Comments: 39 pages, 19 figures; post-print version; added Journal reference

Journal ref: Journal of the Mechanics and Physics of Solids (2023) 105201

Configurational forces on elastic structures

Authors: Davide Bigoni, Federico Bosi, Francesco Dal Corso, Diego Misseroni

Abstract: The discovery of configurational forces acting on elastic structures and its initial applications are reviewed. Configurational forces are related to the possibility that an elastic structure can change its configuration, thus inducing a variation in the potential energy. This concept has already led to several applications (the elastica arm scale, the dripping of an elastic rod, and the torsional… ▽ More The discovery of configurational forces acting on elastic structures and its initial applications are reviewed. Configurational forces are related to the possibility that an elastic structure can change its configuration, thus inducing a variation in the potential energy. This concept has already led to several applications (the elastica arm scale, the dripping of an elastic rod, and the torsional actuator), has been shown to strongly affect stability, and to be related to limbless locomotion. It is believed that these results will open a new research territory in mechanics. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 12 pages, 8 figures

Journal ref: In 50+ Years of AIMETA (pp. 229-241). Springer, Cham. 2022

Double restabilization and design of force-displacement response of the extensible elastica with movable constraints

Authors: Panagiotis Koutsogiannakis, Davide Bigoni, Francesco Dal Corso

Abstract: A highly deformable rod, modelled as the extensible elastica, is connected to a movable clamp at one end and to a pin sliding along a frictionless curved profile at the other. Bifurcation analysis shows that axial compliance provides a stabilizing effect in compression, but unstabilizing in tension. Moreover, with varying the constraint's curvature at the origin and the axial vs bending rod's stif… ▽ More A highly deformable rod, modelled as the extensible elastica, is connected to a movable clamp at one end and to a pin sliding along a frictionless curved profile at the other. Bifurcation analysis shows that axial compliance provides a stabilizing effect in compression, but unstabilizing in tension. Moreover, with varying the constraint's curvature at the origin and the axial vs bending rod's stiffness, in addition to possible buckling in tension, the structure displays none, two, or even four bifurcation loads, the last two associated only to the first buckling mode in compression. Therefore, the straight configuration may lose and recover stability one or two times, thus evidencing single and double restabilization, a feature never observed before. By means of the closed-form solution for the extensible elastica, the quasi-static behaviour of the structure is analytically described under large rotations and axial strain. The presented solution is exploited, together with an {\it ad hoc} developed optimization algorithm, to design the shape of the constraint's profile necessary to obtain a desired force-displacement curve, so to realize a force-limiter or a mechanical device capable of delivering a complex force response upon application of a continuous displacement in both positive and negative direction. △ Less

Submitted 9 September, 2022; v1 submitted 10 August, 2022; originally announced August 2022.

Comments: 35 pages, 19 figures, Published in European Journal of Mechanics - A/Solids (2022)

Journal ref: European Journal of Mechanics - A/Solids, 2022, 104745

Buckling of thin-walled cylinders from three dimensional nonlinear elasticity

Authors: Roberta Springhetti, Gabriel Rossetto, Davide Bigoni

Abstract: The famous bifurcation analysis performed by Flügge on compressed thin-walled cylinders is based on a series of simplifying assumptions, which allow to obtain the bifurcation landscape, together with explicit expressions for limit behaviours: surface instability, wrinkling, and Euler rod buckling. The most severe assumption introduced by Flügge is the use of an incremental constitutive equation, w… ▽ More The famous bifurcation analysis performed by Flügge on compressed thin-walled cylinders is based on a series of simplifying assumptions, which allow to obtain the bifurcation landscape, together with explicit expressions for limit behaviours: surface instability, wrinkling, and Euler rod buckling. The most severe assumption introduced by Flügge is the use of an incremental constitutive equation, which does not follow from any nonlinear hyperelastic constitutive law. This is a strong limitation for the applicability of the theory, which becomes questionable when is utilized for a material characterized by a different constitutive equation, such as for instance a Mooney-Rivlin material. We re-derive the entire Flügge's formulation, thus obtaining a framework where any constitutive equation fits. The use of two different nonlinear hyperelastic constitutive equations, referred to compressible materials, leads to incremental equations, which reduce to those derived by Flügge under suitable simplifications. His results are confirmed, together with all the limit equations, now rigorously obtained, and his theory is extended. This extension of the theory of buckling of thin shells allows for computationally efficient determination of bifurcation landscapes for nonlinear constitutive laws, which may for instance be used to model biomechanics of arteries, or soft pneumatic robot arms. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Crack trajectories in materials containing voids via phase-field modelling

Authors: Riccardo Cavuoto, Pietro Lenarda, Diego Misseroni, Marco Paggi, Davide Bigoni

Abstract: Fracture growth in a material is strongly influenced by the presence of inhomogeneities, which deviate crack trajectories from rectilinearity and deeply affect failure. Increasing crack tortuosity is connected to enhancement of fracture toughness, while often a crack may even be stopped when it impinges a void, which releases the stress concentration. Therefore, the determination of crack trajecto… ▽ More Fracture growth in a material is strongly influenced by the presence of inhomogeneities, which deviate crack trajectories from rectilinearity and deeply affect failure. Increasing crack tortuosity is connected to enhancement of fracture toughness, while often a crack may even be stopped when it impinges a void, which releases the stress concentration. Therefore, the determination of crack trajectories is important in the design against failure of materials and mechanical pieces. The recently developed phase-field approach (AT1 and AT2 models), based on a variational approach to damage localization, is believed to be particularly suited to describe complex crack trajectories. This belief is examined through a comparison between simulations and photoelastic experiments on PMMA plates, which have been designed in a new way, to highlight the effects of notches and circular holes on fracture propagation. The latter is shown to initiate from a notch and to be strongly attracted by voids. When a void is hit, fracture is arrested, unless the void contains a notch on its internal surface, from which a new crack nucleates and propagates. Different mechanical models are tested where fracture initiates and grows (i.) under Mode I compact tension, (ii.) four-point bending and (iii.) a tensile stress indirectly generated during compression of samples containing a circular hole. The experiments show that the fracture propagation may be designed to develop in different tortuous paths, involving multiple arrests and secondary nucleation. Simulations performed with an ad hoc implemented version of the AT1 and AT2 phase-field methods (equipped with spectral decomposition, in which a crack is simulated as a highly localized zone of damage accumulation) are shown to be in close agreement with experiments and therefore confirm the validity of the approach and its potentialities for mechanical design. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: In press on the International Journal of Solids and Structures

MSC Class: 74R10

Elastic disk with isoperimetric Cosserat coating

Authors: Matteo Gaibotti, Davide Bigoni, Sofia G. Mogilevskaya

Abstract: A circular elastic disk is coated with an elastic beam, absorbing shear and normal forces without deformation and linearly reacting to a bending moment with a change in curvature. The inexstensibility of the elastic beam introduces an isoperimetric constraint, so that the length of the initial circumference of the disk is constrained to remain fixed during the loading of the disk/coating system. T… ▽ More A circular elastic disk is coated with an elastic beam, absorbing shear and normal forces without deformation and linearly reacting to a bending moment with a change in curvature. The inexstensibility of the elastic beam introduces an isoperimetric constraint, so that the length of the initial circumference of the disk is constrained to remain fixed during the loading of the disk/coating system. The mechanical model for this system is formulated, solved for general loading, and particularized to the case of two equal and opposite traction distributions, each applied on a small boundary segment (thus modelling indentation of a coated fiber). The stress fields, obtained via complex potentials, are shown to evidence a nice correspondence with photoelastic experiments, ad hoc designed and performed. The presented results are useful for the design of coated fibers at the micro and nano scales. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: 35 pages, 15 figures

Buckling vs unilateral constraint for a multistable metamaterial element

Authors: Nikolin Hima, Davide Bigoni, Francesco Dal Corso

Abstract: A structural element is designed and investigated, forming the basis for the development of an elastic multistable metamaterial. The leitmotif of the structural design is the implementation of a strut characterized by a bifurcation occurring at either vanishing tensile or compressive load. It is shown that buckling at null load leads to a mechanical equivalence with a unilateral constraint formula… ▽ More A structural element is designed and investigated, forming the basis for the development of an elastic multistable metamaterial. The leitmotif of the structural design is the implementation of a strut characterized by a bifurcation occurring at either vanishing tensile or compressive load. It is shown that buckling at null load leads to a mechanical equivalence with a unilateral constraint formulation, introducing shocks in dynamics. Towards a future analysis of the latter, the nonlinear quasi-static response is investigated, showing the multistable character of the structure, which may appear as bistable or tetrastable. △ Less

Submitted 18 July, 2022; v1 submitted 4 May, 2022; originally announced May 2022.

Comments: 22 pages, 7 figures

Tensile material instabilities in elastic beam lattices lead to a bounded stability domain

Authors: G. Bordiga, D. Bigoni, A. Piccolroaz

Abstract: Homogenization of the incremental response of grids made up of preloaded elastic rods leads to homogeneous effective continua which may suffer macroscopic instability, occurring at the same time in both the grid and the effective continuum. This instability corresponds to the loss of ellipticity in the effective material and the formation of localized responses as, for instance, shear bands. Using… ▽ More Homogenization of the incremental response of grids made up of preloaded elastic rods leads to homogeneous effective continua which may suffer macroscopic instability, occurring at the same time in both the grid and the effective continuum. This instability corresponds to the loss of ellipticity in the effective material and the formation of localized responses as, for instance, shear bands. Using lattice models of elastic rods, loss of ellipticity has always been found to occur for stress states involving compression of the rods, as usually these structural elements buckle only under compression. In this way, the locus of material stability for the effective solid is unbounded in tension, i.e. the material is always stable for a tensile prestress. A rigorous application of homogenization theory is proposed to show that the inclusion of sliders (constraints imposing axial and rotational continuity, but allowing shear jumps) in the grid of rods leads to loss of ellipticity in tension, so that the locus for material instability becomes bounded. This result explains (i.) how to design elastic materials passible of localization of deformation and shear banding for all radial stress paths; (ii.) how for all these paths a material may fail by developing strain localization and without involving cracking. △ Less

Submitted 2 June, 2022; v1 submitted 3 January, 2022; originally announced January 2022.

Comments: 24 pages, 9 figures. arXiv admin note: text overlap with arXiv:2001.02096

Journal ref: Phil. Trans. R. Soc. A 380 (2022): 20210388

A way to hypo-elastic artificial materials without a strain potential and displaying flutter instability

Authors: G. Bordiga, A. Piccolroaz, D. Bigoni

Abstract: Cauchy-elastic solids include hyper-elasticity and a subset of elastic materials for which the stress does not follow from a scalar strain potential. More in general, hypo-elastic materials are only defined incrementally and comprise Cauchy-elasticity. Infringement of the hyper-elastic 'dogma' is so far unattempted and normally believed to be impossible, as it apparently violates thermodynamics, b… ▽ More Cauchy-elastic solids include hyper-elasticity and a subset of elastic materials for which the stress does not follow from a scalar strain potential. More in general, hypo-elastic materials are only defined incrementally and comprise Cauchy-elasticity. Infringement of the hyper-elastic 'dogma' is so far unattempted and normally believed to be impossible, as it apparently violates thermodynamics, because energy may be produced in closed strain cycles. Contrary to this belief, we show that non-hyper-elastic behavior is possible and we indicate the way to a practical realization of this new concept. In particular, a design paradigm is established for artificial materials where follower forces, so far ignored in homogenization schemes, are introduced as loads prestressing an elastic two-dimensional grid made up of linear elastic rods (reacting to elongation, flexure and shear). A rigorous application of Floquet-Bloch wave asymptotics yields an unsymmetric acoustic tensor governing the \textit{incremental} dynamics of the effective material. The latter is therefore the incremental response of a hypo-elastic solid, which does not follow from a strain potential and thus does not belong to hyper-elasticity. Through the externally applied follower forces (which could be originated via interaction with a fluid, or gas, or by application of Coulomb friction, or non-holonomic constraints), the artificial material may adsorb/release energy from/to the environment and therefore produce energy in a closed strain loop, without violating any rule of thermodynamics. The solid is also shown to display flutter, a material instability corresponding to a Hopf bifurcation, which was advocated as possible in plastic solids, but never experimentally found and so far believed to be impossible in elasticity. △ Less

Submitted 3 January, 2022; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: 27 pages, 13 figures

Journal ref: J. Mech. Phys. Solids, 158 (2022), Article 104665

Flutter instability and Ziegler destabilization paradox for elastic rods subject to non-holonomic constraints

Authors: Alessandro Cazzolli, Francesco Dal Corso, Davide Bigoni

Abstract: Two types of non-holonomic constraints (imposing a prescription on velocity) are analyzed, connected to an end of a (visco)elastic rod, straight in its undeformed configuration. The equations governing the nonlinear dynamics are obtained and then linearized near the trivial equilibrium configuration. The two constraints are shown to lead to the same equations governing the linearized dynamics of t… ▽ More Two types of non-holonomic constraints (imposing a prescription on velocity) are analyzed, connected to an end of a (visco)elastic rod, straight in its undeformed configuration. The equations governing the nonlinear dynamics are obtained and then linearized near the trivial equilibrium configuration. The two constraints are shown to lead to the same equations governing the linearized dynamics of the Beck (or Pfluger) column in one case and of the Reut column in the other. Therefore, although the structural systems are fully conservative (when viscosity is set to zero), they exhibit flutter and divergence instability. In addition, the Ziegler's destabilization paradox is found when dissipation sources are introduced. It follows that these features are proven to be not only a consequence of 'unrealistic non-conservative loads' (as often stated in the literature), rather, the models proposed by Beck, Reut, and Ziegler can exactly describe the linearized dynamics of structures subject to non-holonomic constraints, which are made now fully accessible to experiments. △ Less

Submitted 24 June, 2021; v1 submitted 22 June, 2021; originally announced June 2021.

Comments: 10 pages, 4 figures

Journal ref: Published in Journal of Applied Mechanics (2021), 88(3): 031003

Dynamics of prestressed elastic lattices: homogenization, instabilities, and strain localization

Authors: Giovanni Bordiga, Luigi Cabras, Andrea Piccolroaz, Davide Bigoni

Abstract: A lattice of elastic Rayleigh rods organized in a parallelepiped geometry can be axially loaded up to an arbitrary amount without distortion and then be subject to incremental time-harmonic dynamic motion. At certain threshold levels of axial load, the grillage manifests instabilities and displays non-trivial axial and flexural incremental vibrations. Floquet-Bloch wave asymptotics is used to homo… ▽ More A lattice of elastic Rayleigh rods organized in a parallelepiped geometry can be axially loaded up to an arbitrary amount without distortion and then be subject to incremental time-harmonic dynamic motion. At certain threshold levels of axial load, the grillage manifests instabilities and displays non-trivial axial and flexural incremental vibrations. Floquet-Bloch wave asymptotics is used to homogenize the in-plane mechanical response, so to obtain an equivalent prestressed elastic solid subject to incremental time-harmonic vibration, which includes, as a particular case, the incremental quasi-static response. The equivalent elastic solid is obtained from its acoustic tensor, directly derived from homogenization and shown to be independent of the rods' rotational inertia. Loss of strong ellipticity in the equivalent continuum coincides with macro-bifurcation in the lattice, while micro-bifurcation remains undetected in the continuum and corresponds to a vibration of vanishing frequency of the lowest dispersion branch of the lattice, occurring at finite wavelength. A perturbative approach based on dynamic Green's function is applied to both the lattice and its equivalent continuum. This shows that only macro-instability corresponds to localization of incremental strain, while micro-instabilities occur in modes that spread throughout the whole lattice with an 'explosive' character. Extremely localized mechanical responses are found both in the lattice and in the solid, with the advantage that the former can be easily realized, for instance via 3D printing. In this way, features such as shear band inclination, or the emergence of a single shear band, or competition between micro and macro instabilities become all designable features. The presented results pave the way for the design of architected materials to be used in applications involving extreme deformations. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 49 pages, 31 figures. arXiv admin note: substantial text overlap with arXiv:2001.02151, arXiv:2001.02096

Journal ref: Journal of the Mechanics and Physics of Solids (2021),146, 104198

Non-holonomic constraints inducing flutter insability in structures under conservative loadings

Authors: Alessandro Cazzolli, Francesco Dal Corso, Davide Bigoni

Abstract: Non-conservative loads of the follower type are usually believed to be the source of dynamic instabilities such as flutter and divergence. It is shown that these instabilities (including Hopf bifurcation, flutter, divergence, and destabilizing effects connected to dissipation phenomena) can be obtained in structural systems loaded by conservative forces, as a consequence of the application of non-… ▽ More Non-conservative loads of the follower type are usually believed to be the source of dynamic instabilities such as flutter and divergence. It is shown that these instabilities (including Hopf bifurcation, flutter, divergence, and destabilizing effects connected to dissipation phenomena) can be obtained in structural systems loaded by conservative forces, as a consequence of the application of non-holonomic constraints. These constraints may be realized through a `perfect skate' (or a non-sliding wheel), or, more in general, through the slipless contact between two circular rigid cylinders, one of which is free of rotating about its axis. The motion of the structure produced by these dynamic instabilities may reach a limit cycle, a feature that can be exploited for soft robotics applications, especially for the realization of limbless locomotion. △ Less

Submitted 4 January, 2024; v1 submitted 22 May, 2020; originally announced May 2020.

Journal ref: Journal of the Mechanics and Physics of Solids, Volume 138, May 2020, 103919

Structures loaded with a force acting along a fixed straight line, or the "Reut's column problem"

Authors: Davide Bigoni, Diego Misseroni

Abstract: An elastic double pendulum subject to a force acting along a fixed straight line, the so-called "Reut's column problem", is a structure exhibiting flutter and divergence instability, which was never realized in practice and thus debated whether to represent reality or mere speculation. It is shown, both theoretically and experimentally, how to obtain the Reut's loading by exploiting the contact wi… ▽ More An elastic double pendulum subject to a force acting along a fixed straight line, the so-called "Reut's column problem", is a structure exhibiting flutter and divergence instability, which was never realized in practice and thus debated whether to represent reality or mere speculation. It is shown, both theoretically and experimentally, how to obtain the Reut's loading by exploiting the contact with friction of a rigid blade against a freely-rotating cylindrical constraint, which moves axially at constant speed, an action recalling that of a bow's hair on a violin string. With this experimental set-up, flutter and divergence instabilities, as well as the detrimental effect of viscosity on critical loads, are documented indisputably, thus bringing an end to a long debate. This result opens a new research area, with perspective applications to mechanical actuators, high-precision cutting tools, or energy harvesting devices. △ Less

Submitted 8 February, 2020; originally announced February 2020.

Comments: 18 pages, 13 figures

Journal ref: Journal of the Mechanics and Physics of Solids 134 (2020): 103741

Omnidirectional flexural invisibility of multiple interacting voids in vibrating elastic plates

Authors: D. Misseroni, A. B. Movchan, D. Bigoni

Abstract: In elasticity, the design of a cloaking for an inclusion or a void to leave a vibrational field unperturbed by its presence, so to achieve its invisibility, is a thoroughly analyzed, but still unchallenged, mechanical problem. The 'cloaking transformation' concept, originally developed in electromagnetism and optics, is not directly applicable to elastic waves, displaying a complex vectorial natur… ▽ More In elasticity, the design of a cloaking for an inclusion or a void to leave a vibrational field unperturbed by its presence, so to achieve its invisibility, is a thoroughly analyzed, but still unchallenged, mechanical problem. The 'cloaking transformation' concept, originally developed in electromagnetism and optics, is not directly applicable to elastic waves, displaying a complex vectorial nature. Consequently, all examples of elastic cloaking presented so far involve complex design and thick coating skins. These cloakings often work only for problems of unidirectional propagation, within narrow ranges of frequency, and considering only one colaked object. Here, a new method based on the concept of reinforcement, achieved via elastic stiffening and mass redistribution, is introduced to cloak multiple voids in an elastic plate. This simple technique produces invisibility of the voids to flexural waves within an extremely broad range of frequencies and thus surpassing in many aspects all existing cloaking techniques. The proposed design principle is applicable in mechanical problems ranging from the micro-scale to the scale of civil engineering. For instance, our results show how to design a perforated load-bearing building wall, vibrating during an earthquake exactly as the same wall, but unperforated, a new finding for seismic protection. △ Less

Submitted 7 February, 2020; originally announced February 2020.

Comments: 16 pages, 5 figures

Journal ref: Proceedings of the Royal Society A 475.2229 (2019): 20190283

Bifurcation of elastic solids with sliding interfaces

Authors: D. Bigoni, N. Bordignon, A. Piccolroaz, S. Stupkiewicz

Abstract: Lubricated sliding contact between soft solids is an interesting topic in biomechanics and for the design of small-scale engineering devices. As a model of this mechanical set-up, two elastic nonlinear solids are considered jointed through a frictionless and bilateral surface, so that continuity of the normal component of the Cauchy traction holds across the surface, but the tangential component i… ▽ More Lubricated sliding contact between soft solids is an interesting topic in biomechanics and for the design of small-scale engineering devices. As a model of this mechanical set-up, two elastic nonlinear solids are considered jointed through a frictionless and bilateral surface, so that continuity of the normal component of the Cauchy traction holds across the surface, but the tangential component is null. Moreover, the displacement can develop only in a way that the bodies in contact do neither detach, nor overlap. Surprisingly, this finite strain problem has not been correctly formulated until now, so that this formulation is the objective of the present paper. The incremental equations are shown to be non-trivial and different from previously (and erroneously) employed conditions. In particular, an exclusion condition for bifurcation is derived to show that previous formulations based on frictionless contact or 'spring-type' interfacial conditions are not able to predict bifurcations in tension, while experiments -- one of which, ad hoc designed, is reported -- show that these bifurcations are a reality and become possible when the correct sliding interface model is used. The presented results introduce a methodology for the determination of bifurcations and instabilities occurring during lubricated sliding between soft bodies in contact. △ Less

Submitted 15 January, 2020; originally announced January 2020.

Comments: 25 pages, 12 figures

Incremental constitutive tensors and strain localization for prestressed elastic lattices: Part II -- incremental dynamics

Authors: G. Bordiga, L. Cabras, A. Piccolroaz, D. Bigoni

Abstract: Floquet-Bloch wave asymptotics is used to homogenize the in-plane mechanical response of a periodic grillage of elastic Rayleigh rods, possessing a distributed mass density, together with rotational inertia. The grid is subject to incremental time-harmonic dynamic motion, superimposed to a given state of axial loading of arbitrary magnitude. In contrast to the quasi-static energy match (addressed… ▽ More Floquet-Bloch wave asymptotics is used to homogenize the in-plane mechanical response of a periodic grillage of elastic Rayleigh rods, possessing a distributed mass density, together with rotational inertia. The grid is subject to incremental time-harmonic dynamic motion, superimposed to a given state of axial loading of arbitrary magnitude. In contrast to the quasi-static energy match (addressed in Part I of the present study), the vibrational properties of the lattice are directly represented by the acoustic tensor of the equivalent solid (without passing through the constitutive tensor). The acoustic tensor is shown to be independent of the rods' rotational inertia and allows the analysis of strong ellipticity of the equivalent continuum, evidencing coincidence with macro-bifurcation in the lattice. On the other hand, micro-bifurcation corresponds to a vibration of vanishing frequency of the lowest dispersion branch of the lattice, occurring at finite wavelength. Dynamic homogenization reveals the structure of the acoustic branches close to ellipticity loss and the analysis of forced vibrations (both in physical space and Fourier space) shows low-frequency wave localizations. A comparison is presented between strain localization occurring near ellipticity loss and forced vibration of the lattice, both corresponding to the application of a concentrated pulsating force. The comparison shows that the homogenization technique allows an almost perfect representation of the lattice, with the exception of cases where micro-bifurcation occurs, which is shown to leave the equivalent solid unaffected. Therefore, the presented results pave the way for the design of architected cellular materials to be used in applications where extreme deformations are involved. △ Less

Submitted 25 August, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: 32 pages, 19 figures

Journal ref: Journal of the Mechanics and Physics of Solids (2021),146, 104198

Incremental constitutive tensors and strain localization for prestressed elastic lattices: Part I -- quasi-static response

Authors: G. Bordiga, L. Cabras, A. Piccolroaz, D. Bigoni

Abstract: A lattice of elastic rods organized in a parallelepiped geometry can be axially loaded up to an arbitrary amount without distortion and then be subject to incremental displacements. Using quasi-static homogenization theory, this lattice can be made equivalent to a prestressed elastic solid subject to incremental deformation, in such a way to obtain extremely localized mechanical responses. These r… ▽ More A lattice of elastic rods organized in a parallelepiped geometry can be axially loaded up to an arbitrary amount without distortion and then be subject to incremental displacements. Using quasi-static homogenization theory, this lattice can be made equivalent to a prestressed elastic solid subject to incremental deformation, in such a way to obtain extremely localized mechanical responses. These responses can be analyzed with reference to a mechanical model which can, in principle, be realized, so that features such as for instance shear bands inclination, or emergence of a single shear band, or competition between micro (occurring in the lattice but not in the equivalent solid) and macro (present in both the lattice and the equivalent continuum) instabilities become all designable features. The analysis of localizations is performed using a Green's function-based perturbative approach to highlight the correspondence between micromechanics of the composite and homogenized response of the equivalent solid. The presented results, limited to quasi-static behaviour, provide a new understanding of strain localization in a continuum and open new possibilities for the realization and experimentation of materials exhibiting these extreme mechanical behaviours. Dynamic homogenization and vibrational localization are deferred to Part II of this study. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 30 pages, 16 figures

Configurational forces and nonlinear structural dynamics

Authors: Costanza Armanini, Francesco Dal Corso, Diego Misseroni, Davide Bigoni

Abstract: Configurational, or Eshelby-like, forces are shown to strongly influence the nonlinear dynamics of an elastic rod constrained with a frictionless sliding sleeve at one end and with an attached mass at the other end. The configurational force, generated at the sliding sleeve constraint and proportional to the square of the bending moment realized there, has been so far investigated only under quasi… ▽ More Configurational, or Eshelby-like, forces are shown to strongly influence the nonlinear dynamics of an elastic rod constrained with a frictionless sliding sleeve at one end and with an attached mass at the other end. The configurational force, generated at the sliding sleeve constraint and proportional to the square of the bending moment realized there, has been so far investigated only under quasi-static setting and is now confirmed (through a variational argument) to be present within a dynamic framework. The deep influence of configurational forces on the dynamics is shown both theoretically (through the development of a dynamic nonlinear model in which the rod is treated as a nonlinear spring, obeying the Euler elastica, with negligible inertia) and experimentally (through a specifically designed experimental set-up). During the nonlinear dynamics, the elastic rod may slip alternatively in and out from the sliding sleeve, becoming a sort of nonlinear oscillator displaying a motion eventually ending with the rod completely injected into or completely ejected from the sleeve. The present results may find applications in the dynamics of compliant and extensible devices, for instance, to guide the movement of a retractable and flexible robot arm. △ Less

Submitted 5 September, 2019; originally announced September 2019.

Journal ref: Journal of the Mechanics and Physics of Solids (2019)

Nested Bloch waves in elastic structures with configurational forces

Authors: Francesco Dal Corso, Domenico Tallarico, Natalia V. Movchan, Alexander B. Movchan, Davide Bigoni

Abstract: Small axial and flexural oscillations are analyzed for a periodic and infinite structure, constrained by sliding sleeves and composed of elastic beams. A nested Bloch-Floquet technique is introduced to treat the non-linear coupling between longitudinal and transverse displacements induced by the configurational forces generated at the sliding sleeve ends. The action of configurational forces is sh… ▽ More Small axial and flexural oscillations are analyzed for a periodic and infinite structure, constrained by sliding sleeves and composed of elastic beams. A nested Bloch-Floquet technique is introduced to treat the non-linear coupling between longitudinal and transverse displacements induced by the configurational forces generated at the sliding sleeve ends. The action of configurational forces is shown to play an important role from two perspectives. First, the band gap structure for purely longitudinal vibration is broken so that axial propagation may occur at frequencies that are forbidden in the absence of a transverse oscillation and, second, a flexural oscillation may induce axial resonance, a situation in which the longitudinal vibrations tend to become unbounded. The presented results disclose the possibility of exploiting configurational forces in the design of mechanical devices towards longitudinal actuation from flexural vibrations of small amplitude at given frequency. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: 20 pages, 5 figures

Journal ref: Philosophical Transactions of the Royal Society A (2019) 377: 20190101

Identification of second-gradient elastic materials from planar hexagonal lattices. Part II: Mechanical characteristics and model validation

Authors: G. Rizzi, F. Dal Corso, D. Veber, D. Bigoni

Abstract: Positive definiteness and symmetry of the constitutive tensors describing a second-gradient elastic (SGE) material, which is energetically equivalent to a hexagonal planar lattice made up of axially deformable bars, are analyzed by exploiting the closed form-expressions obtained in part I of the present study in the \lq condensed' form. It is shown that, while the first-order approximation leads t… ▽ More Positive definiteness and symmetry of the constitutive tensors describing a second-gradient elastic (SGE) material, which is energetically equivalent to a hexagonal planar lattice made up of axially deformable bars, are analyzed by exploiting the closed form-expressions obtained in part I of the present study in the \lq condensed' form. It is shown that, while the first-order approximation leads to an isotropic Cauchy material, a second-order identification procedure provides an equivalent model exhibiting non-locality, non-centrosymmetry, and anisotropy. The derivation of the constitutive properties for the SGE from those of the \lq condensed' one (obtained by considering a quadratic remote displacement which generates stress states satisfying equilibrium) is presented. Comparisons between the mechanical responses of the periodic lattice and of the equivalent SGE material under simple shear and uniaxial strain show the efficacy of the proposed identification procedure and therefore validate the proposed constitutive model. This model reveals that, at higher-order, a lattice material can be made equivalent to a second-gradient elastic material exhibiting an internal length, a finding which is now open for applications in micromechanics. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: 25 pages, 11 figures

Journal ref: International Journal of Solids and Structures, 2019, 176-177, 19-35

Identification of second-gradient elastic materials from planar hexagonal lattices. Part I: Analytical derivation of equivalent constitutive tensors

Authors: G. Rizzi, F. Dal Corso, D. Veber, D. Bigoni

Abstract: A second-gradient elastic (SGE) material is identified as the homogeneous solid equivalent to a periodic planar lattice characterized by a hexagonal unit cell, which is made up of three different linear elastic bars ordered in a way that the hexagonal symmetry is preserved and hinged at each node, so that the lattice bars are subject to pure axial strain while bending is excluded. Closed form-expr… ▽ More A second-gradient elastic (SGE) material is identified as the homogeneous solid equivalent to a periodic planar lattice characterized by a hexagonal unit cell, which is made up of three different linear elastic bars ordered in a way that the hexagonal symmetry is preserved and hinged at each node, so that the lattice bars are subject to pure axial strain while bending is excluded. Closed form-expressions for the identified non-local constitutive parameters are obtained by imposing the elastic energy equivalence between the lattice and the continuum solid, under remote displacement conditions having a dominant quadratic component. In order to generate equilibrated stresses, in the absence of body forces, the applied remote displacement has to be constrained, thus leading to the identification in a \lq condensed' form of a higher-order solid, so that imposition of further constraints becomes necessary to fully quantify the equivalent continuum. The identified SGE material reduces to an equivalent Cauchy material only in the limit of vanishing side length of hexagonal unit cell. The analysis of positive definiteness and symmetry of the equivalent constitutive tensors, the derivation of the second-gradient elastic properties from those of the higher-order solid in the \lq condensed' definition, and a numerical validation of the identification scheme are deferred to Part II of this study. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: 23 pages, 5 figures

Journal ref: International Journal of Solids and Structures, 2019, 176-177, 1-18

Prestress tuning of negative refraction and wave channeling from flexural sources

Authors: G. Bordiga, L. Cabras, A. Piccolroaz, D. Bigoni

Abstract: The quest for wave channeling and manipulation has driven a strong research effort on topological and architected materials, capable of propagating localized electromagnetical or mechanical signals. With reference to an elastic structural grid, which elements can sustain both axial and flexural deformations, it is shown that material interfaces can be created with structural properties tuned by pr… ▽ More The quest for wave channeling and manipulation has driven a strong research effort on topological and architected materials, capable of propagating localized electromagnetical or mechanical signals. With reference to an elastic structural grid, which elements can sustain both axial and flexural deformations, it is shown that material interfaces can be created with structural properties tuned by prestress states to achieve total reflection, negative refraction, and strongly localized signal channeling. The achievement of a flat lens and topologically localized modes is demonstrated and tunability of the system allows these properties to hold for a broad range of wavelengths. An ingredient to obtain these effects is the use, suggested here and never attempted before, of concentrated pulsating moments. The important aspect of the proposed method is that states of prestress can be easily removed or changed to tune with continuity the propagational characteristics of the medium, so that a new use of vibration channeling and manipulation is envisaged for elastic materials. △ Less

Submitted 10 January, 2019; originally announced January 2019.

Comments: 10 pages, 5 figures

Free and forced wave propagation in a Rayleigh-beam grid: flat bands, Dirac cones, and vibration localization vs isotropization

Authors: G. Bordiga, L. Cabras, D. Bigoni, A. Piccolroaz

Abstract: In-plane wave propagation in a periodic rectangular grid beam structure, which includes rotational inertia (so-called 'Rayleigh beams'), is analyzed both with a Floquet-Bloch exact formulation for free oscillations and with a numerical treatment (developed with PML absorbing boundary conditions) for forced vibrations (including Fourier representation and energy flux evaluations), induced by a conc… ▽ More In-plane wave propagation in a periodic rectangular grid beam structure, which includes rotational inertia (so-called 'Rayleigh beams'), is analyzed both with a Floquet-Bloch exact formulation for free oscillations and with a numerical treatment (developed with PML absorbing boundary conditions) for forced vibrations (including Fourier representation and energy flux evaluations), induced by a concentrated force or moment. A complex interplay is observed between axial and flexural vibrations (not found in the common idealization of out-of-plane motion), giving rise to several forms of vibration localization: 'X-', 'cross-' and 'star-' shaped, and channel propagation. These localizations are triggered by several factors, including rotational inertia and slenderness of the beams and the type of forcing source (concentrated force or moment). Although the considered grid of beams introduces an orthotropy in the mechanical response, a surprising 'isotropization' of the vibration is observed at special frequencies. Moreover, rotational inertia is shown to 'sharpen' degeneracies related to Dirac cones (which become more pronounced when the aspect ratio of the grid is increased), while the slenderness can be tuned to achieve a perfectly flat band in the dispersion diagram. The obtained results can be exploited in the realization of metamaterials designed to control wave propagation. △ Less

Submitted 17 July, 2018; v1 submitted 25 June, 2018; originally announced June 2018.

Comments: 25 pages, 20 figures

Flutter and divergence instability in the Pfluger column: experimental evidence of the Ziegler destabilization paradox

Authors: Davide Bigoni, Oleg N. Kirillov, Diego Misseroni, Giovanni Noselli, Mirko Tommasini

Abstract: Flutter instability in elastic structures subject to follower load, the most important cases being the famous Beck's and Pfluger's columns (two elastic rods in a cantilever configuration, with an additional concentrated mass at the end of the rod in the latter case), have attracted, and still attract, a thorough research interest. In this field, the most important issue is the validation of the mo… ▽ More Flutter instability in elastic structures subject to follower load, the most important cases being the famous Beck's and Pfluger's columns (two elastic rods in a cantilever configuration, with an additional concentrated mass at the end of the rod in the latter case), have attracted, and still attract, a thorough research interest. In this field, the most important issue is the validation of the model itself of follower force, a nonconservative action which was harshly criticized and never realized in practice for structures with diffused elasticity. An experimental setup to introduce follower tangential forces at the end of an elastic rod was designed, realized, validated, and tested, in which the follower action is produced by exploiting Coulomb friction on an element (a freely-rotating wheel) in sliding contact against a flat surface (realized by a conveyor belt). It is therefore shown that follower forces can be realized in practice and the first experimental evidence is given for both the flutter and divergence instabilities occurring in the Pfluger's column. In particular, load thresholds for the two instabilities are measured and the detrimental effect of dissipation on the critical load for flutter is experimentally demonstrated, while a slight increase in load is found for the divergence instability. The presented approach to follower forces discloses new horizons for testing self-oscillating structures and for exploring and documenting dynamic instabilities possible when nonconservative loads are applied. △ Less

Submitted 19 April, 2018; originally announced April 2018.

Comments: 28 pages, 12 figures

Journal ref: Journal of the Mechanics and Physics of Solids (2018) 116: 91-116

Identification of higher-order continua equivalent to a Cauchy elastic composite

Authors: A. Bacigalupo, M. Paggi, F. Dal Corso, D. Bigoni

Abstract: A heterogeneous Cauchy elastic material may display micromechanical effects that can be modeled in a homogeneous equivalent material through the introduction of higher-order elastic continua. Asymptotic homogenization techniques provide an elegant and rigorous route to the evaluation of equivalent higher-order materials, but are often of difficult and awkward practical implementation. On the other… ▽ More A heterogeneous Cauchy elastic material may display micromechanical effects that can be modeled in a homogeneous equivalent material through the introduction of higher-order elastic continua. Asymptotic homogenization techniques provide an elegant and rigorous route to the evaluation of equivalent higher-order materials, but are often of difficult and awkward practical implementation. On the other hand, identification techniques, though relying on simplifying assumptions, are of straightforward use. A novel strategy for the identification of equivalent second-gradient Mindlin solids is proposed in an attempt to combine the accuracy of asymptotic techniques with the simplicity of identification approaches. Following the asymptotic homogenization scheme, the overall behaviour is defined via perturbation functions, which (differently from the asymptotic scheme) are evaluated on a finite domain obtained as the periodic repetition of cells and subject to quadratic displacement boundary conditions. As a consequence, the periodicity of the perturbation function is satisfied only in an approximate sense, nevertheless results from the proposed identification algorithm are shown to be reasonably accurate. △ Less

Submitted 19 April, 2018; originally announced April 2018.

The dynamics of a shear band

Authors: Diana Giarola, Domenico Capuani, Davide Bigoni

Abstract: A shear band of finite length, formed inside a ductile material at a certain stage of a con- tinued homogeneous strain, provides a dynamic perturbation to an incident wave field, which strongly influences the dynamics of the material and affects its path to failure. The investigation of this perturbation is presented for a ductile metal, with reference to the incremental mechanics of a material ob… ▽ More A shear band of finite length, formed inside a ductile material at a certain stage of a con- tinued homogeneous strain, provides a dynamic perturbation to an incident wave field, which strongly influences the dynamics of the material and affects its path to failure. The investigation of this perturbation is presented for a ductile metal, with reference to the incremental mechanics of a material obeying the J 2-deformation theory of plasticity (a special form of prestressed, elastic, anisotropic, and incompressible solid). The treatment originates from the derivation of integral representations relating the incremental mechan- ical fields at every point of the medium to the incremental displacement jump across the shear band faces, generated by an impinging wave. The boundary integral equations (under the plane strain assumption) are numerically approached through a collocation technique, which keeps into account the singularity at the shear band tips and permits the analysis of an incident wave impinging a shear band. It is shown that the presence of the shear band induces a resonance, visible in the incremental displacement field and in the stress intensity factor at the shear band tips, which promotes shear band growth. Moreover, the waves scattered by the shear band are shown to generate a fine texture of vibrations, par- allel to the shear band line and propagating at a long distance from it, but leaving a sort of conical shadow zone, which emanates from the tips of the shear band. △ Less

Submitted 4 January, 2018; originally announced January 2018.

Detecting singular weak-dissipation limit for flutter onset in reversible systems

Authors: Davide Bigoni, Oleg N. Kirillov, Diego Misseroni, Giovanni Noselli, Mirko Tommasini

Abstract: A `flutter machine' is introduced for the investigation of a singular interface between the classical and reversible Hopf bifurcations that is theoretically predicted to be generic in nonconservative reversible systems with vanishing dissipation. In particular, such a singular interface exists for the Pflüger viscoelastic column moving in a resistive medium, which is proven by means of the perturb… ▽ More A `flutter machine' is introduced for the investigation of a singular interface between the classical and reversible Hopf bifurcations that is theoretically predicted to be generic in nonconservative reversible systems with vanishing dissipation. In particular, such a singular interface exists for the Pflüger viscoelastic column moving in a resistive medium, which is proven by means of the perturbation theory of multiple eigenvalues with the Jordan block. The laboratory setup, consisting of a cantilevered viscoelastic rod loaded by a positional force with non-zero curl produced by dry friction, demonstrates high sensitivity of the classical Hopf bifurcation onset {to the ratio between} the weak air drag and Kelvin-Voigt damping in the Pflüger column. Thus, the Whitney umbrella singularity is experimentally confirmed, responsible for discontinuities accompanying dissipation-induced instabilities in a broad range of physical contexts. △ Less

Submitted 20 January, 2018; v1 submitted 7 June, 2017; originally announced June 2017.

Comments: 17 pages, 10 figures, typos corrected, refs added. Accepted to PRE

Journal ref: Phys. Rev. E 97, 023003 (2018)

The destabilizing effect of external damping: Singular flutter boundary for the Pfluger column with vanishing external dissipation

Authors: Mirko Tommasini, Oleg N. Kirillov, Diego Misseroni, Davide Bigoni

Abstract: Elastic structures loaded by nonconservative positional forces are prone to instabilities induced by dissipation: it is well-known in fact that internal viscous damping destabilizes the marginally stable Ziegler's pendulum and Pfluger column (of which the Beck's column is a special case), two structures loaded by a tangential follower force. The result is the so-called 'destabilization paradox', w… ▽ More Elastic structures loaded by nonconservative positional forces are prone to instabilities induced by dissipation: it is well-known in fact that internal viscous damping destabilizes the marginally stable Ziegler's pendulum and Pfluger column (of which the Beck's column is a special case), two structures loaded by a tangential follower force. The result is the so-called 'destabilization paradox', where the critical force for flutter instability decreases by an order of magnitude when the coefficient of internal damping becomes infinitesimally small. Until now external damping, such as that related to air drag, is believed to provide only a stabilizing effect, as one would intuitively expect. Contrary to this belief, it will be shown that the effect of external damping is qualitatively the same as the effect of internal damping, yielding a pronounced destabilization paradox. Previous results relative to destabilization by external damping of the Ziegler's and Pfluger's elastic structures are corrected in a definitive way leading to a new understanding of the destabilizating role played by viscous terms. △ Less

Submitted 1 October, 2016; originally announced November 2016.

Comments: 19 pages, 6 figures

Journal ref: Journal of the Mechanics and Physics of Solids (2016) 91: 204-215

Isotoxal star-shaped polygonal voids and rigid inclusions in nonuniform antiplane shear fields. Part II: Singularities, annihilation and invisibility

Authors: Francesco Dal Corso, Summer Shahzad, Davide Bigoni

Abstract: Notch stress intensity factors and stress intensity factors are obtained analytically for isotoxal star-shaped polygonal voids and rigid inclusions (and also for the corresponding limit cases of star-shaped cracks and stiffeners), when loaded through remote inhomogeneous (self-equilibrated, polynomial) antiplane shear stress in an infinite linear elastic matrix. Usually these solutions show stress… ▽ More Notch stress intensity factors and stress intensity factors are obtained analytically for isotoxal star-shaped polygonal voids and rigid inclusions (and also for the corresponding limit cases of star-shaped cracks and stiffeners), when loaded through remote inhomogeneous (self-equilibrated, polynomial) antiplane shear stress in an infinite linear elastic matrix. Usually these solutions show stress singularities at the inclusion corners. It is shown that an infinite set of geometries and loading conditions exist for which not only the singularity is absent, but the stress vanishes ('annihilates') at the corners. Thus the material, which even without the inclusion corners would have a finite stress, remains unstressed at these points in spite of the applied remote load. Moreover, similar conditions are determined in which a star-shaped crack or stiffener leaves the ambient stress completely unperturbed, thus reaching a condition of 'quasi-static invisibility'. Stress annihilation and invisibility define optimal loading modes for the overall strength of a composite and are useful for designing ultra-resistant materials. △ Less

Submitted 24 March, 2016; originally announced May 2016.

Comments: 21 pages, 11 figures

Journal ref: International Journal of Solids and Structures 85-86 (2016), 76-88

Isotoxal star-shaped polygonal voids and rigid inclusions in nonuniform antiplane shear fields. Part I: Formulation and full-field solution

Authors: Francesco Dal Corso, Summer Shahzad, Davide Bigoni

Abstract: An infinite class of nonuniform antiplane shear fields is considered for a linear elastic isotropic space and (non-intersecting) isotoxal star-shaped polygonal voids and rigid inclusions perturbing these fields are solved. Through the use of the complex potential technique together with the generalized binomial and the multinomial theorems, full-field closed-form solutions are obtained in the conf… ▽ More An infinite class of nonuniform antiplane shear fields is considered for a linear elastic isotropic space and (non-intersecting) isotoxal star-shaped polygonal voids and rigid inclusions perturbing these fields are solved. Through the use of the complex potential technique together with the generalized binomial and the multinomial theorems, full-field closed-form solutions are obtained in the conformal plane. The particular (and important) cases of star-shaped cracks and rigid-line inclusions (stiffeners) are also derived. Except for special cases (addressed in Part II), the obtained solutions show singularities at the inclusion corners and at the crack and stiffener ends, where the stress blows-up to infinity, and is therefore detrimental to strength. It is for this reason that the closed-form determination of the stress field near a sharp inclusion or void is crucial for the design of ultra-resistant composites. △ Less

Submitted 24 March, 2016; originally announced April 2016.

Comments: 19 pages, 6 figures, 1 table

Journal ref: International Journal of Solids and Structures 85-86 (2016), 67-75

A Stabilised Nodal Spectral Element Method for Fully Nonlinear Water Waves

Authors: Allan Peter Engsig-Karup, Claes Eskilsson, Daniele Bigoni

Abstract: We present an arbitrary-order spectral element method for general-purpose simulation of non-overturning water waves, described by fully nonlinear potential theory. The method can be viewed as a high-order extension of the classical finite element method proposed by Cai et al (1998) \cite{CaiEtAl1998}, although the numerical implementation differs greatly. Features of the proposed spectral element… ▽ More We present an arbitrary-order spectral element method for general-purpose simulation of non-overturning water waves, described by fully nonlinear potential theory. The method can be viewed as a high-order extension of the classical finite element method proposed by Cai et al (1998) \cite{CaiEtAl1998}, although the numerical implementation differs greatly. Features of the proposed spectral element method include: nodal Lagrange basis functions, a general quadrature-free approach and gradient recovery using global $L^2$ projections. The quartic nonlinear terms present in the Zakharov form of the free surface conditions can cause severe aliasing problems and consequently numerical instability for marginally resolved or very steep waves. We show how the scheme can be stabilised through a combination of over-integration of the Galerkin projections and a mild spectral filtering on a per element basis. This effectively removes any aliasing driven instabilities while retaining the high-order accuracy of the numerical scheme. The additional computational cost of the over-integration is found insignificant compared to the cost of solving the Laplace problem. The model is applied to several benchmark cases in two dimensions. The results confirm the high order accuracy of the model (exponential convergence), and demonstrate the potential for accuracy and speedup. The results of numerical experiments are in excellent agreement with both analytical and experimental results for strongly nonlinear and irregular dispersive wave propagation. The benefit of using a high-order -- possibly adapted -- spatial discretization for accurate water wave propagation over long times and distances is particularly attractive for marine hydrodynamics applications. △ Less

Submitted 29 April, 2016; v1 submitted 8 December, 2015; originally announced December 2015.

Comments: Accepted for publication in Journal of Computational Physics April 29, 2016

ACM Class: G.1

An Elastica Arm Scale

Authors: F. Bosi, D. Misseroni, F. Dal Corso, D. Bigoni

Abstract: The concept of 'deformable arm scale' (completely different from a traditional rigid arm balance) is theoretically introduced and experimentally validated. The idea is not intuitive, but is the result of nonlinear equilibrium kinematics of rods inducing configurational forces, so that deflection of the arms becomes necessary for the equilibrium, which would be impossible for a rigid system. In par… ▽ More The concept of 'deformable arm scale' (completely different from a traditional rigid arm balance) is theoretically introduced and experimentally validated. The idea is not intuitive, but is the result of nonlinear equilibrium kinematics of rods inducing configurational forces, so that deflection of the arms becomes necessary for the equilibrium, which would be impossible for a rigid system. In particular, the rigid arms of usual scales are replaced by a flexible elastic lamina, free of sliding in a frictionless and inclined sliding sleeve, which can reach a unique equilibrium configuration when two vertical dead loads are applied. Prototypes realized to demonstrate the feasibility of the system show a high accuracy in the measure of load within a certain range of use. It is finally shown that the presented results are strongly related to snaking of confined beams, with implications on locomotion of serpents, plumbing, and smart oil drilling. △ Less

Submitted 18 September, 2015; originally announced September 2015.

Comments: 11 pages, 5 figures

Journal ref: Proceedings of the Royal Society A 2014 (470): 20140232

The deformation of an elastic rod with a clamp sliding along a smooth and curved profile

Authors: D. Misseroni, G. Noselli, D. Zaccaria, D. Bigoni

Abstract: The design of compliant mechanisms is crucial in several technologies and relies on the availability of solutions for nonlinear structural problems. One of these solutions is given and experimentally validated in the present article for a compliant mechanism moving along a smooth curved profile. In particular, a deformable elastic rod is held by two clamps, one at each end. The first clamp is cons… ▽ More The design of compliant mechanisms is crucial in several technologies and relies on the availability of solutions for nonlinear structural problems. One of these solutions is given and experimentally validated in the present article for a compliant mechanism moving along a smooth curved profile. In particular, a deformable elastic rod is held by two clamps, one at each end. The first clamp is constrained to slide without friction along a curved profile, while the second clamp moves in a straight line transmitting its motion through the elastic rod to the first clamp. For this system it is shown that the clamp sliding on the profile imposes nontrivial boundary conditions (derived via a variational and an asymptotic approach), which strongly influence buckling and nonlinear structural behaviour. Investigation of this behaviour shows that a compliant mechanism can be designed, which gives an almost neutral response in compression. This behavior could easily be exploited to make a force limiting device. Finally a proof-of-concept device was constructed and tested showing that the analyzed mechanical system can be realized in practice and it behaves tightly to the model, so that it can now be used in the design of machines that use compliant mechanisms. △ Less

Submitted 18 September, 2015; originally announced September 2015.

Comments: 15 pages, 9 figures

Journal ref: International Journal of Solids and Structures 2015 (69-70) 491-497

Efficient uncertainty quantification of a fully nonlinear and dispersive water wave model with random inputs

Authors: Daniele Bigoni, Allan P. Engsig-Karup, Claes Eskilsson

Abstract: A major challenge in next-generation industrial applications is to improve numerical analysis by quantifying uncertainties in predictions. In this work we present a formulation of a fully nonlinear and dispersive potential flow water wave model with random inputs for the probabilistic description of the evolution of waves. The model is analyzed using random sampling techniques and non-intrusive me… ▽ More A major challenge in next-generation industrial applications is to improve numerical analysis by quantifying uncertainties in predictions. In this work we present a formulation of a fully nonlinear and dispersive potential flow water wave model with random inputs for the probabilistic description of the evolution of waves. The model is analyzed using random sampling techniques and non-intrusive methods based on generalized Polynomial Chaos (PC). These methods allow to accurately and efficiently estimate the probability distribution of the solution and require only the computation of the solution in different points in the parameter space, allowing for the reuse of existing simulation software. The choice of the applied methods is driven by the number of uncertain input parameters and by the fact that finding the solution of the considered model is computationally intensive. We revisit experimental benchmarks often used for validation of deterministic water wave models. Based on numerical experiments and assumed uncertainties in boundary data, our analysis reveals that some of the known discrepancies from deterministic simulation in comparison with experimental measurements could be partially explained by the variability in the model input. We finally present a synthetic experiment studying the variance based sensitivity of the wave load on an off-shore structure to a number of input uncertainties. In the numerical examples presented the PC methods have exhibited fast convergence, suggesting that the problem is amenable to being analyzed with such methods. △ Less

Submitted 14 April, 2016; v1 submitted 23 October, 2014; originally announced October 2014.

MSC Class: 65C99; 76B15