Télécharger la plaquette Cast3M


fr en

Télécharger ttri3f.eso

Retour à la liste

Numérotation des lignes : 

ttri3f

  1. C TTRI3F    SOURCE    OF166741  23/12/04    21:15:24     11800          
  2.  
  3. =======================================================================
  4. C=                            T T R I 3 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 TRI3.                                 =
  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 TTRI3F (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 )
  27.       PARAMETER ( O1=1.D0 )
  28.       PARAMETER ( O1s2=0.5D0 )
  29.       PARAMETER ( O1s3=0.333333333333333333333333333333333333D0 )
  30.       PARAMETER ( O66=0.666666666666666666666666666666666667D0 )
  31.       PARAMETER ( O16=0.166666666666666666666666666666666667D0 )
  32.  
  33.       NBNO = 3
  34.       IF      (ISUPTH.EQ.1 .OR. ISUPTH.EQ.2) THEN
  35.         NBPGAU = 3
  36.       ELSE IF (ISUPTH.EQ.3) THEN
  37.         NBPGAU = 1
  38.       ELSE
  39.         CALL ERREUR(5)
  40.         RETURN
  41.       ENDIF
  42.       SEGINI,MINTE
  43.  
  44. C COORDONNEES DES 3 NOEUDS
  45.       IF (ISUPTH.EQ.1) THEN
  46.         QSIGAU(1) = O0
  47.         QSIGAU(2) = O1
  48.         QSIGAU(3) = O0
  49.         ETAGAU(1) = O0
  50.         ETAGAU(2) = O0
  51.         ETAGAU(3) = O1
  52.         POIGAU(1) = O16
  53.         POIGAU(2) = O16
  54.         POIGAU(3) = O16
  55.  
  56. C COORDONNEES DES 3 POINTS DE GAUSS
  57.       ELSE IF (ISUPTH.EQ.2) THEN
  58.         QSIGAU(1) = O16
  59.         QSIGAU(2) = O16
  60.         QSIGAU(3) = O66
  61.         ETAGAU(1) = O16
  62.         ETAGAU(2) = O66
  63.         ETAGAU(3) = O16
  64.         POIGAU(1) = O16
  65.         POIGAU(2) = O16
  66.         POIGAU(3) = O16
  67.  
  68. C COORDONNEES DU CENTRE DE GRAVITE
  69.       ELSE IF (ISUPTH.EQ.3) THEN
  70.         QSIGAU(1) = O1s3
  71.         ETAGAU(1) = O1s3
  72.         POIGAU(1) = O1s2
  73.       ENDIF
  74.  
  75. * FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  76.       DO I = 1, NBPGAU
  77.         XI = QSIGAU(I)
  78.         YI = ETAGAU(I)
  79.  
  80.         SHPTOT(1,1,I) = O1-XI-YI
  81.         SHPTOT(1,2,I) = XI
  82.         SHPTOT(1,3,I) = YI
  83.  
  84.         SHPTOT(2,1,I) = -O1
  85.         SHPTOT(2,2,I) =  O1
  86.         SHPTOT(2,3,I) =  O0
  87.  
  88.         SHPTOT(3,1,I) = -O1
  89.         SHPTOT(3,2,I) =  O0
  90.         SHPTOT(3,3,I) =  O1
  91.  
  92.       ENDDO
  93.  
  94.       IPINTE = MINTE
  95.  
  96. C      RETURN
  97.       END
  98.  
  99.  
  100.  

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