Regularized forms of the traction and tangent derivative boundary integral equations of elasticity are derived for the case of closed regions. The hypersingular and strongly singular integrals of the displacement gradient representation are regularized independently, through identities of the fundamental solution and its various derivatives, before the boundary integral equations are formed. Besides the displacements and the tractions, only the tangential derivatives of the displacements evaluated at the singular point appear in the regularized equations making them well suited for numerical treatment. The regularization of the hypersingular integrals demands that the displacement components have Holder continuous first derivatives at the singular point. Consistent with this requirement, the regularization of the strongly singular integrals is effective if the tractions and the unit vectors normal and tangent to the surface are continuous at that location.;Higher order elements for two and three dimensional elastostatic problems are implemented through the coincident collocation of regularised forms of the displacement and the tangent derivative equations. The nodal values of the displacements, the fractions and their tangential derivatives are used as the degrees of freedom associated with the functional representation of the boundary variables. The tangential derivatives of the displacements and the tractions at the functional nodes are directly recovered from the boundary solution with comparable accuracy as the primative variables. Hence, the nodal values of the stress components are directly obtained through Hooke's law and need not be determined in a post processing manner. Several numerical examples demonstrate the advantages of the higher order elements versus the conventional ones. In two dimensions, four degrees of freedom per node Hermitian elements are used for functional interpolation only on those portions of the boundary where the gradients are high and quadratic Lagrangian elements are employed for the remaining parts of the modelled region. In three dimensions, nine degrees of freedom per node, incomplete quartic elements are employed for the approximation of the displacements and the tractions. Finally, the methodology presented here is general and can be extended to other problems amenable to a boundary integral formulation.

The direct multi-region boundary element method (BEM) is developed to simulate the interaction of ultrasonic waves with submerged finite elastic bodies. This technique is further specialized to model fluid-solid half-space problems with fluid-filled indentations and air bubbles on the solid surface. These problems usually occur in ultrasonic immersion testing and are of particular interest to the NDE community. The problems are formulated in the Fourier (frequency) domain and are described by the three dimensional acoustic and elastodynamic boundary integral equations (BIE) with pressure and displacement serving as primary variables. The techniques developed are general and may be used with any kind of incident wave, although plane waves are used in all numerical experiments. The equations governing the acoustic region are first converted mathematically to equations like those of an elastodynamic region and the resulting two elastodynamic regions are coupled and solved for the interface displacements. Tractions, pressures, and pressure gradients are subsequently computed using these displacements and the interface conditions. Numerical results are obtained for some typical problems and justified using reciprocity relations and by comparison with solutions available for the half-space elastodynamic problem.