Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools


commit:2008:11_18

Commit 2008-11-18

  • Forces d'inertie de rotation

but du commit

  • Pour une série de cas où la rotation d'un dispositif flexible (type aube) est opposé à un dispositif plus rigide (genre carter même si il est déformable), on a l'impression que changer de référentiel, l'observateur tournant avec l'aube et voyant tourner le carter permettrait de réduire le temps de calcul en ne calculant que la vibration du dispositif flexible libérée de sa rotation rigide.
  • Méthode pouvant aussi intéresser Romain dans le cadre du TubeSpinning
  • Pour ce faire, il faut calculer (efficacement si possible) les forces liées à l'accélération centrifuge pour tous les schémas

Modifs

  • Calcul des forces d'inertie à l'aide de la matrice de masses consistante
    • Lorsque un élément se déplace réellement, on calcule les forces d'inertie par Forces = MatriceMasseConsistante * Accélération. Pourquoi faire autrement ici ?
    • Pré calcul de l'axe, vitesse de rotation, point de passage de l'axe dans l'élémentShcut (eviter de faire des get & getPos inutiles)
    • Unification avec les accélération linéaires (gravité) dans le vecteur des forces externes (encore une fois pourquoi faire différement)
    • Suppression du GEN_ROT_FORC (vu que c'est dorénavant une force externe)
    • Matrice de raideur analytique : K = MatriceMassesConsistante * dAcceleration/dX
  • Nettoyage des schémas
  • Auparavant, les forces d'inertie en rotation étaient uniquement accessibles en quasiStatique ou dans la phase d'initialBalancing en rotation des schémas dynamiques (IQSI). Qui plus est, les forces d'inertie remplissaient vitesse & accélérations nodale durant tout le calcul (inutile et dangereux).
  • Dorénavant, la définition d'une rotation dans le propelem génère en permanence les forces d'inertie sauf dans le cas où IQSI = 1 en dynamique auquel cas on refait l'inital balancing (qui pourrait presque être défini comme un schéma QS classique (à condition de définir le facteur de rotation) ⇒ un pas de + vers l'intégration des TimeIntegration dans les stages) suivit d'une phase de “transition vers le schéma dynamique” (remplissage des vitesses & accélérations (fillRotationVA), transfert force externes ⇒ inertie, mise à 0 du facteur de rotation pour qu'on ne tienne plus en compte la rotation…)
  • Maintenant que j'y pense le restart ne fonctionera pas en IQSI1 (rotation factor n'étant pas sauvé dans le fac!!!)⇒ to be done
  • Divers
    • Ajout de Component Product par des Matr2/Matr3 dans Vect2/Vect3

Futur Works - ToDo list

  • restart en IQSI1 (sauver le rotationFactor dans le fac)
  • élargir la base de tests (newac / disques / …)
  • Ajout d'un test dans le LineSearch lorsque le résidu initial LS est > 0 ⇒ rester_ls = false (apparait lorsque l'on a des negativeJacobian ou des DSSPivotNul)

Tests

  • added apps.imp.aube2DFixCarterRigidArcRot
  • Correction ou ajout d'extractors (ca peu servir pour savoir si un test marche encore ou non)

Projet

  • Up To Date

Fichiers ajoutés/supprimés

Code

A 
R 

Test

A oo_meta\apps\imp\aube2DFixCarterRigidArcRot.py  

Luc PAPELEUX 2008/11/18

commit/2008/11_18.txt · Last modified: 2016/03/30 15:23 (external edit)