Télécharger cotrac.eso

Retour à la liste

Numérotation des lignes :

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

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