Télécharger ccotra.eso

Retour à la liste

Numérotation des lignes :

  1. C CCOTRA SOURCE PV 17/12/08 21:15:28 9660
  2. SUBROUTINE CCOTRA(WRK52,WRK2,NCOURB,WRK53)
  3. C COTRAC SOURCE AM 97/12/23 21:41:23 3015
  4. c SUBROUTINE COTRAC(WRK0,WRK2,NCOURB,KERRE)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. -INC SMEVOLL
  8. -INC SMLREEL
  9. -INC DECHE
  10. SEGMENT/WRK2/(TRAC(LTRAC)*D)
  11. DATA NCOMAX/130/
  12. KERRE=0
  13. MEVOLL=nint(XMAT(5))
  14. IF(MEVOLL.EQ.0) THEN
  15. KERRE=37
  16. RETURN
  17. ENDIF
  18. SEGACT MEVOLL
  19. KEVOLL=IEVOLL(1)
  20. SEGACT KEVOLL
  21. MLREEL=IPROGX
  22. MLREE1=IPROGY
  23. * on ne desactive pas les segments pour reduire la contention sur esope en //
  24. *** SEGDES KEVOLL
  25. SEGACT MLREEL,MLREE1
  26. NBPOIX=PROG(/1)
  27. NBPOIY=MLREE1.PROG(/1)
  28. IF(NBPOIX.GT.NCOMAX+1) KERRE=31
  29. IF(NBPOIX.LT.3) KERRE=32
  30. IF(KERRE.EQ.0) GO TO 693
  31. 692 CONTINUE
  32. C
  33. C ERREUR
  34. C
  35. * on ne desactive pas les segments pour reduire la contention sur esope en //
  36. *** SEGDES MLREEL,MLREE1
  37. *** SEGDES MEVOLL
  38. GO TO 777
  39. 693 CONTINUE
  40. C
  41. C TESTS SUR LES POINTS DE LA COURBE ( A COMPLETER )
  42. C
  43. YUNG=XMAT(1)
  44. IF(YUNG.EQ.0.D0) THEN
  45. KERRE=34
  46. GO TO 692
  47. ENDIF
  48. PSIG=MLREE1.PROG(1)
  49. PEPS=PROG(1)
  50. IF(PSIG.NE.0.D0.OR.PEPS.NE.0.D0) THEN
  51. KERRE=35
  52. GO TO 692
  53. ENDIF
  54. C
  55. C VERIF DE LA LIMITE ELASTIQUE NON NULLE
  56. C
  57. PSIG =MLREE1.PROG(2)
  58. IF(PSIG.EQ.0.D0) THEN
  59. KERRE=30
  60. GO TO 692
  61. ENDIF
  62. PEPS=PROG(2)
  63. IF(PEPS.EQ.0.D0) THEN
  64. KERRE=35
  65. GO TO 692
  66. ENDIF
  67. PENTE=PSIG/PEPS
  68. RA=ABS(PENTE-YUNG)/YUNG
  69. IF(RA.GT.1.D-2) THEN
  70. KERRE=36
  71. GO TO 692
  72. ENDIF
  73. C
  74. C VERIF DE LA PENTE
  75. C
  76. DO 100 I=3,NBPOIX
  77. I1=I-1
  78. DEPS=PROG(I)-PROG(I1)
  79. IF(DEPS.EQ.0.D0) THEN
  80. KERRE=33
  81. GO TO 692
  82. ENDIF
  83. PENTE=(MLREE1.PROG(I)-MLREE1.PROG(I1))/DEPS
  84. IF(PENTE.GE.YUNG) THEN
  85. KERRE=33
  86. GO TO 692
  87. ENDIF
  88. 100 CONTINUE
  89. NCOURB=1
  90. TRAC(1)=MLREE1.PROG(2)
  91. TRAC(2)=0.D0
  92. DO 695 I=3,NBPOIX
  93. PEPS=PROG(I)
  94. PSIG=MLREE1.PROG(I)
  95. NCOURB=NCOURB+1
  96. TRAC(2*NCOURB-1)=PSIG
  97. TRAC(2*NCOURB)=PEPS-PSIG/YUNG
  98. 695 CONTINUE
  99. C
  100. * on ne desactive pas les segments pour reduire la contention sur esope en //
  101. *** SEGDES MLREEL,MLREE1
  102. *** SEGDES MEVOLL
  103. RETURN
  104. C
  105. 777 CONTINUE
  106. RETURN
  107. END
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  

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