Чисельно досліджено напружено-деформований стан шаруватих пластин і оболонок із композиційних матеріалів. У рамках дискретного розгляду шарів систематизовано основні варіанти математичних моделей і побудовано відповідні функціональні енергії з урахуванням міжшарових пошкоджень. Виявлено особливості варіаційних постановок задач. На основі змішаної форми методу скінченних елементів розроблено розрахункові схеми і спеціальні алгоритми побудови чисельних розв'язків, запропоновано об'єктно-орієнтований підхід до програмної реалізації. Аналіз отриманих розв'язків низки задач згину дво- і тришарових пластин і циліндричних оболонок з ідеальним і неідеальним міжшаровим контактом показує перспективність даного підходу до оцінки контактної міцності шаруватих структур на стадії їх проектування.