Télécharger la plaquette Cast3M


fr en

Télécharger matp1.eso

Retour à la liste

Numérotation des lignes : 

matp1

  1. C MATP1     SOURCE    CB215821  25/04/23    21:15:29     12247          
  2.       SUBROUTINE MATP1(NBMAIL,IPGEOM,IPMAHY,IPRIG1,COEF,DELTAT,IPCK,
  3.      S                 IPRIG2)
  4. C-----------------------------------------------------------------------
  5. C Calcul de la matrice en trace de charge dans le cas d'une
  6. C formulation mixte hybride pour les equations de DARCY.
  7. C-----------------------------------------------------------------------
  8. C
  9. C---------------------------
  10. C Parametres Entree/Sortie :
  11. C---------------------------
  12. C
  13. C E/  NBMAIL  : Nombre de zones élémentaires de MMODEL
  14. C E/  IPGEOM  : Pointeur de l'objet maillage
  15. C E/  IPMAHY  : Segment contenant le pointeur vers le meleme des
  16. C               connectivites elements/faces pour les zones du MMODEL
  17. C               ou on a defini DARCY.
  18. C E/  IPRIG1  : RIGIDITE de sous type 'DARCY'
  19. C E/  COEF    : Parametre de discretisation temporelle (theta-methode)
  20. C E/  DELTAT  : Pas de discretisation temporelle
  21. C E/  IPCK    : MCHAML donnant pour chaque element Ck|K|
  22. C  /S IPRIG2  : RIGIDITE de sous type 'HYBTP'
  23. C
  24. C----------------------
  25. C Variables en COMMON :
  26. C----------------------
  27. C
  28. C E/  IFOMOD    : Voir CCOPTIO
  29. C E/  NOMDD(20) : Voir CCHAMP
  30. C E/  NOMDU(20) : Voir CCHAMP
  31. C
  32. C
  33. C-----------------------------------------------------------------------
  34. C
  35. C Langage : ESOPE + FORTRAN77
  36. C
  37. C Auteurs : 08/93 F.DABBENE - Cas permanent
  38. C           09/94 X.NOUVELLON - Extension au cas transitoire
  39. C
  40. C-----------------------------------------------------------------------
  41.       IMPLICIT INTEGER(I-N)
  42.       IMPLICIT REAL*8 (A-H,O-Z)
  43. *
  44.  
  45. -INC PPARAM
  46. -INC CCOPTIO
  47. -INC CCHAMP
  48. -INC SMRIGID
  49. -INC SMCOORD
  50. -INC SMCHAML
  51. -INC SMELEME
  52. *
  53.       SEGMENT IPMAHY
  54.          INTEGER MAHYBR(NSOUS)
  55.       ENDSEGMENT
  56. *
  57. * Initialisations
  58. *
  59.       RI1    = IPRIG1
  60.       MCHELM = IPCK
  61. *
  62. * Création du chapeau de l'objet RIGIDITE
  63. *
  64.       NRIGE  = 6
  65.       NRIGEL = NBMAIL
  66.       SEGINI MRIGID
  67.       IPRIG2 = MRIGID
  68.       ICHOLE = 0
  69.       IMGEO1 = 0
  70.       IMGEO2 = 0
  71.       IFORIG = IFOUR
  72.       ISUPEQ = 0
  73.       MTYMAT = 'HYBTP   '
  74.       NBGEOR = 0
  75.       DO 10 IA=1,NBMAIL
  76.          IRIGEL(4,IA) = 0
  77.          COERIG(IA)   = 1.D0
  78.  10   CONTINUE
  79. *
  80. *= BOUCLE SUR LES MAILLAGES ELEMENTAIRES,ZONE IMAIL
  81. *
  82.       SEGACT IPMAHY
  83.       IF (IPCK.NE.0) SEGACT MCHELM
  84.       SEGACT RI1
  85.       DO 30 IMAIL=1,NBMAIL
  86. *
  87. * Recuperation de l'objet maillage ELTFA pour la zone IMAIL
  88. *
  89.          IMAHYB = MAHYBR(IMAIL)
  90.          IF (IMAHYB.EQ.0) GOTO 30
  91. *
  92. * Récupération du MELVAL de la sous zone pour le MCHAML
  93. *
  94.          IF (IPCK.NE.0) THEN
  95.             MCHAML = ICHAML(IMAIL)
  96.             SEGACT MCHAML
  97.             MELVAL = IELVAL(1)
  98.          ELSE
  99.             MELVAL = 0
  100.          ENDIF
  101. *
  102. * Récupération des matrices masses hybrides pour la sous zone
  103. *
  104.          xMATR1 = RI1.IRIGEL(4,IMAIL)
  105. *
  106. * Création du segment DESCRIPTEUR pour la sous zone IMAIL
  107. *
  108.          MELEME = IMAHYB
  109.          SEGACT MELEME
  110.          NBDDL  = NUM(/1)
  111.          NBELEM = NUM(/2)
  112. *
  113.          NLIGRP = NBDDL
  114.          NLIGRD = NBDDL
  115.          SEGINI DESCR
  116.          DO 20 IB=1,NLIGRP
  117.             LISINC(IB) = NOMDD(20)
  118.             LISDUA(IB) = NOMDU(20)
  119.             NOELEP(IB) = IB
  120.             NOELED(IB) = IB
  121.  20      CONTINUE
  122.          IDESCR = DESCR
  123.          SEGDES DESCR
  124. *
  125. * Création du segment IMATRI
  126. * Initialisation du tableau IRIGEL pour la sous zone IMAIL
  127. *
  128.          NELRIG = NBELEM
  129.          SEGINI,xMATRI=xmatr1
  130.          IRIGEL(1,IMAIL) = IMAHYB
  131.          IRIGEL(2,IMAIL) = 0
  132.          IRIGEL(3,IMAIL) = IDESCR
  133.          IRIGEL(4,IMAIL) = xMATRI
  134.          IRIGEL(5,IMAIL) = 0
  135.          IRIGEL(6,IMAIL) = 0
  136.          SEGDES xMATRI
  137. *
  138. * Calcul de RIGIDITE resultat pour la sous zone
  139. *
  140.          CALL MATP2(xmatr1,COEF,DELTAT,MELVAL,xMATRI)
  141.  30   CONTINUE
  142.       SEGDES RI1
  143.       SEGDES IPMAHY
  144.       SEGDES MRIGID
  145.       END
  146.  
  147.  
  148.  
  149.  
  150.  

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