Sources Esope | Cast3M

rotb
C ROTB      SOURCE    CHAT      05/01/13    03:05:22     5004
      SUBROUTINE ROTB(B,NSTRS,COSA,SINA)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
***********************************************************************
*
*    TRANSFORMATION DE LA MATRICE B(8,18) POUR QUE LES INCONUES SOIT
*    EXPRIMEES DANS LE REPERE LOCAL
*
*   ENTREE :
*    B = MATRICE RELIANT LES DEFORMATIONS AUX DEPLACEMENTS
*  COSA
*  SINA =COS-DIRECTEURS DU 1ER AXE D'ORTHO/AXES LOCAUX
*   SORTIE :
*   B  = B T(TRANSPOSE)
***************************************************************************
      DIMENSION B(NSTRS,*),BB(3,2)
*
*    1ER BLOC
*
      DO 10 J=1,3
        JJ1=1+(J-1)*6
        JJ2=JJ1+1
        DO 11 I=1,3
         BB(I,1)=B(I,JJ1)*COSA+B(I,JJ2)*(-SINA)
         BB(I,2)=B(I,JJ1)*SINA+B(I,JJ2)*COSA
 11     CONTINUE
*
          DO 12 I=1,3
           B(I,JJ1)=BB(I,1)
           B(I,JJ2)=BB(I,2)
 12       CONTINUE
 10   CONTINUE
*
*    2EME BLOC
*
      DO 20 J=1,3
        JJ4=4+(J-1)*6
        JJ5=JJ4+1
        DO 30 I=4,6
         I1=I-3
         BB(I1,1)=B(I,JJ4)*COSA+B(I,JJ5)*(-SINA)
         BB(I1,2)=B(I,JJ4)*SINA+B(I,JJ5)*COSA
 30     CONTINUE
        DO 31 I=4,6
          I1=I-3
          B(I,JJ4)=BB(I1,1)
          B(I,JJ5)=BB(I1,2)
 31     CONTINUE
*
        DO 40 I=7,8
         I1=I-6
         BB(I1,1)=B(I,JJ4)*COSA+B(I,JJ5)*(-SINA)
         BB(I1,2)=B(I,JJ4)*SINA+B(I,JJ5)*COSA
 40     CONTINUE
        DO 41 I=7,8
         I1=I-6
         B(I,JJ4)=BB(I1,1)
         B(I,JJ5)=BB(I1,2)
 41     CONTINUE
 20   CONTINUE
*
      RETURN
      END