Решить задачу методом конечных элементов

Решить задачу теории упругости с односторонними граничными условиями методом конечных элементов.   

   Дана область в виде единичного квадрата. Сверху и снизу область закреплена, по бокам к центру на тело действует сила, в результате чего в теле образуется трещина. Необходимо численно решить задачу с использованием схемы двойственности, основанной на модифицированном функционале Лагранжа. Получить графические результаты. По возможности использовать я.п. Python 3+, но не обязательно. Постановка краевой и вариационной задачи имеется. Так же имеются некоторые другие выкладки формул. Литературой и статьями помогу, которые имеются. Так же буду стараться помогать со всем остальным.   

   Что требуется:   

1. Предоставить вывод формул и нахождение данных. (В любом виде, скриншеты/фотографии  для понимания)   

2. Программный код   

3. Решение в графическом виде   

4. Стабильная коммуникация в ходе всего процесса.