• Pilotage en force de matrices de contact rigides:
  • Création de la classe RigidToolMesher permettant d'associer un meshpoint à une matrice. Ce meshpoint est utile si on veut appliquer un élém de masse à la matrice rigide.
  • Vérification du code pour le dynamique.
  • Adaptation des tests:
    • forceDriven: il passe en dynamique et la matrice rebondit sur la matière (passage par une étape où la force sur la matrice est nulle).
    • ddrawing: on avait oublié Paul? il revient avec un test de deep drawing piloté en force. C'était un de mes premiers tests sur Metafor Fortran et je retrouve exactement les mêmes résultats. J'ai passé le test en dynamique pour le stabiliser. Il est ainsi possible d'augmenter la pénalité et éviter les pertes de contact problématiques en quasistatique.

• Ajout d'un test de l'import Gmsh:
  • J'ai ajouté un beau calcul (cfr screenshots) pour tester le module d'importation gmsh. Pour rappel celui-ci est très incomplet (importation de tértraèdres uniquement - cfr. projet cerveau).

Je veux ajouter une masse à ma matrice:

Je “maille” ma matrice (je lui ajoute un meshpoint):

mesher = RigidToolMesher(facset(201))
mesher.execute()

J'ajoute un FieldAplicator relatif à un élément de masse:

app = FieldApplicator(10)
app.push(201, SIDE_ID)
intset.copy(app)

materset.define (3, ConstantMassMaterial)
mater = materset[3]
mater.put(MASS_DENSITY, 1.0e-5)

prp3 = ElementProperties(Mass2DElement)
prp3.put (MATERIAL, 3)
intset(10).addProperty(prp3)

A mtGeo/RigidToolMesher.cpp
A mtGeo/RigidToolMesher.h

— Romain BOMAN 2007/01/10 17:45