getfem-users
[Top][All Lists]
Advanced

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

Re: imposing sliding tangential conditions on part of the surface of a 3


From: Konstantinos Poulios
Subject: Re: imposing sliding tangential conditions on part of the surface of a 3D mesh
Date: Mon, 20 Dec 2021 21:39:52 +0100

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 <AJLesage@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.

reply via email to

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