I was wondering if it is possible to solve the following problem in FreeFEM
in the space
The variational formulation is
I tried using Morley elements to try and solve this in 2D, but I get an error when imposing the Dirichlet boundary conditions on u only. It says that all elements of the vectorial finite element should be in the Dirichlet boundary condition.
Is there a way of solving the above problem in FreeFEM ? Thank you for your advice.