Télécharger la plaquette Cast3M


fr en

Télécharger ella07.eso

Retour à la liste

Numérotation des lignes : 

ella07

  1. C ELLA07    SOURCE    CHAT      05/01/12    23:33:48     5004
  2.       SUBROUTINE ELLA07(XA,YA,ZA,XB,YB,ZB,X1,Y1,Z1,X2,Y2,Z2,IOUI)
  3. C
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8 (A-H,O-Z)
  6. C
  7. C   OPERATEUR ELFE LAPLACE ACOU
  8. C
  9. C   REGARDE SI LE SEGMENT <1-2> APPARTIENT AU SEGMENT <AB>
  10. C
  11. C   PARAMETRES :
  12. C
  13. C   A ( XA,YA,ZA )   B ( XB,YB,ZB )   1 ( X1,Y1,Z1 )   2 ( X2,Y2,Z2 )
  14. C
  15. C
  16. C   SORTIES :
  17. C
  18. C   IOUI    : = 1 SI LA CONDITION EST VERIFIEE
  19. C             = 0 SI LA CONDITION N'EST PAS VERIFIEE
  20. C
  21. C   AUTEURS : SAINT-DIZIER ET GORCY
  22. C   DATE    : 24 JANVIER 1991
  23. C
  24.       IOUI = 0
  25. C
  26. C -- VECTEUR AB ET LONGUEUR DE LE POUTRE
  27. C
  28.       XBA = XB - XA
  29.       YBA = YB - YA
  30.       ZBA = ZB - ZA
  31.       XL = SQRT ( XBA*XBA + YBA*YBA + ZBA*ZBA)
  32.       EPS0 = 1.D-15
  33.       IF (XL.LT.EPS0) RETURN
  34. C
  35. C-- VECTEURS AM1, AM2, BM1, BM2
  36. C
  37.       X1A = X1 - XA
  38.       Y1A = Y1 - YA
  39.       Z1A = Z1 - ZA
  40.       XL1A = SQRT ( X1A*X1A + Y1A*Y1A + Z1A*Z1A )
  41. C
  42.       X2A = X2 - XA
  43.       Y2A = Y2 - YA
  44.       Z2A = Z2 - ZA
  45.       XL2A = SQRT ( X2A*X2A + Y2A*Y2A + Z2A*Z2A )
  46. C
  47.       X1B = X1 - XB
  48.       Y1B = Y1 - YB
  49.       Z1B = Z1 - ZB
  50.       XL1B = SQRT ( X1B*X1B + Y1B*Y1B + Z1B*Z1B )
  51. C
  52.       X2B = X2 - XB
  53.       Y2B = Y2 - YB
  54.       Z2B = Z2 - ZB
  55.       XL2B = SQRT ( X2B*X2B + Y2B*Y2B + Z2B*Z2B )
  56. C
  57. C -- PRODUITS VECTORIELS AM1^AB ET AM2^AB ET NORMES ASSOCIEES
  58. C
  59.       XVECT1 = ZBA*Y1A - YBA*Z1A
  60.       YVECT1 = XBA*Z1A - ZBA*X1A
  61.       ZVECT1 = YBA*X1A - XBA*Y1A
  62.       XNORM1 = SQRT ( XVECT1*XVECT1 + YVECT1*YVECT1 + ZVECT1*ZVECT1 )
  63. C
  64.       XVECT2 = ZBA*Y2A - YBA*Z2A
  65.       YVECT2 = XBA*Z2A - ZBA*X2A
  66.       ZVECT2 = YBA*X2A - XBA*Y2A
  67.       XNORM2 = SQRT ( XVECT2*XVECT2 + YVECT2*YVECT2 + ZVECT2*ZVECT2 )
  68. C
  69. C -- PRODUITS SCALAIRES AM1.BM1 ET AM2.BM2
  70. C
  71.       XP1 = X1A*X1B + Y1A*Y1B + Z1A*Z1B
  72.       XP2 = X2A*X2B + Y2A*Y2B + Z2A*Z2B
  73. C
  74.         IF (((XL1A.LT.EPS0).AND.(XP2.GT.0.)).OR.
  75.      *      ((XL1B.LT.EPS0).AND.(XP2.GT.0.))) RETURN
  76. C
  77.         IF (((XL2A.LT.EPS0).AND.(XP1.GT.0.D0)).OR.
  78.      *      ((XL2B.LT.EPS0).AND.(XP1.GT.0.D0))) RETURN
  79. C
  80.         IF (((XL1A.LT.EPS0).AND.(XP2.LT.0.D0)).OR.
  81.      *      ((XL1B.LT.EPS0).AND.(XP2.LT.0.D0))) THEN
  82.           IOUI = 1
  83.           RETURN
  84.         ENDIF
  85. C
  86.         IF (((XL2A.LT.EPS0).AND.(XP1.LT.0.D0)).OR.
  87.      *      ((XL2B.LT.EPS0).AND.(XP1.LT.0.D0))) THEN
  88.           IOUI = 1
  89.           RETURN
  90.         ENDIF
  91. C
  92. C -- NORMALISATION (REND LES CRITERES INDEPENDANTS DES DIMENSIONS
  93. C     DE LA POUTRE)
  94. C
  95.       XNORM1 = XNORM1 / (XL*XL1A)
  96.       XNORM2 = XNORM2 / (XL*XL2A)
  97.       XP1 = XP1 / (XL*XL)
  98.       XP2 = XP2 / (XL*XL)
  99. C
  100. C -- CRITERES GEOMETRIQUES
  101. C
  102.       EPS1 = 8.7D-3
  103.       EPS2 = 1.D-3
  104. C
  105. C -- IOUI VAUT 1 SI LES VECTEURS SONT COLINEAIRES (M1 ET M2 APPARTIENNEN
  106. C -- A LA DROITE (AB)), ET SI LES PRODUITS SCALAIRES SONT <0 (APPARTE-
  107. C -- NANCE AU SEGMENT AB)
  108. C
  109.       IF (XNORM1.LT.EPS1)THEN
  110.          IF (XNORM2.LT.EPS1) THEN
  111.             IF (XP1.LT.EPS2) THEN
  112.                IF (XP2.LT.EPS2) IOUI=1
  113.             ENDIF
  114.          ENDIF
  115.       ENDIF
  116. C
  117.       END
  118.  
  119.  
  120.  

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