Procedures | Cast3M

   1 : * @DEFPL    PROCEDUR  CHAT      02/10/21    21:19:14     4462           
   2 : DEBP @DEFPL MMM*MMODEL MOD*TABLE MAT*MCHAML EPS*MCHAML SIGI*MCHAML
   3 :            EPSP*MCHAML VAR*MCHAML ERR*FLOTTANT E23D*MCHAML V23D*MCHAML ;
   4 : 
   5 : *-----------------------------------------------------------------------
   6 : * @DEFPL : calcul de l'increment de deformations plastiques et des
   7 : *         variables internes
   8 : *-----------------------------------------------------------------------
   9 : 
  10 :   DP  = tabl ;
  11 :   VV  = tabl ;
  12 :   SS  = tabl ;
  13 :   MAI = tabl ;
  14 :   NB = dime MOD ;
  15 : 
  16 :   si ((vale DIME) ega 3) ;
  17 :     P0 = 0. 0. 0. ; P1 = 0. 0. 1. ; ANGL = (360. / NB) ;
  18 :   finsi ;
  19 : 
  20 : *** boucle sur les colonnes du maillage 2D ou sur les portions du 
  21 : *** disque 3D
  22 : 
  23 :   repeter CAMEM NB ;
  24 :     EPP = redu EPS  (MOD . &CAMEM) ;
  25 :     SII = redu SIGI (MOD . &CAMEM) ;
  26 :     MMM = redu MAT  (MOD . &CAMEM) ;
  27 :     MAI . &CAMEM = extr (MOD . &CAMEM) 'MAIL' ;
  28 : 
  29 :     si (&CAMEM EGA 1) ;
  30 : 
  31 : *** initialisation des deformations plastiques et variables internes de
  32 : *** la premiere colonne (2D) ou de la premiere part (3D) du maillage
  33 : 
  34 :       DP . 1 = E23D ;
  35 :       VV . 1 = V23D ;
  36 :       N_EPSP = DP . 1 ;
  37 :       N_VAR  = VV . 1 ;
  38 : 
  39 : *** pour les colonnes et parts suivantes du maillage : translation ou 
  40 : *** rotation des valeurs trouvees pour la colonne ou la part precedente
  41 : 
  42 :     sinon ;
  43 : 
  44 : *** cas 2D : translation des champs
  45 :       si ((vale DIME) ega 2) ;
  46 :         SSCAM VVCAM DPCAM = 
  47 :           (SS . (&CAMEM - 1)) (VV . (&CAMEM - 1)) (DP . (&CAMEM - 1))
  48 :           dedu tran (MAI . (&CAMEM - 1)) (MAI . &CAMEM) ;
  49 :       finsi ;
  50 : 
  51 : *** cas 3D : rotation des champs
  52 :       si ((vale DIME) ega 3) ;
  53 :         SSCAM VVCAM DPCAM = 
  54 :           (SS . (&CAMEM - 1)) (VV . (&CAMEM - 1)) (DP . (&CAMEM - 1))
  55 :           dedu ANGL P0 P1 rota (MAI . (&CAMEM - 1)) (MAI . &CAMEM) ;
  56 :       finsi ;
  57 : 
  58 :       SS . &CAMEM = SSCAM ;
  59 :       VV . &CAMEM = VVCAM ;
  60 :       DP . &CAMEM = DPCAM ;
  61 :       SS . &CAMEM = SII - (SS . &CAMEM) ;
  62 : lm2 = extr MOD . &CAMEM deformat ;
  63 : lm1 = extr DP . &CAMEM comp ;
  64 : def_ut = changer comp DP . &CAMEM lm1 lm2 ;
  65 : def_ut = changer type def_ut 'DEFORMATIONS' ;
  66 :       SSS = SII + ((MOD . &CAMEM)
  67 : *                   (hook (MOD . &CAMEM) MMM) * (EPP - (DP . &CAMEM))) ;
  68 :                    (hook (MOD . &CAMEM) MMM) * (EPP - def_ut)) ;
  69 :       CC = tabl ;
  70 : 
  71 : *** calcul du critere de plasticite
  72 : 
  73 :       CC = @CRIPL (MOD . &CAMEM) SSS (VV . &CAMEM) MMM ERR ;
  74 : 
  75 : *** si elastique : def. plas. = def. plas. elements precedents
  76 : *** si plastique : def. plas. = def. plas. elements precedents 
  77 : ***                             + increment calcule avec ECOU
  78 : 
  79 :       si (CC . PL) ;
  80 : 
  81 : 
  82 : *        SSSS VAIN DEFP = ecou
  83 : *             (MOD . &CAMEM) (SS . &CAMEM) (VV . &CAMEM) EPP  MMM ;
  84 : *
  85 : com_var = 'EXTR' ('ZERO'  MOD . &CAMEM   'VARINTER') 'COMP' ;
  86 : com_dei = 'EXTR'      DP . &CAMEM  'COMP' ;
  87 :    cho2= comp  MOD . &CAMEM ( SS . &CAMEM et VV . &CAMEM ) 
  88 :           ( EPP  et MMM);
  89 :  DEFP = (changer ('EXCO' cho2 com_dei) 'TYPE' 'DEFORMATIONS'  ) 
  90 :           - dp .&camem ;
  91 : vain = 'CHANGER' ('EXCO' cho2 com_var) 'TYPE' 'VARIABLES INTERNES' ;
  92 : 
  93 :         DP . &CAMEM = (DP . &CAMEM) + DEFP ;
  94 :         VV . &CAMEM = VAIN ;
  95 :       finsi ;
  96 :       N_EPSP = N_EPSP et (DP . &CAMEM) ;
  97 :       N_VAR  = N_VAR  et (VV . &CAMEM) ;
  98 :     finsi ;
  99 : lm2 = extr MOD . &CAMEM deformat ;
 100 : lm1 = extr DP . &CAMEM comp ;
 101 : def_ut = changer comp DP . &CAMEM lm1 lm2 ;
 102 : def_ut = changer type def_ut 'DEFORMATIONS' ;
 103 : *    SS . &CAMEM = (MOD . &CAMEM)
 104 : *                  (hook (MOD . &CAMEM) MMM) * (DP . &CAMEM) ;
 105 :     SS . &CAMEM = (MOD . &CAMEM)
 106 :                   (hook (MOD . &CAMEM) MMM) * def_ut ;
 107 :   fin CAMEM ;
 108 : 
 109 : *** cas 2D : translation des champs de la derniere colonne sur la
 110 : ***          premiere, pour l'initialisation a l'iteration suivante
 111 : 
 112 :   si ((vale DIME) ega 2) ;
 113 :     E23D V23D = (DP . NB) (VV . NB)
 114 :                 dedu tran (MAI . NB) (MAI . 1) ;
 115 : 
 116 : *** cas 3D : rotation des champs de la derniere part sur la premiere,
 117 : ***          pour l'initialisation a l'iteration suivante
 118 : 
 119 :   sinon ;
 120 :     E23D V23D = (DP . NB) (VV . NB)
 121 :                 dedu ANGL P0 P1 rota (MAI . NB) (MAI . 1) ;
 122 :   finsi ;
 123 : 
 124 : FINP N_VAR N_EPSP E23D V23D ;
 125 :  
 126 : 
 127 : 
 128 : 
 129 :