Esope source files

   1 : C ASKSIG    SOURCE    CHAT      07/10/22    21:15:09     5921           
   2 :       SUBROUTINE ASKSIG(S,B,XK)
   3 :       IMPLICIT REAL*8(A-H,O-Z)
   4 :       IMPLICIT INTEGER (I-N)
   5 : 
   6 : *      SUBROUTINE ASSE_KSIGMA(S,B,K)
   7 : *      REAL *8 S(6),B(3,8),K(8,8)
   8 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC      
   9 : C ANTOINE LEGAY LE 13/02/2001
  10 : C CALCUL ET ASSEMBLE LA MATRICE 8*8 K_SIGMA ---> K
  11 : C RQ: ELLE EST DIAGONALE PAR BLOCS DE 8*8, EGAUX
  12 : C B : DERIVEES DES FONCTIONS DE FORMES
  13 : C S : CONTRAINTES
  14 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
  15 : 
  16 :       DIMENSION  S(6),B(3,8),XK(8,8)    
  17 :       DO J=1,8
  18 :          DO I=1,8
  19 :             XK(I,J)=0.
  20 :          ENDDO
  21 :       ENDDO
  22 : C      WRITE(6,*)'*************** SIGMA DANS ASSE-KSIGMA*************'
  23 : C      WRITE(6,*)(S(I),I=1,6)
  24 : C      WRITE(6,*)'*******************************'
  25 : C      WRITE(6,*)'*************** B ****************'
  26 : C      WRITE(6,*)(B(1,I),I=1,8)
  27 : C      WRITE(6,*)(B(2,I),I=1,8)
  28 : C      WRITE(6,*)(B(3,I),I=1,8)
  29 : C      WRITE(6,*)'*******************************'
  30 : C
  31 :       DO J=1,8
  32 :          DO I=1,8
  33 :             XK(I,J) =  S(1)*B(1,I)*B(1,J)
  34 :      !              + S(2)*B(2,I)*B(2,J)
  35 :      !              + S(3)*B(3,I)*B(3,J)
  36 :      !              + S(4)*(B(1,I)*B(2,J)+B(2,I)*B(1,J))
  37 :      !              + S(6)*(B(1,I)*B(3,J)+B(3,I)*B(1,J))
  38 :      !              + S(5)*(B(3,I)*B(2,J)+B(2,I)*B(3,J))
  39 : C            WRITE(6,*)'K IJ=',K(I,J),I,J
  40 :          ENDDO
  41 :       ENDDO
  42 : C
  43 :       END
  44 : 
  45 :