Télécharger epthis.eso

Retour à la liste

Numérotation des lignes :

  1. C EPTHIS SOURCE FD218221 15/09/15 21:15:00 8626
  2.  
  3. C=======================================================================
  4. C= E P T H I S =
  5. C= ----------- =
  6. C= =
  7. C= Fonction : =
  8. C= ---------- =
  9. C= Calcul des deformations d'origine THERMIQUE dans le cas des =
  10. C= materiaux ISOTROPES. Sousprogramme appele par EPTHP (epthp.eso). =
  11. C= =
  12. C= Parametres : (E)=Entree (S)=Sortie =
  13. C= ------------ =
  14. C= MFR (E) Numero de la FORMULATION utilisee =
  15. C= MELE (E) Numero de l'element fini dans NOMTP (cf. CCHAMP) =
  16. C= VALMAT (E) Tableau des caracteristiques du materiau =
  17. C= NDEFS (E) Nombre de composantes de deformations =
  18. C= TEMP (E) Temperature au point de Gauss considere =
  19. C= THIF (E) | Temperatures des differentes couches =
  20. C= THM (E) | dans les cas des elements coques =
  21. C= THSU (E) | (couches INFErieure, MOYEnne, SUPErieure) =
  22. C= E3 (E) Excentrement dans la cas d'elements coques =
  23. C= EPAIST (E) Epaisseur dans le cas d'elements coques =
  24. C= RES (S) Tableau des valeurs des deformation thermiques =
  25. C= KERRE (S) Indicateur d'erreur si non nul =
  26. C=======================================================================
  27.  
  28. SUBROUTINE EPTHIS (MFR,MELE,VALMAT,NDEFS,TEMP,THIF,THM,THSU,
  29. . E3,EPAIST,RES,NPINT,KERRE)
  30.  
  31. IMPLICIT INTEGER(I-N)
  32. IMPLICIT REAL*8 (A-H,O-Z)
  33.  
  34. -INC CCOPTIO
  35. -INC CCREEL
  36.  
  37. DIMENSION VALMAT(*),RES(*)
  38.  
  39. KERRE=0
  40.  
  41. C 1 - Elements MASSIFS et POREUX
  42. C ================================
  43. IF (MFR.EQ.1.OR.MFR.EQ.33.OR.MFR.EQ.31.OR.MFR.EQ.63) THEN
  44. ALPDT=VALMAT(1)*TEMP
  45. C write(6,*) 'EPTHIS',MFR,ALPDT,VALMAT(1),TEMP
  46. C= 1.1 - Bidimensionnel PLAN Contraintes Planes
  47. IF (IFOUR.EQ.-2) THEN
  48. RES(1)=ALPDT
  49. RES(2)=ALPDT
  50. RES(3)=XZero
  51. RES(4)=XZero
  52. C= 1.2 - Bidimensionnel PLAN (DP/DPGE) et Axisymetrie
  53. ELSE IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) THEN
  54. RES(1)=ALPDT
  55. RES(2)=ALPDT
  56. RES(3)=ALPDT
  57. RES(4)=XZero
  58. C= 1.3 - Fourier et Tridimensionnel
  59. ELSE IF (IFOUR.EQ.1.OR.IFOUR.EQ.2) THEN
  60. RES(1)=ALPDT
  61. RES(2)=ALPDT
  62. RES(3)=ALPDT
  63. RES(4)=XZero
  64. RES(5)=XZero
  65. RES(6)=XZero
  66. C= 1.4 - Unidimensionnel (1D)
  67. ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.15.AND.MFR.EQ.1) THEN
  68. RES(1)=ALPDT
  69. IF (IFOUR.EQ.6) THEN
  70. RES(2)=XZero
  71. RES(3)=XZero
  72. ELSE IF (IFOUR.EQ.4.OR.IFOUR.EQ.8) THEN
  73. RES(2)=ALPDT
  74. RES(3)=XZero
  75. ELSE IF (IFOUR.EQ.5.OR.IFOUR.EQ.10.OR.IFOUR.EQ.13) THEN
  76. RES(2)=XZero
  77. RES(3)=ALPDT
  78. ELSE
  79. RES(2)=ALPDT
  80. RES(3)=ALPDT
  81. ENDIF
  82. C= 1.5 - Cas non prevus
  83. ELSE
  84. KERRE=19
  85. ENDIF
  86.  
  87. C 2 - Elements JOINTS 2D PLAN : JOI2 et JOI3
  88. C ============================================
  89. ELSE IF (MELE.EQ.85.OR.MELE.EQ.86) THEN
  90. ALPN=VALMAT(1)
  91. IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1
  92. & .OR.IFOUR.EQ.0) THEN
  93. RES(1)=XZero
  94. RES(2)=ALPN*TEMP
  95. ELSE
  96. KERRE=19
  97. ENDIF
  98.  
  99. C 3 - Element JOINT 3D : JOI4 ET JOT3
  100. C ===================================
  101. ELSE IF (MELE.EQ.88.OR.MELE.EQ.87) THEN
  102. ALPN=VALMAT(1)
  103. IF (IFOUR.EQ.2) THEN
  104. RES(1)=XZero
  105. RES(2)=XZero
  106. RES(3)=ALPN*TEMP
  107. ELSE
  108. KERRE=19
  109. ENDIF
  110.  
  111. C 4 - Elements COQUES MINCES sauf COQ2
  112. C ======================================
  113. ELSE IF (MELE.EQ.27.OR.MELE.EQ.28.OR.MELE.EQ.49.OR.
  114. . MELE.EQ.93) THEN
  115. ALPH=VALMAT(1)
  116. IF (NPINT.EQ.0) THEN
  117. THMO=(THSU+THIF+4.D0*THM)/6.D0
  118. THDI=THSU-THIF
  119. RES(1)=ALPH*THMO
  120. RES(2)=RES(1)
  121. IF (IFOUR.GT.0) THEN
  122. RES(3)=XZero
  123. RES(4)=ALPH*THDI/EPAIST
  124. RES(5)=RES(4)
  125. RES(6)=XZero
  126. ELSE
  127. RES(3)=ALPH*THDI/EPAIST
  128. RES(4)=RES(3)
  129. ENDIF
  130. ELSE
  131. IF (IFOUR.EQ.2) THEN
  132. RES(1)=ALPH*TEMP
  133. RES(2)=RES(1)
  134. RES(3)=XZero
  135. RES(4)=XZero
  136. ELSE
  137. KERRE=19
  138. ENDIF
  139. ENDIF
  140.  
  141. C 5 - Element COQ2
  142. C ==================
  143. ELSE IF (MELE.EQ.44) THEN
  144. ALPH=VALMAT(1)
  145. THMO=(THSU+THIF+4.D0*THM)/6.D0
  146. THDI=THSU-THIF
  147. IF (IFOUR.GT.0) THEN
  148. RES(1)=ALPH*THMO
  149. RES(2)=RES(1)
  150. RES(3)=XZero
  151. RES(4)=ALPH*THDI/EPAIST
  152. RES(5)=RES(4)
  153. RES(6)=XZero
  154. ELSE IF (IFOUR.EQ.-2) THEN
  155. RES(1)=ALPH*THMO
  156. RES(2)=XZero
  157. RES(3)=ALPH*THDI/EPAIST
  158. RES(4)=XZero
  159. ELSE
  160. RES(1)=ALPH*THMO
  161. RES(2)=RES(1)
  162. RES(3)=ALPH*THDI/EPAIST
  163. RES(4)=RES(3)
  164. ENDIF
  165.  
  166. C 6 - Elements COQUES EPAISSES
  167. C ==============================
  168. ELSE IF (MELE.EQ.41.OR.MELE.EQ.56) THEN
  169. ALPH=VALMAT(1)
  170. THG=(0.5*E3)*((E3+1.D0)*THSU+(E3-1.D0)*THIF)+
  171. . (1.D0-E3*E3)*THM
  172. RES(1)=ALPH*THG
  173. RES(2)=ALPH*THG
  174. RES(3)=XZero
  175. RES(4)=XZero
  176. RES(5)=XZero
  177.  
  178. C 7 - Elements POUTRES, TUYAUX et TIMO
  179. C ======================================
  180. ELSE IF (MELE.EQ.29.OR.MELE.EQ.42.OR.MELE.EQ.84) THEN
  181. RES(1)=VALMAT(1)*TEMP
  182.  
  183. C 8 - Elements BARRES, CERCES et POI1
  184. C =====================================
  185. ELSE IF (MELE.EQ.46.OR.MELE.EQ.95.OR.MELE.EQ.45.OR.
  186. . MELE.EQ.123.OR.MELE.EQ.124) THEN
  187. IF ((MELE.EQ.45.AND.IFOUR.NE.-3).OR.
  188. . (MELE.EQ.95.AND.IFOUR.NE.0)) THEN
  189. KERRE=86
  190. ELSE
  191. RES(1)=VALMAT(1)*TEMP
  192. ENDIF
  193.  
  194. C 9 - Autres cas non prevus
  195. C ===========================
  196. ELSE
  197. KERRE=19
  198. ENDIF
  199.  
  200. RETURN
  201. END
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  

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