ON STEADY-STATE SOLUTIONS OF A WAVE EQUATION BY SOLVING A DELAY DIFFERENTIAL EQUATION WITH AN INCREMENTAL HARMONIC BALANCE METHOD

For wave propagation in periodic media with strong nonlinearity, steady-state solutions can be obtained by solving a corresponding nonlinear delay differential equation (DDE). Based on the periodicity, the steady-state response of a repeated particle or segment in the media contains the complete information of solutions for the wave equation. Considering the motion of the selected particle or segment as a variable, motions of its adjacent particles or segments can be described by the same variable function with different phases, which are delayed variables. Thus, the governing equation for wave propagation can be converted to a nonlinear DDE with multiple delays. A modified incremental harmonic balance (IHB) method is presented here to solve the nonlinear DDE by introducing a delay matrix operator, where a direct approach is used to efficiently and automatically construct the Jacobian matrix for the nonlinear residual in the IHB method. This method is presented by solving an example of a one-dimensional monatomic chain under a nonlinear Hertzian contact law. Results are well matched with those in previous work, while calculation time and derivation effort are significantly reduced. Also there is no additional derivation required to solve new wave systems with different governing equations.