Télécharger bfdkt.eso

Retour à la liste

Numérotation des lignes :

  1. C BFDKT SOURCE CHAT 05/01/12 21:39:09 5004
  2. SUBROUTINE BFDKT(QSI,ETA,GEOM,B)
  3. C=======================================================================
  4. C
  5. C EVALUATION DE LA MATRICE B(3,9) AU POINT QSI ETA
  6. C POUR L'ELEMENT DKT
  7. C
  8. C CODE EBERSOLT JANVIER 84
  9. C=======================================================================
  10. IMPLICIT INTEGER(I-N)
  11. IMPLICIT REAL*8 (A-H,O-Z)
  12. DIMENSION BXP(9),BYP(9),BXE(9),BYE(9)
  13. DIMENSION GEOM(*),B(3,*)
  14. C
  15. DATA UN/1.0D0/,DEUX/2.0D0/,QUATRE/4.0D0/,SIX/6.0D0/
  16. C
  17. CL4=GEOM(1)
  18. CL5=GEOM(2)
  19. CL6=GEOM(3)
  20. SL4=GEOM(4)
  21. SL5=GEOM(5)
  22. SL6=GEOM(6)
  23. D4= GEOM(7)
  24. D5= GEOM(8)
  25. D6= GEOM(9)
  26. SU4=GEOM(10)
  27. SU5=GEOM(11)
  28. SU6=GEOM(12)
  29. VJ11=GEOM(13)
  30. VJ12=GEOM(14)
  31. VJ21=GEOM(15)
  32. VJ22=GEOM(16)
  33. DETJ=GEOM(17)
  34. PS2=UN-DEUX*QSI
  35. ET2=UN-DEUX*ETA
  36. SU5E=SU5*ET2
  37. C---- DERIVEE DE BETAX(I) /QSI
  38. BXP(1)=CL6*PS2+ETA*(CL5-CL6)
  39. BXP(2)=-ETA*D5
  40. BXP(3)=-QUATRE+SIX*(QSI+ETA)-ETA*SU5
  41. BXP(4)=-CL6*PS2+ETA*(CL4+CL6)
  42. BXP(5)=ETA*D4
  43. BXP(6)=-DEUX+SIX*QSI+ETA*SU4
  44. BXP(7)=-ETA*(CL5+CL4)
  45. BXP(8)=ETA*(D4-D5)
  46. BXP(9)=-ETA*(SU5-SU4)
  47. C---- DERIVEE DE BETAY(I) /QSI
  48. BYP(1)=ETA*SL5
  49. BYP(2)=UN-ETA*SU5
  50. BYP(3)=-BXP(2)
  51. BYP(4)=ETA*SL4
  52. BYP(5)=-UN+ETA*SU4
  53. BYP(6)=-BXP(5)
  54. BYP(7)=-ETA*(SL5+SL4)
  55. BYP(8)=BXP(9)
  56. BYP(9)=-BXP(8)
  57. C---- DERIVEE DE BETAX(I) /ETA
  58. BXE(1)=-CL5*ET2-QSI*(CL6-CL5)
  59. BXE(2)=D5*ET2-QSI*D5
  60. BXE(3)=-QUATRE+SIX*(QSI+ETA)+SU5E-QSI*SU5
  61. BXE(4)=QSI*(CL4+CL6)
  62. BXE(5)=QSI*D4
  63. BXE(6)=QSI*SU4
  64. BXE(7)=CL5*ET2-QSI*(CL4+CL5)
  65. BXE(8)=D5*ET2+QSI*(D4-D5)
  66. BXE(9)=-DEUX+SIX*ETA+SU5E+QSI*(SU4-SU5)
  67. C---- DERIVEE DE BETAY(I) /ETA
  68. BYE(1)=-SL5*ET2+QSI*SL5
  69. BYE(2)=UN+SU5E-QSI*SU5
  70. BYE(3)=-BXE(2)
  71. BYE(4)=QSI*SL4
  72. BYE(5)=BXE(6)
  73. BYE(6)=-BXE(5)
  74. BYE(7)=SL5*ET2-QSI*(SL4+SL5)
  75. BYE(8)=-UN+SU5E+QSI*(SU4-SU5)
  76. BYE(9)=-BXE(8)
  77. C---- CALCUL DE B
  78. DO 10 I=1,9
  79. B(1,I)=(VJ11*BXP(I)+VJ12*BXE(I))
  80. B(2,I)=(VJ21*BYP(I)+VJ22*BYE(I))
  81. B(3,I)=(VJ11*BYP(I)+VJ12*BYE(I)+VJ21*BXP(I)+VJ22*BXE(I))
  82. 10 CONTINUE
  83. C
  84. C
  85. RETURN
  86. END
  87.  
  88.  

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