Télécharger bmlbst.eso

Retour à la liste

Numérotation des lignes :

  1. C BMLBST SOURCE CHAT 05/01/12 21:41:28 5004
  2. C BMLBST SOURCE BROC 85/07/01 23:05:03
  3. SUBROUTINE BMLBST(BLX,BLY,BLZ,BLT,DJAC,VML22,IFOU,LRE,NDDL,REL)
  4. C=======================================================================
  5. C
  6. C INTEGRE LE PRODUIT BLT ML BL
  7. C POUR LES ELEMENTS LIQUIDES
  8. C
  9. C ROUTINE FORTRAN PUR
  10. C=======================================================================
  11. C INPUT
  12. C BLX =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /X
  13. C BLY =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /Y
  14. C BLZ =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /Z
  15. C BLT =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /TETA
  16. C DJAC=JACOBIEN
  17. C VML22=-(COEFPI**2)/RHO
  18. C IFOU=IFOUR DE CCOPTIO
  19. C LRE =NOMBRE DE COLONNES DES MATRICES BL
  20. C NDDL=NOMBRE DE DEGRES DE LIBERTE
  21. C OUTPUT
  22. C REL(LRE,LRE) =MATRICE MASSE
  23. C=======================================================================
  24. IMPLICIT INTEGER(I-N)
  25. IMPLICIT REAL*8(A-H,O-Z)
  26. DIMENSION BLX(NDDL,1),BLY(NDDL,1),BLZ(NDDL,1),BLT(NDDL,1)
  27. DIMENSION REL(LRE,1)
  28. IFR=IFOU+4
  29. GOTO (666,10,10,10,30,40) ,IFR
  30. C
  31. C ELEMENTS LIQUIDES 2DIM CONT PLANES OU DEFO PLANES OU AXIS
  32. C
  33. 10 CONTINUE
  34. DO 100 IX=1,LRE
  35. DO 100 IY=1,IX
  36. REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+
  37. 1 BLY(2,IY)*BLY(2,IX))
  38. REL(IY,IX)=REL(IX,IY)
  39. 100 CONTINUE
  40. GOTO 666
  41. C
  42. C ELEMENTS LIQUIDES BIDIM FOURIER
  43. C
  44. 30 CONTINUE
  45. DO 300 IX=1,LRE
  46. DO 300 IY=1,IX
  47. REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+
  48. 1 BLY(2,IY)*BLY(2,IX)+
  49. 2 BLT(2,IY)*BLT(2,IX))
  50. REL(IY,IX)=REL(IX,IY)
  51. 300 CONTINUE
  52. GOTO 666
  53. C
  54. C ELEMENTS LIQUIDES TRIDIM
  55. C
  56. 40 CONTINUE
  57. DO 400 IX=1,LRE
  58. DO 400 IY=1,IX
  59. REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+
  60. 1 BLY(2,IY)*BLY(2,IX)+
  61. 2 BLZ(2,IY)*BLZ(2,IX))
  62. REL(IY,IX)=REL(IX,IY)
  63. 400 CONTINUE
  64. GOTO 666
  65. C
  66. 666 CONTINUE
  67. RETURN
  68. END
  69.  
  70.  

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