В статье анализируется алгоритм решения задач об упругом и упругопластическом контактировании твердых деформируемых тел с учетом наличия в зоне контакта дополнительных сред с различными свойствами. Это могут быть нелинейные свойства поверхностных шероховатых слоев и пленок. Для учета влияния таких слоев в конечноэлементной схеме используются контактные конечные элементы, имеющие нелинейные свойства. Описан алгоритм, разработанный для решения задач с различными типами нелинейностей. На примере экспериментальных исследований, подробно описанных в работах [1, 2], выполнена идентификация нелинейных свойств контактных конечных элементов. На ряде примеров демонстрируется работоспособность алгоритма.
