Source author record

Davide Bigoni

Davide Bigoni appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2024arXiv

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

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.

preprint2023arXiv

Fusion of two stable elastic structures resulting in an unstable system

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.

preprint2022arXiv

Buckling of thin-walled cylinders from three dimensional nonlinear elasticity

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.

preprint2022arXiv

Buckling vs unilateral constraint for a multistable metamaterial element

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.

preprint2022arXiv

Crack trajectories in materials containing voids via phase-field modelling

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Elastic disk with isoperimetric Cosserat coating

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2016arXiv

Application of Tomographic Reconstruction Techniques for Density Analysis of Green Bodies

Progress in the manufacturing of ceramics, but also of sintered metals, strongly relies on the evaluation of the density distribution in green bodies. This evaluation is crucial from many points of view, including the calibration of constitutive models for \textit{in-silico} simulation of densification processes. To this end, X-ray tomography and other techniques are possible but can be unmanageable for some institutions. Therefore, a destructive method is introduced in the present article to measure the density field of a green body sample using a CNC mill, an analytical balance, and analysis techniques from the field of computational tomography. A virtual experiment is presented where the method is used to reconstruct a simulated green body density field and is found to satisfactorily correspond to the original solution. The green body density field of a truncated cylinder made of alumina powder is evaluated using this method and the reconstructed field is presented.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Stress channelling in extreme couple-stress materials Part I: Strong ellipticity, wave propagation, ellipticity, and discontinuity relations

Materials with extreme mechanical anisotropy are designed to work near a material instability threshold where they display stress channelling and strain localization, effects that can be exploited in several technologies. Extreme couple stress solids are introduced and for the first time systematically analyzed in terms of several material instability criteria: positive-definiteness of the strain energy (implying uniqueness of the mixed b.v.p.), strong ellipticity (implying uniqueness of the b.v.p. with prescribed kinematics on the whole boundary), plane wave propagation, ellipticity, and the emergence of discontinuity surfaces. Several new and unexpected features are highlighted: (i.) Ellipticity is mainly dictated by the 'Cosserat part' of the elasticity and (ii.) its failure is shown to be related to the emergence of discontinuity surfaces; (iii.) Ellipticity and wave propagation are not interdependent conditions (so that it is possible for waves not to propagate when the material is still in the elliptic range and, in very special cases, for waves to propagate when ellipticity does not hold). The proof that loss of ellipticity induces stress channelling, folding and faulting of an elastic Cosserat continuum (and the related derivation of the infinite-body Green's function under antiplane strain conditions) is deferred to Part II of this study.

preprint2015arXiv

Stress channelling in extreme couple-stress materials Part II: Localized folding vs faulting of a continuum in single and cross geometries

The antiplane strain Green's functions for an applied concentrated force and moment are obtained for Cosserat elastic solids with extreme anisotropy, which can be tailored to bring the material in a state close to an instability threshold such as failure of ellipticity. It is shown that the wave propagation condition (and not ellipticity) governs the behavior of the antiplane strain Green's functions. These Green's functions are used as perturbing agents to demonstrate in an extreme material the emergence of localized (single and cross) stress channelling and the emergence of antiplane localized folding (or creasing, or weak elastostatic shock) and faulting (or elastostatic shock) of a Cosserat continuum, phenomena which remain excluded for a Cauchy elastic material. During folding some components of the displacement gradient suffer a finite jump, whereas during faulting the displacement itself displays a finite discontinuity.

preprint2014arXiv

Eshelby-like forces acting on elastic structures: theoretical and experimental proof

The Eshelbian (or configurational) force is the main concept of a celebrated theoretical framework associated with the motion of dislocations and, more in general, defects in solids. In a similar vein, in an elastic structure where a (smooth and bilateral) constraint can move and release energy, a force driving the configuration is generated, which therefore is called by analogy 'Eshelby-like' or 'configurational'. This force (generated by a specific movable constraint) is derived both via variational calculus and, independently, through an asymptotic approach. Its action on the elastic structure is counterintuitive, but is fully substantiated and experimentally measured on a model structure that we have designed, realized and tested. These findings open a totally new perspective in the mechanics of deformable mechanisms, with possible broad applications, even at the nanoscale.

preprint2014arXiv

Stress concentration near stiff inclusions: validation of rigid inclusion model and boundary layers by means of photoelasticity

Photoelasticity is employed to investigate the stress state near stiff rectangular and rhombohedral inclusions embedded in a 'soft' elastic plate. Results show that the singular stress field predicted by the linear elastic solution for the rigid inclusion model can be generated in reality, with great accuracy, within a material. In particular, experiments: (i.) agree with the fact that the singularity is lower for obtuse than for acute inclusion angles; (ii.) show that the singularity is stronger in Mode II than in Mode I (differently from a notch); (iii.) validate the model of rigid quadrilateral inclusion; (iv.) for thin inclusions, show the presence of boundary layers deeply influencing the stress field, so that the limit case of rigid line inclusion is obtained in strong dependence on the inclusion's shape. The introduced experimental methodology opens the possibility of enhancing the design of thin reinforcements and of analyzing complex situations involving interaction between inclusions and defects.

preprint2013arXiv

Anisotropic effective higher-order response of heterogeneous Cauchy elastic materials

The homogenization results obtained by Bacca et al. (Homogenization of heterogeneous Cauchy-elastic materials leads to Mindlin second-gradient elasticity. Part I: Closed form expression for the effective higher-order constitutive tensor. arXiv:1305.2365 Submitted, 2013), to define effective second-gradient elastic materials from heterogeneous Cauchy elastic solids, are extended here to the case of phases having non-isotropic tensors of inertia. It is shown that the nonlocal constitutive tensor for the homogenized material depends on both the inertia properties of the RVE and the difference between the effective and the matrix local elastic tensors. Results show that: (i.) a composite material can be designed to result locally isotropic but nonlocally orthotropic; (ii.) orthotropic nonlocal effects are introduced when a dilute distribution of aligned elliptical holes and, in the limit case, of cracks is homogenized.

preprint2010arXiv

A dynamical interpretation of flutter instability in a continuous medium

Flutter instability in an infinite medium is a form of material instability corresponding to the occurrence of complex conjugate squares of the acceleration wave velocities. Although its occurrence is known to be possible in elastoplastic materials with nonassociative flow law and to correspond to some dynamically growing disturbance, its mechanical meaning has to date still eluded a precise interpretation. This is provided here by constructing the infinite-body, time-harmonic Green's function for the loading branch of an elastoplastic material in flutter conditions. Used as a perturbation, it reveals that flutter corresponds to a spatially blowing-up disturbance, exhibiting well-defined directional properties, determined by the wave directions for which the eigenvalues become complex conjugate. Flutter is shown to be connected to the formation of localized deformations, a dynamical phenomenon sharing geometrical similarities with the well-known mechanism of shear banding occurring under quasi-static loading. Flutter may occur much earlier than shear banding in a process of continued plastic deformation.

preprint2010arXiv

An elastoplastic framework for granular materials becoming cohesive through mechanical densification. Part I - small strain formulation

Mechanical densification of granular bodies is a process in which a loose material becomes increasingly cohesive as the applied pressure increases. A constitutive description of this process faces the formidable problem that granular and dense materials have completely different mechanical behaviours (nonlinear elastic properties, yield limit, plastic flow and hardening laws), which must both be, in a sense, included in the formulation. A treatment of this problem is provided here, so that a new phenomenological, elastoplastic constitutive model is formulated, calibrated by experimental data, implemented and tested, that is capable of describing the transition between granular and fully dense states of a given material. The formulation involves a novel use of elastoplastic coupling to describe the dependence of cohesion and elastic properties on the plastic strain. The treatment falls within small strain theory, which is thought to be appropriate in several situations; however, a generalization of the model to large strain is provided in Part II of this paper.

preprint2010arXiv

An elastoplastic framework for granular materials becoming cohesive through mechanical densification. Part II - the formulation of elastoplastic coupling at large strain

The two key phenomena occurring in the process of ceramic powder compaction are the progressive gain in cohesion and the increase of elastic stiffness, both related to the development of plastic deformation. The latter effect is an example of `elastoplastic coupling', in which the plastic flow affects the elastic properties of the material, and has been so far considered only within the framework of small strain assumption (mainly to describe elastic degradation in rock-like materials), so that it remains completely unexplored for large strain. Therefore, a new finite strain generalization of elastoplastic coupling theory is given to describe the mechanical behaviour of materials evolving from a granular to a dense state. The correct account of elastoplastic coupling and of the specific characteristics of materials evolving from a loose to a dense state (for instance, nonlinear --or linear-- dependence of the elastic part of the deformation on the forming pressure in the granular --or dense-- state) makes the use of existing large strain formulations awkward, if even possible. Therfore, first, we have resorted to a very general setting allowing general transformations between work-conjugate stress and strain measures; second, we have introduced the multiplicative decomposition of the deformation gradient and, third, employing isotropy and hyperelasticity of elastic response, we have obtained a relation between the Biot stress and its `total' and `plastic' work-conjugate strain measure. This is a key result, since it allows an immediate achievement of the rate elastoplastic constitutive equations. Knowing the general form of these equations, all the specific laws governing the behaviour of ceramic powders are finally introduced as generalizations of the small strain counterparts given in Part I of this paper.

preprint2010arXiv

Yield criteria for quasibrittle and frictional materials

A new yield/damage function is proposed for modelling the inelastic behaviour of a broad class of pressure-sensitive, frictional, ductile and brittle-cohesive materials. The yield function allows the possibility of describing a transition between the shape of a yield surface typical of a class of materials to that typical of another class of materals. This is a fundamental key to model the behaviour of materials which become cohesive during hardening (so that the shape of the yield surface evolves from that typical of a granular material to that typical of a dense material), or which decrease cohesion due to damage accumulation. The proposed yield function is shown to agree with a variety of experimental data relative to soil, concrete, rock, metallic and composite powders, metallic foams, porous metals, and polymers. The yield function represents a single, convex and smooth surface in stress space approaching as limit situations well-known criteria and the extreme limits of convexity in the deviatoric plane. The yield function is therefore a generalization of several criteria, including von Mises, Drucker-Prager, Tresca, modified Tresca, Coulomb-Mohr, modified Cam-clay, and --concerning the deviatoric section-- Rankine and Ottosen. Convexity of the function is proved by developing two general propositions relating convexity of the yield surface to convexity of the corresponding function. These propositions are general and therefore may be employed to generate other convex yield functions.