Télécharger la plaquette Cast3M


fr en

Télécharger tseg2f.eso

Retour à la liste

Numérotation des lignes : 

tseg2f

  1. C TSEG2F    SOURCE    OF166741  23/12/04    21:15:19     11800          
  2.  
  3. C=======================================================================
  4. C=                            T S E G 2 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 SEG2.                                 =
  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 TSEG2F (ISUPTH,IPINTE)
  20.  
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24. -INC SMINTE
  25.  
  26.       PARAMETER ( XZero = 0.D0 , O2=2.D0 , O1=1.D0 , O5=0.5D0 )
  27.       PARAMETER ( O577=0.577350269189625764509148780501957D0 )
  28.  
  29.       NBNO = 2
  30.       IF      (ISUPTH.EQ.1) THEN
  31.         NBPGAU = 2
  32.       ELSE IF (ISUPTH.EQ.2) THEN
  33.         NBPGAU = 2
  34.       ELSE IF (ISUPTH.EQ.3) THEN
  35.         NBPGAU = 1
  36.       ELSE
  37.         CALL ERREUR(5)
  38.         RETURN
  39.       ENDIF
  40.       SEGINI,MINTE
  41.  
  42. C COORDONNEES ET POIDS DES 2 NOEUDS
  43.       IF (ISUPTH.EQ.1) THEN
  44.         QSIGAU(1) = -O1
  45.         QSIGAU(2) = +O1
  46.         POIGAU(1) =  O1
  47.         POIGAU(2) =  O1
  48.  
  49. C COORDONNEES ET POIDS DES 2 POINTS DE GAUSS
  50.       ELSE IF (ISUPTH.EQ.2) THEN
  51.         QSIGAU(1) = -O577
  52.         QSIGAU(2) = +O577
  53.         POIGAU(1) =  O1
  54.         POIGAU(2) =  O1
  55.  
  56. C COORDONNEES ET POIDS DU CENTRE DE GRAVITE
  57.       ELSE IF (ISUPTH.EQ.3) THEN
  58.         QSIGAU(1) = XZero
  59.         POIGAU(1) = O2
  60.       ENDIF
  61.  
  62. * FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  63.       DO I = 1, NBPGAU
  64.         XI = QSIGAU(I)
  65.         SHPTOT(1,1,I) = (O1-XI)*O5
  66.         SHPTOT(1,2,I) = (O1+XI)*O5
  67.  
  68.         SHPTOT(2,1,I) = -O5
  69.         SHPTOT(2,2,I) =  O5
  70.       ENDDO
  71.  
  72.       IPINTE = MINTE
  73.  
  74. c      RETURN
  75.       END
  76.  
  77.  
  78.  

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