El objetivo principal de este trabajo es desarrollar un nuevo método de elementos de contorno (BEM) para analizar la propagación de ondas ultrasónicas en medios porosos viscoelásticos anisotrópicos de tres temperaturas. Debido a los problemas de la fuerte no linealidad de la propagación de ondas ultrasónicas en medios porosos de tres temperaturas, las soluciones analíticas o numéricas a los problemas bajo consideración son siempre desafiantes, lo que hace necesario el desarrollo de nuevas técnicas computacionales. Como resultado, utilizamos un nuevo modelo BEM para resolver tales problemas. Se obtiene un procedimiento de paso de tiempo basado en el método lineal de pasos múltiples después de resolver la ecuación integral de contorno discretizada con la regla de cuadratura. El cálculo de una doble integral es necesario para obtener soluciones fundamentales, pero esto aumenta el tiempo total de cálculo de BEM. Nuestra técnica BEM propuesta se utiliza para resolver el problema actual y mejorar la eficiencia de la formulación. Los resultados numéricos se grafican para demostrar los efectos de la viscosidad y la anisotropía en las ondas de estrés ultrasónicas no lineales en medios porosos de tres temperaturas. La validez, precisión y eficiencia de la metodología propuesta se demuestran al comparar los resultados obtenidos con una solución correspondiente obtenida del método de diferencias finitas (FDM).