Télécharger la plaquette Cast3M


fr en

Télécharger tte10.eso

Retour à la liste

Numérotation des lignes : 

tte10

  1. C TTE10     SOURCE    OF166741  23/12/04    21:15:22     11800          
  2.  
  3. C=======================================================================
  4. C=                            T T E 1 0                                =
  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 TE10.                                 =
  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 TTE10 (ISUPTH,IPINTE)
  20.  
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24. -INC SMINTE
  25.  
  26.       PARAMETER ( XZER=0.D0 , DEUX=2.D0 , UN=1.D0 , QUATRE=4.D0 )
  27.       PARAMETER ( UNDEMI=.5D0 )
  28.       PARAMETER ( UNQUA=.25D0 )
  29.       PARAMETER ( UNSIX=1.D0/6.D0 )
  30.       PARAMETER ( OHUIT=-.80D0 )
  31.       PARAMETER ( OQUCI= .45D0 )
  32.       PARAMETER ( UN30=1.D0/30.D0 )
  33.  
  34.       NBNO = 10
  35.       IF      (ISUPTH.EQ.1) THEN
  36.         NBPGAU = 10
  37.       ELSE IF (ISUPTH.EQ.2) THEN
  38.         NBPGAU = 5
  39.       ELSE IF (ISUPTH.EQ.3) THEN
  40.         NBPGAU = 1
  41.       ELSE
  42.         RETURN
  43.       ENDIF
  44.       SEGINI,MINTE
  45.  
  46. C COORDONNEES ET POIDS DES 10 NOEUDS
  47.       IF (ISUPTH.EQ.1) THEN
  48.         QSIGAU( 2) = UNDEMI
  49.         QSIGAU( 3) = UN
  50.         QSIGAU( 4) = UNDEMI
  51.         QSIGAU( 8) = UNDEMI
  52.         ETAGAU( 4) = UNDEMI
  53.         ETAGAU( 5) = UN
  54.         ETAGAU( 6) = UNDEMI
  55.         ETAGAU( 9) = UNDEMI
  56.         DZEGAU( 7) = UNDEMI
  57.         DZEGAU( 8) = UNDEMI
  58.         DZEGAU( 9) = UNDEMI
  59.         DZEGAU(10) = UN
  60.         DO I = 1, NBPGAU
  61.           POIGAU(I) = UN30
  62.         ENDDO
  63.  
  64. C COORDONNEES ET POIDS DES 5 POINTS DE GAUSS
  65.       ELSE IF (ISUPTH.EQ.2) THEN
  66.  
  67.         QSIGAU(1) = UNQUA
  68.         QSIGAU(2) = UNSIX
  69.         QSIGAU(3) = UNSIX
  70.         QSIGAU(4) = UNSIX
  71.         QSIGAU(5) = UNDEMI
  72.  
  73.         ETAGAU(1) = UNQUA
  74.         ETAGAU(2) = UNSIX
  75.         ETAGAU(3) = UNSIX
  76.         ETAGAU(4) = UNDEMI
  77.         ETAGAU(5) = UNSIX
  78.  
  79.         DZEGAU(1) = UNQUA
  80.         DZEGAU(2) = UNSIX
  81.         DZEGAU(3) = UNDEMI
  82.         DZEGAU(4) = UNSIX
  83.         DZEGAU(5) = UNSIX
  84.  
  85.         POIGAU(1) = OHUIT*UNSIX
  86.         POIGAU(2) = OQUCI*UNSIX
  87.         POIGAU(3) = OQUCI*UNSIX
  88.         POIGAU(4) = OQUCI*UNSIX
  89.         POIGAU(5) = OQUCI*UNSIX
  90.  
  91. C COORDONNEES ET POIDS DU CENTRE DE GRAVITE
  92.       ELSE IF (ISUPTH.EQ.3) THEN
  93.         QSIGAU(1) = XZer
  94.         ETAGAU(1) = XZer
  95.         DZEGAU(1) = XZer
  96.         POIGAU(1) = UN/(DEUX+UN)
  97.       ENDIF
  98.  
  99. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  100.       DO I = 1, NBPGAU
  101.         XI = QSIGAU(I)
  102.         YI = ETAGAU(I)
  103.         ZI = DZEGAU(I)
  104.  
  105.         AUX = UN-XI-YI-ZI
  106.         SHPTOT(1, 1,I) = AUX*(DEUX*AUX-UN)
  107.         SHPTOT(1, 2,I) = QUATRE*XI*AUX
  108.         SHPTOT(1, 3,I) = XI*(DEUX*XI-UN)
  109.         SHPTOT(1, 4,I) = QUATRE*XI*YI
  110.         SHPTOT(1, 5,I) = YI*(DEUX*YI-UN)
  111.         SHPTOT(1, 6,I) = QUATRE*YI*AUX
  112.         SHPTOT(1, 7,I) = QUATRE*ZI*AUX
  113.         SHPTOT(1, 8,I) = QUATRE*XI*ZI
  114.         SHPTOT(1, 9,I) = QUATRE*YI*ZI
  115.         SHPTOT(1,10,I) = ZI*(DEUX*ZI-UN)
  116.  
  117.         SHPTOT(2, 1,I) = UN-QUATRE*AUX
  118.         SHPTOT(2, 2,I) = QUATRE*(AUX-XI)
  119.         SHPTOT(2, 3,I) = QUATRE*XI-UN
  120.         SHPTOT(2, 4,I) = QUATRE*YI
  121.         SHPTOT(2, 5,I) = XZER
  122.         SHPTOT(2, 6,I) = -QUATRE*YI
  123.         SHPTOT(2, 7,I) = -QUATRE*ZI
  124.         SHPTOT(2, 8,I) = QUATRE*ZI
  125.         SHPTOT(2, 9,I) = XZER
  126.         SHPTOT(2,10,I) = XZER
  127.  
  128.         SHPTOT(3, 1,I) = UN-QUATRE*AUX
  129.         SHPTOT(3, 2,I) = -QUATRE*XI
  130.         SHPTOT(3, 3,I) = XZER
  131.         SHPTOT(3, 4,I) = QUATRE*XI
  132.         SHPTOT(3, 5,I) = QUATRE*YI-UN
  133.         SHPTOT(3, 6,I) = QUATRE*(AUX-YI)
  134.         SHPTOT(3, 7,I) = -QUATRE*ZI
  135.         SHPTOT(3, 8,I) = XZER
  136.         SHPTOT(3, 9,I) = QUATRE*ZI
  137.         SHPTOT(3,10,I) = XZER
  138.  
  139.         SHPTOT(4, 1,I) =  UN-QUATRE*AUX
  140.         SHPTOT(4, 2,I) = -QUATRE*XI
  141.         SHPTOT(4, 3,I) =  XZER
  142.         SHPTOT(4, 4,I) =  XZER
  143.         SHPTOT(4, 5,I) =  XZER
  144.         SHPTOT(4, 6,I) = -QUATRE*YI
  145.         SHPTOT(4, 7,I) = QUATRE*(AUX-ZI)
  146.         SHPTOT(4, 8,I) = QUATRE*XI
  147.         SHPTOT(4, 9,I) = QUATRE*YI
  148.         SHPTOT(4,10,I) = QUATRE*ZI-UN
  149.  
  150.       ENDDO
  151.  
  152.       IPINTE = MINTE
  153.  
  154. c      RETURN
  155.       END
  156.  
  157.  
  158.  

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