Esope source files

   1 : C CRIT1     SOURCE    CHAT      05/01/12    22:30:38     5004           
   2 :       FUNCTION CRIT1(XM,XP,THETA,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     CALCULE LA VALEUR DE LA FONCTION DE CHARGE
   9 : C     
  10 : C     ENTREES:XM      VALEUR DU MOMENT
  11 : C             XP      VALEUR DE L EFFORT NORMAL
  12 : C             THETA   ANGLE DE LA FISSURE
  13 : C     
  14 : C     SORTIE:CRIT1    VALEUR DE LA FONCTION
  15 : C=======================================================================
  16 :       CRI=2.D0*XP0*(1.D0-THETA/XPI)*XM
  17 :       CRI=CRI+XM0*SIN(THETA)*XP
  18 :       IF(CRI.GE.XZERO) THEN
  19 :          GO TO 1
  20 :       ELSE
  21 :          GO TO 2
  22 :       ENDIF
  23 :  1    CONTINUE
  24 :       CRIT1=SIN(THETA)/2.D0
  25 :       CRIT1=CRIT1-COS((THETA+XPI*XP/XP0)/2.D0)
  26 :       CRIT1=CRIT1+XM/XM0
  27 :       RETURN
  28 :  2    CONTINUE
  29 :       CRIT1=SIN(THETA)/2.D0
  30 :       CRIT1=CRIT1-COS((THETA-XPI*XP/XP0)/2.D0)
  31 :       CRIT1=CRIT1-XM/XM0
  32 :       RETURN
  33 :       END
  34 : 
  35 : 
  36 :  
  37 :  
  38 :