Esope source files

   1 : C TRANSK    SOURCE    CHAT      05/01/13    03:45:07     5004           
   2 :       SUBROUTINE TRANSK(VKE,B,LRE,NBNO,ICOM)                                    
   3 : C=====================================================                          
   4 : C                                                                               
   5 : C     TRANSFORME VKE(LRE,LRE) DU REPERE LOCAL VKE SYMETRIQUE                    
   6 : C     AU REPERE GENERAL A L AIDE DE LA MATRICE DE PASSAGE                       
   7 : C     B(3,3)                                                                    
   8 : C     NBNO NOMBRE DE NOEUDS DE L ELEMENT 3 NOEUDS DKT                           
   9 : C                                        4 NOEUDS DKQ                           
  10 : C     LRE NOMBRE DE DDL DE L ELEMENT                                            
  11 : C     6 DDL PAR NOEUD                                                           
  12 : C  ENTREES                                                                      
  13 : C    VKE(LRE,LRE)  = MATRICE DE RAIDEUR APRES CHGMNT DE REPERE                  
  14 : C    B(3,3)   = MATRICE DE PASSAGE                                              
  15 : C    LRE      = NOMBRE  DE D.D.L. DANS LA MATRICE DE RIGIDITE                   
  16 : C    NBNO     = NOMBRE DE NOEUDS ( LRE = NBNO*6                                 
  17 : C    ICOM     = 0 SI LES D.D.L. DE FLEXION SONT DECOUPLES DES                   
  18 : C                        D.D.L. DE MEMBRANE  1 SINON                            
  19 : C  SORTIES                                                                      
  20 : C    VKE(LRE,LRE)  = MATRICE DE RAIDEUR APRES CHGMNT DE REPERE                  
  21 : C                                                                               
  22 : C       EBERSOLT AOUT 85                                                        
  23 : C=====================================================                          
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  26 :       DIMENSION VKE(LRE,*),B(3,*)                                               
  27 :       DIMENSION AA(3,3),A(3,3)                                                  
  28 :       DIMENSION IMLT(36)                                                        
  29 :       DATA IMLT/1,2,1,4,3,1,2,4,2,1,4,3,4,3,1,2,4,2,4,2,1,                      
  30 :      1    4,3,4,3,4,3,1,2,4,2,4,2,4,2,1/                                        
  31 :       NB=NBNO+NBNO                                                              
  32 :       IJ=1                                                                      
  33 :       JTROIS=0                                                                  
  34 :       DO 100 J=1,NB                                                             
  35 :       ITROIS=0                                                                  
  36 :       DO 200 I=1,J                                                              
  37 :       ICH=IMLT(IJ)                                                              
  38 :       IF(ICOM.EQ.1) ICH=4                                                       
  39 : C     TRANSVASE UN BLOC 3 3 DE VKE DANS LA MATRICE DE TRAVAIL A                 
  40 :       DO 300 IA=1,3                                                             
  41 :       DO 300 IB=1,3                                                             
  42 :       A(IA,IB)=VKE(ITROIS+IA,JTROIS+IB)                                         
  43 :   300 CONTINUE                                                                  
  44 : C     SELECTION DU TYPE DE CHANGEMENT DE REPERE                                 
  45 : C     TRANSFORME A EN AA A L AIDE DE B                                          
  46 :       GOTO (310,320,330,340) ICH                                                
  46 :   310 CALL VMULT1(A,AA,B)                                                       
  47 :       GOTO 350                                                                  
  47 :   320 CALL VMULT2(A,AA,B)                                                       
  48 :       GOTO 350                                                                  
  48 :   330 CALL VMULT3(A,AA,B)                                                       
  49 :       GOTO 350                                                                  
  49 :   340 CALL VMULT4(A,AA,B)                                                       
  50 :   350 CONTINUE                                                                  
  51 : C     TRANSVASE LA MATRICE AA DANS LE BLOC 3 3 DE VKE                           
  52 : C     ET SYMMETRISE LE BLOC TRANSPOSE DE VKE CORRESPONDANT                      
  53 :       DO 400 IA=1,3                                                             
  54 :       DO 400 IB=1,3                                                             
  55 :       VKE(ITROIS+IA,JTROIS+IB)=AA(IA,IB)                                        
  56 :       VKE(JTROIS+IB,ITROIS+IA)=AA(IA,IB)                                        
  57 :   400 CONTINUE                                                                  
  58 :       IJ=IJ+1                                                                   
  59 :   200 ITROIS=ITROIS+3                                                           
  60 :   100 JTROIS=JTROIS+3                                                           
  61 :       RETURN                                                                    
  62 :       END                                                                       
  63 :