Télécharger dktshp.eso

Retour à la liste

Numérotation des lignes :

  1. C DKTSHP SOURCE CHAT 05/02/25 21:15:12 5035
  2. SUBROUTINE DKTSHP(IGAU,XEL,SHP,DJAC)
  3. C=======================================================================
  4. C EBERSOLT OCT 85 FONCTIONS DE FORME POUR W DU DKT
  5. C ENTREES
  6. C IGAU = NUMERO DU POINT DE GAUSS
  7. C XEL(3,3) = COORDONNEES LOCALES DE L ELEMENT
  8. C SORTIES
  9. C SHP(6,9) = W W,X W,Y EN CE POINT DE GAUSS
  10. C DJAC = JACOBIEN AU POINT DE GAUSS
  11. C NE PAS CONFONDRE AVEC BFDKT DONNANT THET X,X THET Y,X
  12. C THET Y,X THET Y,Y
  13. C
  14. C======================================================================
  15. IMPLICIT INTEGER(I-N)
  16. IMPLICIT REAL*8(A-H,O-Z)
  17. DIMENSION SHP(6,*),XEL(3,*)
  18. DATA UNDEMI,UNQUA,TRDEMI/.5D0,.25D0,1.5D0/
  19. DATA UNHUIT,TRQUA/.125D0,.75D0/
  20. C
  21. CALL ZERO(SHP,6,9)
  22. C
  23. X21=XEL(1,2)-XEL(1,1)
  24. Y31=XEL(2,3)-XEL(2,1)
  25. DJAC=Y31*X21
  26. C
  27. GOTO (100,200,300) ,IGAU
  28. 100 CONTINUE
  29. DX=XEL(1,2)-XEL(1,1)
  30. DY=XEL(2,2)-XEL(2,1)
  31. IDEB=0
  32. IFIN=3
  33. GOTO 500
  34. 200 CONTINUE
  35. DX=XEL(1,3)-XEL(1,1)
  36. DY=XEL(2,3)-XEL(2,1)
  37. IDEB=0
  38. IFIN=6
  39. GOTO 500
  40. 300 CONTINUE
  41. DX=XEL(1,3)-XEL(1,2)
  42. DY=XEL(2,3)-XEL(2,2)
  43. IDEB=3
  44. IFIN=6
  45. GOTO 500
  46. 500 CONTINUE
  47. XLONG= DX*DX + DY*DY
  48. XLON2= SQRT(XLONG)
  49. CC=DX/XLON2
  50. SS=DY/XLON2
  51. C
  52. C
  53. SHP(1,IDEB+1)= UNDEMI
  54. SHP(2,IDEB+1)=-TRDEMI*CC/XLON2
  55. SHP(3,IDEB+1)=-TRDEMI*SS/XLON2
  56. C
  57. SHP(1,IDEB+2)=-UNHUIT*XLON2*CC
  58. SHP(2,IDEB+2)= -TRQUA*CC*SS
  59. SHP(3,IDEB+2)= - SS*SS*UNQUA+CC*CC*UNDEMI
  60. C
  61. SHP(1,IDEB+3)=-UNHUIT*XLON2*SS
  62. SHP(2,IDEB+3)= CC*CC*UNQUA-SS*SS*UNDEMI
  63. SHP(3,IDEB+3)= TRQUA*CC*SS
  64. C
  65.  
  66. DO 510 IA=1,3
  67. DO 510 IB=1,3
  68. SHP(IA,IFIN+IB)=SHP(IA,IDEB+IB)
  69. 510 CONTINUE
  70. C
  71. SHP(2,IFIN+1)=-SHP(2,IDEB+1)
  72. SHP(3,IFIN+1)=-SHP(3,IDEB+1)
  73. C
  74.  
  75. SHP(1,IFIN+2)=-SHP(1,IDEB+2)
  76. SHP(1,IFIN+3)=-SHP(1,IDEB+3)
  77. C
  78. RETURN
  79. END
  80.  
  81.  

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