Télécharger la plaquette Cast3M


fr en

Télécharger cubebb.eso

Retour à la liste

Numérotation des lignes : 

cubebb

  1. C CUBEBB    SOURCE    CHAT      07/10/22    21:15:33     5921
  2.       SUBROUTINE CUBEBB(X,B,VOL)
  3.       implicit real*8(A-H,O-Z)
  4.       implicit integer (I-N)
  5.  
  6. *      SUBROUTINE CUBE_BB(X,B,VOL)
  7. *
  8. *
  9. *      IMPLICIT NONE
  10. *
  11. *
  12. *
  13. *      REAL *8 UNS12
  14. *
  15. *      REAL *8 X(24)
  16. *
  17. *      REAL *8 B(3,8),VOL
  18. *
  19. *
  20. *
  21. *      INTEGER I,J
  22.  
  23.       DIMENSION X(24),B(3,8)
  24.       UNS12=1.D0/12.D0
  25. *
  26. ******* CALCUL DES MATRICES B_BARRE
  27. *
  28.       B(1,1)= X(5)*(X(18)-X(9)+X(15)-X(12))+X(8)*(X(6)-X(12))+
  29.      & X(11)*(X(9)-X(24)+X(6)-X(15))+X(14)*(X(24)-X(18)+X(12)-X(6))+
  30.      & X(17)*(X(15)-X(6))+X(23)*(X(12)-X(15))
  31. *
  32.       B(1,2)= X(8)*(X(21)-X(12)+X(18)-X(3))+X(11)*(X(9)-X(3))+
  33.      & X(2)*(X(12)-X(15)+X(9)-X(18))+X(17)*(X(15)-X(21)+X(3)-X(9))+
  34.      & X(20)*(X(18)-X(9))+X(14)*(X(3)-X(18))
  35. *
  36.       B(1,3)= X(11)*(X(24)-X(3)+X(21)-X(6))+X(2)*(X(12)-X(6))+
  37.      & X(5)*(X(3)-X(18)+X(12)-X(21))+X(20)*(X(18)-X(24)+X(6)-X(12))+
  38.      & X(23)*(X(21)-X(12))+X(17)*(X(6)-X(21))
  39. *
  40.       B(1,4)= X(2)*(X(15)-X(6)+X(24)-X(9))+X(5)*(X(3)-X(9))+
  41.      & X(8)*(X(6)-X(21)+X(3)-X(24))+X(23)*(X(21)-X(15)+X(9)-X(3))+
  42.      & X(14)*(X(24)-X(3))+X(20)*(X(9)-X(24))
  43. *
  44.       B(1,5)= X(23)*(X(12)-X(21)+X(3)-X(18))+X(20)*(X(24)-X(18))+
  45.      & X(17)*(X(21)-X(6)+X(24)-X(3))+X(2)*(X(6)-X(12)+X(18)-X(24))+
  46.      & X(11)*(X(3)-X(24))+X(5)*(X(18)-X(3))
  47. *
  48.       B(1,6)= X(14)*(X(3)-X(24)+X(6)-X(21))+X(23)*(X(15)-X(21))+
  49.      & X(20)*(X(24)-X(9)+X(15)-X(6))+X(5)*(X(9)-X(3)+X(21)-X(15))+
  50.      & X(2)*(X(6)-X(15))+X(8)*(X(21)-X(6))
  51. *
  52.       B(1,7)= X(17)*(X(6)-X(15)+X(9)-X(24))+X(14)*(X(18)-X(24))+
  53.      & X(23)*(X(15)-X(12)+X(18)-X(9))+X(8)*(X(12)-X(6)+X(24)-X(18))+
  54.      & X(5)*(X(9)-X(18))+X(11)*(X(24)-X(9))
  55. *
  56.       B(1,8)= X(20)*(X(9)-X(18)+X(12)-X(15))+X(17)*(X(21)-X(15))+
  57.      & X(14)*(X(18)-X(3)+X(21)-X(12))+X(11)*(X(3)-X(9)+X(15)-X(21))+
  58.      & X(8)*(X(12)-X(21))+X(2)*(X(15)-X(12))
  59. *
  60. *
  61. *
  62.       B(2,1)= X(6)*(X(16)-X(7)+X(13)-X(10))+X(9)*(X(4)-X(10))+
  63.      & X(12)*(X(7)-X(22)+X(4)-X(13))+X(15)*(X(22)-X(16)+X(10)-X(4))+
  64.      & X(18)*(X(13)-X(4))+X(24)*(X(10)-X(13))
  65. *
  66.       B(2,2)= X(9)*(X(19)-X(10)+X(16)-X(1))+X(12)*(X(7)-X(1))+
  67.      & X(3)*(X(10)-X(13)+X(7)-X(16))+X(18)*(X(13)-X(19)+X(1)-X(7))+
  68.      & X(21)*(X(16)-X(7))+X(15)*(X(1)-X(16))
  69. *
  70.       B(2,3)= X(12)*(X(22)-X(1)+X(19)-X(4))+X(3)*(X(10)-X(4))+
  71.      & X(6)*(X(1)-X(16)+X(10)-X(19))+X(21)*(X(16)-X(22)+X(4)-X(10))+
  72.      & X(24)*(X(19)-X(10))+X(18)*(X(4)-X(19))
  73. *
  74.       B(2,4)= X(3)*(X(13)-X(4)+X(22)-X(7))+X(6)*(X(1)-X(7))+
  75.      & X(9)*(X(4)-X(19)+X(1)-X(22))+X(24)*(X(19)-X(13)+X(7)-X(1))+
  76.      & X(15)*(X(22)-X(1))+X(21)*(X(7)-X(22))
  77. *
  78.       B(2,5)= X(24)*(X(10)-X(19)+X(1)-X(16))+X(21)*(X(22)-X(16))+
  79.      & X(18)*(X(19)-X(4)+X(22)-X(1))+X(3)*(X(4)-X(10)+X(16)-X(22))+
  80.      & X(12)*(X(1)-X(22))+X(6)*(X(16)-X(1))
  81. *
  82.       B(2,6)= X(15)*(X(1)-X(22)+X(4)-X(19))+X(24)*(X(13)-X(19))+
  83.      & X(21)*(X(22)-X(7)+X(13)-X(4))+X(6)*(X(7)-X(1)+X(19)-X(13))+
  84.      & X(3)*(X(4)-X(13))+X(9)*(X(19)-X(4))
  85. *
  86.       B(2,7)= X(18)*(X(4)-X(13)+X(7)-X(22))+X(15)*(X(16)-X(22))+
  87.      & X(24)*(X(13)-X(10)+X(16)-X(7))+X(9)*(X(10)-X(4)+X(22)-X(16))+
  88.      & X(6)*(X(7)-X(16))+X(12)*(X(22)-X(7))
  89. *
  90.       B(2,8)= X(21)*(X(7)-X(16)+X(10)-X(13))+X(18)*(X(19)-X(13))+
  91.      & X(15)*(X(16)-X(1)+X(19)-X(10))+X(12)*(X(1)-X(7)+X(13)-X(19))+
  92.      & X(9)*(X(10)-X(19))+X(3)*(X(13)-X(10))
  93. *
  94. *
  95. *
  96.       B(3,1)= X(4)*(X(17)-X(8)+X(14)-X(11))+X(7)*(X(5)-X(11))+
  97.      & X(10)*(X(8)-X(23)+X(5)-X(14))+X(13)*(X(23)-X(17)+X(11)-X(5))+
  98.      & X(16)*(X(14)-X(5))+X(22)*(X(11)-X(14))
  99. *
  100.       B(3,2)= X(7)*(X(20)-X(11)+X(17)-X(2))+X(10)*(X(8)-X(2))+
  101.      & X(1)*(X(11)-X(14)+X(8)-X(17))+X(16)*(X(14)-X(20)+X(2)-X(8))+
  102.      & X(19)*(X(17)-X(8))+X(13)*(X(2)-X(17))
  103. *
  104.       B(3,3)= X(10)*(X(23)-X(2)+X(20)-X(5))+X(1)*(X(11)-X(5))+
  105.      & X(4)*(X(2)-X(17)+X(11)-X(20))+X(19)*(X(17)-X(23)+X(5)-X(11))+
  106.      & X(22)*(X(20)-X(11))+X(16)*(X(5)-X(20))
  107. *
  108.       B(3,4)= X(1)*(X(14)-X(5)+X(23)-X(8))+X(4)*(X(2)-X(8))+
  109.      & X(7)*(X(5)-X(20)+X(2)-X(23))+X(22)*(X(20)-X(14)+X(8)-X(2))+
  110.      & X(13)*(X(23)-X(2))+X(19)*(X(8)-X(23))
  111. *
  112.       B(3,5)= X(22)*(X(11)-X(20)+X(2)-X(17))+X(19)*(X(23)-X(17))+
  113.      & X(16)*(X(20)-X(5)+X(23)-X(2))+X(1)*(X(5)-X(11)+X(17)-X(23))+
  114.      & X(10)*(X(2)-X(23))+X(4)*(X(17)-X(2))
  115. *
  116.       B(3,6)= X(13)*(X(2)-X(23)+X(5)-X(20))+X(22)*(X(14)-X(20))+
  117.      & X(19)*(X(23)-X(8)+X(14)-X(5))+X(4)*(X(8)-X(2)+X(20)-X(14))+
  118.      & X(1)*(X(5)-X(14))+X(7)*(X(20)-X(5))
  119. *
  120.       B(3,7)= X(16)*(X(5)-X(14)+X(8)-X(23))+X(13)*(X(17)-X(23))+
  121.      & X(22)*(X(14)-X(11)+X(17)-X(8))+X(7)*(X(11)-X(5)+X(23)-X(17))+
  122.      & X(4)*(X(8)-X(17))+X(10)*(X(23)-X(8))
  123. *
  124.       B(3,8)= X(19)*(X(8)-X(17)+X(11)-X(14))+X(16)*(X(20)-X(14))+
  125.      & X(13)*(X(17)-X(2)+X(20)-X(11))+X(10)*(X(2)-X(8)+X(14)-X(20))+
  126.      & X(7)*(X(11)-X(20))+X(1)*(X(14)-X(11))
  127. *
  128. *
  129. ******* CALCUL DU VOLUME DE L'ELEMENT
  130. *
  131.       VOL= UNS12*(B(1,1)*X(1)+B(1,2)*X(4)+B(1,3)*X(7)+B(1,4)*X(10)+
  132.      &  B(1,5)*X(13)+B(1,6)*X(16)+B(1,7)*X(19)+B(1,8)*X(22))
  133. *
  134. ******* CALCUL DES MATRICES B
  135. *
  136.       DO J=1,3
  137.          DO I=1,8
  138.               B(J,I)= UNS12*B(J,I)/VOL
  139.          END DO
  140.       END DO
  141. *
  142. *
  143.       RETURN
  144.       END
  145.  
  146.  
  147.  

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