Dear Edouard,
Your second picture is missing. Can you send the link again ?
Best regards,
Yves
Le 11/05/2021 à 16:11, Edouard Oudet a écrit :
Dear getfem users,
Using the python interface of Getfem++, I obtained surprising results
related to the representation of degrees of freedom for a singular
mesh associated to a square with two circular holes. One of the hole
has a very small radius see picture
https://www-ljk.imag.fr/membres/Edouard.Oudet/download/singmesh.png
This mesh can be downloaded at
https://www-ljk.imag.fr/membres/Edouard.Oudet/download/savemesh.mesh
Using the following instructions
filemesh = '/tmp/savemesh.mesh'
tmesh = getfem.Mesh("load", filemesh)
tmeshfem = getfem.MeshFem(tmesh, 1)
tmeshfem.set_fem(getfem.Fem('FEM_PK(2, 4)'))
pdof = getfem.MeshFem.get(tmeshfem, 'basic dof nodes')
and plotting these pdof points I obtained the following picture around
the small circle
https://www-ljk.imag.fr/membres/Edouard.Oudet/download/savemesh.mesh
As it can be observed some dof are missing.. I do not get this problem
for a degree smaller than 4.
Any suggestions are welcome,
best,
Edouard.
--
Yves Renard ([email protected]) tel : (33) 04.72.43.87.08
INSA-Lyon
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE
http://math.univ-lyon1.fr/~renard
---------
