Télécharger la plaquette Cast3M


fr en

Télécharger @defpl.procedur

Retour à la liste

Numérotation des lignes : 

  1. * @DEFPL    PROCEDUR  JK148537  23/08/21    21:15:12     11723          
  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 et (EPP*0.)) 
  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.  
  130.  
  131.  

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