Télécharger la plaquette Cast3M


fr en

Télécharger ichaxr.eso

Retour à la liste

Numérotation des lignes : 

ichaxr

  1. C ICHAXR    SOURCE    CB215821  25/04/23    21:15:23     12247          
  2.  
  3.       SUBROUTINE ICHAXR (IPMODL,ICHELX,MCHELX)
  4.  
  5. C-----------------------------------------------------------------------
  6. C  SOUS-PROGRAMME DE RECUPERATION DU MCHAML D'ENRICHISSEMENT CONTENU
  7. C  DANS LE MODELE ASSOCIE A DES ELEMENTS DE TYPE XFEM/RUPTURE
  8. C-----------------------------------------------------------------------
  9. C  IPMODL (E)  Modele elementaire de REFERENCE (complet) s'appyant sur
  10. C              des elements de type XFEM (formulation MFR=63), contenant
  11. C              en particulier le MCHAML d'ENRICHISSEMENT
  12. C              Segment IMODEL : ACTIF en Entree/Sortie & non modifie
  13. C  ICHELX (S)  Indice dans le tableau IVAMOD du modele IPMODL du MCHAML
  14. C              d'ENRICHISSEMENT, =0 si MCHAML n'existe pas
  15. C  MCHELX (S)  MCHAML d'ENRICHISSEMENT contenu dans le modele
  16. C              = 0 s'il n'est pas present dans le modele
  17. C              Segment MCHELM : ACTIF en Sortie  (si present)
  18. C-----------------------------------------------------------------------
  19.  
  20.       IMPLICIT INTEGER (I-N)
  21.       IMPLICIT REAL*8 (A-H,O-Z)
  22.  
  23.  
  24. -INC PPARAM
  25. -INC CCOPTIO
  26.  
  27. -INC SMCHAML
  28. -INC SMCOORD
  29. -INC SMMODEL
  30.  
  31.       POINTEUR MCHEX1.MCHELM
  32.  
  33.       IMODEL = IPMODL
  34.       ICHELX = 0
  35.       MCHELX = 0
  36.  
  37. C*      SEGACT,IMODEL
  38. C*      IF (IPMODL.EQ.0) THEN
  39. C*        WRITE(IOIMP,*) '#####  SOUS-PROGRAMME ICHAXR #####'
  40. C*        CALL ERREUR(21)
  41. C*        RETURN
  42. C*      ENDIF
  43.  
  44. C- Recuperation du MCHAML d enrichissement dans IMODEL
  45.       NOBMOD = IMODEL.IVAMOD(/1)
  46.       DO i = 1, NOBMOD
  47.         IF (IMODEL.TYMODE(i) .EQ. 'MCHAML') THEN
  48.           MCHELM = IMODEL.IVAMOD(i)
  49.           SEGACT,MCHELM
  50.           IF (TITCHE .EQ. 'ENRICHIS') THEN
  51.             ICHELX = i
  52.             MCHELX = MCHELM
  53.             GOTO 10
  54.           ENDIF
  55.           SEGDES,MCHELM
  56.         ENDIF
  57.         IF (IMODEL.TYMODE(i) .EQ. 'MCHAMLL') THEN
  58.  
  59.           MCHAML = IMODEL.IVAMOD(i)
  60.  
  61.           SEGACT,MCHAML
  62. C++++  INITIALISATION d'un nouveau chamelem chapeau vierge
  63.         L1 = 8
  64. C nombre de sous champs
  65.         N1 = 1
  66. C nombre de composantes
  67.         N2 = 3  
  68.         N3 = 6
  69.  
  70.         SEGINI, MCHEX1
  71.         MCHEX1.TITCHE = 'ENRICHIS'
  72.         MCHEX1.IFOCHE = IFOUR
  73.         MCHEX1.IMACHE(1)= IMAMOD
  74.         MCHEX1.CONCHE(1)= IMODEL.CONMOD
  75.         MCHEX1.INFCHE(1,2) = 0
  76.         MCHEX1.INFCHE(1,3) = NIFOUR
  77.         MCHEX1.INFCHE(1,6) = 1
  78.         ICHELX = i
  79.  
  80.  
  81.         MCHEX1.ICHAML(1)=MCHAML
  82.  
  83.         MCHELX = MCHEX1
  84.         ICHELX = i
  85.             GOTO 10
  86.         ENDIF
  87.       ENDDO
  88. c      IF (ICHELX .EQ. 0) THEN
  89. c        WRITE(IOIMP,*) 'Il n y a pas de MCHAML d ENRICHISsement ',
  90. c     &                 'dans le modele'
  91. c        RETURN
  92. c      ENDIF
  93.  
  94.  10   CONTINUE
  95. C*      SEGDES,IMODEL
  96.  
  97.       RETURN
  98.       END
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  

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