Télécharger dechat.eso

Retour à la liste

Numérotation des lignes :

  1. C DECHAT SOURCE CB215821 16/04/21 21:16:20 8920
  2. C DECHAT SOURCE INSL 24/10/96
  3. SUBROUTINE DECHAT(EPSRX,STRNX,STRNRX,SIGMRX,TANG,S1X,EDT,EDC,
  4. 1 EPST,DEFR,IFISU,IPLA,EQSTR,RTM,EPSC,IREFE,EX,RB,ALPHA,EPSU,
  5. 2 EPO1,ICAL,IBB1,IGAU1,IDIR,PASDT)
  6. C
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8(A-H,O-Z)
  9. CHARACTER*8 PASDT
  10. C
  11. C CE SOUS-PROGRAMME GERE LA REFERMETURE DES FISSURES.
  12. C CALCUL DE LA NOUVELLE PENTE EDT (RESTAURATION PROGRESSIVE
  13. C DE LA RAIDEUR DE LA FISSURE
  14. C
  15. IREFE=0
  16. RBT=ALPHA*RB
  17. REFERM=3.D0*EPSU
  18. IF(ABS(EDC).LT.0.1D-06) EDC=0.D0
  19. IF(EDC.NE.0.D0) REFERM=REFERM+DEFR
  20. IF(STRNRX.GT.EPSRX) RETURN
  21. IF(STRNRX.GT.REFERM) THEN
  22. IREFE=1
  23. S1X=0.D0
  24. RETURN
  25. ENDIF
  26. IFISU=0
  27. IF(EDC.NE.0.D0) EPST=DEFR-RBT/EDC
  28. C
  29. C REFERMETURE POUR UN POINT INITIALLEMENT TENDU
  30. C
  31. C CALCUL DU MODULE EDT LORS DE LA REFERMETURE
  32. C
  33. C
  34. IF(STRNRX.GT.EPST) THEN
  35. IF(EPSRX.LE.REFERM)THEN
  36. EDT=(RBT+RTM)/(EPSRX-EPST)
  37. IF(IPLA.EQ.2.AND.EQSTR.LT.RBT) EDT=EQSTR/(EPSRX-EPST)
  38. S1X=EDT*STRNX+RTM
  39. ELSE
  40. EDT=RBT/(REFERM-EPST)
  41. IF(IPLA.EQ.2.AND.EQSTR.LT.RBT) EDT=EQSTR/(REFERM-EPST)
  42. S1X=EDT*(STRNRX-REFERM)
  43. ENDIF
  44. TANG=EDT
  45. ELSE
  46. IF(STRNRX.LT.EPSC) THEN
  47. EPEQ=ABS(STRNRX)
  48. CALL CDCY(EPEQ,SEQ,IPLA,TANG,EX,RB,EMAX,EPO1,ICAL)
  49. C
  50. C ACTUALISATION DE EQSTR
  51. C
  52. EQSTR=SEQ
  53. S1X=-SEQ
  54. ELSE
  55. C
  56. C ON EST SUR LA PENTE EDC
  57. C
  58. S1X=EDC*(STRNRX-DEFR)
  59. TANG=EDC
  60. ENDIF
  61. ENDIF
  62. IF(EDT.LT.0.D0) THEN
  63. WRITE(*,*) ' !!!! ATTENTION DANS DECHART EDT < 0 APPELE PAR'
  64. & ,PASDT
  65. WRITE(*,94) IBB1,IGAU1,IDIR
  66. 94 FORMAT(' ELEME=',I3,'POINT=',I2,' IDIR=',I2)
  67. WRITE(*,*)'RTM=',RTM,'RBT=',RBT,'IPLA=',IPLA,'IFISU=',IFISU
  68. WRITE(*,*) 'EPSRX= ',EPSRX,' SIGMRX= ',SIGMRX
  69. WRITE(*,*) 'STRNRX= ',STRNRX,' EPSC= ',EPSC
  70. WRITE(*,*) 'EQSTR= ',EQSTR,' EPST= ',EPST,' DEFR=',DEFR
  71. WRITE(*,*) 'ERREUR EDT<0 ! EDT=',EDT,' EDC=',EDC
  72. STOP
  73. ENDIF
  74. RETURN
  75. END
  76.  
  77.  
  78.  
  79.  
  80.  

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