Search

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).

International audience; A system composed of two cubic nonlinear oscillators with close natural frequencies, and thus displaying a 1:1 internal resonance, is studied both theoretically and experimentally, with a special emphasis on the free oscillations and the backbone curves. The instability regions of uncoupled solutions are derived and the bifurcation scenario as a function of the parameters of the problem is established, showing in an exhaustive manner all possible solutions. The backbone curves are then experimentally measured on a circular plate, where the asymmetric modes are known to display companion configurations with close eigenfrequencies. A control system based on a Phase-Locked Loop (PLL) is used to measure the backbone curves and also the frequency response function in the forced and damped case, including unstable branches. The model is used for a complete identification of the unknown parameters and an excellent comparison is drawn out between theoretical prediction and measurements.

International audience; A unimorph piezoelectric cantilever equipped with an Acoustic Black Hole (ABH) termination is designed for broadband energy harvesting. The ABH termination, with its tapered region, induces a focusing of the flexural vibrations which can be used to increase the efficiency of an energy harvesting device. A modal-based analytical model is presented, providing an explicit form of the electro-mechanical coupling for each beam eigenmode. Closed-form expressions for the coupled mechanical response and electrical outputs are obtained, allowing one to draw out a complete parametric study to optimize the device. The optimization procedure is conducted following two steps: first, optimal location and dimensions of a single piezoelectric patch are achieved by maximizing the modal electro-mechanical coupling factor (MEMCF) for each structural mode. Thanks to the proposed analytical approach, it is clearly shown that by putting the piezoelectric patch at the maximum of the strain field in the tapered termination, and by adjusting its length in accordance with the focalization created by the ABH effect, the ABH cantilever produces much higher MEMCFs over a wide frequency range and thus outperforms those of a uniform beam. Second, optimization of the shunted circuit is comprehensively performed for a circuit with only resistance, or both resistance and inductance, in series or in parallel. Analytical results show that the key design rule resides in matching the time scale of the circuit with that of the forcing frequency. Addition of the inductance allows enhancing the performance, but on a narrow frequency band. Finally, broadband advantages can be further obtained by considering multiple piezoelectric patches, in which the optimum is obtained when the shunted circuit in each patch is tuned targeting an eigenmode of the ABH beam.

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.

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., Comment: 64 pages, 8 figures

International audience; An experimental demonstration of the broadband passive damping capacity of a vibro-impact acoustic black hole (VI-ABH) is reported. A VI-ABH is an adaptation of the classical ABH design consisting of a beam with a tapered edge of decreasing thickness creating an acoustic black hole (ABH), complemented by contact points on which the beam impacts during its vibration. The contact nonlinearity creates a rapid and efficient transfer of vibrational energy from the low-frequency range, where the ABH is known to be ineffective, to the high-frequency range, thus improving the global passive vibration mitigation characteristics. The optimal design of a VI-ABH follows the rule of locating the contact points at local maxima of the low-frequency modes. Experiments clearly demonstrate the gain in performance, both in forced and free vibrations.

International audience; Nonlinear vibrations of free-edge shallow spherical shells with large amplitudes are investigated, with the aim of predicting the type of nonlinearity (hardening/softening behaviour) for each mode of the shell, as a function of the radius R of curvature of the shell, from the plate case (R → ∞) to the limit of non-shallow shell. Two different models (based on von Kármán's assumptions or on full numerical finite element approach), and two different methods (normal form and modal derivatives) are contrasted.

Challenging applications and the increasing request of higher performances of MEMS devices are focusing the attention of the MEMS community on the understanding and modelling of complex dynamic phenomena. In this work, we experimentally detect a mechanical frequency comb arising in a quad-mass MEMS structure electrostatically actuated through comb fingers electrodes and we propose a simplified model able to predict such nonlinear dynamic behaviour. In particular, we hypothesize a 1:3 internal resonance between two in-plane modes of the structure and exhibit the presence of a Neimark-Sacker bifurcation responsible of frequency comb phenomena.

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., Comment: 34 pages, 10 figures

Non-intrusive methods have been used since two decades to derive reduced-order models for geometrically nonlinear structures, with a particular emphasis on the so-called STiffness Evaluation Procedure (STEP), relying on the static application of prescribed displacements in a finite-element context. We show that a particularly slow convergence of the modal expansion is observed when applying the method with 3D elements, because of nonlinear couplings occurring with very high frequency modes involving 3D thickness deformations. Focusing on the case of flat structures, we first show by computing all the modes of the structure that a converged solution can be exhibited by using either static condensation or normal form theory. We then show that static modal derivatives provide the same solution with fewer calculations. Finally, we propose a modified STEP, where the prescribed displacements are imposed solely on specific degrees of freedom of the structure, and show that this adjustment also provides efficiently a converged solution., Comment: 6 tables, 14 figures, 27 pages

A comparison between two methods to derive reduced-order models (ROM) for geometrically nonlinear structures is proposed. The implicit condensation and expansion (ICE) method relies on a series of applied static loadings. From this set, a stress manifold is constructed for building the ROM. On the other hand, nonlinear normal modes rely on invariant manifold theory in order to keep the key property of invariance for the reduced subspaces. When the model coefficients are fully known, the ICE method reduces to a static condensation. However, in the framework of finite element discretization, getting all these coefficients is generally too computationally expensive. The stress manifold is shown to tend to the invariant manifold only when a slow/fast decomposition between master and slave coordinates can be assumed. Another key problem in using the ICE method is related to the fitting procedure when a large number of modes need to be taken into account. A simplified procedure, relying on normal form theory and identification of only resonant monomial terms in the nonlinear stiffness, is proposed and contrasted with the current method. All the findings are illustrated on beams and plates examples.

The direct computation of the third-order normal form for a geometrically nonlinear structure discretised with the finite element (FE) method, is detailed. The procedure allows to define a nonlinear mapping in order to derive accurate reduced-order models (ROM) relying on invariant manifold theory. The proposed reduction strategy is direct and simulation free, in the sense that it allows to pass from physical coordinates (FE nodes) to normal coordinates, describing the dynamics in an invariant-based span of the phase space. The number of master modes for the ROM is not a priori limited since a complete change of coordinate is proposed. The underlying theory ensures the quality of the predictions thanks to the invariance property of the reduced subspace, together with their curvatures in phase space that accounts for the nonresonant nonlinear couplings. The method is applied to a beam discretised with 3D elements and shows its ability in recovering internal resonance at high energy. Then a fan blade model is investigated and the correct prediction given by the ROMs are assessed and discussed. A method is proposed to approximate an aggregate value for the damping, that takes into account the damping coefficients of all the slave modes, and also using the Rayleigh damping model as input. Frequency-response curves for the beam and the blades are then exhibited, showing the accuracy of the proposed method., 34 pages, 10 figures, 2 tables, submitted to CMAME

The objective of this contribution is to compare two methods proposed recently in order to build efficient reduced-order models for geometrically nonlinear structures. The first method relies on the normal form theory that allows one to obtain a nonlinear change of coordinates for expressing the reduced-order dynamics in an invariant-based span of the phase space. The second method is the modal derivative (MD) approach, and more specifically the quadratic manifold defined in order to derive a second-order nonlinear change of coordinates. Both methods share a common point of view, willing to introduce a nonlinear mapping to better define a reduced order model that could take more properly into account the nonlinear restoring forces. However the calculation methods are different and the quadratic manifold approach has not the invariance property embedded in its definition. Modal derivatives and static modal derivatives are investigated, and their distinctive features in the treatment of the quadratic nonlinearity is underlined. Assuming a slow/fast decomposition allows understanding how the three methods tend to share equivalent properties. While they give proper estimations for flat symmetric structures having a specific shape of nonlinearities and a clear slow/fast decomposition between flexural and in-plane modes, the treatment of the quadratic nonlinearity makes the predictions different in the case of curved structures such as arches and shells. In the more general case, normal form approach appears preferable since it allows correct predictions of a number of important nonlinear features, including for example the hardening/softening behaviour, whatever the relationships between slave and master coordinates are., 43 pages, 12 figures, published in Nonlinear Dynamics

International audience; Nonlinear flexural vibrations of slender beams holding both an Acoustic Black Hole termination and a contact non-linearity are numerically studied. The Acoustic Black Hole (ABH) effect is a passive vibration mitigation technique, which has shown attractive properties above a given cut-on frequency. In this contribution, a vibro-impact acoustic black hole (VI-ABH) is introduced, the contact nonlinearity being used as a mean to transfer energy from low to high frequencies. A numerical model of a VI-ABH is derived from an Euler-Bernoulli beam. The contact law is handled with a penalization approach, the visco-elastic layer with a Ross-Kerwin-Ungard model and the problem is solved with a modal approach combined with an energy-conserving time integration scheme. Numerical results show that the VI-ABH brings about important modifications, and changes the nature of more traditional black holes, by redistributing all the vibrational energy. It can lead to a strong decrease of the resonance magnitude at low frequencies. Under steady state noise excitation, parametric studies are realised in the cases of a single contact, a grid of contacts and bilateral contacts layouts, in order to find some optimal designs. Transient dynamics is also studied through the analysis of displacement signal envelope and energy decay time. All the numerical results constantly show that the combination the ABH effect and an energy transfer provided by contact nonlinearity leads to very attractive mitigation template including low frequencies.

International audience; Acoustic Black Hole effect (ABH) describes a passive vibration mitigation technique used to damp out theflexural vibrations of thin structures such as beams. It relies on a wedge profile where the thickness is gradually decreasing, sothat wave velocities slows down, potentially decaying to zero without reflection. The device generally shows very interestingdamping properties in the mid and high-frequency range, however its efficiency in the low-frequency range still remainslimited. In this contribution, the inclusion of a nonlinearity as a way to improve the low-frequency performance of an ABHis investigated. Indeed, the nonlinearity may be used as a vector to transfer energy from the low to the high frequencyrange, where the damping properties are much more significant. Two different mechanisms are tested, based respectively ongeometric and contact nonlinearity.

International audience; Geometrically nonlinear vibrations of thin plates and shells with variable thickness are investigated numerically with the purpose of synthesizing the sound of cymbals. In cymbal making, taper refers to the gradual change in thickness from the centre to the rim, and is known to be a key feature that determines the tone of the instrument. It is generally used in conjunction with shape variations in order to enable the cymbal to play a bell-like sound when hit near its centre, or a crash sound when struck close to the edge. The von Kármán equations for thin plates with thickness and shape variations are derived, and a numerical method combining a Rayleigh-Ritz approach together with a Störmer-Verlet scheme for advancing the problem in time, is detailed. One main advantage of the method is its ability to implement easily any frequency-dependent loss mechanism which is a key property for sound synthesis. Also, the accuracy of the computation of the nonlinear restoring force is especially preserved. The method is employed to synthesize the sounds of cymbal-like instruments. The impact of taper is addressed and the relative effects of both thickness and shape variations, are contrasted.

International audience; A nonlinear magnetic vibration absorber is presented and used to control vibration of a three-storey structure. A distinctive feature of the absorber concerns its versatility for tuning the linear and nonlinear stiffness coefficients, depending on simple geometric design parameters such as the distance between fixed magnets and the moving one. In particular, the absorber can be set either as a nonlinear tuned vibration absorber, a nonlinear energy sink, or a bistable tuned vibration absorber, according to whether the linear stiffness term is positive, vanishing, or negative. The response of the primary structure and the vibration mitigation are investigated in the cases of impulsive shock, free vibration with imposed initial displacement, and single frequency excitation. Significant reductions of the primary structure vibrations are obtained for the three cases investigated, showing the ability of using a vibration absorber only relying on magnetic forces for passive control. The detailed comparisons of the absorbers performance show that, in this case study, no general guidelines can be easily deduced for selecting one of the three tunings for a nonlinear absorber. Depending on the excitation, the vibratory levels, and the frequency content of the excitation, the three configurations show advantages and drawbacks that are discussed.

International audience; Collisions in musical string instruments play a fundamental role in explaining the sound production in various instruments such as sitars, tanpuras and electric basses. Contacts occuring during the vibration provide a nonlinear effect which shapes a specific tone due to energy transfers and enriches the hearing experience. As such, they must be carefully simulated for the purpose of physically-based sound synthesis. Most of the numerical methods presented in the literature rely on a compliant modeling of the contact force between the string and the obstacle. In this contribution, numerical methods from nonsmooth contact dynamics are used to integrate the problem in time. A Moreau-Jean time-stepping scheme is combined with an exact scheme for phases with no contact, thus controlling the numerical dispersion. Results for a two-point bridge mimicking a tanpura and an electric bass are presented, showing the ability of the method to deal efficiently with such problems while invoking, as compared to a compliant approach, less modelling parameters and a reduced computational burden.

International audience; For particular playing techniques such as "pop" or "slap" in the electric bass guitar, the string collides with frets, producing a percussive sound used in dierent music styles. The string/frets contacts introduce a non-linearity which is investigated both numerically and experimentally in this paper. A physical model, based on a modal description of the string, is implemented with an unconditionally stable scheme. Simulations including a string/structure coupling and the two polarisations of the string are confronted to controlled experiments, showing a good agreement for increasing amplitudes of initial conditions. A parametric study is then conducted numerically in order to highlight the inuence of physical parameters on the transient behaviour and raises questions related to tuning and playing issues.

International audience; This article is concerned with the vibration of a stiff linear string in the presence of a rigid obstacle. A numerical method for unilateral and arbitrary-shaped obstacles is developed, based on a modal approach in order to take into account the frequency dependence of losses in strings. The contact force of the barrier interaction is treated using a penalty approach, while a conservative scheme is derived for time integration, in order to ensure long-term numerical stability. In this way, the linear behaviour of the string when not in contact with the barrier can be controlled via a mode by mode fitting, so that the model is particularly well suited for comparisons with experiments. An experimental configuration is used with a point obstacle either centered or near an extremity of the string. In this latter case, such a pointwise obstruction approximates the end condition found in the tanpura, an Indian stringed instrument. The second polarisation of the string is also analysed and included in the model. Numerical results are compared against experiments, showing good accuracy over a long time scale.

International audience; Acoustic Black Hole effect (ABH) is a passive vibration damping technique without added mass based on flexural waves properties in thin structures with variable thickness. A common implementation is a plate edge where the thickness is locally reduced with a power law profile and covered with a viscoelastic layer. The plate displacement in the small thickness region is large and easily exceeds the plate thickness. This is the origin of geometric nonlinearity which can generate couplings between linear eigenmodes of the structure and induce energy transfer between low and high frequency regimes. This phenomenon may be used to increase the efficiency of the ABH treatment in the low frequency regime where it is usually inefficient. An experimental investigation evidenced that usual ABH implementation gives rise to measurable geometric nonlinearity and typical nonlinear phenomena. In particular, strongly nonlinear regime and wave turbulence are reported. The nonlinear ABH beam is then modeled as a von Kármán plate with variable thickness. The model is solved numerically by using a modal method combined with an energy-conserving time integration scheme. The effects of both the thickness profile and the damping layer are then investigated in order to improve the damping properties of an ABH beam. It is found that a compromise between the two effects can lead to an important gain of efficiency in the low frequency range.

International audience; The influence of a hysteretic damper on the airfoil flutter instability is investigated. In particular, its effect on the post-critical limit cycle oscillations (LCOs) is emphasized. For that purpose, an aeroelastic model including large amplitude motions and dynamic stall phenomenon, is considered for a rigid flat plate having two degrees of freedom in pitch and plunge motions. The hysteretic behaviour is modeled thanks to a generalized Bouc-Wen formulation. A parametric study of aeroelastic as well as hysteresis model parameters, allows one to draw a complete picture of the bifurcation scenario, highlighting the capacity of the hysteretic damper in precluding the occurrence of stall. The special case of shape memory alloy (SMA) springs is then used numerically and experimentally for controlling the flutter oscillations of a flat plate. The study reveals the ability of the SMA springs to drastically reduce the amplitudes of the LCOs caused by the flutter instability.

A phenomenological model describing the time-frequency dependence of the power spectrum of thin plates vibrating in a wave turbulence regime, is introduced. The model equation contains as basic solutions the Rayleigh-Jeans equipartition of energy, as well as the Kolmogorov-Zakharov spectrum of wave turbulence. In the Wave Turbulence Theory framework, the model is used to investigate the self-similar, non-stationary solutions of forced and free turbulent vibrations. Frequency-dependent damping laws can easily be accounted for. Their effects on the characteristics of the stationary spectra of turbulence are then investigated. Thanks to this analysis, self-similar universal solutions are given, relating the power spectrum to both the injected power and the damping law., Comment: 22 pages, 9 figures

International audience; The present work studies the double polarisation phenomenon observed in vibrating piano strings. From the experimental viewpoint , it is known that when a string is given an initial displacement in one transverse direction (e.g. hammer excitation in the vertical plane), the second transverse displacement (e.g. in the horizontal plane) is also excited after a few milliseconds and the amplitude can be of similar order to the first transverse displacement. This phenomenon contributes to a characteristic piano sound feature called the "double decay". The purpose of this study is to investigate the role of nonlinearities in inducing double polarisations. The nonlinear vibrations of the strings are studied with a two-degrees-of-freedom (dofs) system extracted from the Kirchhoff-Carrier string equations. The method of multiple scales is used to study the free vibrations of two po-larisations having nearly equal eigenfrequencies and thus presenting a 1:1 internal resonance. For an imperfect string with slightly different eigenfrequencies between the two polarisa-tions, it is found out that depending on the energy of the excita-tion, an uncoupled transverse mode can develop into a coupled mode where there is energy exchange between the two transverse polarisations. The coupled mode is stable and the string oscillates in an elliptic path. Numerical experiments are also carried out, confirming the findings of the analytical approach.

The authors are grateful to Bertrand David (Telecom-ParisTech) for computing the code allowing the STFT filtering procedure used in Section 5.1. The filter has been designed in the framework of the PAFI project (Plateforme d’Aide la facture Instrumentale, www.pafi.fr) which is also thanked.; The vibrations and sounds produced by two notes of a double second steelpan are investigated, the main objective being to quantify the nonlinear energy exchanges occurring between vibration modes that are responsible of the peculiar sound of the instrument. A modal analysis first reveals the particular tuning of the modes and the systematic occurence of degenerate modes, from the second one, this feature being a consequence of the tuning and the mode localization. Forced vibrations experiments are then performed to follow precisely the energy exchange between harmonics of the vibration and thus quantify properly the mode couplings. In particular, it is found that energy exchanges are numerous, resulting in complicated frequency response curves even for very small levels of vibration amplitude. Simple models displaying1:2:2 and 1:2:4 internal resonance are then fitted to the measurements, allowing to identify the values of the nonlinear quadratic coupling coefficients resulting from the geometric nonlinearity. The identified 1:2:4 model is finally used to recover the time domain variations of an impacted note in normal playing condition, resulting in an excellent agreement for the temporal behaviour of the first four harmonics.; International audience; The vibrations and sounds produced by two notes of a double second steelpan are investigated, the main objective being to quantify the nonlinear energy exchanges occurring between vibration modes that are responsible of the peculiar sound of the instrument. A modal analysis first reveals the particular tuning of the modes and the systematic occurence of degenerate modes, from the second one, this feature being a consequence of the tuning and the mode localization. Forced vibrations experiments are then performed to follow precisely the energy exchange between harmonics of the vibration and thus quantify properly the mode couplings. In particular, it is found that energy exchanges are numerous, resulting in complicated frequency response curves even for very small levels of vibration amplitude. Simple models displaying1:2:2 and 1:2:4 internal resonance are then fitted to the measurements, allowing to identify the values of the nonlinear quadratic coupling coefficients resulting from the geometric nonlinearity. The identified 1:2:4 model is finally used to recover the time domain variations of an impacted note in normal playing condition, resulting in an excellent agreement for the temporal behaviour of the first four harmonics.

International audience; This paper presents a modal, time-domain scheme for the nonlinear vibrations of perfect and imperfect plates. The scheme can take into account a large number of degrees-of-freedom and is energy-conserving. The targeted application is the sound synthesis of cymbals and gong-like musical instruments, which are known for displaying a strongly nonlinear vibrating behaviour. This behaviour is typical of a wave turbulence regime, in which the wide-band spectrum of excited modes is observable in the form of an energy cascade. The modal method is selected for its versatility in handling complex damping laws that can be implemented easily by selecting appropriate damping values in each one of the modal equations. In the first part of the paper, the modal method is explained in its generality, and it will be seen that the method is valid for plates with arbitrary geometry and boundary conditions as long as the eigenmodes are known. Secondly, a time-integration, energy-conserving scheme for perfect and imperfect plates is presented, and implementation comments are given in order to treat efficiently the high-dimensionality of the resulting dynamical system. The scheme is run with appropriate parameters in order to produce sound samples. A simple impact law is considered for the exci-tation, whereas the flexibility of the method is highlighted by showing simulations for free-edge circular plates and simply-supported rectangular plates, together with various damping laws.

International audience; A passive control of aeroelastic instabilities on a two-degrees-of-freedom (dofs) system is considered here using shape memory alloys (SMA) springs in their pseudo-elastic regime. SMA present a solid-solid phase change that allow them to face strong deformations (∼10%) ; in the pseudo-elastic regime, an hysteresis loop appears in the stress-strain relationship which in turn gives rise to an important amount of dissipated energy. This property makes the SMA a natural candidate for mitigating undesired vibrations in a passive manner. A 2-dofs system is here used to model the classical flutter instability of a wing section in a uniform flow. The SMA spring is selected to act on the pitch in order to dissipate energy of the predominant motion. A simple phenomenological model for the SMA hysteresis loop is introduced, allowing for a quantitative study of the important parameters to optimize in view of an experimental design. Thanks to a simple phenomenological model for the SMA hysteresis loop, a quantitative numerical study is performed in order to exhibit the best tuning of the material parameters for controlling the flutter instability.

International audience; A passive control of aeroelastic instabilities on a two-degrees-of-freedom (dofs) system is considered here using shape memory alloys (SMA) springs in their pseudo-elastic regime. SMA present a solid-solid phase change that allow them to face strong deformations (∼ 10%); in the pseudo-elastic regime, an hysteresis loop appears in the stress-strain relationship which in turn gives rise to an important amount of dissipated energy. This property makes the SMA a natural candidate for damping undesired vibrations in a passive manner. A 2-dofs system is here used to model the classical flutter instability of a wing section in an uniform flow. The SMA spring is selected on the pitch mode in order to dissipate energy of the predominant motion. A simple model for the SMA hysteresis loop is introduced, allowing for a quantitative study of the important parameters to optimize in view of an experimental design.

International audience; We propose a direct method for computing modal coupling coefficients - due to geometrically nonlinear effects - for thin shells vibrating at large amplitude and discretized by a finite element (FE) procedure. These coupling coefficients arise when considering a discrete expansion of the unknown displacement onto the eigenmodes of the linear operator. The evolution problem is thus projected onto the eigenmodes basis and expressed as an assembly of oscillators with quadratic and cubic nonlinearities. The nonlinear coupling coefficients are directly derived from the finite element formulation, with specificities pertaining to the shell elements considered, namely, here elements of the ''Mixed Interpolation of Tensorial Components'' family (MITC). Therefore, the computation of coupling coefficients, combined with an adequate selection of the significant eigenmodes, allows the derivation of effective reduced-order models for computing - with a continuation procedure - the stable and unstable vibratory states of any vibrating shell, up to large amplitudes. The procedure is illustrated on a hyperbolic paraboloid panel. Bifurcation diagrams in free and forced vibrations are obtained. Comparisons with direct time simulations of the full FE model are given. Finally, the computed coefficients are used for a maximal reduction based on asymptotic nonlinear normal modes (NNMs), and we find that the most important part of the dynamics can be predicted with a single oscillator equation.

These lecture notes are related to the CISM course on ”Modal Analysis of nonlinear Mechanical systems”, held at Udine, Italy, from June 25 to 29, 2012. The key concept at the core of all the lessons given during this week is the notion of Nonlinear Normal Mode (NNM), a theoretical tool allowing one to extend, through some well-chosen assumptions and limitations, the linear modes of vibratory systems, to nonlinear regimes. More precisely concerning these notes, they are intended to show the explicit link between Normal Form theory and NNMs, for the specific case of vibratory systems displaying polynomial type nonlinearities. After a brief introduction reviewing the main concepts for deriving the normal form for a given dynamical system, the relationship between normal form theory and nonlinear normal modes (NNMs) will be the core of the developments. Once the main results presented, application of NNMs to vibration problem where geometric nonlinearity is present, will be highlighted. In particular, the developments of reduced-order models based on NNMs expressed asymptotically with the formalism of real normal form, will be deeply presented.

Nonlinear vibrations of thin rectangular plates are considered, using the von kármán equations in order to take into account the effect of geometric nonlinearities. Solutions are derived through an expansion over the linear eigenmodes of the system for both the transverse displacement and the Airy stress function, resulting in a series of coupled oscillators with cubic nonlinearities, where the coupling coefficients are functions of the linear eigenmodes. A general strategy for the calculation of these coefficients is outlined, and the particular case of a simply supported plate with movable edges is thoroughly investigated. To this extent, a numerical method based on a new series expansion is derived to compute the Airy stress function modes, for which an analytical solution is not available. It is shown that this strategy allows the calculation of the nonlinear coupling coefficients with arbitrary precision, and several numerical examples are provided. Symmetry properties are derived to speed up the calculation process and to allow a substantial reduction in memory requirements. Finally, analysis by continuation allows an investigation of the nonlinear dynamics of the first two modes both in the free and forced regimes. Modal interactions through internal resonances are highlighted, and their activation in the forced case is discussed, allowing to compare the nonlinear normal modes (NNMs) of the undamped system with the observable periodic orbits of the forced and damped structure. © 2013 Springer-Verlag Wien.

The effect of damping in the wave turbulence regime for thin vibrating plates is studied. An experimental method, allowing measurements of dissipation in the system at all scales, is first introduced. Practical experimental devices for increasing the dissipation are used. The main observable consequence of increasing the damping is a significant modification in the slope of the power spectral density, so that the observed power laws are not in a pure inertial regime. However, the system still displays a turbulent behavior with a cut-off frequency that is determined by the injected power which does not depend on damping. By using the measured damping power-law in numerical simulations, similar conclusions are drawn out., Comment: 6 pages, 6 figures

International audience; Numerical computation of Nonlinear Normal Modes (NNMs) for conservative vibratory systems is addressed, with the aim of deriving accurate reduced-order models up to large amplitudes. A numerical method is developed, based on the center manifold approach for NNMs, which uses an interpretation of the equations as a transport problem, coupled to a periodicity condition for ensuring manifold's continuity. Systematic comparisons are drawn with other numerical methods, and especially with continuation of periodic orbits, taken as reference solutions. Three di erent mechanical systems, displaying peculiar characteristics allowing for a general view of the performance of the methods for vibratory systems, are selected. Numerical results show that invariant manifolds encounter folding points at large amplitude, generically (but not only) due to internal resonances. These folding points involve an intrinsic limitation to reduced-order models based on the center manifold and on the idea of a functional relationship between slave and master coordinates. Below that amplitude limit, numerical methods are able to produce reduced-order models allowing for a precise prediction of the backbone curve.

International audience; The non-linear dynamic thermomechanical behaviour of superelastic shape memory alloys is investigated. To this end, the Zaki-Moumni model, initially developed for quasi-static loading cases, is extended to simulate the uniaxial forced oscillations of a shape memory alloy device. First, the influence of loading rate is accounted for by considering the thermomechanical coupling in the behaviour of NiTi shape memory alloy. Comparisons between simulations and experimental results show good agreement. Then, the forced response of a shape memory alloy device is investigated at resonance. Both isothermal and non-isothermal conditions are studied, as well as non-symmetric tensile-compressive restoring force. In the case of large values of forcing amplitudes, simulation results show that the dynamic response is prone to jumps, bifurcations and chaotic solutions.

International audience; A thin plate, excited by a harmonic external forcing of increasing amplitude, shows transitions from a periodic response to a chaotic state of wave turbulence. By analogy with the transition to turbulence observed in fluid mechanics as the Reynolds number is increased, a generic transition scenario for thin vibrating plates, first experimentally observed, is here numerically studied. The von Kármán equations for thin plates, which include geometric non-linear effects, are used to model large amplitude vibrations, and an energy-conserving finite difference scheme is employed for discretisation. The transition scenario involves two bifurcations separating three distinct regimes. The first regime is the periodic, weakly non-linear response. The second is a quasiperiodic state where energy is exchanged between internally resonant modes. It is observed only when specific internal resonance relationships are fulfilled between the eigenfrequencies of the structure and the forcing frequency; otherwise a direct transition to the last turbulent state is observed. This third, or turbulent, regime is characterized by a broadband Fourier spectrum and a cascade of energy from large to small wavelengths. For perfect plates including cubic non-linearity, only third-order internal resonances are likely to exist. For imperfect plates displaying quadratic nonlinearity, the energy exchanges and the quasiperiodic states are favored and thus are more easily obtained. Finally, the turbulent regime is characterized in the light of available theoretical results from wave turbulence theory.

International audience; Perturbation methods are routinely used in all fields of applied mathematics where analytical solutions for nonlinear dynamical systems are searched. Among them, normal form theory provides a reliable method for systematically simplifying dynamical systems via nonlinear change of coordinates, and is also used in a mechanical context to define Nonlinear Normal Modes (NNMs). The main recognized drawback of perturbation methods is the absence of a criterion establishing their range of validity in terms of amplitude. In this paper, we propose a method to obtain upper bounds for amplitudes of changes of variables in normal form transformations. The criterion is tested on simple mechanical systems with one and two degrees-of-freedom, and for complex as well as real normal form. Its behavior with increasing order in the normal transform is established, and comparisons are drawn between exact solutions and normal form computations for increasing levels of amplitudes. The results clearly establish that the criterion gives an upper bound for validity limit of normal transforms.

International audience; This paper investigates experimentally the quasi-static and dynamic torsional behaviour of shape memory alloys wires under cyclic loading. A specifically designed torsional pendulum made of a Ni–Ti wire is described. Results on the quasi-static behaviour of the wire obtained using this setup are presented, giving an overall view of the damping capacity of the material as function of the amplitude of the loading (imposed torsional angle), the frequency and the temperature. The dynamical behaviour is then presented through measured frequency response function between forcing angle at the top of the pendulum and the difference between top and bottom rotation angles in the vicinity of the first eigenfrequency of the wire, i.e. in the range [0.3 Hz, 1 Hz]. The softening-type non-linearity and its subsequent jump phenomenon, predicted theorically by the decrease of the effective stiffness when martensite transformation starts is clearly evidenced and analysed.

International audience; The transition from periodic to chaotic vibrations in free-edge, perfect and imperfect circular plates, is numerically studied. A pointwise harmonic forcing with constant frequency and increasing amplitude is applied to observe the bifurcation scenario. The von Kármán equations for thin plates, including geometric non-linearity, are used to model the large-amplitude vibrations. A Galerkin approach based on the eigenmodes of the perfect plate allows discretizing the model. The resulting ordinary-differential equations are numerically integrated. Bifurcation diagrams of Poincaré maps, Lyapunov exponents and Fourier spectra analysis reveal the transitions and the energy exchange between modes. The transition to chaotic vibration is studied in the frequency range of the first eigenfrequencies. The complete bifurcation diagram and the critical forces needed to attain the chaotic regime are especially addressed. For perfect plates, it is found that a direct transition from periodic to chaotic vibrations is at hand. For imperfect plates displaying specific internal resonance relationships, the energy is first exchanged between resonant modes before the chaotic regime. Finally, the nature of the chaotic regime, where a high-dimensional chaos is numerically found, is questioned within the framework of wave turbulence. These numerical findings confirm a number of experimental observations made on shells, where the generic route to chaos displays a quasiperiodic regime before the chaotic state, where the modes, sharing internal resonance relationship with the excitation frequency, appear in the response.

International audience; A numerical study of the transition from periodic to chaotic motions in forced vibrations of circular plates, is proposed. A pointwise harmonic forcing of constant excitation frequency Ω and increasing values of the amplitude is considered. Perfect and imperfect circular plates with a free edge are studied within the von Kármán assumptions for large displacements (geometric non-linearity). The transition scenario is observed for different excitation frequencies in the range of the first eigenfrequencies of the plate. For perfect plate with no specific internal resonance relationships , a direct transition to chaos is at hand. For imperfect plate tuned so as to fulfill specific internal resonance relations, a coupling between internally resonant modes is first observed. The chaotic regime shows an attractor of large dimension, and thus is studied within the framework of wave turbulence.

International audience; A vibrating plate is set into a chaotic state of wave turbulence by a forcing having periodic and random components. Both components are weighted in order to explore continuously intermediate forcing from the periodic to the random one, but keeping constant its rms value. The transverse velocity of the plate is measured at the application point of the force. It is found that whatever the detail of the forcing is, the velocity spectra exhibit a universal cascade for frequencies larger than the forcing frequency range. In contrast, the velocity spectra strongly depend on the nature of the forcing within the range of forcing frequencies. The coherence function is used to extract the contribution of the velocity fluctuations that display a linear relationship with the forcing. The nonlinear contribution to the velocity fluctuations is found to be almost constant, about 55% of the total velocity fluctuations whatever the nature of the forcing from random to periodic. On the other hand, the nonlinear contribution to the fluctuations of the injected power depends on the nature of the forcing; it is significantly larger for the periodic forcing (60%) and decreases continuously as the randomness is increased, reaching a value of 40% for the pure random forcing. For all the cases of intermediate forcing from random to periodic, a simple model of the velocity response recovers in a fairly good agreement the probability density function of the injected power. The consequence of the existence of a linear-response component is discussed in the context of the fluctuation-dissipation theorem validation in experiments of out-of-equilibrium systems. © 2010 The American Physical Society.

The nonlinear interaction of waves in a driven medium may lead to wave turbulence, a state such that energy is transferred from large to small length scales. Here, wave turbulence is observed in experiments on a vibrating plate. The frequency power spectra of the normal velocity of the plate may be rescaled on a single curve, with power-law behaviors that are incompatible with the weak turbulence theory of Düring et al. [Phys. Rev. Lett. 97, 025503 (2006)10.1103/PhysRevLett.97.025503]. Alternative scenarios are suggested to account for this discrepancy -- in particular the occurrence of wave breaking at high frequencies. Finally, the statistics of velocity increments do not display an intermittent behavior.

International audience; Non-linear normal modes (NNMs) are used in order to derive reduced-order models for large amplitude, geometrically non-linear vibrations of thin shells. The main objective of the paper is to compare the accuracy of different truncations, using linear and non-linear modes, in order to predict the response of shells structures subjected to harmonic excitation. For an exhaustive comparison, three different shell problems have been selected: (i) a doubly curved shallow shell, simply supported on a rectangular base; (ii) a circular cylindrical panel with simply supported, in-plane free edges; and (iii) a simply supported, closed circular cylindrical shell. In each case, the models are derived by using refined shell theories for expressing the strain-displacement relationship. As a consequence, in-plane inertia is retained in the formulation. Reduction to one or two NNMs shows perfect results for vibration amplitude lower or equal to the thickness of the shell in the three cases, and this limitation is extended to two times the thickness for two of the selected models. © 2008 Elsevier B.V. All rights reserved.

19 pages; International audience; The effect of geometric imperfections and viscous damping on the type of nonlinearity (i.e., the hardening or softening behaviour) of circular plates and shallow spherical shells with free edge is here investigated. The Von Kármán large-deflection theory is used to derive the continuous models. Then, nonlinear normal modes (NNMs) are used for predicting with accuracy the coefficient, the sign of which determines the hardening or softening behaviour of the structure. The effect of geometric imperfections, unavoidable in real systems, is studied by adding a static initial component in the deflection of a circular plate. Axisymmetric as well as asymmetric imperfections are investigated, and their effect on the type of nonlinearity of the modes of an imperfect plate is documented. Transitions from hardening to softening behaviour are predicted quantitatively for imperfections having the shapes of eigenmodes of a perfect plate. The role of 2:1 internal resonance in this process is underlined. When damping is included in the calculation, it is found that the softening behaviour is generally favoured, but its effect remains limited.

International audience; A vibrating plate is set into a chaotic state of wave turbulence by either a periodic or a random local forcing. Correlations between the forcing and the local velocity response of the plate at the forcing point are studied. Statistical models with fairly good agreement with the experiments are proposed for each forcing. Both distributions of injected power have a logarithmic cusp for zero power, while the tails are Gaussian for the periodic driving and exponential for the random one. The distributions of injected work over long time intervals are investigated in the framework of the fluctuation theorem, also known as the Gallavotti-Cohen theorem. It appears that the conclusions of the theorem are verified only for the periodic, deterministic forcing. Using independent estimates of the phase space contraction, this result is discussed in the light of available theoretical framework.

International audience; The aim of the present paper is to compare two different methods available for reducing the complicated dynamics exhibited by large amplitude, geometrically nonlinear vibrations of a thin shell. The two methods are: the proper orthogonal decomposition (POD), and an asymptotic approximation of the nonlinear normal modes (NNMs) of the system. The structure used to perform comparisons is a water-filled, simply supported circular cylindrical shell subjected to harmonic excitation in the spectral neighbourhood of the fundamental natural frequency. A reference solution is obtained by discretizing the partial differential equations (PDEs) of motion with a Galerkin expansion containing 16 eigenmodes. The POD model is built by using responses computed with the Galerkin model; the NNM model is built by using the discretized equations of motion obtained with the Galerkin method, and taking into account also the transformation of damping terms. Both the POD and NNMs allow to reduce significantly the dimension of the original Galerkin model. The computed nonlinear responses are compared in order to verify the accuracy and the limits of these two methods. For vibration amplitudes equal to 1.5 times the shell thickness, the two methods give very close results to the original Galerkin model. By increasing the excitation and vibration amplitude, significant differences are observed and discussed. The response is investigated also for a fixed excitation frequency by using the excitation amplitude as bifurcation parameter for a wide range of variation. Bifurcation diagrams of Poincaré maps obtained from direct time integration and calculation of the maximum Lyapunov exponent have been used to characterize the system. © 2007 Elsevier Ltd. All rights reserved.

International audience; This study is devoted to the experimental validation of a theoretical model of large amplitude vibrations of thin spherical shells described in a previous study by the same authors. A modal analysis of the structure is first detailed. Then, a specific mode coupling due to a 1:1:2 internal resonance between an axisymmetric mode and two companion asymmetric modes is especially addressed. The structure is forced with a simple-harmonic signal of frequency close to the natural frequency of the axisymmetric mode. The experimental setup, which allows precise measurements of the vibration amplitudes of the three involved modes, is presented. Experimental frequency response curves showing the amplitude of the modes as functions of the driving frequency are compared to the theoretical ones. A good qualitative agreement is obtained with the predictions given by in the model. Some quantitative discrepancies are observed and discussed, and improvements of the model are proposed.