Sources Esope | Cast3M

tufcri
C TUFCRI    SOURCE    CHAT      05/01/13    03:55:02     5004
      SUBROUTINE TUFCRI(SIG0,VAR0,XMAT,XCAR,
     .     KERRE,CRITER)
C----------------------------------------------------------------------
C          TUYAU FISSURE  -  CRITERE
C
C ENTREES
C     SIG0(NSTRS) = CONTRAINTES INITIALES
C     VAR0(NVARI) =  VARIABLES INTERNES   DEBUT
C     XMAT(NCOMAT)=  COMPOSANTES DE MATERIAU
C     XCAR(ICARA) =  CARACTERISTIQUES
C  SORTIES
C     CRITER  = LE CRITERE
C     KERRE   = 0    TOUT OK
C       30   LIMITE ELASTIQUE NULLE
C       34   MODULE D YOUNG NUL
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C    Include contenant quelques constantes dont XPI :
-INC CCREEL
      DIMENSION SIG0(*),VAR0(*),XMAT(*),XCAR(*)
C
      CRITER=0.D0
      KERRE=0
      YOUN=XMAT(1)
      IF(YOUN.EQ.0.D0) THEN
         KERRE=34
         RETURN
      ENDIF
      SIGMA0=XMAT(5)
      IF(SIGMA0.EQ.0.D0) THEN
         KERRE=30
         RETURN
      ENDIF
      EPAI=XCAR(2)
      RAYOM=XCAR(1)-EPAI/2.D0
      THETA=VAR0(3)*XPI/180.D0/2.D0
      XM=-SIG0(6)
      XP=SIG0(1)
      XM0=4.D0*SIGMA0*RAYOM*RAYOM*EPAI
      XP0=2.D0*SIGMA0*XPI*RAYOM*EPAI
C
      CRI=2.D0*XP0*(1.D0-THETA/XPI)*XM
      CRI=CRI+XM0*SIN(THETA)*XP
      IF(CRI.GE.0.D0) THEN
         GO TO 1
      ELSE
         GO TO 2
      ENDIF
 1    CRITER=SIN(THETA)/2.D0
      CRITER=CRITER-COS((THETA+XPI*XP/XP0)/2.D0)
      CRITER=CRITER+XM/XM0
      RETURN
 2    CRITER=SIN(THETA)/2.D0
      CRITER=CRITER-COS((THETA-XPI*XP/XP0)/2.D0)
      CRITER=CRITER-XM/XM0
      RETURN
      END