I am trying to implement the tutorial of COMET-FENICSX of linear elasticity problem in my computer and save results in vtk format. The results do get saved but when I open it using Paraview, I see no results but the image below
@Gaurav_Verma will add a link to the original tutorial, and highlight what modifications (if any) he made to the tutorial.
Please note that you are seeing some results (a stress field), except that it is all numerically zero. Without the code, the best guess is that you modified forcing terms or boundary conditions in a way that the displacement field ends up being identically zero.
Dear Francesco, this is the result I obtained by exactly executing the code provided in the tutorial. I tried multiple ways but with my limited experience and knowledge, I was not able to attain the results provided in the tutoriuals with the same code.
you must link the page containing the tutorial, otherwise we may not be able to come up with it on our own.
Dear Francesco, please see the link of the tutorial page below.
Can confirm that I obtain zero solution as well.
I suspect you’ll have to change the definition of L_form as follows:
ds = ufl.Measure("ds", domain=domain, subdomain_data=facets)
L_form = ufl.dot(T * n, u_) * ds((2, 3))
so that the integration measure ds is aware of the facet markers. cc @bleyerj who may want to fix that upstream.
2 Likes
Dear Francesco, thanks for your assistance. The code is working nicely now.
Thank you for spotting this. This has now been fixed