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.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34. C
  35. DIMENSION XNTH(LRN,LRN),XNTB(LRN,LRN),XNTT(LRN)
  36. PARAMETER(XZER=0.D0,UN=1.D0)
  37. C
  38. CALL ZERO(XNTH,LRN,LRN)
  39. CALL ZERO(XNTB,LRN,LRN)
  40. CALL ZERO(XNTT,1,LRN)
  41. C
  42. IF(MELE.EQ.108.OR.MELE.EQ.185.OR.MELE.EQ.188)THEN
  43. *
  44. * CAS DES ELEMENTS JOP3, JOQ3, JOR3
  45. *
  46. XNTH(1,1)=UN
  47. XNTH(1,2)=UN
  48. XNTH(1,3)=XZER
  49. XNTH(1,4)=XZER
  50. XNTH(1,5)=-UN
  51. XNTH(1,6)=-UN
  52. *
  53. XNTH(2,1)=UN
  54. XNTH(2,2)=UN
  55. XNTH(2,3)=XZER
  56. XNTH(2,4)=XZER
  57. XNTH(2,5)=-UN
  58. XNTH(2,6)=-UN
  59. *
  60. XNTH(5,1)=-UN
  61. XNTH(5,2)=-UN
  62. XNTH(5,3)=XZER
  63. XNTH(5,4)=XZER
  64. XNTH(5,5)=UN
  65. XNTH(5,6)=UN
  66. *
  67. XNTH(6,1)=-UN
  68. XNTH(6,2)=-UN
  69. XNTH(6,3)=XZER
  70. XNTH(6,4)=XZER
  71. XNTH(6,5)=UN
  72. XNTH(6,6)=UN
  73. ***
  74. XNTB(3,1)=XZER
  75. XNTB(3,2)=XZER
  76. XNTB(3,3)=UN
  77. XNTB(3,4)=UN
  78. XNTB(3,5)=-UN
  79. XNTB(3,6)=-UN
  80. *
  81. XNTB(4,1)=XZER
  82. XNTB(4,2)=XZER
  83. XNTB(4,3)=UN
  84. XNTB(4,4)=UN
  85. XNTB(4,5)=-UN
  86. XNTB(4,6)=-UN
  87. *
  88. XNTB(5,1)=XZER
  89. XNTB(5,2)=XZER
  90. XNTB(5,3)=-UN
  91. XNTB(5,4)=-UN
  92. XNTB(5,5)=UN
  93. XNTB(5,6)=UN
  94. *
  95. XNTB(6,1)=XZER
  96. XNTB(6,2)=XZER
  97. XNTB(6,3)=-UN
  98. XNTB(6,4)=-UN
  99. XNTB(6,5)=UN
  100. XNTB(6,6)=UN
  101. ***
  102. XNTT(1)=XZER
  103. XNTT(2)=XZER
  104. XNTT(3)=XZER
  105. XNTT(4)=XZER
  106. XNTT(5)=UN
  107. XNTT(6)=UN
  108. ELSE IF(MELE.EQ.109.OR.MELE.EQ.186.OR.MELE.EQ.189) THEN
  109. *
  110. * CAS DES ELEMENTS JOP6, JOQ6, JOR6
  111. *
  112. XNTH(1,1)=UN
  113. XNTH(1,2)=UN
  114. XNTH(1,3)=UN
  115. XNTH(1,4)=XZER
  116. XNTH(1,5)=XZER
  117. XNTH(1,6)=XZER
  118. XNTH(1,7)=-UN
  119. XNTH(1,8)=-UN
  120. XNTH(1,9)=-UN
  121. *
  122. XNTH(7,1)=-UN
  123. XNTH(7,2)=-UN
  124. XNTH(7,3)=-UN
  125. XNTH(7,4)=XZER
  126. XNTH(7,5)=XZER
  127. XNTH(7,6)=XZER
  128. XNTH(7,7)=UN
  129. XNTH(7,8)=UN
  130. XNTH(7,9)=UN
  131. *
  132. DO 1 M=1,7,3
  133. DO 2 K=1,2
  134. DO 3 L=1,9
  135. XNTH(K+M,L)=XNTH(M,L)
  136. 3 CONTINUE
  137. 2 CONTINUE
  138. 1 CONTINUE
  139. ***
  140. XNTB(4,1)=XZER
  141. XNTB(4,2)=XZER
  142. XNTB(4,3)=XZER
  143. XNTB(4,4)=UN
  144. XNTB(4,5)=UN
  145. XNTB(4,6)=UN
  146. XNTB(4,7)=-UN
  147. XNTB(4,8)=-UN
  148. XNTB(4,9)=-UN
  149. *
  150. XNTB(7,1)=XZER
  151. XNTB(7,2)=XZER
  152. XNTB(7,3)=XZER
  153. XNTB(7,4)=-UN
  154. XNTB(7,5)=-UN
  155. XNTB(7,6)=-UN
  156. XNTB(7,7)=UN
  157. XNTB(7,8)=UN
  158. XNTB(7,9)=UN
  159. *
  160. DO 7 M=1,7,3
  161. DO 8 K=1,2
  162. DO 9 L=1,9
  163. XNTB(K+M,L)=XNTB(M,L)
  164. 9 CONTINUE
  165. 8 CONTINUE
  166. 7 CONTINUE
  167. ***
  168. XNTT(1)=XZER
  169. XNTT(2)=XZER
  170. XNTT(3)=XZER
  171. XNTT(4)=XZER
  172. XNTT(5)=XZER
  173. XNTT(6)=XZER
  174. XNTT(7)=UN
  175. XNTT(8)=UN
  176. XNTT(9)=UN
  177. ELSE IF(MELE.EQ.110.OR.MELE.EQ.187.OR.MELE.EQ.190)THEN
  178. *
  179. * CAS DES ELEMENTS JOP8, JOQ8, JOR8
  180. *
  181. XNTH(1,1) =UN
  182. XNTH(1,2) =UN
  183. XNTH(1,3) =UN
  184. XNTH(1,4) =UN
  185. XNTH(1,5) =XZER
  186. XNTH(1,6) =XZER
  187. XNTH(1,7) =XZER
  188. XNTH(1,8) =XZER
  189. XNTH(1,9) =-UN
  190. XNTH(1,10)=-UN
  191. XNTH(1,11)=-UN
  192. XNTH(1,12)=-UN
  193. *
  194. XNTH(9,1) =-UN
  195. XNTH(9,2) =-UN
  196. XNTH(9,3) =-UN
  197. XNTH(9,4) =-UN
  198. XNTH(9,5) =XZER
  199. XNTH(9,6) =XZER
  200. XNTH(9,7) =XZER
  201. XNTH(9,8) =XZER
  202. XNTH(9,9) =UN
  203. XNTH(9,10)=UN
  204. XNTH(9,11)=UN
  205. XNTH(9,12)=UN
  206. *
  207. DO 4 M=1,9,4
  208. DO 5 K=1,3
  209. DO 6 L=1,12
  210. XNTH(K+M,L)=XNTH(M,L)
  211. 6 CONTINUE
  212. 5 CONTINUE
  213. 4 CONTINUE
  214. ***
  215. XNTB(5,1) =XZER
  216. XNTB(5,2) =XZER
  217. XNTB(5,3) =XZER
  218. XNTB(5,4) =XZER
  219. XNTB(5,5) =UN
  220. XNTB(5,6) =UN
  221. XNTB(5,7) =UN
  222. XNTB(5,8) =UN
  223. XNTB(5,9) =-UN
  224. XNTB(5,10)=-UN
  225. XNTB(5,11)=-UN
  226. XNTB(5,12)=-UN
  227. *
  228. XNTB(9,1) =XZER
  229. XNTB(9,2) =XZER
  230. XNTB(9,3) =XZER
  231. XNTB(9,4) =XZER
  232. XNTB(9,5) =-UN
  233. XNTB(9,6) =-UN
  234. XNTB(9,7) =-UN
  235. XNTB(9,8) =-UN
  236. XNTB(9,9) =UN
  237. XNTB(9,10)=UN
  238. XNTB(9,11)=UN
  239. XNTB(9,12)=UN
  240. *
  241. DO 10 M=1,9,4
  242. DO 11 K=1,3
  243. DO 12 L=1,12
  244. XNTB(K+M,L)=XNTB(M,L)
  245. 12 CONTINUE
  246. 11 CONTINUE
  247. 10 CONTINUE
  248. ***
  249. XNTT(1) =XZER
  250. XNTT(2) =XZER
  251. XNTT(3) =XZER
  252. XNTT(4) =XZER
  253. XNTT(5) =XZER
  254. XNTT(6) =XZER
  255. XNTT(7) =XZER
  256. XNTT(8) =XZER
  257. XNTT(9) =UN
  258. XNTT(10)=UN
  259. XNTT(11)=UN
  260. XNTT(12)=UN
  261. ELSE
  262. CALL ERREUR(5)
  263. END IF
  264. C
  265. RETURN
  266. END
  267.  
  268.  
  269.  
  270.  
  271.  

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