Télécharger la plaquette Cast3M


fr en

Télécharger theths.eso

Retour à la liste

Numérotation des lignes : 

theths

  1. C THETHS    SOURCE    BP208322  15/06/22    21:23:23     8543
  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 PPARAM
  22. -INC CCOPTIO
  23. -INC CCHAMP
  24. *
  25.       DIMENSION RES(NSTRS)
  26.       DIMENSION VALMAT(*)
  27. *
  28.       KERRE=0
  29. *
  30. *
  31. *     TEST SUR LA FORMULATION
  32. *
  33. *____________________________________________________________________
  34. *
  35. *     SECTEUR DE CALCUL POUR LES ELEMENTS JOINTS 2D JOI2 ET JOI3
  36. *____________________________________________________________________
  37. *
  38.       IF (MELE.EQ.85.OR.MELE.EQ.86) THEN
  39. *
  40.          ZKN = VALMAT(2)
  41.          ALPN= VALMAT(3)
  42. C
  43.          IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) THEN
  44.             RES(1)=0.D0
  45.             RES(2)=ALPN * ZKN * TEMP
  46.          ELSE
  47.           KERRE=19
  48.           GOTO 99
  49.          ENDIF
  50. *____________________________________________________________________
  51. *
  52. *     SECTEUR DE CALCUL POUR LES ELEMENTS JOINTS 3D JOI4
  53. *____________________________________________________________________
  54. *
  55.       ELSE IF (MELE.EQ.88) THEN
  56.           ZKN = VALMAT(3)
  57.           ALPN= VALMAT(4)
  58.          IF (IFOUR.EQ.2) THEN
  59.             RES(1)=0.D0
  60.             RES(2)=0.D0
  61.             RES(3)=ALPN * ZKN * TEMP
  62.          ELSE
  63.           KERRE=19
  64.           GOTO 99
  65.          ENDIF
  66. *____________________________________________________________________
  67. *
  68. *               ELEMENTS POUTRES ET TUYAUX
  69. *____________________________________________________________________
  70. *
  71.       ELSEIF(MELE.EQ.29.OR.MELE.EQ.42)THEN
  72.             YOU=VALMAT(1)
  73.             ALPH=VALMAT(2)
  74.             AUX = YOU*ALPH*SD
  75. *
  76.             RES(1)= AUX*TEMP
  77. *
  78. *____________________________________________________________________
  79. *
  80. *               ELEMENTS BARRES CERCES ET POINT
  81. *____________________________________________________________________
  82. *
  83.       ELSEIF(MELE.EQ.46.OR.MELE.EQ.95.OR.MELE.EQ.45)THEN
  84. *
  85.        IF(MELE.EQ.45.AND.IFOUR.NE.-3) THEN
  86.           KERRE=86
  87.           GO TO 99
  88.        ENDIF
  89.        IF(MELE.EQ.95.AND.IFOUR.NE.0) THEN
  90.           KERRE=86
  91.           GO TO 99
  92.        ENDIF
  93. *
  94.        YOU=VALMAT(1)
  95.        ALPH=VALMAT(2)
  96.        AUX = YOU*ALPH*SD
  97.        RES(1)= AUX*TEMP
  98. *
  99. *____________________________________________________________________
  100. *
  101. *  CAS DE LA POUTRE TIMO - MODELE SECTION
  102. *____________________________________________________________________
  103. *
  104.       ELSE IF(MELE.EQ.84) THEN
  105.          RES(1)=SD * TEMP
  106.          RES(2)=0.D0
  107.          RES(3)=0.D0
  108.          RES(4)=0.D0
  109.          RES(5)=0.D0
  110.          RES(6)=0.D0
  111.       ELSE
  112. *
  113. *   CAS NON TRAITES
  114. *
  115.        KERRE=19
  116.       ENDIF
  117. *
  118.  99   CONTINUE
  119.       RETURN
  120.       END
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  

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