Télécharger @defpl.procedur

Retour à la liste

Numérotation des lignes :

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

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