Télécharger la plaquette Cast3M


fr en

Télécharger shp1ce.eso

Retour à la liste

Numérotation des lignes : 

shp1ce

  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.  
  66.  

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