Dear getfem users,
I've been trying to implement an elastostatic problem with contacts using
the python wrapper for GetFEM v5.3. At first, I tried the contact bricks
with non-matching meshes without any issues. However, since my problem
considers large deformation, this approximation to contacts is not
satisfying.
I then tried the large sliding large deformation contact brick as following:
md.add_filtered_fem_variable('lambda_n', mfu, CONTACT)
contact_bid = md.add_integral_large_sliding_contact_brick_raytracing('r',
release_dist)
md.add_slave_contact_boundary_to_large_sliding_contact_brick(contact_bid,
mim, CONTACT, 'u', 'lambda_n')
md.add_master_contact_boundary_to_large_sliding_contact_brick(contact_bid,
mim_c, CONTACT_C, 'u_c')
where u, u_c, CONTACT and CONTACT_C represent the displacements and contact
boundaries for the slave and master meshes respectively.
Even though I get no errors at solving, the result appears to be incorrect.
After deformation, penetration between slave and master meshes is
significant, and changing the release distance for the raytracing
transformation does not seem to help (here I used three times the mean
element size). Also, I get a bunch of these in the console:
Level 1 Warning in getfem_contact_and_friction_common.cc, line 47: Inverted
element !-0.20579
I did read the note in the documentation saying that the brick is not
working. However, it does not specify for which version so I wanted to be
sure.
Is this brick working properly and if not, is there any version of GetFEM
in the repository that has one working version ?
Thank you in advance,
Pablo
PD. Keep up the good work, this is an amazing library !
