Re: [Libmesh-users] Adding a force to a boundary through an additional node

2021-06-27 Thread Vikram Garg
Direct modification of the matrix (in general) to apply such a constraint will need consideration of the basis type and order, and how these map to the dofs. I am not sure this will be that easy to do, but I might be wrong. Vikram Garg vikramvgarg.github.io/ On Sun, Jun 27, 2021 at 9:00 PM Rena

Re: [Libmesh-users] Adding a force to a boundary through an additional node

2021-06-27 Thread edgar
On 2021-06-28 00:18, Renato Poli wrote: Hi, thanks for the answer: @edgar: it is a FORCE (Neumann) while DISPLACEMENTS (primary variables) must be tied together ... Please let me know if it is clear. Totally. Thanks. I wish I could do more :) . Thanks, Renato _

Re: [Libmesh-users] Adding a force to a boundary through an additional node

2021-06-27 Thread Renato Poli
Hi, thanks for the answer: @edgar: it is a FORCE (Neumann) while DISPLACEMENTS (primary variables) must be tied together. (I am trying to solve Mandel's problem in poroelasticity). Please let me know if it is clear. @Vikram: I did not get exactly the use of the SCALAR variable. It is a single val

Re: [Libmesh-users] Adding a force to a boundary through an additional node

2021-06-27 Thread Vikram Garg
You might need a Lagrange multiplier formulation that sets a state variable in an entire subdomain or boundary to a single unknown scalar. Example 3 in systems of equations might help, it shows how SCALAR variables can be used. You could incorporate the scalar variable into the weak form via a pena

Re: [Libmesh-users] Adding a force to a boundary through an additional node

2021-06-27 Thread edgar
On 2021-06-27 23:01, Renato Poli wrote: I'd like to add a force to a whole boundary. The constraint is that the whole boundary must have the same displacement. ---8<--- snip For the sake of disambiguation: do you want to impose displacement or force? The displacement of a boundary is not only

[Libmesh-users] Adding a force to a boundary through an additional node

2021-06-27 Thread Renato Poli
Hi, I'd like to add a force to a whole boundary. The constraint is that the whole boundary must have the same displacement. It seems desirable to add the force to an extra node, and then constraint the whole boundary to this note through add_constrain_row. Is that so? Do you see any better appro