Télécharger kreper.eso

Retour à la liste

Numérotation des lignes :

  1. C KREPER SOURCE CB215821 16/04/21 21:17:37 8920
  2. SUBROUTINE KREPER(NES,NSOM,U1,A1,HCM)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. DIMENSION U1(NES),A1(NES,NSOM),HCM(NES,NES)
  6. DIMENSION R(3,3),T1(3),T2(3)
  7. C
  8. C Calcul des facteurs de forme en 3D
  9. C calcul de la matrice de passage du repere H.C au repere global
  10. C
  11. C NES : dim espace
  12. C NSOM : nombre de sommets
  13. C U1 : vecteur normal unitaire
  14. C A1 : sommets du triangle
  15. C HCM : matrice de passage
  16. C le repere local est orthonorme =(U1,A1.12,produit vectoriel)
  17. C
  18. DO 1 I=1,NES
  19. R(I,1)=U1(I)
  20. T1(I)=A1(I,2)-A1(I,1)
  21. 1 CONTINUE
  22.  
  23. S=0.D0
  24. DO 2 I=1,NES
  25. S = S+T1(I)*T1(I)
  26. 2 CONTINUE
  27. S = SQRT(S)
  28.  
  29. IF(S.LT.1D-10) THEN
  30. WRITE(6,*) ' >KREPER longueur cote 1 nulle'
  31. ELSE
  32. DO 3 I=1,NES
  33. T1(I)=T1(I)/S
  34. 3 CONTINUE
  35.  
  36. C produit vectoriel
  37. CALL KPVEC(U1,T1,T2)
  38.  
  39. DO 4 I=1,NES
  40. R(I,2)=T1(I)
  41. R(I,3)=T2(I)
  42. 4 CONTINUE
  43.  
  44. C matrice inverse
  45. CALL KPINV3(R,HCM)
  46.  
  47. C write(6,*) ' >KREPER : matrice R '
  48. C CALL UTPRIN(R,NES,NES)
  49. C write(6,*) ' >KREPER : matrice inverse '
  50. C CALL UTPRIN(HCM,NES,NES)
  51. ENDIF
  52.  
  53. RETURN
  54. END
  55.  
  56.  
  57.  
  58.  
  59.  

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