Esope source files

   1 : C DKTHSH    SOURCE    CHAT      05/02/25    21:15:10     5035           
   2 :       SUBROUTINE DKTHSH(tabuv,tabw,tabrot,DJAC,XSTRS,REL,epai)                                     
   3 : C---------------------------------------------------------------------          
   4 : C   CALCULE LA MATRICE  KSIGMA   POUR L ELEMENT DKT                             
   5 : C                                                                               
   6 : C      ENTREES                                                                  
   7 : C         tabuv(6,3)    = FONCTIONS U U,X U,Y V V,X V,Y
   8 : C         tabw(6,9)   = FONCTIONS W W,X  W,Y
   9 : C         tabrot(4,9) = FONCTIONS BETAx,x BETAx,y BETA y,x BETAy,y                           
  10 : C         DJAC        = JACOBIEN                                                
  11 : C         XSTRS(6)    = CONTRAINTES AU POINT DE GAUSS                           
  12 : C         REL (18,18) = MATRICE KSIGMA                                          
  13 : C         epai        = Epaisseur
  14 : C      SORTIE                                                                   
  15 : C         REL (18,18) = MATRICE KSIGMA                                          
  16 : C                                                                               
  17 : C   RE-ECRIT  LE 2/4/90                                                         
  18 : C---------------------------------------------------------------------          
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  21 :       DIMENSION REL(18,*),XSTRS(*),tabuv(6,*),tabw(6,9),tabrot(4,9)                                     
  22 :       DATA UNTIER/.33333333333333333D0/ 
  23 :       DATA UNDOUZ/.08333333333333333D0/
  24 :       xNXX=XSTRS(1)                                                            
  25 :       xNYY=XSTRS(2)                                                            
  26 :       xNXY=XSTRS(3)
  27 :       xNXX2=XSTRS(1) * UNDOUZ * epai * epai                                                           
  28 :       xNYY2=XSTRS(2) * UNDOUZ * epai * epai                                                           
  29 :       xNXY2=XSTRS(3) * UNDOUZ * epai * epai                                                           
  30 :                                                                               
  31 :       DO 16  J=1,3                                                              
  32 :       JJ=(J-1)*6+1                                                              
  33 :       DO 17  I=1,3                                                              
  34 :       II=(I-1)*6+1                                                              
  35 :         S= xNXX*tabuv(2,I)*tabuv(2,J) +                                            
  36 :      .     xNYY*tabuv(3,I)*tabuv(3,J) +                                            
  37 :      .     xNXY*(tabuv(2,I)*tabuv(3,J)+tabuv(3,I)*tabuv(2,J))                          
  38 :       S=S*DJAC                                                                  
  39 :       REL(II,JJ)    =REL(II,JJ)    +S                                           
  40 :       REL(II+1,JJ+1)=REL(II+1,JJ+1)+S                                                         
  41 :    17 CONTINUE                                                                  
  42 :    16 CONTINUE                                                                  
  43 : 
  44 :       DO 20 L=1,3
  45 :       DO 21 M=1,3
  46 :       DO 22 I=1,3
  47 :       DO 23 K=1,3
  48 :          S= xNXX*tabw(2,3*L-3+I)*tabw(2,3*M-3+K) +
  49 :      .      xNYY*tabw(3,3*L-3+I)*tabw(3,3*M-3+K) +
  50 :      .      xNXY*tabw(2,3*L-3+I)*tabw(3,3*M-3+K) +
  51 :      .      xNXY*tabw(3,3*L-3+I)*tabw(2,3*M-3+K) +
  52 :      .      xNXX2*tabrot(1,3*L-3+I)*tabrot(1,3*M-3+K) +
  53 :      .      xNXX2*tabrot(3,3*L-3+I)*tabrot(3,3*M-3+K) +
  54 :      .      xNYY2*tabrot(2,3*L-3+I)*tabrot(2,3*M-3+K) +
  55 :      .      xNYY2*tabrot(4,3*L-3+I)*tabrot(4,3*M-3+K) +     
  56 :      .      xNXY2*tabrot(1,3*L-3+I)*tabrot(2,3*M-3+K) +
  57 :      .      xNXY2*tabrot(2,3*L-3+I)*tabrot(1,3*M-3+K) +
  58 :      .      xNXY2*tabrot(3,3*L-3+I)*tabrot(4,3*M-3+K) +
  59 :      .      xNXY2*tabrot(4,3*L-3+I)*tabrot(3,3*M-3+K)
  60 :          S=S*DJAC
  61 :          REL(6*L-4+I,6*M-4+K) = REL(6*L-4+I,6*M-4+K) + S 
  62 :    23 CONTINUE
  63 :    22 CONTINUE
  64 :    21 CONTINUE
  65 :    20 CONTINUE
  66 : 
  67 :       RETURN                                                                    
  68 :       END                                                                       
  69 :