1 : C TYPFAC SOURCE BP208322 10/07/12 21:15:14 6711 2 : SUBROUTINE TYPFAC(MELM,NBNN,IRET) 3 : C======================================================================= 4 : C 5 : C DONNE LA FORMULATION D UNE FACE DE ELEMENT MASSIF MELM 6 : C ENTREE 7 : C MELM = NUMERO DE L ELEMENT MASSIF DANS NOMTP 8 : C NBNN = NOMBRE DE NOEUDS DE LA FACE 9 : C SORTIE 10 : C IRET = NUMERO DE LA FACE DANS NOMTP 11 : C 0 SI OPERATION IMPOSSIBLE 12 : C 13 : C====================================================================== 14 : IMPLICIT INTEGER(I-N) 15 : IF (MELM.EQ.4.OR.MELM.EQ.5.OR.MELM.EQ.8.OR.MELM.EQ.9.OR. 16 : . MELM.EQ.69.OR.MELM.EQ.70.OR.MELM.EQ.183.or. 17 : . MELM.eq.263) THEN 18 : C 19 : C ELEMENTS MASSIFS BIDIMENSIONNELS LINEAIRES 20 : C 21 : IRET=2 22 : ELSE IF (MELM.GE.111.AND.MELM.LE.122) THEN 23 : C 24 : C ELEMENTS MASSIFS BIDIMENSIONNELS POLYGONAUX 25 : C 26 : IRET=2 27 : ELSE IF (MELM.EQ.6.OR.MELM.EQ.7.OR.MELM.EQ.10.OR.MELM.EQ.11. 28 : . OR.MELM.EQ.79.OR.MELM.EQ.80.OR.MELM.EQ.71.OR.MELM.EQ.72. 29 : . OR.MELM.EQ.184) THEN 30 : C 31 : C ELEMENTS MASSIFS BIDIMENSIONNELS QUADRATIQUES 32 : C 33 : IRET=3 34 : ELSE IF (MELM.EQ.14.OR.MELM.EQ.16.OR.MELM.EQ.23.OR.MELM.EQ.25.OR. 35 : . MELM.EQ.73.OR.MELM.EQ.264)THEN 36 : CTY . MELM.EQ.73)THEN 37 : C 38 : C ELEMENTS MASSIFS TRIDIMENSIONNELS LINEAIRES 39 : C 40 : IF (NBNN.EQ.3) THEN 41 : IRET=31 42 : ELSE 43 : IRET=32 44 : ENDIF 45 : ELSE IF (MELM.EQ.15.OR.MELM.EQ.17.OR.MELM.EQ.24.OR.MELM.EQ.26. 46 : . OR.MELM.EQ.81.OR.MELM.EQ.82.OR.MELM.EQ.83.OR.MELM.EQ.76) 47 : . THEN 48 : C 49 : C ELEMENTS MASSIFS TRIDIMENSIONNELS QUADRATIQUES 50 : C 51 : IF (NBNN.EQ.6) THEN 52 : IRET=33 53 : ELSE 54 : IRET=34 55 : ENDIF 56 : C 57 : C elements fluides 2d et 3d 58 : C 59 : else if (melm.eq.35 .or. melm.eq.36) then 60 : iret = 97 61 : else if (melm.eq.38. or. melm.eq.40) then 62 : if (nbnn.eq.3) then 63 : iret = 35 64 : else 65 : iret = 36 66 : endif 67 : else if (melm.eq.37) then 68 : iret = 36 69 : else if (melm.eq.39) then 70 : iret = 35 71 : 72 : C= Elements MASSIFs UNIDIMENSIONNELs (1D) 73 : ELSE IF (MELM.EQ.193.OR.MELM.EQ.194) THEN 74 : IRET=45 75 : NBNN=1 76 : 77 : ELSE 78 : IRET=0 79 : ENDIF 80 : RETURN 81 : END 82 : 83 : 84 : 85 : 86 : 87 : 88 : 89 : 90 : 91 :
© Cast3M 2003 - All rights reserved.
Disclaimer