Esope source files

   1 : C MFDKT     SOURCE    CHAT      05/01/13    01:44:32     5004           
   2 :       SUBROUTINE MFDKT(QSI,ETA,XEL,VNM)                                         
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C     EVALUATION DE LA  MATRICE DE MASSE DKT  D.D.L. FLEXION                    
   6 : C                                                                               
   7 : C  ENTREES                                                                      
   8 : C   QSI ETA COORDONNEES REDUITES DU POINT DE GAUSS                              
   9 : C   XEL(3,3)=COORDONNEES LOCALES DE L ELEMENT                                   
  10 : C  SORTIES                                                                      
  11 : C   VNM(9)  =FONCTION DE FORME POUR LA MASSE                                    
  12 : C                 EBERSOLT AOUT 85                                              
  13 : C=======================================================================        
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  16 :       DIMENSION XEL(3,*),VNM(*)                                                 
  17 :       DATA UN,DEUX,TROIS/1.D0,2.D0,3.D0/                                        
  18 :       DATA XZER/0.D0/                                                           
  19 : C                                                                               
  20 : C     MISE A ZERO  DE BB                                                        
  21 : C                                                                               
  22 :       DO 10 IA=1,9                                                              
  23 :       VNM(IA)=XZER                                                              
  24 :   10  CONTINUE                                                                  
  25 : C                                                                               
  26 : C         CALCUL DES FONCTIONS D INTERPOLATION                                  
  27 : C                                                                               
  28 :       AL=UN-QSI-ETA                                                             
  29 :       A=QSI*ETA*AL                                                              
  30 : C                                                                               
  31 :       VNM(1)=(AL**2)*(TROIS-DEUX*AL)+DEUX*A                                     
  32 :       VNM(2)=(AL**2)*QSI+(A/DEUX)                                               
  33 :       VNM(3)=(AL**2)*ETA+(A/DEUX)                                               
  34 :       VNM(4)=(TROIS-DEUX*QSI)*(QSI**2)+DEUX*A                                   
  35 :       VNM(5)=(QSI-UN)*(QSI**2)-A                                                
  36 :       VNM(6)=(QSI**2)*ETA+(A/DEUX)                                              
  37 :       VNM(7)=(TROIS-DEUX*ETA)*(ETA**2)+DEUX*A                                   
  38 :       VNM(8)=(ETA**2)*QSI+(A/DEUX)                                              
  39 :       VNM(9)=(ETA-UN)*(ETA**2)-A                                                
  40 : C                                                                               
  41 : C         VARIABLES DANS LE REPERE X,Y                                          
  42 : C                                                                               
  43 :       X21=XEL(1,2)-XEL(1,1)                                                     
  44 :       X13=XEL(1,1)-XEL(1,3)                                                     
  45 :       Y13=XEL(2,1)-XEL(2,3)                                                     
  46 : C                                                                               
  47 :       C     =VNM(2)*X21-VNM(3)*X13                                              
  48 :       VNM(2)=-VNM(3)*Y13                                                        
  49 :       VNM(3)=-C                                                                 
  50 :       C     =VNM(5)*X21-VNM(6)*X13                                              
  51 :       VNM(5)=-VNM(6)*Y13                                                        
  52 :       VNM(6)=-C                                                                 
  53 :       C     =VNM(8)*X21-VNM(9)*X13                                              
  54 :       VNM(8)=-VNM(9)*Y13                                                        
  55 :       VNM(9)=-C                                                                 
  56 : C                                                                               
  57 :       RETURN                                                                    
  58 :       END                                                                       
  59 :