Dear Fabrice,
I don't really think such a brick exists (unless Yves corrects me).
But you can obtain the desired effect if you deform your mesh during
the assembly of you pressure Neumann term.
For that you need to define your problem not on getfem::mesh, but
getfem::deformable_mesh.
Subsequently, you can place a single line in the assembly of
normal_source_term_brick
getfem::temporary_mesh_deformator d (your_mesh, mesh_fem_of_displacement,
displacementVector);
So, this way your pressure term will be assembled on the current
configuration.
However, the above trick is not fully consistent, of course, as it won't
add the corresponding stiffness
contribution that comes from the changing geometry. So, your "follower"
pressure will be added in
explicit sense. If you want fully implicit procedure, I think you'll need
to create a new brick.
Best regards,
Andriy
On 23 October 2013 11:31, Fabrice Silva <[email protected]> wrote:
> Hi,
> Dealing with nonlinear elasticitya and lagrangian formulation, I need to
> apply an external pressure to a moving boundary. This is, in french
> words, usually called 'pression suiveuse' (see [1,2] for example), i.e.
> it takes into account the rotation of the normals of the faces.
> Is there such a brick in getfem ? If so, how to use it ?
>
> Best regards,
>
> [1] Code Aster documentation
> www.code-aster.org/V2/doc/default/fr/man_r/r3/r3.03.07.pdf
> [2] Fortin and Garon, page 305
> www.giref.ulaval.ca/~afortin/cours_elements_finis/documents/notes_elements_finis.pdf
>
> --
> Fabrice
>
>
> _______________________________________________
> Getfem-users mailing list
> [email protected]
> https://mail.gna.org/listinfo/getfem-users
>
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
