1 : C EPTHIS SOURCE AM 12/06/07 21:15:11 7390 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 et TUYAUX 179 : C ================================ 180 : ELSE IF (MELE.EQ.29.OR.MELE.EQ.42) 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 :
© Cast3M 2003 - All rights reserved.
Disclaimer