Télécharger la plaquette Cast3M


fr en

Télécharger interf.eso

Retour à la liste

Numérotation des lignes : 

interf

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

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