Télécharger geradk.eso

Retour à la liste

Numérotation des lignes :

  1. C GERADK SOURCE CHAT 05/01/13 00:17:55 5004
  2. SUBROUTINE GERADK(XEL,GEOM)
  3. C |====================================================================|
  4. C | |
  5. C | CALCULE DES CONSTANTES D UN TRIANGLE POUR L ELEMENT DKT |
  6. C | |
  7. C |==== ENTREES |
  8. C | XEL(3,3)= COORDONNEES DE L ELEMENT |
  9. C |==== SORTIES |
  10. C | GEOM(17)=DES CONSTANTES POUR LE DKT |
  11. C | CODE SUO X.Z. |
  12. C |====================================================================|
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8(A-H,O-Z)
  15. DIMENSION XEL(3,*),GEOM(*)
  16. DATA TRSQUA/0.75D0/,UNQUA/0.25D0/,UNDEMI/0.5D0/
  17. X12=XEL(1,1)-XEL(1,2)
  18. X23=XEL(1,2)-XEL(1,3)
  19. X31=XEL(1,3)-XEL(1,1)
  20. Y12=XEL(2,1)-XEL(2,2)
  21. Y23=XEL(2,2)-XEL(2,3)
  22. Y31=XEL(2,3)-XEL(2,1)
  23. XSQ12=X12*X12
  24. XSQ23=X23*X23
  25. XSQ31=X31*X31
  26. YSQ12=Y12*Y12
  27. YSQ23=Y23*Y23
  28. YSQ31=Y31*Y31
  29. DIST12=XSQ12+YSQ12
  30. DIST23=XSQ23+YSQ23
  31. DIST31=XSQ31+YSQ31
  32. A4=-X23/DIST23
  33. A5=-X31/DIST31
  34. A6=-X12/DIST12
  35. D4=-Y23/DIST23
  36. D5=-Y31/DIST31
  37. D6=-Y12/DIST12
  38. B4=TRSQUA*X23*Y23/DIST23
  39. B5=TRSQUA*X31*Y31/DIST31
  40. B6=TRSQUA*X12*Y12/DIST12
  41. C4=(UNQUA*XSQ23-UNDEMI*YSQ23)/DIST23
  42. C5=(UNQUA*XSQ31-UNDEMI*YSQ31)/DIST31
  43. C6=(UNQUA*XSQ12-UNDEMI*YSQ12)/DIST12
  44. E4=(UNQUA*YSQ23-UNDEMI*XSQ23)/DIST23
  45. E5=(UNQUA*YSQ31-UNDEMI*XSQ31)/DIST31
  46. E6=(UNQUA*YSQ12-UNDEMI*XSQ12)/DIST12
  47. GEOM(1)=A4
  48. GEOM(2)=A5
  49. GEOM(3)=A6
  50. GEOM(4)=B4
  51. GEOM(5)=B5
  52. GEOM(6)=B6
  53. GEOM(7)=C4
  54. GEOM(8)=C5
  55. GEOM(9)=C6
  56. GEOM(10)=D4
  57. GEOM(11)=D5
  58. GEOM(12)=D6
  59. GEOM(13)=E4
  60. GEOM(14)=E5
  61. GEOM(15)=E6
  62. GEOM(16)=1.D0
  63. GEOM(17)=1.D0
  64. RETURN
  65. END
  66.  
  67.  

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