Télécharger la plaquette Cast3M


fr en

Télécharger transi.eso

Retour à la liste

Numérotation des lignes : 

transi

  1. C TRANSI    SOURCE    CHAT      05/01/13    03:45:00     5004
  2. C     INTERACTION    SOURCE   NECHNECH   00/03/15
  3.       SUBROUTINE TRANSI(ITHHER,T0,TF,NSTRS,HI,BETJEF,NECH2)
  4. C
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7.       DIMENSION HI(4,4)
  8. *      COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,
  9. *     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
  10. *     &                ITER,ISIM,IBB,IGAU1,IZON
  11. *      COMMON /DBETNECH2/ATR,GTR,ALPH
  12.       SEGMENT BETJEF
  13.        REAL*8 AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
  14.      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  15.        INTEGER ICT,ICC,IMOD,IVIS,ITER,
  16.      &                ISIM,IBB,IGAU,IZON
  17.       ENDSEGMENT
  18.       SEGMENT NECH2
  19.        REAL*8 ATR,GTR,ALPH
  20.       ENDSEGMENT
  21. C
  22. *--------------------------------------------------------------------*
  23. *                                                                    *
  24. *    CALCUL DE  LA MATRICE DU FLUAGE TRANSITOIRE                     *
  25. *                    [HI]                                            *
  26. *                                                                    *
  27. *--------------------------------------------------------------------*
  28. C
  29. C       INITIALISATION
  30.         CALL ZERO(HI,4,4)
  31. C
  32. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  33. C                 CALCUL DE MATRICE [HI]
  34. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  35. C      CAS D'UN CALCUL THERMO-MECANIQUE
  36. C
  37.        IF(ITHHER.NE.0)THEN
  38.         DTETA=TF-T0
  39.         COEF=(DTETA*ALPH*ATR)/RB
  40. C
  41.         IF(DTETA.GE.0.D0)THEN
  42.          IF((IMOD.EQ.2).OR.(IMOD.EQ.4))THEN
  43.           HI(1,1)=COEF
  44.           HI(1,2)=0.D0
  45.           HI(1,3)=0.D0
  46.           HI(2,1)=0.D0
  47.           HI(2,2)=HI(1,1)
  48.           HI(2,3)=0.D0
  49.           HI(3,1)=0.D0
  50.           HI(3,2)=0.D0
  51.           HI(3,3)=HI(1,1)
  52.           HI(4,4)=0.D0
  53.          ENDIF
  54. C
  55.          IF((IMOD.EQ.1).OR.(IMOD.EQ.3))THEN
  56.           HI(1,1)=COEF
  57.           HI(1,2)=0.D0
  58.           HI(2,1)=0.D0
  59.           HI(2,2)=HI(1,1)
  60.           HI(3,3)=0.D0
  61.          ENDIF
  62.         ELSE
  63.          IF((IMOD.EQ.2).OR.(IMOD.EQ.4))THEN
  64.           HI(1,1)=0.D0
  65.           HI(2,2)=HI(1,1)
  66.           HI(3,3)=HI(1,1)
  67.           HI(4,4)=HI(1,1)
  68.          ENDIF
  69.          IF((IMOD.EQ.1).OR.(IMOD.EQ.3))THEN
  70.           HI(1,1)=0.D0
  71.           HI(2,2)=HI(1,1)
  72.           HI(3,3)=HI(1,1)
  73.          ENDIF
  74.         ENDIF
  75.       ENDIF
  76. C
  77. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  78. C       CAS D'UN CALCUL MECANIQUE
  79. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  80. C
  81.        IF(ITHHER.EQ.0)THEN
  82.         IF((IMOD.EQ.2).OR.(IMOD.EQ.4))THEN
  83.          HI(1,1)=0.D0
  84.          HI(2,2)=HI(1,1)
  85.          HI(3,3)=HI(1,1)
  86.          HI(4,4)=HI(1,1)
  87.        ENDIF
  88.        IF((IMOD.EQ.1).OR.(IMOD.EQ.3))THEN
  89.         HI(1,1)=0.D0
  90.         HI(2,2)=HI(1,1)
  91.         HI(3,3)=HI(1,1)
  92.        ENDIF
  93.       ENDIF
  94. C
  95.       RETURN
  96.       END
  97.  
  98.  
  99.  

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