Розвинено математичну модель для аналізу контактної взаємодії двох електропружних трансверсально-ізотропних півпросторів під час стискання за наявності жорсткого дископодібного включення сталої товщини між ними. Модель базується на використанні зв’язаної системи рівнянь електропружності для кожного з п’єзоелектричних півпросторів. Аналітичний розв’язок задачі отримано за допомогою загального подання розв’язків рівнянь електропружності на основі гармонічних функцій, зведення задачі до розгляду інтегрального рівняння і розкладу шуканої функції за малим параметром.