Sources Esope | Cast3M

asksig
C ASKSIG    SOURCE    CHAT      07/10/22    21:15:09     5921
      SUBROUTINE ASKSIG(S,B,XK)
      IMPLICIT REAL*8(A-H,O-Z)
      IMPLICIT INTEGER (I-N)
 
*      SUBROUTINE ASSE_KSIGMA(S,B,K)
*      REAL *8 S(6),B(3,8),K(8,8)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C ANTOINE LEGAY LE 13/02/2001
C CALCUL ET ASSEMBLE LA MATRICE 8*8 K_SIGMA ---> K
C RQ: ELLE EST DIAGONALE PAR BLOCS DE 8*8, EGAUX
C B : DERIVEES DES FONCTIONS DE FORMES
C S : CONTRAINTES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
 
      DIMENSION  S(6),B(3,8),XK(8,8)
      DO J=1,8
         DO I=1,8
            XK(I,J)=0.
         ENDDO
      ENDDO
C      WRITE(6,*)'*************** SIGMA DANS ASSE-KSIGMA*************'
C      WRITE(6,*)(S(I),I=1,6)
C      WRITE(6,*)'*******************************'
C      WRITE(6,*)'*************** B ****************'
C      WRITE(6,*)(B(1,I),I=1,8)
C      WRITE(6,*)(B(2,I),I=1,8)
C      WRITE(6,*)(B(3,I),I=1,8)
C      WRITE(6,*)'*******************************'
C
      DO J=1,8
         DO I=1,8
            XK(I,J) =  S(1)*B(1,I)*B(1,J)
     !              + S(2)*B(2,I)*B(2,J)
     !              + S(3)*B(3,I)*B(3,J)
     !              + S(4)*(B(1,I)*B(2,J)+B(2,I)*B(1,J))
     !              + S(6)*(B(1,I)*B(3,J)+B(3,I)*B(1,J))
     !              + S(5)*(B(3,I)*B(2,J)+B(2,I)*B(3,J))
C            WRITE(6,*)'K IJ=',K(I,J),I,J
         ENDDO
      ENDDO
C
      END