Télécharger blmast.eso

Retour à la liste

Numérotation des lignes :

  1. C BLMAST SOURCE CHAT 05/01/12 21:40:59 5004
  2. SUBROUTINE BLMAST(IGAU,MFR,NBNO,LRE,IFOU,NDDL,NN,XEL,
  3. 1 SHPTOT,SHP,BLX,BLY,BLZ,BLT,DJAC)
  4. C=======================================================================
  5. C
  6. C CALCULE LES MATRICES BL POUR LES ELEMENTS LIQUIDES
  7. C
  8. C ROUTINE FORTRAN PUR
  9. C=======================================================================
  10. C INPUT
  11. C IGAU=NUMERO DU POINT DE GAUSS
  12. C MFR =NUMERO DE LA FORMULATION
  13. C NBNO=NOMBRE DE NOEUDS
  14. C LRE =NOMBRE DE COLONNES DE LA MATRICE B
  15. C IFOU=IFOUR DE CCOPTIO
  16. C NDDL=NOMBRE DE DEGRES DE LIBERTE
  17. C NN =NUMERO DU MODE DE FOURIER
  18. C XEL =COORDONNEES LOCALES DE L ELEMENT
  19. C SHPTOT(6,NBNO,NBGAU)=FONCTIONS DE FORMES ET DERIVEES
  20. C ZONE DE TRAVAIL
  21. C SHP(6,NBNO)=TABLEAU DE TRAVAIL
  22. C OUTPUT
  23. C DJAC=JACOBIEN
  24. C BLX = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / X
  25. C BLY = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / Y
  26. C BLZ = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / Z
  27. C BLT = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / TETA
  28. C=======================================================================
  29. IMPLICIT INTEGER(I-N)
  30. IMPLICIT REAL*8(A-H,O-Z)
  31. DIMENSION XEL(3,1),SHP(6,1),SHPTOT(6,NBNO,1)
  32. DIMENSION BLX(NDDL,1),BLY(NDDL,1),BLZ(NDDL,1),BLT(NDDL,1)
  33. CALL ZERO(BLX,NDDL,LRE)
  34. CALL ZERO(BLY,NDDL,LRE)
  35. CALL ZERO(BLZ,NDDL,LRE)
  36. CALL ZERO(BLT,NDDL,LRE)
  37. IFR=IFOU+4
  38. DIM3=1.D0
  39. RR=0.D0
  40. GOTO (666,10,10,20,30,40) ,IFR
  41. C
  42. C ELEMENTS LIQUIDES BIDIM CONTRAINTES PLANES OU DEFRMTNS PLANES
  43. C
  44. 10 CONTINUE
  45. DO 101 NP=1,NBNO
  46. SHP(2,NP)=SHPTOT(2,NP,IGAU)
  47. SHP(3,NP)=SHPTOT(3,NP,IGAU)
  48. 101 CONTINUE
  49. CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC)
  50. K=1
  51. DO 102 NP=1,NBNO
  52. BLX(2,K+1)=SHP(2,NP)
  53. BLY(2,K+1)=SHP(3,NP)
  54. 102 K=K+2
  55. GOTO 666
  56. C
  57. C ELEMENTS LIQUIDES BIDIM AXISYMETRIQUE
  58. C
  59. 20 CONTINUE
  60. DO 201 NP=1,NBNO
  61. SHP(1,NP)=SHPTOT(1,NP,IGAU)
  62. SHP(2,NP)=SHPTOT(2,NP,IGAU)
  63. SHP(3,NP)=SHPTOT(3,NP,IGAU)
  64. 201 CONTINUE
  65. CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC)
  66. K=1
  67. DO 202 NP=1,NBNO
  68. BLX(2,K+1)=SHP(2,NP)
  69. BLY(2,K+1)=SHP(3,NP)
  70. 202 K=K+2
  71. GOTO 666
  72. C
  73. C ELEMENTS LIQUIDES BIDIM FOURIER
  74. C
  75. 30 CONTINUE
  76. DO 301 NP=1,NBNO
  77. SHP(1,NP)=SHPTOT(1,NP,IGAU)
  78. SHP(2,NP)=SHPTOT(2,NP,IGAU)
  79. SHP(3,NP)=SHPTOT(3,NP,IGAU)
  80. 301 CONTINUE
  81. CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC)
  82. K=1
  83. DO 302 NP=1,NBNO
  84. BLX(2,K+1)= SHP(2,NP)
  85. BLY(2,K+1)= SHP(3,NP)
  86. BLT(2,K+1)= -SHP(1,NP)*NN/RR
  87. 302 K=K+2
  88. GOTO 666
  89. C
  90. C ELEMENTS LIQUIDES TRIDIM
  91. C
  92. 40 CONTINUE
  93. DO 401 NP=1,NBNO
  94. SHP(2,NP)=SHPTOT(2,NP,IGAU)
  95. SHP(3,NP)=SHPTOT(3,NP,IGAU)
  96. SHP(4,NP)=SHPTOT(4,NP,IGAU)
  97. 401 CONTINUE
  98. CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,3,DIM3,RR,DJAC)
  99. K=1
  100. DO 402 NP=1,NBNO
  101. BLX(2,K+1)=SHP(2,NP)
  102. BLY(2,K+1)=SHP(3,NP)
  103. BLZ(2,K+1)=SHP(4,NP)
  104. 402 K=K+2
  105. GOTO 666
  106. C
  107. 666 CONTINUE
  108. RETURN
  109. END
  110.  
  111.  
  112.  
  113.  

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