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

Решить задачу теории упругости с односторонними граничными условиями методом конечных элементов.    
   Дана область в виде единичного квадрата. Сверху и снизу  и с лева область закреплена, но к центру, с права, на тело действует сила, в результате чего образуется трещина. Необходимо численно решить задачу с использованием схемы двойственности, основанной на функционале Лагранжа. Получить графические результаты. По возможности использовать я.п. Python 3+, но не обязательно. Постановка краевой и вариационной задачи имеется. Так же имеются некоторые другие выкладки формул. Литературой и статьями помогу, которые имеются. Так же буду стараться помогать со всем остальным.    
   Что требуется:    
1. Предоставить вывод формул и нахождение данных. (В любом виде, скриншеты/фотографии  для понимания)    
2. Программный код    
3. Решение в графическом виде    
4. Стабильная коммуникация в ходе всего процесса.