Télécharger kchrep.eso

Retour à la liste

Numérotation des lignes :

  1. C KCHREP SOURCE CHAT 05/01/13 00:52:50 5004
  2. SUBROUTINE KCHREP(NES,NSOM,O,A,U,AH,UH,H)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. DIMENSION O(NES),A(NES,NSOM),U(NES+1)
  6. DIMENSION AH(NES,NSOM),UH(NES+1)
  7. DIMENSION H(NES,NES),B(3,3)
  8. C
  9. C Calcul des facteurs de forme en 3D
  10. C calcul des coordonnées d'un point A dans le repere local du HC
  11. C
  12. C NES : dim espace
  13. C NSOM : nombre de sommets
  14. C O : point de vue
  15. C A : sommets du triangle
  16. C U : vecteur normal unitaire
  17. C AH,UH : idem repere local du HC
  18. C H : matrice de passage durepere HC au repere global
  19. C
  20.  
  21. DO 2 L=1,NSOM
  22. DO 20 I=1,NES
  23. B(I,L)=A(I,L)-O(I)
  24. 20 CONTINUE
  25. 2 CONTINUE
  26.  
  27. DO 1 L=1,NSOM
  28. AH(1,L) = B(1,L)*H(1,1) + B(2,L)*H(1,2) + B(3,L)*H(1,3)
  29. AH(2,L) = B(1,L)*H(2,1) + B(2,L)*H(2,2) + B(3,L)*H(2,3)
  30. AH(3,L) = B(1,L)*H(3,1) + B(2,L)*H(3,2) + B(3,L)*H(3,3)
  31. 1 CONTINUE
  32.  
  33. UH(1) = U(1)*H(1,1) + U(2)*H(1,2) + U(3)*H(1,3)
  34. UH(2) = U(1)*H(2,1) + U(2)*H(2,2) + U(3)*H(2,3)
  35. UH(3) = U(1)*H(3,1) + U(2)*H(3,2) + U(3)*H(3,3)
  36.  
  37. UH(4) = -(UH(1)*AH(1,1) + UH(2)*AH(2,1) + UH(3)*AH(3,1))
  38.  
  39. RETURN
  40. END
  41.  
  42.  
  43.  

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