Télécharger la plaquette Cast3M


fr en

Télécharger tshape.eso

Retour à la liste

Numérotation des lignes : 

  1. C TSHAPE    SOURCE    CHAT      12/06/07    21:16:15     7389
  2.  
  3. C=======================================================================
  4. C=                            T S H A P E                              =
  5. C=                            -----------                              =
  6. C=                                                                     =
  7. C=  Fonction :                                                         =
  8. C=  ----------                                                         =
  9. C=  Calcul des fonctions de forme et de leurs derivees soit aux noeuds =
  10. C=  soit aux points de Gauss d'un element fini MASSIF a integration    =
  11. C=  numerique (branchement vers les sous-programmes adequats).         =
  12. C=                                                                     =
  13. C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  14. C=  ------------                                                       =
  15. C=   NEF     (E)  Numero de l'ELEMENT FINI dans NOMTP (cf. CCHAMP)     =
  16. C=   POINTS  (E)  Chaine de caracteres indiquant si l'on souhaite les  =
  17. C=                valeurs aux noeuds (='NOEUD') ou aux points de Gauss =
  18. C=                (='GAUSS') de l'element fini considere               =
  19. C=   IPINTE  (S)  Pointeur sur le segment MINTE                        =
  20. C=                                                                     =
  21. C=  Remarque : Il s'agit des derivees par rapport aux coordonnees de   =
  22. C=  ---------- l'element de reference (Eta,Qsi,Dzeta).                 =
  23. C=                                                                     =
  24. C=  Denis ROBERT, le 6 janvier 1988. P. DOWLATYARI, mars 1992          =
  25. C=======================================================================
  26.  
  27.       SUBROUTINE TSHAPE (NEF,POINTS,IPINTE)
  28.  
  29.       IMPLICIT INTEGER(I-N)
  30.       IMPLICIT REAL*8 (A-H,O-Z)
  31.  
  32. -INC CCOPTIO
  33.  
  34.       CHARACTER*(*) POINTS
  35.  
  36. C  1 - Verification :  Les valeurs ne peuvent etre donnees qu'aux
  37. C ==================== noeuds ou aux points d'integration (Gauss).
  38.       IF (POINTS(1:5).NE.'GAUSS'.AND.POINTS(1:5).NE.'NOEUD') THEN
  39.         CALL ERREUR(19)
  40.         RETURN
  41.       ENDIF
  42.  
  43. C  2 - Appel a la procedure suivant le type de l'ELEMENT FINI
  44. C ============================================================
  45. C  2.1 - Elements SEG2 RAC2 COQ2 BARRe T1D2
  46. C =====
  47.       IF (NEF.EQ.2.OR.NEF.EQ.12.OR.NEF.EQ.44.OR.NEF.EQ.46.OR.
  48.      .    NEF.EQ.191.or.nef.eq.269) THEN
  49.         CALL TSEG2F(POINTS,IPINTE)
  50. C =====
  51. C  2.2 - Elements SEG3 RAC3 T1D3
  52. C =====
  53.       ELSE IF (NEF.EQ.3.OR.NEF.EQ.13.OR.NEF.EQ.192.or.nef.eq.270) THEN
  54.         CALL TSEG3(POINTS,IPINTE)
  55. C =====
  56. C  2.3 - Elements TRI3 LIA3 COQ3
  57. C =====
  58.       ELSE IF (NEF.EQ.4.OR.NEF.EQ.18.OR.NEF.EQ.27) THEN
  59.         CALL TTRI3F(POINTS,IPINTE)
  60. C =====
  61. C  2.4 - Elements TRI6 LIA6
  62. C =====
  63.       ELSE IF (NEF.EQ.6.OR.NEF.EQ.20) THEN
  64.         CALL TTRI6(POINTS,IPINTE)
  65. C =====
  66. C  2.5 - Elements QUA4 LIA4
  67. C =====
  68.       ELSE IF (NEF.EQ.8.OR.NEF.EQ.19) THEN
  69.         CALL TQUA4(POINTS,IPINTE)
  70. C =====
  71. C  2.6 - Elements QUA8 LIA8
  72. C =====
  73.       ELSE IF (NEF.EQ.10.OR.NEF.EQ.21) THEN
  74.         CALL TQUA8(POINTS,IPINTE)
  75. C =====
  76. C  2.7 - Element CUB8
  77. C =====
  78.       ELSE IF (NEF.EQ.14) THEN
  79.         CALL TCUB8(POINTS,IPINTE)
  80. C =====
  81. C  2.8 - Element CU20
  82. C =====
  83.       ELSE IF (NEF.EQ.15) THEN
  84.         CALL TCU20(POINTS,IPINTE)
  85. C =====
  86. C  2.9 - Element PRI6
  87. C =====
  88.       ELSE IF (NEF.EQ.16) THEN
  89.         CALL TPRI6(POINTS,IPINTE)
  90. C ======
  91. C  2.10 - Element PR15
  92. C ======
  93.       ELSE IF (NEF.EQ.17) THEN
  94.         CALL TPR15(POINTS,IPINTE)
  95. C ======
  96. C  2.11 - Element TET4
  97. C ======
  98.       ELSE IF (NEF.EQ.23) THEN
  99.         CALL TTET4F(POINTS,IPINTE)
  100. C ======
  101. C  2.12 - Element TE10
  102. C ======
  103.       ELSE IF (NEF.EQ.24) THEN
  104.         CALL TTE10(POINTS,IPINTE)
  105. C ======
  106. C  2.13 - Element PYR5
  107. C ======
  108.       ELSE IF (NEF.EQ.25) THEN
  109.         CALL TPYR5(POINTS,IPINTE)
  110. C ======
  111. C  2.14 - Element PY13
  112. C ======
  113.       ELSE IF (NEF.EQ.26) THEN
  114.         CALL TPY13(POINTS,IPINTE)
  115. C ======
  116. C  2.15 - Element COQ4
  117. C ======
  118.       ELSE IF (NEF.EQ.49) THEN
  119.         CALL TCOQ4(POINTS,IPINTE)
  120. C ======
  121. C  2.16 - Element COQ8
  122. C ======
  123.       ELSE IF (NEF.EQ.41) THEN
  124.         CALL TCOQ8(POINTS,IPINTE)
  125. C ======
  126. C  2.17 - Element COQ6
  127. C ======
  128.       ELSE IF (NEF.EQ.56) THEN
  129.         CALL TCOQ6(POINTS,IPINTE)
  130. C ======
  131. C  2.18 - Element POI1
  132. C ======
  133.       ELSE IF ((NEF.EQ.45).OR.(NEF.EQ.1)) THEN
  134.         CALL TPOI1(POINTS,IPINTE)
  135. C ======
  136. C  2.19 - ERREUR : Element fini non implemente
  137. C ======
  138.       ELSE
  139.         CALL ERREUR(19)
  140.       ENDIF
  141.  
  142.       RETURN
  143.       END
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  

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