Télécharger la plaquette Cast3M


fr en

Télécharger iddein.eso

Retour à la liste

Numérotation des lignes : 

iddein

  1. C IDDEIN    SOURCE    OF166741  24/05/06    21:15:17     11082          
  2.  
  3. C=======================================================================
  4. C=   RECHERCHE DES NOMS DE COMPOSANTES DE DEFORMATIONS INELASTIQUES    =
  5. C=   --------------------------------------------------------------    =
  6. C=                                                                     =
  7. C=  Entrees :                                                          =
  8. C=    IPMODL    pointeur ACTIF sur objet IMODEL                        =
  9. C=    IFOU      valeur de IFOUR de CCOPTIO                             =
  10. C=                                                                     =
  11. C=  Sorties :                                                          =
  12. C=    IPNOMC    pointeur de type NOMID sur les listes de noms de       =
  13. C=              composantes OBLigatoires et FACultatives               =
  14. C=    NBROBL    nombre de composantes OBLigatoires                     =
  15. C=    NBRFAC    nombre de composantes FACultatives                     =
  16. C=======================================================================
  17.  
  18.       SUBROUTINE IDDEIN (IPMODL,IFOU,IPNOMC,NBROBL,NBRFAC)
  19.  
  20.       IMPLICIT INTEGER(I-N)
  21.       IMPLICIT REAL*8(A-H,O-Z)
  22.  
  23. -INC PPARAM
  24. -INC CCOPTIO
  25. C==DEB= FORMULATION HHO == Include specifique ==========================
  26. -INC CCHHOPA
  27. C==FIN= FORMULATION HHO ================================================
  28.  
  29. -INC SMMODEL
  30.  
  31.       IPNOMC=0
  32.       NBROBL=0
  33.       NBRFAC=0
  34.       NOMID =0
  35.  
  36.       IMODEL=IPMODL
  37.  
  38.       MELE=NEFMOD
  39.       MFR=NUMMFR(MELE)
  40.  
  41.       NFOR=FORMOD(/2)
  42.       CALL PLACE(FORMOD,NFOR,ITHHY,'THERMOHYDRIQUE')
  43.       CALL PLACE(FORMOD,NFOR,ITHER,'THERMIQUE')
  44.       CALL PLACE(FORMOD,NFOR,IMAGN,'MAGNETODYNAMIQUE')
  45.       CALL PLACE(FORMOD,NFOR,IELEC,'ELECTROSTATIQUE')
  46.       CALL PLACE(FORMOD,NFOR,IDIFF,'DIFFUSION')
  47.       CALL PLACE(FORMOD,NFOR,IMET,'METALLURGIE')
  48.  
  49.       mfr2=0
  50.       IF (ITHHY.NE.0) mfr2=65
  51.       if (ITHER.NE.0) mfr2=29
  52.       if (IMAGN.NE.0) mfr2=69
  53.       IF (IELEC.NE.0) mfr2=71
  54.       IF (IDIFF.NE.0) mfr2=73
  55.       if (IMET.NE.0) mfr2=29
  56.       IF (mfr2.NE.0) GOTO 100
  57.  
  58. C  Coque integree ou pas ?
  59.       NPINT=0
  60.       IF (INFMOD(/1).NE.0) NPINT=INFMOD(1)
  61. *
  62. *     -1/ MACRO ELEMENT
  63. *
  64.       IF (MFR.EQ.61)THEN
  65. C
  66.          IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN
  67. *
  68. *        0-A/ CONTRAINTES PLANES - DEFORMATIONS PLANES
  69. *
  70.             NBROBL=5
  71.             SEGINI NOMID
  72.             LESOBL(1)='EIS1'
  73.             LESOBL(2)='CIZ1'
  74.             LESOBL(3)='GIXY'
  75.             LESOBL(4)='EIS2'
  76.             LESOBL(5)='CIZ2'
  77.          ENDIF
  78.  
  79. C  1 - Element JOINT
  80. C ===================
  81.       ELSE IF (MFR.EQ.35) THEN
  82. C =====
  83. C  1.1 - Tridimensionnel
  84. C =====
  85.         IF (IFOU.EQ.2) THEN
  86.           NBROBL=3
  87.           SEGINI,NOMID
  88.           LESOBL(1)='FRS1'
  89.           LESOBL(2)='FRS2'
  90.           LESOBL(3)='FRN '
  91. C =====
  92. C  1.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie
  93. C =====
  94.         ELSE IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR.
  95.      .           (IFOU.EQ.0)) THEN
  96.           NBROBL=2
  97.           SEGINI,NOMID
  98.           LESOBL(1)='FRSN'
  99.           LESOBL(2)='FRN '
  100.         ENDIF
  101.  
  102. C  2 - Element JOINT CISAILLEMENT (2D)
  103. C      Bidimensionnel PLAN (CP/DP/DPGE)
  104. C ======================================
  105.       ELSE IF (MFR.EQ.53) THEN
  106.         IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1)) THEN
  107.           NBROBL=2
  108.           SEGINI,NOMID
  109.           LESOBL(1)='FRS1'
  110.           LESOBL(2)='FRS2'
  111.         ENDIF
  112.  
  113. C  3 - Element JOINT GENERALISE
  114. C ==============================
  115.       ELSE IF (MFR.EQ.55) THEN
  116. C =====
  117. C  3.1 - Bidimensionnel PLAN (CP/DP/DPGE)
  118. C =====
  119.         IF ((IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR.(IFOU.EQ.-3)) THEN
  120.           NBROBL=4
  121.           SEGINI,NOMID
  122.           LESOBL(1)='FRSN'
  123.           LESOBL(2)='FRN '
  124.           LESOBL(3)='EITT'
  125.           LESOBL(4)='EIZZ'
  126. C =====
  127. C  3.2 - Tridimensionnel
  128. C =====
  129.         ELSE IF (IFOU.EQ.2) THEN
  130.           NBROBL=6
  131.           SEGINI,NOMID
  132.           LESOBL(1)='FRS1'
  133.           LESOBL(2)='FRS2'
  134.           LESOBL(3)='FRN '
  135.           LESOBL(4)='EITT'
  136.           LESOBL(5)='EIUU'
  137.           LESOBL(6)='GITU'
  138.         ENDIF
  139.  
  140. C  4 - Elements MASSIFS + FORMULATION HHO
  141. C ======================
  142.       ELSE IF (MFR.EQ.1.OR.MFR.EQ.31.or.MFR.EQ.63.or.
  143.      &         MFR.EQ.HHO_MFR_ELEMENT) THEN
  144. C =====
  145. C  4.1 - Bidimensionnel PLAN (CP/DP/DPGE)
  146. C =====
  147.         IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  148.           NBROBL=4
  149.           SEGINI,NOMID
  150.           LESOBL(1)='EIXX'
  151.           LESOBL(2)='EIYY'
  152.           LESOBL(3)='EIZZ'
  153.           LESOBL(4)='GIXY'
  154. C =====
  155. C  4.2 - Axisymetrie
  156. C =====
  157.         ELSE IF (IFOU.EQ.0) THEN
  158.           NBROBL=4
  159.           SEGINI,NOMID
  160.           LESOBL(1)='EIRR'
  161.           LESOBL(2)='EIZZ'
  162.           LESOBL(3)='EITT'
  163.           LESOBL(4)='GIRZ'
  164. C =====
  165. C  4.3 - Fourier
  166. C =====
  167.         ELSE IF (IFOU.EQ.1) THEN
  168.           NBROBL=6
  169.           SEGINI,NOMID
  170.           LESOBL(1)='EIRR'
  171.           LESOBL(2)='EIZZ'
  172.           LESOBL(3)='EITT'
  173.           LESOBL(4)='GIRZ'
  174.           LESOBL(5)='GIRT'
  175.           LESOBL(6)='GIZT'
  176. C =====
  177. C  4.4 - Tridimensionnel
  178. C =====
  179.         ELSE IF (IFOU.EQ.2) THEN
  180.           NBROBL=6
  181.           SEGINI,NOMID
  182.           LESOBL(1)='EIXX'
  183.           LESOBL(2)='EIYY'
  184.           LESOBL(3)='EIZZ'
  185.           LESOBL(4)='GIXY'
  186.           LESOBL(5)='GIXZ'
  187.           LESOBL(6)='GIYZ'
  188. C =====
  189. C  4.5 - Unidimensionnel (1D)
  190. C =====
  191.         ELSE IF (IFOU.GE.3.AND.IFOU.LE.15.AND.MFR.EQ.1) THEN
  192.           NBROBL=3
  193.           SEGINI,NOMID
  194.           IF (IFOU.LE.11) THEN
  195.             LESOBL(1)='EIXX'
  196.             LESOBL(2)='EIYY'
  197.             LESOBL(3)='EIZZ'
  198.           ELSE IF (IFOU.GE.12.AND.IFOU.LE.14) THEN
  199.             LESOBL(1)='EIRR'
  200.             LESOBL(2)='EIZZ'
  201.             LESOBL(3)='EITT'
  202.           ELSE IF (IFOU.EQ.15) THEN
  203.             LESOBL(1)='EIRR'
  204. C  Voir TC pour LESOBL(2)='EIFF'
  205.             LESOBL(2)='EIZZ'
  206.             LESOBL(3)='EITT'
  207.           ENDIF
  208.         ENDIF
  209.  
  210. C  5 - Elements COQUES
  211. C =====================
  212.       ELSE IF (MFR.EQ.3.OR.MFR.EQ.39) THEN
  213.         IF (NPINT.EQ.0) THEN
  214. C =====
  215. C  5.1 - Tridimensionnel et Fourier
  216. C =====
  217.           IF (IFOU.EQ.2.OR.IFOU.EQ.1) THEN
  218.             NBROBL=6
  219.             SEGINI,NOMID
  220.             LESOBL(1)='EISS'
  221.             LESOBL(2)='EITT'
  222.             LESOBL(3)='GIST'
  223.             LESOBL(4)='RISS'
  224.             LESOBL(5)='RITT'
  225.             LESOBL(6)='RIST'
  226. C =====
  227. C  5.2 - Axisymetrie
  228. C =====
  229.           ELSE IF (IFOU.EQ.0) THEN
  230.             NBROBL=4
  231.             SEGINI,NOMID
  232.             LESOBL(1)='EISS'
  233.             LESOBL(2)='EITT'
  234.             LESOBL(3)='RISS'
  235.             LESOBL(4)='RITT'
  236. C =====
  237. C  5.3 - Bidimensionnel PLAN (CP/DP/DPGE)
  238. C =====
  239.           ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
  240.             NBROBL=4
  241.             SEGINI,NOMID
  242.             LESOBL(1)='EISS'
  243.             LESOBL(2)='EIZZ'
  244.             LESOBL(3)='RISS'
  245.             LESOBL(4)='RIZZ'
  246.           ENDIF
  247.         ELSE
  248. C =====
  249. C  5.4 - Tridimensionnel - Cas particulier
  250. C =====
  251.           IF (MFR.EQ.3.AND.IFOU.EQ.2) THEN
  252.             NBROBL=4
  253.             SEGINI,NOMID
  254.             LESOBL(1)='EISS'
  255.             LESOBL(2)='EITT'
  256.             LESOBL(3)='EINN'
  257.             LESOBL(4)='GIST'
  258.           ENDIF
  259.         ENDIF
  260.  
  261. C  6 - Elements de COQUE EPAISSE
  262. C      Tridimensionnel, axisymetrie et Fourier
  263. C =============================================
  264.       ELSE IF (MFR.EQ.5) THEN
  265.         NBROBL=5
  266.         SEGINI,NOMID
  267.         LESOBL(1)='EISS'
  268.         LESOBL(2)='EITT'
  269.         LESOBL(3)='GIST'
  270.         LESOBL(4)='GISN'
  271.         LESOBL(5)='GITN'
  272.  
  273. C  7 - LIA2 : Element 3D de liaison a 2 noeuds
  274. C =============================================
  275.       ELSE IF (MFR.EQ.51) THEN
  276.         NBROBL=6
  277.         SEGINI,NOMID
  278.         LESOBL(1)='FEPS'
  279.         LESOBL(2)='FGXY'
  280.         LESOBL(3)='FGXZ'
  281.         LESOBL(4)='FFCX'
  282.         LESOBL(5)='FFCY'
  283.         LESOBL(6)='FFCZ'
  284. C  7bis - JOI1: Element 3D de liaison a 2 noeuds
  285. C =============================================
  286.       ELSE IF (MFR.EQ.75) THEN
  287.        IF (IFOU.EQ.2) THEN
  288.         NBROBL=6
  289.         SEGINI,NOMID
  290.         LESOBL(1)='FEXX '
  291.         LESOBL(2)='FGXY '
  292.         LESOBL(3)='FGXZ '
  293.         LESOBL(4)='FCXX  '
  294.         LESOBL(5)='FCXY  '
  295.         LESOBL(6)='FCXZ  '
  296.        ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  297.         NBROBL=3
  298.         SEGINI,NOMID
  299.         LESOBL(1)='FEXX '
  300.         LESOBL(2)='FGXY '
  301.         LESOBL(3)='FCXZ '
  302.        ENDIF
  303.  
  304. C  8 - BAEX : BARRE EXCENTREE
  305. C ============================
  306.       ELSE IF (MFR.EQ.49) THEN
  307.         NBROBL=1
  308.         SEGINI,NOMID
  309.         LESOBL(1)='EISS'
  310.  
  311. C  9 - Elements POUTRE et TUYAU
  312. C ==============================
  313.       ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN
  314. C =====
  315. C  9.1 - Tridimensionnel
  316. C =====
  317.         IF (IFOU.EQ.2) THEN
  318.           NBROBL=6
  319.           SEGINI,NOMID
  320.           LESOBL(1)='FEPS'
  321.           LESOBL(2)='FGXY'
  322.           LESOBL(3)='FGXZ'
  323.           LESOBL(4)='FFCX'
  324.           LESOBL(5)='FFCY'
  325.           LESOBL(6)='FFCZ'
  326. C =====
  327. C  9.2 - Bidimensionnel PLAN (CP/DP/DPGE)
  328. C =====
  329.         ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  330.           NBROBL=3
  331.           SEGINI,NOMID
  332.           LESOBL(1)='FEPS'
  333.           LESOBL(2)='FGXY'
  334.           LESOBL(3)='FFCZ'
  335.         ENDIF
  336.  
  337. C  10 - Element de COQUE en CISAILLEMENT TRANSVERSE
  338. C ==================================================
  339.       ELSE IF (MFR.EQ.9) THEN
  340.         IF (NPINT.EQ.0) THEN
  341. C ======
  342. C  10.1 - Tridimensionnel
  343. C ======
  344.           IF (IFOU.EQ.2) THEN
  345.             NBROBL=8
  346.             SEGINI,NOMID
  347.             LESOBL(1)='EISS'
  348.             LESOBL(2)='EITT'
  349.             LESOBL(3)='GIST'
  350.             LESOBL(4)='RISS'
  351.             LESOBL(5)='RITT'
  352.             LESOBL(6)='RIST'
  353.             LESOBL(7)='GISN'
  354.             LESOBL(8)='GITN'
  355.           ENDIF
  356.         ENDIF
  357.  
  358. C  11 - Elements LINESPRING (tridimensionnel)
  359. C ============================================
  360.       ELSE IF (MFR.EQ.15) THEN
  361.         IF (IFOU.EQ.2) THEN
  362.           NBROBL=6
  363.           SEGINI,NOMID
  364.           LESOBL(1)='EIZZ'
  365.           LESOBL(2)='GIXZ'
  366.           LESOBL(3)='GIYZ'
  367.           LESOBL(4)='RIXX'
  368.           LESOBL(5)='RIZZ'
  369.           LESOBL(6)='FDJP'
  370.         ENDIF
  371.  
  372. C  12 - Elements de MEMBRANE
  373. C ===========================
  374.       ELSE IF (MFR.EQ.25) THEN
  375.         NBROBL=3
  376.         SEGINI,NOMID
  377.         LESOBL(1)='EISS'
  378.         LESOBL(2)='EITT'
  379.         LESOBL(3)='GIST'
  380.  
  381. C  13 - Formulation UNIAXIALE
  382. C =============================
  383.       ELSE IF (MFR.EQ.27) THEN
  384.         NBROBL=1
  385.         SEGINI,NOMID
  386.         LESOBL(1)='EISS'
  387.  
  388. C  14 - Formulation THERMIQUE
  389. C ============================
  390.       ELSE IF (MFR.EQ.29) THEN
  391. C ======
  392. C  14.1 - Tridimensionnel
  393. C ======
  394.         IF (IFOU.EQ.2) THEN
  395.           NBROBL=3
  396.           SEGINI,NOMID
  397.           LESOBL(1)='FRDX'
  398.           LESOBL(2)='FRDY'
  399.           LESOBL(3)='FRDZ'
  400. C ======
  401. C  14.2 - Bidimensionnel PLAN (CP/DP/DPGE)
  402. C ======
  403.         ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
  404.           NBROBL=2
  405.           SEGINI,NOMID
  406.           LESOBL(1)='FRDX'
  407.           LESOBL(2)='FRDY'
  408. C ======
  409. C  14.3 - Axisymetrie et Fourier
  410. C ======
  411.         ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN
  412.           NBROBL=2
  413.           SEGINI,NOMID
  414.           LESOBL(1)='FRDX'
  415.           LESOBL(2)='FRDY'
  416. C*OF      LESOBL(1)='FRDR'
  417. C*OF      LESOBL(2)='FRDZ'
  418. C ======
  419. C  14.4 - Unidimensionnel (1D)
  420. C ======
  421.         ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN
  422.           NBROBL=1
  423.           SEGINI,NOMID
  424. C*OF      IF (IFOU.LE.11) THEN
  425.             LESOBL(1)='FRDX'
  426. C*OF      ELSE
  427. C*OF        LESOBL(1)='FRDR'
  428. C*OF      ENDIF
  429.         ENDIF
  430.  
  431. C  15 - Formulation POREUX
  432. C =========================
  433.       ELSE IF (MFR.EQ.33) THEN
  434. C =====
  435. C  15.A - Elements POREUX
  436. C ======
  437.         IF (MELE.GE.79.AND.MELE.LE.83) THEN
  438. C ========
  439. C  15.A.1 - Bidimensionnel PLAN (CP/DP)
  440. C ========
  441.           IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  442.             NBROBL=5
  443.             SEGINI,NOMID
  444.             LESOBL(1)='EIXX'
  445.             LESOBL(2)='EIYY'
  446.             LESOBL(3)='EIZZ'
  447.             LESOBL(4)='GIXY'
  448.             LESOBL(5)='PIDE'
  449. C ========
  450. C  15.A.2 - Axisymetrie
  451. C ========
  452.           ELSE IF (IFOU.EQ.0) THEN
  453.             NBROBL=5
  454.             SEGINI,NOMID
  455.             LESOBL(1)='EIRR'
  456.             LESOBL(2)='EIZZ'
  457.             LESOBL(3)='EITT'
  458.             LESOBL(4)='GIRZ'
  459.             LESOBL(5)='PIDE'
  460. C ========
  461. C  15.A.3 - Fourier
  462. C ========
  463.           ELSE IF (IFOU.EQ.1) THEN
  464.             NBROBL=7
  465.             SEGINI,NOMID
  466.             LESOBL(1)='EIRR'
  467.             LESOBL(2)='EIZZ'
  468.             LESOBL(3)='EITT'
  469.             LESOBL(4)='GIRZ'
  470.             LESOBL(5)='GIRT'
  471.             LESOBL(6)='GIZT'
  472.             LESOBL(7)='PIDE'
  473. C ========
  474. C  15.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  475. C ========
  476.           ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  477.             NBROBL=7
  478.             SEGINI,NOMID
  479.             LESOBL(1)='EIXX'
  480.             LESOBL(2)='EIYY'
  481.             LESOBL(3)='EIZZ'
  482.             LESOBL(4)='GIXY'
  483.             LESOBL(5)='GIXZ'
  484.             LESOBL(6)='GIYZ'
  485.             LESOBL(7)='PIDE'
  486.          ENDIF
  487. C ======
  488. C  15.B - Elements JOINT POREUX
  489. C ======
  490.         ELSE IF (MELE.GE.108.AND.MELE.LE.110) THEN
  491. C ========
  492. C  15.B.1 - Tridimensionnel
  493. C ========
  494.           IF (IFOU.EQ.2) THEN
  495.             NBROBL=4
  496.             SEGINI,NOMID
  497.             LESOBL(1)='FRS1'
  498.             LESOBL(2)='FRS2'
  499.             LESOBL(3)='FRN '
  500.             LESOBL(4)='PIDE'
  501. C ========
  502. C  15.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
  503. C ========
  504.           ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  505.             NBROBL=3
  506.             SEGINI,NOMID
  507.             LESOBL(1)='FRSN'
  508.             LESOBL(2)='FRN '
  509.             LESOBL(3)='PIDE'
  510.           ENDIF
  511.         ENDIF
  512.  
  513. C  16 - Formulation POREUX type Q
  514. C ================================
  515.       ELSE IF (MFR.EQ.57) THEN
  516. C =====
  517. C  16.A - Elements POREUX Q
  518. C ======
  519.         IF (MELE.GE.173.AND.MELE.LE.177) THEN
  520. C ========
  521. C  16.A.1 - Bidimensionnel PLAN (CP/DP)
  522. C ========
  523.           IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  524.             NBROBL=6
  525.             SEGINI,NOMID
  526.             LESOBL(1)='EIXX'
  527.             LESOBL(2)='EIYY'
  528.             LESOBL(3)='EIZZ'
  529.             LESOBL(4)='GIXY'
  530.             LESOBL(5)='PIDE'
  531.             LESOBL(6)='QIDE'
  532. C ========
  533. C  16.A.2 - Axisymetrie
  534. C ========
  535.           ELSE IF (IFOU.EQ.0) THEN
  536.             NBROBL=6
  537.             SEGINI,NOMID
  538.             LESOBL(1)='EIRR'
  539.             LESOBL(2)='EIZZ'
  540.             LESOBL(3)='EITT'
  541.             LESOBL(4)='GIRZ'
  542.             LESOBL(5)='PIDE'
  543.             LESOBL(6)='QIDE'
  544. C ========
  545. C  16.A.3 - Fourier
  546. C ========
  547.           ELSE IF (IFOU.EQ.1) THEN
  548.             NBROBL=8
  549.             SEGINI,NOMID
  550.             LESOBL(1)='EIRR'
  551.             LESOBL(2)='EIZZ'
  552.             LESOBL(3)='EITT'
  553.             LESOBL(4)='GIRZ'
  554.             LESOBL(5)='GIRT'
  555.             LESOBL(6)='GIZT'
  556.             LESOBL(7)='PIDE'
  557.             LESOBL(8)='QIDE'
  558. C ========
  559. C  16.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  560. C ========
  561.           ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  562.             NBROBL=8
  563.             SEGINI,NOMID
  564.             LESOBL(1)='EIXX'
  565.             LESOBL(2)='EIYY'
  566.             LESOBL(3)='EIZZ'
  567.             LESOBL(4)='GIXY'
  568.             LESOBL(5)='GIXZ'
  569.             LESOBL(6)='GIYZ'
  570.             LESOBL(7)='PIDE'
  571.             LESOBL(8)='QIDE'
  572.           ENDIF
  573. C ======
  574. C  16.B - Elements JOINT POREUX Q
  575. C ======
  576.         ELSE IF (MELE.GE.185.AND.MELE.LE.187) THEN
  577. C ========
  578. C  16.B.1 - Tridimensionnel
  579. C ========
  580.           IF (IFOU.EQ.2) THEN
  581.             NBROBL=5
  582.             SEGINI,NOMID
  583.             LESOBL(1)='FRS1'
  584.             LESOBL(2)='FRS2'
  585.             LESOBL(3)='FRN '
  586.             LESOBL(4)='PIDE'
  587.             LESOBL(5)='QIDE'
  588. C ========
  589. C  16.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
  590. C ========
  591.           ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  592.             NBROBL=4
  593.             SEGINI,NOMID
  594.             LESOBL(1)='FRSN'
  595.             LESOBL(2)='FRN '
  596.             LESOBL(3)='PIDE'
  597.             LESOBL(4)='QIDE'
  598.           ENDIF
  599.         ENDIF
  600.  
  601. C  17 - Formulation POREUX type R
  602. C ================================
  603.       ELSE IF (MFR.EQ.59) THEN
  604. C =====
  605. C  17.A - Elements POREUX R
  606. C ======
  607.         IF (MELE.GE.178.AND.MELE.LE.182) THEN
  608. C ========
  609. C  17.A.1 - Bidimensionnel PLAN (CP/DP)
  610. C ========
  611.           IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  612.             NBROBL=7
  613.             SEGINI,NOMID
  614.             LESOBL(1)='EIXX'
  615.             LESOBL(2)='EIYY'
  616.             LESOBL(3)='EIZZ'
  617.             LESOBL(4)='GIXY'
  618.             LESOBL(5)='PIDE'
  619.             LESOBL(6)='QIDE'
  620.             LESOBL(7)='TIDE'
  621. C ========
  622. C  17.A.2 - Axisymetrie
  623. C ========
  624.           ELSE IF (IFOU.EQ.0) THEN
  625.             NBROBL=7
  626.             SEGINI,NOMID
  627.             LESOBL(1)='EIRR'
  628.             LESOBL(2)='EIZZ'
  629.             LESOBL(3)='EITT'
  630.             LESOBL(4)='GIRZ'
  631.             LESOBL(5)='PIDE'
  632.             LESOBL(6)='QIDE'
  633.             LESOBL(7)='TIDE'
  634. C ========
  635. C  17.A.3 - Fourier
  636. C ========
  637.           ELSE IF (IFOU.EQ.1) THEN
  638.             NBROBL=9
  639.             SEGINI,NOMID
  640.             LESOBL(1)='EIRR'
  641.             LESOBL(2)='EIZZ'
  642.             LESOBL(3)='EITT'
  643.             LESOBL(4)='GIRZ'
  644.             LESOBL(5)='GIRT'
  645.             LESOBL(6)='GIZT'
  646.             LESOBL(7)='PIDE'
  647.             LESOBL(8)='QIDE'
  648.             LESOBL(9)='TIDE'
  649. C ========
  650. C  17.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  651. C ========
  652.           ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  653.             NBROBL=9
  654.             SEGINI,NOMID
  655.             LESOBL(1)='EIXX'
  656.             LESOBL(2)='EIYY'
  657.             LESOBL(3)='EIZZ'
  658.             LESOBL(4)='GIXY'
  659.             LESOBL(5)='GIXZ'
  660.             LESOBL(6)='GIYZ'
  661.             LESOBL(7)='PIDE'
  662.             LESOBL(8)='QIDE'
  663.             LESOBL(9)='TIDE'
  664.           ENDIF
  665. C ======
  666. C  17.B - Elements JOINT POREUX R
  667. C ======
  668.         ELSE IF (MELE.GE.188.AND.MELE.LE.190) THEN
  669. C ========
  670. C  17.B.1 - Tridimensionnel
  671. C ========
  672.           IF (IFOU.EQ.2) THEN
  673.             NBROBL=6
  674.             SEGINI,NOMID
  675.             LESOBL(1)='FRS1'
  676.             LESOBL(2)='FRS2'
  677.             LESOBL(3)='FRN '
  678.             LESOBL(4)='PIDE'
  679.             LESOBL(5)='QIDE'
  680.             LESOBL(6)='TIDE'
  681. C ========
  682. C  17.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
  683. C ========
  684.           ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  685.             NBROBL=5
  686.             SEGINI,NOMID
  687.             LESOBL(1)='FRSN'
  688.             LESOBL(2)='FRN '
  689.             LESOBL(3)='PIDE'
  690.             LESOBL(4)='QIDE'
  691.             LESOBL(5)='TIDE'
  692.           ENDIF
  693.         ENDIF
  694.  
  695. C  18 - Elements TUYAU FISSURE (tridimensionnel)
  696. C ===============================================
  697.       ELSE IF (MFR.EQ.17) THEN
  698.         IF (IFOU.EQ.2) THEN
  699.           NBROBL=8
  700.           SEGINI,NOMID
  701.           LESOBL(1)='FEPS'
  702.           LESOBL(2)='FGXY'
  703.           LESOBL(3)='FGXZ'
  704.           LESOBL(4)='FFCX'
  705.           LESOBL(5)='FFCY'
  706.           LESOBL(6)='FFCZ'
  707.           LESOBL(7)='FPS7'
  708.           LESOBL(8)='FPS8'
  709.         ENDIF
  710.  
  711. C  19 - Elements de SECTION/FIBRE
  712. C ================================
  713.       ELSE IF (MFR.EQ.47) THEN
  714. C ======
  715. C  19.1 - Bidimensionnel PLAN (CP/DP/DPGE)
  716. C ======
  717.         IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
  718.           NBROBL=2
  719.           SEGINI,NOMID
  720.           LESOBL(1)='EIXX'
  721.           LESOBL(2)='GIXY'
  722. C*OF      LESOBL(3)='GIXZ'
  723. C ======
  724. C  19.2 - Tridimensionnel
  725. C ======
  726.         ELSE IF (IFOU.EQ.2) THEN
  727.           NBROBL=3
  728.           SEGINI,NOMID
  729.           LESOBL(1)='EIXX'
  730.           LESOBL(2)='GIXY'
  731.           LESOBL(3)='GIXZ'
  732.         ENDIF
  733.  
  734. C  20 - Elements de zones cohesives
  735. C ==================================
  736.       ELSE IF (MFR.EQ.77) THEN
  737. C =====
  738. C  20.1 - Tridimensionnel
  739. C =====
  740.         IF (IFOU.GT.0) THEN
  741.           NBROBL=3
  742.           SEGINI,NOMID
  743.           LESOBL(1)='FRS1'
  744.           LESOBL(2)='FRS2'
  745.           LESOBL(3)='FRN '
  746. C =====
  747. C  20.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie
  748. C =====
  749.         ELSE IF (IFOU.LE.0) THEN
  750.           NBROBL=2
  751.           SEGINI,NOMID
  752.           LESOBL(1)='FRSN'
  753.           LESOBL(2)='FRN '
  754.         ENDIF
  755. c
  756. c    element coaxial COS2 (3D pour liaison acier-beton)
  757. c
  758.       ELSEIF ( MFR.EQ.78) then
  759.          NBROBL= 3
  760.          SEGINI NOMID
  761.          LESOBL(1)='GINT'
  762.          LESOBL(2)='GINN'
  763.          LESOBL(3)='GINM'
  764.       ENDIF
  765. C ===========================
  766. C  Par DEFAUT : segment VIDE
  767. C ===========================
  768.  100  continue
  769.       IF (NOMID.EQ.0) SEGINI,NOMID
  770.  
  771.       SEGACT,NOMID*NOMOD
  772.       IPNOMC=NOMID
  773.  
  774.       RETURN
  775.       END
  776.  
  777.  
  778.  

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