Esope source files

   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 :