Esope source files

   1 : C CANG      SOURCE    CHAT      05/01/12    21:48:25     5004           
   2 :       SUBROUTINE CANG(TENS,ANG,YOUN)
   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 :       DIMENSION TENS(3)
   8 : C     
   9 : C-------------------------------------------
  10 : C     CETTE ROUTINE PERMET DE CALCULER     |
  11 : C     L ANGLE DE FISSURATION               |
  12 : C-------------------------------------------
  13 : C     
  14 :       ANG=0.D0
  15 :       SIREF=YOUN*1.D-4
  16 :       B=TENS(3)/SIREF
  17 :       A=TENS(1)-TENS(2)
  18 :       IF(ABS(B).LE.1.D-10) GOTO 10
  19 :       IF(A.EQ.0.D0) THEN
  20 :          X=1.D0
  21 :          IF(B.LT.0.D0) X=-1.D0
  22 :          ANG=X*XPI/4.D0
  23 :          GOTO 20
  24 :       ENDIF
  25 :       ANG=1.D0/2.D0*ATAN(2.D0*TENS(3)/A)
  26 :  10   IF(A.LT.0.D0) THEN
  27 :          ANG=ANG+XPI/2.D0
  28 :       ENDIF
  29 :       IF(ANG.GT.(XPI/2.D0)) THEN
  30 :          ANG=ANG-XPI
  31 :       ENDIF
  32 :  20   RETURN
  33 :       END
  34 : 
  35 : 
  36 :  
  37 :  
  38 :