Procedures | Cast3M

   1 : * @CRIPL    PROCEDUR  CHAT      05/02/04    21:15:02     5023           
   2 : DEBP @CRIPL MOD*MMODEL SIG*MCHAML VAR*MCHAML MAT*MCHAML ER*FLOTTANT ;
   3 : 
   4 : *-----------------------------------------------------------------------
   5 : * @CRIPL : calcul du critere de plasticite : f/sy = (J2 - sy)/sy < ER
   6 : *         - . PL  : vrai si f/sy > ER (solution plastique)
   7 : *                   faux si f/sy < ER (solution elastique)
   8 : *         - . NPL : nb de points de Gauss ou f/sy > ER
   9 : *         - . CR  : champ par elements, 1 si f/sy > ER
  10 : *                                       0 si f/sy < ER
  11 : *-----------------------------------------------------------------------
  12 : 
  13 :   CRIPL = tabl ;
  14 :   CRIPL . PL = FAUX ;
  15 :   FPL = crit MOD SIG VAR MAT ;
  16 :   SSS = chan 'STRESSES' MOD ((exco 'SIGY' MAT 'SCAL') 
  17 :                'CHAN' 'TYPE' 'SCALAIRE');
  18 :   ERR = ER * SSS ;
  19 :   CRIPL . CR  = FPL masq egsupe ERR ;
  20 :   CRIPL . NPL = FPL masq egsupe somme ERR ;
  21 : 
  22 :   si ((CRIPL . NPL) > 0.) ;
  23 :     CRIPL . PL = VRAI ;
  24 :   finsi ;
  25 : 
  26 : FINP CRIPL ;
  27 :  
  28 :