Esope source files

   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 :