У цій роботі використано чисельний метод для дослідження еволюції пошкоджень двовимірних анізотропних пластин. Напружено-деформована поведінка анізотропних матеріалів в умовах плоских напружень описана в поєднанні з критерієм розрахунку ефективного напруження Хілла. Для визначення еволюції полів деформацій і пошкоджень використано методи скінченних елементів та інкрементного навантаження. Чисельне моделювання було інтегровано з Abaqus для моделювання, Matlab для визначення параметрів і ParaView для візуалізації. Валідація в умовах одновісного розтягування пластин з вуглепластику з центральними отворами показала, що пошкодження починається від краю отворів і поширюється всередину, що ефективно ілюструється контурами ізопошкоджень. Ці результати забезпечують теоретичну основу для оптимізації проектування та оцінки безпеки анізотропних пластин для інженерних застосувань.