(Deleted) Cylinder Mesh and B.C. for an Elastic problem
3 months ago by
I'm currently trying to write with FeniCS a Linear Elastic problem over a Cylindrical domain (Very thin Cylinder, with base radius 400 times the height, here i used 80 to be faster), problems is that i'm very new to this and i don't know many things that are going on.
So far i wrote as follow:
from fenics import *
from dolfin import *
from mshr import *
# Define 3D geometry
Cylinder = Cylinder(dolfin.Point(0, 0, 0.5), dolfin.Point(0, 0, -0.5), 40, 40)
generator = CSGCGALMeshGenerator3D()
generator.parameters["mesh_resolution"] = 15.0;
generator.parameters["cell_size"] = 1.0;
# Mesh generator
mesh = generator.generate(CSGCGALDomain3D(Cylinder))
dolfin.plot(mesh, "3D mesh")
vtkfile = File("cylinder.pvd")
vtkfile << (mesh)
few questions: Is correct to write the Cylinder as a geometry and then mesh it? If yes, i used this Mesh Generator but i don't know practically what's the "mesh resolution" and "cell size".. So i get a bad mesh and overall long to compute, i tried changing the number a few times but still they are pretty bad, if someone knows how to correct please could him/she inform me.
i was also wondering wha't the most efficent way to mesh a 3d domain like this? Is there any way automate the mesh in order to have a good mesh e.g. all the elements with the same size and connectivities all well made?
Thank you everyone,
Community: FEniCS Project
Please login to add an answer/comment or follow this question.
The thread is closed. No new answer/comment may be added.