[Libmesh-users] Discontinuous base with springs

[Renato P. <re...@gm...>]

Hi

I need to represent fractures in my elastic system. For that I would
like to use a discontinuous base (L2_LAGRANGE variable?) and tie
overlapping dofs with a penalty function, representing springs
(probably weakly).

Is there any example to follow?
How can I efficiently find the dofs to tie?
I should then use the "add_constraint_row", right?

I found another application in the list that uses the
"find_nodal_neighbors" function.
Is that the right thing to do?

Thanks
Renato