Numerical investigation of periodic solutions for a suspension bridge model

Numerical methods are discussed for the computation of periodic solution paths for a suspension bridge model(UNFORMATTED TABLE OR EQUATION FOLLOWS){dollar}{dollar}{bsol}cases{lcub}&{lcub}{bsol}it u{rcub}{bsol}sb{lcub}{bsol}it tt{rcub} + {bsol}delta{lcub}{bsol}it u{rcub}{bsol}sb{lcub}{bsol}it t{rcub} + EIu{bsol}sb{lcub}xxxx{rcub} + Ku{bsol}sp+ = W + {bsol}lambda f(x,t,{bsol}mu){bsol}cr{bsol}cr& {lcub}{bsol}it u{rcub}(0, {lcub}{bsol}it t{rcub}) = {lcub}{bsol}it u{rcub}(1, {lcub}{bsol}it t{rcub}) = {lcub}{bsol}it u{rcub}{bsol}sb{lcub}{bsol}it xx{rcub}(0, {lcub}{bsol}it t{rcub}) = {lcub}{bsol}it u{rcub}{bsol}sb{lcub}{bsol}it xx{rcub}(1,{lcub}{bsol}it t{rcub}) = 0,{bsol}cr{rcub}{dollar}{dollar}(TABLE/EQUATION ENDS)where {dollar}{bsol}delta, E, I, K,{dollar} and {dollar}W{dollar} are constants, {dollar}{bsol}lambda{dollar} and {dollar}{bsol}mu{dollar} are parameters, and {dollar}f{dollar} is a periodic function in {dollar}t{dollar}.; The shooting method with Newton's and quasi-Newton's solver, and pseudoarclength continuation method around simple turning points are treated. The initial value problem solver that we adopted in the shooting method is implicit in the linear part, and explicit in the nonlinear part and forcing terms. We call it two-step Crank-Nicolson scheme.; Multiple periodic solutions for the suspension bridge equations are found, with their stability investigated. The unconditional stabilities of the Crank-Nicolson implicit scheme for solving suspension bridge equations in both infinite domain and finite domain are proved.