Télécharger la plaquette Cast3M


fr en

Télécharger asskp2.eso

Retour à la liste

Numérotation des lignes : 

asskp2

  1. C ASSKP2    SOURCE    CHAT      07/10/22    21:15:11     5921
  2.       SUBROUTINE ASSKP2(IFACE,BKP,XN,tKP)
  3.       IMPLICIT REAL*8(A-H,O-Z)
  4.       IMPLICIT INTEGER (I-N)
  5.  
  6.  
  7. C ASSEMBLAGE DE KP POUR LE SHB8
  8. *      REAL *8 BKP(3,4),XN(3),tkp(24,24),XKP(24),YKP(24),ZKP(24),B(3,8)
  9. *      REAL *8 KPTMP(24,24)
  10. *      INTEGER I,J,IFACE
  11. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  12.  
  13.       DIMENSION BKP(3,4),XN(3),tkp(24,24),XKP(24),YKP(24),ZKP(24),B(3,8)
  14.       dimension tkptmp(24,24)
  15.       IF(IFACE.NE.1.AND.IFACE.NE.2) THEN
  16.          WRITE(6,*)'******** FACE NON DEFINIE DANS ASSE_KP2 *******'
  17.          CALL ERREUR(5)
  18.          STOP
  19.       ENDIF
  20.       DO I=1,3
  21.          DO J=1,8
  22.             B(I,J)=0.
  23.          ENDDO
  24.       ENDDO
  25.       IF(IFACE.EQ.1)THEN
  26. CCCCCCCCCCCCCCC PRESSION SUR FACE INFERIEURE
  27.          DO I=1,3
  28.             DO J=1,4
  29.                B(I,J)=BKP(I,J)
  30.             ENDDO
  31.          ENDDO
  32.       ENDIF
  33.       IF(IFACE.EQ.2)THEN
  34. CCCCCCCCCCCCCCC PRESSION SUR FACE SUPERIEURE
  35.          DO I=1,3
  36.             DO J=5,8
  37.                B(I,J)=BKP(I,J-4)
  38.             ENDDO
  39.          ENDDO
  40.       ENDIF
  41.       DO I=1,24
  42.          XKP(I)=0.
  43.          YKP(I)=0.
  44.          ZKP(I)=0.
  45.       ENDDO
  46.       DO I=9,16
  47.          XKP(I)=B(2,I-8)*XN(1)-B(1,I-8)*XN(2)
  48.       ENDDO
  49.       DO I=17,24
  50.          XKP(I)=B(3,I-16)*XN(1)-B(1,I-16)*XN(3)
  51.       ENDDO
  52.       DO I=1,8
  53.          YKP(I)=B(1,I)*XN(2)-B(2,I)*XN(1)
  54.       ENDDO
  55.       DO I=17,24
  56.          YKP(I)=B(3,I-16)*XN(2)-B(2,I-16)*XN(3)
  57.       ENDDO
  58.       DO I=1,8
  59.          ZKP(I)=B(1,I)*XN(3)-B(3,I)*XN(1)
  60.       ENDDO
  61.       DO I=9,16
  62.          ZKP(I)=B(2,I-8)*XN(3)-B(3,I-8)*XN(2)
  63.       ENDDO
  64.       DO I=1,24
  65.          DO J=1,24
  66.             tkptmp(J,I)=0.
  67.          ENDDO
  68.       ENDDO
  69.       DO I=1,24
  70.          DO J=1,8
  71.             tkptmp(J,I)=XKP(I)
  72.          ENDDO
  73.          DO J=9,16
  74.             tkptmp(J,I)=YKP(I)
  75.          ENDDO
  76.          DO J=17,24
  77.             tkptmp(J,I)=ZKP(I)
  78.          ENDDO
  79.       ENDDO
  80.       DO I=1,24
  81.          DO J=1,24
  82.             tkp(I,J)=(tkptmp(I,J)+tkptmp(J,I))/2.
  83. C            KP(I,J)=KPTMP(I,J)
  84.          ENDDO
  85.       ENDDO
  86.       RETURN
  87.       END
  88.  
  89.  
  90.  

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