Esope source files

   1 : C ROTB      SOURCE    CHAT      05/01/13    03:05:22     5004           
   2 :       SUBROUTINE ROTB(B,NSTRS,COSA,SINA)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : ***********************************************************************
   6 : *
   7 : *    TRANSFORMATION DE LA MATRICE B(8,18) POUR QUE LES INCONUES SOIT
   8 : *    EXPRIMEES DANS LE REPERE LOCAL  
   9 : *
  10 : *   ENTREE :
  11 : *    B = MATRICE RELIANT LES DEFORMATIONS AUX DEPLACEMENTS
  12 : *  COSA 
  13 : *  SINA =COS-DIRECTEURS DU 1ER AXE D'ORTHO/AXES LOCAUX    
  14 : *   SORTIE :
  15 : *   B  = B T(TRANSPOSE)
  16 : ***************************************************************************
  17 :       DIMENSION B(NSTRS,*),BB(3,2)
  18 : *
  19 : *    1ER BLOC
  20 : *
  21 :       DO 10 J=1,3
  22 :         JJ1=1+(J-1)*6
  23 :         JJ2=JJ1+1
  24 :         DO 11 I=1,3
  25 :          BB(I,1)=B(I,JJ1)*COSA+B(I,JJ2)*(-SINA)    
  26 :          BB(I,2)=B(I,JJ1)*SINA+B(I,JJ2)*COSA
  27 :  11     CONTINUE
  28 : *
  29 :           DO 12 I=1,3
  30 :            B(I,JJ1)=BB(I,1)
  31 :            B(I,JJ2)=BB(I,2) 
  32 :  12       CONTINUE
  33 :  10   CONTINUE 
  34 : *
  35 : *    2EME BLOC
  36 : *
  37 :       DO 20 J=1,3
  38 :         JJ4=4+(J-1)*6
  39 :         JJ5=JJ4+1
  40 :         DO 30 I=4,6
  41 :          I1=I-3
  42 :          BB(I1,1)=B(I,JJ4)*COSA+B(I,JJ5)*(-SINA)    
  43 :          BB(I1,2)=B(I,JJ4)*SINA+B(I,JJ5)*COSA
  44 :  30     CONTINUE
  45 :         DO 31 I=4,6
  46 :           I1=I-3
  47 :           B(I,JJ4)=BB(I1,1)
  48 :           B(I,JJ5)=BB(I1,2)
  49 :  31     CONTINUE 
  50 : *
  51 :         DO 40 I=7,8
  52 :          I1=I-6 
  53 :          BB(I1,1)=B(I,JJ4)*COSA+B(I,JJ5)*(-SINA)    
  54 :          BB(I1,2)=B(I,JJ4)*SINA+B(I,JJ5)*COSA
  55 :  40     CONTINUE
  56 :         DO 41 I=7,8
  57 :          I1=I-6 
  58 :          B(I,JJ4)=BB(I1,1)
  59 :          B(I,JJ5)=BB(I1,2)
  60 :  41     CONTINUE
  61 :  20   CONTINUE
  62 : *
  63 :       RETURN
  64 :       END 
  65 : 
  66 :