DOLFIN

symmetric boundary condition

Asked by Sungick Kim

In 3D problem,
if I want to fix x and y-direction displacements as 0, and let z-direction free to move,
how can I achieve this?
Thank you.

Question information

Language:
English Edit question
Status:
Solved
For:
DOLFIN Edit question
Assignee:
No assignee Edit question
Solved by:
Sungick Kim
Solved:
Last query:
Last reply:
Revision history for this message
Sungick Kim (sungick) said : 		#1

It was solved by Krishna Garikipati's help.

V = VectorFunctionSpace(mesh, "CG", 1)
xplane = compile_subdomains(" x[0] < 0 +5*DOLFIN_EPS && on_boundary")
displacement = Constant(0.0)
bc_m1 = DirichletBC(V.sub(0), displacement, xplane)