Using exported XDMF mesh file from FreeCAD on FEniCS

[dokken]

There are many post on this throughout the forum (simply search for Measure or XDMFFile) and you Get results such as:

Gmsh mesh volume (and surface) is worth 0

I am facing another mesh problem, which I used to solve by using an old mesh format with gmsh (msh2), but this time this won’t fix the problem. I have the following mesh.geo: SetFactory("OpenCASCADE"); Box(1) = {0, 0, 0, 1e-2, 5e-3, 1e-4}; //+ Box(2) = {2.5e-3, 2.5e-3, 1e-4, 5e-3, 1e-4, 1e-4}; //+ BooleanFragments{ Volume{1}; Delete; }{ Volume{2}; Delete; } //+ Physical Volume("the_box") = {1}; //+ Physical Volume("thin box") = {2}; //+ Physical Surface("end_left") = {7}; //+ Physical Surface(…

and

Different results for the same mesh created differently (Hyperelasticity)

I claim that the meshes are not the same, see attached picture: [cube_comparison] } I would also recommend you not to use xml to load meshes, as it has been deprecated for a while. The suggested way of converting a mesh from gmsh is by using meshio. Also, to make sure that your visualization is P2 data is saved, use xdmf.write_checkpoint. Here is a script doing all of the suggested actions: from dolfin import * import meshio import numpy as np class Left(SubDomain): def inside(self, x,…