Parallelized computation of quasi-periodic solutions for finite element problems: A Fourier series expansion-based shooting method

TL;DR

This paper introduces a parallelized Fourier series expansion-based shooting method (FSE-Shooting) for efficiently computing quasi-periodic solutions in high-dimensional nonlinear finite element systems, with applications to stability analysis.

Contribution

The paper presents a novel FSE-Shooting method that enables parallelized computation of quasi-periodic solutions with multiple frequencies in complex PDE-governed systems.

Findings

Efficient computation of quasi-periodic solutions in high-dimensional systems.

Versatility demonstrated on systems with up to 1872 DOFs.

Capability to compute Lyapunov exponents for stability assessment.

Abstract

High-dimensional nonlinear mechanical systems admit quasi-periodic solutions that are essential for the understanding of the dynamical systems. These quasi-periodic solutions stay on some invariant tori governed by complex PDEs in hyper-time. Here, we propose a Fourier series expansion-based shooting method (FSE-Shooting) for the parallelized computation of quasi-periodic solution with $d$ base frequencies ($d \ge 2$). We represent the associated $d$-torus as a collection of trajectories initialized at a ($d-1$)-torus. We drive a set of ODEs that hold for any of these trajectories. We also derive a set of boundary conditions that couple the initial and terminal states of these trajectories and then formulate a set of nonlinear algebraic equations via the coupling conditions. We use Fourier series expansion to parameterize the ($d-1$)-torus and shooting method to iterate the Fourier coefficients associated with initial torus such that the coupling conditions are satisfied. In particular, the terminal points of these trajectories are parallelized computed via Newmark integration, where the time points and Fourier coefficients are transformed to each other by alternating Frequency-Time method. A straightforward phase condition is devised to track the quasi-periodic solutions with priori unknown base frequencies. Additionally, the by-product of the FSE-Shooting can be also directly used to compute the Lyapunov exponents to assess the stabilities of quasi-periodic solutions. The results of three finite element systems show the efficiency and versatility of FSE-Shooting in high-dimensional nonlinear dynamical systems, including a three-dimensional shell structure with $1872$ DOFs.