Esope source files

   1 : C SHP1CE    SOURCE    MAGN      05/02/18    21:16:23     5031           
   2 :       SUBROUTINE SHP1CE(QSI,ETA,DZE,MELE,SHP,IRET)
   3 : C=======================================================================
   4 : C
   5 : C     FONCTIONS DE FORME ET LEUR DERIVEES  CENTREP1
   6 : C
   7 : C=======================================================================
   8 :       IMPLICIT REAL*8(A-H,O-Z)
   9 :       PARAMETER(UNDEMI=.5D0,DEUX=2.D0)
  10 :       PARAMETER(TROIS=3.D0,QUATRE=4.D0)
  11 :       DIMENSION SHP(6,*)
  12 : C=======================================================================
  13 : C     QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS
  14 : C
  15 : C     SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE FORME DU ND INO
  16 : C                            SHP(2,INO) DERIVEES PR RPPRT A QSI
  17 : C                            SHP(3,INO) DERIVEES PR RPPRT A ETA
  18 : C                            SHP(4,INO) DERIVEES PR RPPRT A DZE
  19 : C                            SHP(5,INO) DERIVEES PR RPPRT A 4
  20 : C                               COORDONNEES BARYCENTRIQUES
  21 : C     CAS DE DIMENSION 2
  22 : C               3 NOEUDS TRIANGLE
  23 : C=======================================================================
  24 :       IRET=1
  25 : 
  26 : C    CAS QUAF/CENTREP1 pour un element geometrique de type TRI7
  27 : C    CAS MACRO/CENTREP1 pour un element geometrique de type TRI7
  28 :       IF(MELE.EQ.238.OR.MELE.EQ.245.OR.MELE.EQ.252
  29 :      &    .OR.MELE.EQ.224.OR.MELE.EQ.231) THEN
  30 :       
  31 :       SHP(1,1)=(5.0D0/3.0D0)-DEUX*QSI-DEUX*ETA
  32 :       SHP(1,2)=DEUX*QSI-1.0D0/3.0D0
  33 :       SHP(1,3)=DEUX*ETA-1.0D0/3.0D0
  34 :       
  35 :       SHP(2,1)=-DEUX
  36 :       SHP(2,2)=DEUX
  37 :       SHP(2,3)=0.0D0
  38 :       
  39 :       SHP(3,1)=-DEUX
  40 :       SHP(3,2)=0.0D0
  41 :       SHP(3,3)=DEUX
  42 :       
  43 : C    CAS QUAF/CENTREP1 pour un element geometrique de type QUA9
  44 : C ET CAS MACRO/CENTREP1 pour un element geometrique de type QUA9
  45 :       ELSEIF(MELE.EQ.239.OR.MELE.EQ.246.OR.MELE.EQ.253
  46 :      &   .OR.MELE.EQ.225.OR.MELE.EQ.232) THEN
  47 :       
  48 :       SHP(1,1)=DEUX/TROIS*(ETA+QSI+UNDEMI)
  49 :       SHP(1,2)=DEUX/TROIS*(ETA-(DEUX*QSI)+UNDEMI)
  50 :       SHP(1,3)=DEUX/TROIS*(QSI-(DEUX*ETA)+UNDEMI)
  51 :       
  52 :       SHP(2,1)=DEUX/TROIS
  53 :       SHP(2,2)=-QUATRE/TROIS
  54 :       SHP(2,3)=DEUX/TROIS
  55 :       
  56 :       SHP(3,1)=DEUX/TROIS
  57 :       SHP(3,2)=DEUX/TROIS
  58 :       SHP(3,3)=-QUATRE/TROIS
  59 :       
  60 :       ENDIF
  61 :       
  62 :       RETURN
  63 :       END
  64 : 
  65 :