Télécharger @statio.notice

Retour à la liste

Afficher cette notice en

Numérotation des lignes :
$$$$ @STATIO  NOTICE  CHAT      11/09/12    21:18:12     7124           
                                             DATE     11/09/12



L'algorithme stationnaire a ete cree et developpe par DANG VAN K.,
MAITOURNAM H. et NGUYEN Q. S. du Laboratoire de Mecanique des Solides
a l'Ecole Polytechnique.

Pour en connaitre le principe, voir l'article paru dans "Journal of
the Mechanics and Physics of Solids", ecrit par DANG VAN K. et
MAITOURNAM H., intitule : "Steady-state flow in classical elastoplas-
ticity : applications to repeated rolling and sliding contact" (1993,
41(11), pp. 1691-1710).

Cette procedure a ete programmee en langage Gibiane par DRAGON M., en
these au Laboratoire de Mecanique des Solides.



  Procedure @STATIO                               Voir aussi :
    -----------------  
    @STATIO  TAB1 ;


    TAB1.  BLOCAGES_MECANIQUES
           CARACTERISTIQUES
           CHARGEMENT
           CONTRAINTES
           CONTRAINTES_PLASTIQUES
           CONVERGENCE
           CRITERE_PLASTICITE
           DEFORMATIONS
           DEFORMATIONS_PLASTIQUES
           DEPLACEMENTS
           EP2D
           EP3D
           FORCES_PLASTIQUES
           MAXITERATION
           MODELE
           MODELE_TABLE
           PRECISION
           VARIABLES_INTERNES
           VA2D
           VA3D


           
    Objet :
    _______

Cette procedure permet de calculer les deformations plastiques d'une 
structure soumise a un chargement mobile, dans l'etat stationnaire atteint
apres un grand nombre de cycles de chargement. 

On se place dans le repere lie au chargement, et l'etat stationnaire est 
determine directement. Exemples de telles structures : rail soumis au
passage repete d'une roue (2D), disque en rotation soumis au contact d'un 
pion (3D) ...


   Commentaire :
   _____________

Les indices de l'objet TAB1 sont des mots (a ecrire en toutes lettres,
et en majuscules s'ils sont mis entre cotes).

TAB1 contient les parametres du calcul a definir en entree d'une part, 
et les grandeurs mecaniques determinees au cours de la procedure d'autre
part.

     -  Liste des parametres a definir en entree :

 BLOCAGES_MECANIQUES : (type RIGIDITE) blocages mecaniques

 CARACTERISTIQUES    : (type MCHAML, sous-type CARACTERISTIQUES) champ de 
                       caracteristiques materielles, materiau elasto-plas-
                       tique

 CHARGEMENT          : (type CHPOINT) chargement defini sous la forme de
                       de forces nodales

 MAXITERATION        : (type ENTIER) nombre maximal d'iterations

 MODELE              : (type MMODEL) objet modele s'appuyant sur le mail-
                       lage entier de la structure, elasto-plastique 

 MODELE_TABLE        : (type TABLE) table indicee par des entiers de 1 a n,
                       contenant les n colonnes du maillage 2D ou les n 
                       parts (type MMODEL) du maillage 3D axisymetrique
                       (voir Remarques)

 PRECISION           : (type FLOTTANT) precision utilisee dans le calcul du
                       critere de plasticite, et dans le test de stationna-
                       rite (par exemple 1.e-3)


     -  Liste des parametres utilises au cours du calcul :

 Ces objets de type TABLE, contiennent sous les indices i designant les ite-
 rations, pour chacune de ces iterations (jusqu'a la convergence ou jusqu'a
 TAB1 . MAXITERATION) :

 CONTRAINTES             : (type MCHAML, sous-type CONTRAINTES) contraintes
                           obtenues par le calcul elastique

 CONTRAINTES_PLASTIQUES  : (type MCHAML, sous-type CONTRAINTES) (L : epsp),
                           ou L est la matrice de Hook et epsp le champ des 
                           deformations plastiques

 CONVERGENCE             : (type LOGIQUE) mis a VRAI lorsque le calcul
                           converge a la derniere iteration

 CRITERE                 : (type TABLE) table resultant du calcul du critere
                           de plasticite contenant les 3 indices suivants :
                           PL  : logique VRAI si la solution est plastique
                           NPL : nombre de points de Gauss ou la solution 
                                 est plastique
                           CR  : (type MCHAML) vaut 1 en chaque point ou le
                                 critere est viole, et 0 sinon

 DEFORMATIONS            : (type MCHAML, sous-type DEFORMATIONS) deformations
                           obtenues par le calcul elastique

 DEFORMATIONS_PLASTIQUES : (type MCHAML, sous-type DEFORMATIONS) deformations
                           plastiques calculees par l'algorithme stationnaire

 DEPLACEMENTS            : (type CHPOINT, sous-type DEPLACEMENTS) deplace-
                           ments obtenues par le calcul elastique

 EP2D                    : (type MCHAML, sous-type DEFORMATIONS) deformations
                           plastiques initiales pour chaque iteration en 2D

 EP3D                    : (type MCHAML, sous-type DEFORMATIONS) deformations
                           plastiques initiales pour chaque iteration en 3D

 FORCES_PLASTIQUES       : (type CHPOINT, sous-type FORCE) forces calculees
                           a partir des CONTRAINTES_PLASTIQUES

 VARIABLES_INTERNES      : (type MCHAML, sous-type VARIABLES INTERNES) vari-
                           ables internes

 VA2D                    : (type MCHAML, sous-type VARIABLES INTERNES) vari-
                           ables internes initiales pour chaque iteration en
                           2D

 VA3D                    : (type MCHAML, sous-type VARIABLES INTERNES) vari-
                           ables internes initiales pour chaque iteration en
                           3D


Remarques :
-----------

Cette procedure, du fait du probleme mecanique et de l'algorithme,
necessite un maillage particulier (voir MODELE_TABLE) :

2D : maillage forme de colonnes d'elements, contenant chacune le meme
     nombre d'elements en hauteur et un seul en largeur (comme un tab-
     leau, rectangle modelisant une portion de rail par exemple). Il peut
     etre realise par translation d'une premiere colonne, dans la direc-
     tion du mouvement de la structure.


       y          V mouvement du chargement
       |                  ------>
       |__x     ______________________________
               |___|       |___|___|      |___|
               |___|       |___|___|      |___|
               |___|       |___|___|      |___|
               |___|       |___|___|      |___|
               |___|       |___|___|      |___|
               |___|       |___|___|      |___|
               |___|       |___|___|      |___|
               |___|       |___|___|      |___|
               |___|_______|___|___|______|___|
                n          (i+1) i          1


3D : maillage axisymetrique, forme de parts (de gateau) toutes identiques,
     realise par exemple par rotation d'une premiere part autour de l'axe 
     Oz. (Camembert modelisant un disque par exemple.)

Un exemple peut etre consulte : sta2d.dgibi.

 
 

© Cast3M 2003 - Tous droits réservés.
Mentions légales