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

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