Esope source files

   1 : C DBDKT2    SOURCE    CHAT      05/02/28    21:15:01     5037           
   2 :       SUBROUTINE DBDKT2(XEL,DD,XDDL,IGAU,XSTRS,SHP1,SHPTOT,BGENE,
   3 :      *  NBNO,LRE,NST)                                 
   4 : C======================================================================         
   5 : C                                                                               
   6 : C  AJOUTE DANS LES CONTRAINTES XSTRS LA PARTIE DHO * EPS2                       
   7 : C                  QUE POUR LES ELEMENTS DKT                                    
   8 : C  ENTREE                                                                       
   9 : C    XEL(3,3)           =COORDONNEES LOCALES  DU DKT                            
  10 : C    DD(6,6)            =MATRICE DE HOOKE                                       
  11 : C    XDDL(18)          =D D L LOCAUX DE LA FORMULATION                          
  12 : C    IGAU               =NUMERO DU POINT DE GAUSS    DE 1 A 3                   
  13 : C    XSTRS(6)           =CONTRAINTES  DHO*(EPS1        )                        
  14 : C  SORTIE                                                                       
  15 : C    XSTRS(6)           =CONTRAINTES DHO*(EPS1 + EPS2 )                         
  16 : C                                                                               
  17 : C   CODE EBERSOLT AOUT  85                                                      
  18 : C=======================================================================        
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  21 :       DIMENSION XEL(3,*),DD(6,*),XSTRS(*),XDDL(*)                               
  22 :       DIMENSION EPS2(6),SHP(6,9),ROT(9)
  23 :       DIMENSION SHP1(6,NBNO),SHPTOT(6,NBNO,*)
  24 :       DIMENSION BGENE(NST,*)                                         
  25 : C                                                                               
  26 : C                                                                               
  27 :       DATA XZER,UNDEMI/0.D0,.5D0/                                               
  28 : C                                                                               
  29 :       DUDX=XZER
  30 :       DUDY=XZER
  31 :       DVDX=XZER
  32 :       DVDY=XZER                                                                          
  33 :       DWDX=XZER                                                                 
  34 :       DWDY=XZER
  35 :       eps2(4)=xzer  
  36 :       eps2(5)=xzer  
  37 :       eps2(6)=xzer                                                                                                                             
  38 : C                                                                               
  38 :       CALL DKTSHP(IGAU,XEL,SHP,DJAC)                                            
  39 : C                                                                               
  40 : C     ON EXTRAIT LES DEGRES W ROT X ROT Y                                       
  41 : C                                                                               
  42 :       IACC= 0                                                                   
  43 :       IADD= 2                                                                   
  44 :       DO 100 IA=1,3                                                             
  45 :       DO 110 IB=1,3                                                             
  46 :       ROT(IACC+IB)=XDDL(IADD+IB)                                                
  47 :  110  CONTINUE                                                                  
  48 :       IACC=IACC+3                                                               
  49 :       IADD=IADD+6                                                               
  50 :  100  CONTINUE                                                                  
  51 : C                                                                               
  52 :       DO 200 IA=1,9                                                             
  53 :       DWDX =SHP(2,IA)*ROT(IA) +DWDX                                             
  54 :       DWDY =SHP(3,IA)*ROT(IA) +DWDY                                             
  55 :   200 CONTINUE                                                                  
  56 : C                                                                               
  57 : C--Calcul des derivees de u et v
  58 :       DO 300 IA=0,2
  59 :       DUDX=DUDX + BGENE(1,1+6*IA)*XDDL(1+6*IA)
  60 :       DUDY=DUDY + BGENE(3,1+6*IA)*XDDL(1+6*IA)
  61 :       DVDX=DVDX + BGENE(3,2+6*IA)*XDDL(2+6*IA)
  62 :       DVDY=DVDY + BGENE(2,2+6*IA)*XDDL(2+6*IA)
  63 :   300 CONTINUE
  64 : C
  65 : C                                               
  66 :       EPS2(1)=0.5D0*(DUDX*DUDX+DVDX*DVDX+DWDX*DWDX)                                      
  67 :       EPS2(2)=0.5D0*(DUDY*DUDY+DVDY*DVDY+DWDY*DWDY)                                      
  68 :       EPS2(3)=DUDX*DUDY+DVDX*DVDY+DWDX*DWDY                                                                                                                                         
  69 : C                                                                               
  70 : C     MULTIPLICATION PAR LA MATRICE DE HOOKE                                    
  71 : C                                                                                   
  72 :       DO 700 IB=1,6                                                             
  73 :       DO 700 IC=1,6                                                             
  74 :       XSTRS(IB)=XSTRS(IB)+ DD(IB,IC)*EPS2(IC)                                   
  75 :   700 CONTINUE                                                           
  76 : C                                                                               
  77 :       RETURN                                                                    
  78 :       END                                                                       
  79 :  
  80 :