Télécharger la plaquette Cast3M


fr en

Télécharger dechat.eso

Retour à la liste

Numérotation des lignes : 

dechat

  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