El objetivo principal de este trabajo es proponer un nuevo método de elementos de contorno (BEM) para modelar la sensibilidad al estrés de problemas de daño termoelasto-plástico no local. La solución numérica de la ecuación de conducción de calor sujeta a una condición no local se describe utilizando un modelo de elementos de contorno. La cantidad total de energía térmica contenida dentro del sólido bajo consideración está especificada por la condición no local. El procedimiento para resolver la ecuación de calor revelará una función de control desconocida que rige la temperatura en una región específica del límite del sólido. El BEM de estrés inicial para estructuras con daño de ablandamiento por deformación se emplea en un programa de elementos de contorno con iteraciones en cada incremento de carga para desarrollar un modelo de plasticidad con deterioro del límite de fluencia. Para evitar las dificultades asociadas con el cálculo numérico de integrales singulares, la técnica de regularización es aplicable a los operadores integrales. Para validar la corrección física y eficiencia de la formulación sugerida, se resuelve un caso numérico.