Periodic Motion and Frequency Energy Plots of Dynamical Systems Coupled with Piecewise Nonlinear Energy Sink

The frequency-energy plots (FEPs) of two-degree-of-freedom linear structures attached to a piecewise nonlinear energy sink (PNES) are generated here and thoroughly investigated. This study provides the FEP analysis of such systems for further understanding to nonlinear targeted energy transfer (TET) by the PNES. The attached PNES to the considered linear dynamical systems incorporates a symmetrical clearance zone of zero stiffness content where the boundaries of the zone are coupled with the linear structure by linear stiffness elements. In addition, linear viscous damping is selected to be continuous during the PNES mass oscillation. The underlying nonlinear dynamical behaviour of the considered structure-PNES systems is investigated by generating the fundamental backbone curves of the FEP and the bifurcated subharmonic resonance branches using numerical continuation methods. Accordingly, interesting dynamical behaviour of the nonlinear normal modes (NNMs) of the structure-PNES system on different backbones and subharmonic resonance branches has been observed. In addition, the imposed wavelet transform frequency spectrums on the FEPs have revealed that the TET takes place in multiple resonance captures where it is dominated by the nonlinear action of the PNES.

Stability of Fixed Points in an Approximate Solution of the Spring-Mass Running Model

We consider a classical spring-mass model of human running which is built upon an inverted elastic pendulum. Based on our previous results concerning asymptotic solutions for large spring constant (or small angle of attack), we construct analytical approximations of solutions in the considered model. The model itself consists of two sets of differential equations - one set describes the motion of the centre of mass of a runner in contact with the ground (support phase), and the second set describes the phase of no contact with the ground (flight phase). By appropriately concatenating asymptotic solutions for the two phases we are able to reduce the dynamics to a onedimensional apex to apex return map. We find sufficient conditions for this map to have a unique stable fixed point. By numerical continuation of fixed points with respect to energy, we find a transcritical bifurcation in our model system.MSC 2020 Classification: 34C20, 34D05, 37N25, 70K20, 70K42, 70K50, 70K60

A numerical approach for detecting switch-like bistability in mass action chemical reaction networks with conservation laws

Background Theoretical analysis of signaling pathways can provide a substantial amount of insight into their function. One particular area of research considers signaling pathways capable of assuming two or more stable states given the same amount of signaling ligand. This phenomenon of bistability can give rise to switch-like behavior, a mechanism that governs cellular decision making. Investigation of whether or not a signaling pathway can confer bistability and switch-like behavior, without knowledge of specific kinetic rate constant values, is a mathematically challenging problem. Recently a technique based on optimization has been introduced, which is capable of finding example parameter values that confer switch-like behavior for a given pathway. Although this approach has made it possible to analyze moderately sized pathways, it is limited to reaction networks that presume a uniterminal structure. It is this limited structure we address by developing a general technique that applies to any mass action reaction network with conservation laws. Results In this paper we developed a generalized method for detecting switch-like bistable behavior in any mass action reaction network with conservation laws. The method involves (1) construction of a constrained optimization problem using the determinant of the Jacobian of the underlying rate equations, (2) minimization of the objective function to search for conditions resulting in a zero eigenvalue, (3) computation of a confidence level that describes if the global minimum has been found and (4) evaluation of optimization values, using either numerical continuation or directly simulating the ODE system, to verify that a bistability region exists. The generalized method has been tested on three motifs known to be capable of bistability. Conclusions We have developed a variation of an optimization-based method for the discovery of bistability, which is not limited to uniterminal chemical reaction networks. Successful completion of the method provides an S-shaped bifurcation diagram, which indicates that the network acts as a bistable switch for the given optimization parameters.

Dynamics and bifurcations of a ratio-dependent predator-prey model

In this paper, we study a ratio-dependent predator-prey model with modied Holling-Tanner formalism, by using dynamical techniques and numerical continuation algorithms implemented in Matcont. We determine codim-1 and 2 bifurcation points and their corresponding normal form coecients. We also compute a curve of limit cycles of the system emanating from a Hopf point.

Pressure-driven wrinkling of soft inner-lined tubes

A simple equation modelling an inextensible elastic lining of an inner-lined tube subject to an imposed pressure difference is derived from a consideration of the idealised elastic properties of the lining and the pressure and soft-substrate forces. Two cases are considered in detail, one with prominent wrinkling and a second one in which wrinkling is absent and only buckling remains. Bifurcation diagrams are computed via numerical continuation for both cases. Wrinkling, buckling, folding, and mixed-mode solutions are found and organised according to system-response measures including tension, in-plane compression, maximum curvature and energy. Approximate wrinkle solutions are constructed using weakly nonlinear theory, in excellent agreement with numerics. Our approach explains how the wavelength of the wrinkles is selected as a function of the parameters in compressed wrinkling systems and shows how localised folds and mixed-mode states form in secondary bifurcations from wrinkled states. Our model aims to capture the wrinkling response of arterial endothelium to blood pressure changes but applies much more broadly.

Unified framework for localized patterns in reaction–diffusion systems; the Gray–Scott and Gierer–Meinhardt cases

A recent study of canonical activator-inhibitor Schnakenberg-like models posed on an infinite line is extended to include models, such as Gray–Scott, with bistability of homogeneous equilibria. A homotopy is studied that takes a Schnakenberg-like glycolysis model to the Gray–Scott model. Numerical continuation is used to understand the complete sequence of transitions to two-parameter bifurcation diagrams within the localized pattern parameter regime as the homotopy parameter varies. Several distinct codimension-two bifurcations are discovered including cusp and quadruple zero points for homogeneous steady states, a degenerate heteroclinic connection and a change in connectedness of the homoclinic snaking structure. The analysis is repeated for the Gierer–Meinhardt system, which lies outside the canonical framework. Similar transitions are found under homotopy between bifurcation diagrams for the case where there is a constant feed in the active field, to it being in the inactive field. Wider implications of the results are discussed for other pattern-formation systems arising as models of natural phenomena. This article is part of the theme issue ‘Recent progress and open frontiers in Turing’s theory of morphogenesis’.

Friction-induced traveling wave coupling in tuned bladed-disks

Flutter is a major constraint on modern turbomachines; as the designs move toward more slender, thinner, and loaded blades, they become more prone to experience high cycle fatigue problems. Dry friction, present at the root attachment for cantilever configurations, is one of the main sources of energy dissipation. It saturates the flutter vibration amplitude growth, producing a limit cycle oscillation whose amplitude depends on the balance between the energy injected and dissipated by the system. Both phenomena, flutter and friction, typically produce a small correction of the purely elastic response of the structure. A large number of elastic cycles is required to notice their effect, which appears as a slow modulation of the oscillation amplitude. Furthermore, even longer time scales appear when multiple traveling waves are aerodynamically unstable and exhibit similar growth rates. All these slow scales make the system time integration very stiff and CPU expensive, bringing some doubts about whether the final solutions are properly converged. In order to avoid these uncertainties, a numerical continuation procedure is applied to analyze the solutions that set in, their traveling wave content, their bifurcations and their stability. The system is modeled using an asymptotic reduced order model and the continuation results are validated against direct time integrations. New final states with multiple traveling wave content are found and analyzed. These solutions have not been obtained before for the case of microslip friction at the blade attachment; only solutions consisting of a single traveling wave have been reported in previous works.

Nonlinear stability analysis and numerical continuation of bifurcations of a rotor supported by floating ring bearings

Floating ring bearings have been widely used, over the last decades, in rotors of automotive turbochargers because of their improved damping behavior and their good emergency-operating capabilities. They also offer a cost-effective design and have good assembly properties. Nevertheless, rotors with floating ring bearings show vibration effects of nonlinear nature induced by self-excited oscillations originating from the bearing oil films (oil whirl/whip phenomena) and may exhibit various nonlinear vibration effects which may cause damage to the rotor. In order to investigate these dynamic phenomena, this paper has developed a nonlinear model of a perfectly balanced rigid rotor supported by two identical floating ring bearings with consideration of their vibration behavior mainly governed by fluid dynamics. The dimensionless hydrodynamic forces of floating ring bearings have been derived based on the short bearing theory and the half Sommerfeld solution. Using the numerical continuation approach, different bifurcations are detected when a control parameter, the journal speed, is varied. Depending on the system’s physical parameters, the rotor can show stable or unstable limit cycles which themselves may collapse beyond a certain rotor speed to exhibit a fold bifurcation. Bifurcation analysis is performed to investigate the occurring instabilities and nonlinear phenomena. Such results explain the instabilities characteristics of the floating ring bearing in high-speed applications. It has also been found that the selection of the bearing modulus plays an important role in the characterization of the rotor stability threshold speed and bifurcation sequences. An understanding of the system’s nonlinear behavior serves as the basis for new and rational criteria for the design and the safe operation of rotating machines.

Unlocking a nose landing gear in different flight conditions: folds, cusps and a swallowtail

This paper investigates the unlocking of a non-conventional nose landing gear mechanism that uses a single lock to fix the landing gear in both its downlocked and uplocked states (as opposed to having two separate locks as in most present nose landing gears in operation today). More specifically, we present a bifurcation analysis of a parameterized mathematical model for this mechanical system that features elastic constraints and takes into account internal and external forces. This formulation makes it possible to employ numerical continuation techniques to determine the robustness of the proposed unlocking strategy with respect to changing aircraft attitude. In this way, we identify as a function of several parameters the steady-state solutions of the system, as well as their bifurcations: fold bifurcations where two steady states coalesce, cusp points on curves of fold bifurcations, and a swallowtail bifurcation that generates two cusp points. Our results are presented as surfaces of steady states, joined by curves of fold bifurcations, over the plane of retraction actuator force and unlock actuator force, where we consider four scenarios of the aircraft: level flight; steep climb; steep descent; intermediate descent. A crucial cusp point is found to exist irrespective of aircraft attitude: it corresponds to the mechanism being at overcentre, which is a position that creates a mechanical singularity with respect to the effect of forces applied by the actuators. Furthermore, two cusps on a key fold locus are unfolded in a (codimension-three) swallowtail bifurcation as the aircraft attitude is changed: physical factors that create these bifurcations are presented. A practical outcome of this research is the realization that the design of this and other types of landing gear mechanism should be undertaken by considering the effects of forces over considerable ranges, with a special focus on the overcentre position, to ensure a smooth retraction occurs. More generally, continuation methods are shown to be a valuable tool for determining the overall geometric structure of steady states of mechanisms subject to (external) forces.

Continuation and Bifurcation in Nonlinear PDEs – Algorithms, Applications, and Experiments

Numerical continuation and bifurcation methods can be used to explore the set of steady and time–periodic solutions of parameter dependent nonlinear ODEs or PDEs. For PDEs, a basic idea is to first convert the PDE into a system of algebraic equations or ODEs via a spatial discretization. However, the large class of possible PDE bifurcation problems makes developing a general and user–friendly software a challenge, and the often needed large number of degrees of freedom, and the typically large set of solutions, often require adapted methods. Here we review some of these methods, and illustrate the approach by application of the package to some advanced pattern formation problems, including the interaction of Hopf and Turing modes, patterns on disks, and an experimental setting of dead core pattern formation.