Hello,
I am comparing the results of 3D linear elasticity solutions with the corresponding analytical approximations. In the case of 3D cylinders (quadratic and circular) cross section that goes well. But in the case of “fat” cylinders (or plates) I do not succeed.
The homogeneous Dirchlet condition at the lateral boundary of the 3D plate domain is not found… " *** Warning: Found no facets matching domain for boundary condition."
The corresponding following code worked for the cylinder examples (where Dirchlet was at the longitudinal ends)
h = 0.5
r =5
mesh = Cylinder(Point(0, 0, 0), Point(0, 0, h), r, r)
mesh = generate_mesh(mesh, 100)
SyntaxEditor Code Snippet
tol = 0.001
def clamped_boundary(x, on_boundary):
return on_boundary and abs(r-np.sqrt(pow(x[0],2)+pow(x[1],2))) < tol
bc = DirichletBC(V, Constant((0, 0, 0)), clamped_boundary)
How can I the homogeneous Dirichlet conditions at the lateral boundary of the 3D plate domain?
Thank you!
Bezalel