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    PASCAL    21/05/06    21:15:02     10998          
  2.       SUBROUTINE CCOTRA(WRK52,WRK2,NCOURB,WRK53)
  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.       IMPLICIT INTEGER(I-N)
  14.       IMPLICIT REAL*8(A-H,O-Z)
  15. -INC SMEVOLL
  16. -INC SMLREEL
  17. -INC DECHE
  18. -INC CCREEL
  19.  
  20. -INC PPARAM
  21. -INC CCOPTIO
  22.       SEGMENT/WRK2/(TRAC(LTRAC)*D)
  23.       DATA NCOMAX/130/
  24.       KERRE=0
  25.       MEVOLL=nint(XMAT(5))
  26.  
  27.       SEGACT MEVOLL
  28.       KEVOLL=IEVOLL(1)
  29.       SEGACT KEVOLL
  30.       MLREEL=IPROGX
  31.       MLREE1=IPROGY
  32.  
  33.       SEGACT MLREEL,MLREE1
  34.       NBPOIX=PROG(/1)
  35.       NBPOIY=MLREE1.PROG(/1)
  36.       IF(NBPOIX.GT.NCOMAX+1) KERRE=31
  37.       IF(KERRE.NE.0) RETURN
  38.       IF(NBPOIX.LT.2) THEN
  39.         CALL ERREUR(1100) 
  40.         RETURN
  41.       ENDIF 
  42. C
  43. C   VERIFICATION DU MODULE DE YOUNG
  44. C
  45.       YUNG=XMAT(1)
  46.       IF(YUNG.LE.0.D0) THEN
  47.         REAERR(1)=YUNG
  48.         CALL ERREUR(1090)
  49.         KERRE=34
  50.         RETURN
  51.       ENDIF
  52. C
  53. C  VERIF DE LA LIMITE ELASTIQUE NON NULLE
  54. C
  55.       PSIG =MLREE1.PROG(1)
  56.       IF(PSIG.LT.0.D0) THEN
  57.         REAERR(1)=PSIG
  58.         CALL ERREUR(1091)
  59.         KERRE=30
  60.         RETURN
  61.       ENDIF
  62.       PEPS=PROG(1)
  63.       IF(ABS(PEPS).GT.XZPREC) THEN
  64.         REAERR(1)=PEPS
  65.         CALL ERREUR(1099)
  66.         KERRE=35
  67.         RETURN
  68.       ENDIF
  69. C
  70. C  VERIF DE LA PENTE
  71. C
  72.       DO 100 I=2,NBPOIX
  73.         I1=I-1
  74.         DEPS=PROG(I)-PROG(I1)
  75.         IF(DEPS.LE.0.D0) THEN
  76.           CALL ERREUR(1101)
  77.           KERRE=33
  78.           RETURN
  79.         ENDIF
  80.         PENTE=(MLREE1.PROG(I)-MLREE1.PROG(I1))/DEPS
  81.         IF(ABS(PENTE-YUNG).LT.(1.D-3*YUNG)) THEN
  82. C          REAERR(1)=PENTE
  83. C          REAERR(2)=YUNG
  84.           CALL ERREUR(1102)
  85.           KERRE=33
  86.           RETURN
  87.         ENDIF
  88.  100  CONTINUE
  89. C
  90. C  DONNEE DE LA COURBE SOUS FORME D'UNE LISTE 
  91. C
  92.       NCOURB = 1
  93.       TRAC(1)= MLREE1.PROG(1)
  94.       TRAC(2)= 0.D0
  95.       DO 695 I=2,NBPOIX
  96.         PEPS = PROG(I)
  97.         PSIG = MLREE1.PROG(I)
  98.         NCOURB = NCOURB+1
  99.         TRAC(2*NCOURB-1) = PSIG
  100.         TRAC(2*NCOURB) = PEPS
  101.   695 CONTINUE
  102.  
  103.       END
  104.  
  105.  
  106.  
  107.  
  108.  

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