Esope source files

   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
  38 :         CALL ERREUR(19)
  39 :         RETURN
  40 :       ENDIF
  41 : 
  42 : C  2 - Appel a la procedure suivant le type de l'ELEMENT FINI
  43 : C ============================================================
  44 : C  2.1 - Elements SEG2 RAC2 COQ2 BARRe T1D2
  45 : C =====
  46 :       IF (NEF.EQ.2.OR.NEF.EQ.12.OR.NEF.EQ.44.OR.NEF.EQ.46.OR.
  47 :      .    NEF.EQ.191.or.nef.eq.269) THEN
  47 :         CALL TSEG2F(POINTS,IPINTE)
  48 : C =====
  49 : C  2.2 - Elements SEG3 RAC3 T1D3
  50 : C =====
  51 :       ELSE IF (NEF.EQ.3.OR.NEF.EQ.13.OR.NEF.EQ.192.or.nef.eq.270) THEN
  51 :         CALL TSEG3(POINTS,IPINTE)
  52 : C =====
  53 : C  2.3 - Elements TRI3 LIA3 COQ3
  54 : C =====
  55 :       ELSE IF (NEF.EQ.4.OR.NEF.EQ.18.OR.NEF.EQ.27) THEN
  55 :         CALL TTRI3F(POINTS,IPINTE)
  56 : C =====
  57 : C  2.4 - Elements TRI6 LIA6
  58 : C =====
  59 :       ELSE IF (NEF.EQ.6.OR.NEF.EQ.20) THEN
  59 :         CALL TTRI6(POINTS,IPINTE)
  60 : C =====
  61 : C  2.5 - Elements QUA4 LIA4
  62 : C =====
  63 :       ELSE IF (NEF.EQ.8.OR.NEF.EQ.19) THEN
  63 :         CALL TQUA4(POINTS,IPINTE)
  64 : C =====
  65 : C  2.6 - Elements QUA8 LIA8
  66 : C =====
  67 :       ELSE IF (NEF.EQ.10.OR.NEF.EQ.21) THEN
  67 :         CALL TQUA8(POINTS,IPINTE)
  68 : C =====
  69 : C  2.7 - Element CUB8
  70 : C =====
  71 :       ELSE IF (NEF.EQ.14) THEN
  71 :         CALL TCUB8(POINTS,IPINTE)
  72 : C =====
  73 : C  2.8 - Element CU20
  74 : C =====
  75 :       ELSE IF (NEF.EQ.15) THEN
  75 :         CALL TCU20(POINTS,IPINTE)
  76 : C =====
  77 : C  2.9 - Element PRI6
  78 : C =====
  79 :       ELSE IF (NEF.EQ.16) THEN
  79 :         CALL TPRI6(POINTS,IPINTE)
  80 : C ======
  81 : C  2.10 - Element PR15
  82 : C ======
  83 :       ELSE IF (NEF.EQ.17) THEN
  83 :         CALL TPR15(POINTS,IPINTE)
  84 : C ======
  85 : C  2.11 - Element TET4
  86 : C ======
  87 :       ELSE IF (NEF.EQ.23) THEN
  87 :         CALL TTET4F(POINTS,IPINTE)
  88 : C ======
  89 : C  2.12 - Element TE10
  90 : C ======
  91 :       ELSE IF (NEF.EQ.24) THEN
  91 :         CALL TTE10(POINTS,IPINTE)
  92 : C ======
  93 : C  2.13 - Element PYR5
  94 : C ======
  95 :       ELSE IF (NEF.EQ.25) THEN
  95 :         CALL TPYR5(POINTS,IPINTE)
  96 : C ======
  97 : C  2.14 - Element PY13
  98 : C ======
  99 :       ELSE IF (NEF.EQ.26) THEN
  99 :         CALL TPY13(POINTS,IPINTE)
 100 : C ======
 101 : C  2.15 - Element COQ4
 102 : C ======
 103 :       ELSE IF (NEF.EQ.49) THEN
 103 :         CALL TCOQ4(POINTS,IPINTE)
 104 : C ======
 105 : C  2.16 - Element COQ8
 106 : C ======
 107 :       ELSE IF (NEF.EQ.41) THEN
 107 :         CALL TCOQ8(POINTS,IPINTE)
 108 : C ======
 109 : C  2.17 - Element COQ6
 110 : C ======
 111 :       ELSE IF (NEF.EQ.56) THEN
 111 :         CALL TCOQ6(POINTS,IPINTE)
 112 : C ======
 113 : C  2.18 - Element POI1
 114 : C ======
 115 :       ELSE IF ((NEF.EQ.45).OR.(NEF.EQ.1)) THEN
 115 :         CALL TPOI1(POINTS,IPINTE)
 116 : C ======
 117 : C  2.19 - ERREUR : Element fini non implemente
 118 : C ======
 119 :       ELSE
 119 :         CALL ERREUR(19)
 120 :       ENDIF
 121 : 
 122 :       RETURN
 123 :       END
 124 : 
 125 :  
 126 :  
 127 :  
 128 :  
 129 :