Se propone un método de elementos finitos equilibrados basado en fuerzas de interfaz, combinado con la teoría del análisis dual, para lograr una solución de tensiones de alta precisión en problemas de contacto con deslizamiento limitado, así como una estimación rigurosa del error a posteriori. Primero, en el marco del deslizamiento limitado, se deriva el principio de mínima energía residual para el problema de contacto y, mediante la introducción de restricciones de contacto (la fuerza de contacto como presión), se transforma en un problema de programación cuadrática; posteriormente, combinando la técnica de macroelementos, se construye un elemento equilibrado basado en fuerzas de interfaz e incorpora en el principio de mínima energía residual del problema de contacto, obteniendo finalmente un campo de tensiones equilibrado de alta precisión que satisface estrictamente las ecuaciones de equilibrio. Combinando dicho campo de tensiones equilibrado con el campo de desplazamientos coordinado de los elementos finitos, se deriva y calcula el error de las relaciones constitutivas en el problema de contacto, demostrando que este error es un límite superior estricto del error de discretización. Finalmente, el método propuesto se valida numéricamente mediante un ejemplo típico de contacto multicuásico. Los resultados muestran que el método garantiza la precisión de la solución, la convergencia y la efectividad de la estimación del error de discretización.

problema de contacto con deslizamiento limitado;elementos finitos equilibrados;principio de energía residual;algoritmo de búsqueda local;análisis dual