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  24/04/12    21:16:18     11897          
  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 SMMODEL
  29.  
  30.       POINTEUR MCHEX1.MCHELM
  31.  
  32.       IMODEL = IPMODL
  33.       ICHELX = 0
  34.       MCHELX = 0
  35.  
  36. C*      SEGACT,IMODEL
  37. C*      IF (IPMODL.EQ.0) THEN
  38. C*        WRITE(IOIMP,*) '#####  SOUS-PROGRAMME ICHAXR #####'
  39. C*        CALL ERREUR(21)
  40. C*        RETURN
  41. C*      ENDIF
  42.  
  43. C- Recuperation du MCHAML d enrichissement dans IMODEL
  44.       NOBMOD = IMODEL.IVAMOD(/1)
  45.       DO i = 1, NOBMOD
  46.         IF (IMODEL.TYMODE(i) .EQ. 'MCHAML') THEN
  47.           MCHELM = IMODEL.IVAMOD(i)
  48.           SEGACT,MCHELM
  49.           IF (TITCHE .EQ. 'ENRICHIS') THEN
  50.             ICHELX = i
  51.             MCHELX = MCHELM
  52.             GOTO 10
  53.           ENDIF
  54.           SEGDES,MCHELM
  55.         ENDIF
  56.         IF (IMODEL.TYMODE(i) .EQ. 'MCHAMLL') THEN
  57.  
  58.           MCHAML = IMODEL.IVAMOD(i)
  59.  
  60.           SEGACT,MCHAML
  61. C++++  INITIALISATION d'un nouveau chamelem chapeau vierge
  62.         L1 = 8
  63. C nombre de sous champs
  64.         N1 = 1
  65. C nombre de composantes
  66.         N2 = 3  
  67.         N3 = 6
  68.  
  69.         SEGINI, MCHEX1
  70.         MCHEX1.TITCHE = 'ENRICHIS'
  71.         MCHEX1.IFOCHE = IFOUR
  72.         MCHEX1.IMACHE(1)= IMAMOD
  73.         MCHEX1.CONCHE(1)= IMODEL.CONMOD
  74.         MCHEX1.INFCHE(1,2) = 0
  75.         MCHEX1.INFCHE(1,3) = NIFOUR
  76.         MCHEX1.INFCHE(1,6) = 1
  77.         ICHELX = i
  78.  
  79.  
  80.         MCHEX1.ICHAML(1)=MCHAML
  81.  
  82.         MCHELX = MCHEX1
  83.         ICHELX = i
  84.             GOTO 10
  85.         ENDIF
  86.       ENDDO
  87. c      IF (ICHELX .EQ. 0) THEN
  88. c        WRITE(IOIMP,*) 'Il n y a pas de MCHAML d ENRICHISsement ',
  89. c     &                 'dans le modele'
  90. c        RETURN
  91. c      ENDIF
  92.  
  93.  10   CONTINUE
  94. C*      SEGDES,IMODEL
  95.  
  96.       RETURN
  97.       END
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  

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