Télécharger la plaquette Cast3M


fr en

Télécharger tqua8.eso

Retour à la liste

Numérotation des lignes : 

tqua8

  1. C TQUA8     SOURCE    OF166741  24/12/09    21:15:02     12093          
  2.  
  3. C=======================================================================
  4. C=                            T Q U A 8                                =
  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 massif QUA8.                                  =
  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 TQUA8 (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, O25 = 0.25D0 )
  27.       PARAMETER ( O4 = 4.D0 )
  28.       PARAMETER ( XI  = 0.774596669241483377035853079956479D0 )
  29.       PARAMETER ( XI2 = 0.555555555555555555555555555555555D0 )
  30.       PARAMETER ( ETE = 0.888888888888888888888888888888888D0 )
  31.       PARAMETER ( XI3 = XI2*XI2 )
  32.       PARAMETER ( ET2 = ETE*ETE )
  33.       PARAMETER ( XIET = XI2*ETE )
  34.  
  35.       NBNO = 8                                                                 
  36.       IF      (ISUPTH.EQ.1) THEN
  37.         NBPGAU = 8
  38.       ELSE IF (ISUPTH.EQ.2) THEN
  39.         NBPGAU = 9
  40.       ELSE IF (ISUPTH.EQ.3) THEN
  41.         NBPGAU = 1
  42.       ELSE
  43.         CALL ERREUR(5)
  44.         RETURN
  45.       ENDIF
  46.       SEGINI,MINTE
  47.  
  48. C POIDS ET COORDONNEES DES 8 NOEUDS
  49.       IF (ISUPTH.EQ.1) THEN
  50.         QSIGAU(1) =  O1
  51.         QSIGAU(2) =  O0
  52.         QSIGAU(3) = -O1
  53.         QSIGAU(4) = -O1
  54.         QSIGAU(5) = -O1
  55.         QSIGAU(6) =  O0
  56.         QSIGAU(7) =  O1
  57.         QSIGAU(8) =  O1
  58.  
  59.         ETAGAU(1) =  O1
  60.         ETAGAU(2) =  O1
  61.         ETAGAU(3) =  O1
  62.         ETAGAU(4) =  O0
  63.         ETAGAU(5) = -O1
  64.         ETAGAU(6) = -O1
  65.         ETAGAU(7) = -O1
  66.         ETAGAU(8) =  O0
  67.  
  68.         POIGAU(1) = O5
  69.         POIGAU(2) = O5
  70.         POIGAU(3) = O5
  71.         POIGAU(4) = O5
  72.         POIGAU(5) = O5
  73.         POIGAU(6) = O5
  74.         POIGAU(7) = O5
  75.         POIGAU(8) = O5
  76.  
  77. C POIDS ET COORDONNEES DES 9 POINTS DE GAUSS
  78.       ELSE IF (ISUPTH.EQ.2) THEN
  79.         QSIGAU(1) = -XI
  80.         QSIGAU(2) =  O0
  81.         QSIGAU(3) =  XI
  82.         QSIGAU(4) = -XI
  83.         QSIGAU(5) =  O0
  84.         QSIGAU(6) =  XI
  85.         QSIGAU(7) = -XI
  86.         QSIGAU(8) =  O0
  87.         QSIGAU(9) =  XI
  88.  
  89.         ETAGAU(1) = -XI
  90.         ETAGAU(2) = -XI
  91.         ETAGAU(3) = -XI
  92.         ETAGAU(4) =  O0
  93.         ETAGAU(5) =  O0
  94.         ETAGAU(6) =  O0
  95.         ETAGAU(7) =  XI
  96.         ETAGAU(8) =  XI
  97.         ETAGAU(9) =  XI
  98.  
  99.         POIGAU(1) = XI3
  100.         POIGAU(2) = XIET
  101.         POIGAU(3) = XI3
  102.         POIGAU(4) = XIET
  103.         POIGAU(5) = ET2
  104.         POIGAU(6) = XIET
  105.         POIGAU(7) = XI3
  106.         POIGAU(8) = XIET
  107.         POIGAU(9) = XI3
  108.  
  109. C POIDS ET COORDONNEES DU CENTRE DE GRAVITE
  110.       ELSE IF (ISUPTH.EQ.3) THEN
  111.         QSIGAU(1) = O0
  112.         ETAGAU(1) = O0
  113.         POIGAU(1) = O4
  114.       ENDIF
  115.  
  116. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  117.       DO I = 1, NBPGAU
  118.         X = QSIGAU(I)
  119.         Y = ETAGAU(I)
  120.  
  121.         SHPTOT(1,1,I) = -O25*(O1+X+Y)*(O1-X)*(O1-Y)
  122.         SHPTOT(1,2,I) = O5*(O1-X*X)*(O1-Y)
  123.         SHPTOT(1,3,I) = -O25*(O1-X+Y)*(O1+X)*(O1-Y)
  124.         SHPTOT(1,4,I) = O5*(O1+X)*(O1-Y*Y)
  125.         SHPTOT(1,5,I) = -O25*(O1-X-Y)*(O1+X)*(O1+Y)
  126.         SHPTOT(1,6,I) = O5*(O1-X*X)*(O1+Y)
  127.         SHPTOT(1,7,I) = -O25*(O1+X-Y)*(O1-X)*(O1+Y)
  128.         SHPTOT(1,8,I) = O5*(O1-X)*(O1-Y*Y)
  129.  
  130.         SHPTOT(2,1,I) = O25*(O1-Y)*(X+X+Y)
  131.         SHPTOT(2,2,I) = -X*(O1-Y)
  132.         SHPTOT(2,3,I) = O25*(O1-Y)*(X+X-Y)
  133.         SHPTOT(2,4,I) = O5*(O1-Y*Y)
  134.         SHPTOT(2,5,I) = O25*(O1+Y)*(X+X+Y)
  135.         SHPTOT(2,6,I) = -X*(O1+Y)
  136.         SHPTOT(2,7,I) = O25*(O1+Y)*(X+X-Y)
  137.         SHPTOT(2,8,I) = -O5*(O1-Y*Y)
  138.  
  139.         SHPTOT(3,1,I) = O25*(O1-X)*(X+Y+Y)
  140.         SHPTOT(3,2,I) = -O5*(O1-X*X)
  141.         SHPTOT(3,3,I) = O25*(O1+X)*(-X+Y+Y)
  142.         SHPTOT(3,4,I) = -Y*(O1+X)
  143.         SHPTOT(3,5,I) = O25*(O1+X)*(X+Y+Y)
  144.         SHPTOT(3,6,I) = O5*(O1-X*X)
  145.         SHPTOT(3,7,I) = O25*(O1-X)*(Y+Y-X)
  146.         SHPTOT(3,8,I) = -Y*(O1-X)
  147.       ENDDO
  148.  
  149.       IPINTE = MINTE
  150.  
  151. c      RETURN
  152.       END
  153.  
  154.  
  155.  

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