Télécharger asksig.eso

Retour à la liste

Numérotation des lignes :

  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.  
  46.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales