Institut des Sciences de la Mécanique et Applications Industrielles

This paper aims at reviewing nonlinear methods for model order reduction of structures with geometric nonlinearity, with a special emphasis on the techniques based on invariant manifold theory. Nonlinear methods differ from linear based techniques by their use of a nonlinear mapping instead of adding new vectors to enlarge the projection basis. Invariant manifolds have been first introduced in vibration theory within the context of nonlinear normal modes (NNMs) and have been initially computed from the modal basis, using either a graph representation or a normal form approach to compute mappings and reduced dynamics. These developments are first recalled following a historical perspective, where the main applications were first oriented toward structural models that can be expressed thanks to partial differential equations (PDE). They are then replaced in the more general context of the parametrisation of invariant manifold that allows unifying the approaches. Then the specific case of structures discretized with the finite element method is addressed. Implicit condensation, giving rise to a projection onto a stress manifold, and modal derivatives, used in the framework of the quadratic manifold, are first reviewed. Finally, recent developments allowing direct computation of reduced-order models (ROMs) relying on invariant manifolds theory are detailed. Applicative examples are shown and the extension of the methods to deal with further complications are reviewed. Finally, open problems and future directions are highlighted.

Summary In this contribution, we propose a dynamic gradient damage model as a phase‐field approach for studying brutal fracture phenomena in quasi‐brittle materials under impact‐type loading conditions. Several existing approaches to account for the tension–compression asymmetry of fracture behavior of materials are reviewed. A better understanding of these models is provided through a uniaxial traction experiment. We then give an efficient numerical implementation of the model in an explicit dynamics context. Simulations results obtained with parallel computing are discussed both from a computational and physical point of view. Different damage constitutive laws and tension–compression asymmetry formulations are compared with respect to their aptitude to approximate brittle fracture. Copyright © 2016 John Wiley & Sons, Ltd.

Abstract This paper investigates model-order reduction methods for geometrically nonlinear structures. The parametrisation method of invariant manifolds is used and adapted to the case of mechanical systems in oscillatory form expressed in the physical basis, so that the technique is directly applicable to mechanical problems discretised by the finite element method. Two nonlinear mappings, respectively related to displacement and velocity, are introduced, and the link between the two is made explicit at arbitrary order of expansion, under the assumption that the damping matrix is diagonalised by the conservative linear eigenvectors. The same development is performed on the reduced-order dynamics which is computed at generic order following different styles of parametrisation. More specifically, three different styles are introduced and commented: the graph style, the complex normal form style and the real normal form style. These developments allow making better connections with earlier works using these parametrisation methods. The technique is then applied to three different examples. A clamped-clamped arch with increasing curvature is first used to show an example of a system with a softening behaviour turning to hardening at larger amplitudes, which can be replicated with a single mode reduction. Secondly, the case of a cantilever beam is investigated. It is shown that invariant manifold of the first mode shows a folding point at large amplitudes. This exemplifies the failure of the graph style due to the folding point on a real structure, whereas the normal form style is able to pass over the folding. Finally, a MEMS (Micro Electro Mechanical System) micromirror undergoing large rotations is used to show the importance of using high-order expansions on an industrial example.

Abstract Dimensionality reduction in mechanical vibratory systems poses challenges for distributed structures including geometric nonlinearities, mainly because of the lack of invariance of the linear subspaces. A reduction method based on direct normal form computation for large finite element (FE) models is here detailed. The main advantage resides in operating directly from the physical space, hence avoiding the computation of the complete eigenfunctions spectrum. Explicit solutions are given, thus enabling a fully non-intrusive version of the reduction method. The reduced dynamics is obtained from the normal form of the geometrically nonlinear mechanical problem, free of non-resonant monomials, and truncated to the selected master coordinates, thus making a direct link with the parametrisation of invariant manifolds. The method is fully expressed with a complex-valued formalism by detailing the homological equations in a systematic manner, and the link with real-valued expressions is established. A special emphasis is put on the treatment of second-order internal resonances and the specific case of a 1:2 resonance is made explicit. Finally, applications to large-scale models of micro-electro-mechanical structures featuring 1:2 and 1:3 resonances are reported, along with considerations on computational efficiency.

We inspect the propagation of shear polarized surface waves akin to Love waves through a forest of trees of the same height atop a guiding layer on a soil substrate. An asymptotic analysis shows that the forest behaves like an infinitely anisotropic wedge with effective boundary conditions. We discover that the foliage of trees brings a radical change in the nature of the dispersion relation of these surface waves, which behave like spoof plasmons in the limit of a vanishing guiding layer, and like Love waves in the limit of trees with a vanishing height. When we consider a forest with trees of increasing or decreasing height, this hybrid ``spoof Love wave'' is either trapped within the trees or converted into a downward propagating bulk (shear) wave. These mechanisms of wave trapping and wave conversion appear to be robust with respect to perturbations of height or position of trees in the metawedge and with respect to three-dimensional effects such as regarding a potential change of elastic wave polarization.

Wake modes of a three-dimensional blunt-based body near a wall are investigated at a Reynolds number states in the horizontal (respectively vertical) direction, while bistable dynamics persists for the symmetric LR and TB configurations. The shape of periodic modes follows the changes in wake static orientation. The transition between the two regimes is governed by both the total injected flow rate and the location of the injection.

International audience

Abstract Quasi-periodic solutions can arise in assemblies of nonlinear oscillators as a consequence of Neimark-Sacker bifurcations. In this work, the appearance of Neimark-Sacker bifurcations is investigated analytically and numerically in the specific case of a system of two coupled oscillators featuring a 1:2 internal resonance. More specifically, the locus of Neimark-Sacker points is analytically derived and its evolution with respect to the system parameters is highlighted. The backbone curves, solution of the conservative system, are first investigated, showing in particular the existence of two families of periodic orbits, denoted as parabolic modes. The behaviour of these modes, when the detuning between the eigenfrequencies of the system is varied, is underlined. The non-vanishing limit value, at the origin of one solution family, allows explaining the appearance of isolated solutions for the damped-forced system. The results are then applied to a Micro-Electro-Mechanical System-like shallow arch structure, to show how the analytical expression of the Neimark-Sacker boundary curve can be used for rapid prediction of the appearance of quasiperiodic regime, and thus frequency combs, in Micro-Electro-Mechanical System dynamics.

Abstract This paper is devoted to a detailed analysis of the appearance of frequency combs in the dynamics of a micro-electro-mechanical systems (MEMS) resonator featuring 1:2 internal resonance. To that purpose, both experiments and numerical predictions are reported and analysed to predict and follow the appearance of the phononic frequency comb arising as a quasi-periodic regime between two Neimark-Sacker bifurcations. Numerical predictions are based on a reduced-order model built thanks to an implicit condensation method, where both mechanical nonlinearities and electrostatic forces are taken into account. The reduced order model is able to predict a priori, i.e. without the need of experimental calibration of parameters, and in real time, i.e. by solving one or two degrees-of-freedom system of equations, the nonlinear behaviour of the MEMS resonator. Numerical predictions show a good agreement with experiments under different operating conditions, thus proving the great potentiality of the proposed simulation tool. In particular, the bifurcation points and frequency content of the frequency comb are carefully predicted by the model, and the main features of the periodic and quasi-periodic regimes are given with accuracy, underlining that the complex dynamics of such MEMS device is effectively driven by the characteristics of the 1:2 internal resonance.

The aim of this contribution is to present numerical comparisons of model-order reduction methods for geometrically nonlinear structures in the general framework of finite element (FE) procedures. Three different methods are compared: the implicit condensation and expansion (ICE), the quadratic manifold computed from modal derivatives (MD), and the direct normal form (DNF) procedure, the latter expressing the reduced dynamics in an invariant-based span of the phase space. The methods are first presented in order to underline their common points and differences, highlighting in particular that ICE and MD use reduction subspaces that are not invariant. A simple analytical example is then used in order to analyze how the different treatments of quadratic nonlinearities by the three methods can affect the predictions. Finally, three beam examples are used to emphasize the ability of the methods to handle curvature (on a curved beam), 1:1 internal resonance (on a clamped-clamped beam with two polarizations), and inertia nonlinearity (on a cantilever beam).

In the context of turbulent flows laden with inertial particles, the accurate estimation of preferential concentration is particularly relevant. We have recently proposed to use Voronoï diagrams to estimate concentration fields from two-dimensional imaging techniques implemented around wind tunnel experiments. Due to various experimental biases, the relevance of such an analysis becomes questionable. In this paper, we show the robustness of the Voronoï analysis with respect to the three more important identified biases possibly present in such experiments.

We propose a novel nonlinear manifold learning from snapshot data and demonstrate its superiority over proper orthogonal decomposition (POD) for shedding-dominated shear flows. Key enablers are isometric feature mapping, Isomap, as encoder and, $K$ -nearest neighbours ( $K$ NN) algorithm as decoder. The proposed technique is applied to numerical and experimental datasets including the fluidic pinball, a swirling jet and the wake behind a couple of tandem cylinders. Analysing the fluidic pinball, the manifold is able to describe the pitchfork bifurcation and the chaotic regime with only three feature coordinates. These coordinates are linked to the vortex-shedding phases and the force coefficients. The manifold coordinates of the swirling jet are comparable to the POD mode amplitudes, yet allow for a more distinct and less noise-sensitive manifold identification. A similar observation is made for the wake of two tandem cylinders. The tandem cylinders are aligned and located at a streamwise distance which corresponds to the transition between the single bluff body and the reattachment regimes of vortex shedding. Isomap unveils these two shedding regimes while the Lissajous plot of the first two POD mode amplitudes features a single circle. The reconstruction error of the manifold model is small compared with the fluctuation level, indicating that the low embedding dimensions contain the coherent structure dynamics. The proposed Isomap– $K$ NN manifold learner is expected to be of great importance in estimation, dynamic modelling and control for a large range of configurations with dominant coherent structures.

Micro-Electro-Mechanical Systems revolutionized the consumer market for their small dimensions, high performances and low costs. In recent years, the evolution of the Internet of Things is posing new challenges to MEMS designers that have to deal with complex multiphysics systems experiencing highly nonlinear dynamic responses. To be able to simulate a priori and in real-time the behavior of such systems it is thus becoming mandatory to understand the sources of nonlinearities and avoid them when harmful or exploit them for the design of innovative devices. In this work, we present the first numerical tool able to estimate a priori and in real-time the complex nonlinear responses of MEMS devices without resorting to simplified theories. Moreover, the proposed tool predicts different working conditions without the need of ad-hoc calibration procedures. It consists in a nonlinear Model Order Reduction Technique based on the Implicit Static Condensation that allows to condense the high fidelity FEM models into few degrees of freedom, thus greatly speeding-up the solution phase and improving the design process of MEMS devices. In particular, the 1:2 internal resonance experienced in a MEMS gyroscope test-structure fabricated with a commercial process is numerically investigated and an excellent agreement with experiments is found.

Summary The objective of studies presented in this paper is to analyse the spatial incoherency of seismic ground motions using signals from a velocimeter dense array located on a rock site, recording the aftershock sequence of the two M6 Kefalonia earthquakes that occurred in January/February 2014 (Kefalonia island, Greece). The analyses are carried out with both horizontal and vertical components of velocigrams for small separation distances of stations (<100 m). The coherencies of seismic ground motions identified from strong motion windows are compared with those identified from coda parts of signals. It is realized that there is no significant difference between the coherencies estimated from those two parts of signals. The influence of earthquake event number on the result of coherencies and the dispersions of coherencies estimated from different earthquake events are presented. Finally, coherencies estimated from the dense array are compared with several coherency models proposed and widely used in the literature. The possibility of modifying some parameters of those existing coherency models to fit with in situ coherencies are discussed and presented. Copyright © 2017 John Wiley & Sons, Ltd.

Hippocampal place cells and entorhinal grid cells are thought to form a representation of space by integrating internal and external sensory cues. Experimental data show that different subsets of place cells are controlled by vision, self-motion or a combination of both. Moreover, recent studies in environments with a high degree of visual aliasing suggest that a continuous interaction between place cells and grid cells can result in a deformation of hexagonal grids or in a progressive loss of visual cue control over grid fields. The computational nature of such a bidirectional interaction remains unclear. In this work we present a neural network model of the dynamic interaction between place cells and grid cells within the entorhinal-hippocampal processing loop. The model was tested in two recent experimental paradigms involving environments with visually similar compartments that provided conflicting evidence about visual cue control over self-motion-based spatial codes. Analysis of the model behavior suggests that the strength of entorhinal-hippocampal dynamical loop is the key parameter governing differential cue control in multi-compartment environments. Moreover, construction of separate spatial representations of visually identical compartments required a progressive weakening of visual cue control over place fields in favor of self-motion based mechanisms. More generally our results suggest a functional segregation between plastic and dynamic processes in hippocampal processing.

We propose a self-supervised cluster-based hierarchical reduced-order modelling methodology to model and analyse the complex dynamics arising from a sequence of bifurcations for a two-dimensional incompressible flow of the fluidic pinball. The hierarchy is guided by a triple decomposition separating a slowly varying base flow, dominant shedding and secondary flow structures. All these flow components are kinematically resolved by a hierarchy of clusters. The transition dynamics between these clusters is described by a directed network, called the cluster-based hierarchical network model (HiCNM). Three consecutive Reynolds number regimes for different dynamics are considered: (i) periodic shedding at $Re=80$ , (ii) quasi-periodic shedding at $Re=105$ and (iii) chaotic shedding at $Re=130$ , involving three unstable fixed points, three limit cycles, two quasi-periodic attractors and a chaotic attractor. The HiCNM enables identification of the dynamics between multiple invariant sets in a self-supervised manner. Both the global trends and the local structures during the transition are well resolved by a moderate number of hierarchical clusters. The proposed HiCNM provides a visual representation of transient and post-transient, multi-frequency, multi-attractor behaviour and may automate the identification and analysis of complex dynamics with multiple scales and multiple invariant sets.

Purpose – The purpose of this paper is to present a new bond slip model for reinforced concrete structures. It consists in an interface element (3D) which represents the interface between concrete (modeled in 3D) and steel, modeled using 1D truss elements. Design/methodology/approach – The formulation of the interface element is presented and verified through a comparison with an analytical solution on an academic case. Finally, the model is compared with experimental results on a reinforced concrete tie. Findings – Contrary to the classical perfect or “no-slip” relation which supposes the same displacement between steel and concrete, the proposed model is able to reproduce both global (force-displacement curve) and local (crack openings) results. Originality/value – The proposed approach, applicable to large-scale computations, represents a valuable alternative to the no-slip relation hypothesis to correctly capture the crack properties of reinforced concrete structures.

The direct parametrisation method for invariant manifolds is a nonlinear reduction technique which derives nonlinear mappings and reduced-order dynamics that describe the evolution of dynamical systems along a low-dimensional invariant-based span of the phase space. It can be directly applied to finite element problems. When the development is performed using an arbitrary order asymptotic expansion, it provides an efficient reduced-order modeling strategy for geometrically nonlinear structures. It is here applied to the case of rotating structures featuring centrifugal effect. A rotating cantilever beam with large amplitude vibrations is first selected in order to highlight the main features of the method. Numerical results show that the method provides accurate reduced-order models (ROMs) for any rotation speed and vibration amplitude of interest with a single master mode, thus offering remarkable reduction in the computational burden. The hardening/softening transition of the fundamental flexural mode with increasing rotation speed is then investigated in detail and a ROM parametrised with respect to rotation speed and forcing frequencies is detailed. The method is then applied to a twisted plate model representative of a fan blade, showing how the technique can handle more complex structures. Hardening/softening transition is also investigated as well as interpolation of ROMs, highlighting the efficacy of the method.

This paper presents a novel derivation of the direct parametrisation method for invariant manifolds able to build simulation-free reduced-order models for nonlinear piezoelectric structures, with a particular emphasis on applications to Micro Electro Mechanical Systems. The constitutive model adopted accounts for the hysteretic and electrostrictive response of the piezoelectric material by resorting to the Landau-Devonshire theory of ferroelectrics. Results are validated with full-order simulations operated with a harmonic balance finite element method to highlight the reliability of the proposed reduction procedure. Numerical results show a remarkable gain in terms of computing time as a result of the dimensionality reduction process over low dimensional invariant sets. Results are also compared with experimental data to highlight the remarkable benefits of the proposed model order reduction technique.

We devise and evaluate numerically Hybrid High-Order (HHO) methods for hyperelastic materials undergoing finite deformations. The HHO methods use as discrete unknowns piecewise polynomials of order $$k\ge 1$$ on the mesh skeleton, together with cell-based polynomials that can be eliminated locally by static condensation. The discrete problem is written as the minimization of a broken nonlinear elastic energy where a local reconstruction of the displacement gradient is used. Two HHO methods are considered: a stabilized method where the gradient is reconstructed as a tensor-valued polynomial of order k and a stabilization is added to the discrete energy functional, and an unstabilized method which reconstructs a stable higher-order gradient and circumvents the need for stabilization. Both methods satisfy the principle of virtual work locally with equilibrated tractions. We present a numerical study of the two HHO methods on test cases with known solution and on more challenging three-dimensional test cases including finite deformations with strong shear layers and cavitating voids. We assess the computational efficiency of both methods, and we compare our results to those obtained with an industrial software using conforming finite elements and to results from the literature. The two HHO methods exhibit robust behavior in the quasi-incompressible regime.