Télécharger la plaquette Cast3M


fr en

Télécharger tufjda.eso

Retour à la liste

Numérotation des lignes : 

tufjda

  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
  36.          CALL ERREUR(519)
  37. C
  38. C   ERREUR
  39. C
  40.          GO TO 777
  41.       ENDIF
  42. C
  43. C STOCKAGE DES DEUX LISTES DE REELS J ET DELTAA
  44. C
  45.       DO 10 I=1,NBPOIX
  46.          DELTAA(I)=PROG(I)
  47.          XJLIS(I)=MLREE1.PROG(I)
  48.  10   CONTINUE
  49.       XJ1C=XJLIS(1)
  50.       XJTOT=XJP+((XK1**2)/YOUN)
  51.       IF(XJTOT.LT.XJ1C)THEN
  52. C
  53. C               PAS DE PROPAGATION
  54. C
  55.          THETA=THETA0
  56.       ELSE
  57. C
  58. C               PROPAGATION
  59. C
  60.          IBI=0
  61.          CALL TRACTI(DA,XJTOT,DELTAA,XJLIS,NBPOIX,2,IBI)
  62.          IF(IBI.EQ.1) THEN
  63.            KERRE=75
  64.          ENDIF
  65.          THETA=THETA0+((DA/RAYOM)/2.D0)
  66.       ENDIF
  67. C
  68.  777  CONTINUE
  69.       SEGDES MLREEL,MLREE1
  70.       SEGDES MEVOLL
  71.       RETURN
  72.       END
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales