Esope source files

   1 : C DIAXYZ    SOURCE    CHAT      05/01/12    22:50:02     5004           
   2 :       SUBROUTINE DIAXYZ(ITYP,SIGT,SIGEL,VECPRO,IBO,A,B,C,B1,B2)                 
   3 : C                                                                               
   4 : C         EFFECTUE LA TRANSFORMATION DE SIGAMA DIAGONALISE EN                   
   5 : C         SIGAM GLOBAL                                                          
   6 : C                                                                               
   7 : C         ENTREE  SIGT                                                          
   8 : C         SORTIE SIGEL                                                          
   9 : C             ITYP TYPE DE L'ELEMENT                                            
  10 : C                                                                               
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  13 :       DIMENSION SIGT(*),SIGEL(*),VECPRO(*)                                      
  14 : C                                                                               
  15 :       GO TO (70,71,71,73,74,75,76,77,78,79,80,81,71),ITYP                       
  16 : C                                                                               
  17 : C   MODIF DU 05-01-96  
  18 : C   POUR ITYP=13 (COQUES EPAISSES) ON PASSAIT ICI (AU 83)  
  19 : C   DESORMAIS ON VA AU 71 COMME LES CALCULS EN CONTRAINTE PLANE 
  20 : C 
  21 : C   83 CONTINUE                                                               
  22 : C               COQUES AVEC CISAILLEMENTS TRANSVERSES                           
  23 : C      SIGEL(1)=B*(SIGT(1)-SIGT(2))                     
  24 : C      SIGEL(2)=B*(SIGT(1)+SIGT(2))            
  25 : C      SIGEL(3)=SIGT(3)                     
  26 : C      SIGEL(4)=SIGT(4)                                 
  27 : C      SIGEL(5)=SIGT(5)                               
  28 : C      SIGEL(6)=SIGT(6)                     
  29 : C      GO TO 740                           
  30 :   81  CONTINUE                                                                  
  31 : C                      TUYAUX                                                   
  32 :   80  CONTINUE                                                                  
  33 : C                      POUTRES                                                  
  34 :       DO 810 I=1,6                                                              
  35 : 810   SIGEL(I)=SIGT(I)                                                          
  36 :       GO TO 740                                                                 
  37 :    79 CONTINUE                                                                  
  38 :       GO TO 740                                                                 
  39 :    78 CONTINUE                                                                  
  40 :       GO TO 740                                                                 
  41 :    77 SIGEL(1)=SIGT(1)                                                          
  42 :       GO TO 740                                                                 
  43 :    76 DO 760 IA=1,4                                                             
  44 :   760 SIGEL(IA)=SIGT(IA)                                                        
  45 :       GO TO 740                                                                 
  46 :    75 CONTINUE                                                                  
  47 :       SIGEL(1)=B*(SIGT(1)-SIGT(2))                                              
  48 :       SIGEL(2)=B*(SIGT(1)+SIGT(2))                                              
  49 :       SIGEL(4)=SIGT(4)                                                          
  50 :       SIGEL(3)=SIGT(3)                                                          
  51 :       GO TO 740                                                                 
  52 :    74 CONTINUE                                                                  
  53 :       DO 7401 IB=1,3                                                            
  54 :       SIGEL(IB)=0.                                                              
  55 :       DO 7402 IC=1,3                                                            
  56 :       IP=3*(IC-1)+IB                                                            
  57 :  7402 SIGEL(IB)=SIGEL(IB)+VECPRO(IP)*SIGT(IC)                                   
  58 :  7401 SIGEL(3+IB)=SIGT(3+IB)                                                    
  59 :       GO TO 740                                                                 
  60 :    71 CONTINUE                                                                  
  61 :       B1=B                                                                      
  62 :       B2=B                                                                      
  63 :       DO 710 I=1,IBO                                                            
  64 :       IP=(I-1)*3                                                                
  65 :       SIGEL(IP+1)=B1*SIGT(IP+1)-B2*SIGT(IP+2)                                   
  66 :       SIGEL(IP+2)=B2*SIGT(IP+1)+B1*SIGT(IP+2)                                   
  67 :   710 SIGEL(3+IP)=SIGT(3+IP)                                                    
  68 :       GO TO 740                                                                 
  69 :    73 CONTINUE                                                                  
  70 :       SIGEL(3)=SIGT(3)                                                          
  71 :       GO TO 740                                                                 
  72 :    70 CONTINUE                                                                  
  73 :       SIGT(1)=SIGT(1)*A                                                         
  74 :       SIGT(2)=SIGT(2)*B                                                         
  75 :       SIGT(3)=SIGT(3)*C                                                         
  76 :       SIGEL(1)=SIGT(1)+SIGT(2)+SIGT(3)                                          
  77 :       SIGEL(2)=SIGT(1)-SIGT(2)+SIGT(3)                                          
  78 :       SIGEL(3)=SIGT(1)-2.*SIGT(3)                                               
  79 :       SIGEL(4)=SIGT(4)                                                          
  80 :       SIGEL(5)=SIGT(5)                                                          
  81 :       SIGEL(6)=SIGT(6)                                                          
  82 :   740 CONTINUE                                                                  
  83 :       RETURN                                                                    
  84 :       END                                                                       
  85 :  
  86 :