Esope source files

   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
  29 :         CALL ZERO(HI,4,4)
  30 : C
  31 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  32 : C                 CALCUL DE MATRICE [HI]
  33 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  34 : C      CAS D'UN CALCUL THERMO-MECANIQUE
  35 : C
  36 :        IF(ITHHER.NE.0)THEN
  37 :         DTETA=TF-T0
  38 :         COEF=(DTETA*ALPH*ATR)/RB
  39 : C        
  40 :         IF(DTETA.GE.0.D0)THEN
  41 :          IF((IMOD.EQ.2).OR.(IMOD.EQ.4))THEN
  42 :           HI(1,1)=COEF
  43 :           HI(1,2)=0.D0
  44 :           HI(1,3)=0.D0
  45 :           HI(2,1)=0.D0
  46 :           HI(2,2)=HI(1,1)
  47 :           HI(2,3)=0.D0
  48 :           HI(3,1)=0.D0
  49 :           HI(3,2)=0.D0
  50 :           HI(3,3)=HI(1,1)
  51 :           HI(4,4)=0.D0
  52 :          ENDIF
  53 : C
  54 :          IF((IMOD.EQ.1).OR.(IMOD.EQ.3))THEN
  55 :           HI(1,1)=COEF
  56 :           HI(1,2)=0.D0
  57 :           HI(2,1)=0.D0
  58 :           HI(2,2)=HI(1,1)
  59 :           HI(3,3)=0.D0
  60 :          ENDIF
  61 :         ELSE
  62 :          IF((IMOD.EQ.2).OR.(IMOD.EQ.4))THEN
  63 :           HI(1,1)=0.D0
  64 :           HI(2,2)=HI(1,1)
  65 :           HI(3,3)=HI(1,1)
  66 :           HI(4,4)=HI(1,1)
  67 :          ENDIF
  68 :          IF((IMOD.EQ.1).OR.(IMOD.EQ.3))THEN
  69 :           HI(1,1)=0.D0
  70 :           HI(2,2)=HI(1,1)
  71 :           HI(3,3)=HI(1,1)
  72 :          ENDIF 
  73 :         ENDIF
  74 :       ENDIF
  75 : C
  76 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  77 : C       CAS D'UN CALCUL MECANIQUE
  78 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  79 : C
  80 :        IF(ITHHER.EQ.0)THEN
  81 :         IF((IMOD.EQ.2).OR.(IMOD.EQ.4))THEN
  82 :          HI(1,1)=0.D0
  83 :          HI(2,2)=HI(1,1)
  84 :          HI(3,3)=HI(1,1)
  85 :          HI(4,4)=HI(1,1)
  86 :        ENDIF
  87 :        IF((IMOD.EQ.1).OR.(IMOD.EQ.3))THEN
  88 :         HI(1,1)=0.D0
  89 :         HI(2,2)=HI(1,1)
  90 :         HI(3,3)=HI(1,1)
  91 :        ENDIF 
  92 :       ENDIF 
  93 : C      
  94 :       RETURN
  95 :       END 
  96 :  
  97 :