Table of Contents
Commit 2015-01-21
Etat plan défo : épaisseur hors plan variable
Il est possible de faire varier l'épaisseur hors plan associée à une iteraction. Les commandes ci-dessous montrent comment utiliser cette fonctionnalité avec des éléments 2D.
prp1 = ElementProperties(Volume2DElement) prp1.put (CAUCHYMECHVOLINTMETH,VES_CMVIM_SRIPR) prp1.put (STIFFMETHOD, STIFF_ANALYTIC) prp1.put (THICKNESS, p['hzz'] ) prp1.depend(THICKNESS, fctEvol, Field1D(TM,RE))
L'aire nodale est calculée comme la somme des demi-longueurs des segments incidents aux noeuds de contact multipliée par l'épaisseur hors plan des éléments volumiques.
Validation
Pour valider la modification, j'ai fait varier l'épaisseur hors plan d'un carré de 1 mm de coté (cfr. cas hVar_2D_… dans la liste ci-dessous) et j'ai comparé les résultats obtenus (pression, contrainte sigmazz, déplacements d'un noeud) avec ceux obtenus lors d'une simulation de traction sur un cube de 1 mm de côté (cfr. cas hVar_ref3D_…).
Les cas tests ajoutés dans la batterie sont listés ci-dessous. Les éléments SRI - SRIPR - EAS sont testés.
[a] apps/gps/hVar_2D_EAS.py [a] apps/gps/hVar_2D_SRI.py [a] apps/gps/hVar_2D_SRIPR.py [a] apps/gps/hVar_ref3D_EAS.py [a] apps/gps/hVar_ref3D_SRI.py [a] apps/gps/hVar_ref3D_SRIPR.py