How to compute invariant manifolds and their reduced dynamics in high-dimensional finite element models

Invariant manifolds are important constructs for the quantitative and qualitative understanding of nonlinear phenomena in dynamical systems. In nonlinear damped mechanical systems, for instance, spectral submanifolds have emerged as useful tools for the computation of forced response curves, backbone curves, detached resonance curves (isolas) via exact reduced-order models. For conservative nonlinear mechanical systems, Lyapunov subcenter manifolds and their reduced dynamics provide a way to identify nonlinear amplitude–frequency relationships in the form of conservative backbone curves. Despite these powerful predictions offered by invariant manifolds, their use has largely been limited to low-dimensional academic examples. This is because several challenges render their computation unfeasible for realistic engineering structures described by finite element models. In this work, we address these computational challenges and develop methods for computing invariant manifolds and their reduced dynamics in very high-dimensional nonlinear systems arising from spatial discretization of the governing partial differential equations. We illustrate our computational algorithms on finite element models of mechanical structures that range from a simple beam containing tens of degrees of freedom to an aircraft wing containing more than a hundred–thousand degrees of freedom.

Semi-Automated Procedure to Estimate Nonlinear Kinematic Hardening Model to Simulate the Nonlinear Dynamic Properties of Soil and Rock

The strain-dependent nonlinear properties of ground materials, such as shear modulus degradation (G/Gmax) and damping, are of significant importance in seismic-related analyses. However, the ABAQUS program lacks a comprehensive procedure to estimate parameters for a built-in model. In this study, a nonlinear kinematic hardening (NKH) model with three back-stress values was used, which allows better fitting to the backbone curves compared to the simplified nonlinear kinematic hardening (SNKH) model previously proposed. Instead of modeling in ABAQUS, a semi-automated procedure was implemented in MATLAB, which can predict shear stress–shear strain hysteretic loops, to find the fitting parameters to the target G/Gmax and/or damping curves. The procedure was applied for three soil and two rock samples, and the results indicate a good match between model and target backbone curves, which proves the application of the procedure and the NKH model in simulating the nonlinear properties of ground materials.

Backbone curves, Neimark-Sacker boundaries and appearance of quasi-periodicity in nonlinear oscillators: application to 1:2 internal resonance and frequency combs in MEMS

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.

An Accurate Numerical Model Simulating Hysteretic Behavior of Reinforced Concrete Columns Irrespective of Types of Loading Protocols

In older reinforced concrete (RC) buildings, columns are fragile elements that can induce collapse of entire buildings during earthquakes. An accurate assessment of the seismic vulnerability of RC buildings using nonlinear response history analyses requires an accurate numerical model. The peak-oriented hysteretic rule is often used in existing numerical models to simulate the hysteretic behavior of RC members, with predefined backbone curves and cyclic deterioration. A monotonic backbone curve is commonly constructed from a cyclic envelope. Because cyclic envelope varies according to loading protocols, particularly in a softening branch, it is difficult to obtain a unique backbone curve irrespective of loading protocols. In addition, cyclic deterioration parameters irrespective of loading protocols cannot be found because these parameters are estimated with respect to the backbone curves. Modeling parameters of existing numerical models can also vary with respect to loading protocol. The objective of this study is to propose a loading protocol-independent numerical model that does not require estimates of modeling parameters specifically tuned for a certain loading protocol. The accuracy of the proposed model is verified by comparing the simulated and measured cyclic curves of different sets of identical RC column specimens under various loading protocols.