El modelado numérico es importante para explorar los procesos fundamentales que ocurren en la roca y para evaluar el rendimiento real de las estructuras construidas sobre y en el sistema de masa rocosa, y así apoyar el diseño de problemas de ingeniería en roca. Estimar la estabilidad de los sistemas de fundación de masa rocosa basándose únicamente en un enfoque teórico es una tarea complicada si existe superposición de sus modos de colapso potenciales. Este artículo aplica el análisis de límites por elementos finitos para evaluar la capacidad de carga de zapatas múltiples igualmente espaciadas que descansan sobre una masa rocosa que obedece al criterio de falla de Hoek-Brown no lineal modificado. Las soluciones numéricas se expresan en términos del factor de eficiencia que depende del espaciado entre las zapatas, así como de las propiedades de la masa rocosa. Además, se cuantifican los efectos de la sobrecarga superficial y la rugosidad de las zapatas. Se evalúa el espaciado máximo en el que el efecto interferente de las zapatas adyacentes desaparece y se propone una expresión algebraica para aproximar el espaciado máximo. Se examinan los mecanismos de falla para algunos casos de masa rocosa bajo múltiples zapatas.