Télécharger la plaquette Cast3M


fr en

Télécharger cukskp.eso

Retour à la liste

Numérotation des lignes : 

cukskp

  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