Télécharger cubebb.eso

Retour à la liste

Numérotation des lignes :

  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