Harmonic balance-boundary element and continuation methods for steady-state wave scattering by interior and surface-breaking cracks with contact acoustic nonlinearity
A numerical method for nonlinear steady-state wave scattering due to a crack with contact acoustic nonlinearity (CAN) is developed in this article. The in-plane scattering problem is addressed here. The two systems studied are composed of an unbounded elastic solid that includes an interior crack and an elastic half plane that includes a surface-breaking crack. A time-harmonic plane wave is incident, and clapping and rubbing motions on the crack faces are induced as nonlinear phenomena. A harmonic balance-boundary element method (HB-BEM) is utilized to solve the nonlinear steady-state problem. The system of equations obtained by the HB-BEM consists of an elastostatic boundary integral equation (BIE) and frequency-domain elastodynamic BIEs of several different frequencies, which are coupled by nonlinear boundary conditions. To solve the system numerically, the nonlinear boundary conditions are regularized. In the regularized boundary conditions, traction is expressed by a smooth single-valued function with respect to crack-opening displacement. The solution of the steady-state system is tracked with a parameter continuation method to investigate bifurcations. The parameters that vary in the continuation method are the incident frequency and static field variables. To determine the stability of the steady-state solutions, a stability analysis method based on Hill’s infinite determinant is proposed.
The stable solutions obtained by the proposed method agree with the conventional time-domain solutions after enough time has elapsed. A
-order subharmonic resonance is observed with bifurcations for particular parameters. The interaction between the local resonant properties due to the crack and the CAN causes subharmonic resonance.