Télécharger @statio.procedur

Retour à la liste

Numérotation des lignes :

  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.  
  55.  

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