Télécharger la plaquette Cast3M


fr en

Télécharger ttet4f.eso

Retour à la liste

Numérotation des lignes : 

ttet4f

  1. C TTET4F    SOURCE    OF166741  23/12/04    21:15:23     11800          
  2.  
  3. C=======================================================================
  4. C=                            T T E T 4 F                              =
  5. C=                            -----------                              =
  6. C=                                                                     =
  7. C=  Fonction :                                                         =
  8. C=  ----------                                                         =
  9. C=   Acquisition des caracteristiques d'integration aux noeuds, aux    =
  10. C=   points d'integration (Gauss) ou au centre de gravite pour         =
  11. C=   l'element thermique support TET4.                                 =
  12. C=                                                                     =
  13. C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  14. C=  ------------                                                       =
  15. C=   ISUPTH  (E)   Indique le support d'integration demande            =
  16. C=   IPINTE  (S)   Pointeur sur le segment MINTE associe (ACTIF en S)  =
  17. C=======================================================================
  18.  
  19.       SUBROUTINE TTET4F (ISUPTH,IPINTE)
  20.  
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24. -INC SMINTE
  25.  
  26.       PARAMETER ( UN=1.D0 )
  27.       PARAMETER ( XZER=0.D0 )
  28.       PARAMETER ( UNDEMI=.5D0 )
  29.       PARAMETER ( UNQUA=.25D0 )
  30.       PARAMETER ( UNSIX=1.D0/6.D0 )
  31.       PARAMETER ( OHUIT=-.80D0 )
  32.       PARAMETER ( OQUCI= .45D0 )
  33.       PARAMETER ( UNDOU=1.D0/12.D0 )
  34.  
  35.       NBNO = 4
  36.       IF      (ISUPTH.EQ.1) THEN
  37.         NBPGAU = 4
  38.       ELSE IF (ISUPTH.EQ.2) THEN
  39.         NBPGAU = 5
  40.       ELSE IF (ISUPTH.EQ.3) THEN
  41.         NBPGAU = 1
  42.       ELSE
  43.         RETURN
  44.       ENDIF
  45.       SEGINI,MINTE
  46.  
  47. C COORDONNEES ET POIDS DES 4 NOEUDS
  48.       IF (ISUPTH.EQ.1) THEN
  49.         QSIGAU(1) = XZER
  50.         QSIGAU(2) = UN
  51.         QSIGAU(3) = XZER
  52.         QSIGAU(4) = XZER
  53.         ETAGAU(1) = XZER
  54.         ETAGAU(2) = XZER
  55.         ETAGAU(3) = UN
  56.         ETAGAU(4) = XZER
  57.         DZEGAU(1) = XZER
  58.         DZEGAU(2) = XZER
  59.         DZEGAU(3) = XZER
  60.         DZEGAU(4) = UN
  61.         POIGAU(1) = UNDOU
  62.         POIGAU(2) = UNDOU
  63.         POIGAU(3) = UNDOU
  64.         POIGAU(4) = UNDOU
  65.  
  66. C COORDONNEES ET POIDS DES 5 POINTS DE GAUSS
  67.       ELSE IF (ISUPTH.EQ.2) THEN
  68.  
  69.         QSIGAU(1) = UNQUA
  70.         QSIGAU(2) = UNSIX
  71.         QSIGAU(3) = UNSIX
  72.         QSIGAU(4) = UNSIX
  73.         QSIGAU(5) = UNDEMI
  74.  
  75.         ETAGAU(1) = UNQUA
  76.         ETAGAU(2) = UNSIX
  77.         ETAGAU(3) = UNSIX
  78.         ETAGAU(4) = UNDEMI
  79.         ETAGAU(5) = UNSIX
  80.  
  81.         DZEGAU(1) = UNQUA
  82.         DZEGAU(2) = UNSIX
  83.         DZEGAU(3) = UNDEMI
  84.         DZEGAU(4) = UNSIX
  85.         DZEGAU(5) = UNSIX
  86.  
  87.         POIGAU(1) = OHUIT*UNSIX
  88.         POIGAU(2) = OQUCI*UNSIX
  89.         POIGAU(3) = OQUCI*UNSIX
  90.         POIGAU(4) = OQUCI*UNSIX
  91.         POIGAU(5) = OQUCI*UNSIX
  92.  
  93. C COORDONNEES ET POIDS DU CENTRE DE GRAVITE
  94.       ELSE IF (ISUPTH.EQ.3) THEN
  95.         QSIGAU(1) = XZer
  96.         ETAGAU(1) = XZer
  97.         DZEGAU(1) = XZer
  98.         POIGAU(1) = UN/(UN+UN+UN)
  99.  
  100.       ENDIF
  101.  
  102. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  103.       DO I = 1, NBPGAU
  104.         XI = QSIGAU(I)
  105.         YI = ETAGAU(I)
  106.         ZI = DZEGAU(I)
  107.  
  108.         AUX = UN-XI-YI-ZI
  109.         SHPTOT(1,1,I) = AUX
  110.         SHPTOT(1,2,I) = XI
  111.         SHPTOT(1,3,I) = YI
  112.         SHPTOT(1,4,I) = ZI
  113.  
  114.         SHPTOT(2,1,I) = -UN
  115.         SHPTOT(2,2,I) =  UN
  116.         SHPTOT(2,3,I) =  XZER
  117.         SHPTOT(2,4,I) =  XZER
  118.  
  119.         SHPTOT(3,1,I) = -UN
  120.         SHPTOT(3,2,I) =  XZER
  121.         SHPTOT(3,3,I) =  UN
  122.         SHPTOT(3,4,I) =  XZER
  123.  
  124.         SHPTOT(4,1,I) = -UN
  125.         SHPTOT(4,2,I) =  XZER
  126.         SHPTOT(4,3,I) =  XZER
  127.         SHPTOT(4,4,I) =  UN
  128.  
  129.       ENDDO
  130.  
  131.       IPINTE = MINTE
  132.  
  133. c      RETURN
  134.       END
  135.  
  136.  
  137.  

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