Esope source files

   1 : C TUFIC1    SOURCE    CHAT      05/01/13    03:55:18     5004           
   2 :       SUBROUTINE TUFIC1(XM,XP,THETA,IRETOU,PRECIS,XM0,XP0)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C    Include contenant quelques constantes dont XPI :
   6 : -INC CCREEL
   7 : C=======================================================================
   8 : C      POUR SAVOIR SI ON EST ELASTIQUE OU PAS
   9 : C     
  10 : C         ENTREES:XM,XP      ETAT DE CONTRAINTES
  11 : C                 THETA    ANGLE DE LA FISSURE
  12 : C     
  13 : C         SORTIE:IRETOU    0 SI ELASTIQUE,1 SINON
  14 : C=======================================================================
  15 :       IRETOU=0
  16 : C     
  17 : C      CALCUL DE LA VALEUR DE LA FONCTION DE CHARGE EN OUVERTURE
  18 : C     
  19 :       CRIT=SIN(THETA)/2.D0
  20 :       CRIT=CRIT-COS((THETA+XPI*XP/XP0)/2.D0)
  21 :       CRIT=CRIT+XM/XM0
  22 :       IF(CRIT.GE.PRECIS) THEN
  23 :          IRETOU=1
  24 :          RETURN
  25 :       ENDIF
  26 : C     
  27 : C      CALCUL DE LA VALEUR DE LA FONCTION DE CHARGE EN FERMETURE
  28 : C     
  29 :       CRIT=SIN(THETA)/2.D0
  30 :       CRIT=CRIT-COS((THETA-XPI*XP/XP0)/2.D0)
  31 :       CRIT=CRIT-XM/XM0
  32 :       IF(CRIT.GE.PRECIS) THEN
  33 :          IRETOU=1
  34 :          RETURN
  35 :       ENDIF
  36 :       RETURN
  37 :       END
  38 : 
  39 : 
  40 :  
  41 :  
  42 :