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