getfem-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

## Re: dof with singular mesh

 From: Edouard Oudet Subject: Re: dof with singular mesh Date: Tue, 11 May 2021 17:57:09 +0200 User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:78.0) Gecko/20100101 Thunderbird/78.7.1

```oups, sorry

https://www-ljk.imag.fr/membres/Edouard.Oudet/download/dof_singular.png

Best

On 5/11/21 5:46 PM, Yves Renard wrote:
```
```
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.

```
```
```
```
--
Edouard Oudet : http://www-ljk.imag.fr/membres/Edouard.Oudet/
IMAG - Bureau 164
700 avenue Centrale
38400 Saint Martin d'Hères
+33 (0)4 57 42 17 71  (office LJK)

```

reply via email to

 [Prev in Thread] Current Thread [Next in Thread]