Télécharger la plaquette Cast3M


fr en

Télécharger tcoq4.eso

Retour à la liste

Numérotation des lignes : 

tcoq4

  1. C TCOQ4     SOURCE    OF166741  23/12/04    21:15:08     11800          
  2.  
  3. C=======================================================================
  4. C=                            T C O Q 4                                =
  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 COQ4.                           =
  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 TCOQ4 (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)
  27.       PARAMETER ( O1=1.D0 ,O4=4.D0, O25=0.25D0 )
  28.       PARAMETER ( X=0.577350269189625764509148780501957D0 )
  29.  
  30.       NBNO = 4
  31.       IF      (ISUPTH.EQ.1) THEN
  32.         NBPGAU = 4
  33.       ELSE IF (ISUPTH.EQ.2) THEN
  34.         NBPGAU = 8
  35.       ELSE IF (ISUPTH.EQ.3) THEN
  36.         NBPGAU = 1
  37.       ELSE
  38.         CALL ERREUR(5)
  39.         RETURN
  40.       ENDIF
  41.       SEGINI,MINTE
  42.  
  43. C COORDONNEES DES 4 NOEUDS
  44.       IF (ISUPTH.EQ.1) THEN
  45.         QSIGAU(1) = -O1
  46.         QSIGAU(2) =  O1
  47.         QSIGAU(3) =  O1
  48.         QSIGAU(4) = -O1
  49.         ETAGAU(1) = -O1
  50.         ETAGAU(2) = -O1
  51.         ETAGAU(3) =  O1
  52.         ETAGAU(4) =  O1
  53.         POIGAU(1) = O1
  54.         POIGAU(2) = O1
  55.         POIGAU(3) = O1
  56.         POIGAU(4) = O1
  57.  
  58. C COORDONNEES DES 8 POINTS DE GAUSS
  59. C 2 NAPPES DE 4 POINTS DE GAUSS
  60.       ELSE IF (ISUPTH.EQ.2) THEN
  61.         QSIGAU(1) = -X
  62.         QSIGAU(2) =  X
  63.         QSIGAU(3) =  X
  64.         QSIGAU(4) = -X
  65.         QSIGAU(5) = -X
  66.         QSIGAU(6) =  X
  67.         QSIGAU(7) =  X
  68.         QSIGAU(8) = -X
  69.         ETAGAU(1) = -X
  70.         ETAGAU(2) = -X
  71.         ETAGAU(3) =  X
  72.         ETAGAU(4) =  X
  73.         ETAGAU(5) = -X
  74.         ETAGAU(6) = -X
  75.         ETAGAU(7) =  X
  76.         ETAGAU(8) =  X
  77.         DZEGAU(1) = -X
  78.         DZEGAU(2) = -X
  79.         DZEGAU(3) = -X
  80.         DZEGAU(4) = -X
  81.         DZEGAU(5) =  X
  82.         DZEGAU(6) =  X
  83.         DZEGAU(7) =  X
  84.         DZEGAU(8) =  X
  85.         DO i = 1, NBPGAU
  86.           POIGAU(i) = O1
  87.         ENDDO
  88.  
  89. C COORDONNEES DU CENTRE DE GRAVITE
  90.       ELSE IF (ISUPTH.EQ.3) THEN
  91.         QSIGAU(1) = XZero
  92.         ETAGAU(1) = XZero
  93.         POIGAU(1) = O4
  94.       ENDIF
  95.  
  96. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  97.       DO I = 1, NBPGAU
  98.         XI = QSIGAU(I)
  99.         YI = ETAGAU(I)
  100.  
  101.         SHPTOT(1,1,I) =  O25*(O1-XI)*(O1-YI)
  102.         SHPTOT(1,2,I) =  O25*(O1+XI)*(O1-YI)
  103.         SHPTOT(1,3,I) =  O25*(O1+XI)*(O1+YI)
  104.         SHPTOT(1,4,I) =  O25*(O1-XI)*(O1+YI)
  105. C
  106.         SHPTOT(2,1,I) = -O25*(O1-YI)
  107.         SHPTOT(2,2,I) =  O25*(O1-YI)
  108.         SHPTOT(2,3,I) =  O25*(O1+YI)
  109.         SHPTOT(2,4,I) = -O25*(O1+YI)
  110. C
  111.         SHPTOT(3,1,I) = -O25*(O1-XI)
  112.         SHPTOT(3,2,I) = -O25*(O1+XI)
  113.         SHPTOT(3,3,I) =  O25*(O1+XI)
  114.         SHPTOT(3,4,I) =  O25*(O1-XI)
  115.       ENDDO
  116.  
  117.       IPINTE = MINTE
  118.  
  119. c      RETURN
  120.       END
  121.  
  122.  
  123.  

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