Metafor

Calcul des vitesses et accélérations des degrés de liberté fixés

Suppression de oo_nda/tpi

• dans le cadre des aubes tournantes, cherchant à comprendre la différence de comportement vibratoire entre un modèle aube tournante et un modèle où les forces centrifuges sont ajoutées à un modèle dont l'aube est fixe, une piste était liée au schéma de Chung & Hulbert (à travers la pondération pas de temps 0 / pas de temps 1 des forces inertielles (via alpha_m) & extérieures (via alpha_f).
• une propriété booléenne des éléments (propelem) INERTIAINFEXT a été ajoutée telle que par défaut, les forces de gravité et centrifuge soient ajoutées dans le vecteur FExt (INERTIAINFEXT = true) et que l'on puisse les ajouter dans FInert (INERTIAINFEXT = false)
• tests :
  • apps/imp/aube3dFix1 : INERTIAINFEXT = true
  • apps/imp/aube3dFix2 : INERTIAINFEXT = false
  • apps/imp/aube3dFix2Num : INERTIAINFEXT = false & matrice de raideur numérique
  • apps/imp/aube2DFix1CarterRigidArcRot: INERTIAINFEXT = true
  • apps/imp/aube2DFix2CarterRigidArcRot: INERTIAINFEXT = false

• Lors de ces implémentations, j'ai repéré que les vitesses et accélérations des degrés de libertés fixés sont calculés non pas en utilisant les formules du schéma (voir thèse JPP pg VIII.25 : formules VIII.4.13 VIII.4.14 VIII.4.15), mais bien par un Euler arrière (pour des raisons de stabilité dans le cadre d'un schéma de Newmark :Thèse JPP pg VIII.28 et suivantes, formules VIII.5.8)
• Cette formulation permet d'éviter des oscillations sur des calculs fait via Newmark (cad jamais avec Metafor), mais entraîne les erreurs le bruit important à chaque division de pas de temps : ci dessous les contraintes et forces d'une “aube” en rotation avec variation (imposée/importante) du pas de temps en 0.014sec (+ variations du pas de temps lié à des accidents de convergence)

• Calcul des Vitesses & Accelerations DDL's fixés Euler Arrière :

• Calcul des Vitesses & Accelerations DDL's fixés Consistant schéma :

• Pour l'instant seul le schéma AlphaGeneralizedTimeIntegration est adapté
• le choix entre le calcul consistant avec le schéma des vitesses et accélérations se fait via la commande
  • mim.setConstrainedDofsVAComputationMethod(VACOMPUTATION_SCHEMECONSISTANT) (default)
  • mim.setConstrainedDofsVAComputationMethod(VACOMPUTATION_LINEAR)

• apps/imp/massSpring3d2NewmarkVAConsistant : démontre le problème de conditions limites en newmark
• apps/imp/massSpring3d2 : ancien test en Euler Arrière prend 6500 fois moins de pas de temps
• apps/imp/aube3DRotVAConsistant : Schéma consistant
• apps/imp/aube3DRotVALinear : Schéma Euler Arrière

• le module tpi n'étant qu'un gros copier-coller des éléments volumiques
• l'objectif était de mettre un espace de développement à disposition de Arnaud Collet
• vu qu'il n'a finalement pas développé sa thèse dans Metafor, le module n'a plus de raison d'être …

• stats.py :
  • Correction du serveur svn

 
A oo_meta\mtElements\boundaries\ShearContactValueExtractor.h/cpp
A oo_meta\mtElements\boundaries\PressureContactValueExtractor.h/cpp
A oo_meta\mtMaterials\volumesKevinVoigtViscoElastHypoGpState.h/cpp
A oo_meta\mtMaterials\volumesKevinVoigtViscoElastHypoMaterial.h/cpp
R 

 
A oo_meta/apps/monosMaterials2/KevinVoigtViscoElast.py
R  

— Luc Papeleux 2013/11/22