El objetivo principal del presente trabajo es encontrar una solución analítica aproximada para la ecuación diferencial no lineal del oscilador de impacto vibratorio bajo la influencia de la actuación electromagnética cerca de la resonancia primaria. El desencadenante del régimen de impacto vibratorio se debe al contacto hertziano. El método de funciones auxiliares óptimas (OAFM) se utiliza para dar una solución aproximada analítica del problema. Se estudian completamente las influencias de la carga normal estática y la actuación electromagnética cerca de la resonancia primaria. Las principales novedades del procedimiento propuesto son la presencia de algunas nuevas funciones auxiliares adecuadas, la introducción de los parámetros de control de convergencia, la construcción original de la inicial y de la primera iteración, y la libertad para elegir el método para determinar los valores óptimos de los parámetros de control de convergencia. Todo esto condujo a una solución analítica explícita y precisa, que es otra novedad propuesta en el artículo. Esta técnica es muy precisa, simple, efectiva y fácil de aplicar utilizando solo la primera iteración. Un segundo objetivo fue realizar un análisis de estabilidad del modelo utilizando el método de múltiples escalas y los valores propios de la matriz jacobiana.