Télécharger la plaquette Cast3M


fr en

Télécharger ccotra.eso

Retour à la liste

Numérotation des lignes : 

ccotra

  1. C CCOTRA    SOURCE    OF166741  25/09/30    21:15:03     12371          
  2.  
  3. C-----------------------------------------------------------------------
  4. C
  5. C     CETTE SUBROUTINE RENVOIE LA RELATION CONTRAINTE-DEFORMATION
  6. C     DANS LE SEGMENT WRK2 A PARTIR DE LA COURBE D'ECROUISSAGE,
  7. C     I.E. LA COURBE CONTRAINTE VS DEFORMATION PLASTIQUE.
  8. C
  9. C     LE PREMIER POINT DE LA COURBE FOURNIT SIGY,
  10. C     LA DEFORMATION PLASTIQUE ASSOCIEE EST NULLE.
  11. C
  12. C-----------------------------------------------------------------------
  13.       SUBROUTINE CCOTRA(WRK52,WRK2,NCOURB,WRK53)
  14.  
  15.       IMPLICIT INTEGER(I-N)
  16.       IMPLICIT REAL*8(A-H,O-Z)
  17.  
  18. -INC PPARAM
  19. -INC CCOPTIO
  20. -INC CCREEL
  21.  
  22. -INC SMEVOLL
  23. -INC SMLREEL
  24. -INC DECHE
  25.  
  26.       SEGMENT/WRK2/(TRAC(LTRAC)*D)
  27.       DATA NCOMAX/130/
  28.  
  29.       KERRE = 0
  30.  
  31.       MEVOLL = nint(XMAT(5))
  32.       SEGACT MEVOLL
  33.       KEVOLL = mevoll.IEVOLL(1)
  34. c//      SEGDES,MEVOLL
  35.  
  36.       SEGACT KEVOLL
  37.       MLREEL = kevoll.IPROGX
  38.       MLREE1 = kevoll.IPROGY
  39. c//      SEGDES,KEVOLL
  40.  
  41.       SEGACT MLREEL,MLREE1
  42.       NBPOIX = mlreel.PROG(/1)
  43.       NBPOIY = MLREE1.PROG(/1)
  44.       IF (NBPOIX.GT.NCOMAX+1) THEN
  45.         KERRE=31
  46.         RETURN
  47.       ELSE IF (NBPOIX.LT.2) THEN
  48.         CALL ERREUR(1100) 
  49.         KERRE=1100
  50.         RETURN
  51.       ENDIF 
  52. c*C
  53. c*C   VERIFICATION DU MODULE DE YOUNG
  54. c*C
  55. c*      YUNG=XMAT(1)
  56. c*      IF (YUNG.LE.0.D0) THEN
  57. c*        REAERR(1)=YUNG
  58. c*        CALL ERREUR(1090)
  59. c*        KERRE=34
  60. c*        RETURN
  61. c*      ENDIF
  62. C
  63. C  VERIFICATIONS PREMIER POINT DE LA COURBE :
  64. C  DEFORMATION PLASTIQUE NULLE & LIMITE ELASTIQUE NON NULLE
  65. C
  66.       PEPS = mlreel.PROG(1)
  67.       PSIG = MLREE1.PROG(1)
  68.       IF (PSIG.LE.0.D0) THEN
  69.         REAERR(1)=PSIG
  70.         CALL ERREUR(1091)
  71.         KERRE=30
  72.         RETURN
  73.       ENDIF
  74.       IF (ABS(PEPS).GT.XZPREC) THEN
  75.         REAERR(1)=PEPS
  76.         CALL ERREUR(1099)
  77.         KERRE=35
  78.         RETURN
  79.       ENDIF
  80. C
  81. C  DONNEE DE LA COURBE SOUS FORME D'UNE LISTE 
  82. c*C  VERIFICATION DE LA PENTE : test supprime
  83. C
  84.       NCOURB = 1
  85.       TRAC(1) = PSIG
  86.       TRAC(2) = 0.D0
  87.  
  88.       DO I = 2, NBPOIX
  89.         I1=I-1
  90.         PEPS = mlreel.PROG(I)
  91.         PSIG = mlree1.PROG(I)
  92.  
  93.         DEPS = PEPS - mlreel.PROG(I1)
  94.         IF (DEPS.LE.0.D0) THEN
  95.           CALL ERREUR(1101)
  96.           KERRE=33
  97.           RETURN
  98.         ENDIF
  99. c*        PENTE=(PSIG-MLREE1.PROG(I1))/DEPS
  100. c*        IF(ABS(PENTE-YUNG).LT.(1.D-3*YUNG)) THEN
  101. c*          REAERR(1)=PENTE
  102. c*          REAERR(2)=YUNG
  103. c*          CALL ERREUR(1102)
  104. c*          KERRE=33
  105. c*          RETURN
  106. c*        ENDIF
  107.  
  108.         NCOURB = NCOURB+1
  109.         TRAC(2*NCOURB-1) = PSIG
  110.         TRAC(2*NCOURB)   = PEPS
  111.       ENDDO
  112.  
  113. c      RETURN
  114.       END
  115.  
  116.  
  117.  

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