1 : C THETHS SOURCE AM 08/02/14 21:31:46 6050 2 : SUBROUTINE THETHS(MELE,VALMAT,NSTRS,TEMP,SD,RES,KERRE) 3 : ******************************************************************* 4 : * APPELE PAR THETAP :CONTRAINTES THERMIQUES POUR LES MATERIAUX 5 : * HOMOGENEISES ET SECTION 6 : * 7 : * ENTREES : 8 : * MELE = NUMERO DE L'ELEMENT FINI 9 : * VALMAT = TABLEAU DE MATERIAU 10 : * TEMP = TEMPERATURES 11 : * SD = SECTION 12 : * SORTIES : 13 : * RES = TABLEAU DE CONTRAINTES 14 : * KERRE= INDICATEUR D'ERREUR 15 : * 16 : * P DOWLATYARI MARS 92 17 : *********************************************************************** 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : * 21 : -INC CCOPTIO 22 : -INC CCHAMP 23 : * 24 : DIMENSION RES(NSTRS) 25 : DIMENSION VALMAT(*) 26 : * 27 : KERRE=0 28 : * 29 : * 30 : * TEST SUR LA FORMULATION 31 : * 32 : *____________________________________________________________________ 33 : * 34 : * SECTEUR DE CALCUL POUR LES ELEMENTS JOINTS 2D JOI2 ET JOI3 35 : *____________________________________________________________________ 36 : * 37 : IF (MELE.EQ.85.OR.MELE.EQ.86) THEN 38 : * 39 : ZKN = VALMAT(2) 40 : ALPN= VALMAT(3) 41 : C 42 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) THEN 43 : RES(1)=0.D0 44 : RES(2)=ALPN * ZKN * TEMP 45 : ELSE 46 : KERRE=19 47 : GOTO 99 48 : ENDIF 49 : *____________________________________________________________________ 50 : * 51 : * SECTEUR DE CALCUL POUR LES ELEMENTS JOINTS 3D JOI4 52 : *____________________________________________________________________ 53 : * 54 : ELSE IF (MELE.EQ.88) THEN 55 : ZKN = VALMAT(3) 56 : ALPN= VALMAT(4) 57 : IF (IFOUR.EQ.2) THEN 58 : RES(1)=0.D0 59 : RES(2)=0.D0 60 : RES(3)=ALPN * ZKN * TEMP 61 : ELSE 62 : KERRE=19 63 : GOTO 99 64 : ENDIF 65 : *____________________________________________________________________ 66 : * 67 : * ELEMENTS POUTRES ET TUYAUX 68 : *____________________________________________________________________ 69 : * 70 : ELSEIF(MELE.EQ.29.OR.MELE.EQ.42)THEN 71 : YOU=VALMAT(1) 72 : ALPH=VALMAT(2) 73 : AUX = YOU*ALPH*SD 74 : * 75 : RES(1)= AUX*TEMP 76 : * 77 : *____________________________________________________________________ 78 : * 79 : * ELEMENTS BARRES CERCES ET POINT 80 : *____________________________________________________________________ 81 : * 82 : ELSEIF(MELE.EQ.46.OR.MELE.EQ.95.OR.MELE.EQ.45)THEN 83 : * 84 : IF(MELE.EQ.45.AND.IFOUR.NE.-3) THEN 85 : KERRE=86 86 : GO TO 99 87 : ENDIF 88 : IF(MELE.EQ.95.AND.IFOUR.NE.0) THEN 89 : KERRE=86 90 : GO TO 99 91 : ENDIF 92 : * 93 : YOU=VALMAT(1) 94 : ALPH=VALMAT(2) 95 : AUX = YOU*ALPH*SD 96 : RES(1)= AUX*TEMP 97 : * 98 : *____________________________________________________________________ 99 : * 100 : * CAS DE LA POUTRE TIMO - MODELE SECTION 101 : *____________________________________________________________________ 102 : * 103 : ELSE IF(MELE.EQ.84) THEN 104 : RES(1)=SD * TEMP 105 : RES(2)=0.D0 106 : RES(3)=0.D0 107 : RES(4)=0.D0 108 : RES(5)=0.D0 109 : RES(6)=0.D0 110 : ELSE 111 : * 112 : * CAS NON TRAITES 113 : * 114 : KERRE=19 115 : ENDIF 116 : * 117 : 99 CONTINUE 118 : RETURN 119 : END 120 : 121 : 122 : 123 :
© Cast3M 2003 - All rights reserved.
Disclaimer