Esope source files

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