Télécharger asskp2.eso

Retour à la liste

Numérotation des lignes :

  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