Dear all
To describe a brain geometrical space contrains, I have three meshes part in my
python script
meshb=gf.Mesh('import','gid','brainr25tet4h5.GiD.msh')
meshtf=gf.Mesh('import','gid','tentoriumfalxr25h5.GiD.msh')
meshh=gf.Mesh('import','gid','headr25tet4h5.GiD.msh')
meshb and mesh are 3d and meshtf is a fine membrane described by a triangle
mesh surface
mesh and meshtf are attached at several places (nodes in common), how do I
ensure that the fem takes this boundary condition in account?
To set up that the brain is sliding without friction inside the head, I wrote
the following python lines
fbb = meshb.outer_faces() # Boundary of the brain
fbh = meshh.outer_faces() # Boundary of the head
HEAD_BOUND = 1;
BRAIN_BOUND = 2;
# selection of faces for contact condition
meshh.set_region(HEAD_BOUND, fbh)
meshb.set_region(BRAIN_BOUND, fbb)
md.add_penalized_contact_between_nonmatching_meshes_brick(mim3h, uh, ub,
datanamer, HEAD_BOUND, BRAIN_BOUND, 1, lambda1_n)
Is it a good option? How to choose the parameters datanamer and lambda1?
I would like the brain to slide along the whole meshtf. How to define a set of
faces that contains the entire triangles mesh?
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.
