A computational framework is proposed to path follow the periodic solutions of nonlinear spatially continuous systems and more general coupled multiphysics problems represented by systems of partial differential equations with time-dependent excitations. The set of PDEs is cast in first order differential form (in time) u˙ = f(u,s,t;c) where u(s,t) is the vector collecting all state variables including the velocities/time rates, s is a space coordinate (here, one-dimensional systems are considered without lack of generality for the space dependence) and t denotes time. The vector field f depends, in general, not only on the classical state variables (such as positions and velocities) but also on the space gradients of the leading unknowns. The space gradients are introduced as part of the state variables. This is justified by the mathematical and computational requirements on the continuity in space up to the proper differential order of the space gradients associated with the unknown position vector field. The path following procedure employs, for the computation of the periodic solutions, only the evaluation of the vector field f. This part of the path following procedure within the proposed combined scheme was formerly implemented by Dankowicz and coworkers in a MATLAB software package called COCO. The here proposed procedure seeks to discretize the space dependence of the variables using finite elements based on Lagrangian polynomials which leads to a discrete form of the vector field f. A concurrent bifurcation analysis is carried out by calculating the eigenvalues of the monodromy matrix. A hinged-hinged nonlinear beam subject to a primary-resonance harmonic transverse load or to a parametric-resonance horizontal end displacement is considered as a case study. Some primary-resonance frequency-response curves are calculated along with their stability to assess the convergence of the discretization scheme. The frequency-response curves are shown to be in close agreement with those calculated by direct integration of the PDEs through the FE software called COMSOL Multiphysics. Besides primary-resonance direct forcing conditions, also parametric forcing causing the principal parametric resonance of the lowest two bending modes is considered through construction of the associated transition curves. The proposed approach integrates algorithms from the finite element and bifurcation domains thus enabling an accurate and effective unfolding of the bifurcation and post-bifurcation scenarios of nonautonomous PDEs with the underlying structures.
Simulation of lid-driven cavity with top and bottom moving boundary conditions using implicit finite difference method and staggered grid
