Esope source files

   1 : C TUFCRI    SOURCE    CHAT      05/01/13    03:55:02     5004           
   2 :       SUBROUTINE TUFCRI(SIG0,VAR0,XMAT,XCAR,
   3 :      .     KERRE,CRITER)
   4 : C----------------------------------------------------------------------
   5 : C          TUYAU FISSURE  -  CRITERE
   6 : C     
   7 : C ENTREES
   8 : C     SIG0(NSTRS) = CONTRAINTES INITIALES
   9 : C     VAR0(NVARI) =  VARIABLES INTERNES   DEBUT
  10 : C     XMAT(NCOMAT)=  COMPOSANTES DE MATERIAU
  11 : C     XCAR(ICARA) =  CARACTERISTIQUES
  12 : C  SORTIES
  13 : C     CRITER  = LE CRITERE
  14 : C     KERRE   = 0    TOUT OK
  15 : C       30   LIMITE ELASTIQUE NULLE
  16 : C       34   MODULE D YOUNG NUL
  17 : C-----------------------------------------------------------------------
  18 :       IMPLICIT INTEGER(I-N)
  19 :       IMPLICIT REAL*8(A-H,O-Z)
  20 : C    Include contenant quelques constantes dont XPI :
  21 : -INC CCREEL
  22 :       DIMENSION SIG0(*),VAR0(*),XMAT(*),XCAR(*)
  23 : C     
  24 :       CRITER=0.D0
  25 :       KERRE=0
  26 :       YOUN=XMAT(1)
  27 :       IF(YOUN.EQ.0.D0) THEN
  28 :          KERRE=34
  29 :          RETURN
  30 :       ENDIF
  31 :       SIGMA0=XMAT(5)
  32 :       IF(SIGMA0.EQ.0.D0) THEN
  33 :          KERRE=30
  34 :          RETURN
  35 :       ENDIF
  36 :       EPAI=XCAR(2)
  37 :       RAYOM=XCAR(1)-EPAI/2.D0
  38 :       THETA=VAR0(3)*XPI/180.D0/2.D0
  39 :       XM=-SIG0(6)
  40 :       XP=SIG0(1)
  41 :       XM0=4.D0*SIGMA0*RAYOM*RAYOM*EPAI
  42 :       XP0=2.D0*SIGMA0*XPI*RAYOM*EPAI
  43 : C     
  44 :       CRI=2.D0*XP0*(1.D0-THETA/XPI)*XM
  45 :       CRI=CRI+XM0*SIN(THETA)*XP
  46 :       IF(CRI.GE.0.D0) THEN
  47 :          GO TO 1
  48 :       ELSE
  49 :          GO TO 2
  50 :       ENDIF
  51 :  1    CRITER=SIN(THETA)/2.D0
  52 :       CRITER=CRITER-COS((THETA+XPI*XP/XP0)/2.D0)
  53 :       CRITER=CRITER+XM/XM0
  54 :       RETURN
  55 :  2    CRITER=SIN(THETA)/2.D0
  56 :       CRITER=CRITER-COS((THETA-XPI*XP/XP0)/2.D0)
  57 :       CRITER=CRITER-XM/XM0
  58 :       RETURN
  59 :       END
  60 : 
  61 : 
  62 :  
  63 :  
  64 :