[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] Theory behind mdbrick_Coulomb_friction.set_augmented_
Re: [Getfem-users] Theory behind mdbrick_Coulomb_friction.set_augmented_matrix
Wed, 18 Nov 2009 14:47:49 +0100
The mdbrick_Coulomb_friction brick contains a certain number of version due to
the fact that it has been adapted to different situations (unstationnary
problems, Barbosa-Hughes augmentation ...).
In fact, I intend to rewrite this brick soon for the new brick system.
There is two type of augmentation. There is the augmentation due to the
augmented lagrangian technique used to solve numerically the problem with a
Newton-Raphson algorithm (this is indeed presented in
The second augmentation is the Barbosa-Hughes one presented indeed in
This is not a mandatory augmentation. The advantage is that it adds a
contribution on the multipliers which helps to solve the system and allow to
have an inf-sup condition between the space of multipliers and the space of
displacements (for instance, it allows affine displacements and piecewise
constant multipliers). This augmentation is not fully taken into acocunt by
the brick (in fact there is just the possibility to take into account the
additional term wich modifies a bit the contac/friction condition). This
means that the additional term should be added by hand (as it is actually
done in contrib/static_friction.cc )
My question : do you really need this augmentation in your problem ?
If yes, I can give you some additional explanations.
Besides, you propose an algorithm to produce matrices BN and BT. This would be
a very interesting contribution of course !
My idea for the new contact/friction brick should be to write a basic brick
being essentially the same than the old one (BN, BT, gap, r are not computed
by the brick) and some derived bricks (for instance thanks to a virtual
function) which would be able to compute the matrices BN and BT in different
standard situations (elastic/rigid, elastic/elastic conforming mesh,
elastic/elastic non-conforming mesh) given the contact boundaries and the
space of multiplier (and an optional additional stailization / dynamic part).
In the new brick system, there is already a mecanism to basically ensure that
the discrete inf-sup condition is satisfied (it suppresses the multiplier
exceeding dofs). But of course this do no prevent from locking phenomena. So
I don't know if we have to use it for this contact/brick (the problem is when
there is a contact boundary next to a Dirichlet one).
If you are interested in, we could iterate to build this new brick. This
should be interesting and very usefull !
I also intend to write a more general contact/friction brick for large
deplacement/deformation problems. For the moment, this is very prospective
because I would like to build a robust algorithme. I found the existing
algorithm not really satisfactory. So this is for a longer term.
On samedi 14 novembre 2009, Konstantinos Poulios wrote:
> I am writing an application for gear calculations based on getfem and
> using the "mdbrick_Coulomb_friction" brick. The first results seem
> but I have difficulties to add augmentation to the Lagrange
> multipliers for the contact constraints. I 've already read the
> relevant material from http://math.univ-lyon1.fr/~renard/ -->
> Publications and
> http://www.lms.polytechnique.fr/users/khenous/publications.htm but I
> am not yet in a position to derive the part of the tangent matrix
> corresponding to the "set_augmented_matrix" method myself. Up to now
> all my attempts to derive this matrix yield to a zero matrix. I have
> already studied the contrib/static_friction.cc example but it didn't
> help me to resolve my question.
> I would like to ask for some hints or bibliography on the derivation
> of the tangent matrix of the "mdbrick_Coulomb_friction" brick in the
> case of augmentation. Please take into account that I haven't studied
> Mathematics but Mechanical Engineering and most of the books I have
> read on contact mechanics are from the point of view of an Engineer.
> Actually I understand more or less the theory presented in e.g.
> http://math.univ-lyon1.fr/~renard/papers/2007_HR_contact.pdf and
> http://www.lms.polytechnique.fr/users/khenous/files/coulomb_num.pdf .
> I also understand the implementation in getfem_Coulomb_friction.h
> pretty well. What I miss are some steps between the abstract theory
> and its practical implementation.
> Additionally I would like to offer my help to expand the contact
> mechanics infrastructure in getfem. I have already programmed a quite
> general contact searching algorithm which derives the "BN" and "gap"
> arrays for the "mdbrick_Coulomb_friction" brick. I could also help
> with the documentation of the contact brick.
> Best Regards
> Getfem-users mailing list
Yves Renard (address@hidden) tel : (33) 04.72.43.87.08
Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE