Télécharger intdel.eso

Retour à la liste

Numérotation des lignes :

  1. C INTDEL SOURCE AM 15/08/21 21:15:29 8599
  2. SUBROUTINE INTDEL(XNTH,XNTB,XNTT,LRN,MELE)
  3. C
  4. C------------------------------------------------------------------------
  5. C
  6. C CETTE ROUTINE CONTIENT LES INTEGRALS QUI PRESENTENT DES
  7. C FONCTIONS DELTA DE DIRAC DANS LE CAS DES ELEMENTS JOINTS
  8. C POREUX (BALD)
  9. C CETTE ROUTINE EST APPELEE PAR: PERMAB.ESO, RIGI2.ESO,
  10. C EPSI2.ESO, SIGPOR.ESO.
  11. C
  12. C------------------------------------------------------------------------
  13. C
  14. C ENTREES
  15. C MELE = NUMERO DE L ELEMENT FINI DANS NOMTP
  16. C LRN = DIMENSIONS DES TABLEAUX
  17. C
  18. C SORTIES
  19. C XNTH = TABLEAU QUI SERT POUR LA MATRICE DE PERMEABILITE
  20. C XNTB = TABLEAU QUI SERT POUR LA MATRICE DE PERMEABILITE
  21. C XNTT = TABLEAU QUI SERT POUR LES MATRICES DE PERMEABILITE
  22. C ET DE RIGIDITE
  23. C
  24. C AUTEUR: BALD 20/01/96
  25. C
  26. C------------------------------------------------------------------------
  27. C
  28. IMPLICIT INTEGER(I-N)
  29. IMPLICIT REAL*8(A-H,O-Z)
  30. C
  31. -INC CCOPTIO
  32. C
  33. DIMENSION XNTH(LRN,LRN),XNTB(LRN,LRN),XNTT(LRN)
  34. PARAMETER(XZER=0.D0,UN=1.D0)
  35. C
  36. CALL ZERO(XNTH,LRN,LRN)
  37. CALL ZERO(XNTB,LRN,LRN)
  38. CALL ZERO(XNTT,1,LRN)
  39. C
  40. IF(MELE.EQ.108.OR.MELE.EQ.185.OR.MELE.EQ.188)THEN
  41. *
  42. * CAS DES ELEMENTS JOP3, JOQ3, JOR3
  43. *
  44. XNTH(1,1)=UN
  45. XNTH(1,2)=UN
  46. XNTH(1,3)=XZER
  47. XNTH(1,4)=XZER
  48. XNTH(1,5)=-UN
  49. XNTH(1,6)=-UN
  50. *
  51. XNTH(2,1)=UN
  52. XNTH(2,2)=UN
  53. XNTH(2,3)=XZER
  54. XNTH(2,4)=XZER
  55. XNTH(2,5)=-UN
  56. XNTH(2,6)=-UN
  57. *
  58. XNTH(5,1)=-UN
  59. XNTH(5,2)=-UN
  60. XNTH(5,3)=XZER
  61. XNTH(5,4)=XZER
  62. XNTH(5,5)=UN
  63. XNTH(5,6)=UN
  64. *
  65. XNTH(6,1)=-UN
  66. XNTH(6,2)=-UN
  67. XNTH(6,3)=XZER
  68. XNTH(6,4)=XZER
  69. XNTH(6,5)=UN
  70. XNTH(6,6)=UN
  71. ***
  72. XNTB(3,1)=XZER
  73. XNTB(3,2)=XZER
  74. XNTB(3,3)=UN
  75. XNTB(3,4)=UN
  76. XNTB(3,5)=-UN
  77. XNTB(3,6)=-UN
  78. *
  79. XNTB(4,1)=XZER
  80. XNTB(4,2)=XZER
  81. XNTB(4,3)=UN
  82. XNTB(4,4)=UN
  83. XNTB(4,5)=-UN
  84. XNTB(4,6)=-UN
  85. *
  86. XNTB(5,1)=XZER
  87. XNTB(5,2)=XZER
  88. XNTB(5,3)=-UN
  89. XNTB(5,4)=-UN
  90. XNTB(5,5)=UN
  91. XNTB(5,6)=UN
  92. *
  93. XNTB(6,1)=XZER
  94. XNTB(6,2)=XZER
  95. XNTB(6,3)=-UN
  96. XNTB(6,4)=-UN
  97. XNTB(6,5)=UN
  98. XNTB(6,6)=UN
  99. ***
  100. XNTT(1)=XZER
  101. XNTT(2)=XZER
  102. XNTT(3)=XZER
  103. XNTT(4)=XZER
  104. XNTT(5)=UN
  105. XNTT(6)=UN
  106. ELSE IF(MELE.EQ.109.OR.MELE.EQ.186.OR.MELE.EQ.189) THEN
  107. *
  108. * CAS DES ELEMENTS JOP6, JOQ6, JOR6
  109. *
  110. XNTH(1,1)=UN
  111. XNTH(1,2)=UN
  112. XNTH(1,3)=UN
  113. XNTH(1,4)=XZER
  114. XNTH(1,5)=XZER
  115. XNTH(1,6)=XZER
  116. XNTH(1,7)=-UN
  117. XNTH(1,8)=-UN
  118. XNTH(1,9)=-UN
  119. *
  120. XNTH(7,1)=-UN
  121. XNTH(7,2)=-UN
  122. XNTH(7,3)=-UN
  123. XNTH(7,4)=XZER
  124. XNTH(7,5)=XZER
  125. XNTH(7,6)=XZER
  126. XNTH(7,7)=UN
  127. XNTH(7,8)=UN
  128. XNTH(7,9)=UN
  129. *
  130. DO 1 M=1,7,3
  131. DO 2 K=1,2
  132. DO 3 L=1,9
  133. XNTH(K+M,L)=XNTH(M,L)
  134. 3 CONTINUE
  135. 2 CONTINUE
  136. 1 CONTINUE
  137. ***
  138. XNTB(4,1)=XZER
  139. XNTB(4,2)=XZER
  140. XNTB(4,3)=XZER
  141. XNTB(4,4)=UN
  142. XNTB(4,5)=UN
  143. XNTB(4,6)=UN
  144. XNTB(4,7)=-UN
  145. XNTB(4,8)=-UN
  146. XNTB(4,9)=-UN
  147. *
  148. XNTB(7,1)=XZER
  149. XNTB(7,2)=XZER
  150. XNTB(7,3)=XZER
  151. XNTB(7,4)=-UN
  152. XNTB(7,5)=-UN
  153. XNTB(7,6)=-UN
  154. XNTB(7,7)=UN
  155. XNTB(7,8)=UN
  156. XNTB(7,9)=UN
  157. *
  158. DO 7 M=1,7,3
  159. DO 8 K=1,2
  160. DO 9 L=1,9
  161. XNTB(K+M,L)=XNTB(M,L)
  162. 9 CONTINUE
  163. 8 CONTINUE
  164. 7 CONTINUE
  165. ***
  166. XNTT(1)=XZER
  167. XNTT(2)=XZER
  168. XNTT(3)=XZER
  169. XNTT(4)=XZER
  170. XNTT(5)=XZER
  171. XNTT(6)=XZER
  172. XNTT(7)=UN
  173. XNTT(8)=UN
  174. XNTT(9)=UN
  175. ELSE IF(MELE.EQ.110.OR.MELE.EQ.187.OR.MELE.EQ.190)THEN
  176. *
  177. * CAS DES ELEMENTS JOP8, JOQ8, JOR8
  178. *
  179. XNTH(1,1) =UN
  180. XNTH(1,2) =UN
  181. XNTH(1,3) =UN
  182. XNTH(1,4) =UN
  183. XNTH(1,5) =XZER
  184. XNTH(1,6) =XZER
  185. XNTH(1,7) =XZER
  186. XNTH(1,8) =XZER
  187. XNTH(1,9) =-UN
  188. XNTH(1,10)=-UN
  189. XNTH(1,11)=-UN
  190. XNTH(1,12)=-UN
  191. *
  192. XNTH(9,1) =-UN
  193. XNTH(9,2) =-UN
  194. XNTH(9,3) =-UN
  195. XNTH(9,4) =-UN
  196. XNTH(9,5) =XZER
  197. XNTH(9,6) =XZER
  198. XNTH(9,7) =XZER
  199. XNTH(9,8) =XZER
  200. XNTH(9,9) =UN
  201. XNTH(9,10)=UN
  202. XNTH(9,11)=UN
  203. XNTH(9,12)=UN
  204. *
  205. DO 4 M=1,9,4
  206. DO 5 K=1,3
  207. DO 6 L=1,12
  208. XNTH(K+M,L)=XNTH(M,L)
  209. 6 CONTINUE
  210. 5 CONTINUE
  211. 4 CONTINUE
  212. ***
  213. XNTB(5,1) =XZER
  214. XNTB(5,2) =XZER
  215. XNTB(5,3) =XZER
  216. XNTB(5,4) =XZER
  217. XNTB(5,5) =UN
  218. XNTB(5,6) =UN
  219. XNTB(5,7) =UN
  220. XNTB(5,8) =UN
  221. XNTB(5,9) =-UN
  222. XNTB(5,10)=-UN
  223. XNTB(5,11)=-UN
  224. XNTB(5,12)=-UN
  225. *
  226. XNTB(9,1) =XZER
  227. XNTB(9,2) =XZER
  228. XNTB(9,3) =XZER
  229. XNTB(9,4) =XZER
  230. XNTB(9,5) =-UN
  231. XNTB(9,6) =-UN
  232. XNTB(9,7) =-UN
  233. XNTB(9,8) =-UN
  234. XNTB(9,9) =UN
  235. XNTB(9,10)=UN
  236. XNTB(9,11)=UN
  237. XNTB(9,12)=UN
  238. *
  239. DO 10 M=1,9,4
  240. DO 11 K=1,3
  241. DO 12 L=1,12
  242. XNTB(K+M,L)=XNTB(M,L)
  243. 12 CONTINUE
  244. 11 CONTINUE
  245. 10 CONTINUE
  246. ***
  247. XNTT(1) =XZER
  248. XNTT(2) =XZER
  249. XNTT(3) =XZER
  250. XNTT(4) =XZER
  251. XNTT(5) =XZER
  252. XNTT(6) =XZER
  253. XNTT(7) =XZER
  254. XNTT(8) =XZER
  255. XNTT(9) =UN
  256. XNTT(10)=UN
  257. XNTT(11)=UN
  258. XNTT(12)=UN
  259. ELSE
  260. CALL ERREUR(5)
  261. END IF
  262. C
  263. RETURN
  264. END
  265.  
  266.  
  267.  
  268.  
  269.  

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