Dear Anne-Cecile,
You are basically looking for something like:
mfu = gf.MeshFem(m, 3)
mfu.set_classical_fem(disp_fem_order)
mfmult = gf.MeshFem(m, 1)
mfmult.set_classical_fem(mult_fem_order)
md.add_fem_variable("u", mfu)
md.add_filtered_fem_variable("multn", mfmult, BRAINS_BOUND)
md.add_linear_term(mim, "(u.Normal)*multn", BRAINS_BOUND)
The last line is basically equivalent to:
md.add_linear_term(mim, "(u.Normal)*Test_multn + (Test_u.Normal)*multn",
BRAINS_BOUND)
The value for mult_fem_order can be the same as for disp_fem_order.
Best regards
Kostas
On Mon, Dec 20, 2021 at 9:07 PM Lesage,Anne Cecile J <
ajles...@mdanderson.org> wrote:
> Dear all
>
>
>
> In a python script, I defined a sliding boundary for a 3D tetra mesh
>
>
>
> meshb.set_region(BRAINS_BOUND, fbbs)
>
>
>
> How do I impose tangential sliding conditions for this surface U_n = 0
> with u being the displacement field?
>
>
>
> Thank you
>
> Regards
>
> Anne-Cecile
>
>
>
>
>
>
>
>
>
>
> The information contained in this e-mail message may be privileged,
> confidential, and/or protected from disclosure. This e-mail message may
> contain protected health information (PHI); dissemination of PHI should
> comply with applicable federal and state laws. If you are not the intended
> recipient, or an authorized representative of the intended recipient, any
> further review, disclosure, use, dissemination, distribution, or copying of
> this message or any attachment (or the information contained therein) is
> strictly prohibited. If you think that you have received this e-mail
> message in error, please notify the sender by return e-mail and delete all
> references to it and its contents from your systems.
>
