La difusión de estructuras compuestas exige métodos numéricos eficientes para simular el comportamiento dinámico de laminados elásticos con delaminaciones de interfaz con caras interactivas. Se propone un método avanzado de ecuación integral de contorno que emplea la transformada de Hankel de las matrices de Green para modelar la dispersión de ondas y analizar las frecuencias propias de delaminaciones circulares de interfaz parcialmente cerradas entre medios disímiles. Se introduce un caso más general de delaminación circular parcialmente cerrada utilizando condiciones de frontera de resorte con distribución de rigidez de resorte no uniforme. La desconocida apertura de la grieta se expande como una serie de Fourier con respecto a la coordenada angular y en términos de polinomios asociados de Legendre del primer tipo a través de la coordenada radial. El problema se descompone en un sistema de ecuaciones integrales de contorno y se resuelve utilizando el método de Bubnov-Galerkin. El método de ecuación integral de contorno se compara con el método sin malla y los trabajos publicados para un espacio homogéneo con una grieta abierta circular. Se demuestran los resultados del análisis numérico que muestran la eficiencia y la convergencia del método. El método propuesto podría ser útil para la identificación de daños empleando la información sobre las frecuencias propias estimadas experimentalmente. Además, se puede extender para composites multicapa con contacto imperfecto entre subcapas y múltiples delaminaciones circulares.