Hello Community,
I want to solve following contact problem known as the signorini problem:
-div sigma (u) = f_i on Omega
u = 0 on boundary_1
sigma_ti(u) = 0 on boundary_2
under the complementary conditions
G sigma_n(u) = 0
G >= 0
sigma_n(U) >= 0
all on boundary_2.
G is a function R^2 -> R describing the distance between a rigid surface and Omega and sigma stands for the cauchy stress tensor.
I want to solve this by transforming the complementary conditions via the Fischer Burmeister Function.
My Problem is how to implement this function into the solving process. Or is there a more easier way to implement complementary boundary conditions into FEniCS ?
every hint for a solution is appreciated.
thanks
Roman
