Télécharger la plaquette Cast3M


fr en

Télécharger intdel.eso

Retour à la liste

Numérotation des lignes : 

intdel

  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