Télécharger cukskp.eso

Retour à la liste

Numérotation des lignes :

  1. C CUKSKP SOURCE CHAT 07/10/22 21:15:40 5921
  2. SUBROUTINE CUKSKP(IFACE,BKSI)
  3. implicit real*8(A-H,O-Z)
  4. implicit integer (I-N)
  5.  
  6.  
  7. * SUBROUTINE CUB_KSI_KP(IFACE,BKSI)
  8. C ELEMENT SHB8
  9. C DEDIEE AU CALCUL DE KP: POUR LES 4 NOEUDS DE LA FACE PRESSION
  10. C
  11. C ------------------------------------------------------------------
  12. C
  13. C DERIVEE P/R VARIABLES CANONIQUE AUX POINTS D INTEGRATION
  14. C POUR LES ELEMENTS CUBES A 8 NOEUDS
  15. C H. BUNG 03-93
  16. C CORRIGE PAR A. LEGAY POUR 4 NOEUDS
  17. C ------------------------------------------------------------------
  18. C
  19. C
  20. C ENTREES
  21. C IFACE : =1 FACE INFERIEURE
  22. C =2 FACE SUPERIEURE
  23. C SORTIES :
  24. C BKSI(3,8) : LES DERIVEES
  25. C
  26. * IMPLICIT NONE
  27. C
  28. C VARIABLES GLOBALES
  29. C
  30. * INTEGER IFACE
  31. * REAL *8 BKSI(3,4)
  32. DIMENSION BKSI(3,4)
  33. C
  34. C VARIABLES LOCALES
  35. C
  36. * REAL *8 XPOINT(4,3),ZI,UNS8
  37. DIMENSION XPOINT(4,3)
  38. UNS8=0.125D0
  39. C
  40. XPOINT(1,1)=-1.
  41. XPOINT(1,2)=-1.
  42. XPOINT(2,1)=1.
  43. XPOINT(2,2)=-1.
  44. XPOINT(3,1)=1.
  45. XPOINT(3,2)=1.
  46. XPOINT(4,1)=-1.
  47. XPOINT(4,2)=1.
  48. IF(IFACE.EQ.1)THEN
  49. XPOINT(1,3)=-1.
  50. XPOINT(2,3)=-1.
  51. XPOINT(3,3)=-1.
  52. XPOINT(4,3)=-1.
  53. ZI=-1.
  54. ENDIF
  55. IF(IFACE.EQ.2)THEN
  56. XPOINT(1,3)=1.
  57. XPOINT(2,3)=1.
  58. XPOINT(3,3)=1.
  59. XPOINT(4,3)=1.
  60. ZI=+1.
  61. ENDIF
  62. IF(IFACE.NE.1.AND.IFACE.NE.2) THEN
  63. WRITE(6,*)'******** FACE NON DEFINIE DANS CUB_KSI_KP *******'
  64. STOP
  65. ENDIF
  66. C
  67. BKSI(1,1) = UNS8*(1+XPOINT(1,3)*ZI)*XPOINT(1,1)
  68. BKSI(1,2) = UNS8*(1+XPOINT(2,3)*ZI)*XPOINT(2,1)
  69. BKSI(1,3) = UNS8*(1+XPOINT(3,3)*ZI)*XPOINT(3,1)
  70. BKSI(1,4) = UNS8*(1+XPOINT(4,3)*ZI)*XPOINT(4,1)
  71. C
  72. BKSI(2,1) = UNS8*(1+XPOINT(1,3)*ZI)*XPOINT(1,2)
  73. BKSI(2,2) = UNS8*(1+XPOINT(2,3)*ZI)*XPOINT(2,2)
  74. BKSI(2,3) = UNS8*(1+XPOINT(3,3)*ZI)*XPOINT(3,2)
  75. BKSI(2,4) = UNS8*(1+XPOINT(4,3)*ZI)*XPOINT(4,2)
  76. C
  77. BKSI(3,1) = UNS8*XPOINT(1,3)
  78. BKSI(3,2) = UNS8*XPOINT(2,3)
  79. BKSI(3,3) = UNS8*XPOINT(3,3)
  80. BKSI(3,4) = UNS8*XPOINT(4,3)
  81. RETURN
  82. END
  83.  
  84.  
  85.  
  86.  

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