Télécharger la plaquette Cast3M


fr en

Télécharger @MODTRI.procedur

Retour à la liste

Numérotation des lignes : 

  1. * @MODTRI   PROCEDUR  BR232186  14/06/03    21:15:01     8061           
  2. DEBP @MODTRI TABPARA * TABLE DATA * TABLE;
  3. **********************************************
  4. * PROCEDURE @MODTRI.PROCEDUR
  5. **********************************************
  6. * Entree :      TABPARA (type TABLE)
  7. *               TABLE DE PARAMETRES
  8. *
  9. *               DATA (type TABLE)
  10. *               TABLE DE VARIABLES INTERNES
  11. *               INITIALES
  12. **********************************************
  13. * Sortie :      TAB3 (type TABLE)
  14. *               TABLE DE VARIABLES INTERNES
  15. *               A CONVERGENCE
  16. **********************************************
  17. * Commentaire : Loi de comportement en 
  18. *               variables generalisees
  19. **********************************************
  20. * Developpeur : COMPUTATIONAL ENGINEERING 
  21. *               AND STRUCTURES SAS
  22. *               contact@c-e-s.fr
  23. **********************************************
  24.  
  25. * Parametres materiaux
  26. KB1 = TABPARA . 1;
  27. KB  = TABPARA . 1;
  28. KB2 = TABPARA . 2;
  29. KB3 = TABPARA . 3;
  30. X0  = TABPARA . 4;
  31. X1  = TABPARA . 5;
  32.  
  33. * Variables internes
  34. DPLUS = DATA . 1;
  35. DMOIN = DATA . 2;
  36. DJ    = DATA . 3;
  37. MAXDP = DATA . 4;
  38. MAXDM = DATA . 5;
  39.  
  40. * Calcul des seuils et de l endommagement
  41. SI (DJ >EG 0.0);
  42.  
  43. SI (DJ > MAXDP);
  44. MAXDP = DJ;
  45. FINSI;
  46.  
  47. SI (MAXDP > X0);
  48. TOTO  = ((KB1 * X0) + (KB2 * (MAXDP - X0)));
  49. DPLUS = (1.0 - (TOTO / (MAXDP * KB1)));
  50. FINSI;
  51.  
  52. SI (MAXDP > X1);
  53. TOTO  = (((KB1 * X0) + (KB2 * (X1 - X0)))
  54.  + (KB3 * (MAXDP - X1)));
  55. DPLUS = (1.0 - (TOTO / (MAXDP * KB1)));
  56. FINSI;
  57.  
  58. SINON;
  59.  
  60. SI (DJ < MAXDM);
  61. MAXDM = DJ;
  62. FINSI;
  63.  
  64. SI (MAXDM < (-1.0 * X0));
  65. TOTO  = ((KB1 * X0) + (KB2 * (-1.0 * (MAXDM + X0))));
  66. DMOIN = 1.0 - (TOTO / ((-1.0 * MAXDM) * KB1));
  67. FINSI;
  68.  
  69. SI (MAXDM < (-1.0 * X1));
  70. TOTO = (KB1 * X0); 
  71. TOTO = TOTO + (KB2 * (X1 - X0));
  72. TOTO = TOTO + (KB3 * (-1.0 * (MAXDM + X1)));
  73. TOTO = TOTO / (-1.0 * MAXDM);
  74. TOTO = TOTO / KB1;
  75. DMOIN = 1.0 - TOTO; 
  76. FINSI;
  77. FINSI; 
  78.  
  79. SI (DJ >EG 0.0);
  80. FO1 = KB * (DJ * (1.0 + (-1.0 * DPLUS)));
  81. FF  = DPLUS;
  82. SINON;
  83. FO1 = KB * (DJ * (1.0 + (-1.0 * DMOIN)));
  84. FF  = DMOIN;
  85. FINSI;
  86.  
  87. * Rendu des variables en sortie de la procedure
  88. TAB3 = TABLE;
  89. TAB3 . 1 = FO1;
  90. TAB3 . 2 = DPLUS;
  91. TAB3 . 3 = DMOIN;
  92. TAB3 . 4 = MAXDP;
  93. TAB3 . 5 = MAXDM;
  94. TAB3 . 6 = FF;
  95.  
  96. RESPRO TAB3;
  97.  
  98. FINP;
  99.  
  100.  

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