Télécharger la plaquette Cast3M


fr en

Télécharger tcoq6.eso

Retour à la liste

Numérotation des lignes : 

tcoq6

  1. C TCOQ6     SOURCE    OF166741  23/12/04    21:15:09     11800          
  2.  
  3. C=======================================================================
  4. C=                            T C O Q 6                                =
  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 coque COQ6.                           =
  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 TCOQ6 (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, O2=2.D0, O3=3.D0, O4=4.D0, O5=0.5D0, O6=6.D0,
  28.      &            O9=9.D0, O12=12.D0, O21=21.D0, O80=80.D0, O155=155.D0,
  29.      &            O2400=2400.D0 )
  30.       PARAMETER ( X=0.577350269189625764509148780501957D0 )
  31.  
  32.       NBNO = 6
  33.       IF      (ISUPTH.EQ.1) THEN
  34.         NBPGAU = 6
  35.       ELSE IF (ISUPTH.EQ.2) THEN
  36.         NBPGAU = 14
  37.       ELSE IF (ISUPTH.EQ.3) THEN
  38.         NBPGAU = 1
  39.       ELSE
  40.         RETURN
  41.       ENDIF
  42.       SEGINI,MINTE
  43.  
  44. C COORDONNEES DES 6 NOEUDS
  45.       IF (ISUPTH.EQ.1) THEN
  46.         QSIGAU(1) = O0
  47.         QSIGAU(2) = O5
  48.         QSIGAU(3) = O1
  49.         QSIGAU(4) = O5
  50.         QSIGAU(5) = O0
  51.         QSIGAU(6) = O0
  52.         ETAGAU(1) = O0
  53.         ETAGAU(2) = O0
  54.         ETAGAU(3) = O0
  55.         ETAGAU(4) = O5
  56.         ETAGAU(5) = O1
  57.         ETAGAU(6) = O5
  58.         POIGAU(1) = O1/O12
  59.         POIGAU(2) = O1/O12
  60.         POIGAU(3) = O1/O12
  61.         POIGAU(4) = O1/O12
  62.         POIGAU(5) = O1/O12
  63.         POIGAU(6) = O1/O12
  64.  
  65. C COORDONNEES DES 14 POINTS DE GAUSS
  66. C 2 NAPPES DE 7 POINTS DE GAUSS
  67.       ELSE IF (ISUPTH.EQ.2) THEN
  68.         rac15 = SQRT(15.D0)
  69.         QSIGAU(1) = O1/O3
  70.         QSIGAU(2) = (O9-O2*rac15)/O21
  71.         QSIGAU(3) = (O6+rac15)/O21
  72.         QSIGAU(4) = (O6+rac15)/O21
  73.         QSIGAU(5) = (O9+O2*rac15)/O21
  74.         QSIGAU(6) = (O6-rac15)/O21
  75.         QSIGAU(7) = (O6-rac15)/O21
  76.         ETAGAU(1) = O1/O3
  77.         ETAGAU(2) = (O6+rac15)/O21
  78.         ETAGAU(3) = (O9-O2*rac15)/O21
  79.         ETAGAU(4) = (O6+rac15)/O21
  80.         ETAGAU(5) = (O6-rac15)/O21
  81.         ETAGAU(6) = (O9+O2*rac15)/O21
  82.         ETAGAU(7) = (O6-rac15)/O21
  83.         POIGAU(1) = O9/O80
  84.         POIGAU(2) = (O155+rac15)/O2400
  85.         POIGAU(3) = (O155+rac15)/O2400
  86.         POIGAU(4) = (O155+rac15)/O2400
  87.         POIGAU(5) = (O155-rac15)/O2400
  88.         POIGAU(6) = (O155-rac15)/O2400
  89.         POIGAU(7) = (O155-rac15)/O2400
  90.         DO I = 1, 7
  91.           QSIGAU(I+7) = QSIGAU(I)
  92.           ETAGAU(I+7) = ETAGAU(I)
  93.           DZEGAU(I)   = -X
  94.           DZEGAU(I+7) =  X
  95.           POIGAU(I+7) = POIGAU(I)
  96.         ENDDO
  97.  
  98. C COORDONNEES DU CENTRE DE GRAVITE
  99.       ELSE IF (ISUPTH.EQ.3) THEN
  100.         QSIGAU(1) = O1/O3
  101.         ETAGAU(1) = O1/O3
  102.         POIGAU(1) = O5
  103.       ENDIF
  104.  
  105. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  106.       DO I = 1, NBPGAU
  107.         XI = QSIGAU(I)
  108.         YI = ETAGAU(I)
  109.  
  110.         SHPTOT(1,1,I) = (O1-XI-YI)*(O1-O2*(XI+YI))
  111.         SHPTOT(1,2,I) = O4*XI*(O1-XI-YI)
  112.         SHPTOT(1,3,I) = XI*(XI+XI-O1)
  113.         SHPTOT(1,4,I) = O4*XI*YI
  114.         SHPTOT(1,5,I) = YI*(YI+YI-O1)
  115.         SHPTOT(1,6,I) = O4*YI*(O1-XI-YI)
  116.  
  117.         SHPTOT(2,1,I) = O4*(XI+YI)-O3
  118.         SHPTOT(2,2,I) = O4*(O1-XI-XI-YI)
  119.         SHPTOT(2,3,I) = O4*XI-O1
  120.         SHPTOT(2,4,I) = O4*YI
  121.         SHPTOT(2,5,I) = O0
  122.         SHPTOT(2,6,I) = -O4*YI
  123.  
  124.         SHPTOT(3,1,I) = O4*(XI+YI)-O3
  125.         SHPTOT(3,2,I) = -O4*XI
  126.         SHPTOT(3,3,I) = O0
  127.         SHPTOT(3,4,I) = O4*XI
  128.         SHPTOT(3,5,I) = O4*YI-O1
  129.         SHPTOT(3,6,I) = O4*(O1-XI-YI-YI)
  130.  
  131.       ENDDO
  132.  
  133.       IPINTE = MINTE
  134.  
  135. c      RETURN
  136.       END
  137.  
  138.  
  139.  

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