-
Notifications
You must be signed in to change notification settings - Fork 15
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
FEM for 3D Tet mesh #46
Comments
Hi! We would definitely be interested. A very easy place to start incorporating tet FEM to gpytoolbox is all the After that, the obvious next step would be Another eventual issue will be writing a reader and writer of tetrahedral meshes (our current Let us know if we can help! Edit: Oh, and I forgot to answer two of your questions: yes, we would like everything to be numpy/scipy-based (like our triangle mesh FEM implementations). And no, we definitely do not need it next week :-) |
Thanks for the fast reply! Some more questions to restrictions of the FEM
|
Indeed, there is currently no I don't think we have any powerful reason to restrict ourselves to linear, tetrahedral elements, but for the fact that they make life much easier. Thus if I were writing this, I would start by doing linear tets, then quadratic tets, then other elements. As long as (like you say) the user can easily know what they're calling and it has a reasonable default, I have no opposition to the code being richer and would be happy to merge it. I have no strong opinions on any of the other points you raise |
Ok thanks, then I will write the methods to make them easily extendable and annotate where it can be done. |
Sounds good! Looking forward to it :-) |
@sgsellan @odedstein I would suggest to use the gmsh Does this sound reasonable to you? |
I like the gmsh .msh format. I think it's a reasonble choice. |
Just an update for you, im currently working on the gmsh importer (works fine already). I hadn't much spare time recently. |
Are you interested in me implementing the FEM for 3D Tets?
I already have a working version for every element type using gmsh.
You probably want to stay with numpy / scipy only, correct?
If you don't need it next week, I can adapt it (for linear elastic tet elements) to only use numpy/scipy 😀.
Can you give me some hints where to find existing code/functionality for 3D Tets in gpytoolbox?
Do you use voigt notation for stresses/straits?
Or diagonal notation for stresses/straits?
(only relevant for anisotropic materials)
The text was updated successfully, but these errors were encountered: