Télécharger interf.eso

Retour à la liste

Numérotation des lignes :

  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