Télécharger ccotr4.eso

Retour à la liste

Numérotation des lignes :

  1. C CCOTR4 SOURCE PASCAL 19/07/26 21:15:01 10198
  2. SUBROUTINE CCOTR4(WRK52,WRK2,NCOURB,WRK53)
  3. C CCOTRA SOURCE CB215821 18/12/04 21:15:10 10020
  4. C COTRAC SOURCE AM 97/12/23 21:41:23 3015
  5. c SUBROUTINE COTRAC(WRK0,WRK2,NCOURB,KERRE)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. -INC SMEVOLL
  9. -INC SMLREEL
  10. -INC DECHE
  11. -INC CCREEL
  12. -INC CCOPTIO
  13. SEGMENT/WRK2/(TRAC(LTRAC)*D)
  14. DATA NCOMAX/130/
  15. KERRE=0
  16. MEVOLL=nint(XMAT(5))
  17.  
  18. SEGACT MEVOLL
  19. KEVOLL=IEVOLL(1)
  20. SEGACT KEVOLL
  21. MLREEL=IPROGX
  22. MLREE1=IPROGY
  23.  
  24. SEGACT MLREEL,MLREE1
  25. NBPOIX=PROG(/1)
  26. NBPOIY=MLREE1.PROG(/1)
  27. IF(NBPOIX.GT.NCOMAX+1) KERRE=31
  28. IF(NBPOIX.LT.3) KERRE=32
  29. IF(KERRE.NE.0) RETURN
  30. C
  31. C TESTS SUR LES POINTS DE LA COURBE ( A COMPLETER )
  32. C
  33. YUNG=XMAT(1)
  34. IF(YUNG .LE. 0.D0) THEN
  35. REAERR(1)=YUNG
  36. CALL ERREUR(1090)
  37. KERRE=34
  38. RETURN
  39. ENDIF
  40. PSIG=MLREE1.PROG(1)
  41. PEPS=PROG(1)
  42. IF(PSIG.NE.0.D0 .OR. PEPS.NE.0.D0) THEN
  43. REAERR(1)=PEPS
  44. REAERR(2)=PSIG
  45. CALL ERREUR(1092)
  46. KERRE=35
  47. RETURN
  48. ENDIF
  49. C
  50. C VERIF DE LA LIMITE ELASTIQUE NON NULLE
  51. C
  52. PSIG =MLREE1.PROG(2)
  53. IF(PSIG.EQ.0.D0) THEN
  54. REAERR(1)=PSIG
  55. CALL ERREUR(1091)
  56. KERRE=30
  57. RETURN
  58. ENDIF
  59. PEPS=PROG(2)
  60. IF(PEPS.EQ.0.D0) THEN
  61. REAERR(1)=PEPS
  62. CALL ERREUR(1089)
  63. KERRE=35
  64. RETURN
  65. ENDIF
  66. PENTE=PSIG/PEPS
  67. RA=ABS(PENTE-YUNG)/YUNG
  68. IF(RA.GT.1.D-3) THEN
  69. REAERR(1)=PENTE
  70. REAERR(2)=YUNG
  71. CALL ERREUR(1088)
  72. KERRE=36
  73. RETURN
  74. ENDIF
  75. C
  76. C VERIF DE LA PENTE
  77. C
  78. DO 100 I=3,NBPOIX
  79. I1=I-1
  80. DEPS=PROG(I)-PROG(I1)
  81. IF(DEPS.LE.0.D0) THEN
  82. CALL ERREUR(1087)
  83. KERRE=33
  84. RETURN
  85. ENDIF
  86. PENTE=(MLREE1.PROG(I)-MLREE1.PROG(I1))/DEPS
  87. IF(PENTE.GT.1.01*YUNG) THEN
  88. REAERR(1)=PENTE
  89. REAERR(2)=YUNG
  90. CALL ERREUR(1088)
  91. KERRE=33
  92. RETURN
  93. ENDIF
  94. 100 CONTINUE
  95. NCOURB=1
  96. TRAC(1)=MLREE1.PROG(2)
  97. TRAC(2)=0.D0
  98. DO 695 I=3,NBPOIX
  99. PEPS=PROG(I)
  100. PSIG=MLREE1.PROG(I)
  101. NCOURB=NCOURB+1
  102. TRAC(2*NCOURB-1)=PSIG
  103. TRAC(2*NCOURB)=PEPS-PSIG/YUNG
  104. 695 CONTINUE
  105.  
  106. END
  107.  
  108.  
  109.  

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