Procedures | Cast3M

   1 : * @STATIO   PROCEDUR  AM        98/12/30    21:15:23     3421           
   2 : DEBP @STATIO ttt*TABLE ;
   3 : 
   4 : *-----------------------------------------------------------------------
   5 : * @STATIO : teste si le probleme est elastique et gere les iterations
   6 : *          - si la solution est elastique : arret du programme
   7 : *          - sinon : lancement de la boucle des iterations (ALGSTA)
   8 : *            jusqu'a atteindre le nombre maximal, ou bien jusqu'a la 
   9 : *            convergence (. 'CONVERGENCE' . i = vrai)
  10 : *-----------------------------------------------------------------------
  11 : 
  12 : *** initialisation des tableaux de calcul
  13 : 
  14 :   ttt = @INITIA ttt ;
  15 : 
  16 : *** calcul du critere de plasticite pour le probleme a l'iteration 0
  17 : 
  18 :   ttt . 'CRITERE_PLASTICITE' . 0 = @CRIPL (ttt . 'MODELE')
  19 :         (ttt . 'CONTRAINTES' . 0) (ttt . 'VARIABLES_INTERNES' . 0)
  20 :         (ttt . 'CARACTERISTIQUES') 0. ;
  21 : 
  22 : *** test du critere de plasticite
  23 : 
  24 :   si (ttt . 'CRITERE_PLASTICITE' . 0 . PL) ;
  25 : 
  26 : *** probleme plastique : calcul stationnaire
  27 : 
  28 :   mess ' ' ;
  29 :   mess '***                CALCUL STATIONNAIRE                ***' ;
  30 :   mess ' ' ;
  31 : 
  32 :   repeter INCC (ttt . 'MAXITERATION') ;
  33 :     i = &INCC ;
  34 :     ttt = @ALGSTA i ttt ;
  35 :     si (ttt . 'CONVERGENCE' . i) ;
  36 :       quit INCC ;
  37 :     finsi ;
  38 :     si (i ega (ttt . 'MAXITERATION')) ;
  39 :       mess ' ' ; mess ' ' ;
  40 :       mess '*********************************************************';
  41 :       mess '           LE NB MAXI D ITERATIONS EST ATTEINT' ;
  42 :       mess '*********************************************************';
  43 :     finsi ;
  44 :   fin INCC ;
  45 : 
  46 : *** probleme elastique : arret du programme
  47 : 
  48 :   sinon ;
  49 :     mess '**** LA STRUCTURE NE PLASTIFIE PAS ****' ;
  50 :     mess '****        ARRET DU CALCUL        ****' ;
  51 :   finsi ;
  52 : 
  53 : FINP ;
  54 :