Télécharger la plaquette Cast3M


fr en

Télécharger coupla.eso

Retour à la liste

Numérotation des lignes : 

coupla

  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