Télécharger la plaquette Cast3M


fr en

Télécharger shrot3.eso

Retour à la liste

Numérotation des lignes : 

shrot3

  1. C SHROT3    SOURCE    CHAT      05/01/13    03:15:49     5004
  2.       SUBROUTINE SHROT3(QSI,ETA,DZE,IELE,SHP,IRET)
  3. C=======================================================================
  4. C
  5. C     FONCTIONS DE FORME VECTORIELLES POUR L'ELEMENT ROT3 :
  6. C       ELLES SONT LES ROTATIONNELS SURFACIQUES DES FONCTIONS
  7. C       DE LAGRANGE DE L'ELEMENT TRI3 :
  8. C         Wi = Grad(PHIi) x n ou n est la normale exterieure au plan de
  9. C                                l'element
  10. C         Wi = Grad(PHIi) tourne de -pi/2
  11. C         Wi est parallele au cote oppose au sommet i dans le plan de
  12. C                                l'element
  13. C
  14. C=======================================================================
  15.       IMPLICIT INTEGER(I-N)
  16.       IMPLICIT REAL*8(A-H,O-Z)
  17.       PARAMETER(XZER=0.D0,UNDEMI=.5D0,UN=1.D0,DEUX=2.D0)
  18.       PARAMETER(TROIS=3.D0,QUATRE=4.D0)
  19.       DIMENSION SHP(6,*)
  20. C=======================================================================
  21. C     QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS
  22. C
  23. C     SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE LAGRANGE DU ND INO
  24. C                            SHP(2,INO) COORDONNEE X DE LA FONCTION DU ND INO
  25. C                            SHP(3,INO) COORDONNEE Y DE LA FONCTION DU ND INO
  26. C                            SHP(4,INO) COORDONNEE Z DE LA FONCTION DU ND INO
  27. C=======================================================================
  28.       IRET=1
  29. C
  30. C    66 INDIQUE QUE L ELEMENT PAS ENCORE IMPLEMENTE
  31. C
  32. C           P1 S2 S3 T3 T4 T6 T7 Q4 Q5 Q8 Q9 R2 R3 C8 C20 P6 P15
  33.       GOTO (66,66,66, 1,66,66,66,66,66,66,66,66,66,66, 66,66, 66,
  34. C        L3 L4 L6 L8 MU T4 T10 P5 P13 ATT SUP RP3 LP6 LP8
  35.      &   66,66,66,66,66,66, 66,66, 66, 66, 66, 66, 66, 66) IELE
  36. C
  37.   66  IRET=0
  38.       GOTO 666
  39. C
  40. C
  41. C     FONCTIONS DE FORME BIDIMENSIONNELLES
  42. C
  43.    1  CONTINUE
  44. C
  45. C     FONCTIONS DE FORME TRIANGLE A 3 NOEUDS
  46. C      NOEUD1 = (0,0)
  47. C      NOEUD2 = (1,0)
  48. C      NOEUD3 = (0,1)
  49. C
  50.       AUX=UN-QSI-ETA
  51.       SHP(1,1)=AUX
  52.       SHP(2,1)=-UN
  53.       SHP(3,1)=UN
  54.       SHP(1,2)=QSI
  55.       SHP(2,2)=XZER
  56.       SHP(3,2)=-UN
  57.       SHP(1,3)=ETA
  58.       SHP(2,3)=UN
  59.       SHP(3,3)=XZER
  60. C
  61.  666  CONTINUE
  62.       RETURN
  63.       END
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70.  

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