Esope source files

   1 : C DBCO33    SOURCE    CHAT      11/01/31    21:15:01     6862           
   2 :       SUBROUTINE DBCO33(XE,DD,XDDL,EPS2,SHP,BPSS,XEL,XLOC,XSTRS)                
   3 : C======================================================================         
   4 : C                                                                               
   5 : C  AJOUTE DANS LES CONTRAINTES XSTRS LA PARTIE DHO * EPS2                       
   6 : C                  QUE POUR LES ELEMENTS COQ3                                   
   7 : C  ENTREE                                                                       
   8 : C    XE(3,3)           =COORDONNEES GLOBALES DU COQ3                            
   9 : C    DD(6,6)            =MATRICE DE HOOKE                                       
  10 : C    XDDL(18)          =D.D.L.GLOBAUX DE LA FORMULATION                         
  11 : C    XSTRS(6)           =CONTRAINTES  DHO*(EPS1        )                        
  12 : C  TRAVAIL                                                                      
  13 : C    EPS2(3)            = EPSILON 2                                             
  14 : C    SHP(3,9)           = FONCTIONS DE FORME                                    
  15 : C    XEL(3,3)           = COORDONNEES LOCALES                                   
  16 : C    XLOC(18)         = D.D.L. LOCAUX                                           
  17 : C  SORTIE                                                                       
  18 : C    XSTRS(6)           =CONTRAINTES DHO*(EPS1 + EPS2 )                         
  19 : C                                                                               
  20 : C   CODE EBERSOLT AOUT  85                                                      
  21 : C=======================================================================        
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  24 :       DIMENSION XE(3,*),DD(6,*),XSTRS(*),XDDL(*)                                
  25 :       DIMENSION EPS2(*),SHP(3,*),BPSS(3,*)                                      
  26 :       DIMENSION XLOC(*),XEL(3,*)                                                
  27 : C                                                                               
  28 : C                                                                               
  29 :       DATA XZER,UNDEMI,UNTIER/0.D0,.5D0,.333333333333333333D0/                  
  30 : C                                                                               
  31 :       DWDX=XZER                                                                 
  32 :       DWDY=XZER                                                                 
  33 : C                                                                             
  33 :       CALL VPAST(XE,BPSS)                                                       
  33 :       CALL VCORLC(XE,XEL,BPSS) 
  34 : *      WRITE(6,*) ' BPSS', ((BPSS( IOU,IYU),IOU=1,3),IYU=1,3) 
  35 : *      WRITE(6,*) 'XEL' , ((XEL(iou,iyu),iou=1,3), iyu=1,3)                                               
  35 :       CALL MATVEC(XDDL,XLOC,BPSS,6)                                             
  35 :       CALL COQ3SH(UNTIER,UNTIER,XEL,SHP)   
  36 : *      write(6,*) ((shp(iou,iyu), iou=1,3),iyu=1,9) 
  37 : *      write(6,*) (xloc(iou),iou=1,18)                                   
  38 :       DO 200 IA=1,3                                                             
  39 :       IA1 = 2 +(IA-1)*6                                                         
  40 :       IA2 =    (IA-1)*3                                                         
  41 :       DO 200 IB=1,3                                                             
  42 :       DWDX =SHP(2,IA2+IB)*XLOC(IA1+IB) +DWDX                                    
  43 :       DWDY =SHP(3,IA2+IB)*XLOC(IA1+IB) +DWDY                                    
  44 :   200 CONTINUE                                                                  
  45 : C                                                                               
  46 :       EPS2(1)= DWDX*DWDX*UNDEMI                                                 
  47 :       EPS2(2)= DWDY*DWDY*UNDEMI                                                 
  48 :       EPS2(3)= DWDX*DWDY                                                        
  49 : **      WRITE(6,888) EPS2(1),EPS2(2),EPS2(3),DWDX,DWDY                            
  50 : *  888 FORMAT(4X,' DEBUGG EPS2 ',5E12.5)                                         
  51 : C                                                                               
  52 : C     MULTIPLICATION PAR LA MATRICE DE HOOKE                                    
  53 : C                                                                               
  54 :       DO 300 IB=1,3                                                             
  55 :       DO 300 IC=1,3                                                             
  56 :       XSTRS(IB)=XSTRS(IB)+ DD(IB,IC)*EPS2(IC)                                   
  57 :   300 CONTINUE                                                                  
  58 : C                                                                               
  59 :       RETURN                                                                    
  60 :       END                                                                       
  61 :  
  62 :