Télécharger cucalb.eso

Retour à la liste

Numérotation des lignes :

  1. C CUCALB SOURCE CHAT 07/10/22 21:15:37 5921
  2. SUBROUTINE CUCALB(BKSIP,XNOE,B,AJAC)
  3. IMPLICIT REAL*8(A-H,O-Z)
  4. IMPLICIT INTEGER (I-N)
  5.  
  6. *
  7. * ------------------------------------------------------------------
  8. * CALCUL DE LA MATRICE B AU POINT D INTEGRATION DONNNE
  9. *
  10. * H. BUNG 06-98
  11. * ------------------------------------------------------------------
  12. *
  13. * ENTREES :
  14. * XNOE(24) : COORDONNEES DES NOEUDS
  15. * IPINT : NO DU POINT D INTEGRATION
  16. * MOT : NOM DE L ELEMENT FINI
  17. * SORTIE
  18. * B(NBLIB,NBN) : MATRICE B
  19. * AJAC : JACOBIEN
  20. *
  21. *
  22. DIMENSION XNOE(24),BKSIP(3,8),B(3,8)
  23.  
  24. * REAL *8 XNOE(24),BKSIP(3,8)
  25. * REAL *8 B(3,8),AJAC
  26. *
  27. *---- VARIABLES LOCALES
  28. *
  29. * REAL *8 DJ(3,3),UJ(3,3)
  30. * INTEGER LRET,I,J,K,NBN
  31. DIMENSION DJ(3,3),UJ(3,3)
  32. *
  33. * REAL *8 ZER, UN
  34. *
  35. NBN=8
  36. ZER=0D0
  37. UN=1.D0
  38. *
  39. *--- DJ = BKSIP * TRANSPOSE(XNOE)
  40. *
  41. DO I=1,3
  42. DO J=1,3
  43. DJ(I,J) = 0.
  44. END DO
  45. END DO
  46. DO I=1,3
  47. DO J=1,3
  48. DO K=1,NBN
  49. C DJ(J,I)=DJ(J,I)+BKSIP(J,K)*XNOE(I,K)
  50. DJ(J,I)=DJ(J,I)+BKSIP(J,K)*XNOE((K-1)*3+I)
  51. END DO
  52. END DO
  53. END DO
  54. *
  55. *----- UJ(J,I) MATRICE INVERSE DE DJ(J,I)
  56. *
  57.  
  58. CALL INV33(DJ,UJ,AJAC,LRET)
  59. C TEST SI ELEMENT TROP DEFORME: CROISEMENT
  60. AJAC=ABS(AJAC)
  61. *
  62. *----- MATRICE ( B ) = UJ * BKSIP
  63. *
  64. DO I=1,3
  65. DO J=1,NBN
  66. B(I,J)=0.
  67. END DO
  68. END DO
  69. DO K=1,3
  70. DO J=1,3
  71. DO I=1,NBN
  72. B(J,I)=B(J,I)+UJ(J,K)*BKSIP(K,I)
  73. END DO
  74. END DO
  75. END DO
  76. *
  77. END
  78.  
  79.  
  80.  

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