Télécharger tracti.eso

Retour à la liste

Numérotation des lignes :

  1. C TRACTI SOURCE PV 18/03/21 21:15:06 9783
  2. SUBROUTINE TRACTI (RESU,EPST,SIG,EPS,NCOURB,ICLE,IBI)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. DIMENSION SIG(*),EPS(*)
  6. C
  7. C ICLE = 1 ON CALCULE LA PENTE A LA COURBE DE TRACTION
  8. C ICLE = 2 ON CALCULE LA CONTRAINTE
  9. C IBI NON NUL : ERREUR : ON DEPASSE LE BOUT DE LA COURBE
  10.  
  11. C
  12. C ON INTERPOLE LINEAIREMENT
  13. C
  14. IBI=0
  15. *
  16. * Mars 2018
  17. * On autorise mmaintenant a depasser l'extremite de la courbe de traction
  18. * mais on teste que la contrainte est positive
  19. *
  20. *
  21. * TEST POUR VOIR SI ON DEPASSE EN EPSILON
  22. *
  23. IF(EPST.GT.1.1D0*EPS(NCOURB)) THEN
  24. ** IBI=1
  25. ** RETURN
  26. ENDIF
  27. *
  28. DO 11 I=2,NCOURB
  29. IP=I
  30. IF(EPS(I).GT.EPST) GO TO 13
  31. 11 CONTINUE
  32. 13 CONTINUE
  33. IP1=IP-1
  34. PENTE=(SIG(IP)-SIG(IP1))/(EPS(IP)-EPS(IP1))
  35. GO TO (1,2),ICLE
  36. RETURN
  37. 2 RESU=SIG(IP1)+PENTE*(EPST-EPS(IP1))
  38. if (resu.lt.0.d0) IBI=1
  39. RETURN
  40. END
  41.  
  42.  
  43.  
  44.  

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