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

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

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

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

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

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

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

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