Télécharger la plaquette Cast3M


fr en

Télécharger facref.eso

Retour à la liste

Numérotation des lignes : 

facref

  1. C FACREF    SOURCE    CHAT      05/01/12    23:55:49     5004
  2.       SUBROUTINE FACREF(X,ITYP,IFAC)
  3. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  4. C
  5. C      RECHERCHE IFAC NUMERO DE LA FACE OU SE TROUVE LE POINT DE
  6. C       COORDONNEES X DANS L'ELEMENT DE REFERENCE
  7. C
  8. C       ITYP TYPE DE L ELEMENT
  9. C
  10. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  11.       IMPLICIT INTEGER(I-N)
  12.       IMPLICIT REAL*8 (A-H,O-Z)
  13.       DIMENSION X(*)
  14.        IFAC=0
  15.       IF(ITYP.EQ.4.OR.ITYP.EQ.6.OR.ITYP.EQ.7)THEN
  16. C               L ELEMENT DE REFERENCE EST UN TRIANGLE
  17.            IF(X(2).EQ.0.D0)THEN
  18.                IFAC=1
  19.            ELSEIF((X(1)+X(2)).EQ.1.D0)THEN
  20.                IFAC=2
  21.            ELSEIF(X(1).EQ.0.D0)THEN
  22.                IFAC=3
  23.            ENDIF
  24.       ELSEIF(ITYP.EQ.8.OR.ITYP.EQ.11)THEN
  25. C              L ELEMENT DE REFERENCE EST UN CARRE
  26.            IF(X(2).EQ.-1.D0)THEN
  27.                IFAC=1
  28.            ELSEIF(X(1).EQ.1.D0)THEN
  29.                IFAC=2
  30.            ELSEIF(X(2).EQ.1.D0)THEN
  31.                IFAC=3
  32.            ELSEIF(X(1).EQ.-1.D0)THEN
  33.                IFAC=4
  34.            ENDIF
  35.       ELSEIF(ITYP.EQ.14)THEN
  36. C              L ELEMENT DE REFERENCE EST UN CUBE
  37.            IF(X(3).EQ.-1.D0)THEN
  38.                IFAC=1
  39.            ELSEIF(X(3).EQ.1.D0)THEN
  40.                IFAC=2
  41.            ELSEIF(X(2).EQ.-1.D0)THEN
  42.                IFAC=3
  43.            ELSEIF(X(1).EQ.1.D0)THEN
  44.                IFAC=4
  45.            ELSEIF(X(2).EQ.1.D0)THEN
  46.                IFAC=5
  47.            ELSEIF(X(1).EQ.-1.D0)THEN
  48.                IFAC=3
  49.            ENDIF
  50.       ELSEIF(ITYP.EQ.16)THEN
  51. C              L ELEMENT DE REFERENCE EST UN PRISME
  52.            IF(X(3).EQ.-1.D0)THEN
  53.                IFAC=1
  54.            ELSEIF(X(3).EQ.1.D0)THEN
  55.                IFAC=2
  56.            ELSEIF(X(2).EQ.0.D0)THEN
  57.                IFAC=3
  58.            ELSEIF((X(1)+X(2)).EQ.1.D0)THEN
  59.                IFAC=4
  60.            ELSEIF(X(1).EQ.0.D0)THEN
  61.                IFAC=5
  62.            ENDIF
  63.       ELSEIF(ITYP.EQ.23)THEN
  64. C              L ELEMENT DE REFERENCE EST TETRAEDRE
  65.            IF(X(2).EQ.0.D0)THEN
  66.                IFAC=1
  67.            ELSEIF(X(3).EQ.0.D0)THEN
  68.                IFAC=2
  69.            ELSEIF((X(1)+X(2)+X(3)).EQ.1.D0)THEN
  70.                IFAC=3
  71.            ELSEIF(X(1).EQ.0.D0)THEN
  72.                IFAC=4
  73.            ENDIF
  74.       ELSE
  75.            CALL ERREUR(16)
  76.       ENDIF
  77.       RETURN
  78.       END
  79.  
  80.  
  81.  

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