Nonlinear forced vibration analysis of postbuckled beams

A numerical solution methodology is proposed herein to investigate the nonlinear forced vibrations of Euler–Bernoulli beams with different boundary conditions around the buckled configurations. By introducing a set of differential and integral matrix operators, the nonlinear integro-differential equation that governs the buckling of beams is discretized and then solved using the pseudo-arc-length method. The discretized governing equation of free vibration around the buckled configurations is also solved as an eigenvalue problem after imposing the boundary conditions and some complicated matrix manipulations. To study forced and nonlinear vibrations that take place around a buckled configuration, a Galerkin-based numerical method is applied to reduce the partial integro-differential equation into a time-varying ordinary differential equation of Duffing type. The Duffing equation is then discretized using time differential matrix operators, which are defined based on the derivatives of a periodic base function. Finally, for any given magnitude of axial load, the pseudo -arc-length method is used to obtain the nonlinear frequencies of buckled beams. The effects of axial load on the free vibration, nonlinear, and forced vibrations of beams in both prebuckling and postbuckling domains for the lowest three vibration modes are analyzed. This study shows that the nonlinear response of beams subjected to periodic excitation is complex in the postbuckling domain. For example, the type of boundary conditions significantly affects the nonlinear response of the postbuckled beams.