Télécharger la plaquette Cast3M


fr en

Télécharger tpri6.eso

Retour à la liste

Numérotation des lignes : 

tpri6

  1. C TPRI6     SOURCE    OF166741  23/12/04    21:15:15     11800          
  2.  
  3. C=======================================================================
  4. C=                            T P R I 6                                =
  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 PRI6.                                 =
  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 TPRI6 (ISUPTH,IPINTE)
  20.  
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24. -INC SMINTE
  25.  
  26.       PARAMETER ( O0=0.D0, O1=1.D0, O5=0.5D0 )
  27.       PARAMETER ( W=1.D0/6.D0 )
  28.       PARAMETER ( Z=0.577350269189625764509148780501957D0 )
  29.  
  30.       NBNO = 6
  31.       IF      (ISUPTH.EQ.1) THEN
  32.         NBPGAU = 6
  33.       ELSE IF (ISUPTH.EQ.2) THEN
  34.         NBPGAU = 6
  35.       ELSE IF (ISUPTH.EQ.3) THEN
  36.         NBPGAU = 1
  37.       ELSE
  38.         CALL ERREUR(5)
  39.         RETURN
  40.       ENDIF
  41.       SEGINI,MINTE
  42.  
  43. C COORDONNEES INTRINSEQUES DES 6 NOEUDS
  44.       IF      (ISUPTH.EQ.1) THEN
  45.         QSIGAU(1) = O0
  46.         QSIGAU(2) = O1
  47.         QSIGAU(3) = O0
  48.         QSIGAU(4) = O0
  49.         QSIGAU(5) = O1
  50.         QSIGAU(6) = O0
  51.         ETAGAU(1) = O0
  52.         ETAGAU(2) = O0
  53.         ETAGAU(3) = O1
  54.         ETAGAU(4) = O0
  55.         ETAGAU(5) = O0
  56.         ETAGAU(6) = O1
  57.         DZEGAU(1) = -O1
  58.         DZEGAU(2) = -O1
  59.         DZEGAU(3) = -O1
  60.         DZEGAU(4) =  O1
  61.         DZEGAU(5) =  O1
  62.         DZEGAU(6) =  O1
  63. C POIDS
  64.         DO I = 1, NBPGAU
  65.           POIGAU(I) = W
  66.         ENDDO
  67.  
  68. C COORDONNEES ET POIDS DES 6 POINTS DE GAUSS
  69.       ELSE IF (ISUPTH.EQ.2) THEN
  70.         QSIGAU(1) = O5
  71.         QSIGAU(2) = O5
  72.         QSIGAU(3) = O0
  73.         QSIGAU(4) = O5
  74.         QSIGAU(5) = O5
  75.         QSIGAU(6) = O0
  76.         ETAGAU(1) = O0
  77.         ETAGAU(2) = O5
  78.         ETAGAU(3) = O5
  79.         ETAGAU(4) = O0
  80.         ETAGAU(5) = O5
  81.         ETAGAU(6) = O5
  82.         DZEGAU(1) = -Z
  83.         DZEGAU(2) = -Z
  84.         DZEGAU(3) = -Z
  85.         DZEGAU(4) =  Z
  86.         DZEGAU(5) =  Z
  87.         DZEGAU(6) =  Z
  88. C POIDS
  89.         DO I = 1, NBPGAU
  90.           POIGAU(I) = W
  91.         ENDDO
  92.  
  93. C COORDONNEES DU CENTRE DE GRAVITE
  94.       ELSE IF (ISUPTH.EQ.3) THEN
  95.         QSIGAU(1) = O0
  96.         ETAGAU(1) = O0
  97.         DZEGAU(1) = O0
  98.         POIGAU(1) = O1
  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. C- FONCTIONS DE FORME DU PRISME
  108.         SHPTOT(1,1,I) = O5*(O1-XI-YI)*(O1-ZI)
  109.         SHPTOT(1,2,I) = O5*    XI    *(O1-ZI)
  110.         SHPTOT(1,3,I) = O5*    YI    *(O1-ZI)
  111.         SHPTOT(1,4,I) = O5*(O1-XI-YI)*(O1+ZI)
  112.         SHPTOT(1,5,I) = O5*    XI    *(O1+ZI)
  113.         SHPTOT(1,6,I) = O5*    YI    *(O1+ZI)
  114. C- DERIVEES EN XI DES FONCTIONS DE FORME
  115.         SHPTOT(2,1,I) = -O5*(O1-ZI)
  116.         SHPTOT(2,2,I) =  O5*(O1-ZI)
  117.         SHPTOT(2,3,I) =  O0
  118.         SHPTOT(2,4,I) = -O5*(O1+ZI)
  119.         SHPTOT(2,5,I) =  O5*(O1+ZI)
  120.         SHPTOT(2,6,I) =  O0
  121. C- DERIVEES EN YI DES FONCTIONS DE FORME
  122.         SHPTOT(3,1,I)=-O5*(O1-ZI)
  123.         SHPTOT(3,2,I)= O0
  124.         SHPTOT(3,3,I)= O5*(O1-ZI)
  125.         SHPTOT(3,4,I)=-O5*(O1+ZI)
  126.         SHPTOT(3,5,I)= O0
  127.         SHPTOT(3,6,I)= O5*(O1+ZI)
  128. C- DERIVEES EN ZI DES FONCTIONS DE FORME
  129.         SHPTOT(4,1,I) = -O5*(O1-XI-YI)
  130.         SHPTOT(4,2,I) = -O5*XI
  131.         SHPTOT(4,3,I) = -O5*YI
  132.         SHPTOT(4,4,I) = O5*(O1-XI-YI)
  133.         SHPTOT(4,5,I) = O5*XI
  134.         SHPTOT(4,6,I) = O5*YI
  135.       ENDDO
  136.  
  137.       IPINTE = MINTE
  138.  
  139. c      RETURN
  140.       END
  141.  
  142.  
  143.  

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