Télécharger la plaquette Cast3M


fr en

Télécharger ecro.procedur

Retour à la liste

Numérotation des lignes : 

  1. * ECRO      PROCEDUR  SP204843  24/07/19    21:15:02     11967          
  2. *----------------------------------------------------------------------*
  3. * Procedure ECRO : conversion d'une courbe de traction convennelle     *
  4. * en courbe d'ecrouissage selon hypothese de calcul.                   *
  5. *                                                                      *
  6. * Syntaxe : voir notice.                                               *
  7. *                                                                      *
  8. *----------------------------------------------------------------------*
  9.  
  10. debp ECRO mot1*mot evol1*evolution ;
  11.  
  12. *----------------------- Verification EVOLUTION -----------------------*
  13.  
  14. * Verification evolution debute a l'origine :
  15. leps1    = extr evol1 absc ;
  16. lsig1    = extr evol1 ordo ;
  17. x1       = extr leps1 1 ;
  18. y1       = extr lsig1 1 ;
  19.  
  20. tol1     = (vale prec) * 10. ;
  21. si ((x1 neg 0. tol1) ou (y1 neg 0. tol1)) ;
  22.   ERREUR ': le 1er point de l''evolution n''est pas a l''origine.' ;
  23. fins ;
  24.  
  25. * Verification lecture du mot-cle :
  26. ipaslu1  = vrai ;
  27.  
  28. *---------------------- Cacul du module de Young ----------------------*
  29.  
  30. x2       = extr leps1 2 ;
  31. y2       = extr lsig1 2 ;
  32. dx       = x2 - x1 ;
  33. dY       = y2 - y1 ;
  34. si (dx ega 0. tol1) ;
  35.   ERREUR ': la pente initiale de la courbe est verticale.' ;
  36. sino ;
  37.   ym1      = dy / dx ;
  38. fins ;
  39.  
  40. *----------------------------- Option HPP -----------------------------*
  41.  
  42. si (ega mot1 'HPP') ;
  43.   ipaslu1  = faux ;
  44.  
  45.   lepe1    = lsig1 / ym1 ;
  46.   lepp1    = (leps1 - lepe1) enle 1 ;
  47.   lsmp1    = lsig1 enle 1 ;
  48.   ecro1    = evol roug manu epeq lepp1 smeq lsmp1 ;
  49.  
  50. fins ;
  51.  
  52. *---------------------------- Option TOTAL ----------------------------*
  53.  
  54. si (ega mot1 'TOTAL' 4) ;
  55.   ipaslu1  = faux ;
  56.  
  57.   epsy1    = x2 ;
  58.   sigy1    = y2 ;
  59.   sigy1    = sigy1 / (1. + epsy1) ;
  60.   epsy1    = 0.5 * epsy1 * epsy1 + epsy1 ;
  61.   ym1      = sigy1 / epsy1 ;
  62.   lsig1    = lsig1 / (1. + leps1) ;
  63.   leps1    = 0.5 * leps1 * leps1 + leps1 ;
  64.   lepe1    = lsig1 / ym1 ;
  65.   lepp1    = (leps1 - lepe1) enle 1 ;
  66.   lsmp1    = lsig1 enle 1 ;
  67.   ecro1    = evol roug manu epeq lepp1 smeq lsmp1 ;
  68.  
  69. fins ;
  70.  
  71. *------------------------- Option REACTUALISE -------------------------*
  72.  
  73. si (ega mot1 'REACTUALISE' 4) ;
  74.   ipaslu1  = faux ;
  75.  
  76.   epsy1    = x2 ;
  77.   sigy1    = y2 ;
  78.   sigy1    = sigy1 * (1. + epsy1) ;
  79.   epsy1    = log (1. + epsy1) ;
  80.   ym1      = sigy1 / epsy1 ;
  81.   lsig1    = lsig1 * (1. + leps1) ;
  82.   leps1    = log (1. + leps1) ;
  83.   lepe1    = lsig1 / ym1 ;
  84.   lepp1    = (leps1 - lepe1) enle 1 ;
  85.   lsmp1    = lsig1 enle 1 ;
  86.   ecro1    = evol roug manu epeq lepp1 smeq lsmp1 ;
  87.  
  88. fins ;
  89.  
  90. *-------------------- Pas le bon mot-cle => erreur --------------------*
  91.  
  92. si ipaslu1 ;
  93.   ERREUR ': on attent le mot-cle HPP, TOTAL ou REACTUALISE.' ;
  94. fins ;
  95.  
  96. finp ecro1 ym1 ;
  97.  
  98.  
  99.  

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