Télécharger coupla.eso

Retour à la liste

Numérotation des lignes :

  1. C COUPLA SOURCE CB215821 16/04/21 21:16:06 8920
  2. C COUPLA SOURCE INSL 24/10/96
  3. SUBROUTINE COUPLA(EPS1,COURBE,IPCOUR,YOUN,SIGX,PENTE)
  4. C
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. DIMENSION COURBE(IPCOUR)
  8. C
  9. C **************************************
  10. C * CAS ELASTOPLASTIQUE ECROUISSABLE *
  11. C * COURBE INTRODUITE POINT PAR POINT *
  12. C * CALCUL DE LA PENTE DE LA COURBE *
  13. C **************************************
  14. C
  15. A1=1.D0
  16. IF(EPS1.LT.0.D0) A1=-1.D0
  17. EPSX=ABS(EPS1)
  18. SIGX1=ABS(SIGX)
  19. C
  20. IF(EPSX.LE.COURBE(4).AND.SIGX1.LT.COURBE(3)) THEN
  21. C --------------------------------------------
  22. C * PREMIER TRONCON DE LA COURBE DE TRACTION *
  23. C --------------------------------------------
  24. PENTE=YOUN
  25. SIGX=A1*YOUN*EPSX
  26. EPPX=0.D0
  27. ELSE
  28. C
  29. DO 1 I=6,IPCOUR,2
  30. C --------------------------------------------------------
  31. C * RECHERCHE DU POINT DE LA COURBE CORRESPONDANT A SET *
  32. C --------------------------------------------------------
  33. IF(EPSX.LE.COURBE(I)) THEN
  34. A=COURBE(I-3)
  35. B=COURBE(I-2)
  36. C=COURBE(I-1)
  37. D=COURBE(I)
  38. GOTO 3
  39. ENDIF
  40. 1 CONTINUE
  41. C
  42. PENTE=0.D0
  43. SIGX=0.D0
  44. EPPX=EPS1-SIGX/YOUN
  45. GOTO 4
  46. C ------------------------------------------------------------
  47. C * PENTEC = PENTE DE LA COURBE DE TRACTION (BP= A-B*PENTEC) *
  48. C ------------------------------------------------------------
  49. 3 CONTINUE
  50. PENTE=(C-A)/(D-B)
  51. SIGX=A1*(A+(PENTE*(EPSX-B)))
  52. EPPX=EPS1-SIGX/YOUN
  53. ENDIF
  54. 4 CONTINUE
  55. C
  56. RETURN
  57. END
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  

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