В статье представлен метод, позволяющий учитывать контактное взаимодействие между поверхностями на основе метода конечных элементов. Метод основан на итерационных циклах: определения статусов контактных элементов, обновления контактных усилий для выполнения условия нулевого проникновения всех активных контактных элементов.