1 : C THETIS SOURCE AM 12/06/07 21:16:09 7390 2 : 3 : C======================================================================= 4 : C= T H E T I S = 5 : C= ----------- = 6 : C= = 7 : C= Fonction : = 8 : C= ---------- = 9 : C= Calcul des contraintes 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= NSTRS (E) Nombre de composantes de contraintes = 18 : C= NPINT (E) Nombre de points d'integration pour les coques = 19 : C= TEMP (E) Temperature au point de Gauss considere = 20 : C= THIF (E) | Temperatures des differentes couches = 21 : C= THM (E) | dans les cas des elements coques = 22 : C= THSU (E) | (couches INFErieure, MOYEnne, SUPErieure) = 23 : C= E3 (E) Excentrement dans la cas d'elements coques = 24 : C= SD (E) Section dans le cas d'elements barres = 25 : C= EPAIST (E) Epaisseur dans le cas d'elements coques = 26 : C= RES (S) Tableau des valeurs des contraintes thermiques = 27 : C= KERRE (S) Indicateur d'erreur si non nul = 28 : C======================================================================= 29 : 30 : SUBROUTINE THETIS (MFR,MELE,VALMAT,NSTRS,NPINT,TEMP,THIF, 31 : . THM,THSU,E3,SD,EPAIST,RES,KERRE) 32 : 33 : IMPLICIT INTEGER(I-N) 34 : IMPLICIT REAL*8 (A-H,O-Z) 35 : 36 : -INC CCOPTIO 37 : -INC CCHAMP 38 : 39 : DIMENSION VALMAT(*),RES(*) 40 : 41 : KERRE=0 42 : 43 : C 1 - Elements MASSIFS et POREUX 44 : C ================================ 45 : IF (MFR.EQ.1.OR.MFR.EQ.33.OR.MFR.EQ.31.OR.MFR.EQ.63) THEN 46 : YOU=VALMAT(1) 47 : XNU=VALMAT(2) 48 : ALPDT=VALMAT(3)*TEMP 49 : C= 1.1 - Bidimensionnel PLAN Contraintes Planes 50 : IF (IFOUR.EQ.-2) THEN 51 : AUX=YOU*ALPDT/(1.D0-XNU) 52 : RES(1)=AUX 53 : RES(2)=AUX 54 : RES(3)=0. 55 : RES(4)=0. 56 : C= 1.2 - Bidimensionnel PLAN (DP/DPGE) et Axisymetrie 57 : ELSE IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) THEN 58 : AUX=YOU*ALPDT/(1.D0-2.D0*XNU) 59 : RES(1)=AUX 60 : RES(2)=AUX 61 : RES(3)=AUX 62 : RES(4)=0.D0 63 : C= 1.3 - Fourier et Tridimensionnel 64 : ELSE IF (IFOUR.EQ.1.OR.IFOUR.EQ.2) THEN 65 : AUX=YOU*ALPDT/(1.D0-2.D0*XNU) 66 : RES(1)=AUX 67 : RES(2)=AUX 68 : RES(3)=AUX 69 : RES(4)=0. 70 : RES(5)=0. 71 : RES(6)=0. 72 : C= 1.4 - Unidimensionnel (1D) 73 : ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.15.AND.MFR.EQ.1) THEN 74 : IF (IFOUR.EQ.6) THEN 75 : RES(1)=YOU*ALPDT 76 : RES(2)=0. 77 : RES(3)=0. 78 : ELSE IF (IFOUR.EQ.4.OR.IFOUR.EQ.8) THEN 79 : AUX=YOU*ALPDT/(1.-XNU) 80 : RES(1)=AUX 81 : RES(2)=AUX 82 : RES(3)=0. 83 : ELSE IF (IFOUR.EQ.5.OR.IFOUR.EQ.10.OR.IFOUR.EQ.13) THEN 84 : AUX=YOU*ALPDT/(1.-XNU) 85 : RES(1)=AUX 86 : RES(2)=0. 87 : RES(3)=AUX 88 : ELSE 89 : AUX=YOU*ALPDT/(1.-2.*XNU) 90 : RES(1)=AUX 91 : RES(2)=AUX 92 : RES(3)=AUX 93 : ENDIF 94 : C= 1.5 - Cas non prevus 95 : ELSE 96 : KERRE=19 97 : ENDIF 98 : 99 : C 2 - Elements JOINTS 2D PLAN : JOI2 et JOI3 100 : C ============================================ 101 : ELSE IF (MELE.EQ.85.OR.MELE.EQ.86) THEN 102 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) THEN 103 : ALPN=VALMAT(3)*VALMAT(2)*TEMP 104 : RES(1)=0. 105 : RES(2)=ALPN 106 : ELSE 107 : KERRE=19 108 : ENDIF 109 : 110 : C 3 - Element JOINT 3D : JOT4, JOI4 111 : C =================================== 112 : ELSE IF (MELE.EQ.87.OR.MELE.EQ.88) THEN 113 : IF (IFOUR.EQ.2) THEN 114 : ALPN=VALMAT(3)*VALMAT(2)*TEMP 115 : RES(1)=0. 116 : RES(2)=0. 117 : RES(3)=ALPN 118 : ELSE 119 : KERRE=19 120 : ENDIF 121 : 122 : C 4 - Elements COQUES MINCES sauf COQ2 123 : C ====================================== 124 : C+mdj !!!! Ajouter le cas de coque integre (NPINT.NE.0) : 125 : C on recupere un champs par elememt de temperatures aux 126 : C points de GAUSS cette temperature unique remplacera THMO 127 : ELSE IF (MELE.EQ.27.OR.MELE.EQ.28.OR.MELE.EQ.49.OR. 128 : . MELE.EQ.93) THEN 129 : YOU=VALMAT(1) 130 : XNU=VALMAT(2) 131 : ALPH=VALMAT(3) 132 : AUX=YOU*ALPH/(1.D0-XNU) 133 : EXCEN =0. 134 : EP2=EPAIST*EPAIST/12.D0 135 : THMO=(THSU+THIF+4.D0*THM)/6.D0 136 : THDI=THSU-THIF 137 : IF (NPINT.EQ.0) THEN 138 : RES(1)=AUX*THMO*EPAIST 139 : RES(2)=RES(1) 140 : IF (IFOUR.GT.0) THEN 141 : RES(3)=0. 142 : RES(4)=AUX*(THMO*EPAIST*EXCEN+THDI*EP2) 143 : RES(5)=RES(4) 144 : RES(6)=0. 145 : C** ELSE IF (IFOUR.LE.0) THEN 146 : ELSE 147 : RES(3)=AUX*(THMO*EPAIST*EXCEN+THDI*EP2) 148 : RES(4)=RES(3) 149 : ENDIF 150 : ELSE 151 : KERRE=19 152 : ENDIF 153 : 154 : C 5 - Element COQ2 155 : C ================== 156 : ELSE IF (MELE.EQ.44) THEN 157 : YOU=VALMAT(1) 158 : XNU=VALMAT(2) 159 : ALPH=VALMAT(3) 160 : EP2=EPAIST*EPAIST/12.D0 161 : THMO=(THSU+THIF+4.D0*THM)/6.D0 162 : THDI=THSU-THIF 163 : EXCEN =0. 164 : IF (IFOUR.GT.0) THEN 165 : AUX=YOU*ALPH/(1.D0-XNU) 166 : RES(1)=AUX*THMO*EPAIST 167 : RES(2)=RES(1) 168 : RES(3)=0. 169 : RES(4)=AUX*(THMO*EPAIST*EXCEN+THDI*EP2) 170 : RES(5)=RES(4) 171 : RES(6)=0. 172 : ELSE IF (IFOUR.EQ.-2) THEN 173 : AUX=YOU*ALPH 174 : RES(1)=AUX*THMO*EPAIST 175 : RES(2)=0. 176 : RES(3)=AUX*(THMO*EPAIST*EXCEN+THDI*EP2) 177 : RES(4)=0. 178 : ELSE 179 : AUX=YOU*ALPH/(1.D0-XNU) 180 : RES(1)=AUX*THMO*EPAIST 181 : RES(2)=RES(1) 182 : RES(3)=AUX*(THMO*EPAIST*EXCEN+THDI*EP2) 183 : RES(4)=RES(3) 184 : ENDIF 185 : 186 : C 6 - Elements COQUES EPAISSES 187 : C ============================== 188 : ELSE IF (MELE.EQ.41.OR.MELE.EQ.56) THEN 189 : YOU=VALMAT(1) 190 : XNU=VALMAT(2) 191 : ALPH=VALMAT(3) 192 : AUX=YOU*ALPH/(1.-XNU) 193 : THG=0.5*E3*((E3+1.D0)*THSU+(E3-1.D0)*THIF)+(1.-E3*E3)*THM 194 : RES(1)=AUX*THG 195 : RES(2)=RES(1) 196 : RES(3)=0. 197 : RES(4)=0. 198 : RES(5)=0. 199 : 200 : C 7 - Elements POUTRES et TUYAUX 201 : C ================================ 202 : ELSE IF (MELE.EQ.29.OR.MELE.EQ.42.OR.MELE.EQ.84) THEN 203 : YOU=VALMAT(1) 204 : ALPH=VALMAT(2) 205 : RES(1)=YOU*ALPH*SD*TEMP 206 : 207 : C 8 - Elements BARRES, CERCES et POI1 208 : C ===================================== 209 : ELSE IF (MELE.EQ. 46.OR.MELE.EQ. 95.OR.MELE.EQ.45.OR. 210 : . MELE.EQ.123.OR.MELE.EQ.124) THEN 211 : IF (MELE.EQ.45.AND.IFOUR.NE.-3) THEN 212 : KERRE=86 213 : ELSE IF (MELE.EQ.95.AND.IFOUR.NE.0) THEN 214 : KERRE=86 215 : ELSE 216 : YOU=VALMAT(1) 217 : ALPH=VALMAT(2) 218 : RES(1)=YOU*ALPH*SD*TEMP 219 : ENDIF 220 : 221 : C 9 - Autres cas non prevus 222 : C =========================== 223 : ELSE 224 : KERRE=19 225 : ENDIF 226 : 227 : RETURN 228 : END 229 : 230 : 231 : 232 : 233 : 234 :
© Cast3M 2003 - All rights reserved.
Disclaimer