Télécharger iddefo.eso

Retour à la liste

iddefo
C IDDEFO    SOURCE    OF166741  24/10/07    21:15:24     12016          
 
C=======================================================================
C=       RECHERCHE DES NOMS DE COMPOSANTES DE DEFORMATIONS             =
C=       -------------------------------------------------             =
C=                                                                     =
C=  Entrees :                                                          =
C=    IPMODL    pointeur ACTIF sur objet IMODEL                        =
C=    IFOU      valeur de IFOUR de CCOPTIO                             =
C=                                                                     =
C=  Sorties :                                                          =
C=    IPNOMC    pointeur de type NOMID sur les listes de noms de       =
C=              composantes OBLigatoires et FACultatives               =
C=    NBROBL    nombre de composantes OBLigatoires                     =
C=    NBRFAC    nombre de composantes FACultatives                     =
C=======================================================================
 
      SUBROUTINE IDDEFO (IPMODL,IFOU,IPNOMC,NBROBL,NBRFAC)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
C==DEB= FORMULATION HHO == Include specifique ==========================
-INC CCHHOPA
C==FIN= FORMULATION HHO ================================================
 
-INC SMMODEL
      POINTEUR nomid1.nomid
 
      NBROBL=0
      NBRFAC=0
      NOMID=0
 
      IMODEL=IPMODL
      MELE=NEFMOD
      MFR=NUMMFR(MELE)
 
      IF (FORMOD(1).eq.'MECANIQUE') THEN
        DO imat = 1 , matmod(/2)
          if (matmod(imat).eq.'IMPEDANCE') then
            MFR2 = infele(13)
            if (mfr2.eq.26.or.mfr2.eq.28) mfr = mfr2
 
C           Recuperation du NOMID des composantes primales
            NOMID1=IMODEL.LNOMID(1)
            nbrobl=NOMID1.LESOBL(/2)
            nbrfac=NOMID1.LESFAC(/2)
            segini,NOMID
            do io = 1,nbrobl
              LESOBL(io)='E'//NOMID1.LESOBL(io)(1:3)
            enddo
            do io = 1,nbrfac
              LESFAC(io)='E'//NOMID1.LESFAC(io)(1:3)
            enddo
 
            IPNOMC=NOMID
 
            RETURN
          endif
        ENDDO
      ENDIF
      IF (FORMOD(1).EQ.'ELECTROSTATIQUE ') MFR = 71
      IF (FORMOD(1).EQ.'DIFFUSION       ') MFR = 73
 
      NFOR = FORMOD(/2)
      CALL PLACE(FORMOD,NFOR,ITHHY,'THERMOHYDRIQUE  ')
      CALL PLACE(FORMOD,NFOR,ITHER,'THERMIQUE       ')
      CALL PLACE(FORMOD,NFOR,IMAGN,'MAGNETODYNAMIQUE')
      CALL PLACE(FORMOD,NFOR,IMETA,'METALLURGIE     ')
 
      mfr2 = 0
      IF (ITHHY.EQ.1) mfr2=65
      IF (ITHER.eq.1) mfr2=29
      IF (IMAGN.eq.1) mfr2=69
      IF (IMETA.eq.1) mfr2=70
      IF (mfr2.ne.0) go to 100
 
C  Coque integree ou pas ?
      NPINT=INFMOD(1)
*
*     -1/ MACRO ELEMENT
*
      IF (MFR.EQ.61)THEN
         IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN
*
*        0-A/ CONTRAINTES PLANES - DEFORMATIONS PLANES
*
            NBROBL=5
            SEGINI NOMID
            LESOBL(1)='EPS1'
            LESOBL(2)='CZ1 '
            LESOBL(3)='GXY '
            LESOBL(4)='EPS2'
            LESOBL(5)='CZ2 '
         ENDIF
c
c    element coaxial COS2 (3D pour liaison acier-beton)
c
      ELSE IF ( MFR.EQ.78) THEN
           NBROBL=3
           SEGINI NOMID
           LESOBL(1)='GLIT'
           LESOBL(2)='GLIN'
           LESOBL(3)='GLIM'
c
C  1 - Element JOINT
C ===================
      ELSE IF (MFR.EQ.35) THEN
C =====
C  1.1 - Tridimensionnel
C =====
        IF (IFOU.EQ.2) THEN
          NBROBL=3
          SEGINI,NOMID
          LESOBL(1)='DRS1'
          LESOBL(2)='DRS2'
          LESOBL(3)='DRN '
C =====
C  1.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie
C =====
        ELSE IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR.
     .           (IFOU.EQ.0)) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='DRSN'
          LESOBL(2)='DRN '
        ENDIF
 
C  2 - Element JOINT CISAILLEMENT (2D)
C      Bidimensionnel PLAN (CP/DP/DPGE)
C ======================================
      ELSE IF (MFR.EQ.53) THEN
        IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1)) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='DRS1'
          LESOBL(2)='DRS2'
        ENDIF
 
C  3 - Element JOINT GENERALISE
C ==============================
      ELSE IF (MFR.EQ.55) THEN
C =====
C  3.1 - Bidimensionnel PLAN (CP/DP/DPGE)
C =====
        IF ((IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR.(IFOU.EQ.-3)) THEN
          NBROBL=4
          SEGINI,NOMID
          LESOBL(1)='DRSN'
          LESOBL(2)='DRN '
          LESOBL(3)='EPTT'
          LESOBL(4)='EPZZ'
C =====
C  3.2 - Tridimensionnel
C =====
        ELSE IF (IFOU.EQ.2) THEN
          NBROBL=6
          SEGINI,NOMID
          LESOBL(1)='DRS1'
          LESOBL(2)='DRS2'
          LESOBL(3)='DRN '
          LESOBL(4)='EPTT'
          LESOBL(5)='EPUU'
          LESOBL(6)='GATU'
        ENDIF
 
C  4 - Elements MASSIFS/INCOMPRESSIBLES/XFEM/HHO
C ===============================================
      ELSE IF (MFR.EQ.1.OR.MFR.EQ.31.OR.MFR.EQ.63.OR.
     &         MFR.EQ.HHO_MFR_ELEMENT) THEN
C =====
C  4.1 - Bidimensionnel PLAN (CP/DP/DPGE)
C =====
        IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
          NBROBL=4
          SEGINI,NOMID
          LESOBL(1)='EPXX'
          LESOBL(2)='EPYY'
          LESOBL(3)='EPZZ'
          LESOBL(4)='GAXY'
C =====
C  4.2 - Axisymetrie
C =====
        ELSE IF (IFOU.EQ.0) THEN
          NBROBL=4
          SEGINI,NOMID
          LESOBL(1)='EPRR'
          LESOBL(2)='EPZZ'
          LESOBL(3)='EPTT'
          LESOBL(4)='GARZ'
C =====
C  4.3 - Fourier
C =====
        ELSE IF (IFOU.EQ.1) THEN
          NBROBL=6
          SEGINI,NOMID
          LESOBL(1)='EPRR'
          LESOBL(2)='EPZZ'
          LESOBL(3)='EPTT'
          LESOBL(4)='GARZ'
          LESOBL(5)='GART'
          LESOBL(6)='GAZT'
C =====
C  4.4 - Tridimensionnel
C =====
        ELSE IF (IFOU.EQ.2) THEN
          NBROBL=6
          SEGINI,NOMID
          LESOBL(1)='EPXX'
          LESOBL(2)='EPYY'
          LESOBL(3)='EPZZ'
          LESOBL(4)='GAXY'
          LESOBL(5)='GAXZ'
          LESOBL(6)='GAYZ'
C =====
C  4.5 - Unidimensionnel (1D)
C =====
        ELSE IF (IFOU.GE.3.AND.IFOU.LE.15.AND.MFR.EQ.1) THEN
          NBROBL=3
          SEGINI,NOMID
          IF (IFOU.LE.11) THEN
            LESOBL(1)='EPXX'
            LESOBL(2)='EPYY'
            LESOBL(3)='EPZZ'
          ELSE IF (IFOU.GE.12.AND.IFOU.LE.14) THEN
            LESOBL(1)='EPRR'
            LESOBL(2)='EPZZ'
            LESOBL(3)='EPTT'
          ELSE IF (IFOU.EQ.15) THEN
            LESOBL(1)='EPRR'
C  Voir TC pour LESOBL(2)='EPFF' ?
            LESOBL(2)='EPZZ'
            LESOBL(3)='EPTT'
          ENDIF
        ENDIF
 
C  5 - Elements COQUES
C =====================
      ELSE IF (MFR.EQ.3.OR.MFR.EQ.39) THEN
        IF (NPINT.EQ.0) THEN
C =====
C  5.1 - Tridimensionnel et Fourier
C =====
          IF (IFOU.EQ.2.OR.IFOU.EQ.1) THEN
            NBROBL=6
            SEGINI,NOMID
            LESOBL(1)='EPSS'
            LESOBL(2)='EPTT'
            LESOBL(3)='GAST'
            LESOBL(4)='RTSS'
            LESOBL(5)='RTTT'
            LESOBL(6)='RTST'
C =====
C  5.2 - Axisymetrie
C =====
          ELSE IF (IFOU.EQ.0) THEN
            NBROBL=4
            SEGINI,NOMID
            LESOBL(1)='EPSS'
            LESOBL(2)='EPTT'
            LESOBL(3)='RTSS'
            LESOBL(4)='RTTT'
C =====
C  5.3 - Bidimensionnel PLAN (CP/DP/DPGE)
C =====
          ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
            NBROBL=4
            SEGINI,NOMID
            LESOBL(1)='EPSS'
            LESOBL(2)='EPZZ'
            LESOBL(3)='RTSS'
            LESOBL(4)='RTZZ'
          ENDIF
        ELSE
C =====
C  5.4 - Tridimensionnel - Cas particulier
C =====
          IF (MFR.EQ.3.AND.IFOU.EQ.2) THEN
            NBROBL=4
            SEGINI,NOMID
            LESOBL(1)='EPSS'
            LESOBL(2)='EPTT'
            LESOBL(3)='EPNN'
            LESOBL(4)='GAST'
          ENDIF
        ENDIF
 
C  6 - Elements de COQUE EPAISSE
C      Tridimensionnel, axisymetrie et Fourier
C =============================================
      ELSE IF (MFR.EQ.5) THEN
        NBROBL=5
        SEGINI,NOMID
        LESOBL(1)='EPSS'
        LESOBL(2)='EPTT'
        LESOBL(3)='GAST'
        LESOBL(4)='GASN'
        LESOBL(5)='GATN'
 
C  7 - LIA2: Element 3D de liaison a 2 noeuds
C =============================================
      ELSE IF (MFR.EQ.51) THEN
        NBROBL=6
        SEGINI,NOMID
        LESOBL(1)='EPS '
        LESOBL(2)='GXY '
        LESOBL(3)='GXZ '
        LESOBL(4)='CX  '
        LESOBL(5)='CY  '
        LESOBL(6)='CZ  '
C  7bis - JOI1: Element 3D de liaison a 2 noeuds
C =============================================
      ELSE IF (MFR.EQ.75) THEN
       IF (IFOU.EQ.2) THEN
        NBROBL=6
        SEGINI,NOMID
        LESOBL(1)='EXX '
        LESOBL(2)='GXY '
        LESOBL(3)='GXZ '
        LESOBL(4)='CXX  '
        LESOBL(5)='CXY  '
        LESOBL(6)='CXZ  '
       ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
        NBROBL=3
        SEGINI,NOMID
        LESOBL(1)='EXX '
        LESOBL(2)='GXY '
        LESOBL(3)='CXZ '
       ENDIF
 
C  8 - BAEX : BARRE EXCENTREE
C ============================
      ELSE IF (MFR.EQ.49) THEN
        NBROBL=1
        SEGINI,NOMID
        LESOBL(1)='EPSS'
 
C  9 - Elements POUTRE et TUYAU
C ==============================
      ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN
C =====
C  9.1 - Tridimensionnel
C =====
        IF (IFOU.EQ.2) THEN
          NBROBL=6
          SEGINI,NOMID
          LESOBL(1)='EPS '
          LESOBL(2)='GXY '
          LESOBL(3)='GXZ '
          LESOBL(4)='CX  '
          LESOBL(5)='CY  '
          LESOBL(6)='CZ  '
C =====
C  9.2 - Bidimensionnel PLAN (CP/DP/DPGE)
C =====
        ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
          NBROBL=3
          SEGINI,NOMID
          LESOBL(1)='EPS '
          LESOBL(2)='GXY '
          LESOBL(3)='CZ  '
        ENDIF
 
C  10 - Element de COQUE en CISAILLEMENT TRANSVERSE
C ==================================================
      ELSE IF (MFR.EQ.9) THEN
        IF (NPINT.EQ.0) THEN
C ======
C  10.1 - Tridimensionnel
C ======
          IF (IFOU.EQ.2) THEN
            NBROBL=8
            SEGINI,NOMID
            LESOBL(1)='EPSS'
            LESOBL(2)='EPTT'
            LESOBL(3)='GAST'
            LESOBL(4)='RTSS'
            LESOBL(5)='RTTT'
            LESOBL(6)='RTST'
            LESOBL(7)='GASN'
            LESOBL(8)='GATN'
          ENDIF
        ENDIF
 
C  11 - Elements LINESPRING (tridimensionnel)
C ============================================
      ELSE IF (MFR.EQ.15) THEN
        IF (IFOU.EQ.2) THEN
          NBROBL=6
          SEGINI,NOMID
          LESOBL(1)='EPZZ'
          LESOBL(2)='GAXZ'
          LESOBL(3)='GAYZ'
          LESOBL(4)='RTXX'
          LESOBL(5)='RTZZ'
          LESOBL(6)='DJP '
        ENDIF
 
C  12 - Elements de MEMBRANE
C ===========================
      ELSE IF (MFR.EQ.25) THEN
        NBROBL=3
        SEGINI,NOMID
        LESOBL(1)='EPSS'
        LESOBL(2)='EPTT'
        LESOBL(3)='GAST'
 
C  13 - Formulation UNIAXIALE
C =============================
      ELSE IF (MFR.EQ.27) THEN
        NBROBL=1
        SEGINI,NOMID
        LESOBL(1)='EPSS'
 
C  14 - Formulation THERMIQUE
C ============================
      ELSE IF (MFR.EQ.29) THEN
C ======
C  14.1 - Tridimensionnel
C ======
        IF (IFOU.EQ.2) THEN
          NBROBL=3
          SEGINI,NOMID
          LESOBL(1)='GRDX'
          LESOBL(2)='GRDY'
          LESOBL(3)='GRDZ'
C ======
C  14.2 - Bidimensionnel PLAN (CP/DP/DPGE)
C ======
        ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='GRDX'
          LESOBL(2)='GRDY'
C ======
C  14.3 - Axisymetrie et Fourier
C ======
        ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='GRDX'
          LESOBL(2)='GRDY'
C*OF      LESOBL(1)='GRDR'
C*OF      LESOBL(2)='GRDZ'
C ======
C  14.4 - Unidimensionnel (1D)
C ======
        ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN
          NBROBL=1
          SEGINI,NOMID
C*OF      IF (IFOU.LE.11) THEN
            LESOBL(1)='GRDX'
C*OF      ELSE
C*OF        LESOBL(1)='GRDR'
C*OF      ENDIF
        ENDIF
 
C  15 - Formulation POREUX
C =========================
      ELSE IF (MFR.EQ.33) THEN
C =====
C  15.A - Elements POREUX
C ======
        IF (MELE.GE.79.AND.MELE.LE.83) THEN
C ========
C  15.A.1 - Bidimensionnel PLAN (CP/DP)
C ========
          IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
            NBROBL=5
            SEGINI,NOMID
            LESOBL(1)='EPXX'
            LESOBL(2)='EPYY'
            LESOBL(3)='EPZZ'
            LESOBL(4)='GAXY'
            LESOBL(5)='PRES'
C ========
C  15.A.2 - Axisymetrie
C ========
          ELSE IF (IFOU.EQ.0) THEN
            NBROBL=5
            SEGINI,NOMID
            LESOBL(1)='EPRR'
            LESOBL(2)='EPZZ'
            LESOBL(3)='EPTT'
            LESOBL(4)='GARZ'
            LESOBL(5)='PRES'
C ========
C  15.A.3 - Fourier
C ========
          ELSE IF (IFOU.EQ.1) THEN
            NBROBL=7
            SEGINI,NOMID
            LESOBL(1)='EPRR'
            LESOBL(2)='EPZZ'
            LESOBL(3)='EPTT'
            LESOBL(4)='GARZ'
            LESOBL(5)='GART'
            LESOBL(6)='GAZT'
            LESOBL(7)='PRES'
C ========
C  15.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
C ========
          ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
            NBROBL=7
            SEGINI,NOMID
            LESOBL(1)='EPXX'
            LESOBL(2)='EPYY'
            LESOBL(3)='EPZZ'
            LESOBL(4)='GAXY'
            LESOBL(5)='GAXZ'
            LESOBL(6)='GAYZ'
            LESOBL(7)='PRES'
          ENDIF
C ======
C  15.B - Elements JOINT POREUX
C ======
        ELSE IF (MELE.GE.108.AND.MELE.LE.110) THEN
C ========
C  15.B.1 - Tridimensionnel
C ========
          IF (IFOU.EQ.2) THEN
            NBROBL=4
            SEGINI,NOMID
            LESOBL(1)='DRS1'
            LESOBL(2)='DRS2'
            LESOBL(3)='DRN '
            LESOBL(4)='PRES'
C ========
C  15.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
C ========
          ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
            NBROBL=3
            SEGINI,NOMID
            LESOBL(1)='DRSN'
            LESOBL(2)='DRN '
            LESOBL(3)='PRES'
          ENDIF
        ENDIF
 
C  16 - Formulation POREUX type Q
C ================================
      ELSE IF (MFR.EQ.57) THEN
C =====
C  16.A - Elements POREUX Q
C ======
        IF (MELE.GE.173.AND.MELE.LE.177) THEN
C ========
C  16.A.1 - Bidimensionnel PLAN (CP/DP)
C ========
          IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
            NBROBL=6
            SEGINI,NOMID
            LESOBL(1)='EPXX'
            LESOBL(2)='EPYY'
            LESOBL(3)='EPZZ'
            LESOBL(4)='GAXY'
            LESOBL(5)='PRES'
            LESOBL(6)='PREQ'
C ========
C  16.A.2 - Axisymetrie
C ========
          ELSE IF (IFOU.EQ.0) THEN
            NBROBL=6
            SEGINI,NOMID
            LESOBL(1)='EPRR'
            LESOBL(2)='EPZZ'
            LESOBL(3)='EPTT'
            LESOBL(4)='GARZ'
            LESOBL(5)='PRES'
            LESOBL(6)='PREQ'
C ========
C  16.A.3 - Fourier
C ========
          ELSE IF (IFOU.EQ.1) THEN
            NBROBL=8
            SEGINI,NOMID
            LESOBL(1)='EPRR'
            LESOBL(2)='EPZZ'
            LESOBL(3)='EPTT'
            LESOBL(4)='GARZ'
            LESOBL(5)='GART'
            LESOBL(6)='GAZT'
            LESOBL(7)='PRES'
            LESOBL(8)='PREQ'
C ========
C  16.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
C ========
          ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
            NBROBL=8
            SEGINI,NOMID
            LESOBL(1)='EPXX'
            LESOBL(2)='EPYY'
            LESOBL(3)='EPZZ'
            LESOBL(4)='GAXY'
            LESOBL(5)='GAXZ'
            LESOBL(6)='GAYZ'
            LESOBL(7)='PRES'
            LESOBL(8)='PREQ'
          ENDIF
C ======
C  16.B - Elements JOINT POREUX Q
C ======
        ELSE IF (MELE.GE.185.AND.MELE.LE.187) THEN
C ========
C  16.B.1 - Tridimensionnel
C ========
          IF (IFOU.EQ.2) THEN
            NBROBL=5
            SEGINI,NOMID
            LESOBL(1)='DRS1'
            LESOBL(2)='DRS2'
            LESOBL(3)='DRN '
            LESOBL(4)='PRES'
            LESOBL(5)='PREQ'
C ========
C  16.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
C ========
          ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
            NBROBL=4
            SEGINI,NOMID
            LESOBL(1)='DRSN'
            LESOBL(2)='DRN '
            LESOBL(3)='PRES'
            LESOBL(4)='PREQ'
          ENDIF
        ENDIF
 
C  17 - Formulation POREUX type R
C ================================
      ELSE IF (MFR.EQ.59) THEN
C =====
C  17.A - Elements POREUX R
C ======
        IF (MELE.GE.178.AND.MELE.LE.182) THEN
C ========
C  17.A.1 - Bidimensionnel PLAN (CP/DP)
C ========
          IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
            NBROBL=7
            SEGINI NOMID
            LESOBL(1)='EPXX'
            LESOBL(2)='EPYY'
            LESOBL(3)='EPZZ'
            LESOBL(4)='GAXY'
            LESOBL(5)='PRES'
            LESOBL(6)='PREQ'
            LESOBL(7)='TPT '
C ========
C  17.A.2 - Axisymetrie
C ========
          ELSE IF (IFOU.EQ.0) THEN
            NBROBL=7
            SEGINI NOMID
            LESOBL(1)='EPRR'
            LESOBL(2)='EPZZ'
            LESOBL(3)='EPTT'
            LESOBL(4)='GARZ'
            LESOBL(5)='PRES'
            LESOBL(6)='PREQ'
            LESOBL(7)='TPT '
C ========
C  17.A.3 - Fourier
C ========
          ELSE IF (IFOU.EQ.1) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='EPRR'
            LESOBL(2)='EPZZ'
            LESOBL(3)='EPTT'
            LESOBL(4)='GARZ'
            LESOBL(5)='GART'
            LESOBL(6)='GAZT'
            LESOBL(7)='PRES'
            LESOBL(8)='PREQ'
            LESOBL(9)='TPT '
C ========
C  17.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
C ========
          ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
            NBROBL=9
            SEGINI,NOMID
            LESOBL(1)='EPXX'
            LESOBL(2)='EPYY'
            LESOBL(3)='EPZZ'
            LESOBL(4)='GAXY'
            LESOBL(5)='GAXZ'
            LESOBL(6)='GAYZ'
            LESOBL(7)='PRES'
            LESOBL(8)='PREQ'
            LESOBL(9)='TPT '
          ENDIF
C ======
C  17.B - Elements JOINT POREUX R
C ======
        ELSE IF (MELE.GE.188.AND.MELE.LE.190) THEN
C ========
C  17.B.1 - Tridimensionnel
C ========
          IF (IFOU.EQ.2) THEN
            NBROBL=6
            SEGINI,NOMID
            LESOBL(1)='DRS1'
            LESOBL(2)='DRS2'
            LESOBL(3)='DRN '
            LESOBL(4)='PRES'
            LESOBL(5)='PREQ'
            LESOBL(6)='TPT '
C ========
C  17.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
C ========
          ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
            NBROBL=5
            SEGINI,NOMID
            LESOBL(1)='DRSN'
            LESOBL(2)='DRN '
            LESOBL(3)='PRES'
            LESOBL(4)='PREQ'
            LESOBL(5)='TPT '
          ENDIF
        ENDIF
 
C  18 - Elements TUYAU FISSURE (tridimensionnel)
C ===============================================
      ELSE IF (MFR.EQ.17) THEN
        IF (IFOU.EQ.2) THEN
          NBROBL=8
          SEGINI,NOMID
          LESOBL(1)='EPS '
          LESOBL(2)='GXY '
          LESOBL(3)='GXZ '
          LESOBL(4)='CX  '
          LESOBL(5)='CY  '
          LESOBL(6)='CZ  '
          LESOBL(7)='EPS7'
          LESOBL(8)='EPS8'
        ENDIF
 
C  19 - Elements de SECTION/FIBRE
C ================================
      ELSE IF (MFR.EQ.47) THEN
C ======
C  19.1 - Bidimensionnel PLAN (CP/DP/DPGE)
C ======
        IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='EPXX'
          LESOBL(2)='GAXY'
C ======
C  19.2 - Tridimensionnel
C ======
        ELSE IF (IFOU.EQ.2) THEN
          NBROBL=3
          SEGINI,NOMID
          LESOBL(1)='EPXX'
          LESOBL(2)='GAXY'
          LESOBL(3)='GAXZ'
        ENDIF
 
C  20 - Elements de zones cohesives
C ==================================
      ELSE IF (MFR.EQ.77) THEN
C =====
C  20.1 - Tridimensionnel
C =====
        IF (IFOU.GT.0) THEN
          NBROBL=3
          SEGINI,NOMID
          LESOBL(1)='DRS1'
          LESOBL(2)='DRS2'
          LESOBL(3)='DRN '
C =====
C  20.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie
C =====
        ELSE IF (IFOU.LE.0) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='DRSN'
          LESOBL(2)='DRN '
        ENDIF
 
C  21 - Formulation ELECTROSTATIQUE (base MASSIF)
C ==================================
      ELSE IF (MFR.EQ.71) THEN
C ======
C  21.1 - Tridimensionnel
C ======
        IF (IFOU.EQ.2) THEN
          NBROBL=3
          SEGINI,NOMID
          LESOBL(1)='EELX    '
          LESOBL(2)='EELY    '
          LESOBL(3)='EELZ    '
C ======
C  21.2 - Bidimensionnel PLAN (CP/DP/DPGE)
C ======
        ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='EELX    '
          LESOBL(2)='EELY    '
C ======
C  21.3 - Axisymetrie et Fourier
C ======
        ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN
          NBROBL=2
          SEGINI,NOMID
          LESOBL(1)='EELR    '
          LESOBL(2)='EELZ    '
C ======
C  21.4 - Unidimensionnel (1D)
C ======
        ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN
          NBROBL=1
          SEGINI,NOMID
          IF (IFOU.LE.11) THEN
            LESOBL(1)='EELX    '
          ELSE
            LESOBL(1)='EELR    '
          ENDIF
        ENDIF
 
C  22 - Formulation DIFFUSION (Nomid VIDE)
C ========================================
      ELSE IF (MFR.EQ.73) THEN
 
C*DEB = Merci de conserver ces lignes ! (Developpements futurs)
C*XC ======
C*XC  22.0 - BARRE
C*XC ======
C*X       IF (MELE.EQ.46) THEN
C*X          NBROBL=1
C*X          SEGINI,NOMID
C*X          LESOBL(1)='GS      '
C*X       ELSE
C*XC ======
C*XC  22.1 - Tridimensionnel
C*XC ======
C*X        IF (IFOU.EQ.2) THEN
C*X          NBROBL=3
C*X          SEGINI,NOMID
C*X          LESOBL(1)='GX      '
C*X          LESOBL(2)='GY      '
C*X          LESOBL(3)='GZ      '
C*XC ======
C*XC  22.2 - Bidimensionnel PLAN (CP/DP/DPGE)
C*XC ======
C*X        ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
C*X          NBROBL=2
C*X          SEGINI,NOMID
C*X          LESOBL(1)='GX      '
C*X          LESOBL(2)='GY      '
C*XC ======
C*XC  22.3 - Axisymetrie et Fourier
C*XC ======
C*X        ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN
C*X          NBROBL=2
C*X          SEGINI,NOMID
C*X          LESOBL(1)='GR      '
C*X          LESOBL(2)='GZ      '
C*XC ======
C*XC  22.4 - Unidimensionnel (1D)
C*XC ======
C*X        ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN
C*X          NBROBL=1
C*X          SEGINI,NOMID
C*X          IF (IFOU.LE.11) THEN
C*X            LESOBL(1)='GX      '
C*X          ELSE
C*X            LESOBL(1)='GR      '
C*X          ENDIF
C*X        ENDIF
C*X       ENDIF
C*XC* Ajout du nom de l'INCOnnue (DDL) contenu dans TYMODE(1)
C*X       DO i = 1, NBROBL
C*X         LESOBL(i)(3:8) = imodel.TYMODE(1)(1:6)
C*X       ENDDO
C*FIN = Merci de conserver ces lignes ! (Developpements futurs)
 
      ENDIF
 
C  Par DEFAUT : (Nomid VIDE)
C ===========================
  100 continue
      IF (NOMID.EQ.0) THEN
        SEGINI,NOMID
      ELSE
        if (ifomod.eq.6) then
          nbrfa0 = nbrfac
          NBRFAC = NBROBL + (nbrfa0*2)
          segadj,nomid
          do imo = 1,nbrobl
            lesfac(nbrfa0 + imo)(2:4) = lesobl(imo)(1:3)
            lesfac(nbrfa0 + imo)(1:1) = 'I'
          enddo
          do imo = 1,nbrfa0
            lesfac(nbrfa0+nbrobl+imo)(2:4) = lesfac(imo)(1:3)
            lesfac(nbrfa0+nbrobl+imo)(1:1) = 'I'
          enddo
        endif
      ENDIF
 
      SEGACT,NOMID*NOMOD
      IPNOMC=NOMID
 
c      RETURN
      END
 
 
 
 

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