Esope source files

   1 : C TUFJDA    SOURCE    CHAT      05/01/13    03:56:29     5004           
   2 :       SUBROUTINE TUFJDA(THETA,XJP,RAYOM,JDA,THETA0,XK1,YOUN,KERRE)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 :       DIMENSION DELTAA(130),XJLIS(130)
   6 : C     
   7 : C-----------------------------------------------------------------
   8 : C              PROPAGATION TUYAU FISSURE
   9 : C  CE SOUS-PROGRAMME EST APPELE PAR TUFPLA
  10 : C     
  11 : C  ENTREES:  XJP          VALEUR DE J PLASTIQUE AU POINT CONSIDERE
  12 : C            XJE          VALEUR DE J ELASTIQUE AU POINT CONSIDERE
  13 : C            RAYOM       RAYON MOYEN DU TUYAU
  14 : C            XJDA         COURBE (J-DA) DU MATERIAU
  15 : C            THETA0      ANGLE DE FISSURE INITIAL
  16 : C     
  17 : C  SORTIES:  THETA       ANGLE DE LA FISSURE DANS L'ETAT DE CHARGE
  18 : C                        CONSIDERE
  19 : C-----------------------------------------------------------------
  20 : C     
  21 : -INC SMEVOLL
  22 : -INC SMLREEL
  23 : C     
  24 :       DATA NCOMAX/130/
  25 :       KERRE=0
  26 :       MEVOLL=JDA
  27 :       SEGACT MEVOLL
  28 :       KEVOLL=IEVOLL(1)
  29 :       SEGACT KEVOLL
  30 :       MLREEL=IPROGX
  31 :       MLREE1=IPROGY
  32 :       SEGDES KEVOLL
  33 :       SEGACT MLREEL,MLREE1
  34 :       NBPOIX=PROG(/1)
  35 :       IF(NBPOIX.GT.NCOMAX+1) THEN
  35 :          CALL ERREUR(519)
  36 : C     
  37 : C   ERREUR
  38 : C     
  39 :          GO TO 777
  40 :       ENDIF
  41 : C     
  42 : C STOCKAGE DES DEUX LISTES DE REELS J ET DELTAA
  43 : C     
  44 :       DO 10 I=1,NBPOIX
  45 :          DELTAA(I)=PROG(I)
  46 :          XJLIS(I)=MLREE1.PROG(I)
  47 :  10   CONTINUE
  48 :       XJ1C=XJLIS(1)
  49 :       XJTOT=XJP+((XK1**2)/YOUN)
  50 :       IF(XJTOT.LT.XJ1C)THEN
  51 : C     
  52 : C               PAS DE PROPAGATION
  53 : C     
  54 :          THETA=THETA0
  55 :       ELSE
  56 : C     
  57 : C               PROPAGATION
  58 : C     
  59 :          IBI=0
  59 :          CALL TRACTI(DA,XJTOT,DELTAA,XJLIS,NBPOIX,2,IBI)
  60 :          IF(IBI.EQ.1) THEN
  61 :            KERRE=75
  62 :          ENDIF
  63 :          THETA=THETA0+((DA/RAYOM)/2.D0)
  64 :       ENDIF
  65 : C     
  66 :  777  CONTINUE
  67 :       SEGDES MLREEL,MLREE1
  68 :       SEGDES MEVOLL
  69 :       RETURN
  70 :       END
  71 : 
  72 : 
  73 : 
  74 : 
  75 : 
  76 :  
  77 :  
  78 :