Télécharger la plaquette Cast3M


fr en

Télécharger tseg3.eso

Retour à la liste

Numérotation des lignes : 

tseg3

  1. C TSEG3     SOURCE    OF166741  23/12/04    21:15:20     11800          
  2.  
  3. C=======================================================================
  4. C=                            T S E G 3                                =
  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 SEG3.                                 =
  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 TSEG3 (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 , O2=2.D0 , O1=1.D0 , O1s2=0.5D0 )
  27.       PARAMETER ( O2s3=0.666666666666666666666666666666667D0 )
  28.       PARAMETER ( O555=0.555555555555555555555555555555556D0 )
  29.       PARAMETER ( O774=0.774596669241483377035853079956479D0 )
  30.       PARAMETER ( O888=0.888888888888888888888888888888889D0 )
  31.  
  32.       NBNO = 3
  33.       IF      (ISUPTH.EQ.1) THEN
  34.         NBPGAU = 3
  35.       ELSE IF (ISUPTH.EQ.2) THEN
  36.         NBPGAU = 3
  37.       ELSE IF (ISUPTH.EQ.3) THEN
  38.         NBPGAU = 1
  39.       ELSE
  40.         CALL ERREUR(5)
  41.         RETURN
  42.       ENDIF
  43.       SEGINI,MINTE
  44.  
  45. * COORDONNEES ET POIDS DES 3 NOEUDS
  46.       IF (ISUPTH.EQ.1) THEN
  47.         QSIGAU(1) = -O1
  48.         QSIGAU(2) =  O0
  49.         QSIGAU(3) = +O1
  50.         POIGAU(1) =  O2s3
  51.         POIGAU(2) =  O2s3
  52.         POIGAU(3) =  O2s3
  53.  
  54. * COORDONNEES ET POIDS DES 3 POINTS DE GAUSS
  55.       ELSE IF (ISUPTH.EQ.2) THEN
  56.         QSIGAU(1) = -O774
  57.         QSIGAU(2) =  O0
  58.         QSIGAU(3) = +O774
  59.         POIGAU(1) =  O555
  60.         POIGAU(2) =  O888
  61.         POIGAU(3) =  O555
  62.  
  63. * COORDONNEES ET POIDS DU CENTRE DE GRAVITE
  64.       ELSE IF (ISUPTH.EQ.3) THEN
  65.         QSIGAU(1) = O0
  66.         POIGAU(1) = O2
  67.       ENDIF
  68.  
  69. * FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  70.       DO I = 1, NBPGAU
  71.         XI = QSIGAU(I)
  72.  
  73.         SHPTOT(1,1,I) = O1s2*XI*(XI-O1)
  74.         SHPTOT(1,2,I) = O1-XI*XI
  75.         SHPTOT(1,3,I) = O1s2*XI*(XI+O1)
  76.  
  77.         SHPTOT(2,1,I) = XI-O1s2
  78.         SHPTOT(2,2,I) = -O2*XI
  79.         SHPTOT(2,3,I) = XI+O1s2
  80.       ENDDO
  81.  
  82.       IPINTE = MINTE
  83.  
  84. C      RETURN
  85.       END
  86.  
  87.  
  88.  

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