Hello,
I have a unit square as my domain and I am solving the equations of linear elasticity on it. I want to impose no y-displacement on the lower edge and I want to fix the node located at the middle of the lower edge, that is at (0.5,0). I checked and I have a node at that point.
I am doing
def lower_boundary(x,on_boundary):
tol = 1E-15
return (abs(x[1]) < tol)
def lower_boundary_fixed_point(x,on_boundary):
tol=1E-15
return (abs(x[1]) < tol) and (abs(x[0]-0.5)<tol)
bc_lower=DirichletBC(V.sub(1),Constant(0),lower_boundary)
bc_lower_fixed_point=DirichletBC(V.sub(0),Constant(0),lower_boundary_fixed_point)
From bc_lower_fixed_point I am getting
*** Warning: Found no facets matching domain for boundary condition.
Could someone tell me what I am doing wrong?