Télécharger la plaquette Cast3M


fr en

Télécharger idvar1.eso

Retour à la liste

Numérotation des lignes : 

idvar1

  1. C IDVAR1    SOURCE    FANDEUR   10/02/01    21:15:36     6603
  2.  
  3. C=======================================================================
  4. *       RECHERCHE DES NOMS DE COMPOSANTES DES VARIABLES INTERNES
  5. *       --------------------------------------------------------
  6. *
  7. *  ENTREES :
  8. *  ---------
  9. *     MATEPL  NUMERO DU MATERIAU PLASTIQUE
  10. *     MFR     NUMERO DE FORMULATION
  11. *     IFOUR   VALEUR DE IFOU DANS CCOPTIO
  12. *     NPINT   NOMBRE DE POINTS D'INTEGRATION DANS L'EPAISSEUR
  13. *             POUR LES ELEMENTS DE COQUE INTEGRES
  14. *
  15. *  SORTIES :
  16. *  ---------
  17. *     IPCOMP  POINTEUR SUR LES LISTES DES COMPOSANTES
  18. *             OBLIGATOIRES ET FACULTATIFS
  19. *     NBROBL  NOMBRE DES COMPOSANTES OBLIGATOIRES
  20. *     NBRFAC  NOMBRE DES COMPOSANTES FACULTATIF
  21. C=======================================================================
  22.  
  23.       SUBROUTINE IDVAR1 (MATEPL,MFR,IFOUR,NPINT,IPCOMP,NBROBL,NBRFAC)
  24.  
  25.       IMPLICIT INTEGER(I-N)
  26.       IMPLICIT REAL*8 (A-H,O-Z)
  27.  
  28.       SEGMENT NOMID
  29.         CHARACTER*8 LESOBL(NBROBL),LESFAC(NBRFAC)
  30.       ENDSEGMENT
  31.  
  32.       NBROBL=0
  33.       NBRFAC=0
  34.       NOMID=0
  35.  
  36. C  MODELES DE FLUAGE : NORTON (MATEPL=19)  et  POLYNOMIAL (MATEPL=21)
  37. C ====================================================================
  38.  
  39. *  ELEMENT TUYAU FISSURE
  40. *************************
  41.       IF (MFR.EQ.17) THEN
  42.         NBROBL=7
  43.         SEGINI,NOMID
  44.         LESOBL(1)='EPLU '
  45.         LESOBL(2)='EMIN'
  46.         LESOBL(3)='EEFF '
  47.         LESOBL(4)='EPSH'
  48.         LESOBL(5)='LL  '
  49.         LESOBL(6)='EPSE'
  50.         LESOBL(7)='JPOI'
  51.  
  52. *  ELEMENT MASSIF
  53. ******************
  54.       ELSE IF (MFR.EQ.1.OR.MFR.EQ.31.OR.MFR.EQ.33) THEN
  55. * CAS CONTRAINTES PLANES OU DEFORMATIONS PLANES
  56. * OU DEFORMATIONS PLANES GENERALISEES
  57.         IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
  58.           NBROBL=12
  59.           SEGINI,NOMID
  60.           LESOBL( 1)='PXX '
  61.           LESOBL( 2)='PYY '
  62.           LESOBL( 3)='PZZ '
  63.           LESOBL( 4)='PXY '
  64.           LESOBL( 5)='VNXX'
  65.           LESOBL( 6)='VNYY'
  66.           LESOBL( 7)='VNZZ'
  67.           LESOBL( 8)='VNXY'
  68.           LESOBL( 9)='EPCE'
  69.           LESOBL(10)='EPSH'
  70.           LESOBL(11)='LL  '
  71.           LESOBL(12)='EPSE'
  72. * CAS AXISYMETRIQUE
  73.         ELSE IF (IFOUR.EQ.0)  THEN
  74.           NBROBL=12
  75.           SEGINI,NOMID
  76.           LESOBL( 1)='PRR '
  77.           LESOBL( 2)='PZZ '
  78.           LESOBL( 3)='PTT '
  79.           LESOBL( 4)='PRZ '
  80.           LESOBL( 5)='VNRR'
  81.           LESOBL( 6)='VNZZ'
  82.           LESOBL( 7)='VNTT'
  83.           LESOBL( 8)='VNRZ'
  84.           LESOBL( 9)='EPCE'
  85.           LESOBL(10)='EPSH'
  86.           LESOBL(11)='LL  '
  87.           LESOBL(12)='EPSE'
  88. * CAS SERIE DE FOURIER
  89.         ELSE IF (IFOUR.EQ.1) THEN
  90.           NBROBL=16
  91.           SEGINI,NOMID
  92.           LESOBL( 1)='PRR '
  93.           LESOBL( 2)='PZZ '
  94.           LESOBL( 3)='PTT '
  95.           LESOBL( 4)='PRZ '
  96.           LESOBL( 5)='PRT '
  97.           LESOBL( 6)='PZT '
  98.           LESOBL( 7)='VNRR'
  99.           LESOBL( 8)='VNZZ'
  100.           LESOBL( 9)='VNTT'
  101.           LESOBL(10)='VNRZ'
  102.           LESOBL(11)='VNRT'
  103.           LESOBL(12)='VNZT'
  104.           LESOBL(13)='EPCE'
  105.           LESOBL(14)='EPSH'
  106.           LESOBL(15)='LL  '
  107.           LESOBL(16)='EPSE'
  108. * CAS TRIDIMENSIONNEL
  109.         ELSE IF (IFOUR.EQ.2) THEN
  110.           NBROBL=16
  111.           SEGINI,NOMID
  112.           LESOBL( 1)='PXX '
  113.           LESOBL( 2)='PYY '
  114.           LESOBL( 3)='PZZ '
  115.           LESOBL( 4)='PXY '
  116.           LESOBL( 5)='PYZ '
  117.           LESOBL( 6)='PZX '
  118.           LESOBL( 7)='VNXX'
  119.           LESOBL( 8)='VNYY'
  120.           LESOBL( 9)='VNZZ'
  121.           LESOBL(10)='VNXY'
  122.           LESOBL(11)='VNYZ'
  123.           LESOBL(12)='VNZX'
  124.           LESOBL(13)='EPCE'
  125.           LESOBL(14)='EPSH'
  126.           LESOBL(15)='LL  '
  127.           LESOBL(16)='EPSE'
  128. * CAS UNIDIMENSIONNEL (1D) PLAN
  129.         ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN
  130.           NBROBL=10
  131.           SEGINI,NOMID
  132.           LESOBL( 1)='PXX '
  133.           LESOBL( 2)='PYY '
  134.           LESOBL( 3)='PZZ '
  135.           LESOBL( 4)='VNXX'
  136.           LESOBL( 5)='VNYY'
  137.           LESOBL( 6)='VNZZ'
  138.           LESOBL( 7)='EPCE'
  139.           LESOBL( 8)='EPSH'
  140.           LESOBL( 9)='LL  '
  141.           LESOBL(10)='EPSE'
  142. * CAS UNIDIMENSIONNEL (1D) AXIS ou SPHE
  143.         ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN
  144.           NBROBL=10
  145.           SEGINI,NOMID
  146.           LESOBL( 1)='PRR '
  147.           LESOBL( 2)='PZZ '
  148.           LESOBL( 3)='PTT '
  149.           LESOBL( 4)='VNRR'
  150.           LESOBL( 5)='VNZZ'
  151.           LESOBL( 6)='VNTT'
  152.           LESOBL( 7)='EPCE'
  153.           LESOBL( 8)='EPSH'
  154.           LESOBL( 9)='LL  '
  155.           LESOBL(10)='EPSE'
  156.         ENDIF
  157.  
  158. *  ELEMENT COQUE
  159. *****************
  160.       ELSE IF (MFR.EQ.3) THEN
  161.         IF (NPINT.EQ.0) THEN
  162. * CAS TRIDIMENSIONNEL ou FOURIER
  163.           IF (IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN
  164.             NBROBL=20
  165.             SEGINI,NOMID
  166.             LESOBL( 1)='PN11'
  167.             LESOBL( 2)='PN22'
  168.             LESOBL( 3)='PN33'
  169.             LESOBL( 4)='PN12'
  170.             LESOBL( 5)='PM11'
  171.             LESOBL( 6)='PM22'
  172.             LESOBL( 7)='PM33'
  173.             LESOBL( 8)='PM12'
  174.             LESOBL( 9)='NN11'
  175.             LESOBL(10)='NN22'
  176.             LESOBL(11)='NN33'
  177.             LESOBL(12)='NN12'
  178.             LESOBL(13)='NM11'
  179.             LESOBL(14)='NM22'
  180.             LESOBL(15)='NM33'
  181.             LESOBL(16)='NM12'
  182.             LESOBL(17)='EPCE'
  183.             LESOBL(18)='EPSH'
  184.             LESOBL(19)='LL  '
  185.             LESOBL(20)='EPSE'
  186. * CAS AXISYMETRIQUE
  187.           ELSE IF (IFOUR.EQ.0) THEN
  188.             NBROBL=16
  189.             SEGINI,NOMID
  190.             LESOBL( 1)='PN11'
  191.             LESOBL( 2)='PN22'
  192.             LESOBL( 3)='PN33'
  193.             LESOBL( 4)='PM11'
  194.             LESOBL( 5)='PM22'
  195.             LESOBL( 6)='PM33'
  196.             LESOBL( 7)='NN11'
  197.             LESOBL( 8)='NN22'
  198.             LESOBL( 9)='NN33'
  199.             LESOBL(10)='NM11'
  200.             LESOBL(11)='NM22'
  201.             LESOBL(12)='NM33'
  202.             LESOBL(13)='EPCE'
  203.             LESOBL(14)='EPSH'
  204.             LESOBL(15)='LL  '
  205.             LESOBL(16)='EPSE'
  206. * CAS CONTRAINTES PLANES ou DEFORMATIONS PLANES ou
  207. *     DEFORMATIONS PLANES GENERALISEES
  208.           ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-3) THEN
  209.             NBROBL=16
  210.             SEGINI,NOMID
  211.             LESOBL( 1)='PN11'
  212.             LESOBL( 2)='PN22'
  213.             LESOBL( 3)='PN33'
  214.             LESOBL( 4)='PM11'
  215.             LESOBL( 5)='PM22'
  216.             LESOBL( 6)='PM33'
  217.             LESOBL( 7)='NN11'
  218.             LESOBL( 8)='NN22'
  219.             LESOBL( 9)='NN33'
  220.             LESOBL(10)='NM11'
  221.             LESOBL(11)='NM22'
  222.             LESOBL(12)='NM33'
  223.             LESOBL(13)='EPCE'
  224.             LESOBL(14)='EPSH'
  225.             LESOBL(15)='LL  '
  226.             LESOBL(16)='EPSE'
  227.           ENDIF
  228.         ELSE
  229.           IF (IFOUR.EQ.2) THEN
  230.             NBROBL=12
  231.             SEGINI,NOMID
  232.             LESOBL( 1)='PMSS'
  233.             LESOBL( 2)='PMTT'
  234.             LESOBL( 3)='PMNN'
  235.             LESOBL( 4)='PMST'
  236.             LESOBL( 5)='NMSS'
  237.             LESOBL( 6)='NMTT '
  238.             LESOBL( 7)='NMNN '
  239.             LESOBL( 8)='NMST '
  240.             LESOBL( 9)='EPCE'
  241.             LESOBL(10)='EPSH'
  242.             LESOBL(11)='LL  '
  243.             LESOBL(12)='EPSE'
  244.           ENDIF
  245.         ENDIF
  246.  
  247. *  COQUE EPAISSE
  248. *****************
  249. *  CAS TRIDIMENSIONNEL
  250.       ELSE IF (MFR.EQ.5) THEN
  251.         IF (IFOUR.EQ.2) THEN
  252.           NBROBL=16
  253.           SEGINI,NOMID
  254.           LESOBL( 1)='PMSS'
  255.           LESOBL( 2)='PMTT'
  256.           LESOBL( 3)='PMNN'
  257.           LESOBL( 4)='PMST'
  258.           LESOBL( 5)='PMSN'
  259.           LESOBL( 6)='PMTN'
  260.           LESOBL( 7)='NMSS'
  261.           LESOBL( 8)='NMTT'
  262.           LESOBL( 9)='NMNN'
  263.           LESOBL(10)='NMST'
  264.           LESOBL(11)='NMSN'
  265.           LESOBL(12)='NMTN'
  266.           LESOBL(13)='EPCE'
  267.           LESOBL(14)='EPSH'
  268.           LESOBL(15)='LL  '
  269.           LESOBL(16)='EPSE'
  270.         ENDIF
  271.  
  272. *  POUTRE ET TUYAU
  273. *******************
  274. * CAS TRIDIMENSIONNEL
  275.       ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN
  276.         IF (IFOUR.EQ.2) THEN
  277.           NBROBL=16
  278.           SEGINI,NOMID
  279.           LESOBL( 1)='PEFX'
  280.           LESOBL( 2)='PEFY'
  281.           LESOBL( 3)='PEFZ'
  282.           LESOBL( 4)='PMOX'
  283.           LESOBL( 5)='PMOY'
  284.           LESOBL( 6)='PMOZ'
  285.           LESOBL( 7)='NEFX'
  286.           LESOBL( 8)='NEFY'
  287.           LESOBL( 9)='NEFZ'
  288.           LESOBL(10)='NMOX'
  289.           LESOBL(11)='EMOY'
  290.           LESOBL(12)='EMOZ'
  291.           LESOBL(13)='EPCE'
  292.           LESOBL(14)='EPSH'
  293.           LESOBL(15)='LL  '
  294.           LESOBL(16)='EPSE'
  295.         ENDIF
  296.  
  297. *  COQUE CISAILLEMENT TRANSVERSE
  298. *********************************
  299. * CAS TRIDIMENSIONNEL
  300.       ELSE IF (MFR.EQ.9) THEN
  301.         IF (NPINT.EQ.0) THEN
  302.           IF (IFOUR.EQ.2) THEN
  303.             NBROBL=20
  304.             SEGINI,NOMID
  305.             LESOBL( 1)='PN11'
  306.             LESOBL( 2)='PN22'
  307.             LESOBL( 3)='PN12'
  308.             LESOBL( 4)='PM11'
  309.             LESOBL( 5)='PM22'
  310.             LESOBL( 6)='PM12'
  311.             LESOBL( 7)='PV1 '
  312.             LESOBL( 8)='PV2 '
  313.             LESOBL( 9)='NN11'
  314.             LESOBL(10)='NN22'
  315.             LESOBL(11)='NN12'
  316.             LESOBL(12)='NM11'
  317.             LESOBL(13)='NM22'
  318.             LESOBL(14)='NM12'
  319.             LESOBL(15)='NV1 '
  320.             LESOBL(16)='NV2 '
  321.             LESOBL(17)='EPCE'
  322.             LESOBL(18)='EPSH'
  323.             LESOBL(19)='LL  '
  324.             LESOBL(20)='EPSE'
  325.           ENDIF
  326.         ENDIF
  327.  
  328. *  FORMULATION MEMBRANE
  329. ************************
  330.       ELSE IF (MFR.EQ.25) THEN
  331.         NBROBL=10
  332.         SEGINI,NOMID
  333.         LESOBL( 1)='PMSS'
  334.         LESOBL( 2)='PMTT'
  335.         LESOBL( 3)='PMST'
  336.         LESOBL( 4)='NMSS'
  337.         LESOBL( 5)='NMTT'
  338.         LESOBL( 6)='NMST'
  339.         LESOBL( 7)='EPCE'
  340.         LESOBL( 8)='EPSH'
  341.         LESOBL( 9)='LL  '
  342.         LESOBL(10)='EPSE'
  343.  
  344. *  FORMULATION UNIAXIALE
  345. *************************
  346.       ELSE IF (MFR.EQ.27) THEN
  347.         NBROBL=6
  348.         SEGINI,NOMID
  349.         LESOBL(1)='PEFX'
  350.         LESOBL(2)='NEFX'
  351.         LESOBL(3)='EPCE'
  352.         LESOBL(4)='EPSH'
  353.         LESOBL(5)='LL  '
  354.         LESOBL(6)='EPSE'
  355.       ENDIF
  356.  
  357.       IPCOMP=NOMID
  358.       SEGDES,NOMID
  359.  
  360.       RETURN
  361.       END
  362.  
  363.  
  364.  
  365.  
  366.  

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